charm decays with photons in the final state and trigger counter ...

charma.uprm.edu

charm decays with photons in the final state and trigger counter ...

79CAlso changed the gap size (to 4 inches)CChanged simulation of V counters and numberingCC Modified: june/19/00c Alexis Paris - Major rewrite to implement a more precise geometryC and efficiency simulation. Changed V numberingCbecause there were really only 12 V counters.C$$INCLUDE ’IOPARAM.INC’$$INCLUDE ’FLATCH.INC’$$INCLUDE ’RGPARAM.INC’$$INCLUDE ’RGCNTRL.INC’$$INCLUDE ’HEADER.INC’$$INCLUDE ’hxvoh_geom.inc’$$INCLUDE ’hxvoh_eff.inc’C Local variablesREAL hlower,hupper,vlower,vupper,x,y,X0INTEGER i,id,nh,nhi,nv,periodlogical firstreal VNEWEFFdata first /.true./if(first)thenhlower = hyc-0.5*hwhupper = hlower+12*hwvlower = vyoff-vlvupper = vyoff+vlfirst = .false.endifC-- Just return if this is a neutral particle&&&IF(ID.EQ. 1.OR.ID.EQ. 4.OR.ID.EQ. 7.OR.ID.EQ.10.OR.ID.EQ.13.OR.ID.EQ.16.OR.ID.EQ.17.OR.ID.EQ.18.OR.ID.EQ.20.OR.ID.EQ.22.OR.ID.EQ.25.OR.ID.EQ.26.OR.ID.EQ.28.OR.ID.EQ.30) THENRETURNEND IFc-- Simulating the HXV run periodsperiod=0if (IRGRUN.ge.5000.and.IRGRUN.le.9000) period = 1if (IRGRUN.ge.9001.and.IRGRUN.le.12000) period = 2if (IRGRUN.ge.12001.and.IRGRUN.le.15000) period = 3if (period.eq.0) returncc--first the horizontal counterscNH=0IF(y.gt.hlower.and.y.lt.hupper) thenNHI=(Y-hlower)/Hw+1IF(NHI.LT.1.OR.NHI.GT.12) THENPRINT *,’ ***WARNING IN HXVSIM INVALID NHI = ’,NHIgoto 10END IFif (x.le.Xhgapn.and.x.gt.(Xhgapn-hl)) thennh=nhi !negative xelseif (x.gt.Xhgapp.and.x.lt.(Xhgapp+hl))thennh=nhi+12 !positive xendifIF(NH.GE.1.and.NH.LE.24) THEN

More magazines by this user
Similar magazines