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*' '/
&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)
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)
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)
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.
RETURN
END
+
+ SUBROUTINE PYPTFS(MODE,PTMAX,PTMIN,PTGEN)
+C DUMMY
+C common /pyqpar/ T0u,tau0u,nfu,ienglu,ianglu
+C common /parimp/ b1,psib1,rb1,rb2,noquen
+ RETURN
+ END