1 AliGenerator* CreateGenerator();
3 void fastGenDPMjet(Int_t nev = 1, char* filename = "galice.root")
6 gSystem->Load("liblhapdf");
7 gSystem->Load("libpythia6");
8 gSystem->Load("libdpmjet.so");
9 gSystem->Load("libTDPMjet.so");
11 AliRunLoader* rl = AliRunLoader::Open("galice.root","FASTRUN","recreate");
12 AliPythiaRndm::SetPythiaRandom(new TRandom());
14 rl->SetCompressionLevel(2);
15 rl->SetNumberOfEventsPerFile(nev);
16 rl->LoadKinematics("RECREATE");
18 gAlice->SetRunLoader(rl);
22 AliStack* stack = rl->Stack();
25 AliHeader* header = rl->GetHeader();
27 // Create and Initialize Generator
28 AliGenerator *gener = CreateGenerator();
30 gener->SetStack(stack);
37 for (iev = 0; iev < nev; iev++) {
39 printf("\n \n Event number %d \n \n", iev);
43 rl->SetEventNumber(iev);
46 // stack->ConnectTree();
51 Int_t npart = stack->GetNprimary();
52 printf("Analyse %d Particles\n", npart);
53 for (Int_t part=0; part<npart; part++) {
54 TParticle *MPart = stack->Particle(part);
55 Int_t pdg = MPart->GetPdgCode();
56 Int_t pis = MPart->GetStatusCode();
57 // printf("Particle %5d %5d\n", part, pdg, pis);
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()
90 gener = new AliGenDPMjet(1);
91 gener->SetProcess(kDpmMb);
92 gener->SetProjectile("P", 1, 1);
93 gener->SetTarget("P", 1, 1);
94 gener->SetEnergyCMS(14000.);