version 1.9, 2005/12/16 07:23:56
|
version 1.18, 2006/02/12 06:18:50
|
|
|
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) |