1 void sim(Int_t nev=100) {
\r
2 // AliLog::SetGlobalLogLevel(AliLog::kError);
\r
5 simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");
\r
6 simu.SetMakeDigitsFromHits("ITS TPC");
\r
8 simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");
\r
11 simu.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd()));
23 TPythia6 *pythia = TPythia6::Instance();
25 Double_t xsection = pythia->GetPARI(1);
26 Int_t ntrials = pythia->GetMSTI(5);
28 TFile *file = new TFile("pyxsec.root","recreate");
29 TTree *tree = new TTree("Xsection","Pythia cross section");
30 TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
31 TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
37 cout << "Pythia cross section: " << xsection
38 << ", number of trials: " << ntrials << endl;