set proper OCD location
[u/mrichter/AliRoot.git] / test / QA / sim.C
CommitLineData
478d5b53 1void 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);
59932366 8// simu.SetDefaultStorage("alien://Folder=/alice/data/2008/LHC08c/OCDB/");
f66339d1 9 simu.SetDefaultStorage("local://$ALICE_ROOT");
478d5b53 10 simu.SetSpecificStorage("EMCAL/*","local://DB");
11
f66339d1 12 simu.SetRunQA("ALL:ALL") ;
51757634 13 // AliQA::SetQARefStorage(Form("%s%s/", AliQA::GetQARefDefaultStorage(), kYear)) ;
14 AliQA::SetQARefStorage("local://$ALICE_ROOT") ;
80336a64 15 // AliQA::SetQARefDataDirName(AliQA::kMONTECARLO) ; //RUN_TYPE
478d5b53 16
17 TStopwatch timer;
18 timer.Start();
19 simu.Run(nev);
20 WriteXsection();
21 timer.Stop();
22 timer.Print();
23}
24
25WriteXsection()
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