fix macro control plots npart eccentricity
[u/mrichter/AliRoot.git] / PWG2 / FLOW / Documentation / examples / runGlauberMCexample.C
CommitLineData
ec852657 1{
2 //load libraries
e06729d3 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+");
ec852657 8
a0278ee7 9 //set the random seed from current time
10 TTimeStamp time;
11 Int_t seed = time->GetSec();
12 gRandom->SetSeed(seed);
13
ba441bc7 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
74d0740e 19 //const char *fname="GlauberMC_PbPb_ntuple.root"; // name output file
ba441bc7 20
21 // run the code to produce an ntuple:
22 // AliGlauberMC::runAndSaveNucleons(10000,"Pb","Pb",72);
23 Double_t mind=0.4;
74d0740e 24 AliGlauberMC::runAndSaveNtuple(nevents,sysA,sysB,signn,mind);
ba441bc7 25
ec852657 26}