Pyquen large angle radiation
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 3 Dec 2011 08:48:56 +0000 (08:48 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 3 Dec 2011 08:48:56 +0000 (08:48 +0000)
Leticia Cunqueiro

PYTHIA6/AliGenPythia.cxx
PYTHIA6/AliGenPythia.h

index 62d334a..b6cc0ef 100644 (file)
@@ -84,6 +84,11 @@ AliGenPythia::AliGenPythia():
     fQuench(0),
     fQhat(0.),
     fLength(0.),
+    fpyquenT(1.),
+    fpyquenTau(0.1),
+    fpyquenNf(0),
+    fpyquenEloss(0),
+    fpyquenAngle(0), 
     fImpact(0.),
     fPtKick(1.),
     fFullEvent(kTRUE),
@@ -182,6 +187,11 @@ AliGenPythia::AliGenPythia(Int_t npart)
      fQuench(kFALSE),
      fQhat(0.),
      fLength(0.),
+     fpyquenT(1.),
+     fpyquenTau(0.1),
+     fpyquenNf(0),
+     fpyquenEloss(0),
+     fpyquenAngle(0), 
      fImpact(0.),
      fPtKick(1.),
      fFullEvent(kTRUE),
@@ -530,6 +540,8 @@ void AliGenPythia::Init()
        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);
index 38e7e73..4b5a253 100644 (file)
@@ -145,7 +145,9 @@ class AliGenPythia : public AliGenMC
     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;}    
@@ -231,6 +233,11 @@ class AliGenPythia : public AliGenMC
     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