c ----------- Start of Geant variables ---------------------------- COMMON/GCCUTS/CUTGAM,CUTELE,CUTNEU,CUTHAD,CUTMUO,BCUTE,BCUTM + ,DCUTE ,DCUTM ,PPCUTM,TOFMAX,GCUTS(5) COMMON/GCFLAG/IDEBUG,IDEMIN,IDEMAX,ITEST,IDRUN,IDEVT,IEORUN + ,IEOTRI,IEVENT,ISWIT(10),IFINIT(20),NEVENT,NRNDM(2) COMMON/GCFLAX/BATCH, NOLOG LOGICAL BATCH, NOLOG COMMON/GCJLOC/NJLOC(2),JTM,JMA,JLOSS,JPROB,JMIXT,JPHOT,JANNI + ,JCOMP,JBREM,JPAIR,JDRAY,JPFIS,JMUNU,JRAYL + ,JMULOF,JCOEF,JRANG C COMMON/GCJLCK/NJLCK(2),JTCKOV,JABSCO,JEFFIC,JINDEX,JCURIN + ,JPOLAR,JTSTRA,JTSTCO,JTSTEN,JTASHO C EQUIVALENCE (JLASTV,JTSTEN) PARAMETER (MAXJMP=30) COMMON/GCJUMP/JUDCAY, JUDIGI, JUDTIM, JUFLD , JUHADR, JUIGET, + JUINME, JUINTI, JUKINE, JUNEAR, JUOUT , JUPHAD, + JUSKIP, JUSTEP, JUSWIM, JUTRAK, JUTREV, JUVIEW, + JUPARA DIMENSION JMPADR(MAXJMP) EQUIVALENCE (JMPADR(1), JUDCAY) COMMON/GCKINE/IKINE,PKINE(10),ITRA,ISTAK,IVERT,IPART,ITRTYP + ,NAPART(5),AMASS,CHARGE,TLIFE,VERT(3),PVERT(4),IPAOLD INTEGER MXGKIN PARAMETER (MXGKIN=100) PARAMETER (MULTRA=50) CHARACTER*4 GNASH, GNNVV, GNVNV COMMON/GCMUTR/NCVOLS,KSHIFT,NSHIFT,ICUBE,NAIN,JJJ, + NIET,IOLDSU,IVOOLD,IWPOIN,IHPOIN,IVECVO(100), + PORGX,PORGY,PORGZ,POX(15),POY(15),POZ(15),GBOOM, + PORMIR(18),PORMAR(18),IPORNT, + ICGP,CLIPMI(6),CLIPMA(6), + ABCD(4),BMIN(6),BMAX(6),CGB(16000),CGB1(16000), + GXMIN(MULTRA),GXMAX(MULTRA),GYMIN(MULTRA), + GYMAX(MULTRA),GZMIN(MULTRA),GZMAX(MULTRA), + GXXXX(MULTRA),GYYYY(MULTRA),GZZZZ(MULTRA) COMMON/GCMUTC/ GNASH(MULTRA),GNNVV(MULTRA),GNVNV(MULTRA) COMMON/GCKING/KCASE,NGKINE,GKIN(5,MXGKIN), + TOFD(MXGKIN),IFLGK(MXGKIN) INTEGER KCASE,NGKINE ,IFLGK,MXPHOT,NGPHOT REAL GKIN,TOFD,XPHOT C PARAMETER (MXPHOT=800) COMMON/GCKIN2/NGPHOT,XPHOT(11,MXPHOT) C COMMON/GCKIN3/GPOS(3,MXGKIN) REAL GPOS COMMON/GCLIST/NHSTA,NGET ,NSAVE,NSETS,NPRIN,NGEOM,NVIEW,NPLOT + ,NSTAT,LHSTA(20),LGET (20),LSAVE(20),LSETS(20),LPRIN(20) + ,LGEOM(20),LVIEW(20),LPLOT(20),LSTAT(20) COMMON/GCMATE/NMAT,NAMATE(5),A,Z,DENS,RADL,ABSL COMMON/GCMULO/SINMUL(101),COSMUL(101),SQRMUL(101),OMCMOL,CHCMOL + ,EKMIN,EKMAX,NEKBIN,NEK1,EKINV,GEKA,GEKB,EKBIN(200),ELOW(200) COMMON/GCMZFO/IOMATE,IOPART,IOTMED,IOSEJD,IOSJDD,IOSJDH,IOSTAK + ,IOMZFO(13) COMMON/GCNUM/NMATE ,NVOLUM,NROTM,NTMED,NTMULT,NTRACK,NPART + ,NSTMAX,NVERTX,NHEAD,NBIT COMMON /GCNUMX/ NALIVE,NTMSTO c COMMON/GCOMIS/ICOMIS,JUINIT,JUGEOM,JUKINE,JUSTEP,JUOUT,JULAST COMMON/GCONST/PI,TWOPI,PIBY2,DEGRAD,RADDEG,CLIGHT,BIG,EMASS COMMON/GCONSX/EMMU,PMASS,AVO COMMON/GCPHYS/IPAIR,SPAIR,SLPAIR,ZINTPA,STEPPA + ,ICOMP,SCOMP,SLCOMP,ZINTCO,STEPCO + ,IPHOT,SPHOT,SLPHOT,ZINTPH,STEPPH + ,IPFIS,SPFIS,SLPFIS,ZINTPF,STEPPF + ,IDRAY,SDRAY,SLDRAY,ZINTDR,STEPDR + ,IANNI,SANNI,SLANNI,ZINTAN,STEPAN + ,IBREM,SBREM,SLBREM,ZINTBR,STEPBR + ,IHADR,SHADR,SLHADR,ZINTHA,STEPHA + ,IMUNU,SMUNU,SLMUNU,ZINTMU,STEPMU + ,IDCAY,SDCAY,SLIFE ,SUMLIF,DPHYS1 + ,ILOSS,SLOSS,SOLOSS,STLOSS,DPHYS2 + ,IMULS,SMULS,SOMULS,STMULS,DPHYS3 + ,IRAYL,SRAYL,SLRAYL,ZINTRA,STEPRA COMMON/GCPHLT/ILABS,SLABS,SLLABS,ZINTLA,STEPLA + ,ISYNC + ,ISTRA ccc GCSCANScan geometry control parameters PARAMETER (MSLIST=32,MAXMDT=3) COMMON/GCSCAN/SCANFL,NPHI,PHIMIN,PHIMAX,NTETA,TETMIN,TETMAX, + MODTET,IPHIMI,IPHIMA,IPHI1,IPHIL,NSLMAX, + NSLIST,ISLIST(MSLIST),VSCAN(3),FACTX0,FACTL, + FACTR,IPHI,ITETA,ISCUR,SX0,SABS,TETMID(MAXMDT), + TETMAD(MAXMDT) + ,SX0S,SX0T,SABSS,SABST,FACTSF + ,DLTPHI,DLTETA,DPHIM1,DTETM1 + ,FCX0M1,FCLLM1,FCRRM1 LOGICAL SCANFL COMMON/GCSCAC/SFIN,SFOUT CHARACTER*80 SFIN,SFOUT ccc GCSETSIdentification of current sensitive detector COMMON/GCSETS/IHSET,IHDET,ISET,IDET,IDTYPE,NVNAME,NUMBV(20) ccc GCSTAKControl variables for parallel tracking PARAMETER (NWSTAK=12,NWINT=11,NWREAL=12,NWTRAC=NWINT+NWREAL+5) COMMON /GCSTAK/ NJTMAX, NJTMIN, NTSTKP, NTSTKS, NDBOOK, NDPUSH, + NJFREE, NJGARB, NJINVO, LINSAV(15), LMXSAV(15) ccc GCTMEDArray of current tracking medium parameters COMMON/GCTMED/NUMED,NATMED(5),ISVOL,IFIELD,FIELDM,TMAXFD,STEMAX + ,DEEMAX,EPSIL,STMIN,CFIELD,PREC,IUPD,ISTPAR,NUMOLD COMMON/GCTLIT/THRIND,PMIN,DP,DNDL,JMIN,ITCKOV,IMCKOV,NPCKOV ccc GCTRAKTrack parameters at the end of the current step PARAMETER (MAXMEC=30) COMMON/GCTRAK/VECT(7),GETOT,GEKIN,VOUT(7),NMEC,LMEC(MAXMEC) + ,NAMEC(MAXMEC),NSTEP ,MAXNST,DESTEP,DESTEL,SAFETY,SLENG + ,STEP ,SNEXT ,SFIELD,TOFG ,GEKRAT,UPWGHT,IGNEXT,INWVOL + ,ISTOP ,IGAUTO,IEKBIN, ILOSL, IMULL,INGOTO,NLDOWN,NLEVIN + ,NLVSAV,ISTORY PARAMETER (MAXME1=30) COMMON/GCTPOL/POLAR(3), NAMEC1(MAXME1) ccc List of mechanisms active in the current step. CHARACTER*4 MEC(MAXMEC),MEC1(MAXME1),DFLT(2) PARAMETER (LEFTM1=MAXME1-9) DATA MEC/'NEXT','MULS','LOSS','FIEL','DCAY','PAIR','COMP','PHOT' + ,'BREM','DRAY','ANNI','HADR','ECOH','EVAP','FISS','ABSO' + ,'ANNH','CAPT','EINC','INHE','MUNU','TOFM','PFIS','SCUT' + ,'RAYL','PARA','PRED','LOOP','NULL','STOP'/ DATA MEC1/'LABS','LREF','SMAX','SCOR','CKOV','REFL','REFR', + 'SYNC','STRA',LEFTM1*' '/ ccc GCUNITDescription of logical units COMMON/GCUNIT/LIN,LOUT,NUNITS,LUNITS(5) INTEGER LIN,LOUT,NUNITS,LUNITS COMMON/GCMAIL/CHMAIL CHARACTER*132 CHMAIL ccc GCVOLUCurrent geometrical information COMMON/GCVOLU/NLEVEL,NAMES(15),NUMBER(15), + LVOLUM(15),LINDEX(15),INFROM,NLEVMX,NLDEV(15),LINMX(15), + GTRAN(3,15),GRMAT(10,15),GONLY(15),GLX(3) COMMON/GCDRAW/NUMNOD,MAXNOD,NUMND1,LEVVER,LEVHOR,MAXV,IPICK, + MLEVV,MLEVH,NWCUT,JNAM,JMOT,JXON,JBRO,JDUP,JSCA,JDVM,JPSM, + JNAM1,JMOT1,JXON1,JBRO1,JDUP1,JSCA1,JULEV,JVLEV, + LOOKTB(16), + GRMAT0(10),GTRAN0(3),IDRNUM,GSIN(41),GCOS(41),SINPSI,COSPSI, + GTHETA,GPHI,GPSI,GU0,GV0,GSCU,GSCV,NGVIEW, + ICUTFL,ICUT,CTHETA,CPHI,DCUT,NSURF,ISURF, + GZUA,GZVA,GZUB,GZVB,GZUC,GZVC,PLTRNX,PLTRNY, + LINATT,LINATP,ITXATT,ITHRZ,IPRJ,DPERS,ITR3D,IPKHIT,IOBJ,LINBUF, + MAXGU,MORGU,MAXGS,MORGS,MAXTU,MORTU,MAXTS,MORTS, + IGU,IGS,ITU,ITS,NKVIEW,IDVIEW, + NOPEN,IGMR,IPIONS,ITRKOP,IHIDEN, + ZZFU,ZZFV,MYISEL, + DDUMMY(15) c common/gcbank/q(ngbank) c ------- End of Geant variables ------------------- c integer nwpawc integer lrecl integer istat real n_tuple(6) integer icycle integer nexits,nstops,nmuons integer numeropart real momentpart character*10 chtags(6) c parameter (ngbank=1000000,nwpawc=1000000) c parameter (lrecl=1024) c common/pawc/h(nwpawc) common/my/istat,n_tuple,icycle,numeropart,momentpart common/my2/nexits(10000),nstops(10000),nmuons(10000) ! num of moth vol exits and num of "stops" common/ENLOSS/Esplane(60),Efplane(60),Ealplane(60) common/STRIP/Estrip(60,24),Estrip2(60,24),Emip(60,24) ! Energy loss at each strip data chtags /'EIR','ESC','EAL','MIP','YY','ZZ'/