From e17f27a2b694b8b46bd22e31836e4e128b1d0b17 Mon Sep 17 00:00:00 2001 From: mverweij Date: Sun, 22 Jun 2014 20:14:24 +0200 Subject: [PATCH 1/1] add pyquen quenching strength --- ANALYSIS/macros/train/AddMCGenQuench.C | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/ANALYSIS/macros/train/AddMCGenQuench.C b/ANALYSIS/macros/train/AddMCGenQuench.C index 4c39a18530a..c1eef73fd06 100644 --- 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; -- 2.39.3