Correction.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 8 Nov 2008 17:05:36 +0000 (17:05 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 8 Nov 2008 17:05:36 +0000 (17:05 +0000)
PYTHIA6/pyquen1_5.F

index e3846aa..06aea33 100644 (file)
@@ -1298,8 +1298,9 @@ C************************** END PARINV *************************************
 * arg=(((dsin(uu)*cosh(uu))**2)+((dcos(uu)*sinh(uu))**2))/(2.d0*uu*uu);   
 * here quark production inside the QGP  
 c AM: better numerical stability      
-      arg     = ((dcos(uu)*cosh(uu))**2)+((dsin(uu)*sinh(uu))**2)  
+      
       if (uu .lt. 5.) then
+         arg     = ((dcos(uu)*cosh(uu))**2)+((dsin(uu)*sinh(uu))**2)  
          xlogarg = dlog(max(arg,1.d-20))
       else
          xlogarg = log(0.25) + 2. * uu