changing run periof
[u/mrichter/AliRoot.git] / test / QA / sim.C
1 void sim(Int_t nev=4) {
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/LHC08d/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   for (Int_t det = 0 ; det < AliQA::kNDET ; det++) {
17     simu.SetQACycles(det, 2) ;
18   }
19   
20   TStopwatch timer;
21   timer.Start();
22   simu.Run(nev);
23   WriteXsection();
24   timer.Stop();
25   timer.Print();
26 }
27
28 WriteXsection()
29 {
30   TPythia6 *pythia = TPythia6::Instance();
31   pythia->Pystat(1);
32   Double_t xsection = pythia->GetPARI(1);
33   Int_t    ntrials  = pythia->GetMSTI(5);
34
35   TTree   *tree   = new TTree("Xsection","Pythia cross section");
36   TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
37   TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
38   tree->Fill();
39
40   TFile *file = new TFile("pyxsec.root","recreate");
41   tree->Write();
42   file->Close();
43
44   cout << "Pythia cross section: " << xsection 
45        << ", number of trials: " << ntrials << endl;
46 }
47