1 AliGenerator* CreateGenerator();
3 void fastGen(Int_t nev = 1, char* filename = "galice.root")
6 gSystem->Load("liblhapdf");
7 gSystem->Load("libpythia6");
8 gSystem->Load("libEGPythia6");
9 gSystem->Load("libAliPythia6");
13 AliRunLoader* rl = AliRunLoader::Open("galice.root","FASTRUN","recreate");
15 rl->SetCompressionLevel(2);
16 rl->SetNumberOfEventsPerFile(nev);
17 rl->LoadKinematics("RECREATE");
19 gAlice->SetRunLoader(rl);
23 AliStack* stack = rl->Stack();
26 AliHeader* header = rl->GetHeader();
28 // Create and Initialize Generator
29 AliGenerator *gener = CreateGenerator();
31 gener->SetStack(stack);
38 for (iev = 0; iev < nev; iev++) {
40 printf("\n \n Event number %d \n \n", iev);
44 rl->SetEventNumber(iev);
47 // stack->ConnectTree();
52 Int_t npart = stack->GetNprimary();
53 printf("Analyse %d Particles\n", npart);
54 for (Int_t part=0; part<npart; part++) {
55 TParticle *MPart = stack->Particle(part);
56 Int_t mpart = MPart->GetPdgCode();
57 printf("Particle %d\n", mpart);
61 header->SetNprimary(stack->GetNprimary());
62 header->SetNtrack(stack->GetNtrack());
66 header->SetStack(stack);
68 rl->WriteKinematics("OVERWRITE");
76 rl->WriteHeader("OVERWRITE");
83 AliGenerator* CreateGenerator()
85 gener = new AliGenPythia(1);
88 // vertex position and smearing
89 gener->SetVertexSmear(kPerEvent);
91 gener->SetStrucFunc(kCTEQ6);
92 // charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
93 gener->SetProcess(kPyJets);
94 // Centre of mass energy
95 gener->SetEnergyCMS(5500.);
96 // Pt transfer of the hard scattering
97 gener->SetPtHard(50.,50.2);
98 // Initialize generator