AliGenerator* CreateGenerator(); void fastGenDPMjet(Int_t nev = 1, char* filename = "galice.root") { // Runloader gSystem->Load("libdpmjet.so"); gSystem->Load("libTDPMjet.so"); AliRunLoader* rl = AliRunLoader::Open("galice.root","FASTRUN","recreate"); rl->SetCompressionLevel(2); rl->SetNumberOfEventsPerFile(nev); rl->LoadKinematics("RECREATE"); rl->MakeTree("E"); gAlice->SetRunLoader(rl); // Create stack rl->MakeStack(); AliStack* stack = rl->Stack(); // Header AliHeader* header = rl->GetHeader(); // // Create and Initialize Generator AliGenerator *gener = CreateGenerator(); gener->Init(); gener->SetStack(stack); // // Event Loop // Int_t iev; for (iev = 0; iev < nev; iev++) { printf("\n \n Event number %d \n \n", iev); // Initialize event header->Reset(0,iev); rl->SetEventNumber(iev); stack->Reset(); rl->MakeTree("K"); // stack->ConnectTree(); // Generate event gener->Generate(); // Analysis Int_t npart = stack->GetNprimary(); printf("Analyse %d Particles\n", npart); for (Int_t part=0; partParticle(part); Int_t pdg = MPart->GetPdgCode(); Int_t pis = MPart->GetStatusCode(); // printf("Particle %5d %5d\n", part, pdg, pis); } // Finish event header->SetNprimary(stack->GetNprimary()); header->SetNtrack(stack->GetNtrack()); // I/O // stack->FinishEvent(); header->SetStack(stack); rl->TreeE()->Fill(); rl->WriteKinematics("OVERWRITE"); } // event loop // // Termination // Generator gener->FinishRun(); // Write file rl->WriteHeader("OVERWRITE"); gener->Write(); rl->Write(); } AliGenerator* CreateGenerator() { // kDpmMb // kDpmMbNonDiffr // kDpmDiffr // kDpmSingleDiffr // kDpmDoubleDiffr gener = new AliGenDPMjet(1); gener->SetProcess(kDpmMb); gener->SetEnergyCMS(14000.); return gener; }