1 C ****COMMON BLOCK FILE FOR HERWIG VERSION 6.5****
3 C ALTERATIONS: Layout completely overhauled for 5.9
6 C New common blocks added for version 6.1:
7 C HWCLUS,HWSUSY,HWRPAR,HWMINB
9 C New variables added for version 6.1:
10 C OMHMIX,ET2MIX,PH3MIX,IOP4JT,NPRFMT,
11 C PRNDEF,PRNTEX,PRNWEB,EFFMIN,GCUTME,
12 C IOP4JT,NPRFMT see HWPRAM
13 C Y4JT,DURHAM see HWHARD
14 C QORQQB,QBORQQ see HWPROP
17 C PPCL,NCL,IDCL see HWCLUS
18 C TANB,ALPHAH,COSBPA,SINBPA,COSBMA,
19 C SINBMA,COSA,SINA,COSB,SINB,COTB,
20 C ZMIXSS,ZMXNSS,ZSGNSS,LFCH,RFCH,
21 C SLFCH,SRFCH,WMXUSS,WMXVSS,WSGNSS,
22 C QMIXSS,LMIXSS,THETAT,THETAB,THETAL,
23 C ATSS,ABSS,ALSS,MUSS,FACTSS,GHWWSS,
24 C GHZZSS,GHDDSS,GHUUSS,GHWHSS,GHSQSS,
25 C XLMNSS,RMMNSS,DMSSM,SENHNC,
26 C SSPARITY,SUSYIN see HWSUSY
27 C LAMDA1,LAMDA2,LAMDA3,HRDCOL,RPARTY,
29 C PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,
30 C PMBM1,PMBM2,PMBP1,PMBP2,PMBP3 see HWMINB
32 C New parameters added for version 6.1:
35 C Parameter NMXRES raised to 500
37 C Scalar variables changed to arrays of size 2:
38 C CLSMR,PSPLT,CLDIR see HWPRAM
42 C New common blocks added for version 6.2:
47 C New common block added for version 6.202:
49 C which contains all other new variables, parameters and
50 C control flags introduced since version 6.1, so that
51 C other common blocks become identical to version 6.1
53 C New parameters added for version 6.2:
54 C VIPWID,DXRCYL,DXZMAX,DXRSPH,LRSUSY see HWIGIN
55 C GRVLAM,EMGRV,GAMGRV see HWHGRV
57 C New control flags added for version 6.2:
60 C IMSSM,IHIGGS see HWUINC
62 C New variable added for version 6.2:
67 C NMXHEP raised to 4000 for LHC studies
71 DOUBLE PRECISION ZERO,ONE,TWO,THREE,FOUR,HALF
72 PARAMETER (ZERO =0.D0, ONE =1.D0, TWO =2.D0,
73 & THREE=3.D0, FOUR=4.D0, HALF=0.5D0)
76 & ACCUR,AFCH,ALPFAC,ALPHEM,ANOMSC,ASFIXD,AVWGT,B1LIM,BETAF,BRFRAC,
77 & BRHIG,BTCLM,CAFAC,CFFAC,CLDKWT,CLMAX,CLPOW,CLQ,CLSMR,CMMOM,COSS,
78 & COSTH,CSPEED,CTHRPW,CTMAX,DECPAR,DECWT,DISF,DKLTM,EBEAM1,EBEAM2,
79 & EMLST,EMMAX,EMMIN,EMPOW,EMSCA,ENHANC,ENSOF,EPOLN,ETAMIX,EVWGT,
80 & EXAG,F0MIX,F1MIX,F2MIX,GAMH,GAMMAX,GAMW,GAMWT,GAMZ,GAMZP,GCOEF,
81 & GEV2NB,GEV2MM,GPOLN,H1MIX,HBAR,HARDST,OMEGA0,PBEAM1,PBEAM2,PDIQK,
82 & PGSMX,PGSPL,PHEP,PHIMIX,PHIPAR,PHOMAS,PIFAC,PLTCUT,PPAR,PPOLN,
83 & PRECO,PRSOF,PSPLT,PTINT,PTMAX,PTMIN,PTPOW,PTRMS,PXRMS,PWT,Q2MAX,
84 & Q2MIN,Q2POW,Q2WWMN,Q2WWMX,QCDL3,QCDL5,QCDLAM,QDIQK,QEV,QFCH,QG,
85 & QLIM,QSPAC,QV,QWT,REPWT,RESN,RHOHEP,RHOPAR,RLTIM,RMASS,RMIN,
86 & RSPIN,SCABI,SINS,SNGWT,SWEIN,SWTEF,SUD,THMAX,TLOUT,TMTOP,TMNISR,
87 & TQWT,VCKM,VFCH,VGCUT,VHEP,VMIN2,VPAR,VPCUT,VQCUT,VTXPIP,VTXQDK,
88 & WBIGST,WGTMAX,WGTSUM,WHMIN,WSQSUM,XFACT,XLMIN,XMIX,XMRCT,XX,
89 & XXMIN,YBMAX,YBMIN,YJMAX,YJMIN,YMIX,YMRCT,YWWMAX,YWWMIN,ZBINM,
90 & ZJMAX,ZMXISR,Y4JT,EFFMIN,PPCL,
91 & TANB,ALPHAH,COSBPA,SINBPA,COSBMA,SINBMA,COSA,SINA,COSB,SINB,COTB,
92 & ZMIXSS,ZMXNSS,ZSGNSS,LFCH,RFCH,SLFCH,SRFCH, WMXUSS,WMXVSS,WSGNSS,
93 & QMIXSS,LMIXSS,THETAT,THETAB,THETAL,ATSS,ABSS,ALSS,MUSS,FACTSS,
94 & GHWWSS,GHZZSS,GHDDSS,GHUUSS,GHWHSS,GHSQSS
95 C--fix by PR 12/7/02 to avoid problems with nag compiler
97 & XLMNSS,RMMNSS,DMSSM,SENHNC,SSPARITY,LAMDA1,LAMDA2,LAMDA3,
98 & PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,PMBM1,PMBM2,PMBP1,PMBP2,PMBP3,
99 & OMHMIX,ET2MIX,PH3MIX,GCUTME
102 & CLDIR,IAPHIG,IBRN,IBSH,ICHRG,ICO,IDCMF,IDHEP,IDHW,IDK,IDKPRD,IDN,
103 & IDPAR,IDPDG,IERROR,IFLAV,IFLMAX,IFLMIN,IHPRO,IMQDK,INHAD,INTER,
104 & IOPDKL,IOPHIG,IOPREM,IPART1,IPART2,IPRINT,IPRO,IPROC,ISLENT,
105 & ISPAC,ISTAT,ISTHEP,ISTPAR,JCOPAR,JDAHEP,JDAPAR,JMOHEP,JMOPAR,
106 & JNHAD,LNEXT,LOCN,LOCQ,LRSUD,LSTRT,LWEVT,LWSUD,MAPQ,MAXER,MAXEV,
107 & MAXFL,MAXPR,MODBOS,MODMAX,MODPDF,NBTRY,NCLDK,NCOLO,NCTRY,NDKYS,
108 & NDTRY,NETRY,NEVHEP,NEVPAR,NFLAV,NGSPL,NHEP,NME,NMODES,NMXCDK,
109 & NMXDKS,NMXHEP,NMXJET,NMXMOD,NMXPAR,NMXQDK,NMXRES,NMXSUD,NPAR,
110 & NPRODS,NQDK,NQEV,NRES,NRN,NSPAC,NSTRU,NSTRY,NSUD,NUMER,NUMERU,
111 & NWGTS,NZBIN,SUDORD,IOP4JT,HRDCOL,NMXCL,NCL,IDCL,NPRFMT,NRECO
114 & AZSOFT,AZSPIN,BGSHAT,BREIT,CLRECO,COLISR,DKPSET,FROST,FSTEVT,
115 & FSTWGT,GENEV,GENSOF,HARDME,HVFCEN,MAXDKL,MIXING,NOSPAC,NOWGT,
116 & PRNDEC,PIPSMR,PRVTX,RSTAB,SOFTME,TMPAR,TPOL,USECMF,VTOCDK,VTORDK,
117 & ZPRIME,RPARTY,COLUPD,PRNDEF,PRNTEX,PRNWEB,DURHAM,SUSYIN,
129 C New standard event common
130 PARAMETER (NMXHEP=4000)
131 COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
132 & JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
134 C Beams, process and number of events
135 COMMON/HWBEAM/IPART1,IPART2
136 COMMON/HWBMCH/PART1,PART2
137 COMMON/HWPROC/EBEAM1,EBEAM2,PBEAM1,PBEAM2,IPROC,MAXEV
139 C Basic parameters (and quantities derived from them)
140 COMMON/HWPRAM/AFCH(16,2),ALPHEM,B1LIM,BETAF,BTCLM,CAFAC,CFFAC,
141 & CLMAX,CLPOW,CLSMR(2),CSPEED,ENSOF,ETAMIX,F0MIX,F1MIX,F2MIX,GAMH,
142 & GAMW,GAMZ,GAMZP,GEV2NB,H1MIX,PDIQK,PGSMX,PGSPL(4),PHIMIX,PIFAC,
143 & PRSOF,PSPLT(2),PTRMS,PXRMS,QCDL3,QCDL5,QCDLAM,QDIQK,QFCH(16),QG,
144 & QSPAC,QV,SCABI,SWEIN,TMTOP,VFCH(16,2),VCKM(3,3),VGCUT,VQCUT,
145 & VPCUT,ZBINM,EFFMIN,OMHMIX,ET2MIX,PH3MIX,GCUTME,
146 & IOPREM,IPRINT,ISPAC,LRSUD,LWSUD,MODPDF(2),NBTRY,NCOLO,NCTRY,
147 & NDTRY,NETRY,NFLAV,NGSPL,NSTRU,NSTRY,NZBIN,IOP4JT(2),NPRFMT,
148 & AZSOFT,AZSPIN,CLDIR(2),HARDME,NOSPAC,PRNDEC,PRVTX,SOFTME,ZPRIME,
149 & PRNDEF,PRNTEX,PRNWEB
151 COMMON/HWPRCH/AUTPDF(2),BDECAY
153 C Parton shower common (same format as /HEPEVT/)
154 PARAMETER (NMXPAR=500)
155 COMMON/HWPART/NEVPAR,NPAR,ISTPAR(NMXPAR),IDPAR(NMXPAR),
156 & JMOPAR(2,NMXPAR),JDAPAR(2,NMXPAR),PPAR(5,NMXPAR),VPAR(4,NMXPAR)
158 C Parton polarization common
159 COMMON/HWPARP/DECPAR(2,NMXPAR),PHIPAR(2,NMXPAR),RHOPAR(2,NMXPAR),
162 C Electroweak boson common
163 PARAMETER (MODMAX=50)
164 COMMON/HWBOSC/ALPFAC,BRHIG(12),ENHANC(12),GAMMAX,RHOHEP(3,NMXHEP),
165 & IOPHIG,MODBOS(MODMAX)
167 C Parton colour common
168 COMMON/HWPARC/JCOPAR(4,NMXPAR)
170 C other HERWIG branching, event and hard subprocess common blocks
171 COMMON/HWBRCH/ANOMSC(2,2),HARDST,PTINT(3,2),XFACT,INHAD,JNHAD,
172 & NSPAC(7),ISLENT,BREIT,FROST,USECMF
174 COMMON/HWEVNT/AVWGT,EVWGT,GAMWT,TLOUT,WBIGST,WGTMAX,WGTSUM,WSQSUM,
175 & IDHW(NMXHEP),IERROR,ISTAT,LWEVT,MAXER,MAXPR,NOWGT,NRN(2),NUMER,
176 & NUMERU,NWGTS,GENSOF
178 COMMON/HWHARD/ASFIXD,CLQ(7,6),COSS,COSTH,CTMAX,DISF(13,2),EMLST,
179 & EMMAX,EMMIN,EMPOW,EMSCA,EPOLN(3),GCOEF(7),GPOLN,OMEGA0,PHOMAS,
180 & PPOLN(3),PTMAX,PTMIN,PTPOW,Q2MAX,Q2MIN,Q2POW,Q2WWMN,Q2WWMX,QLIM,
181 & SINS,THMAX,Y4JT,TMNISR,TQWT,XX(2),XLMIN,XXMIN,YBMAX,YBMIN,YJMAX,
182 & YJMIN,YWWMAX,YWWMIN,WHMIN,ZJMAX,ZMXISR,IAPHIG,IBRN(2),IBSH,
183 & ICO(10),IDCMF,IDN(10),IFLMAX,IFLMIN,IHPRO,IPRO,MAPQ(6),MAXFL,
184 & BGSHAT,COLISR,FSTEVT,FSTWGT,GENEV,HVFCEN,TPOL,DURHAM
186 C Arrays for particle properties (NMXRES = max no of particles defined)
187 PARAMETER(NMXRES=500)
188 COMMON/HWPROP/RLTIM(0:NMXRES),RMASS(0:NMXRES),RSPIN(0:NMXRES),
189 & ICHRG(0:NMXRES),IDPDG(0:NMXRES),IFLAV(0:NMXRES),NRES,
190 & VTOCDK(0:NMXRES),VTORDK(0:NMXRES),
191 & QORQQB(0:NMXRES),QBORQQ(0:NMXRES)
193 COMMON/HWUNAM/RNAME(0:NMXRES),TXNAME(2,0:NMXRES)
195 C Arrays for particle decays (NMXDKS = max total no of decays,
196 C NMXMOD = max no of modes for a particle)
197 PARAMETER(NMXDKS=4000,NMXMOD=200)
198 COMMON/HWUPDT/BRFRAC(NMXDKS),CMMOM(NMXDKS),DKLTM(NMXRES),
199 & IDK(NMXDKS),IDKPRD(5,NMXDKS),LNEXT(NMXDKS),LSTRT(NMXRES),NDKYS,
200 & NME(NMXDKS),NMODES(NMXRES),NPRODS(NMXDKS),DKPSET,RSTAB(0:NMXRES)
202 C Weights used in cluster decays
203 COMMON/HWUWTS/REPWT(0:3,0:4,0:4),SNGWT,DECWT,QWT(3),PWT(12),
206 C Parameters for cluster decays (NMXCDK = max total no of cluster
208 PARAMETER(NMXCDK=4000)
209 COMMON/HWUCLU/CLDKWT(NMXCDK),CTHRPW(12,12),PRECO,RESN(12,12),
210 & RMIN(12,12),LOCN(12,12),NCLDK(NMXCDK),NRECO,CLRECO
212 C Variables controling mixing and vertex information
213 C--VTXPIP should have been a 5-vector, problems with NAG compiler
214 COMMON/HWDIST/EXAG,GEV2MM,HBAR,PLTCUT,VMIN2,VTXPIP(5),XMIX(2),
215 & XMRCT(2),YMIX(2),YMRCT(2),IOPDKL,MAXDKL,MIXING,PIPSMR
217 C Arrays for temporarily storing heavy-b,c-hadrons decaying partonicaly
218 C (NMXQDK = max no such decays in an event)
219 PARAMETER (NMXQDK=20)
220 COMMON/HWQDKS/VTXQDK(4,NMXQDK),IMQDK(NMXQDK),LOCQ(NMXQDK),NQDK
222 C Parameters for Sudakov form factors
223 C (NMXSUD= max no of entries in lookup table)
224 PARAMETER (NMXSUD=1024)
225 COMMON/HWUSUD/ACCUR,QEV(NMXSUD,6),SUD(NMXSUD,6),INTER,NQEV,NSUD,
228 PARAMETER (NMXJET=200)
232 & TANB,ALPHAH,COSBPA,SINBPA,COSBMA,SINBMA,COSA,SINA,COSB,SINB,COTB,
233 & ZMIXSS(4,4),ZMXNSS(4,4),ZSGNSS(4), LFCH(16),RFCH(16),
234 & SLFCH(16,4),SRFCH(16,4), WMXUSS(2,2),WMXVSS(2,2), WSGNSS(2),
235 & QMIXSS(6,2,2),LMIXSS(6,2,2),
236 & THETAT,THETAB,THETAL,ATSS,ABSS,ALSS,MUSS,FACTSS,
237 & GHWWSS(3),GHZZSS(3),GHDDSS(4),GHUUSS(4),GHWHSS(3),
238 & GHSQSS(4,6,2,2),XLMNSS,RMMNSS,DMSSM,SENHNC(24),SSPARITY,SUSYIN
240 C R-Parity violating parameters and colours
241 COMMON /HWRPAR/ LAMDA1(3,3,3),LAMDA2(3,3,3),
242 & LAMDA3(3,3,3),HRDCOL(2,5),RPARTY,COLUPD
244 C Parameters for minimum bias/soft underlying event
246 & PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,PMBM1,PMBM2,PMBP1,PMBP2,PMBP3
248 C Cluster common used by soft event routines
249 PARAMETER (NMXCL=500)
250 COMMON/HWCLUS/PPCL(5,NMXCL),IDCL(NMXCL),NCL
252 C Parameters for resonant graviton production
253 DOUBLE PRECISION GRVLAM,EMGRV,GAMGRV
254 COMMON/HWGRAV/GRVLAM,EMGRV,GAMGRV
256 C Other new parameters for version 6.2
257 DOUBLE PRECISION VIPWID,DXRCYL,DXZMAX,DXRSPH
259 INTEGER IMSSM,IHIGGS,PARITY,LRSUSY
260 COMMON/HW6202/VIPWID(3),DXRCYL,DXZMAX,DXRSPH,WZRFR,FIX4JT,
261 & IMSSM,IHIGGS,PARITY,LRSUSY
263 C New parameters for version 6.203
264 DOUBLE PRECISION ABWGT,ABWSUM,AVABW
265 INTEGER NNEGWT,NNEGEV
267 COMMON/HW6203/ABWGT,ABWSUM,AVABW,NNEGWT,NNEGEV,NEGWTS
269 C New parameters for version 6.3
270 INTEGER IMAXCH,IMAXOP
271 PARAMETER (IMAXCH=20,IMAXOP=40)
272 DOUBLE PRECISION MJJMIN,CHNPRB(IMAXCH)
274 LOGICAL OPTM,CHON(IMAXCH)
275 COMMON/HW6300/MJJMIN,CHNPRB,IOPSTP,IOPSH,OPTM,CHON
276 C New PDF's for version 6.3
277 INTEGER NXMRS,NQMRS,NPMRS
278 PARAMETER(NXMRS=49,NQMRS=37,NPMRS=8)
279 DOUBLE PRECISION FMRS(3,NPMRS,NXMRS,NQMRS+1)
281 C Circe interface for version 6.3
282 INTEGER CIRCOP,CIRCAC,CIRCVR,CIRCRV,CIRCCH
283 COMMON /HWCIRC/CIRCOP,CIRCAC,CIRCVR,CIRCRV,CIRCCH
284 C New parameters and commons for spin correlations
285 C--constants for the arrays
286 INTEGER NMXSPN,NCFMAX
287 PARAMETER(NMXSPN=50,NCFMAX=3)
288 INTEGER NMODE2,NMODE3,NDIAGR,NMODEB,NMODE4
289 PARAMETER(NMODE2=500,NMODE3=500,NDIAGR=8,NMODEB=50,NMODE4=4)
290 C--common block for X --> X gauge boson
291 DOUBLE PRECISION ABMODE(2,NMODEB),BBMODE(2,12,NMODEB),
292 & PBMODE(12,NMODEB),WTBMAX(12,NMODEB)
293 INTEGER IDBPRT(NMODEB),IBMODE(NMODEB),IBDRTP(NMODEB),NBMODE
294 COMMON /HWDSPB/ABMODE,BBMODE,PBMODE,WTBMAX,IDBPRT,IBDRTP,IBMODE,
296 C--common block for two body decays
297 DOUBLE PRECISION A2MODE(2,NMODE2),P2MODE(NMODE2),WT2MAX(NMODE2)
298 INTEGER ID2PRT(NMODE2),I2DRTP(NMODE2),N2MODE
299 COMMON /HWDSP2/A2MODE,P2MODE,WT2MAX,ID2PRT,I2DRTP,N2MODE
300 C--common block for three body decays
301 DOUBLE PRECISION A3MODE(2,NDIAGR,NMODE3),B3MODE(2,NDIAGR,NMODE3),
302 & P3MODE(NMODE3),WT3MAX(NMODE3),SPN3CF(NCFMAX,NCFMAX,NMODE3)
303 INTEGER ID3PRT(NMODE3),I3MODE(NDIAGR,NMODE3),
304 & I3DRTP(NDIAGR,NMODE3),N3MODE,NDI3BY(NMODE3),N3NCFL(NMODE3),
305 & I3DRCF(NDIAGR,NMODE3)
306 COMMON /HWDSP3/A3MODE,B3MODE,P3MODE,WT3MAX,SPN3CF,ID3PRT,I3MODE,
307 & I3DRTP,N3MODE,NDI3BY,N3NCFL,I3DRCF
308 C--common block for four body decays
309 DOUBLE PRECISION A4MODE(2,12,NMODE4),B4MODE(2,12,NMODE4),
310 & P4MODE(12,12,NMODE4),WT4MAX(12,12,NMODE4)
311 INTEGER ID4PRT(NMODE4),I4MODE(2,NMODE4),N4MODE
312 COMMON /HWDSP4/A4MODE,B4MODE,P4MODE,WT4MAX,ID4PRT,I4MODE,N4MODE
313 C--common block for spin correlations in event
314 INTEGER NDECSY,NSEARCH,LRDEC,LWDEC
315 LOGICAL SYSPIN,THREEB,FOURB
317 COMMON /HWDSPN/NDECSY,NSEARCH,LRDEC,LWDEC,SYSPIN,THREEB,
320 INTEGER IDSPN(NMXSPN),JMOSPN(NMXSPN),JDASPN(2,NMXSPN),NSPN,
321 & ISNHEP(NMXHEP),NSNTRY,NCFL(NMXSPN),SPCOPT
322 DOUBLE COMPLEX MESPN(2,2,2,2,NCFMAX,NMXSPN),RHOSPN(2,2,NMXSPN)
323 DOUBLE PRECISION SPNCFC(NCFMAX,NCFMAX,NMXSPN)
324 LOGICAL DECSPN(NMXSPN)
325 COMMON /HWSPIN/MESPN,RHOSPN,SPNCFC,IDSPN,JMOSPN,JDASPN,
326 & NSPN,ISNHEP,NSNTRY,DECSPN,NCFL,SPCOPT
327 INTEGER JAK1,JAK2,ITDKRC,IFPHOT
328 COMMON /HWSTAU/ JAK1,JAK2,ITDKRC,IFPHOT
330 C--common block for Les Houches interface to store information we need
333 PARAMETER (MAXHRP=100)
334 DOUBLE PRECISION LHWGT(MAXHRP),LHWGTS(MAXHRP),LHMXSM,
335 & LHXSCT(MAXHRP),LHXERR(MAXHRP),LHXMAX(MAXHRP)
336 INTEGER LHIWGT(MAXHRP),ITYPLH,LHNEVT(MAXHRP)
337 LOGICAL LHSOFT,LHGLSF
338 COMMON /HWGUPR/LHWGT,LHWGTS,LHXSCT,LHXERR,LHXMAX,LHMXSM,LHIWGT,
339 & LHNEVT,ITYPLH,LHSOFT,LHGLSF
341 C--common block for HERWIG6.5
344 COMMON /HW6500/ PRESPL
346 COMMON /HW6504/ITOPRD
347 C--the only change for 6505 was to increase MODMAX from 5 to 50
348 DOUBLE PRECISION PDFX0,PDFPOW
349 COMMON /HW6506/PDFX0,PDFPOW
350 C--add new variable to prevent infinite loops in HWDFOR/FIV
352 COMMON /HW6510/NDETRY