C ****COMMON BLOCK FILE FOR HERWIG VERSION 6.1**** C C ALTERATIONS: Layout completely overhauled for 5.9 C C C New common blocks added for version 6.1: C HWCLUS,HWSUSY,HWRPAR,HWMINB C C New variables added for version 6.1: C OMHMIX,ET2MIX,PH3MIX,IOP4JT,NPRFMT, C PRNDEF,PRNTEX,PRNWEB,EFFMIN,GCUTME, C IOP4JT,NPRFMT see HWPRAM C Y4JT,DURHAM see HWHARD C QORQQB,QBORQQ see HWPROP C NRECO see HWUCLU C TXNAME see HWUNAM C PPCL,NCL,IDCL see HWCLUS C TANB,ALPHAH,COSBPA,SINBPA,COSBMA, C SINBMA,COSA,SINA,COSB,SINB,COTB, C ZMIXSS,ZMXNSS,ZSGNSS,LFCH,RFCH, C SLFCH,SRFCH,WMXUSS,WMXVSS,WSGNSS, C QMIXSS,LMIXSS,THETAT,THETAB,THETAL, C ATSS,ABSS,ALSS,MUSS,FACTSS,GHWWSS, C GHZZSS,GHDDSS,GHUUSS,GHWHSS,GHSQSS, C XLMNSS,RMMNSS,IMSSM,SENHNC, C SSPARITY,SUSYIN see HWSUSY C LAMDA1,LAMDA2,LAMDA3,HRDCOL,RPARTY, C COLUPD see HWRPAR C PMBN1,PMBN2,PMBN3,PMBK1,PMBK2, C PMBM1,PMBM2,PMBP1,PMBP2,PMBP3 see HWMINB C C New parameters added for version 6.1: C NMXCL C C Parameter NMXRES raised to 500 C C Scalar variables changed to arrays of size 2: C CLSMR,PSPLT,CLDIR see HWPRAM C IMPLICIT NONE DOUBLE PRECISION ZERO,ONE,TWO,THREE,FOUR,HALF PARAMETER (ZERO =0.D0, ONE =1.D0, TWO =2.D0, & THREE=3.D0, FOUR=4.D0, HALF=0.5D0) C DOUBLE PRECISION & ACCUR,AFCH,ALPFAC,ALPHEM,ANOMSC,ASFIXD,AVWGT,B1LIM,BETAF,BRFRAC, & BRHIG,BTCLM,CAFAC,CFFAC,CLDKWT,CLMAX,CLPOW,CLQ,CLSMR,CMMOM,COSS, & COSTH,CSPEED,CTHRPW,CTMAX,DECPAR,DECWT,DISF,DKLTM,EBEAM1,EBEAM2, & EMLST,EMMAX,EMMIN,EMPOW,EMSCA,ENHANC,ENSOF,EPOLN,ETAMIX,EVWGT, & EXAG,F0MIX,F1MIX,F2MIX,GAMH,GAMMAX,GAMW,GAMWT,GAMZ,GAMZP,GCOEF, & GEV2NB,GEV2MM,GPOLN,H1MIX,HBAR,HARDST,OMEGA0,PBEAM1,PBEAM2,PDIQK, & PGSMX,PGSPL,PHEP,PHIMIX,PHIPAR,PHOMAS,PIFAC,PLTCUT,PPAR,PPOLN, & PRECO,PRSOF,PSPLT,PTINT,PTMAX,PTMIN,PTPOW,PTRMS,PXRMS,PWT,Q2MAX, & Q2MIN,Q2POW,Q2WWMN,Q2WWMX,QCDL3,QCDL5,QCDLAM,QDIQK,QEV,QFCH,QG, & QLIM,QSPAC,QV,QWT,REPWT,RESN,RHOHEP,RHOPAR,RLTIM,RMASS,RMIN, & RSPIN,SCABI,SINS,SNGWT,SWEIN,SWTEF,SUD,THMAX,TLOUT,TMTOP,TMNISR, & TQWT,VCKM,VFCH,VGCUT,VHEP,VMIN2,VPAR,VPCUT,VQCUT,VTXPIP,VTXQDK, & WBIGST,WGTMAX,WGTSUM,WHMIN,WSQSUM,XFACT,XLMIN,XMIX,XMRCT,XX, & XXMIN,YBMAX,YBMIN,YJMAX,YJMIN,YMIX,YMRCT,YWWMAX,YWWMIN,ZBINM, & ZJMAX,ZMXISR,Y4JT,EFFMIN,PPCL, & TANB,ALPHAH,COSBPA,SINBPA,COSBMA,SINBMA,COSA,SINA,COSB,SINB,COTB, & ZMIXSS,ZMXNSS,ZSGNSS,LFCH,RFCH,SLFCH,SRFCH, WMXUSS,WMXVSS,WSGNSS, & QMIXSS,LMIXSS,THETAT,THETAB,THETAL,ATSS,ABSS,ALSS,MUSS,FACTSS, & GHWWSS,GHZZSS,GHDDSS,GHUUSS,GHWHSS,GHSQSS, & XLMNSS,RMMNSS,IMSSM,SENHNC,SSPARITY,LAMDA1,LAMDA2,LAMDA3, & PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,PMBM1,PMBM2,PMBP1,PMBP2,PMBP3, & OMHMIX,ET2MIX,PH3MIX,GCUTME C INTEGER & CLDIR,IAPHIG,IBRN,IBSH,ICHRG,ICO,IDCMF,IDHEP,IDHW,IDK,IDKPRD,IDN, & IDPAR,IDPDG,IERROR,IFLAV,IFLMAX,IFLMIN,IHPRO,IMQDK,INHAD,INTER, & IOPDKL,IOPHIG,IOPREM,IPART1,IPART2,IPRINT,IPRO,IPROC,ISLENT, & ISPAC,ISTAT,ISTHEP,ISTPAR,JCOPAR,JDAHEP,JDAPAR,JMOHEP,JMOPAR, & JNHAD,LNEXT,LOCN,LOCQ,LRSUD,LSTRT,LWEVT,LWSUD,MAPQ,MAXER,MAXEV, & MAXFL,MAXPR,MODBOS,MODMAX,MODPDF,NBTRY,NCLDK,NCOLO,NCTRY,NDKYS, & NDTRY,NETRY,NEVHEP,NEVPAR,NFLAV,NGSPL,NHEP,NME,NMODES,NMXCDK, & NMXDKS,NMXHEP,NMXJET,NMXMOD,NMXPAR,NMXQDK,NMXRES,NMXSUD,NPAR, & NPRODS,NQDK,NQEV,NRES,NRN,NSPAC,NSTRU,NSTRY,NSUD,NUMER,NUMERU, & NWGTS,NZBIN,SUDORD,IOP4JT,HRDCOL,NMXCL,NCL,IDCL,NPRFMT,NRECO C LOGICAL & AZSOFT,AZSPIN,BGSHAT,BREIT,CLRECO,COLISR,DKPSET,FROST,FSTEVT, & FSTWGT,GENEV,GENSOF,HARDME,HVFCEN,MAXDKL,MIXING,NOSPAC,NOWGT, & PRNDEC,PIPSMR,PRVTX,RSTAB,SOFTME,TMPAR,TPOL,USECMF,VTOCDK,VTORDK, & ZPRIME,RPARTY,COLUPD,PRNDEF,PRNTEX,PRNWEB,DURHAM,SUSYIN, & QORQQB,QBORQQ C CHARACTER*4 & BDECAY CHARACTER*8 & PART1,PART2,RNAME CHARACTER*20 & AUTPDF CHARACTER*37 & TXNAME C C New standard event common PARAMETER (NMXHEP=2000) COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP), & JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP) C C Beams, process and number of events COMMON/HWBEAM/IPART1,IPART2 COMMON/HWBMCH/PART1,PART2 COMMON/HWPROC/EBEAM1,EBEAM2,PBEAM1,PBEAM2,IPROC,MAXEV C C Basic parameters (and quantities derived from them) COMMON/HWPRAM/AFCH(16,2),ALPHEM,B1LIM,BETAF,BTCLM,CAFAC,CFFAC, & CLMAX,CLPOW,CLSMR(2),CSPEED,ENSOF,ETAMIX,F0MIX,F1MIX,F2MIX,GAMH, & GAMW,GAMZ,GAMZP,GEV2NB,H1MIX,PDIQK,PGSMX,PGSPL(4),PHIMIX,PIFAC, & PRSOF,PSPLT(2),PTRMS,PXRMS,QCDL3,QCDL5,QCDLAM,QDIQK,QFCH(16),QG, & QSPAC,QV,SCABI,SWEIN,TMTOP,VFCH(16,2),VCKM(3,3),VGCUT,VQCUT, & VPCUT,ZBINM,EFFMIN,OMHMIX,ET2MIX,PH3MIX,GCUTME, & IOPREM,IPRINT,ISPAC,LRSUD,LWSUD,MODPDF(2),NBTRY,NCOLO,NCTRY, & NDTRY,NETRY,NFLAV,NGSPL,NSTRU,NSTRY,NZBIN,IOP4JT(2),NPRFMT, & AZSOFT,AZSPIN,CLDIR(2),HARDME,NOSPAC,PRNDEC,PRVTX,SOFTME,ZPRIME, & PRNDEF,PRNTEX,PRNWEB C COMMON/HWPRCH/AUTPDF(2),BDECAY C C Parton shower common (same format as /HEPEVT/) PARAMETER (NMXPAR=500) COMMON/HWPART/NEVPAR,NPAR,ISTPAR(NMXPAR),IDPAR(NMXPAR), & JMOPAR(2,NMXPAR),JDAPAR(2,NMXPAR),PPAR(5,NMXPAR),VPAR(4,NMXPAR) C C Parton polarization common COMMON/HWPARP/DECPAR(2,NMXPAR),PHIPAR(2,NMXPAR),RHOPAR(2,NMXPAR), & TMPAR(NMXPAR) C C Electroweak boson common PARAMETER (MODMAX=5) COMMON/HWBOSC/ALPFAC,BRHIG(12),ENHANC(12),GAMMAX,RHOHEP(3,NMXHEP), & IOPHIG,MODBOS(MODMAX) C C Parton colour common COMMON/HWPARC/JCOPAR(4,NMXPAR) C C other HERWIG branching, event and hard subprocess common blocks COMMON/HWBRCH/ANOMSC(2,2),HARDST,PTINT(3,2),XFACT,INHAD,JNHAD, & NSPAC(7),ISLENT,BREIT,FROST,USECMF C COMMON/HWEVNT/AVWGT,EVWGT,GAMWT,TLOUT,WBIGST,WGTMAX,WGTSUM,WSQSUM, & IDHW(NMXHEP),IERROR,ISTAT,LWEVT,MAXER,MAXPR,NOWGT,NRN(2),NUMER, & NUMERU,NWGTS,GENSOF C COMMON/HWHARD/ASFIXD,CLQ(7,6),COSS,COSTH,CTMAX,DISF(13,2),EMLST, & EMMAX,EMMIN,EMPOW,EMSCA,EPOLN(3),GCOEF(7),GPOLN,OMEGA0,PHOMAS, & PPOLN(3),PTMAX,PTMIN,PTPOW,Q2MAX,Q2MIN,Q2POW,Q2WWMN,Q2WWMX,QLIM, & SINS,THMAX,Y4JT,TMNISR,TQWT,XX(2),XLMIN,XXMIN,YBMAX,YBMIN,YJMAX, & YJMIN,YWWMAX,YWWMIN,WHMIN,ZJMAX,ZMXISR,IAPHIG,IBRN(2),IBSH, & ICO(10),IDCMF,IDN(10),IFLMAX,IFLMIN,IHPRO,IPRO,MAPQ(6),MAXFL, & BGSHAT,COLISR,FSTEVT,FSTWGT,GENEV,HVFCEN,TPOL,DURHAM C C Arrays for particle properties (NMXRES = max no of particles defined) PARAMETER(NMXRES=500) COMMON/HWPROP/RLTIM(0:NMXRES),RMASS(0:NMXRES),RSPIN(0:NMXRES), & ICHRG(0:NMXRES),IDPDG(0:NMXRES),IFLAV(0:NMXRES),NRES, & VTOCDK(0:NMXRES),VTORDK(0:NMXRES), & QORQQB(0:NMXRES),QBORQQ(0:NMXRES) C COMMON/HWUNAM/RNAME(0:NMXRES),TXNAME(2,0:NMXRES) C C Arrays for particle decays (NMXDKS = max total no of decays, C NMXMOD = max no of modes for a particle) PARAMETER(NMXDKS=4000,NMXMOD=200) COMMON/HWUPDT/BRFRAC(NMXDKS),CMMOM(NMXDKS),DKLTM(NMXRES), & IDK(NMXDKS),IDKPRD(5,NMXDKS),LNEXT(NMXDKS),LSTRT(NMXRES),NDKYS, & NME(NMXDKS),NMODES(NMXRES),NPRODS(NMXDKS),DKPSET,RSTAB(0:NMXRES) C C Weights used in cluster decays COMMON/HWUWTS/REPWT(0:3,0:4,0:4),SNGWT,DECWT,QWT(3),PWT(12), & SWTEF(NMXRES) C C Parameters for cluster decays (NMXCDK = max total no of cluster C decay channels) PARAMETER(NMXCDK=4000) COMMON/HWUCLU/CLDKWT(NMXCDK),CTHRPW(12,12),PRECO,RESN(12,12), & RMIN(12,12),LOCN(12,12),NCLDK(NMXCDK),NRECO,CLRECO C C Variables controling mixing and vertex information COMMON/HWDIST/EXAG,GEV2MM,HBAR,PLTCUT,VMIN2,VTXPIP(4),XMIX(2), & XMRCT(2),YMIX(2),YMRCT(2),IOPDKL,MAXDKL,MIXING,PIPSMR C C Arrays for temporarily storing heavy-b,c-hadrons decaying partonicaly C (NMXQDK = max no such decays in an event) PARAMETER (NMXQDK=20) COMMON/HWQDKS/VTXQDK(4,NMXQDK),IMQDK(NMXQDK),LOCQ(NMXQDK),NQDK C C Parameters for Sudakov form factors C (NMXSUD= max no of entries in lookup table) PARAMETER (NMXSUD=1024) COMMON/HWUSUD/ACCUR,QEV(NMXSUD,6),SUD(NMXSUD,6),INTER,NQEV,NSUD, & SUDORD C PARAMETER (NMXJET=200) C C SUSY parameters COMMON/HWSUSY/ & TANB,ALPHAH,COSBPA,SINBPA,COSBMA,SINBMA,COSA,SINA,COSB,SINB,COTB, & ZMIXSS(4,4),ZMXNSS(4,4),ZSGNSS(4), LFCH(16),RFCH(16), & SLFCH(16,4),SRFCH(16,4), WMXUSS(2,2),WMXVSS(2,2), WSGNSS(2), & QMIXSS(6,2,2),LMIXSS(6,2,2), & THETAT,THETAB,THETAL,ATSS,ABSS,ALSS,MUSS,FACTSS, & GHWWSS(3),GHZZSS(3),GHDDSS(4),GHUUSS(4),GHWHSS(3), & GHSQSS(4,6,2,2),XLMNSS,RMMNSS,IMSSM,SENHNC(24),SSPARITY,SUSYIN C C R-Parity violating parameters and colours COMMON /HWRPAR/ LAMDA1(3,3,3),LAMDA2(3,3,3), & LAMDA3(3,3,3),HRDCOL(2,5),RPARTY,COLUPD C C Parameters for minimum bias/soft underlying event COMMON/HWMINB/ & PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,PMBM1,PMBM2,PMBP1,PMBP2,PMBP3 C C Cluster common used by soft event routines PARAMETER (NMXCL=500) COMMON/HWCLUS/PPCL(5,NMXCL),IDCL(NMXCL),NCL