version 1.7, 2005/10/05 16:07:27
|
version 1.9, 2005/12/16 07:23:56
|
|
|
z=0.0 | z=0.0 |
dxdz=tan(phtgt*3.14159265/180.0) | dxdz=tan(phtgt*3.14159265/180.0) |
dydz=tan(thtgt*3.14159265/180.0) | dydz=tan(thtgt*3.14159265/180.0) |
c write(*,*)p_spec,th_spec,dpp |
if(nevent.le.2)write(*,*)'******** mom = ',p_spec,' ******' |
c write(*,*)x,y,z | c write(*,*)x,y,z |
c write(*,*)dxdz,dydz | c write(*,*)dxdz,dydz |
c | c |
|
|
> x_fp, dx_fp, y_fp, dy_fp, m2, | > x_fp, dx_fp, y_fp, dy_fp, m2, |
> ms_flag, wcs_flag, decay_flag, resmult, fry, | > ms_flag, wcs_flag, decay_flag, resmult, fry, |
> ok_spec, pathlen) | > ok_spec, pathlen) |
write(*,*)'Momentum = ',p_spec,' Angle = ',th_spec |
c write(*,*)'Momentum = ',p_spec,' Angle = ',th_spec |
c write(*,*)x,y,z | c write(*,*)x,y,z |
c write(*,*)dxdz,dydz | c write(*,*)dxdz,dydz |
c 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 |
|
|
phfp=atan(dx_fp)/3.14159265*180.0 | phfp=atan(dx_fp)/3.14159265*180.0 |
yfp=y_fp*100.0 | yfp=y_fp*100.0 |
thfp=atan(dy_fp)/3.14159265*180.0*100.0 | thfp=atan(dy_fp)/3.14159265*180.0*100.0 |
write(6,*)'x,phi,y,theta,ok_spec =' |
c write(6,*)'x,phi,y,theta,ok_spec =' |
write(6,*)xfp,phfp,yfp,thfp,ok_spec |
c write(6,*)xfp,phfp,yfp,thfp,ok_spec |
c | c |
do ii=1,20 | do ii=1,20 |
ntuple_array(ii)=0.0 | ntuple_array(ii)=0.0 |
|
|
c phie=-.025+rndm(1)*.050 | c phie=-.025+rndm(1)*.050 |
c thetae=-.009+rndm(2)*.018 | c thetae=-.009+rndm(2)*.018 |
if(pcentral.ge.5000) then | if(pcentral.ge.5000) then |
|
c phie=-.00+rndm(1)*.00 |
|
c thetae=-.00+rndm(2)*.00 |
phie=-.130+rndm(1)*.260 | phie=-.130+rndm(1)*.260 |
thetae=-.065+rndm(2)*.130 | thetae=-.065+rndm(2)*.130 |
else if(pcentral.ge.3000.and.pcentral.lt.5000) then | else if(pcentral.ge.3000.and.pcentral.lt.5000) then |