fix macro control plots npart eccentricity
[u/mrichter/AliRoot.git] / PWG2 / FLOW / Documentation / examples / runGlauberMCexample.C
1 {
2   //load libraries
3  gSystem->Load("libPWG2flowTools");
4  //  gSystem->SetBuildDir("/tmp");
5  //  gROOT->LoadMacro("$ALICE_ROOT/PWG2/FLOW/AliFlowTools/glauberMC/AliGlauberNucleon.cxx+");
6  //  gROOT->LoadMacro("$ALICE_ROOT/PWG2/FLOW/AliFlowTools/glauberMC/AliGlauberNucleus.cxx+");
7  //  gROOT->LoadMacro("$ALICE_ROOT/PWG2/FLOW/AliFlowTools/glauberMC/AliGlauberMC.cxx+");
8
9   //set the random seed from current time
10   TTimeStamp time;
11   Int_t seed = time->GetSec();
12   gRandom->SetSeed(seed);
13
14   Int_t nevents = 10000; // number of events to simulate 
15   // supported systems are e.g. "p", "d", "Si", "Au", "Pb", "U" 
16   Option_t *sysA="Pb"; 
17   Option_t *sysB="Pb";
18   Double_t signn=72; // inelastic nucleon nucleon cross section
19   //const char *fname="GlauberMC_PbPb_ntuple.root"; // name output file
20
21   // run the code to produce an ntuple:
22   //  AliGlauberMC::runAndSaveNucleons(10000,"Pb","Pb",72);
23   Double_t mind=0.4;
24   AliGlauberMC::runAndSaveNtuple(nevents,sysA,sysB,signn,mind);
25
26 }