[u/mrichter/AliRoot.git] / prod / LHC09a14 / sim.C

void sim(Int_t nev=100) {\r
//  AliLog::SetGlobalLogLevel(AliLog::kError);\r
\r
  AliSimulation simu;\r
  simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");\r
  simu.SetMakeDigitsFromHits("ITS TPC");\r
\r
  simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");\r
\r

  simu.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd()));

  TStopwatch timer;\r
  timer.Start();\r
  simu.Run(nev);\r
  WriteXsection();
  timer.Stop();\r
  timer.Print();\r
}\r

WriteXsection()
{
  TPythia6 *pythia = TPythia6::Instance();
  pythia->Pystat(1);
  Double_t xsection = pythia->GetPARI(1);
  Int_t    ntrials  = pythia->GetMSTI(5);

  TFile *file = new TFile("pyxsec.root","recreate");
  TTree   *tree   = new TTree("Xsection","Pythia cross section");
  TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
  TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
  tree->Fill();

  tree->Write();
  file->Close();

  cout << "Pythia cross section: " << xsection 
       << ", number of trials: " << ntrials << endl;
}
Commit	Line	Data
b0491165	1	void sim(Int_t nev=100) {\r
	2	// AliLog::SetGlobalLogLevel(AliLog::kError);\r
	3	\r
	4	AliSimulation simu;\r
	5	simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");\r
	6	simu.SetMakeDigitsFromHits("ITS TPC");\r
	7	\r
	8	simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");\r
	9	\r
	10
	11	simu.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd()));
	12
	13	TStopwatch timer;\r
	14	timer.Start();\r
	15	simu.Run(nev);\r
	16	WriteXsection();
	17	timer.Stop();\r
	18	timer.Print();\r
	19	}\r
	20
	21	WriteXsection()
	22	{
	23	TPythia6 *pythia = TPythia6::Instance();
	24	pythia->Pystat(1);
	25	Double_t xsection = pythia->GetPARI(1);
	26	Int_t ntrials = pythia->GetMSTI(5);
	27
	28	TFile *file = new TFile("pyxsec.root","recreate");
	29	TTree *tree = new TTree("Xsection","Pythia cross section");
	30	TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
	31	TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
	32	tree->Fill();
	33
	34	tree->Write();
	35	file->Close();
	36
	37	cout << "Pythia cross section: " << xsection
	38	<< ", number of trials: " << ntrials << endl;
	39	}