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*' '/
WRITE(MSTU(11),5100)
STOP
ENDIF
-
C...Initial values for some counters.
N=0
MINT(5)=MINT(5)+1
&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 for the Lambda(1520).
C -------------------------------------------------------------------------
- IF (IABS(KF).EQ.3122) THEN
+C IF (IABS(KF).EQ.3122) THEN
C Converting a fraction (0.20) of Lambda0 to Lambda(1520) + c.c.
C This fraction is based on the experimental measurement at ISR
C Bobbink 83, NP B217,11 (1983)
C The region 0.5 < XF < 1.0 has been extrapolated to XF=0
- IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
- ENDIF
+C IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
+C ENDIF
- IF(IABS(KF).EQ.3212) THEN
+C IF(IABS(KF).EQ.3212) THEN
C Converting a fraction (0.20) of Sigma0 to Lambda(1520) + c.c.
C We suppose the same fraction as for Lambda0
- IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
- ENDIF
+C IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
+C ENDIF
- IF (IABS(KF).EQ.3214) THEN
+C IF (IABS(KF).EQ.3214) THEN
C Converting a fraction (0.30) of Sigma0(1385) to Lambda(1520) + c.c.
C This is conservative extimate supposing that the ratio
C scales as (M_Sigma1385/M_Lambda0)^2 ~ 1.5
- IF(PYR(0).LE.0.30) KF=ISIGN(3124,KF)
- ENDIF
+C IF(PYR(0).LE.0.30) KF=ISIGN(3124,KF)
+C ENDIF
RETURN
C...Use tabulated probabilities to select new flavour and hadron.
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.