@@ -500,6 +500,8 @@ SUBROUTINE LUPREP(IP)
500500
501501C...Rearrange parton shower product listing along strings: begin loop.
502502 I1=N
503+ DO 101 J=1,5
504+ 101 DPC(J)=0
503505 DO 130 MQGST=1,2
504506 DO 120 I=MAX(1,IP),N
505507 IF(K(I,1).NE.3) GOTO 120
@@ -3459,6 +3461,11 @@ SUBROUTINE LUSHOW(IP1,IP2,QMAX)
34593461 PT2MIN=MAX(0.5*PARJ(82),1.1*PARJ(81))**2
34603462 ALAMS=PARJ(81)**2
34613463 ALFM=LOG(PT2MIN/ALAMS)
3464+ DO 101 I=1,4
3465+ ISI(I)=0
3466+ IPA(I)=0
3467+ KFLD(I)=0
3468+ 101 CONTINUE
34623469
34633470C...Store positions of shower initiating partons.
34643471 M3JC=0
@@ -9582,6 +9589,8 @@ SUBROUTINE PYMULT(MMUL)
95829589 XI=0.
95839590 YI=0.
95849591 XK=0.5
9592+ XF=1.
9593+ YF=1.
95859594 IIT=0
95869595 130 IF(IIT.EQ.0) THEN
95879596 XK=2.*XK
@@ -9877,6 +9886,9 @@ SUBROUTINE PYMULT(MMUL)
98779886
98789887 IF(RFLAV.LT.PARP(85).AND.NSTR.GE.1) THEN
98799888C....Choose relevant string pieces to place gluons on.
9889+ IST1=0
9890+ IST2=0
9891+ ISTM=0
98809892 DO 210 I=N+1,N+2
98819893 DMIN=1E8
98829894 DO 200 ISTR=1,NSTR
@@ -11248,6 +11260,7 @@ SUBROUTINE PYWIDT(KFLR,RMAS,WDTP,WDTE)
1124811260C...H0 -> gamma + gamma; quark, charged lepton and W loop contributions
1124911261 ETARE=0.
1125011262 ETAIM=0.
11263+ EJ=0.
1125111264 DO 150 J=1,3*MSTP(1)+1
1125211265 IF(J.LE.2*MSTP(1)) THEN
1125311266 EJ=KCHG(J,1)/3.
@@ -11289,7 +11302,9 @@ SUBROUTINE PYWIDT(KFLR,RMAS,WDTP,WDTE)
1128911302 ELSEIF(I.EQ.15) THEN
1129011303C...H0 -> gamma + Z0; quark, charged lepton and W loop contributions
1129111304 ETARE=0.
11292- ETAIM=0.
11305+ ETAIM=0.
11306+ VJ=0.
11307+ EJ=0.
1129311308 DO 160 J=1,3*MSTP(1)+1
1129411309 IF(J.LE.2*MSTP(1)) THEN
1129511310 EJ=KCHG(J,1)/3.
@@ -11304,7 +11319,7 @@ SUBROUTINE PYWIDT(KFLR,RMAS,WDTP,WDTE)
1130411319 VJ=AJ-4.*EJ*XW
1130511320 EPS=(2.*PMAS(10+JL,1)/RMAS)**2
1130611321 EPSP=(2.*PMAS(10+JL,1)/PMAS(23,1))**2
11307- ELSE
11322+ ELSE
1130811323 EPS=(2.*PMAS(24,1)/RMAS)**2
1130911324 EPSP=(2.*PMAS(24,1)/PMAS(23,1))**2
1131011325 ENDIF
0 commit comments