1 AliGenerator* CreateGenerator();
3 void fastGen(Int_t nev = 1, char* filename = "galice.root")
7 AliRunLoader* rl = AliRunLoader::Open("galice.root","FASTRUN","recreate");
9 rl->SetCompressionLevel(2);
10 rl->SetNumberOfEventsPerFile(nev);
11 rl->LoadKinematics("RECREATE");
13 gAlice->SetRunLoader(rl);
17 AliStack* stack = rl->Stack();
20 AliHeader* header = rl->GetHeader();
22 // Create and Initialize Generator
23 AliGenerator *gener = CreateGenerator();
25 gener->SetStack(stack);
32 for (iev = 0; iev < nev; iev++) {
34 printf("\n \n Event number %d \n \n", iev);
38 rl->SetEventNumber(iev);
41 // stack->ConnectTree();
46 Int_t npart = stack->GetNprimary();
47 printf("Analyse %d Particles\n", npart);
48 for (Int_t part=0; part<npart; part++) {
49 TParticle *MPart = stack->Particle(part);
50 Int_t mpart = MPart->GetPdgCode();
51 printf("Particle %d\n", mpart);
55 header->SetNprimary(stack->GetNprimary());
56 header->SetNtrack(stack->GetNtrack());
60 header->SetStack(stack);
62 rl->WriteKinematics("OVERWRITE");
70 rl->WriteHeader("OVERWRITE");
77 AliGenerator* CreateGenerator()
79 gener = new AliGenPythia(1);
82 // vertex position and smearing
83 gener->SetVertexSmear(kPerEvent);
85 gener->SetStrucFunc(kGRVHO);
86 // charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
87 gener->SetProcess(kPyJets);
88 // Centre of mass energy
89 gener->SetEnergyCMS(5500.);
90 // Pt transfer of the hard scattering
91 gener->SetPtHard(50.,50.2);
92 // Initialize generator