1 AliGenerator* CreateGenerator();
3 void fastGenDPMjet(Int_t nev = 1, char* filename = "dpmjet.root")
7 gSystem->Load("liblhapdf");
8 gSystem->Load("libpythia6");
9 gSystem->Load("libdpmjet.so");
10 gSystem->Load("libTDPMjet.so");
13 AliPDG::AddParticlesToPdgDataBase();
15 AliRunLoader* rl = AliRunLoader::Open(filename,"FASTRUN","recreate");
16 AliPythiaRndm::SetPythiaRandom(new TRandom());
18 rl->SetCompressionLevel(2);
19 rl->SetNumberOfEventsPerFile(nev);
20 rl->LoadKinematics("RECREATE");
22 gAlice->SetRunLoader(rl);
26 AliStack* stack = rl->Stack();
29 AliHeader* header = rl->GetHeader();
30 //AliGenDPMjetEventHEader* dpmHeader;
32 // Create and Initialize Generator
33 AliGenDPMjet *gener = CreateGenerator();
34 //AliCollisionGeometry *coll = gener->CollisionGeometry();
36 gener->SetStack(stack);
41 for (iev = 0; iev < nev; iev++) {
43 if(iev%500==0) printf("\n Event number %d \n", iev);
47 rl->SetEventNumber(iev);
54 Int_t npart = stack->GetNprimary();
57 header->SetNprimary(stack->GetNprimary());
58 header->SetNtrack(stack->GetNtrack());
63 header->SetStack(stack);
65 rl->WriteKinematics("OVERWRITE");
71 rl->WriteHeader("OVERWRITE");
78 AliGenerator* CreateGenerator()
80 AliGenDPMjet *gener = new AliGenDPMjet(-1);
83 /*Float_t pEnergy = 4000.;
84 gener->SetProjectile("P", 1, 1);
85 gener->SetTarget("A", 208, 82);
87 gener->SetEnergyCMS(TMath::Sqrt(82./208.) * 2* pEnergy);
88 gener->SetProjectileBeamEnergy(pEnergy); */
91 Float_t pEnergy = 1577.;
92 gener->SetProjectile("A", 208, 82);
93 gener->SetTarget("P", 1, 1);
95 gener->SetEnergyCMS(TMath::Sqrt(208./82.) * 2* pEnergy);
96 gener->SetProjectileBeamEnergy(pEnergy);
99 /* gener->SetProjectile("A", 208, 82);
100 gener->SetTarget("A", 208, 82);
101 gener->SetEnergyCMS(2760.);
103 gener->SetProcess(kDpmMb);
104 gener->SetImpactParameterRange(0., 100.);
105 //gener->SetFragmentProd(kTRUE);*/
107 gener->SetTrackingFlag(0);