Working on the electron cut
[u/mrichter/AliRoot.git] / prod / LHC09a1 / sim.C
1 void sim(Int_t nev=150) {\r
2   AliSimulation simu;\r
3   simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC T0 VZERO");\r
4   simu.SetMakeDigitsFromHits("ITS TPC");\r
5 \r
6   simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");\r
7   \r
8   // No write access to the OCDB => local specific storage\r
9   simu.SetSpecificStorage("GRP/GRP/Data",\r
10                           Form("local://%s/",gSystem->pwd()));\r
11 \r
12   Printf("%s:%d %d events",(char*)__FILE__,__LINE__,nev);\r
13 \r
14   TStopwatch timer;\r
15   timer.Start();\r
16   simu.Run(nev);\r
17   WriteXsection();\r
18   timer.Stop();\r
19   timer.Print();\r
20 }\r
21 \r
22 WriteXsection()\r
23 {\r
24  TPythia6 *pythia = TPythia6::Instance();\r
25  pythia->Pystat(1);\r
26  Double_t xsection = pythia->GetPARI(1);\r
27  UInt_t    ntrials  = pythia->GetMSTI(5);\r
28 \r
29  TFile *file = new TFile("pyxsec.root","recreate");\r
30  TTree   *tree   = new TTree("Xsection","Pythia cross section");\r
31  TBranch *branch = tree->Branch("xsection", &xsection, "X/D");\r
32  TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");\r
33  tree->Fill();\r
34 \r
35 \r
36 \r
37  tree->Write();\r
38  file->Close();\r
39 \r
40  cout << "Pythia cross section: " << xsection\r
41       << ", number of trials: " << ntrials << endl;\r
42 }\r
43 \r