+++ /dev/null
-void sim(Int_t nev=100) {\r
-// AliLog::SetGlobalLogLevel(AliLog::kError);\r
-\r
- AliSimulation simu;\r
- simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");\r
- simu.SetMakeDigitsFromHits("ITS TPC");\r
-\r
- simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");\r
-\r
-
- simu.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd()));
-
- TStopwatch timer;\r
- timer.Start();\r
- simu.Run(nev);\r
- WriteXsection();
- timer.Stop();\r
- timer.Print();\r
-}\r
-
-WriteXsection()
-{
- TPythia6 *pythia = TPythia6::Instance();
- pythia->Pystat(1);
- Double_t xsection = pythia->GetPARI(1);
- Int_t ntrials = pythia->GetMSTI(5);
-
- TFile *file = new TFile("pyxsec.root","recreate");
- TTree *tree = new TTree("Xsection","Pythia cross section");
- TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
- TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
- tree->Fill();
-
- tree->Write();
- file->Close();
-
- cout << "Pythia cross section: " << xsection
- << ", number of trials: " << ntrials << endl;
-}