Fix some coding convention violations.

[u/mrichter/AliRoot.git] / PYTHIA6 / pyquen1_5.F

diff --git a/PYTHIA6/pyquen1_5.F b/PYTHIA6/pyquen1_5.F
index f99514005045265551c5c1a3a0ce4aeccea96432..06aea339e6131df14556b8409a02d7ae58971fc5 100644 (file)
--- a/PYTHIA6/pyquen1_5.F
+++ b/PYTHIA6/pyquen1_5.F
@@ -896,6 +896,7 @@
        bub=4.d0*tauu/TC   
        alf=6.d0*pi/((33.d0-2.d0*nf)*dlog(max(bub,1.d-10)))
        z=pi*4.d0*tauu*tauu*alf*(1.+nf/6.d0)  
+       if (z < 0.) write(6,*) "Warning in PLJETR: z < 0.", z
        bubs=dsqrt(abs(z))/TC 
        alfs=6.d0*pi/((33.d0-2.d0*nf)*dlog(max(bubs,1.d-10))) 
        phmin2=z 
@@ -1296,13 +1297,22 @@ C************************** END PARINV *************************************
 * if  quark production outside the QGP then 
 * arg=(((dsin(uu)*cosh(uu))**2)+((dcos(uu)*sinh(uu))**2))/(2.d0*uu*uu);   
 * here quark production inside the QGP  
-      arg=((dcos(uu)*cosh(uu))**2)+((dsin(uu)*sinh(uu))**2)  
-      write(6,*) "Warning in PLFUN1: z < 0.", z
+c AM: better numerical stability      
+      
+      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
+      endif
+         
+
       gl1=(ygl/(cfac*abs(z)))**0.3333333d0
       gl2=(um/epa)**1.333333d0  
-      dc=1.d0/((1.d0+((gl1*gl2*or)**1.5d0))**2)       ! massive parton    
+      dc=1.d0/((1.d0+((gl1*gl2*or)**1.5d0))**2)        ! massive parton    
 c      dc=1.d0                                         !massless parton 
-      plfun1=dc*3.d0*alfs*ygl*dlog(max(arg,1.d-20))*spinf/(pi*al*or)   
+c      plfun1 = dc*3.d0*alfs*ygl*dlog(max(arg,1.d-20))*spinf/(pi*al*or)   
+      plfun1   = dc*3.d0*alfs*ygl*xlogarg*spinf/(pi*al*or) 
       return 
       end