1 void sim(Int_t nev=1) {
3 simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD PMD T0 ZDC VZERO");
4 simu.SetMakeDigits ("TRD TOF PHOS HMPID EMCAL MUON FMD PMD T0 ZDC VZERO");
5 simu.SetMakeDigitsFromHits("ITS TPC");
6 simu.SetWriteRawData("ALL","raw.root",kTRUE);
7 simu.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
8 simu.SetSpecificStorage("GRP/GRP/Data",
9 Form("local://%s",gSystem->pwd()));
11 simu.SetRunQA("ALL:ALL") ;
12 simu.SetQARefDefaultStorage("local://$ALICE_ROOT/OCDB") ;
14 for (Int_t det = 0 ; det < AliQA::kNDET ; det++) {
15 simu.SetQACycles(det, 2) ;
28 TPythia6 *pythia = TPythia6::Instance();
30 Double_t xsection = pythia->GetPARI(1);
31 Int_t ntrials = pythia->GetMSTI(5);
33 TTree *tree = new TTree("Xsection","Pythia cross section");
34 TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
35 TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
38 TFile *file = new TFile("pyxsec.root","recreate");
42 cout << "Pythia cross section: " << xsection
43 << ", number of trials: " << ntrials << endl;