- adding refined bad channel map

[u/mrichter/AliRoot.git] / PYTHIA6 / pythia6214.f

diff --git a/PYTHIA6/pythia6214.f b/PYTHIA6/pythia6214.f
index 922c17410d74e237fbaa026110e307d37db946ae..6deccbf2e1a11fb229a1536ea9da8c0cd0386ac2 100644 (file)
--- a/PYTHIA6/pythia6214.f
+++ b/PYTHIA6/pythia6214.f
@@ -2563,8 +2563,8 @@ C...Local arrays and character variables.
       CHARACTER CHFRAM*12,CHBEAM*12,CHTARG*12,CHLH(2)*6
  
 C...Interface to PDFLIB.
-      COMMON/W50512/QCDL4,QCDL5
-      SAVE /W50512/
+      COMMON/LW50512/QCDL4,QCDL5
+      SAVE /LW50512/
       DOUBLE PRECISION VALUE(20),QCDL4,QCDL5
       CHARACTER*20 PARM(20)
       DATA VALUE/20*0D0/,PARM/20*' '/
@@ -27692,8 +27692,8 @@ C...Local arrays.
      &XPPI(-6:6),XPPR(-6:6)
  
 C...Interface to PDFLIB.
-      COMMON/W50513/XMIN,XMAX,Q2MIN,Q2MAX
-      SAVE /W50513/
+      COMMON/LW50513/XMIN,XMAX,Q2MIN,Q2MAX
+      SAVE /LW50513/
       DOUBLE PRECISION XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GLU,
      &VALUE(20),XMIN,XMAX,Q2MIN,Q2MAX
       CHARACTER*20 PARM(20)
@@ -28358,8 +28358,8 @@ C...Local arrays.
       DIMENSION XPEL(-25:25),XPGA(-6:6),SXP(0:6)
  
 C...Interface to PDFLIB.
-      COMMON/W50513/XMIN,XMAX,Q2MIN,Q2MAX
-      SAVE /W50513/
+      COMMON/LW50513/XMIN,XMAX,Q2MIN,Q2MAX
+      SAVE /LW50513/
       DOUBLE PRECISION XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GLU,
      &VALUE(20),XMIN,XMAX,Q2MIN,Q2MAX
       CHARACTER*20 PARM(20)
@@ -45412,7 +45412,6 @@ C...Commonblocks.
       SAVE /PYJETS/,/PYDAT1/,/PYDAT2/,/PYDAT3/,/PYINT4/
 C...Local array.
       DIMENSION PS(2,6),IJOIN(100)
- 
 C...Initialize and reset.
       MSTU(24)=0
       IF(MSTU(12).GE.1) CALL PYLIST(0)
@@ -51415,10 +51414,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 +51432,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.