Update AMPT macro:
authormiweber <m.weber@cern.ch>
Tue, 28 Oct 2014 18:25:36 +0000 (19:25 +0100)
committermivanov <marian.ivanov@cern.ch>
Tue, 28 Oct 2014 21:32:50 +0000 (22:32 +0100)
- adding user options: string melting and no ART
- adding Boost to LHC
- adding decayer

ANALYSIS/macros/train/AddMCGenAmpt.C

index 5db128d..e5d7efd 100644 (file)
@@ -1,37 +1,69 @@
-AliGenerator *AddMCGenAmpt()
+AliGenerator *AddMCGenAmpt(
+                          Double_t Energy      = 2760.,   // CM energy 
+                          Double_t bmin        = 0.0,     // minimum impact parameter
+                          Double_t bmax        = 20.0,    // maximum impact parameter
+                          Bool_t stringMelting = kTRUE,   // string melting option 
+                          Bool_t useART        = kTRUE,   // use hadronic rescattering phase (ART)
+                          )
 {
-// User defined generator
+  // User defined generator
 
   gSystem->Load("libampt.so");       
   gSystem->Load("libTAmpt.so");
 
+
   AliGenAmpt *genAMPT = new AliGenAmpt(-1);
+  //=========================================================================
+
+
+  // User settings
+  Int_t Flag_SM    = 4;       // flag for string melting: 1 = default, 4 = String Melting
+  Int_t NTmax      = 150;     // NTMAX: number of timesteps (Default = 150), to turn off ART set it to 3
+  Double_t mu      = 3.2264;  // parton screening mass in fm^(-1) (Default = 3.2264)
+  Double_t alpha_s = 1./3.;   // change mu and alpha_s (Default = 1./3.) to vary scattering cross-section
+                              // mu = 3.2 fm^-1 and alpha_s = 0.33 ==> sigma_{partonic} = 1.5mb
+  if(!stringMelting)
+    Flag_SM = 1;
+  
+  if(!useART)
+    NTmax = 3;
+  //=========================================================================
+
+
+  // Decayer
+  AliDecayer *decayer = new AliDecayerPythia();
+  genAMPT->SetForceDecay( kHadronicD );
+  genAMPT->SetDecayer( decayer );
+  //=========================================================================
 
-  // will be made optional later
-  genAMPT->SetEnergyCMS(2760);
+  // Collision system
+  genAMPT->SetEnergyCMS(Energy);
   genAMPT->SetReferenceFrame("CMS");
   genAMPT->SetProjectile("A", 208, 82);
   genAMPT->SetTarget    ("A", 208, 82);
   genAMPT->SetPtHardMin (2);
   genAMPT->SetImpactParameterRange(0.00,20.00);
-  genAMPT->SetJetQuenching(0); // enable jet quenching
-  genAMPT->SetShadowing(1);    // enable shadowing
-  genAMPT->SetDecaysOff(1);    // neutral pion and heavy particle decays switched off
-  genAMPT->SetSpectators(0);   // track spectators 
-  genAMPT->SetIsoft(4);        // 4=string melting, 1=standard AMPT
-  genAMPT->SetXmu(3.2264);     // parton xsection
-  genAMPT->SetNtMax(150);      // time bins
-  
-  genAMPT->SetAlpha(1./3.);    //alpha =0.333
-  genAMPT->SetStringFrag(0.5,0.9); //string fragmentation parameters
-  genAMPT->SetIpop(1); //enable popcorn mechanism (net-baryon stopping)
-  // This particular choice of gives scattering cross section to be 1.5 mb
+  //=========================================================================
+
+  // options
+  genAMPT->SetJetQuenching(0);     // enable jet quenching
+  genAMPT->SetShadowing(1);        // enable shadowing
+  genAMPT->SetDecaysOff(1);        // neutral pion and heavy particle decays switched off
+  genAMPT->SetSpectators(0);       // track spectators 
+  genAMPT->SetIsoft(Flag_SM);      // 4=string melting, 1=standard AMPT
+  genAMPT->SetXmu(mu);             // parton xsection
+  genAMPT->SetNtMax(NTmax);        // time bins
+  genAMPT->SetAlpha(alpha_s);      // alpha =0.333
+  genAMPT->SetStringFrag(0.5,0.9); // string fragmentation parameters
+  genAMPT->SetIpop(1);             // enable popcorn mechanism (net-baryon stopping)
+  //=========================================================================
 
+  // Boost into LHC lab frame
+  genAMPT->SetBoostLHC(1);
+  
+  // randomize reaction plane
   genAMPT->SetRandomReactionPlane(kTRUE);
 
  return genAMPT;
 
-
-
-
 }