ANALYSIS/macros/train/AddMCGenAmpt.C

   1 AliGenerator *AddMCGenAmpt()
   2 {
   3 // User defined generator
   4
   5   gSystem->Load("libampt.so ");
   6   gSystem->Load("libTAmpt.so");
   7
   8   AliGenAmpt *genAMPT = new AliGenAmpt(-1);
   9
  10   // will be made optional later
  11   genAMPT->SetEnergyCMS(2760);
  12   genAMPT->SetReferenceFrame("CMS");
  13   genAMPT->SetProjectile("A", 208, 82);
  14   genAMPT->SetTarget    ("A", 208, 82);
  15   genAMPT->SetPtHardMin (2);
  16   genAMPT->SetImpactParameterRange(0.00,20.00);
  17   genAMPT->SetJetQuenching(0); // enable jet quenching
  18   genAMPT->SetShadowing(1);    // enable shadowing
  19   genAMPT->SetDecaysOff(1);    // neutral pion and heavy particle decays switched off
  20   genAMPT->SetSpectators(0);   // track spectators
  21   genAMPT->SetIsoft(4);        // 4=string melting, 1=standard AMPT
  22   genAMPT->SetXmu(3.2264);     // parton xsection
  23   genAMPT->SetNtMax(150);      // time bins
  24
  25   genAMPT->SetAlpha(1./3.);    //alpha =0.333
  26   genAMPT->SetStringFrag(0.5,0.9); //string fragmentation parameters
  27   genAMPT->SetIpop(1); //enable popcorn mechanism (net-baryon stopping)
  28   // This particular choice of gives scattering cross section to be 1.5 mb
  29
  30   genAMPT->SetRandomReactionPlane(kTRUE);
  31
  32  return genAMPT;
  33
  34
  35
  36
  37 }