PWGCF/FLOW/Documentation/examples/runGlauberMCexample.C

   1 {
   2   //load libraries
   3  gSystem->Load("libPWGGlauber");
   4  //  gSystem->SetBuildDir("/tmp");
   5  //  gROOT->LoadMacro("$ALICE_ROOT/PWG/Glauber/AliGlauberNucleon.cxx+");
   6  //  gROOT->LoadMacro("$ALICE_ROOT/PWG/Glauber/AliGlauberNucleus.cxx+");
   7  //  gROOT->LoadMacro("$ALICE_ROOT/PWG/Glauber/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 = 1000; // 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=64; // 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 }