|
version 1.4, 2005/09/22 12:34:24
|
version 1.5, 2005/09/22 17:17:49
|
|
|
|
| | |
| xfp=x_fp | xfp=x_fp |
| phfp=atan(dx_fp)/3.14159265*180.0 | phfp=atan(dx_fp)/3.14159265*180.0 |
| yfp=y_fp |
yfp=y_fp*100.0 |
| thfp=atan(dy_fp)/3.14159265*180.0 |
thfp=atan(dy_fp)/3.14159265*180.0*100.0 |
| c | c |
| do ii=1,20 | do ii=1,20 |
| ntuple_array(ii)=0.0 | ntuple_array(ii)=0.0 |
|
|
|
| c what the corresponding theta and phi are for the other arm. | c what the corresponding theta and phi are for the other arm. |
| | |
| if (abs(eang).gt.abs(hang)) then | if (abs(eang).gt.abs(hang)) then |
| phie=-.065+rndm(1)*.130 |
phie=-.025+rndm(1)*.050 |
| thetae=-.030+rndm(2)*.060 |
thetae=-.009+rndm(2)*.018 |
| |
c phie=-.0+rndm(1)*.0 |
| |
c thetae=-.0+rndm(2)*.0 |
| escat=mp/(1.0+mp/e0-cos(fg*eang+thetae)*cos(phie)) | escat=mp/(1.0+mp/e0-cos(fg*eang+thetae)*cos(phie)) |
| pscat=sqrt(e0**2+escat**2- | pscat=sqrt(e0**2+escat**2- |
| $ 2.0*e0*escat*cos(fg*eang+thetae)*cos(phie)) | $ 2.0*e0*escat*cos(fg*eang+thetae)*cos(phie)) |
|
|
|
| else | else |
| c Hadron arm defining acceptance | c Hadron arm defining acceptance |
| 1221 call grndm ( rndm, 3) | 1221 call grndm ( rndm, 3) |
| phie=-.065+rndm(1)*.130 |
phie=-.025+rndm(1)*.050 |
| thetae=-.030+rndm(2)*.060 |
thetae=-.009+rndm(2)*.018 |
| |
c phie=-.0+rndm(1)*.0 |
| |
c thetae=-.0+rndm(2)*.0 |
| escat=mp/(1.0+mp/e0-cos(fg*eang+thetae)*cos(phie)) | escat=mp/(1.0+mp/e0-cos(fg*eang+thetae)*cos(phie)) |
| pscat=sqrt(e0**2+escat**2- | pscat=sqrt(e0**2+escat**2- |
| $ 2.0*e0*escat*cos(fg*eang+thetae)*cos(phie)) | $ 2.0*e0*escat*cos(fg*eang+thetae)*cos(phie)) |
|
|
|
| endif | endif |
| | |
| ptgt=pscat | ptgt=pscat |
| |
c write(*,*)escat,eang,pscat,hang |
| dptgt=(pcentral-pscat)/pcentral | dptgt=(pcentral-pscat)/pcentral |
| c | c |
| c Following statement to fill just the high and low dp bins | c Following statement to fill just the high and low dp bins |
Return to
gukine.f
CVS log
Up to [HallC] / geant_gep / src