HERWIG/HERWIG61.INC

   1 C          ****COMMON BLOCK FILE FOR HERWIG VERSION 6.1****
   2 C
   3 C ALTERATIONS: Layout completely overhauled for 5.9
   4 C
   5 C
   6 C New common blocks added for version 6.1:
   7 C              HWCLUS,HWSUSY,HWRPAR,HWMINB
   8 C
   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
  15 C              NRECO                               see HWUCLU
  16 C              TXNAME                              see HWUNAM
  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,IMSSM,SENHNC,
  26 C              SSPARITY,SUSYIN                     see HWSUSY
  27 C              LAMDA1,LAMDA2,LAMDA3,HRDCOL,RPARTY,
  28 C              COLUPD                              see HWRPAR
  29 C              PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,
  30 C              PMBM1,PMBM2,PMBP1,PMBP2,PMBP3       see HWMINB
  31 C
  32 C New parameters added for version 6.1:
  33 C              NMXCL
  34 C
  35 C Parameter NMXRES raised to 500
  36 C
  37 C Scalar variables changed to arrays of size 2:
  38 C              CLSMR,PSPLT,CLDIR                   see HWPRAM
  39 C
  40       IMPLICIT NONE
  41       DOUBLE PRECISION ZERO,ONE,TWO,THREE,FOUR,HALF
  42       PARAMETER (ZERO =0.D0, ONE =1.D0, TWO =2.D0,
  43      &           THREE=3.D0, FOUR=4.D0, HALF=0.5D0)
  44 C
  45       DOUBLE PRECISION
  46      & ACCUR,AFCH,ALPFAC,ALPHEM,ANOMSC,ASFIXD,AVWGT,B1LIM,BETAF,BRFRAC,
  47      & BRHIG,BTCLM,CAFAC,CFFAC,CLDKWT,CLMAX,CLPOW,CLQ,CLSMR,CMMOM,COSS,
  48      & COSTH,CSPEED,CTHRPW,CTMAX,DECPAR,DECWT,DISF,DKLTM,EBEAM1,EBEAM2,
  49      & EMLST,EMMAX,EMMIN,EMPOW,EMSCA,ENHANC,ENSOF,EPOLN,ETAMIX,EVWGT,
  50      & EXAG,F0MIX,F1MIX,F2MIX,GAMH,GAMMAX,GAMW,GAMWT,GAMZ,GAMZP,GCOEF,
  51      & GEV2NB,GEV2MM,GPOLN,H1MIX,HBAR,HARDST,OMEGA0,PBEAM1,PBEAM2,PDIQK,
  52      & PGSMX,PGSPL,PHEP,PHIMIX,PHIPAR,PHOMAS,PIFAC,PLTCUT,PPAR,PPOLN,
  53      & PRECO,PRSOF,PSPLT,PTINT,PTMAX,PTMIN,PTPOW,PTRMS,PXRMS,PWT,Q2MAX,
  54      & Q2MIN,Q2POW,Q2WWMN,Q2WWMX,QCDL3,QCDL5,QCDLAM,QDIQK,QEV,QFCH,QG,
  55      & QLIM,QSPAC,QV,QWT,REPWT,RESN,RHOHEP,RHOPAR,RLTIM,RMASS,RMIN,
  56      & RSPIN,SCABI,SINS,SNGWT,SWEIN,SWTEF,SUD,THMAX,TLOUT,TMTOP,TMNISR,
  57      & TQWT,VCKM,VFCH,VGCUT,VHEP,VMIN2,VPAR,VPCUT,VQCUT,VTXPIP,VTXQDK,
  58      & WBIGST,WGTMAX,WGTSUM,WHMIN,WSQSUM,XFACT,XLMIN,XMIX,XMRCT,XX,
  59      & XXMIN,YBMAX,YBMIN,YJMAX,YJMIN,YMIX,YMRCT,YWWMAX,YWWMIN,ZBINM,
  60      & ZJMAX,ZMXISR,Y4JT,EFFMIN,PPCL,
  61      & TANB,ALPHAH,COSBPA,SINBPA,COSBMA,SINBMA,COSA,SINA,COSB,SINB,COTB,
  62      & ZMIXSS,ZMXNSS,ZSGNSS,LFCH,RFCH,SLFCH,SRFCH, WMXUSS,WMXVSS,WSGNSS,
  63      & QMIXSS,LMIXSS,THETAT,THETAB,THETAL,ATSS,ABSS,ALSS,MUSS,FACTSS,
  64      & GHWWSS,GHZZSS,GHDDSS,GHUUSS,GHWHSS,GHSQSS,
  65      & XLMNSS,RMMNSS,IMSSM,SENHNC,SSPARITY,LAMDA1,LAMDA2,LAMDA3,
  66      & PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,PMBM1,PMBM2,PMBP1,PMBP2,PMBP3,
  67      & OMHMIX,ET2MIX,PH3MIX,GCUTME
  68 C
  69       INTEGER
  70      & CLDIR,IAPHIG,IBRN,IBSH,ICHRG,ICO,IDCMF,IDHEP,IDHW,IDK,IDKPRD,IDN,
  71      & IDPAR,IDPDG,IERROR,IFLAV,IFLMAX,IFLMIN,IHPRO,IMQDK,INHAD,INTER,
  72      & IOPDKL,IOPHIG,IOPREM,IPART1,IPART2,IPRINT,IPRO,IPROC,ISLENT,
  73      & ISPAC,ISTAT,ISTHEP,ISTPAR,JCOPAR,JDAHEP,JDAPAR,JMOHEP,JMOPAR,
  74      & JNHAD,LNEXT,LOCN,LOCQ,LRSUD,LSTRT,LWEVT,LWSUD,MAPQ,MAXER,MAXEV,
  75      & MAXFL,MAXPR,MODBOS,MODMAX,MODPDF,NBTRY,NCLDK,NCOLO,NCTRY,NDKYS,
  76      & NDTRY,NETRY,NEVHEP,NEVPAR,NFLAV,NGSPL,NHEP,NME,NMODES,NMXCDK,
  77      & NMXDKS,NMXHEP,NMXJET,NMXMOD,NMXPAR,NMXQDK,NMXRES,NMXSUD,NPAR,
  78      & NPRODS,NQDK,NQEV,NRES,NRN,NSPAC,NSTRU,NSTRY,NSUD,NUMER,NUMERU,
  79      & NWGTS,NZBIN,SUDORD,IOP4JT,HRDCOL,NMXCL,NCL,IDCL,NPRFMT,NRECO
  80 C
  81       LOGICAL
  82      & AZSOFT,AZSPIN,BGSHAT,BREIT,CLRECO,COLISR,DKPSET,FROST,FSTEVT,
  83      & FSTWGT,GENEV,GENSOF,HARDME,HVFCEN,MAXDKL,MIXING,NOSPAC,NOWGT,
  84      & PRNDEC,PIPSMR,PRVTX,RSTAB,SOFTME,TMPAR,TPOL,USECMF,VTOCDK,VTORDK,
  85      & ZPRIME,RPARTY,COLUPD,PRNDEF,PRNTEX,PRNWEB,DURHAM,SUSYIN,
  86      & QORQQB,QBORQQ
  87 C
  88       CHARACTER*4
  89      & BDECAY
  90       CHARACTER*8
  91      & PART1,PART2,RNAME
  92       CHARACTER*20
  93      & AUTPDF
  94       CHARACTER*37
  95      & TXNAME
  96 C
  97 C New standard event common
  98       PARAMETER (NMXHEP=2000)
  99       COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
 100      & JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
 101 C
 102 C Beams, process and number of events
 103       COMMON/HWBEAM/IPART1,IPART2
 104       COMMON/HWBMCH/PART1,PART2
 105       COMMON/HWPROC/EBEAM1,EBEAM2,PBEAM1,PBEAM2,IPROC,MAXEV
 106 C
 107 C Basic parameters (and quantities derived from them)
 108       COMMON/HWPRAM/AFCH(16,2),ALPHEM,B1LIM,BETAF,BTCLM,CAFAC,CFFAC,
 109      & CLMAX,CLPOW,CLSMR(2),CSPEED,ENSOF,ETAMIX,F0MIX,F1MIX,F2MIX,GAMH,
 110      & GAMW,GAMZ,GAMZP,GEV2NB,H1MIX,PDIQK,PGSMX,PGSPL(4),PHIMIX,PIFAC,
 111      & PRSOF,PSPLT(2),PTRMS,PXRMS,QCDL3,QCDL5,QCDLAM,QDIQK,QFCH(16),QG,
 112      & QSPAC,QV,SCABI,SWEIN,TMTOP,VFCH(16,2),VCKM(3,3),VGCUT,VQCUT,
 113      & VPCUT,ZBINM,EFFMIN,OMHMIX,ET2MIX,PH3MIX,GCUTME,
 114      & IOPREM,IPRINT,ISPAC,LRSUD,LWSUD,MODPDF(2),NBTRY,NCOLO,NCTRY,
 115      & NDTRY,NETRY,NFLAV,NGSPL,NSTRU,NSTRY,NZBIN,IOP4JT(2),NPRFMT,
 116      & AZSOFT,AZSPIN,CLDIR(2),HARDME,NOSPAC,PRNDEC,PRVTX,SOFTME,ZPRIME,
 117      & PRNDEF,PRNTEX,PRNWEB
 118 C
 119       COMMON/HWPRCH/AUTPDF(2),BDECAY
 120 C
 121 C Parton shower common (same format as /HEPEVT/)
 122       PARAMETER (NMXPAR=500)
 123       COMMON/HWPART/NEVPAR,NPAR,ISTPAR(NMXPAR),IDPAR(NMXPAR),
 124      & JMOPAR(2,NMXPAR),JDAPAR(2,NMXPAR),PPAR(5,NMXPAR),VPAR(4,NMXPAR)
 125 C
 126 C Parton polarization common
 127       COMMON/HWPARP/DECPAR(2,NMXPAR),PHIPAR(2,NMXPAR),RHOPAR(2,NMXPAR),
 128      & TMPAR(NMXPAR)
 129 C
 130 C Electroweak boson common
 131       PARAMETER (MODMAX=5)
 132       COMMON/HWBOSC/ALPFAC,BRHIG(12),ENHANC(12),GAMMAX,RHOHEP(3,NMXHEP),
 133      & IOPHIG,MODBOS(MODMAX)
 134 C
 135 C Parton colour common
 136       COMMON/HWPARC/JCOPAR(4,NMXPAR)
 137 C
 138 C other HERWIG branching, event and hard subprocess common blocks
 139       COMMON/HWBRCH/ANOMSC(2,2),HARDST,PTINT(3,2),XFACT,INHAD,JNHAD,
 140      & NSPAC(7),ISLENT,BREIT,FROST,USECMF
 141 C
 142       COMMON/HWEVNT/AVWGT,EVWGT,GAMWT,TLOUT,WBIGST,WGTMAX,WGTSUM,WSQSUM,
 143      & IDHW(NMXHEP),IERROR,ISTAT,LWEVT,MAXER,MAXPR,NOWGT,NRN(2),NUMER,
 144      & NUMERU,NWGTS,GENSOF
 145 C
 146       COMMON/HWHARD/ASFIXD,CLQ(7,6),COSS,COSTH,CTMAX,DISF(13,2),EMLST,
 147      & EMMAX,EMMIN,EMPOW,EMSCA,EPOLN(3),GCOEF(7),GPOLN,OMEGA0,PHOMAS,
 148      & PPOLN(3),PTMAX,PTMIN,PTPOW,Q2MAX,Q2MIN,Q2POW,Q2WWMN,Q2WWMX,QLIM,
 149      & SINS,THMAX,Y4JT,TMNISR,TQWT,XX(2),XLMIN,XXMIN,YBMAX,YBMIN,YJMAX,
 150      & YJMIN,YWWMAX,YWWMIN,WHMIN,ZJMAX,ZMXISR,IAPHIG,IBRN(2),IBSH,
 151      & ICO(10),IDCMF,IDN(10),IFLMAX,IFLMIN,IHPRO,IPRO,MAPQ(6),MAXFL,
 152      & BGSHAT,COLISR,FSTEVT,FSTWGT,GENEV,HVFCEN,TPOL,DURHAM
 153 C
 154 C Arrays for particle properties (NMXRES = max no of particles defined)
 155       PARAMETER(NMXRES=500)
 156       COMMON/HWPROP/RLTIM(0:NMXRES),RMASS(0:NMXRES),RSPIN(0:NMXRES),
 157      & ICHRG(0:NMXRES),IDPDG(0:NMXRES),IFLAV(0:NMXRES),NRES,
 158      & VTOCDK(0:NMXRES),VTORDK(0:NMXRES),
 159      & QORQQB(0:NMXRES),QBORQQ(0:NMXRES)
 160 C
 161       COMMON/HWUNAM/RNAME(0:NMXRES),TXNAME(2,0:NMXRES)
 162 C
 163 C Arrays for particle decays (NMXDKS = max total no of decays,
 164 C                             NMXMOD = max no of modes for a particle)
 165       PARAMETER(NMXDKS=4000,NMXMOD=200)
 166       COMMON/HWUPDT/BRFRAC(NMXDKS),CMMOM(NMXDKS),DKLTM(NMXRES),
 167      & IDK(NMXDKS),IDKPRD(5,NMXDKS),LNEXT(NMXDKS),LSTRT(NMXRES),NDKYS,
 168      & NME(NMXDKS),NMODES(NMXRES),NPRODS(NMXDKS),DKPSET,RSTAB(0:NMXRES)
 169 C
 170 C Weights used in cluster decays
 171       COMMON/HWUWTS/REPWT(0:3,0:4,0:4),SNGWT,DECWT,QWT(3),PWT(12),
 172      & SWTEF(NMXRES)
 173 C
 174 C Parameters for cluster decays (NMXCDK = max total no of cluster
 175 C                                         decay channels)
 176       PARAMETER(NMXCDK=4000)
 177       COMMON/HWUCLU/CLDKWT(NMXCDK),CTHRPW(12,12),PRECO,RESN(12,12),
 178      & RMIN(12,12),LOCN(12,12),NCLDK(NMXCDK),NRECO,CLRECO
 179 C
 180 C Variables controling mixing and vertex information
 181       COMMON/HWDIST/EXAG,GEV2MM,HBAR,PLTCUT,VMIN2,VTXPIP(4),XMIX(2),
 182      & XMRCT(2),YMIX(2),YMRCT(2),IOPDKL,MAXDKL,MIXING,PIPSMR
 183 C
 184 C Arrays for temporarily storing heavy-b,c-hadrons decaying partonicaly
 185 C (NMXQDK = max no such decays in an event)
 186       PARAMETER (NMXQDK=20)
 187       COMMON/HWQDKS/VTXQDK(4,NMXQDK),IMQDK(NMXQDK),LOCQ(NMXQDK),NQDK
 188 C
 189 C Parameters for Sudakov form factors
 190 C (NMXSUD= max no of entries in lookup table)
 191       PARAMETER (NMXSUD=1024)
 192       COMMON/HWUSUD/ACCUR,QEV(NMXSUD,6),SUD(NMXSUD,6),INTER,NQEV,NSUD,
 193      & SUDORD
 194 C
 195       PARAMETER (NMXJET=200)
 196 C
 197 C SUSY parameters
 198       COMMON/HWSUSY/
 199      & TANB,ALPHAH,COSBPA,SINBPA,COSBMA,SINBMA,COSA,SINA,COSB,SINB,COTB,
 200      & ZMIXSS(4,4),ZMXNSS(4,4),ZSGNSS(4), LFCH(16),RFCH(16),
 201      & SLFCH(16,4),SRFCH(16,4), WMXUSS(2,2),WMXVSS(2,2), WSGNSS(2),
 202      & QMIXSS(6,2,2),LMIXSS(6,2,2),
 203      & THETAT,THETAB,THETAL,ATSS,ABSS,ALSS,MUSS,FACTSS,
 204      & GHWWSS(3),GHZZSS(3),GHDDSS(4),GHUUSS(4),GHWHSS(3),
 205      & GHSQSS(4,6,2,2),XLMNSS,RMMNSS,IMSSM,SENHNC(24),SSPARITY,SUSYIN
 206 C
 207 C R-Parity violating parameters and colours
 208       COMMON /HWRPAR/ LAMDA1(3,3,3),LAMDA2(3,3,3),
 209      &                LAMDA3(3,3,3),HRDCOL(2,5),RPARTY,COLUPD
 210 C
 211 C Parameters for minimum bias/soft underlying event
 212       COMMON/HWMINB/
 213      & PMBN1,PMBN2,PMBN3,PMBK1,PMBK2,PMBM1,PMBM2,PMBP1,PMBP2,PMBP3
 214 C
 215 C Cluster common used by soft event routines
 216       PARAMETER (NMXCL=500)
 217       COMMON/HWCLUS/PPCL(5,NMXCL),IDCL(NMXCL),NCL