PWG2/FLOW/Documentation/examples/runGlauberMCexample.C

   1 {
   2   //load libraries
   3   gSystem->SetBuildDir("/tmp");
   4   gROOT->LoadMacro("$ALICE_ROOT/PWG2/FLOW/AliFlowTools/glauberMC/AliGlauberNucleon.cxx+");
   5   gROOT->LoadMacro("$ALICE_ROOT/PWG2/FLOW/AliFlowTools/glauberMC/AliGlauberNucleus.cxx+");
   6   gROOT->LoadMacro("$ALICE_ROOT/PWG2/FLOW/AliFlowTools/glauberMC/AliGlauberMC.cxx+");
   7
   8
   9   Int_t nevents = 10000; // number of events to simulate
  10   // supported systems are e.g. "p", "d", "Si", "Au", "Pb", "U"
  11   Option_t *sysA="Pb";
  12   Option_t *sysB="Pb";
  13   Double_t signn=72; // inelastic nucleon nucleon cross section
  14   const char *fname="GlauberMC_PbPb_ntuple.root"; // name output file
  15
  16   // run the code to produce an ntuple:
  17   //  AliGlauberMC::runAndSaveNucleons(10000,"Pb","Pb",72);
  18   Double_t mind=0.4;
  19   AliGlauberMC::runAndSaveNtuple(nevents,sysA,sysB,signn,mind,fname);
  20
  21 }