test/QA/sim.C

   1 void sim(Int_t nev=1) {
   2   const char * kYear = "08" ;
   3   AliSimulation simu;
   4   simu.SetMakeSDigits("TRD TOF PHOS HMPID  EMCAL MUON FMD PMD T0 ZDC VZERO");
   5   simu.SetMakeDigits ("TRD TOF PHOS HMPID  EMCAL MUON FMD PMD T0 ZDC VZERO");
   6   simu.SetMakeDigitsFromHits("ITS TPC");
   7   simu.SetWriteRawData("ALL","raw.root",kTRUE);
   8 //  simu.SetDefaultStorage("alien://Folder=/alice/data/2008/LHC08c/OCDB/");
   9   simu.SetDefaultStorage("local://$ALICE_ROOT");
  10   simu.SetSpecificStorage("EMCAL/*","local://DB");
  11
  12   simu.SetRunQA("ALL:ALL") ;
  13  // AliQA::SetQARefStorage(Form("%s%s/", AliQA::GetQARefDefaultStorage(), kYear)) ;
  14   AliQA::SetQARefStorage("local://$ALICE_ROOT") ;
  15   // AliQA::SetQARefDataDirName(AliQA::kMONTECARLO) ; //RUN_TYPE
  16
  17   TStopwatch timer;
  18   timer.Start();
  19   simu.Run(nev);
  20   WriteXsection();
  21   timer.Stop();
  22   timer.Print();
  23 }
  24
  25 WriteXsection()
  26 {
  27   TPythia6 *pythia = TPythia6::Instance();
  28   pythia->Pystat(1);
  29   Double_t xsection = pythia->GetPARI(1);
  30   Int_t    ntrials  = pythia->GetMSTI(5);
  31
  32   TTree   *tree   = new TTree("Xsection","Pythia cross section");
  33   TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
  34   TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
  35   tree->Fill();
  36
  37   TFile *file = new TFile("pyxsec.root","recreate");
  38   tree->Write();
  39   file->Close();
  40
  41   cout << "Pythia cross section: " << xsection
  42        << ", number of trials: " << ntrials << endl;
  43 }
  44