prod/LHC08d8/sim.C

   1 void sim(Int_t nev=100) {\r
   2 //  AliLog::SetGlobalLogLevel(AliLog::kError);\r
   3 \r
   4   AliSimulation simu;\r
   5   simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");\r
   6   simu.SetMakeDigitsFromHits("ITS TPC");\r
   7 \r
   8   simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");\r
   9 \r
  10
  11   simu.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd()));
  12
  13   TStopwatch timer;\r
  14   timer.Start();\r
  15   simu.Run(nev);\r
  16   WriteXsection();
  17   timer.Stop();\r
  18   timer.Print();\r
  19 }\r
  20
  21 WriteXsection()
  22 {
  23   TPythia6 *pythia = TPythia6::Instance();
  24   pythia->Pystat(1);
  25   Double_t xsection = pythia->GetPARI(1);
  26   Int_t    ntrials  = pythia->GetMSTI(5);
  27
  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");
  32   tree->Fill();
  33
  34   tree->Write();
  35   file->Close();
  36
  37   cout << "Pythia cross section: " << xsection
  38        << ", number of trials: " << ntrials << endl;
  39 }