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);
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.
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.