set proper OCD location
[u/mrichter/AliRoot.git] / 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