1 AliGenerator* CreateGenerator();
3 void fastGenDPMjet(Int_t nev = 1, char* filename = "galice.root")
7 gSystem->Load("libdpmjet.so");
8 gSystem->Load("libTDPMjet.so");
10 AliRunLoader* rl = AliRunLoader::Open("galice.root","FASTRUN","recreate");
12 rl->SetCompressionLevel(2);
13 rl->SetNumberOfEventsPerFile(nev);
14 rl->LoadKinematics("RECREATE");
16 gAlice->SetRunLoader(rl);
20 AliStack* stack = rl->Stack();
23 AliHeader* header = rl->GetHeader();
25 // Create and Initialize Generator
26 AliGenerator *gener = CreateGenerator();
28 gener->SetStack(stack);
35 for (iev = 0; iev < nev; iev++) {
37 printf("\n \n Event number %d \n \n", iev);
41 rl->SetEventNumber(iev);
44 // stack->ConnectTree();
49 Int_t npart = stack->GetNprimary();
50 printf("Analyse %d Particles\n", npart);
51 for (Int_t part=0; part<npart; part++) {
52 TParticle *MPart = stack->Particle(part);
53 Int_t pdg = MPart->GetPdgCode();
54 Int_t pis = MPart->GetStatusCode();
55 // printf("Particle %5d %5d\n", part, pdg, pis);
59 header->SetNprimary(stack->GetNprimary());
60 header->SetNtrack(stack->GetNtrack());
64 header->SetStack(stack);
66 rl->WriteKinematics("OVERWRITE");
74 rl->WriteHeader("OVERWRITE");
81 AliGenerator* CreateGenerator()
88 gener = new AliGenDPMjet(1);
89 gener->SetProcess(kDpmMb);
90 gener->SetEnergyCMS(14000.);