X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=TPHIC%2Ftphic17.f;h=e00da09aebd81b68d62494bf186658e230d773d4;hp=8ffe99a36bb4922364c1f32b8b73e89ea91abebd;hb=8bd87c5e04750bb79676819d9cd5318961713160;hpb=0b5dd07122f86c333da3b33467d629f10db57c6a

diff --git a/TPHIC/tphic17.f b/TPHIC/tphic17.f
index 8ffe99a36bb..e00da09aebd 100644
--- a/TPHIC/tphic17.f
+++ b/TPHIC/tphic17.f
@@ -97,7 +97,7 @@ C          THETAB               = sbottom mixing angle
       REAL AMZISS(4)
       EQUIVALENCE (AMZISS(1),AMZ1SS)
       SAVE /SSPAR/
-      REAL pgam1(4),pgam2(4),xpar(10),ypar(6)
+      REAL xpar(10),ypar(6)
       INTEGER ipar(2)
       INTEGER pycomp
       REAL*8 ggrnd
@@ -445,7 +445,7 @@ CDECK  ID>, GGRUN.
       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
       DOUBLE PRECISION PMAS,PARF,VCKM
       SAVE  /PYDAT2/
-      REAL pgam1(4),pgam2(4),pev(4,10)
+      REAL pgam1(4),pgam2(4)
       REAL p3(4), p4(4)
       REAL*8 dbetaz
       REAL*8 ggrnd
@@ -526,7 +526,7 @@ CDECK  ID>, GGRUN.
           e3     = xkin(5)
           e4     = xmgg - e3
         END IF
-        sinthe = sqrt (1 - costhe*costhe)
+        sinthe = sqrt ((1.-costhe)*(1.+costhe))
         p3(1) = p3abs * cos(phi) * sinthe
         p3(2) = p3abs * sin(phi) * sinthe
         p3(3) = p3abs * costhe
@@ -598,19 +598,19 @@ CDECK  ID>, GGEXIT.
      &               gvconst(4,10)
       REAL nc
 *
-      WRITE (6,'(/x,79(''=''))')
-      WRITE (6,'(x,''I'',36x,''TPHIC'',36x,      ''I'')')
-      WRITE (6,'(x,''I'',34x,''Process '',i1,34x,''I'')') iproc
-      WRITE (6,'(x,''I'',34x,''---------'',  34x,''I'')')
-      WRITE (6,'(x,''I'',77x,                    ''I'')')
-      WRITE (6,'(x,''I'',20x,''Events thrown   : '',i8,31x,''I'')')
+      WRITE (6,'(/1x,79(''=''))')
+      WRITE (6,'(1x,''I'',36x,''TPHIC'',36x,      ''I'')')
+      WRITE (6,'(1x,''I'',34x,''Process '',i1,34x,''I'')') iproc
+      WRITE (6,'(1x,''I'',34x,''---------'',  34x,''I'')')
+      WRITE (6,'(1x,''I'',77x,                    ''I'')')
+      WRITE (6,'(1x,''I'',20x,''Events thrown   : '',i8,31x,''I'')')
      &       ntry
-      WRITE (6,'(x,''I'',20x,''Events accepted : '',i8,31x,''I'')')
+      WRITE (6,'(1x,''I'',20x,''Events accepted : '',i8,31x,''I'')')
      &       nevent
