1 AliGenerator* CreateGenerator();
3 void fastGen(Int_t nev = 1, char* filename = "galice.root")
9 TFile* file = new TFile(filename, "recreate");
11 AliStack* stack = new AliStack(10000);
12 stack->MakeTree(0, filename);
15 AliHeader* header = new AliHeader();
17 TTree* treeE = new TTree("TE","Headers");
18 treeE->Branch("Header", "AliHeader", &header, 4000, 0);
21 // Create and Initialize Generator
22 AliGenerator *gener = CreateGenerator();
24 gener->SetStack(stack);
31 for (iev = 0; iev < nev; iev++) {
33 printf("\n \n Event number %d \n \n", iev);
38 stack->BeginEvent(iev);
43 Int_t npart = stack->GetNprimary();
44 printf("Analyse %d Particles\n", npart);
45 for (Int_t part=0; part<npart; part++) {
46 TParticle *MPart = stack->Particle(part);
47 Int_t mpart = MPart->GetPdgCode();
48 printf("Particle %d\n", mpart);
52 header->SetNprimary(stack->GetNprimary());
53 header->SetNtrack(stack->GetNtrack());
57 header->SetStack(stack);
59 (stack->TreeK())->Write(0,TObject::kOverwrite);
66 treeE->Write(0,TObject::kOverwrite);
67 delete treeE; treeE = 0;
76 AliGenerator* CreateGenerator()
78 gener = new AliGenPythia(1);
81 // vertex position and smearing
82 gener->SetVertexSmear(kPerEvent);
84 gener->SetStrucFunc(kGRVHO);
85 // charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
86 gener->SetProcess(kPyJets);
87 // Centre of mass energy
88 gener->SetEnergyCMS(5500.);
89 // Pt transfer of the hard scattering
90 gener->SetPtHard(50.,50.2);
91 // Initialize generator