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
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
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
& 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 ---
* 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)
*
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)
* 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)
*
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)
*
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
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. /
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