|
version 1.8, 2005/10/11 22:54:47
|
version 1.16, 2006/02/05 20:32:28
|
|
|
|
| integer i,j,k,ichoice,ichoice2,iremain | integer i,j,k,ichoice,ichoice2,iremain |
| real*8 rotmat,xyz(3),xyznew(3),termang | real*8 rotmat,xyz(3),xyznew(3),termang |
| integer ic,ii,jj | integer ic,ii,jj |
| integer*4 junk1,ikinsetting |
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(10),eang(10),hang(10),targ_thick(10) |
| real*8 xtgt,ytgt,thtgt,phtgt,ptgt,dptgt | real*8 xtgt,ytgt,thtgt,phtgt,ptgt,dptgt |
|
|
|
| c include 'option.h' | c include 'option.h' |
| c | c |
| common/kincom/rotmat(3,3) | common/kincom/rotmat(3,3) |
| |
c |
| |
nhu1=0 |
| |
nhx1=0 |
| |
nhv1=0 |
| |
nhu2=0 |
| |
nhx2=0 |
| |
nhv2=0 |
| |
nhu3=0 |
| |
nhx3=0 |
| |
nhv3=0 |
| |
nhu4=0 |
| |
nhx4=0 |
| |
nhv4=0 |
| |
theta_front=999.0 |
| |
theta_rear=999.0 |
| |
c |
| call grndm ( rndm , 3 ) | call grndm ( rndm , 3 ) |
| c | c |
| 111 format(a80) | 111 format(a80) |
|
|
|
| c write(*,*)'Getting kinematics setting ...' | c write(*,*)'Getting kinematics setting ...' |
| open(unit=1,file='geant_kinematics.dat',type='UNKNOWN') | open(unit=1,file='geant_kinematics.dat',type='UNKNOWN') |
| read(1,*)ikinsetting | read(1,*)ikinsetting |
| |
read(1,*)isectrack |
| close(unit=1) | close(unit=1) |
| |
if(isectrack.eq.1) then |
| |
sectrack=.true. |
| |
else |
| |
sectrack=.false. |
| |
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,10 |
| read(1,*)e0(i),eang(i),hang(i),targ_thick(i) | read(1,*)e0(i),eang(i),hang(i),targ_thick(i) |
|
|
|
| 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 |
|
|
|
| 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 |
Return to
gukine.f
CVS log
Up to [HallC] / geant_gep / src