-      WRITE (6,'(x,''I'',20x,''Cross section   : '',e10.3,'' +- '','//
+      WRITE (6,'(1x,''I'',20x,''Cross section   : '',e10.3,'' +- '','//
      &      'e10.3,'' nb'',12x ''I'')') xstot, xstote
-      WRITE (6,'(x,''I'',77x,                    ''I'')')
-      WRITE (6,'(x,79(''=''))')
+      WRITE (6,'(1x,''I'',77x,                    ''I'')')
+      WRITE (6,'(1x,79(''=''))')
 *
       IF (iproc .EQ. 1) THEN
 *   ---  Alternative cross section for minimum bias events ---
@@ -826,10 +826,10 @@ CDECK  ID>, GGDECY.
 *         Decay chi_1^+- --> chi_1^0 W^+- in CMS, uniformly
         eb = (amprt(1)**2 + amprt(2)**2 - amprt(4)**2) / (2 * amprt(1))
         ec = (amprt(1)**2 - amprt(2)**2 + amprt(4)**2) / (2 * amprt(1))
-        pb = sqrt (eb*eb - amprt(2)*amprt(2))
+        pb = sqrt ((eb-amprt(2))*(eb+amprt(2)))
         pc = pb
         costhe = 2*ggrnd(0) - 1.
-        sinthe = sqrt (1 - costhe*costhe)
+        sinthe = sqrt((1.-costhe)*(1.+costhe))
         phi = 2*pi*ggrnd(0)
         pcm(1) = pb * sinthe * cos(phi)
         pcm(2) = pb * sinthe * sin(phi)
@@ -878,7 +878,7 @@ CDECK  ID>, GGDECY.
 *
         phi    = 2.*pi * ggrnd(0)
         costhe = 2.*ggrnd(0) - 1.
-        sinthe = sqrt (1. - costhe*costhe)
+        sinthe = sqrt ((1.-costhe)*(1.+costhe))
 *
         pp = xm23 / 2.
         p2(1) = pp * sinthe * cos(phi)
@@ -895,7 +895,7 @@ CDECK  ID>, GGDECY.
 *       3-momenta p2, p3 are along z-axis
 *
         gamma = e23 / xm23
-        gambe = sqrt (gamma*gamma - 1.)
+        gambe = sqrt ((gamma-1.)*(gamma+1.))
 *
         p2z = gamma*p2(3) + gambe*p2(4)
         e2  = gamma*p2(4) + gambe*p2(3)
@@ -911,7 +911,7 @@ CDECK  ID>, GGDECY.
 *
         phi    = 2.*pi * ggrnd(0)
         costhe = 2.*ggrnd(0) - 1.
-        sinthe = sqrt (1. - costhe*costhe)
+        sinthe = sqrt ((1.-costhe)*(1.+costhe))
         cosphi = cos(phi)
         sinphi = sin(phi)
 *
@@ -976,7 +976,7 @@ CDECK  ID>, GGDECY.
         p12 = am1 * p2(4)
 *
         e4 = am1 - e23
-        pp =sqrt (e4*e4 - am4*am4)
+        pp =sqrt ((e4-am4)*(e4+am4))
         p4(1) =  pp * sinthe * cos(phi)
         p4(2) =  pp * sinthe * sin(phi)
         p4(3) = -pp * costhe
@@ -1160,17 +1160,16 @@ C***********************************************************************
       REAL      xmres,xgres
       LOGICAL   lres,lwate
 C
-      REAL       alpha,pi,wlum,weight
+      REAL       alpha,pi
       PARAMETER (alpha=1./137.036,pi=3.14159265)
 C
       REAL ptag1(10),ptag2(10),pgam1(10),pgam2(10),wt
-      REAL amas,wsq,q1sq,q2sq
+      REAL amas,wsq
 C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
-      REAL*4  xmin,xmax,ymin,ymax,gmma,y3,xm,q0,qs,qs1,qs2
-      REAL*8  pa(5),pb(5),p1(5),p2(5),p3(5)
+      REAL*4  y3,xm,q0,qs,qs1,qs2
+      REAL*8  pa(5),pb(5)
 C
-      REAL*8 ph1,ph2,ak1,ak2,qk1,qk2,q1,q2,q3,e1,e2,e3
-      REAL*8 y,y1,y2,phi,pi2,pae,aeq,beq,wps,tm3,dts,wpt
+      REAL*8  pi2
 C
       LOGICAL     lstr
       DATA q0,pi2,lstr /  0.060, 6.28318530718 d0, .false. /
@@ -1185,7 +1184,7 @@ C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 C
       pa(1)= 0D0
       pa(2)= 0D0
-      pa(3)= dsqrt(1D0*(eb*eb-amas*amas))
+      pa(3)= dsqrt(1D0*((eb-amas)*(eb+amas)))
       pa(4)= eb*1D0
       pa(5)= amas*1D0
 C