1 jones 1.1 subroutine gustep
2 c
3 c this subroutine was written by jpsullivan april 21-22, 1993
4 c the tracking realated part is relatively simple -- if the
5 c particle leave the volume called 'targ', throw it away.
6 c it also makes a bunch of histograms
7 c
8 c *keep,gctrak.
9 *-- author :
10 common/gctrak/vect(7),getot,gekin,vout(7),nmec,lmec(30),namec(30)
11 + ,nstep ,maxnst,destep,destel,safety,sleng
12 + ,step ,snext ,sfield
13 + ,tofg ,gekrat,upwght,ignext,inwvol,istop ,idecad,iekbin
14 + , ilosl, imull,ingoto,nldown,nlevin,nlvsav,istory
15 c
16 integer nmec,lmec,namec,nstep ,maxnst,ignext,inwvol,istop
17 + ,idecad,iekbin,ilosl, imull,ingoto,nldown,nlevin
18 + ,nlvsav,istory
19 real vect,getot,gekin,vout,destep,destel,safety,sleng ,step
20 + ,snext,sfield,tofg ,gekrat,upwght
21 c end gctrak
22 jones 1.1 * keep,gcvolu.
23 *-- author :
24 common/gcvolu/nlevel,names(15),number(15),
25 + lvolum(15),lindex(15),infrom,nlevmx,nldev(15),linmx(15),
26 + gtran(3,15),grmat(10,15),gonly(15),glx(3)
27 c
28 integer nlevel,number,lvolum,lindex,infrom,nlevmx,
29 + nldev,linmx
30 character*4 names
31 real gtran,grmat,gonly,glx
32 c end gcvolu
33 c
34 * keep,gcbank.
35 *-- author :
36 integer iq,lq,nzebra,ixstor,ixdiv,ixcons,lmain,lr1
37 integer kwbank,kwwork,iws
38 real gversn,zversn,fendq,ws,q
39 c
40 parameter (kwbank=69000,kwwork=5200)
41 common/gcbank/nzebra,gversn,zversn,ixstor,ixdiv,ixcons,fendq(16)
42 + ,lmain,lr1,ws(kwbank)
43 jones 1.1 dimension iq(2),q(2),lq(8000),iws(2)
44 equivalence (q(1),iq(1),lq(9)),(lq(1),lmain),(iws(1),ws(1))
45 common/gclink/jdigi ,jdraw ,jhead ,jhits ,jkine ,jmate ,jpart
46 + ,jrotm ,jrung ,jset ,jstak ,jgstat,jtmed ,jtrack,jvertx
47 + ,jvolum,jxyz ,jgpar ,jgpar2,jsklt
48 c
49 integer jdigi ,jdraw ,jhead ,jhits ,jkine ,jmate ,jpart
50 + ,jrotm ,jrung ,jset ,jstak ,jgstat,jtmed ,jtrack,jvertx
51 + ,jvolum,jxyz ,jgpar,jgpar2 ,jsklt
52 c
53 * keep,gcking.
54 *-- author :
55 common/gcking/kcase,ngkine,gkin(5,100),tofd(100),iflgk(100)
56 integer kcase,ngkine ,iflgk
57 real gkin,tofd
58 c end gcking
59 c
60 * keep,gckine.
61 *-- author :
62 *-- author :
63 integer ikine,itra,istak,ivert,ipart,itrtyp,napart,ipaold
64 jones 1.1 real pkine,amass,charge,tlife,vert,pvert
65 common/gckine/ikine,pkine(10),itra,istak,ivert,ipart,itrtyp
66 + ,napart(5),amass,charge,tlife,vert(3),pvert(4),ipaold
67 c end gckine
68 c
69 integer ihset,ihdet,iset,idet,idtype,nvname,numbv
70 common/gcsets/ihset,ihdet,iset,idet,idtype,nvname,numbv(20)
71 real x1,y1,z1,lpar,v1,v2,v3,newdist,x1new,y1new,z1new
72 real xstr,ystr,zstr,xstrnew,ystrnew,zstrnew
73 real rotmat2,rotmat3,rotmat4,rotmat1
74 c
75 common/geomstep/rotmat1(3,3),rotmat2(3,3),rotmat3(3,3),
76 & rotmat4(3,3)
77 c
78 include 'fpp_local.h'
79 include 'geant_local.h'
80 c include 'parameter.h'
81 c include 'espace_type.h'
82 c include 'detector.h'
83 c include 'transport.h'
84 c include 'option.h'
85 jones 1.1 c
86 c
87 integer i,make_hist,ioff,ihit,ieffcheck
88 real pt_pi,ppar_pi,arg1,arg2,rapid_pi,pchmb,hits(6)
89 real a,b,c,beta,z,y,straw,ypath,rndm(3),ycompare
|