fQuench(0),
fQhat(0.),
fLength(0.),
+ fpyquenT(1.),
+ fpyquenTau(0.1),
+ fpyquenNf(0),
+ fpyquenEloss(0),
+ fpyquenAngle(0),
fImpact(0.),
fPtKick(1.),
fFullEvent(kTRUE),
fQuench(kFALSE),
fQhat(0.),
fLength(0.),
+ fpyquenT(1.),
+ fpyquenTau(0.1),
+ fpyquenNf(0),
+ fpyquenEloss(0),
+ fpyquenAngle(0),
fImpact(0.),
fPtKick(1.),
fFullEvent(kTRUE),
fPythia->InitQuenching(0., 0.1, 0.6e6, 0);
}
+ if(fQuench ==2){fPythia->SetPyquenParameters(fpyquenT,fpyquenTau,fpyquenNf,fpyquenEloss,fpyquenAngle);}
+
if (fQuench == 3) {
// Nestor's change of the splittings
fPythia->SetPARJ(200, 0.8);
void SetQhat(Float_t qhat) {fQhat = qhat;}
//Set initial medium length.
void SetLength(Float_t length) {fLength = length;}
-
+ //set parameters for pyquen afterburner
+ virtual void SetPyquenPar(Float_t t0=1., Float_t tau0=0.1, Int_t nf=0,Int_t iengl=0, Int_t iangl=3)
+ {fpyquenT = t0; fpyquenTau = tau0; fpyquenNf=nf;fpyquenEloss=iengl;fpyquenAngle=iangl;}
virtual void SetHadronisation(Int_t flag = 1) {fHadronisation = flag;}
virtual void SetPatchOmegaDalitz(Int_t flag = 1) {fPatchOmegaDalitz = flag;}
virtual void SetReadFromFile(const Text_t *filname) {fkFileName = filname; fReadFromFile = 1;}
Int_t fQuench; //Flag for quenching
Float_t fQhat; //Transport coefficient (GeV^2/fm)
Float_t fLength; //Medium length (fm)
+ Float_t fpyquenT; //Pyquen initial temperature
+ Float_t fpyquenTau; //Pyquen initial proper time
+ Int_t fpyquenNf; //Pyquen number of flavours into the game
+ Int_t fpyquenEloss; //Pyquen type of energy loss
+ Int_t fpyquenAngle; //Pyquen radiation angle for gluons
Float_t fImpact; //Impact parameter for quenching simulation (q-pythia)
Float_t fPtKick; //Transverse momentum kick
Bool_t fFullEvent; //!Write Full event if true