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,
115 & AZSOFT,AZSPIN,BGSHAT,BREIT,CLRECO,COLISR,DKPSET,FROST,FSTEVT,
116 & FSTWGT,GENEV,GENSOF,HARDME,HVFCEN,MAXDKL,MIXING,NOSPAC,NOWGT,
117 & PRNDEC,PIPSMR,PRVTX,RSTAB,SOFTME,TMPAR,TPOL,USECMF,VTOCDK,VTORDK,
118 & ZPRIME,RPARTY,COLUPD,PRNDEF,PRNTEX,PRNWEB,DURHAM,SUSYIN,
130 C New standard event common
131 PARAMETER (NMXHEP=4000)
132 COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
133 & JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
135 C Beams, process and number of events
136 COMMON/HWBEAM/IPART1,IPART2
137 COMMON/HWBMCH/PART1,PART2
138 COMMON/HWPROC/EBEAM1,EBEAM2,PBEAM1,PBEAM2,IPROC,MAXEV
140 C Basic parameters (and quantities derived from them)
141 COMMON/HWPRAM/AFCH(16,2),ALPHEM,B1LIM,BETAF,BTCLM,CAFAC,CFFAC,
142 & CLMAX,CLPOW,CLSMR(2),CSPEED,ENSOF,ETAMIX,F0MIX,F1MIX,F2MIX,GAMH,
143 & GAMW,GAMZ,GAMZP,GEV2NB,H1MIX,PDIQK,PGSMX,PGSPL(4),PHIMIX,PIFAC,
144 & PRSOF,PSPLT(2),PTRMS,PXRMS,QCDL3,QCDL5,QCDLAM,QDIQK,QFCH(16),QG,
145 & QSPAC,QV,SCABI,SWEIN,TMTOP,VFCH(16,2),VCKM(3,3),VGCUT,VQCUT,
146 & VPCUT,ZBINM,EFFMIN,OMHMIX,ET2MIX,PH3MIX,GCUTME,
147 & IOPREM,IPRINT,ISPAC,LRSUD,LWSUD,MODPDF(2),NBTRY,NCOLO,NCTRY,
148 & NDTRY,NETRY,NFLAV,NGSPL,NSTRU,NSTRY,NZBIN,IOP4JT(2),NPRFMT,
149 & AZSOFT,AZSPIN,CLDIR(2),HARDME,NOSPAC,PRNDEC,PRVTX,SOFTME,ZPRIME,
150 & PRNDEF,PRNTEX,PRNWEB
152 COMMON/HWPRCH/AUTPDF(2),BDECAY
154 C Parton shower common (same format as /HEPEVT/)
155 PARAMETER (NMXPAR=500)
156 COMMON/HWPART/NEVPAR,NPAR,ISTPAR(NMXPAR),IDPAR(NMXPAR),
157 & JMOPAR(2,NMXPAR),JDAPAR(2,NMXPAR),PPAR(5,NMXPAR),VPAR(4,NMXPAR)
159 C Parton polarization common
160 COMMON/HWPARP/DECPAR(2,NMXPAR),PHIPAR(2,NMXPAR),RHOPAR(2,NMXPAR),
163 C Electroweak boson common
164 PARAMETER (MODMAX=50)
165 COMMON/HWBOSC/ALPFAC,BRHIG(12),ENHANC(12),GAMMAX,RHOHEP(3,NMXHEP),
166 & IOPHIG,MODBOS(MODMAX)
168 C Parton colour common
169 COMMON/HWPARC/JCOPAR(4,NMXPAR)
171 C other HERWIG branching, event and hard subprocess common blocks
172 COMMON/HWBRCH/ANOMSC(2,2),HARDST,PTINT(3,2),XFACT,INHAD,JNHAD,
173 & NSPAC(7),ISLENT,BREIT,FROST,USECMF
175 COMMON/HWEVNT/AVWGT,EVWGT,GAMWT,TLOUT,WBIGST,WGTMAX,WGTSUM,WSQSUM,
176 & IDHW(NMXHEP),IERROR,ISTAT,LWEVT,MAXER,MAXPR,NOWGT,NRN(2),NUMER,
177 & NUMERU,NWGTS,GENSOF,EV1PR,EV2PR
179 COMMON/HWHARD/ASFIXD,CLQ(7,6),COSS,COSTH,CTMAX,DISF(13,2),EMLST,
180 & EMMAX,EMMIN,EMPOW,EMSCA,EPOLN(3),GCOEF(7),GPOLN,OMEGA0,PHOMAS,
181 & PPOLN(3),PTMAX,PTMIN,PTPOW,Q2MAX,Q2MIN,Q2POW,Q2WWMN,Q2WWMX,QLIM,
182 & SINS,THMAX,Y4JT,TMNISR,TQWT,XX(2),XLMIN,XXMIN,YBMAX,YBMIN,YJMAX,
183 & YJMIN,YWWMAX,YWWMIN,WHMIN,ZJMAX,ZMXISR,IAPHIG,IBRN(2),IBSH,
184 & ICO(10),IDCMF,IDN(10),IFLMAX,IFLMIN,IHPRO,IPRO,MAPQ(6),MAXFL,
185 & BGSHAT,COLISR,FSTEVT,FSTWGT,GENEV,HVFCEN,TPOL,DURHAM
187 C Arrays for particle properties (NMXRES = max no of particles defined)
188 PARAMETER(NMXRES=500)
189 COMMON/HWPROP/RLTIM(0:NMXRES),RMASS(0:NMXRES),RSPIN(0:NMXRES),
190 & ICHRG(0:NMXRES),IDPDG(0:NMXRES),IFLAV(0:NMXRES),NRES,
191 & VTOCDK(0:NMXRES),VTORDK(0:NMXRES),
192 & QORQQB(0:NMXRES),QBORQQ(0:NMXRES)
194 COMMON/HWUNAM/RNAME(0:NMXRES),TXNAME(2,0:NMXRES)
196 C Arrays for particle decays (NMXDKS = max total no of decays,
197 C NMXMOD = max no of modes for a particle)
198 PARAMETER(NMXDKS=4000,NMXMOD=200)
199 COMMON/HWUPDT/BRFRAC(NMXDKS),CMMOM(NMXDKS),DKLTM(NMXRES),
200 & IDK(NMXDKS),IDKPRD(5,NMXDKS),LNEXT(NMXDKS),LSTRT(NMXRES),NDKYS,
201 & NME(NMXDKS),NMODES(NMXRES),NPRODS(NMXDKS),DKPSET,RSTAB(0:NMXRES)
203 C Weights used in cluster decays
204 COMMON/HWUWTS/REPWT(0:3,0:4,0:4),SNGWT,DECWT,QWT(3),PWT(12),
207 C Parameters for cluster decays (NMXCDK = max total no of cluster
209 PARAMETER(NMXCDK=4000)
210 COMMON/HWUCLU/CLDKWT(NMXCDK),CTHRPW(12,12),PRECO,RESN(12,12),
211 & RMIN(12,12),LOCN(12,12),NCLDK(NMXCDK),NRECO,CLRECO
213 C Variables controling mixing and vertex information
214 C--VTXPIP should have been a 5-vector, problems with NAG compiler
215 COMMON/HWDIST/EXAG,GEV2MM,HBAR,PLTCUT,VMIN2,VTXPIP(5),XMIX(2),
216 & XMRCT(2),YMIX(2),YMRCT(2),IOPDKL,MAXDKL,MIXING,PIPSMR
218 C Arrays for temporarily storing heavy-b,c-hadrons decaying partonicaly
219 C (NMXQDK = max no such decays in an event)
220 PARAMETER (NMXQDK=20)
221 COMMON/HWQDKS/VTXQDK(4,NMXQDK),IMQDK(NMXQDK),LOCQ(NMXQDK),NQDK
223 C Parameters for Sudakov form factors
224 C (NMXSUD= max no of entries in lookup table)
225 PARAMETER (NMXSUD=1024)
226 COMMON/HWUSUD/ACCUR,QEV(NMXSUD,6),SUD(NMXSUD,6),INTER,NQEV,NSUD,
229 PARAMETER (NMXJET=200)
233 & TANB,ALPHAH,COSBPA,SINBPA,COSBMA,SINBMA,COSA,SINA,COSB,SINB,COTB,
234 & ZMIXSS(4,4),ZMXNSS(4,4),ZSGNSS(4), LFCH(16),RFCH(16),
235 & SLFCH(16,4),SRFCH(16,4), WMXUSS(2,2),WMXVSS(2,2), WSGNSS(2),
236 & QMIXSS(6,2,2),LMIXSS(6,2,2),
237 & THETAT,THETAB,THETAL,ATSS,ABSS,ALSS,MUSS,FACTSS,
238 & GHWWSS(3),GHZZSS(3),GHDDSS(4),GHUUSS(4),GHWHSS(3),
239 & GHSQSS(4,6,2,2),XLMNSS,RMMNSS,DMSSM,SENHNC(24),SSPARITY,SUSYIN
241 C R-Parity violating parameters and colours
242 COMMON /HWRPAR/ LAMDA1(3,3,3),LAMDA2(3,3,3),
243 & LAMDA3(3,3,3),HRDCOL(2,5),RPARTY,COLUPD
245 C Parameters for minimum bias/soft underlying event
247 & PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,PMBM1,PMBM2,PMBP1,PMBP2,PMBP3
249 C Cluster common used by soft event routines
250 PARAMETER (NMXCL=500)
251 COMMON/HWCLUS/PPCL(5,NMXCL),IDCL(NMXCL),NCL
253 C Parameters for resonant graviton production
254 DOUBLE PRECISION GRVLAM,EMGRV,GAMGRV
255 COMMON/HWGRAV/GRVLAM,EMGRV,GAMGRV
257 C Other new parameters for version 6.2
258 DOUBLE PRECISION VIPWID,DXRCYL,DXZMAX,DXRSPH
260 INTEGER IMSSM,IHIGGS,PARITY,LRSUSY
261 COMMON/HW6202/VIPWID(3),DXRCYL,DXZMAX,DXRSPH,WZRFR,FIX4JT,
262 & IMSSM,IHIGGS,PARITY,LRSUSY
264 C New parameters for version 6.203
265 DOUBLE PRECISION ABWGT,ABWSUM,AVABW
266 INTEGER NNEGWT,NNEGEV
268 COMMON/HW6203/ABWGT,ABWSUM,AVABW,NNEGWT,NNEGEV,NEGWTS
270 C New parameters for version 6.3
271 INTEGER IMAXCH,IMAXOP
272 PARAMETER (IMAXCH=20,IMAXOP=40)
273 DOUBLE PRECISION MJJMIN,CHNPRB(IMAXCH)
275 LOGICAL OPTM,CHON(IMAXCH)
276 COMMON/HW6300/MJJMIN,CHNPRB,IOPSTP,IOPSH,OPTM,CHON
277 C New PDF's for version 6.3
278 INTEGER NXMRS,NQMRS,NPMRS
279 PARAMETER(NXMRS=49,NQMRS=37,NPMRS=8)
280 DOUBLE PRECISION FMRS(3,NPMRS,NXMRS,NQMRS+1)
282 C Circe interface for version 6.3
283 INTEGER CIRCOP,CIRCAC,CIRCVR,CIRCRV,CIRCCH
284 COMMON /HWCIRC/CIRCOP,CIRCAC,CIRCVR,CIRCRV,CIRCCH
285 C New parameters and commons for spin correlations
286 C--constants for the arrays
287 INTEGER NMXSPN,NCFMAX
288 PARAMETER(NMXSPN=50,NCFMAX=3)
289 INTEGER NMODE2,NMODE3,NDIAGR,NMODEB,NMODE4
290 PARAMETER(NMODE2=500,NMODE3=500,NDIAGR=8,NMODEB=50,NMODE4=4)
291 C--common block for X --> X gauge boson
292 DOUBLE PRECISION ABMODE(2,NMODEB),BBMODE(2,12,NMODEB),
293 & PBMODE(12,NMODEB),WTBMAX(12,NMODEB)
294 INTEGER IDBPRT(NMODEB),IBMODE(NMODEB),IBDRTP(NMODEB),NBMODE
295 COMMON /HWDSPB/ABMODE,BBMODE,PBMODE,WTBMAX,IDBPRT,IBDRTP,IBMODE,
297 C--common block for two body decays
298 DOUBLE PRECISION A2MODE(2,NMODE2),P2MODE(NMODE2),WT2MAX(NMODE2)
299 INTEGER ID2PRT(NMODE2),I2DRTP(NMODE2),N2MODE
300 COMMON /HWDSP2/A2MODE,P2MODE,WT2MAX,ID2PRT,I2DRTP,N2MODE
301 C--common block for three body decays
302 DOUBLE PRECISION A3MODE(2,NDIAGR,NMODE3),B3MODE(2,NDIAGR,NMODE3),
303 & P3MODE(NMODE3),WT3MAX(NMODE3),SPN3CF(NCFMAX,NCFMAX,NMODE3)
304 INTEGER ID3PRT(NMODE3),I3MODE(NDIAGR,NMODE3),
305 & I3DRTP(NDIAGR,NMODE3),N3MODE,NDI3BY(NMODE3),N3NCFL(NMODE3),
306 & I3DRCF(NDIAGR,NMODE3)
307 COMMON /HWDSP3/A3MODE,B3MODE,P3MODE,WT3MAX,SPN3CF,ID3PRT,I3MODE,
308 & I3DRTP,N3MODE,NDI3BY,N3NCFL,I3DRCF
309 C--common block for four body decays
310 DOUBLE PRECISION A4MODE(2,12,NMODE4),B4MODE(2,12,NMODE4),
311 & P4MODE(12,12,NMODE4),WT4MAX(12,12,NMODE4)
312 INTEGER ID4PRT(NMODE4),I4MODE(2,NMODE4),N4MODE
313 COMMON /HWDSP4/A4MODE,B4MODE,P4MODE,WT4MAX,ID4PRT,I4MODE,N4MODE
314 C--common block for spin correlations in event
315 INTEGER NDECSY,NSEARCH,LRDEC,LWDEC
316 LOGICAL SYSPIN,THREEB,FOURB
318 COMMON /HWDSPN/NDECSY,NSEARCH,LRDEC,LWDEC,SYSPIN,THREEB,
321 INTEGER IDSPN(NMXSPN),JMOSPN(NMXSPN),JDASPN(2,NMXSPN),NSPN,
322 & ISNHEP(NMXHEP),NSNTRY,NCFL(NMXSPN),SPCOPT
323 DOUBLE COMPLEX MESPN(2,2,2,2,NCFMAX,NMXSPN),RHOSPN(2,2,NMXSPN)
324 DOUBLE PRECISION SPNCFC(NCFMAX,NCFMAX,NMXSPN)
325 LOGICAL DECSPN(NMXSPN)
326 COMMON /HWSPIN/MESPN,RHOSPN,SPNCFC,IDSPN,JMOSPN,JDASPN,
327 & NSPN,ISNHEP,NSNTRY,DECSPN,NCFL,SPCOPT
328 INTEGER JAK1,JAK2,ITDKRC,IFPHOT
329 COMMON /HWSTAU/ JAK1,JAK2,ITDKRC,IFPHOT
331 C--common block for Les Houches interface to store information we need
334 PARAMETER (MAXHRP=100)
335 DOUBLE PRECISION LHWGT(MAXHRP),LHWGTS(MAXHRP),LHMXSM,
336 & LHXSCT(MAXHRP),LHXERR(MAXHRP),LHXMAX(MAXHRP)
337 INTEGER LHIWGT(MAXHRP),ITYPLH,LHNEVT(MAXHRP)
338 LOGICAL LHSOFT,LHGLSF
339 COMMON /HWGUPR/LHWGT,LHWGTS,LHXSCT,LHXERR,LHXMAX,LHMXSM,LHIWGT,
340 & LHNEVT,ITYPLH,LHSOFT,LHGLSF
342 C--common block for HERWIG6.5
345 COMMON /HW6500/ PRESPL
347 COMMON /HW6504/ITOPRD
348 C--the only change for 6505 was to increase MODMAX from 5 to 50
349 DOUBLE PRECISION PDFX0,PDFPOW
350 COMMON /HW6506/PDFX0,PDFPOW