Another fix avoiding fpe
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Apr 2011 15:58:16 +0000 (15:58 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Apr 2011 15:58:16 +0000 (15:58 +0000)
PYTHIA6/pythia6.4.25/pythia-6.4.25.f

index b888a24..be6c5d1 100644 (file)
@@ -15309,8 +15309,12 @@ C...Evolve q -> q gamma separately, pick it if larger pT.
 C...  Evolve massive quark creation separately.
         MCRQQ=0
         IF (MQMASS.NE.0) THEN
-          PT2CR=(RMQ2+VINT(18))*(RML**(TPM/(TPL*PYR(0)**(-TML/WN)-TPM)))
-     &         -VINT(18)
+           IF (WN .EQ. 0.) THEN
+              ARG = -1.
+           ELSE
+              ARG = TPM/(TPL*PYR(0)**(-TML/WN)-TPM)
+           ENDIF
+           PT2CR=(RMQ2+VINT(18))*(RML**ARG)-VINT(18)
 C...If massive quark also on opposite side, ensure sufficient remaining 
 C...phase space also for creation of that quark
           TMINQQ = TMIN