|
version 1.14, 2006/01/26 21:52:59
|
version 1.19, 2006/03/22 11:28:52
|
|
|
|
| integer ic,ii,jj | integer ic,ii,jj |
| integer*4 junk1,ikinsetting,isectrack | integer*4 junk1,ikinsetting,isectrack |
| real*8 xfp,yfp,tthfp,tphfp,pfp,junk2,junk3,thfp,phfp | real*8 xfp,yfp,tthfp,tphfp,pfp,junk2,junk3,thfp,phfp |
| real*8 e0(10),eang(10),hang(10),targ_thick(10) |
real*8 e0(12),eang(12),hang(12),targ_thick(12) |
| real*8 xtgt,ytgt,thtgt,phtgt,ptgt,dptgt | real*8 xtgt,ytgt,thtgt,phtgt,ptgt,dptgt |
| real yf,ff,thetaf,dptgt2 |
real yf,ff,thetaf |
| character*80 junkline | character*80 junkline |
| c | c |
| c The following variables are the SIMC variables | c The following variables are the SIMC variables |
|
|
|
| nhu4=0 | nhu4=0 |
| nhx4=0 | nhx4=0 |
| nhv4=0 | nhv4=0 |
| |
theta_front=999.0 |
| |
theta_rear=999.0 |
| c | c |
| call grndm ( rndm , 3 ) | call grndm ( rndm , 3 ) |
| c | c |
|
|
|
| sectrack=.false. | sectrack=.false. |
| endif | endif |
| open(unit=1,file='hdr_gep.dat',status='old') | open(unit=1,file='hdr_gep.dat',status='old') |
| do i=1,10 |
do i=1,12 |
| read(1,*)e0(i),eang(i),hang(i),targ_thick(i) | read(1,*)e0(i),eang(i),hang(i),targ_thick(i) |
| enddo | enddo |
| einc=e0(ikinsetting) | einc=e0(ikinsetting) |
|
|
|
| $ thtgt,phtgt,ptgt,dptgt) | $ thtgt,phtgt,ptgt,dptgt) |
| | |
| c | c |
| p_spec=ptgt |
c At this point, after using polynomials from Lerose, we should |
| |
c have x,y,theta,phi, and dpmom at the focal plane. |
| c | c |
| |
c write(*,*)'Back from kincalc ...' |
| |
c write(*,*)xtgt,ytgt,thtgt,phtgt,dptgt |
| |
c x(1)=tan(phtgt) |
| |
c x(2)=ytgt |
| |
c x(3)=tan(thtgt) |
| |
c x(4)=dptgt |
| |
y_target=ytgt*100.0 |
| |
c |
| |
c Here is where we call the SIMC routines to calculate the focal plane |
| |
c quantities from the target quantities. |
| |
c |
| |
c xfp=xfinal(x,4)*100.0 |
| |
c phfp=atan(thetaf(x,4))/3.14159265*180.0 |
| |
c yfp=yf(x,4)*100.0 |
| |
c thfp=atan(ff(x,4))/3.14159265*180.0 |
| |
p_spec=ptgt |
| |
th_spec=hrsh_ang |
| |
dpp=dptgt*100.00 |
| |
x=xtgt |
| |
y=ytgt |
| |
z=0.0 |
| |
dxdz=tan(phtgt*3.14159265/180.0) |
| |
dydz=tan(thtgt*3.14159265/180.0) |
| if(nevent.le.2)write(*,*)'******** mom = ',p_spec,' ******' | if(nevent.le.2)write(*,*)'******** mom = ',p_spec,' ******' |
| |
c write(*,*)x,y,z |
| |
c write(*,*)dxdz,dydz |
| c | c |
| call grndm(rndm,3) |
c Call to SIMC routine |
| dptgt2=rndm(1)*0.200-0.100 |
|
| xfp=rndm(2)*120.00-60.00 |
|
| phfp=rndm(3)*80.0-40.0 |
|
| c phfp=3.0 |
|
| c phfp=rndm(3)*0.00001-0.000005 |
|
| c |
|
| call grndm(rndm,3) |
|
| yfp=rndm(1)*60.00-30.00 |
|
| thfp=rndm(2)*80.0-40.0 |
|
| c thfp=3.0 |
|
| c thfp=rndm(3)*0.00001-0.000005 |
|
| c | c |
| |
decay_flag=.false. |
| |
wcs_flag=.false. |
| |
ms_flag=.false. |
| |
m2=(938.2796**2) |
| |
call mc_hms (p_spec, th_spec, dpp, x, y, z, dxdz, dydz, |
| |
> x_fp, dx_fp, y_fp, dy_fp, m2, |
| |
> ms_flag, wcs_flag, decay_flag, resmult, fry, |
| |
> ok_spec, pathlen) |
| |
c write(*,*)'Momentum = ',p_spec,' Angle = ',th_spec |
| |
c write(*,*)x,y,z |
| |
c write(*,*)dxdz,dydz |
| |
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 |
| |
|
| |
xfp=x_fp |
| |
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,*)'x,phi,y,theta,ok_spec =' |
| c write(6,*)xfp,phfp,yfp,thfp,ok_spec | c write(6,*)xfp,phfp,yfp,thfp,ok_spec |
| c | c |
|
|
|
| ntuple_array(14)=real(phfp) | ntuple_array(14)=real(phfp) |
| c | c |
| pcentral=ptgt | pcentral=ptgt |
| pfp=dptgt2 |
pfp=dptgt |
| |
c |
| |
c |
| |
c write(*,*)'Focal plane quantities ...' |
| |
c write(*,*)xfp,yfp,thfp,phfp,pcentral,pfp |
| |
|
| |
c xfp=xtgt ! in cm |
| |
c yfp=ytgt ! in cm |
| |
c thfp=thtgt/3.14159265*180.0 ! in degrees |
| |
c phfp=phtgt/3.14159265*180.0 ! in degrees |
| |
c pcentral=ptgt ! in MeV/c |
| |
c pfp=dptgt ! fraction of central momentum |
| c | c |
| c Finally, we convert this information over to a format that GEANT likes. | c Finally, we convert this information over to a format that GEANT likes. |
| c | c |
|
|
|
| tinit=ea-938.2796 | tinit=ea-938.2796 |
| vert(1)=xfp | vert(1)=xfp |
| vert(2)=yfp | vert(2)=yfp |
| vert(3)=160.0 |
vert(3)=0.0 |
| ntbeam = 0.0 | ntbeam = 0.0 |
| nttarg = 0.0 | nttarg = 0.0 |
| ubuf(1) = 0.0 | ubuf(1) = 0.0 |
Return to
gukine.f
CVS log
Up to [HallC] / geant_gep / src