add pyquen quenching strength

diff --git a/ANALYSIS/macros/train/AddMCGenQuench.C b/ANALYSIS/macros/train/AddMCGenQuench.C
index 4c39a18530add1cedc06453721df14b9efee4178..c1eef73fd0690368580d2b101e76552b2aab88b0 100644 (file)
--- a/ANALYSIS/macros/train/AddMCGenQuench.C
+++ b/ANALYSIS/macros/train/AddMCGenQuench.C
@@ -23,16 +23,16 @@
 //If you want k in fm then you have to divide by the squared of hbarc (hbar c=0.197 GeV fm)
 // This gives k=4.4e^6 fm, which is the quench value we set in SetQhat
 
-AliGenerator* AddMCGenQuench(Float_t e_cms = 2760., Double_t ptHardMin = 0., Double_t ptHardMax = 0., Int_t fpythia = 1, Double_t quench=4.4e6) 
+AliGenerator* AddMCGenQuench(Float_t e_cms = 2760., Double_t ptHardMin = 0., Double_t ptHardMax = 0., Int_t fpythia = 1, Double_t quench=4.4e6, ,Int_t ianglepyquen = 2) 
 {
   //Add Pythia generator: pt-hard bin or min bias
 
   gSystem->Load("liblhapdf.so");
  
-  return CreatePythia6Gen(e_cms, ptHardMin, ptHardMax, fpythia, quench);
+  return CreatePythia6Gen(e_cms, ptHardMin, ptHardMax, fpythia, quench, ianglepyquen);
 }
 
-AliGenerator* CreatePythia6Gen(Float_t e_cms, Int_t ptHardMin, Int_t ptHardMax, Int_t fpythia, Double_t quench) {
+AliGenerator* CreatePythia6Gen(Float_t e_cms, Int_t ptHardMin, Int_t ptHardMax, Int_t fpythia, Double_t quench = 4.4e6, Int_t ianglepyquen=2) {
     
   gSystem->Load("libqpythia.so");
   gSystem->Load("libEGPythia6.so");
@@ -64,6 +64,7 @@ AliGenerator* CreatePythia6Gen(Float_t e_cms, Int_t ptHardMin, Int_t ptHardMax,
   if (fpythia == 2){
     genP->SetTune(103); //tune DW, standard choice for Q2 showers
     genP->SetQuench(2);
+    genP->SetPyquenPar(1,0.1,0,0,ianglepyquen);
   }
 
   return genP;