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

void sim(Int_t nev=10) {
  AliSimulation simu;
  //simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL FMD PMD T0 ZDC VZERO MUON HLT");
  //simu.SetMakeDigits ("TRD TOF PHOS HMPID EMCAL FMD PMD T0 ZDC VZERO MUON HLT");
  //simu.SetMakeDigitsFromHits("ITS TPC");
  //simu.SetWriteRawData("ALL","raw.root",kTRUE);
  simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/?cacheFold=/tmp/CDBCache?operateDisconnected=kFALSE");

  // No write access to the OCDB => specific storage
  simu.SetSpecificStorage("GRP/GRP/Data",
                          Form("local://%s",gSystem->pwd()));

  //simu.SetRunQA(":") ; 
  //simu.SetRunHLT("");
  //simu.SetWriteGRPEntry(kFALSE);

  //Merge with HIJING
  //Load necessary libraries and connect to the GRID
  gSystem->Load("libNetx.so") ;
  gSystem->Load("libRAliEn.so");
  TGrid::Connect("alien://") ;
    
  //Feed Grid with collection file
  TGridCollection * collection = (TGridCollection*) TAlienCollection::Open("collection.xml");   
  if (! collection) {
    AliError(Form("collection: %s not found", kXML)) ;
    return kFALSE ;
  }
  TGridResult* result = collection->GetGridResult("",0 ,0);
  TString alienURL = result->GetKey(0, "turl") ;//split all files in collection, only one is sent per job
  cout << "================== " << alienURL << endl ;
  simu.MergeWith(alienURL);

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

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

  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();

  TFile *file = new TFile("pyxsec.root","recreate");
  tree->Write();
  file->Close();

  cout << "Pythia cross section: " << xsection 
       << ", number of trials: " << ntrials << endl;
}
Commit	Line	Data
333b9f18	1	void sim(Int_t nev=10) {
	2	AliSimulation simu;
	3	//simu.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL FMD PMD T0 ZDC VZERO MUON HLT");
	4	//simu.SetMakeDigits ("TRD TOF PHOS HMPID EMCAL FMD PMD T0 ZDC VZERO MUON HLT");
	5	//simu.SetMakeDigitsFromHits("ITS TPC");
	6	//simu.SetWriteRawData("ALL","raw.root",kTRUE);
	7	simu.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/?cacheFold=/tmp/CDBCache?operateDisconnected=kFALSE");
	8
	9	// No write access to the OCDB => specific storage
	10	simu.SetSpecificStorage("GRP/GRP/Data",
	11	Form("local://%s",gSystem->pwd()));
	12
	13	//simu.SetRunQA(":") ;
	14	//simu.SetRunHLT("");
	15	//simu.SetWriteGRPEntry(kFALSE);
	16
	17	//Merge with HIJING
	18	//Load necessary libraries and connect to the GRID
	19	gSystem->Load("libNetx.so") ;
	20	gSystem->Load("libRAliEn.so");
	21	TGrid::Connect("alien://") ;
	22
	23	//Feed Grid with collection file
	24	TGridCollection * collection = (TGridCollection*) TAlienCollection::Open("collection.xml");
	25	if (! collection) {
	26	AliError(Form("collection: %s not found", kXML)) ;
	27	return kFALSE ;
	28	}
	29	TGridResult* result = collection->GetGridResult("",0 ,0);
	30	TString alienURL = result->GetKey(0, "turl") ;//split all files in collection, only one is sent per job
	31	cout << "================== " << alienURL << endl ;
	32	simu.MergeWith(alienURL);
	33
	34	TStopwatch timer;
	35	timer.Start();
	36	simu.Run(nev);
	37	WriteXsection();
	38	timer.Stop();
	39	timer.Print();
	40	}
	41
	42	WriteXsection()
	43	{
	44	TPythia6 *pythia = TPythia6::Instance();
	45	pythia->Pystat(1);
	46	Double_t xsection = pythia->GetPARI(1);
	47	Int_t ntrials = pythia->GetMSTI(5);
	48
	49	TTree *tree = new TTree("Xsection","Pythia cross section");
	50	TBranch *branch = tree->Branch("xsection", &xsection, "X/D");
	51	TBranch *branch = tree->Branch("ntrials" , &ntrials , "X/i");
	52	tree->Fill();
	53
	54	TFile *file = new TFile("pyxsec.root","recreate");
	55	tree->Write();
	56	file->Close();
	57
	58	cout << "Pythia cross section: " << xsection
	59	<< ", number of trials: " << ntrials << endl;
	60	}