C***************************************************************** C* * C* M C R E A C T . C D E * C* * C***************************************************************** C--- C--- Common block containing variables related to MC event generation C--- C--- ================================================== Vertex size an React C--- *** DEFINED BY THE USER through data cards *** C--- C--- VTXSIZ(3) Vertex size (sigma x, sigma y, sigma z) C--- KREAC Reaction type : 1 resonant production (Breit Wigner) C--- -1 non resonant production (gaussian) C--- 0 single particle generation C--- C--- ============================================ Particles involved C--- *** DEFINED BY THE USER through data cards *** C--- C--- NPARTMC Total # of particles C--- LPARTMC(i) array of Hollerith variables containing the names C--- of the particles involved in the reactions C--- PMASS(i) mass of the i-th particle C--- PWIDT(i) Width of the i-th particle C--- PCHAR(i) Charge C--- PCTAU(i) Particle c*Tau (-1. if stable) C--- PKIND(i) Kind = 0 no cal. shower C--- Kind = 1 e.m. shower C--- Kind = 2 hadr shower C--- C--- =================================== Arrays describing the reaction C--- to be generated C--- *** DEFINED BY THE USER through data cards *** C--- C--- NSTEPMC # of step (or cascade decays) trough which the reaction C--- proceeds C--- C--- NBODMC(n) # of particles produced by the n-th step of a cascade decay C--- DECSTEP(m,n) index of particles involved in the n-th step C--- m=1 pointer in LPART to the decaying particle C--- m=2-NBOD(n) pointer in LPART to the daughter particles C--- C--- NMCPUP expected number of Pile-up TOTAL (charged + neutral) C--- multiplicities !!!! REAL !!!! C--- DT_AVER Average Delta T between events C--- C--- NMCTRK Total number(event+pileup) tracks C--- C--- SPGENL(1-2) Single particle generation theta limits (KREAC=0) C--- SPGENL(3-4) Single particle generation phi limits (KREAC=0) C--- SPGENL(5-6) Single particle generation momentum limits (KREAC=0) C--- C--- ================================= NON phase-space generation C--- *** DEFINED BY THE USER through data cards *** C--- C--- USGEN(n) number identifing the user-defined algorithm for event C--- generation at step n (should be provided by a routine C--- USGENEV( UGEN(n) written by the user) C--- GENPAR(i,n) parameters for the event generation at step N C--- used by the routine USGENEV( UGEN(n) ) C--- C--- ================================ Unifor Phase space ============== C--- *** DEFINED BY THE USER through data cards *** C--- C--- MXWEIG(n) Maximum GENBOD weight for this step (used to get uniform C--- Phase space distribution) C--- If 0 is used to determine Max GENBOD weight for the reaction C--- MCWEIG(n) Weight at step N C--- C--- C--- ================= Generated particle parameters C--- (all generated particles) C--- C--- PPLAB(j,i) Px,Py,Pz E,|P| of i-th particle in Lab. Frame C--- PPCM (j,i) " " " " " " in C.M. Frame C--- C--- ============================ Track parameters C--- (all generated particles) C--- C--- TRACKGEN(j,i) j=1,2,3,4 : Theta, phi, E, p C--- j=5,6,7 : x,y,z of production vertex C--- j=8 : decay length (0 for particles not seen C--- in the detector) C--- j=9 : Charge C--- j=10 : Kind (=1: e.m. =2:hadr not e.m.) C--- j=11 : time distance from generated event C--- (<> 0 for pileup) C--- C--- C--- C--- ================================= Common blocks used by GENBOD (CERNLIB W-515) C--- C--- NP # of OUTGOING particles ( 2<= NP <= 18 ) C--- KGENEV 1: No energy dependence of cross section C--- 2: Fermi energy dependence ( NOT USED BY OUR MC) C--- TECM Total energy in C.M. C--- AMASS(i) Mass of i-th outgoing particle C--- MCPCM(1-5,i) Px,Py,Pz,E,P of outgoing particle C--- WT Weight of the event C--- C--- ===================================== Detector crossing C--- C--- Common block containing the information on elements C--- crossed by the track in each volume of the detector C--- and track parameters followed through the detectors C--- Used to fill the detector maps C--- C--- CROSDET( 1, IVOL,IT) E at the exit of the detector C--- 2 th C--- 3 ph C--- 4 p C--- 5 x impact point ( or conversion point) C--- 6 y C--- 7 z C--- 8 x exit point C--- 9 y C--- 10 z C--- 11 L Distance from exit point to decay point C--- 12 Crossed Length inside the detector C--- C============================================================================ INTEGER MCXSTP PARAMETER (MCXSTP=10) INTEGER MCXPRT PARAMETER (MCXPRT=25) INTEGER MCXDCP PARAMETER (MCXDCP=12) INTEGER NPARTMC, LPARTMC(MCXPRT), PCHAR(MCXPRT), PKIND(MCXPRT) +, NSTEPMC, DECSTEP(MCXDCP,MCXSTP) +, NBODMC(MCXSTP), USGEN(MCXSTP), KREAC ,NMCTRK,MCRUN REAL PARTMC(6,MCXPRT), NMCPUP INTEGER IPARTMC(6,MCXPRT) EQUIVALENCE (PARTMC,IPARTMC) REAL PMASS(MCXPRT), PWIDT(MCXPRT), PCTAU(MCXPRT) REAL PPLAB(5,MCXPRT),PPCM(5,MCXPRT),VTXSIZ(3) REAL MXWEIG(MCXSTP),MCWEIG(MCXSTP) REAL TRACKGEN(11,MCXPRT),DT_AVER REAL GENPAR(10,MCXSTP) REAL SPGENL(6) C COMMON /MCREACT/ KREAC,NPARTMC, PARTMC, LPARTMC, MCRUN +, PCHAR, PMASS, PWIDT, PCTAU, PKIND +, NSTEPMC, DECSTEP, VTXSIZ +, NBODMC, MXWEIG, USGEN, MCWEIG, GENPAR +, PPLAB, PPCM, TRACKGEN ,NMCTRK, NMCPUP, DT_AVER COMMON /MCREACT2/ SPGENL C INTEGER NP,KGENEV REAL TECM, AMASS(18), MCPCM(5,18), WT COMMON /GENIN/ NP, TECM, AMASS, KGENEV COMMON /GENOUT/ MCPCM, WT C REAL CROSDET(12,NMXVOL,30) COMMON /DETCROSS/CROSDET