promoted 2 methods to public needed by AMORE

[u/mrichter/AliRoot.git] / prod / LHC09a1 / sim.C

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