Modified splittings for fQuench = 3.

diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx
index 783b434..d48b00c 100644 (file)
--- a/PYTHIA6/AliGenPythia.cxx
+++ b/PYTHIA6/AliGenPythia.cxx
@@ -520,6 +520,8 @@ void AliGenPythia::Init()
     if (fQuench) {
        fPythia->InitQuenching(0., 0.1, 0.6e6, 0);
     }
+    fPythia->SetPARJ(200, 0.0);
+
     if (fQuench == 3) {
        // Nestor's change of the splittings
        fPythia->SetPARJ(200, 0.8);

diff --git a/PYTHIA6/pythia6214.f b/PYTHIA6/pythia6214.f
index 922c174..22f08da 100644 (file)
--- a/PYTHIA6/pythia6214.f
+++ b/PYTHIA6/pythia6214.f
@@ -51415,10 +51415,13 @@ C...Calculate allowed z range.
  
 C...Integral of Altarelli-Parisi z kernel for QCD.
 C...(Includes squark and gluino; with factor N_C/C_F extra for latter).
+      FMED = PARJ(200)
       IF(MSTJ(49).EQ.0.AND.KFL(1).EQ.21) THEN
-        FBR=6D0*LOG((1D0-ZC)/ZC)+MSTJ(45)*0.5D0
+C Nestor
+        FBR=(1.D0+FMED)*6D0*LOG((1D0-ZC)/ZC)+MSTJ(45)*0.5D0
       ELSEIF(MSTJ(49).EQ.0) THEN
-        FBR=(8D0/3D0)*LOG((1D0-ZC)/ZC)
+C Nestor
+        FBR=(1.D0+FMED)*(8D0/3D0)*LOG((1D0-ZC)/ZC)
         IF(IGLUI.EQ.1.AND.IR.GE.31) FBR=FBR*(9D0/4D0)
  
 C...Integral of Altarelli-Parisi z kernel for scalar gluon.
@@ -51430,9 +51433,9 @@ C...Integral of Altarelli-Parisi z kernel for scalar gluon.
  
 C...Integral of Altarelli-Parisi z kernel for Abelian vector gluon.
       ELSEIF(KFL(1).EQ.21) THEN
-        FBR=6D0*MSTJ(45)*(0.5D0-ZC)
+        FBR=(1.D0+FMED)*6D0*MSTJ(45)*(0.5D0-ZC)
       ELSE
-        FBR=2D0*LOG((1D0-ZC)/ZC)
+        FBR=(1.D0+FMED)*2D0*LOG((1D0-ZC)/ZC)
       ENDIF
  
 C...Reset QCD probability for colourless.