ANALYSIS/macros/train/AddMCGenPythia8.C

   1 AliGenerator* AddMCGenPythia8(Float_t e_cms = 2760., Bool_t kCR = kTRUE)
   2 {
   3   // Add Pythia 8 generator:
   4   //    - Color reconnection = ON/OFF
   5
   6   gSystem->Load("liblhapdf.so");
   7
   8   AliGenerator *genP = NULL;
   9   genP = CreatePythia8Gen(e_cms, kCR);
  10
  11   return genP;
  12 }
  13
  14 AliGenerator* CreatePythia8Gen(Float_t e_cms, Bool_t kCR) {
  15
  16    gSystem->Load("libpythia6.so");
  17    gSystem->Load("libEGPythia6.so");
  18    gSystem->Load("libAliPythia6.so");
  19    gSystem->Load("libpythia8.so");
  20    gSystem->Load("libAliPythia8.so");
  21    gSystem->Setenv("PYTHIA8DATA", gSystem->ExpandPathName("$ALICE_ROOT/PYTHIA8/pythia8175/xmldoc"));
  22    gSystem->Setenv("LHAPDF",      gSystem->ExpandPathName("$ALICE_ROOT/LHAPDF"));
  23    gSystem->Setenv("LHAPATH",     gSystem->ExpandPathName("$ALICE_ROOT/LHAPDF/PDFsets"));
  24
  25
  26   AliGenPythiaPlus* gener = new AliGenPythiaPlus(AliPythia8::Instance());
  27
  28   // set process (MB)
  29   gener->SetProcess(kPyMbDefault);
  30
  31   //   Centre of mass energy
  32   gener->SetEnergyCMS(e_cms); // in GeV
  33
  34   // Event list
  35   gener->SetEventListRange(-1, 2);
  36
  37   // color reconnection
  38   (AliPythia8::Instance())->ReadString("Tune:pp = 5");//CR
  39
  40   if(kCR)
  41     (AliPythia8::Instance())->ReadString("BeamRemnants:reconnectColours = on");
  42   else
  43     (AliPythia8::Instance())->ReadString("BeamRemnants:reconnectColours = off");
  44
  45   return gener;
  46 }