geant_gep/src/gukine.f - diff

Return to gukine.f CVS log

Up to [HallC] / geant_gep / src

Diff for /geant_gep/src/gukine.f between version 1.5 and 1.8

version 1.5, 2005/09/22 17:17:49

version 1.8, 2005/10/11 22:54:47

Line 132

> x_fp, dx_fp, y_fp, dy_fp, m2,

> ms_flag, wcs_flag, decay_flag, resmult, fry,

> ok_spec, pathlen)

c write(*,*)p_spec,th_spec,dpp

c write(*,*)'Momentum = ',p_spec,' Angle = ',th_spec

c write(*,*)x,y,z

c write(*,*)dxdz,dydz

write(6,*)x_fp,dx_fp,y_fp,dy_fp,ok_spec

c write(6,*)x_fp,dx_fp,y_fp,dy_fp,ok_spec

c write(*,*)m2,ms_flag,wcs_flag,decay_flag

c write(*,*)resmult,fry,ok_spec,pathlen

Line 143

phfp=atan(dx_fp)/3.14159265*180.0

yfp=y_fp*100.0

thfp=atan(dy_fp)/3.14159265*180.0*100.0

c write(6,*)'x,phi,y,theta,ok_spec ='

c write(6,*)xfp,phfp,yfp,thfp,ok_spec

do ii=1,20

ntuple_array(ii)=0.0

Line 321

Line 323

c what the corresponding theta and phi are for the other arm.

if (abs(eang).gt.abs(hang)) then

phie=-.025+rndm(1)*.050

c phie=-.025+rndm(1)*.050

thetae=-.009+rndm(2)*.018

c thetae=-.009+rndm(2)*.018

c phie=-.0+rndm(1)*.0

if(pcentral.ge.5000) then

c thetae=-.0+rndm(2)*.0

phie=-.130+rndm(1)*.260

thetae=-.065+rndm(2)*.130

else if(pcentral.ge.3000.and.pcentral.lt.5000) then

phie=-.067+rndm(1)*.135

thetae=-.034+rndm(2)*.067

else

phie=-.087+rndm(1)*.174

thetae=-.044+rndm(2)*.087

endif

c phie=-.00+rndm(1)*.00

c thetae=-.00+rndm(2)*.00

escat=mp/(1.0+mp/e0-cos(fg*eang+thetae)*cos(phie))

pscat=sqrt(e0**2+escat**2-

$ 2.0*e0*escat*cos(fg*eang+thetae)*cos(phie))

Line 334

Line 346

else

c Hadron arm defining acceptance

1221 call grndm ( rndm, 3)

phie=-.025+rndm(1)*.050

phie=-.080+rndm(1)*.160