C
C written by Martin W. Gruenewald (91/01/28)
C ----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
LOGICAL QEDRAD
COMMON /PHORAD/ QEDRAD(NMXHEP)
INTEGER N,IHEP,IST,ID,JMO1,JMO2,JDA1,JDA2,I,IP
C Azimuthal correlation functions for Collins' algorithm,
C see I.G.Knowles, Comp. Phys. Comm. 58 (90) 271 for notation.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION Z1,Z2,DOT12,DOT23,DOT31,TR,FN(7),VEC1(2),VEC2(2),
& VEC3(2),VEC(2)
INTEGER IPAR,JPAR
C OF DECAYS)
C NEW VARAIBLE BACK TO ALLOW CODE TO SEARCH DOWN CHAIN
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,IST,ID,JC,KC,JD,JHEP,LHEP,ID2,NTRY,KHEP
LOGICAL BACK
IF (IERROR.NE.0) RETURN
C********MODIFIED 13/11/00 BY BRW TO ALLOW MULTIPLE APPLICATION IN
C*******SAME EVENT (FOR WW AND ZZ) N.B. NO CLEANUP CALLS FOR THESE!
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,X(3),W,WMAX,WSUM,
& X1MIN,X1MAX,X2MIN,X2MAX,QSCALE,GAMFAC,GLUFAC,R(3,3),CS,SN,M(3),
& E(3),LAMBDA,A,B,C,PTSQ,EM,P1(5),P2(5),PVRT(4),EPS,MASDEP
C IF (IOPT.EQ.1) SET UP EVENT RECORD
C IF (IOPT.EQ.2) CLEAN UP EVENT RECORD AFTER SHOWERING
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWBVMC,HWUALF,HWULDO,P1(5),P2(5),P3(5),
& PCMF(5),L(5),R(3,3),Q,XBJ,RN,XPMIN,XPMAX,XP,ZPMIN,ZPMAX,ZP,FAC,
& X1,X2,XTSQ,XT,PTSQ,SIN1,SIN2,W1,W2,CFAC,PDFOLD(13),PDFNEW(13),
SUBROUTINE HWBDYP(IOPT)
C MATRIX ELEMENT CORRECTIONS TO DRELL-YAN PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,PMODK,AZ,CZ,
& T,U,S,EM,TMIN,TMAX,PMOD2,GLUFAC,SMIN,SMAX,SZ,TEST,
& JAC,M(3),W1,W,PMOD3,SCAPR,CPHI,SPHI,SCALE,XI1,XI2,
C DELETES INTERNAL LINES FROM SHOWER, MAKES COLOUR CONNECTION INDEX
C AND COPIES INTO /HEPEVT/ IN COLOUR ORDER.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,ID,IJET,KHEP,IPAR,JPAR,NXPAR,IP,JP
IF (IERROR.NE.0) RETURN
C---SAVE VIRTUAL PARTON DATA
C HWBGEN EVOLVES QCD JETS ACCORDING TO THE METHOD OF
C G.MARCHESINI & B.R.WEBBER, NUCL. PHYS. B238(1984)1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWULDO,HWRGAU,EINHEP,ERTXI,RTXI,XF
INTEGER NTRY,LASHEP,IHEP,NRHEP,ID,IST,JHEP,KPAR,I,J,IRHEP(NMXJET),
& IRST(NMXJET),JPR
C for events read in from the GUPI (Generic User Process Interface)
C event common block
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER MAXNUP
PARAMETER (MAXNUP=500)
INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
C-----------------------------------------------------------------------
C COMBINES JETS WITH REQUIRED KINEMATICS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWULDO,EPS,PTX,PTY,PF,PTINF,PTCON,CN,CP,SP,PP0,
& PM0,ET0,DET,ECM,EMJ,EMP,EMS,DMS,ES,DPF,ALF,AL(2),ET(2),PP(2),
& PT(3),PA(5),PB(5),PC(5),PQ(5),PR(5),PS(5),RR(3,3),RS(3,3),ETC,
C On output PPAR(1-5,*) (if TMPAR(*)), containts:
C (P-trans,Xi or Xilast,P-long,E,M)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUSQR,EXI,PISQ,PJPK,EJEK,PTSQ,Z,ZMIN,ZMAX,
$ EMI,EMJ,EMK,C,NQ,HWBVMC,RHO,POLD,PNEW,EOLD,ENEW,A,B
INTEGER IPAR,JPAR,KPAR,MPAR,I,J,K
C BRANCHES TIMELIKE PARTON KPAR INTO TWO, PUTS PRODUCTS
C INTO NPAR+1 AND NPAR+2, AND INCREASES NPAR BY TWO
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUTAB,HWRUNI,HWULDO,PMOM,
& QNOW,QLST,QKTHR,RN,QQBAR,DQQ,QGTHR,SNOW,QSUD,ZMIN,ZMAX,ZRAT,WMIN,
& QLAM,Z1,Z2,ETEST,ZTEST,ENOW,XI,XIPREV,EPREV,QMAX,QGAM,SLST,SFNL,
C SUBROUTINE TO REPLACE HWBCON IN RPARITY VIOLATING SUSY
C BASED ON HWBCON BY BRW
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,IST,ID,JC,JD,JHEP,IDP,IDM,IDM2,
& RHEP,IST2,ANTC,XHEP,IP,COLP
LOGICAL BVVUSE,BVVHRD,BVDEC1,BVDEC2,COLRD,ACOLRD,BVHRD,BVHRD2,
C-----------------------------------------------------------------------
C--Function to find the right daugther of a decaying gluino
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ID,JHEP,KC,JC
LOGICAL COL,IFGO
C---N.B. WILL NEED MODS WHEN SUSY PARTICLES CAN SHOWER
C-----------------------------------------------------------------------
C--Function to search in the jet for the particle
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER JC,JD,QHEP,LHEP,IHEP,JHEP,IDM,NCOUNT,ID,IP,IDM2,COLP
LOGICAL CON,BVVUSE,FLA,AFLA,BVVHRD
FLA(IP) = (IP.LE.6.OR.(IP.GE.115.AND.IP.LE.120).
C DECPAR(2). RHOPAR(2) gives the jet spin density matrix.
C See I.G. Knowles, Comp. Phys. Comm. 58 (90) 271.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,DMIN,PT,EIKON,EISCR,EINUM,EIDEN1,EIDEN2,
& WT,SPIN,Z1,Z2,TR,PRMAX,CX,SX,CAZ,ROHEP(3),RMAT(3,3),ZERO2(2)
INTEGER IPAR,JPAR,KPAR,LPAR,MPAR,JSTR,LSTR,MSTR
C in hard subprocess, otherwise zero. Assignments based upon
C Comp. Phys. Comm. 58 (1990) 271.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION C,V12,V23,V13,TR,C1,C2,C3,R1(2),R2(2)
INTEGER IST
SAVE R1,R2,V12
C-----------------------------------------------------------------------
C COMPUTES (OR READS) TABLES OF SUDAKOV FORM FACTORS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUGAU,HWBVMC,HWBSUG,HWBSU1,HWBSU2,G1,G2,QRAT,
& QLAM,POWER,AFAC,QMIN,QFAC,QNOW,ZMIN,ZMAX,Q1,QCOLD,VGOLD,VQOLD,
& RMOLD(6),ACOLD,ZLO,ZHI
C THE SECOND ORDER ALPHAS CASE COMES FROM CONVERTING INTEGRAL OVER
C Q^2 INTO ONE OVER ALPHAS, WITH FLAVOUR THRESHOLDS.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBSUL,HWUALF,Z,QRAT,QLAM,U,AL,BL,QNOW,QMIN,
& BET(6),BEP(6),MUMI(6),MUMA(6),ALMI(6),ALMA(6),FINT(6),ALFINT,
& MUMIN,MUMAX,ALMIN,ALMAX
C due to spin (if AZSPIN) using the algorithm of Knowles & Collins.
C Ses Nucl. Phys. B304 (1988) 794 & Comp. Phys. Comm. 58 (1990) 271.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,DMIN,PT,EIKON,EINUM,EIDEN1,EIDEN2,EISCR,
& WT,SPIN,Z1,Z2,PRMAX,CAZ,CX,SX,ROHEP(3),RMAT(3,3),ZERO2(2)
INTEGER INITBR,INTERF,IPAR,JPAR,KPAR,LPAR,MPAR,NTRY,JOLD
C-----------------------------------------------------------------------
SUBROUTINE HWBTOP
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUSQR,X(3),W,
& X3MIN,X3MAX,X1MIN,X1MAX,QSCALE,GLUFAC,R(3,3),M(3),
& E(3),AW,PTSQ,EM,EPS,MASDEP,A,B,C,GAMDEP,LAMBDA,
C-----------------------------------------------------------------------
C VIRTUAL MASS CUTOFF FOR PARTON TYPE ID
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC
INTEGER ID
IF (ID.EQ.13) THEN
C Subroutine to split a baryonic cluster containing two heavy quarks
C Based on HWCCUT
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUPCM,HWRGEN,HWVDOT,EMC,QM1,QM2,QM3,QM4,
& PXY,PCX,PCY,RCM,PCL(5),AX(5),PA(5),PB(5),PC(5),
& VCLUS(4),DQM,EMX,EMY,SKAPPA,RKAPPA,VTMP(4),
C FINDS UNPAIRED PARTONS AFTER BARYON-NUMBER VIOLATION
C MODIFIED FOR RPARITY VIOLATING SUSY
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
COMMON/HWBVIC/NBV,IBV(18)
DOUBLE PRECISION HWRGEN,PDQ(5)
INTEGER NBV,IBV,JBV,KBV,LBV,IHEP,IP1,IP2,IP3,JP1,JP2,JP3,
C-----------------------------------------------------------------------
C Function to find the baryon number violating vertex a parton came from
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER HWCBVT,IP,JP(2),KP,I,J,ID,TYPE,IDM,IDM2,IDM3,IDM4
JP(1) = IP
ID = IDHW(IP)
C-----------------------------------------------------------------------
C Subroutine to correct colour connections after the gluon splitting
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,STFSPT,LHEP,MHEP,RHEP
IF(IERROR.NE.0) RETURN
C--Find the first particle in the event record with status 150
C-----------------------------------------------------------------------
C Cuts into 2 the cluster, momentum PCL, made of partons JHEP & KHEP
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWREXQ,HWUPCM,HWRGEN,HWVDOT,EMC,QM1,QM2,EMX,EMY,
& QM3,PXY,PCX,PCY,RCM,PCL(5),AX(5),PA(5),PB(5),PC(5),SKAPPA,DELTM,
& VSCA,VTMP(4),RKAPPA,VCLUS
C-----------------------------------------------------------------------
C DECAYS CLUSTERS INTO PRIMARY HADRONS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER JCL,KCL,IP,JP,KP,IST,ID1,ID2,ID3
IF (IERROR.NE.0) RETURN
IF (IPRO/10.EQ.9.OR.IPRO/10.EQ.5) THEN
C MHS's cluster vertices
C MODIFIED 16/10/97 BY BRW FOR SUSY PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWULDO,HWVDOT,HWRGEN,HWUPCM,DCL0,DCL(4),DCL1,
& DFAC,DISP1(4),DISP2(4),DMAX,PCL(5),DOT1,DOT2,FAC,VCLUS,SCA1,SCA2,
& EM0,EM1,EM2,PC0,PC1
C SPLITS ANY TIMELIKE GLUONS REMAINING AFTER PERTURBATIVE
C BRANCHING INTO LIGHT (I.E. U OR D) Q-QBAR PAIRS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,PF
INTEGER HWRINT,IHEP,JHEP,KHEP,LHEP,MHEP,ID,J,IST
EXTERNAL HWRGEN,HWRINT
C
C MODIFIED 15/11/99 TO SMEAR POSITIONS OF HADRONS BY 1/(CLUSTER MASS)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRGAU,HWVDOT,EM0,EM1,EM2,EMADU,EMSQ,
& PCMAX,PCM,PTEST,PCQK,PP(5),EMLOW,RMAT(3,3),CT,ST,CX,SX,HPSMR
INTEGER HWRINT,JCL,ID1,ID2,ID3,ID,IR1,IR2,NTRY,IDMIN,IMAX,I,MHEP,
C-----------------------------------------------------------------------
C Computes the width and maximum weight for a two body mode
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IMODE,I
DOUBLE PRECISION A(2),M(3),PCM,E1,E2,HWUPCM,PHS,WGT,MWGT,PCM2,
& M2(3)
C Subroutine to perform the three body decays for spin correlations
C and SUSY three body modes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IMODE,I,J,ID,IDP(4+NDIAGR),ITYPE,NDIA,ID1,ID2,
& DRTYPE(NDIAGR),NTRY,IDSPIN,NCTHRE,DRCF(NDIAGR)
DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,WGT,M342,HWRUNI,
C-----------------------------------------------------------------------
C Subroutine to calculate the matrix element for a given mode
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,P0,P1,P2,P3,P0P,IB,ID,IDP(4+NDIAGR),IDSPIN,NDIA,
& DRTYPE(NDIAGR),NCTHRE,DRCF(NDIAGR)
DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,WGT,FJAC,M342,HWRUNI,
C Subroutine to calculate the helicity amplitudes for the three body
C gauge boson exchange diagram
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,C(2,2),E(2,2),ZI,APP(2,2),APM(2,2),
& AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for the three body
C Higgs boson exchange diagram
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for the three body
C antisfermion exchange diagram
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for the three body
C sfermion exchange diagram
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for the three body
C gauge boson exchange diagram (antiparticle decay)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),
& F0M(2,2,8),F2(2,2,8),PRE,C(2,2),E(2,2),ZI,APP(2,2),APM(2,2),
& AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8),F01(2,2,8,8)
C Subroutine to calculate the helicity amplitudes for the three body
C gauge boson exchange diagram
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),
& F0M(2,2,8),F2(2,2,8),PRE,C(2,2),ZI,APP(2,2),APM(2,2),
& AMP(2,2),AMM(2,2),F1M(2,2,8),F3(2,2,8),F01(2,2,8,8)
C Subroutine to calculate the helicity amplitudes for the three body
C decay fermion --> gravitino fermion antifermion (via gauge boson)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,F1M(2,2,8),F3(2,2,8)
DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,M342,M232,M242,MR,
C Subroutine to calculate the helicity amplitudes for 1st 3 body RPV
C diagram f--> fbar fbar f
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 2nd 3 body RPV
C diagram f --> fbar fbar f
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 3rd 3 body RPV
C diagram f --> fbar fbar f
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 4th 3 body RPV
C diagram f --> f f f
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 5th 3 body RPV
C diagram f --> f f f
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 6th 3 body RPV
C diagram f --> f f f
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 7th 3 body RPV
C diagram f --> fbar fbar fbar
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 8th 3 body RPV
C diagram f --> fbar fbar fbar
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,V1(2,2),V2(2,2),ZI,F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for 9th 3 body RPV
C diagram f --> fbar fbar fbar
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C Subroutine to calculate the helicity amplitudes for the three body
C Higgs boson exchange diagram antifermion decay
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ME(2,2,2,2),S,D,F0(2,2,8),F1(2,2,8),F01(2,2,8,8),
& F0M(2,2,8),F2(2,2,8),PRE,ZI,V1(2,2),V2(2,2),F1M(2,2,8),
& F3(2,2,8)
C-----------------------------------------------------------------------
C Subroutine to perform the four body Higgs decays
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IMODE,I,J,ID,IDP(4+NDIAGR),ITYPE(2),NTRY,ITYPE1,ITYPE2
DOUBLE PRECISION A,B,MS,MWD,M,M2,WGT,HWRUNI,BRW(6),BRZ(12),
& HWUPCM,WMAX,WSUM,WSSUM,MR,PRE,TEMP,HWRGEN,WTMAX,P(5,5)
C Subroutine to calculate the matrix element for a given four body
C decay mode
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,P0,P1,P2,P3,ID,O(2),IDP(4+NDIAGR),II,P4
DOUBLE PRECISION A,B,MS,MWD,M,M2,WGT,HWRUNI,
& M23,PCMA,PCMB(2),HWUPCM,PHS,N(3),HWVDOT,PP,HWULDO,EPS,
C IF BOSON CAME FROM W+1JET, GIVE IT THE CORRECT DECAY CORRELATIONS
C--BRW FIX 20/07/04: ADD FULL DECAY CORRELATIONS FOR W/Z+HIGGS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWULDO,R(3,3),CV,CA,BR,PCM,
& PBOS(5),PMAX,PROB,RRLL,RLLR
INTEGER HWRINT,IBOS,IBOSON,IPAIR,ICMF,IOPT,IHEL,IMOTH,
C CHOOSE DECAY MODE OF BOSON
C IOPT=2 TO RESET COUNTERS, 1 FOR BOSON PAIR, 0 FOR ANY OTHERS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,BRMODE(12,3),CV,CA,BR,BRLST,BRCOM,FACZ,
& FACW
INTEGER HWRINT,IDBOS,IDEC,IDMODE(2,12,3),IFER,IANT,IOPT,I1,I2,
C IDENTICAL TO HWDBOZ BUT REQUIRES DECAY MODE ACCESSIBLE FOR GIVEN
C MASS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,BRMODE(12,3),CV,CA,BR,BRLST,BRCOM,FACZ,
& FACW,MSMODE(12,3),MASS
INTEGER HWRINT,IDBOS,IDEC,IDMODE(2,12,3),IFER,IANT,IOPT,I1,I2,
C Checks line L of decay table is compatible with decay of particle
C IDKY, tidies up the line and sets NPRODS.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION EPS,QS,Q,DM
INTEGER IDKY,L,IFAULT,I,ID,J
LOGICAL IFGO
C-----------------------------------------------------------------------
C INTERFACE TO QQ-CLEO MONTE CARLO (LS 11/12/91)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,IIHEP,NHEPHF,QQLMAT
LOGICAL QQLERR
CHARACTER*8 NAME
C-----------------------------------------------------------------------
C INTERFACE TO EURODEC PACKAGE (LS 10/29/91)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,IIHEP,NHEPHF,IEUPDG,IPDGEU
CHARACTER*8 NAME
C---EURODEC COMMON'S : INITIAL INPUT
C-----------------------------------------------------------------------
C Generates 4-body decay 0->1+2+3+4 using pure phase space
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,P0(5),P1(5),P2(5),P3(5),P4(5),B,C,AA,BB,
& CC,DD,EE,TT,S1,RS1,FF,S2,PP,QQ,RR,P1CM,P234(5),P2CM,P34(5),P3CM
INTEGER NTRY
C-----------------------------------------------------------------------
C Generates 5-body decay 0->1+2+3+4+5 using pure phase space
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,P0(5),P1(5),P2(5),P3(5),P4(5),P5(5),B,C,
& AA,BB,CC,DD,EE,FF,TT,S1,RS1,GG,S2,RS2,HH,S3,PP,QQ,RR,SS,P1CM,
& P2345(5),P2CM,P345(5),P3CM,P45(5),P4CM
C GENERATES DECAYS OF UNSTABLE HADRONS AND LEPTONS
C Modified for TAUOLA interface 16/10/01 PR
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
COMMON/FFS/TB,BT
COMMON/SFF/IT1,IB1,IT2,IB2
DOUBLE PRECISION TB,BT
C CALCULATE THE COMPLEX FUNCTION F OF HHG eq 2.18
C FOR USE IN H-->GAMMGAMM DECAYS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION TAU,FNREAL,FNIMAG,FNLOG,FNSQR
IF (TAU.GT.ONE) THEN
FNREAL=(ASIN(1/SQRT(TAU)))**2
C CALCULATE THE DOUBLE BREIT-WIGNER INTEGRAL
C X=(EMV/EMH)**2 , Y=EMV*GAMV/EMH**2
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWDHGF,X,Y,CHANGE,X1,X2,FAC1,FAC2,TH1,TH2,TH1HI,
& TH1LO,TH2HI,TH2LO,X2MAX,SQFAC
INTEGER NBIN,IBIN1,IBIN2
C B) FOR GAMinp>0 CALCULATE TOTAL HIGGS WIDTH
C FOR EMH=GAMINP. STORE RESULT IN GAMINP.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWDHGF,HWRGEN,HWRUNI,HWUSQR,HWUPCM,GAMINP,EMH,
& EMF,COLFAC,ENF,K1,K0,BET0,BET1,GAM0,GAM1,SCLOG,CFAC,XF,EM,GAMLIM,
& GAM,XW,EMW,XZ,EMZ,YW,YZ,EMI,TAUT,TAUW,WIDHIG,VECDEC,EMB,GAMB,
C MODIFIED TO CALL A NUMBER OF ROUTINES TO DO THE VARIOUS BITS OF
C THE PROCESS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION PW(5)
INTEGER IHEP,IS,ID,IM,KHEP,LHEP,MHEP,NPR,CLSAVE(2),NHEPST
LOGICAL FOUND
C IE to select the decay mode
C was part of HWDHOB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUMBW,HWRGEN,SDKM,RN,BF
INTEGER IST(3),IHEP,ID,IM,I,JHEP,LHEP,MHEP,NPR,MTRY,NTRY,IS
EXTERNAL HWRGEN
C IE generate the kinematics for the decay
C was part of HWDHOB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
COMMON/FFS/TB,BT
COMMON/SFF/IT1,IB1,IT2,IB2
DOUBLE PRECISION TB,BT
C IE setup the colour connections
C was part of HWDHOB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ID,IM,KHEP,LHEP,MHEP,NPR,CLSAVE(2)
IF (IERROR.NE.0) RETURN
C Colour connections
C IE parton-showers with special treatment for top
C was part of HWDHOB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION PW(5),PDW(5,3)
INTEGER IHEP,ID,IM,I,KHEP,LHEP,MHEP,NPR,NTRY,WHEP,SHEP
DOUBLE COMPLEX RHOIN(2,2)
C IE sort out RPV colour connections
C was part of HWDHOB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ID,LHEP,MHEP,IDM,IDM2,THEP,CLSAVE(2)
IF (IERROR.NE.0) RETURN
C--New to correct colour connections in Rslash
C Subroutine to perform the final part of the heavy object decays
C IE sort out any colour connection problems
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,IM,JHEP,ISM,JCM
IF (IERROR.NE.0) RETURN
C Fix any SUSY colour disconnections
C either, meson/baryon spectator model weak decays;
C or, quarkonia -> 2-gluons, q-qbar, 3-gluons, or 2-gluons + photon.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
COMMON/FFS/TB,BT
COMMON/SFF/IT1,IB1,IT2,IB2
DOUBLE PRECISION TB,BT
C-----------------------------------------------------------------------
C Sets the colour connections in Baryon number violating decays
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,MHEP,ID,ID2,IDM2,IDM3,COLCON(2,2,3),FLACON(2,3),JHEP,
& DECAY,COLANT,KHEP,IDM,IDMB,IDMB2,IDMB3,IDMB4,QHEP,IDM4,
& CLSAVE(2),XHEP,I,HWRINT,THEP
C-----------------------------------------------------------------------
C SUBROUTINE TO IMPLEMENT ALL RPARITY DECAY MATRIX ELEMENTS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION SM(6),SW(6),HWULDO,INFCOL,AM, M12SQ,M23SQ,MSGN,
& M13SQ,A(6),B(6),SWEAK,MW,DECMOM(5),TEST(3),EPS,
& M12SQT(6),M23SQT(6),M13SQT(6),LIMIT,M(4),RAND,
C-----------------------------------------------------------------------
C Subroutine to perform decays including spin correlations
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION PW(5)
INTEGER IDEC,IP,IS,IHEP,ID,IM,LHEP,MHEP,NPR,KHEP,CLSAVE(2),NTRY,
& ID1
C-----------------------------------------------------------------------
C Subroutine to check a vertex and decide which branch to treat
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IDEC,I,IPICK(5),IP,HWRINT,P1,P2,P3,P4,P3P,P4P,NPR,P0,P0P,
& P1P,P2P,IF1,IF2,P5,P5P
DOUBLE PRECISION NORM
C including spin correlations
C was part of HWDHOB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,PW(5),HWDPWT,HWDWWT,PCM,HWUPCM
INTEGER IHEP,IM,KHEP,MHEP,NPR,ISN,RHEP
EXTERNAL HWRGEN,HWDPWT,HWDWWT,HWUPCM
C if using TAUOLA selects the spin and uses TAUOLA to perform the
C decay
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IP,IHEP,ID1,ID,NTRY
DOUBLE PRECISION PPOL,HWRGEN,POL
EXTERNAL HWRGEN
C Subroutine to calculate the two body matrix element for spin
C correlations
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IOUT1,IOUT2,IMODE,IDSPIN,ID,I,J,IDP(3),P0,P1,P2,O(2),P0P,
& NTRY
DOUBLE PRECISION XMASS,PLAB,PRW,PCM,PREF(5),P(5,3),PM(5,3),PCMA,
C Master subroutine for three body SUSY and spin ME's
C Uses HWD3ME to generate the momenta etc
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX F0(2,2,8),F1(2,2,8),F1M(2,2,8),F3(2,2,8),
& F0M(2,2,8),F2(2,2,8),RHOIN(2,2),F01(2,2,8,8)
DOUBLE PRECISION A,B,MS,MWD,MA,MB,MA2,MB2,M342,M232,M242,MR,
C IOPT = 1 select decay mode and generate momenta
C IOPT = 2 enter first decays and perform parton shower
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IOPT,ID,IOUT1,IOUT2,IB(2),I,IDF(4),ITYPE(2),IMODE,
& IDP(4+NDIAGR),ID1,ID2,J
DOUBLE PRECISION CV,CA,A,B,MS,MWD,MR,M,M2,P(5,5),PW(5,2),BR
C IOPT = 1 performs decay
C IOPT = 2 write outs final TAUOLA information
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IOPT,IHEP,ID,ITAU,I,IMO,NHEPPO
DOUBLE PRECISION POL
REAL POL1(4)
C-----------------------------------------------------------------------
C DECIDES WHETHER TO DO TOP QUARK DECAY BEFORE HADRONIZATION
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
LOGICAL DECAY
DECAY=RMASS(6).GT.130D0
END
C MATRIX ELEMENT SQUARED FOR
C ((V-A)*TB1+(V+A)*CT1)*((V-A)*TB2+(V+A)*CT2)) WEAK DECAY
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
COMMON/FFS/TB,BT
COMMON/SFF/IT1,IB1,IT2,IB2
DOUBLE PRECISION TB,BT
C Sets STAB=.TRUE. if DKVRTX lies outside the specified region.
C Revised 05/09/00 by BRW to put parameters in common
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION DKVRTX(4),RR
LOGICAL STAB
STAB=.FALSE.
C Modified 28/03/01 by BRW to handle negative weights
C Modified 15/07/02 by PR for Les Houches Accord
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I
DOUBLE PRECISION RNWGT,SPWGT,ERWGT
C--Les Houches Common Block
C GENERATES A PHOTON IN WEIZSACKER-WILLIAMS (WWA=.TRUE.) OR
C ELSE EQUIVALENT PHOTON APPROX FROM INCOMING E+, E-, MU+ OR MU-
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,EGMIN,ZMIN,ZMAX,ZGAM,SS,ZMI,ZMA,
& PPL,PMI,QT2,Q2,QQMIN,QQMAX,S0,A
INTEGER IHEP,IHADIS
C-----------------------------------------------------------------------
C FIND MINIMUM INVARIANT MASS SQUARED NEEDED FOR HARD PROCESS, S0
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION S0,RPM(2)
INTEGER HQ,I
IF (IPRO.EQ.13.OR.IPRO.EQ.14) THEN
C INITIALISES ELEMENTARY PROCESS
C Modified 28/03/01 by BRW to handle negative weights
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRSET,DUMMY,SAFETY
EXTERNAL HWRSET
PARAMETER (SAFETY=1.001)
C-----------------------------------------------------------------------
C GENERATES AN ISR PHOTON FROM INCOMING E+, E-, MU+ OR MU-
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION CIRCKP(2)
COMMON /HWCIR2/CIRCKP
DOUBLE PRECISION HWRGEN,QSQMAX,QSQMIN,A,B,B1,B2,B3,B4,B5,B6,B7,B8,
C-----------------------------------------------------------------------
C SETS UP 2->1 (COLOUR SINGLET) HARD SUBPROCESS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION PA
INTEGER ICMF,I,IBM,IHEP
C---INCOMING LINES
C OTHERWISE, CHOOSES AND LOADS ALL VARIABLES FOR HARD PROCESS
C modifications for Les Houches accord by PR (7/15/02)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION CIRCKP(2)
COMMON /HWCIR2/CIRCKP
DOUBLE PRECISION Z1,Z2,C1,C2,B1,B2,CIRCEE,CIRCGG,RS,MISS,ETA,
C SETS UP 2->2 HARD SUBPROCESS
c BRW change 18/8/04: BW smearing of mass i only if SMRi is true
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUMBW,HWUPCM,PA,PCM
INTEGER ICMF,IBM,I,J,K,IHEP,NTRY
LOGICAL SMR3,SMR4
C
C...First release: 1-APR-1998 by Stefano Moretti
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I
DOUBLE PRECISION P1(0:3),P2(0:3),P3(0:3),P4(0:3)
DOUBLE PRECISION P(0:3)
C I-L are the particles (all outgoing)
C Z1 and Z2 are the decay products of the Z
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ND,I,J,K,L,Z1,Z2
DOUBLE COMPLEX HWH2DD,ZI,S,D,F
PARAMETER(ZI=(0.0D0,1.0D0))
C...First release: 01-APR-1998 by Kosuke Odagiri
C...First modified: 12-APR-1998 by Stefano Moretti
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER MU,IRES,IFL
DOUBLE PRECISION P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
DOUBLE PRECISION EMB,EMT,EMW,EMH,EMH01,EMH02,EMH03
C Subroutine to implement the F function of Eijk and Kliess
C fixed first momenta and all second momenta
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
DOUBLE COMPLEX F(2,2,8),S,D,SIP(2),SJP(2)
INTEGER I,J,NP
C Subroutine to implement the F function of Eijk and Kliess
C fixed second momenta and all first momenta
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
DOUBLE COMPLEX F(2,2,8),S,D,SIP(2),SJP(2)
INTEGER I,J,NP
C Subroutine to implement the F function of Eijk and Kliess
C All first and second momenta
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION P(5),MQ,PM(5),XMASS,PLAB,PRW,PCM,HWULDO,PDOT,EPS
DOUBLE COMPLEX F(2,2,8,8),SIP(2),SJP(2),S,D
INTEGER I,J,NP
C Massless matrix elements for gg-->qqZ and qq-->qqZ
C using the matrix elements given in Nucl. Phys. B262 (1985) 235-242
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IQ,I,J,OZ(2,2),IDZ,P1,P2,P3,P4,IQI,ID(2),K
DOUBLE PRECISION MG(2),MQ(2,5),G(12,2),FLOW(3,3),CQFC,CQIFC,
& CGFC,CGIFC
C-----------------------------------------------------------------------
C Massive matrix elements for gg --> qqbarZ and qqbar --> qqbarZ
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IQ,I,IDZ,P1,P2,PL,PB,PBB,O(2),J,IQI
DOUBLE PRECISION MG(2),MQ(2,5),G(12,2),CQFC,CQIFC,CGFC,CGIFC,
& PTMP(5,10),XMASS,PLAB,PRW,PCM,HWULDO,QBL,QBBL,Q2B,Q1B,Q2BB,
C-----------------------------------------------------------------------
C Phase Space for vector boson plus 2 jets
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION WEIGHT,XMASS,PLAB,PRW,PCM,Y(3),Y35,Y34,Y45,RAND,
& HWRGEN,HWRUNI,M35,M35S,G(IMAXCH),DEM,MT(3),PT(3),MJAC,ETOT,
& STOT,MQ(3),MQ2(3),PS35,HWUPCM,TWOPI2,MT35,PTJ(3),MT2(3),A,C,
C IOPT=1 return the function at M2
C IOPT=2 calculate M2
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IOPT
DOUBLE PRECISION FJAC,MPOW,MMN,MQ2,M2,A1,A01,RPOW,QPOW,HWRGEN,MMX
EXTERNAL HWRGEN
C IOPT=1 return the function at M2
C IOPT=2 calculate M2
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IOPT
DOUBLE PRECISION FJAC,MPOW,A1,A01,RPOW,QPOW,HWRGEN,PT2,
& PPOW,PTMN2,PTMX2,Z
C Subroutine to calculate the spinor products in the notation of
C Kleiss and Strirling S(1) is S and S(2) is T
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION WRN(2),K(5),KK(5),P(5,2),Q1,Q2,EPS,QTI,PTI,
& PT,QT,DPM,DMP,QP,QM,P1,P2,PP,PM
DOUBLE COMPLEX S(2),ZI,Z1,ZT,ZQ,ZQS,ZPS,ZP,ZDPM,ZDMP
C Z1 and Z2 are the decay products of the Z
C P1 is the polarization of the line I,J
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T1,S,D
INTEGER I,J,K,L,Z1,Z2,P1
COMMON/HWHEWS/S(8,8,2),D(8,8)
C P1 is the polarization of the line I,J
C P2 is the polarization of the gluon K
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T2,S,D
INTEGER I,J,K,L,Z1,Z2,P1,P2
DOUBLE PRECISION B(6)
C P1 is the polarization of the line I,J
C P2 is the polarization of the gluon K
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T3,S,D
INTEGER I,J,K,L,Z1,Z2,P1,P2
DOUBLE PRECISION B(6)
C P1 is the polarization of the line I,J
C P2 is the polarization of the line K,L
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T4,AP,AM,S,D
INTEGER I,J,K,L,Z1,Z2,J1,J2,J3,J4,P1,P2
DOUBLE PRECISION B(6)
C P1 is the polarization of the line I,J
C P2 is the polarization of the line K,L
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T5,AP,AM,S,D
INTEGER I,J,K,L,Z1,Z2,J1,J2,J3,J4,P1,P2
DOUBLE PRECISION B(6)
C P2 is the polarization of the gluon K
C P3 is the polarization of the gluon L
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T6,S,D
INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
DOUBLE PRECISION B(6)
C P2 is the polarization of the gluon K
C P3 is the polarization of the gluon L
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T7,S,D
INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
DOUBLE PRECISION B(6)
C P2 is the polarization of the gluon K
C P3 is the polarization of the gluon L
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T8,S,D
INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
DOUBLE PRECISION B(6)
C P2 is the polarization of the gluon K
C P3 is the polarization of the gluon L
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T9,S,D
INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
DOUBLE PRECISION B(6)
C P2 is the polarization of the gluon K
C P3 is the polarization of the gluon L
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWH2T0,S,D
INTEGER I,J,K,L,Z1,Z2,J1,J2,P1,P2,P3
DOUBLE PRECISION B(6)
C
C WARNING: Phase space factor inaccurate for JADE y_cut > 0.14.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER LM,LP,IQK,I,J,IDMN,IDMX,ID1,ID2,IST(4)
DOUBLE PRECISION HWRGEN,HWUALF,HWUAEM,HWULDO,HWH4J1,HWH4J2,
& HWH4J4,HWH4J5,HWH4J6,HWH4J7,QNOW,Q2NOW,QLST,SCUT,PSFAC,FACT,
C Mean EVWGT = cross section in nanoBarn
C
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,Y,Q2,SHAT,Z,PHI,AJACOB,DSIGMA,ME,MP,
& ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,FSIGMA(18),
& SIGSUM,PROB,PRAN,PVRT(4),X
C of the hard 2-->3 subprocess, to match with HERWIG routines of
C jet evolution.
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUECM,HWUPCM,HWUSQR,Y,Q2,SHAT,Z,PHI,AJACOB,
& DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
& PGAMMA(5),SG,MF1,MF2,EP,PP,EL,PL,E1,E2,Q1,COSBET,SINBET,COSTHE,
C Fill the logical vector INSIDE to tag contributing subprocesses
C to the cross-section
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
LOGICAL IFGO
DOUBLE PRECISION HWRUNI,HWRGEN,HWUPCM,Y,Q2,SHAT,Z,PHI,AJACOB,
& DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
C Returns differential cross section DSIGMA in (Y,Q2,ETA,Z,PHI)
C Scale for structure functions and alpha_s selected by BGSHAT
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUALF,HWUAEM,Y,Q2,SHAT,Z,PHI,AJACOB,DSIGMA,
& ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT,
& SFUN(13),ALPHA,LDSIG,DLQ(7),SG,XG,MF1,MF2,MSUM,MDIF,MPRO,FFUN,
C----------------------------------------------------------------------
C DEEP INELASTIC LEPTON-HADRON SCATTERING: MEAN EVWGT = SIGMA IN NB
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,PRAN,PROB,SAMP,SIG,Q2,
& XBJ,Y,W,S,MLEP,MHAD,MLSCAT,YMIN,YMAX,XXMAX,Q2JAC,XXJAC,
& JACOBI,A1,A2,A3,B1,B2,PCM,PCMEP,PCMLW,PCMEQ,PCMLQ,COSPHI,PA,
C Modified 16/01/01 by BRW to implement Peter Richardson's
C fix for bug in lepton mass effects on branching ratio
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUAEM,EPS,C1,C2,C3,EMSQZ,EMGMZ,
& EMSQZP,EMGMZP,CQF(7,6,16),QPOW,RPOW,A01,A1,A02,A2,A03,A3,CRAN,
& EMJ1,EMJ2,EMJ3,EMJAC,FACT,QSQ,HCS,FACTR,RCS,EXTRA,PMAX,PTHETA
C-----------------------------------------------------------------------
C Drell-Yan production with a q qbar pair
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,MAP(12),ORD,IFL,IDP(6),IFLOW,QCFL(2,2),GCFL(2),IDZ,IQ
DOUBLE PRECISION HCS,RCS,MQ(2,5),HWRGEN,G(12,2),DIST(2),MG(2)
LOGICAL FSTCLL,MASS
C AND COS(THETA) IN CENTRE-OF-MASS SYSTEM
C AND TIMES BRANCHING FRACTION IF WW
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWULDO,EMSQ,BETA,S,T,U,TMIN,TMAX,TRAT,
& DSDT,PROB,X,Z(2),ZMIN,ZMAX,PCMIN,PCMAX,PCFAC,PLOGMI,PLOGMA,PTCMF,
& Q,PC,BLOG,EMCMIN,EMCMAX,EMLMIN,EMLMAX,WGT(6),RWGT,CV,CA,BR,QT(2),
C W + GAMMA --> FF'BAR : MEAN EVWGT = CROSS SECTION IN NANOBARN
C BASED ON BOSON GLUON FUSION OF ABBIENDI AND STANCO
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,GMASS,EV(3),RV,Y,Q2,SHAT,Z,PHI,AJACOB,
& DSIGMA,ME,MP,ML,MREMIF(18),MFIN1(18),MFIN2(18),RS,SMA,W2,RSHAT
INTEGER LEP
C-----------------------------------------------------------------------
C COMPUTES DIFFERENTIAL CROSS SECTION DSIGMA IN (Y,Q2,ETA,Z,PHI)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION TMAX,TMIN,A1,A2,B1,B2,I0,I1,I2,I3,I4,I5,MUSQ,
& MDSQ,ETA,Q1,COSTHE,S,G,T,U,C1,C2,D1,D2,F1,F2,COSBET,WPROP,D(4,4),
& C(4,4),QU,QD,QE,QW,PHOTON,EMWSQ,EMSSQ,CFAC,Y,Q2,SHAT,Z,PHI,
C If fermion flavour specified mass effects fully included.
C EVWGT=sig(e+e- --> ffbar) in nb
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWUAEM,Q2NOW,Q2LST,FACTR,
& VF2,VF,CLF(7),PRAN,PQWT,PMAX,PTHETA,SINTH2,CPHI,SPHI,C2PHI,S2PHI,
& PPHI,SINTH,PCM,PP(5),EWGT
C If flavour specified mass effects fully included.
C EVWGT=sig(e^-e^+ --> qqbar g) in nb
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWUALF,HWUAEM,HWULDO,HWDPWT,Q2NOW,Q2LST,
& PHASP,QGMAX,QGMIN,FACTR,QM2,CLF(7),ORDER,PRAN,PQWT,QQG,QBG,SUM,
& RUT,QQLM,QQLP,QBLM,QBLP,DYN1,DYN2,DYN3,DYN4,DYN5,DYN6,XQ2,X2SUM,
C-----------------------------------------------------------------------
C SUSY E+E- -> 2 SLEPTON PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWUAEM,EPS,HCS,RCS,S,PF,QPE,HWUPCM,PCM,
& FACTR,SN2TH,MZ,ME2(2,2,6),EMSC2,HWUMBW,HWRUNI,T,SQPE
INTEGER ID1,ID2,IL,IL1,IL2,I,J,IG,IG1,IHEP,NTRY,IDL,ILP,IDLR(2),
C-----------------------------------------------------------------------
C SUSY E+E- -> 2 GAUGINO PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWUAEM,HCS,RCS,MNU(4),MNU2(4),HWRUNI,
& FACA,M1(4,4),S2W,XC(4),XD(4),MSNU,
& MZ,HWHSS2,U,T,QPE,SQPE,MSL,MSL2,MSR,MSR2,
C-----------------------------------------------------------------------
C SUSY E+E- -> 2 SPARTICLE PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION SAVWT(3),RANWT,HWRGEN
EXTERNAL HWRGEN
SAVE SAVWT
C-----------------------------------------------------------------------
C SUSY E+E- -> 2 SQUARK PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWUAEM,EPS,HCS,RCS,S,PF,QPE,HWUPCM,PCM,
& FACTR,SN2TH,MZ,ME2(2,2,6),EMSC2,HWUMBW,HWRUNI,SQPE
INTEGER ID1,ID2,IQ,IQ1,IQ2,I,J,IHEP,IDL,IDLR(2),IDSQU(2),NTRY
C-----------------------------------------------------------------------
SUBROUTINE HWHEW0(IP,ETOT,XM,PR,WEIGHT,CR)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,ETOT,XM(2),PR(5,2),WEIGHT,CR,XM1,XM2,S,
& D1,PABS,D,CX,C,E,F,SC,G
INTEGER IP,I
C COUPLING (E**8/4/9) MUST BE INCLUDED AS WELL AS COMPENSATION
C FOR ON POLE APPROXIMATION AS DESIRED.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWHEW4,ZH,ZCH,ZD,ZAMP1,ZAMP3,DWW,CWW,BWW,AWW,
& AWWM,AWWP,AMPTEM,ZTWO,ZHALF
DOUBLE PRECISION XW,ZMASS,T3,EQ1,RR,RL,ZM2,AMP2,RKW,COLFAC(4),
C-----------------------------------------------------------------------
C E+E- -> W+W-/Z0Z0 (BASED ON ZOLTAN KUNSZT'S PROGRAM)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ZH,ZCH,ZD
DOUBLE PRECISION HWUAEM,HWRGEN,HWUPCM,ETOT,STOT,FLUXW,GAMM,GIMM,
& WM2,WXMIN,WX1MAX,WX2MAX,FJAC1,FJAC2,WX1,WX2,WMM1,WMM2,XXM,W2BO,
C-----------------------------------------------------------------------
C Hadron-Hadron to WW/WZ/ZZ (BASED ON ZOLTAN KUNSZT'S PROGRAM)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX ZH,ZCH,ZD
DOUBLE PRECISION HWUAEM,HWRGEN,HWUPCM,FLUXW,CSW,XMASS,
& PLAB,PRW,PCM,HWRUNI,P(5,10),AMPWW,DIST(4),MW2,CFAC1,AMP,
C ICH specifies the channel to be used if OPTM=.TRUE.
C This is used in optimising the weights for the different channels
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ICH,IDBOS,ISM(2,IMAXCH),I,J,IB(2),IDRES,IDP,IOPT,IPRC,ID1
DOUBLE PRECISION XMASS,PLAB,PRW,PCM,RAND,HWRGEN,BMS2(2),TJAC,PLM,
& MJAC(2),TWOPI2,SJAC,STOT,THAT,UHAT,TMIN,TMAX,UMIN,UMAX,PS(2),
C IOPT=1 return the function at MBOS2
C IOPT=2 calculate MBOS2
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IDBOZ,ISM,IOPT
DOUBLE PRECISION MBOZ,FJAC,GBOZ,GMBOZ,MPOW,MMIN,
& MBOS2,A1,A2,A01,A02,RPOW,QPOW,HWRGEN,MMAX,EMSQ
C-----------------------------------------------------------------------
C WW cross section in hadron hadron
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HCS,RCS,HWRGEN,DIST(2),CFAC,WAMP(2),S34,S56,KWW2,
& MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),TAUI(2),
& CSW,CFAC1
C-----------------------------------------------------------------------
C ZZ cross section in hadron hadron
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION AMP(2),RCS,HCS,HWRGEN,DIST(2),S34,S56,CFAC,
& MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),TAUI(2),
& CSW,CFAC1
C-----------------------------------------------------------------------
C WZ cross section in hadron hadron
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION AMP(2),HCS,RCS,HWRGEN,W34,DIST(2),S34,S56,CFAC,
& TCS,S12,MW2,MZ2,GMW,GMZ,G(4,2),EE(4),CKM2(12),RF(2),LF(2),
& TAUI(2),CSW,CFAC1
C-----------------------------------------------------------------------
C Subroutine to select t or u for HWHGBP
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IOPT
DOUBLE PRECISION FJAC,TPOW,TMIN,T,A1,A01,RPOW,QPOW,HWRGEN,TMAX,TN,
& TX,MT
C If GAMGRV is zero, it is revaluated during the first run
C MEAN EVWGT = SIGMA IN NB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,EPS,EMSQG,
& EMGMG,S,CC,SS,SS2,M1(16),M2(16),M3,M4,M5(3),M6(3),
& RNGLU,FACT,HCS,FACTR,RCS,A2,A02,QPE,SQPE,RGRV
C Use the GUPI (Generic User Process Interface) event common block
C as the hard process for HERWIG
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
C--Les Houches Common Block
INTEGER MAXPUP
PARAMETER(MAXPUP=100)
C-----------------------------------------------------------------------
C QCD HEAVY FLAVOUR PRODUCTION: MEAN EVWGT = SIGMA IN NB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,RCS,Z1,Z2,ET,EJ,
& QM2,QPE,FACTR,S,T,U,ST,TU,US,TUS,UST,EN,RN,AF,ASTU,
& AUST,CF,CN,CS,CSTU,CSUT,CTSU,CTUS,HCS,UT,SU,GT,DIST,KK,KK2,
C-----------------------------------------------------------------------
C HIGGS + HEAVY QUARK (BOTTOM & TOP) PRODUCTION (2HDM)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN, HWUALF, HWUAEM, EPS, HCS, RCS,
& DIST, SM, DM, QPE, PF, SQPE, EMSC2, FACTR, S, T3, U4,
& SN2TH, ME2(0:4), MW, XWEIN, PT2MIN, PT2, GQH(0:4), G1, RMMIN,
C-----------------------------------------------------------------------
C ASSOCIATE PRODUCTION W+H- FROM QUARK FUSION (2HDM)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,IHEL
DOUBLE PRECISION EMH,EMHWT,RMW,EMW
DOUBLE PRECISION RMH
C-----------------------------------------------------------------------
C PRODUCTION OF MSSM CHARGED HIGGSES FROM B-QUARK+LIGHT-QUARK FUSION
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,K,L,M,N
INTEGER II,JJ,ITMP
INTEGER IFL,IRES
C =2: infinite mass limit.
C Only top loop included. A factor (alpha_s**3*alpha_W) is extracted
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWHIGB,HWHIGC,HWHIGD,HWHIG5,HWHIG1,HWHIG2,BI(4),
& CI(7),DI(3),EPSI,TAMP(7)
DOUBLE PRECISION S,T,U,EMH2,WTQQ,WTQG,WTGQ,WTGG,EMW2,RNGLU,RNQRK,
C One loop scalar integrals, used in HWHIGJ.
C If NOMASS=.TRUE. use a small mass approx. for particle in loop.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWHIGB,HWUCI2,HWULI2,EPSI,PII
DOUBLE PRECISION S,T,EQ2,EH2,RAT
LOGICAL NOMASS
C One loop scalar integrals, used in HWHIGJ.
C If NOMASS=.TRUE. use a small mass approx. for particle in loop.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWHIGC,HWUCI2,HWULI2,EPSI,PII
DOUBLE PRECISION S,T,EQ2,EH2,RAT,COSH
LOGICAL NOMASS
C One loop scalar integrals, used in HWHIGJ.
C If NOMASS=.TRUE. use a small mass approx. for particle in loop.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE COMPLEX HWHIGD,HWUCI2,HWULI2,EPSI,PII,Z1,Z2
DOUBLE PRECISION S,T,EQ2,EH2,DLS,DLT,DLM,RZ12,DL1,DL2,
& ST,ROOT,XP,XM
C--------------------------------------------------------------------------
C LEPTOPRODUCTION OF MS(SM) HIGGSES IN ASSOCIATION WITH HEAVY QUARK PAIRS
C--------------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER JHIGGS
INTEGER I,L,M,N,NN
INTEGER IH,IQ,JQ,IIQ,JJQ
C-----------------------------------------------------------------------
C DRELL-YAN 2 PARTON -> 2 HIGGS PAIR (2HDM)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN, HWUAEM, EPS, HCS, RCS, DIST, S, PF, QPE,
& FACTR, SN2TH, MZ, MW, MNN(2,2), MCC(2), MCN(3), EMSC2, GW2, GZ2,
& GHH(4), XWEIN, S2W, ECM_MAX, X(3), XL(3),
C Adapted from the program of U. Baur and E.W.N. Glover
C See: Nucl. Phys. B339 (1990) 38
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWUAEM,EPS,RCS,EMH,EMHWT,
& EMHTMP,BR,CV,CA,EMH2,ET,EJ,PT,EMT,EMAX,YMAX,YHINF,YHSUP,EXYH,
& YMIN,YJINF,YJSUP,EXYJ,S,T,U,FACT,AMPQQ,AMPQG,AMPGQ,AMPGG,HCS,
C AS SUGGESTED IN M.H.SEYMOUR, PHYS.LETT.B354(1995)409.
C ENDIF
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRUNI,EM,WEIGHT,EMH,DIF,FUN,THETA,T,EMHLST,W0,
& W1,EMM,GAMEM,T0,TMIN,TMAX,THEMIN,THEMAX,ZMIN,ZMAX,Z,F,GAMOFS
INTEGER I
C-----------------------------------------------------------------------
C PRODUCTION OF MSSM HIGGSES IN ASSOCIATION WITH HEAVY QUARK PAIRS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER JHIGGS
INTEGER I,J,K,L,M,N
INTEGER IS,IH,IQ,JQ,IIQ,JJQ,IQMIN,IQMAX,IGG,IQQ
C HIGGS PRODUCTION VIA GLUON OR QUARK FUSION
C MEAN EVWGT = HIGGS PRODN C-S * BRANCHING FRACTION IN NB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUALF,HWHIGT,HWRGEN,HWUSQR,HWUAEM,BRHIGQ,EMH,
& CSFAC(13),EVSUM(13),EMFAC,CV,CA,BR,RWGT,E1,E2,EMQ,GFACTR,RQM(6)
INTEGER IDEC,I,J,ID1,ID2
C WARNING: THIS IS A FACTOR OF 3 GREATER THAN EHLQ'S ETA FUNCTION
C PARITY=+1 FOR SCALAR AND -1 FOR PSEUDOSCALAR
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWHIGT,RATIO,RAT2,EMH,FREAL,FIMAG,ETALOG,AIREAL,
& AIIMAG
INTEGER I,J,K,L
C MSSM NEUTRAL HIGGS PRODUCTION IN ASSOCIATION WITH GAUGE BOSON
C--BRW fix 27/8/04: corrected off-shell gauge boson mass dependence
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,K,L,M,N
INTEGER IV,IDEC
INTEGER ID1,ID2
C HIGGS PRODUCTION VIA W/Z BOSON FUSION
C MEAN EVWGT = HIGGS PRODN C-S * BRANCHING FRACTION IN NB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWULDO,HWRUNI,HWRGEN,HWUAEM,K1MAX2,K1MIN2,K12,
& K2MAX2,K2MIN2,K22,EMW2,EMW,ROOTS,EMH2,EMH,ROOTS2,P1,PHI1,PHI2,
& COSPHI,COSTH1,SINTH1,COSTH2,SINTH2,P2,WEIGHT,TAU,TAULN,CSFAC,
C
C MEAN EVWGT = CROSS-SECTION (IN NB) * HIGGS BRANCHING FRACTION
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUAEM,HWHIGY,HWRUNI,HWRGEN,HWULDO,EMZ,CVE,CAE,
& POL1,POL2,CE1,CE2,CE3,PMAX,EMZ2,S,B,FACTR,EMH,EMFAC,EMH2,A,XP,
& CV,CA,BRHIGQ,BR,X1,X2,FAC1,FAC2,XPP,XPPSQ,COEF,X,XSQ,PROB,C1,C2,
C-----------------------------------------------------------------------
C PRODUCTION OF MSSM HIGGS PAIRS IN L+L- (L=E,MU)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN, HWUAEM, HCS, RCS, S, PF, QPE,
& FACTR, SN2TH, MZ, MNN(2), MCC, EMSC2, GZ2,
& GHH(4), XWEIN, S2W, X(1), XL(1),
C-----------------------------------------------------------------------
C PRODUCTION OF MSSM HIGGSES IN ASSOCIATION WITH B,T-SQUARK PAIRS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
COMMON/SQSQH/JHIGGS,ILBL,JH,IF1MIN,IF1MAX,IF2MIN,IF2MAX
INTEGER JHIGGS,ILBL,JH,IF1MIN,IF1MAX,IF2MIN,IF2MAX
INTEGER I,J,K,L,M,N
C-----------------------------------------------------------------------
C QQD direct photon pair production: mean EVWGT = sigma in nb
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWHPPB,EPS,RCS,ET,EJ,KK,KK2,
& YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,FACTR,RS,S,T,U,CSTU,TQSQ,
& DSTU,HCS
C-----------------------------------------------------------------------
C QCD DIRECT PHOTON + JET PRODUCTION: MEAN EVWGT = SIGMA IN NB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,HWHPPB,EPS,RCS,ET,EJ,KK,KK2,
& YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,FACTR,FACTF,RS,S,T,U,CF,
& AF,CSTU,CTSU,CUST,DSTU,HCS,TQCH
C point-like photon/QCD heavy flavour single excitation, using exact
C massive lightcone kinematics, mean EVWGT = sigma in nb.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,PP1,PP2,QM2,FACTR,
& PT,PJ,PT2,PTM,EXY,T,CC,EXY2,S,U,C,SIGE,HCS,RCS
INTEGER IQ1,IQ2,ID1,ID2,IHAD1,IHAD2
C Point-like photon/gluon heavy flavour pair production, with
C exact lightcone massive kinematics, mean EVWGT = sigma in nb.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRUNI,HWUALF,EPS,PP1,PP2,QM2,FACTR,ET,EJ,ET2,
& EXY,EXY2,S,T,U,C
INTEGER IQ1,IHAD1,IHAD2
C See M. Benayoun, et al., Nucl. Phys. B282 (1987) 653 for details.
C mean EVWGT = sigma in nb
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,PP1,PP2,ET,EJ,EXY,EXY2,
& FACT,FACTR,S,T,U,REDS,DELT(3,3),C1STU,C3STU,HCS,RCS,CMIX,SMIX,
& C1WVFN,FPI,FETA8,FETA1,FRHO,FPHI8,FPHI1,FPI2,FETA2(3),FETAP2(3),
C-----------------------------------------------------------------------
C point-like photon/QCD di-jet production: mean EVWGT = sigma in nb
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,EPS,RCS,PP1,PP2,ET,EJ,
& EXY,EXY2,FACTR,RS,S,T,U,CSTU,CTSU,HCS
INTEGER ID1,ID2,ID3,ID4,IHAD1,IHAD2
C Compton scattering of point-like photon and (anti)quark
C mean EVWGT = sigma in nb
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,EPS,RCS,PP1,PP2,ET,EJ,EXY,EXY2,
& FACTR,S,T,U,CTSU,HCS
INTEGER ID1,ID2,IHAD1,IHAD2
C-----------------------------------------------------------------------
C QCD HARD 2->2 PROCESSES: MEAN EVWGT = SIGMA IN NB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUALF,RS,EPS,HF,RCS,Z1,Z2,ET,EJ,
& FACTR,S,T,U,ST,TU,US,STU,TUS,UST,EN,RN,GFLA,AF,ASTU,ASUT,AUST,
& BF,BSTU,BSUT,BUST,BUTS,CF,CSTU,CSUT,CTSU,CTUS,DF,DSTU,DTSU,DUTS,
C-----------------------------------------------------------------------
C IDENTIFIES HARD SUBPROCESS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER HWRINT,ID3,ID4,IPERM,IHPR,ND3
EXTERNAL HWRINT
IHPRO=IHPR
C HARD PROCESS: GAMGAM --> QQBAR/LLBAR/W+W-
C MEAN EVENT WEIGHT = CROSS-SECTION IN NB AFTER CUTS ON PT
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION RCS,HCS,RS,S,EMSQ,BE,TMIN,TMAX,T,U,FACTR,Q,CFAC,
$ HWRGEN
INTEGER IHAD1,IHAD2,HQ,ID3,ID4,I1,I2
C-----------------------------------------------------------------------
C Subroutine for 2 parton -> 2 parton via UDD resonant squarks
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HCS,S,RCS,HWRGEN,MQ1,MQ2,TAU,LOWTLM,UPPTLM,RTAB,
& SQSH,MATELM,SCF(12),CHANPB(2),HWRUNI,PCM,MIX(12),
& ME(2,3,3,3,3),WD,MS(12),SWD(12),RAND,TAUA,
C Subroutine for 2 parton -> parton SUSY particle via UDD resonant
C squarks.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HCS,S,RCS,HWRGEN,ME(4),CW,MER(6),MZ,TAU,TAUA,
& TAUB,LOWTLM,UPPTLM,HWRUNI,SH,SQSH,SCF(12),MW2,
& LAMC(3),CHANPB(2),PCM,ECM,RAND,MEN(7,6,3,3),
C SUSY E+E- --> SM PARTICLES VIA RPV
C MODIFIED TO INCLUDE BEAM POLARIZATION EFFECTS BY PETER RICHARDSON
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,HWUMBW,HWUAEM,HCS,RCS,FACA,
& S,T,PCM,MQ1,MQ2,SP,TP,TPZ,TPN,TPN2,MSL2(3),MZ,
& MZ2,MSU2(3,2),MWD(3),GL,GR,GLP,GRP,EC,EE,THTMIN,
C IDENTIFY THE REMNANTS OF THE HARD SCATTERING
C AND BREAK THEIR COLOUR CONNECTION IF NECESSARY
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION PCL(5),
$ P1P2,P1SQ,P2SQ,S,M1SQ,M2SQ,TMP1,TMP2,A,B,C,D,PTOT(4),HWULDO
INTEGER IBEAM,ITARG,IHEP,NTEMP,I,ICOL,IANT
C-----------------------------------------------------------------------
C SUSY E+E- RPV PRODUCTION
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
IF(IPROC.GE.800.AND.IPROC.LE.850) THEN
CALL HWHRES
ELSEIF(IPROC.GE.860.AND.IPROC.LT.890) THEN
C SUSY E+E- --> RPV SINGLE SPARTICLE PRODUCTION
C POLARZATION EFFECTS ADDED 5/4/02 BY PETER RICHARDSON
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWUAEM,HWRUNI,HWUPCM,HWUMBW,HCS,RCS,FACA,
& FACB,FACC,FACD,FACE,M1(4,4),M2(2,4),M3(8,2),
& MW,MZ,MSCL(2,2),MSCL2(2,2),MZ2,MSL2,MSR2,MSNU2,
C slepton mass and mass*width added to save statement to
C avoid problems with Linux by Peter Richardson
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HCS,S,RCS,HWRGEN,FAC,ECM,TH,PCM,CFAC,CHANPB,SH,
& TAU,TAUA,TAUB,LOWTLM,UPPTLM,HWRUNI,MSL(12),
& SQSH,MET(2),SCF(12),MIX(12),ME(4,3,3,3,3,2),
C-----------------------------------------------------------------------
C Subroutine for 2 parton -> sparticle + X via LQD
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HCS,A(6,12),B(6,12),S,RCS,HWRGEN,CW,FAC2,EC,ME2,
& MW,G,TAU,TAUA,TAUB,LOWTLM,UPPTLM,HWRUNI,SW,SQSH,LC,
& SH,MSL(12),MSU(12),MST(6),C(2,6,12),D(2,6,12),UH,
C-----------------------------------------------------------------------
C Subroutine for all hadron-hadron Rparity violating processes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
IF(MOD(IPROC,10000).GE.4000.AND.MOD(IPROC,10000).LT.4060) THEN
C--SINGLE SPARTICLE VIA LQD
CALL HWHRLS
C-----------------------------------------------------------------------
C IDENTIDY HARD R-PARITY VIOLATING PROCESS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ID3, ID4, R4, IPERM,TYPE,ID1,ID2,NEUTD1(8),SLEPID(8),
& NEUTD2(13),SQUID(6),SGN,HWUANT,SQUID2(12),SLPID2(12),
& GAGID1(6),GAGID2(8)
C SCATTERS ABOVE PTMIN WITH PROBABILITY 1/(M+1)
C PTJIM IS THE MINIMUM TRANSVERSE MOMENTUM FOR ADDITIONAL SCATTERS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRGET,HWRSET,WGT,PBOOST(5),RBOOST(3,3),
$ WJMAX,PT,PTJIM,DUMMY,HWUPCM
INTEGER IHEP,IBM,ITG,IBMN,ITGN,IBMT,ITGT,I,REPORT,NHARD,
C CALL HWHQCD, AND SWAP BACK
C WGT IS THE OUTPUT EVENT WEIGHT
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION WGT,PTJIM,XMIN,XMAX,XPOW,
$ TMPXMN,TMPXMX,TMPXPW,TMPWGT
LOGICAL FIRST
C TREATS ALL PARTONS ON EQUAL FOOTING WITH HWHSNM(ID1,ID2,S,T)
C PROVIDING THE MATRIX ELEMENT SQUARED FOR PARTON TYPES ID1 AND ID2
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ID1,ID2
DOUBLE PRECISION HWRGEN,HWRUNI,HWHSNM,EPS,RCS,ET,EJ,KK,KK2,
& YJ1INF,YJ1SUP,Z1,YJ2INF,YJ2SUP,Z2,FACT,S,T,U,HCS
C INCLUDED FOR IDENTICAL QUARKS AND LIKEWISE S-CHANNEL CONTRIBUTION
C FOR IDENTICAL QUARK-ANTIQUARK PAIRS.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWHSNM,HWUAEM,HWUALF,S,T,ASQ,AINU,AINS,Y,SOLD,
$ TOLD,QQ(13,13),ZETA3
INTEGER ID1,ID2
C-----------------------------------------------------------------------
C Calculates the spin correlations for the hard process
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ZI,S,D,ME(2,2,2,2,NCFMAX),MED(2,2,2,2),F3(2,2,8),
C section f fbar --> gauge boson --> fermion fermion
C This diagram 1 from DAMTP-2001-83 with opposite sign of P4
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section f fbar ---> fermion fermion via t-channel scalar exchange
C This diagram 2 from DAMTP-2001-83 with opposite sign of P4
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f fbar ---> fermion fermion via u-channel scalar exchange
C This diagram 3 from DAMTP-2001-83 with opposite sign of P4
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,
C section f fbar --> gauge boson --> fermion antifermion
C This diagram 1 from DAMTP-2001-83
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C N.B. a gauge choice has been made to simplify the triple gluon vertex
C This diagram 4 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C N.B. a gauge choice has been made to simplify the triple gluon vertex
C This diagram 5 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C N.B. a gauge choice has been made to simplify the triple gluon vertex
C This diagram 6 from DAMTP-2001-83 with the gauge choice L1=2 L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section quark gluon --> fermion sfermion
C This diagram 7 from DAMTP-2001-83 with the gauge choice L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section antiquark gluon --> fermion antisfermion
C This diagram 10 from DAMTP-2001-83 with the gauge choice L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section quark gluon --> fermion antisfermion (s-channel quark)
C This is diagram 8 from DAMTP-2001-83 with the gauge choice L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section quark gluon --> fermion antisfermion (s-channel quark)
C This is diagram 11 from DAMTP-2001-83 with the gauge choice L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section quark gluon --> fermion antisfermion (s-channel quark)
C This is diagram 9 from DAMTP-2001-83 with the gauge choice L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section quark gluon --> fermion antisfermion (s-channel quark)
C This is diagram 12 from DAMTP-2001-83 with the gauge choice L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C This diagram 4 from DAMTP-2001-83 with opposite helicity for 4
C and gauge choice L1=2 L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
C This diagram 5 from DAMTP-2001-83 with opposite helicity for 4
C and gauge choice L1=2 L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
C This diagram 6 from DAMTP-2001-83 with opposite helicity for 4
C and gauge choice L1=2 L2=1
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI
C section fermion fermion --> fermion fermion (t-channel boson)
C This diagram 13 from DAMTP-2001-83
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
C section fermion antifermion --> fermion antifermion (t-channel boson)
C This diagram 14 from DAMTP-2001-83
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
C section antifermion fermion --> antifermion fermion (t-channel boson)
C This diagram 15 from DAMTP-2001-83
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
C section antifermion fermion --> antifermion fermion (t-channel boson)
C This diagram 16 from DAMTP-2001-83
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE,ZI,
C section f fbar ---> f fbar via s-channel scalar exchange
C This is diagram 1 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f fbar ---> f fbar via t-channel scalar exchange
C This is diagram 2 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f fbar ---> fermion fermion via t-channel scalar exchange
C This is diagram 3 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f fbar ---> f f via s-channel scalar exchange
C This is diagram 4 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f fbar ---> f f via u-channel scalar exchange
C This is diagram 5 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f fbar ---> f f via t-channel scalar exchange
C This is diagram 6 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f f ---> f fbar via s-channel scalar exchange
C This is diagram 7 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f f ---> f fbar via t-channel scalar exchange
C This is diagram 8 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f f ---> f fbar via u-channel scalar exchange
C This is diagram 9 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section fbar fbar ---> f f via s-channel scalar exchange
C This is diagram 10 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section fbar fbar ---> f f via t-channel scalar exchange
C This is diagram 11 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section fbar fbar ---> f f via u-channel scalar exchange
C This is diagram 12 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section f f ---> f f via s-channel scalar exchange
C This is diagram 13 from RPV
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C section fbar fbar ---> fbar fbar via t-channel scalar exchange
C This is diagram 14 from RPV notes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NDIAHD
PARAMETER(NDIAHD=10)
DOUBLE COMPLEX ME(2,2,2,2),S,D,F3(2,2,8),F4(2,2,8),PRE
C-----------------------------------------------------------------------
C LL(BAR) -> GAUGINOS (including beam polarization)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWHSS2, S, T, U, M3, M4, SGN
DOUBLE COMPLEX CLL, CLR, CRL, CRR
HWHSS2 =
C SUSY 2 PARTON -> 2 GAUGINOS PROCESSES (1 - 3)
C -> GAUGINO + SPARTON PROCESSES (4 - 7)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN, HWUALF, HWUAEM, EPS, HCS, RCS, DIST,
& ML(6), ML2(6), MR(6), MR2(6), MCH(2), MCH2(2), MNU(4), MNU2(4),
& MSQK, MG, MG2, SM, DM, DAB, QPE, SGN, PF, SQPE, EMSC2,
C-----------------------------------------------------------------------
C SUSY 2 PARTON -> 2 SLEPTON PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN, HWUAEM, EPS, HCS, RCS, DIST, S, PF, QPE,
& FACTR, SN2TH, MZ, MW, ME2(2,2,6,2), ME2W(2,3), EMSC2, GW2
INTEGER IQ, IQ1, IQ2, ID1, ID2, IL, IL1, IL2, I, J
C-----------------------------------------------------------------------
C SUSY HARD 2 PARTON -> 2 SPARTON PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN, HWUALF, EPS, HCS, RCS, DIST, NC, NC2,
& NC2C, ML2(6), ML4(6), MR2(6), MR4(6), MG2, SM, DM, QPE,
& SQPE, FACTR, AFAC, AF, BONE, CFAC, CFC2, CFC3, CONE,
C-----------------------------------------------------------------------
C SUSY HARD 2 PARTON -> 2 SPARTON/GAUGINO/SLEPTON PROCESSES
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION SAVWT(3),RANWT,HWRGEN,HWRUNI,Z1,Z2,ET,EJ,
& QPE,S,T,U,KK,KK2,YJ1INF,YJ1SUP,YJ2INF,YJ2SUP,SVEMSC
INTEGER ISP
C-----------------------------------------------------------------------
C IDENTIFIES HARD SUSY SUBPROCESS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ID3, R3, ID4, R4, IPERM, IHPR, SSL
PARAMETER (SSL = 400)
IHPRO = 3000 + IHPR
C USES CROSS-SECTIONS OF EHLQ FOR ANNIHILATION AND COMPTON SCATTERING
C IHPRO=0 FOR BOTH, 1 FOR ANNIHILATION, AND 2 FOR COMPTON.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,DISFAC(2,12,2),EMV2,DISMAX,S,T,U,
& SHAT,THAT,UHAT,Z,HWUALF,PT,EMT,GFACTR,SIGANN,SIGCOM(2),CSFAC,ET,
& EJ,YMIN,YMAX,VYMIN,VYMAX,EMAX,CV,CA,BR,EMV,GAMV,HWUAEM,TMIN,TMAX
C Master subroutine for all vector boson + 2 jet processes
C Currently implemented qqbar Z only
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,K,IDBS,IPRC,IDP(6),ORD,IB,ICMF,IHEP,IFLOW,IZ,IBRAD,
& ICOL(5),IDZ,IQ
DOUBLE PRECISION HWRGEN,HWRUNI,XMASS,PLAB,PRW,PCM,HWUAEM,BR,FLUX,
C---DSDCOS HOLDS THE CROSS-SECTIONS FOR THE PROCESSES LISTED ABOVE
C (1-8) ARE WITH B FROM BEAM 1, (9-16) ARE WITH B FROM BEAM 2.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,DSDCOS(16),EMT2,EMT,EMW2,EMW,
& CMFMIN,TAUMIN,TAUMLN,S,T,U,ROOTS,DSMAX
INTEGER HWRINT,IDHWEX(2,16),I
C W+/- PRODUCTION AND DECAY VIA DRELL-YAN PROCESS
C MEAN EVWGT IS SIG(W+/-)*(BRANCHING FRACTION) IN NB
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRUNI,HWUPCM,PRAN,PROB,COEF,CSFAC,EMW,
& FTQK,PTOP,ETOP,EBOT,PMAX,FHAD,FTOT,BRAF,FLEP,TMIN,HWUAEM,TMAX
INTEGER HWRINT,ICH,IC,IL,ID,IDEC,JDEC,IWP(2,16)
C-----------------------------------------------------------------------
SUBROUTINE HWICHK
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
IF(RMASS(1).LT.0.1.OR.RMASS(1).GT.1.0.OR.
& FMRS(1,1,20,1).LT.0.1.OR.FMRS(1,1,20,1).GT.1.0) THEN
STOP 'Block data hwudat not loaded, stop execution'
C When reading in if IME =1 matrix element codes >= 100 are accepted
C 0 are set zero.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IUNIT,IOPT,IME,JUNIT,I,J,K,L,IDKY,ITMP(5),IDUM
CHARACTER*8 CDK(NMXDKS),CDKPRD(5,NMXDKS),CDUM
JUNIT=ABS(IUNIT)
C-----------------------------------------------------------------------
C SETS INPUT PARAMETERS
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION FAC,ANGLE
INTEGER I,J,N,L
CHARACTER*28 TITLE
C Use the GUPI (Generic User Process Interface) run common block
C to initialise HERWIG -- Initialization for Les Houches interface
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER MAXPUP
PARAMETER(MAXPUP=100)
INTEGER IDBMUP,PDFGUP,PDFSUP,IDWTUP,NPRUP,LPRUP
C-----------------------------------------------------------------------
C Subroutine to merge Higgs WW/ZZ decay modes for four body ME
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IH,I,NMODE,J,K
LOGICAL REMOVE
DOUBLE PRECISION BR
C default values
C IOPT = 2 optimises the weights for the process selected
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,IPRC,ICH,IOPT,ISTP,IWGT,IFER,IANT,IGAU,IQRK
LOGICAL CALLED,TEV,LHC
DOUBLE PRECISION CHNPST(IMAXCH,IMAXOP),D(IMAXOP),CHWGTS(IMAXCH),
C Calculates the couplings for the SUSY decays for spin correlations
C and 3/4 body matrix elements
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUALF,PRE,MCHAR(2),QIJPP(4,4),SIJPP(4,4),
& DIJ(2,2),QIJ(2,2),R(4,2),SIJ(2,2)
INTEGER I,J,K,L,IH,IK,IL,IQ
C Initialise all the decay modes for three/four body MEs and spin
C correlations
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,K,NDKYST
C--set the number of two and three body modes to zero
N2MODE = 0
C-----------------------------------------------------------------------
C Initialise the SUSY two body modes for spin correlations
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,IL,IH,L,L1,IM,O(2),II,JJ,III,JJJ,KKK
COMMON /HWSPNC/ SW,CW,TW,E,G,RT,ORT,MW,MZ,GS,AFN,AFG,AFC,OIJ,OIJP,
& OIJPP,HNN,HCC,HNC,HFF,HWW,HZZ,ZAB,HHB
C Initialise the top/SUSY three body decay modes
C gravitino and RPV modes added by Peter Richardson
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,K,L,L1,IL,IQ,IQ1,IQ2,IFR,SIFR,IH,IH1,IM,O(2),II,JJ,
& III,JJJ,KKK
DOUBLE PRECISION SW,CW,TW,E,G,RT,ORT,MW,MZ,GS,AFN(2,12,2,4),
C-----------------------------------------------------------------------
C Initialise the Higgs four body modes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,K,IL,IH,II,JJ
DOUBLE PRECISION COL(2),SW,CW,TW,E,G,RT,ORT,MW,MZ,AFN(2,12,2,4),
& AFG(2,6,2),AFC(2,12,2,2),OIJ(2,4,2),OIJP(2,2,2),OIJPP(2,4,4),
C Reads in SUSY particle properties and decays,
C in format generated by ISAWIG
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,K,IH,IHW,NSSP,NDEC,MDKYS
DOUBLE PRECISION BETAH, WEINCOS,WEINSIN, MW,MZ, RMMAX
DOUBLE PRECISION FTM,FTMUU(4),FTMDD(4),FTMTT(4),FTMBB(4),FTMU,FTMD
C = 8000: CREATES MINIMUM-BIAS EVENT
C SUPPRESSED BY ADDING 10000 TO IPROC
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWREXP,ENFAC,TECM,SECM,SUMM,EMCL,BMP(5),BMR(3,3)
INTEGER HWRINT,NETC,IBT,IDBT,ID1,ID2,ID3,KHEP,LHEP,NTRY,ICMS,
& NPPBAR,MCHT,JCL,JD1,JD2,JD3,ICH,MODC,NCHT,INHEP(2),
C GENERATES CYLINDRICAL PHASE SPACE USING THE METHOD OF JADACH
C RETURNS WITH NCL=0 IF UNSUCCESSFUL
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWREXT,HWRUNG,HWUSQR,TECM,ESS,ALOGS,EPS,SUMX,
& SUMY,PT,PX,PY,PT2,SUMPT2,SUMTM,XIMIN,XIMAX,YY,SUM1,SUM2,SUM3,
& SUM4,EX,FY,DD,DYY,ZZ,E1,TM,SLOP,XI(NMXCL)
C if necessary. The branching ratios of any other IDKTMP decays are
C scaled by (1.-BRTMP)/(1.-BR_OLD)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUPCM,BRTMP,SCALE,EPS
INTEGER IDKTMP,IMETMP,IATMP,IBTMP,ICTMP,IDTMP,IETMP,IDKY,ITMP(5),
& L,I,J,K,JPREV
C-----------------------------------------------------------------------
C Chooses charged multiplicity NCHT at the p-pbar c.m. energy EPPBAR
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWMNBI,HWRGEN,EPPBAR,E0,ALOGS,RK,EK,AVN,SUM,R,
& CUM(500)
INTEGER NCHT,IMAX,I,N
C-----------------------------------------------------------------------
C COMPUTES WEIGHT FOR MINIMUM-BIAS EVENT
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION S,X,Y
INTEGER IDB,IDT,IDBT
IF (IERROR.NE.0) RETURN
C-----------------------------------------------------------------------
C subroutine for radiation in top decays
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP,KK,IPOS,NN,NHEP0,KK1,KK2,JMOH(NMXHEP)
DOUBLE PRECISION HWDPWT
EXTERNAL HWDPWT
C-----------------------------------------------------------------------
C subroutine for radiation in top production
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
C--local variables
INTEGER IMO(10),IFOUND,JMO(2),I,J,K,L,NSTART,NHEPX
C--initialisation
C Gaussian random number, mean A, standard deviation B.
C Generates uncorrelated pairs and throws one of them away.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGAU,HWRGEN,A,B,X,TRASH
INTEGER J
EXTERNAL HWRGEN
*-- Author : F. James, modified by Mike Seymour
*- Split in 3 files by M. Kirsanov. Initial seeds ISEED set in HWUDAT
C-----------------------------------------------------------------------
- FUNCTION HWRGEN(I)
+
+C FUNCTION HWRGEN(I)
C-----------------------------------------------------------------------
C MAIN RANDOM NUMBER GENERATOR
C USES METHOD OF l'Ecuyer, (VIA F.JAMES, COMP PHYS COMM 60(1990)329)
C-----------------------------------------------------------------------
- IMPLICIT NONE
- DOUBLE PRECISION HWRGEN
- COMMON/HWSEED/ISEED(2)
- INTEGER ISEED
- INTEGER I,K,IZ
-C
- K=ISEED(1)/53668
- ISEED(1)=40014*(ISEED(1)-K*53668)-K*12211
- IF (ISEED(1).LT.0) ISEED(1)=ISEED(1)+2147483563
- K=ISEED(2)/52774
- ISEED(2)=40692*(ISEED(2)-K*52774)-K*3791
- IF (ISEED(2).LT.0) ISEED(2)=ISEED(2)+2147483399
- IZ=ISEED(1)-ISEED(2)
- IF (IZ.LT.1) IZ=IZ+2147483562
- HWRGEN=DBLE(IZ)*4.656613001013252D-10
+C IMPLICIT NONE
+C DOUBLE PRECISION HWRGEN
+C COMMON/HWSEED/ISEED(2)
+C INTEGER ISEED
+C INTEGER I,K,IZ
+C
+C K=ISEED(1)/53668
+C ISEED(1)=40014*(ISEED(1)-K*53668)-K*12211
+C IF (ISEED(1).LT.0) ISEED(1)=ISEED(1)+2147483563
+C K=ISEED(2)/52774
+C ISEED(2)=40692*(ISEED(2)-K*52774)-K*3791
+C IF (ISEED(2).LT.0) ISEED(2)=ISEED(2)+2147483399
+C IZ=ISEED(1)-ISEED(2)
+C IF (IZ.LT.1) IZ=IZ+2147483562
+C HWRGEN=DBLE(IZ)*4.656613001013252D-10
C---> (4.656613001013252D-10 = 1.D0/2147483589)
- END
+c END
CDECK ID>, HWRSET.
*CMZ :- -26/04/91 12.42.30 by Federico Carminati
*-- Author : F. James, modified by Mike Seymour
C-----------------------------------------------------------------------
C Generates a random primary IP using a triple Gaussian distribution
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGAU
INTEGER I
EXTERNAL HWRGAU
C-----------------------------------------------------------------------
C DOES BRANCHING OF SPACELIKE PARTON KPAR
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWRGEN,HWRUNI,HWSTAB,HWUALF,HWUTAB,HWSGQQ,
& HWSSUD,XLAST,QNOW,QLST,QP,QMIN,QLAM,QSAV,SMAX,SLST,SNOW,RN,SUDA,
& SUDB,ZZ,ENOW,XI,PMOM,DIST(13),DMIN,X1,X2,REJFAC,OTHXI,OTHZ,QTMP,
C ID1.LT.0 ON RETURN MEANS NO PHASE SPACE
C ID1.EQ.0 ON RETURN FLAGS REJECTED BRANCHINGS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWRGEN,HWUALF,HWUAEM,QP,X,QQ,Z,WQG,WQV,
& WQP,XQV,ZMIN,ZMAX,YMIN,YMAX,DELY,YY,PSUM,EZ,WQN,WR,ZR,WZ,ZZ,AZ,
& PVAL,EY,DIST(13),PROB(13,100),PPHO
C-----------------------------------------------------------------------
C---BRW change 27/8/04: include Frixione's fix to reduce PDFSET calls
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWSGAM,X,SCALE,XOLD,QOLD,XMWN,QSCA,SS,SMIN,S,T,
& TMIN,TMAX,VX,AA,VT,WT,UPV,DNV,SEA,STR,CHM,BTM,TOP,GLU,WX,XQSUM,
& DMIN,TPMIN,TPMAX,DIST(13),G(2),Q0(5),QL(5),F(5),A(6,5),
IF(PARM(1).NE.PARMSAVE.OR.VAL(1).NE.VALSAVE)THEN
PARMSAVE=PARM(1)
VALSAVE=VAL(1)
- CALL PDFSET(PARM,VAL)
+ CALL PDFSET_HERWIG(PARM,VAL)
ENDIF
IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR.
& X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN
IF(PARM(1).NE.PARMSAVE.OR.VAL(1).NE.VALSAVE)THEN
PARMSAVE=PARM(1)
VALSAVE=VAL(1)
- CALL PDFSET(PARM,VAL)
+ CALL PDFSET_HERWIG(PARM,VAL)
ENDIF
IF (X.LT.PDFXMN.AND.PDFWRX(IBEAM,1) .OR.
& X.GT.PDFXMX.AND.PDFWRX(IBEAM,2)) THEN
C GENERATES X VALUES (IF GENEX)
C EVALUATES STRUCTURE FUNCTIONS AND ENFORCES CUTOFFS ON X
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWRUNI,X,QL
INTEGER I,J
LOGICAL GENEX
C CORRECTION TO GLUON STRUCTURE FUNCTION FOR BACKWARD EVOLUTION:
C G->Q-QBAR PART OF FORM FACTOR
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWSGQQ,HWUALF,QSCA,GG
EXTERNAL HWUALF
GG=HWUALF(1,QSCA)**(-ONE/BETAF)
C-----------------------------------------------------------------------
C MRST98 Leading order PDF's central and higher gluon + average
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION X,Q,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU,XMIN,XMAX,
& QSQMIN,QSQMAX,Q2,QQ(NQMRS),XXMRS(NXMRS),G(NPMRS),N0(NPMRS),
& XSAVE,Q2SAVE,XXX,A,B,FAC
C-----------------------------------------------------------------------
C REPLACES SPACELIKE PARTONS BY SPECTATORS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUSQR,EMSQ,EMTR,EPAR,XPAR,QSQ,PCL(5)
INTEGER KHEP,IP,JP,IDH,IDP,ISP,IDSPC,JHEP
EXTERNAL HWUSQR
C-----------------------------------------------------------------------
FUNCTION HWSSUD(I)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWSSUD,HWSGQQ,DMIN,QSCA,XLAST,DIST(13)
INTEGER I,N0,IS,ID
EXTERNAL HWSGQQ
C-----------------------------------------------------------------------
C TRUE FOR VALENCE PARTON ID IN INCOMING HADRON INHAD
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ID,IDHAD
LOGICAL HWSVAL
HWSVAL=.FALSE.
C See R. Kleiss et al.: CERN yellow report 89-08, vol.3 p.129
C Hadronic component from: H. Burkhardt et al.: Z. Phys C43 (89) 497
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUAEM,HWUAER,Q2,EPS,A1,B1,C1,A2,B2,C2,A3,B3,C3,
& A4,B4,C4,AEMPI,EEL2,EMU2,ETAU2,ETOP2,REPIGG,X
LOGICAL FIRST
C WITH 5-FLAVOUR BETA, LAMBDA=QCDL3
C .EQ.3 ONE-LOOP WITH 5-FLAVOUR BETA, LAMBDA=QCDL3
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUALF,SCALE,KAFAC,B3,B4,B5,B6,C3,C4,C5,C6,C35,
& C45,C65,D35,RHO,RAT,RLF,DRH,EPS
INTEGER IOPT,ITN
C-----------------------------------------------------------------------
C Returns the antiparticle of IPART; uses HERWIG numbering
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER HWUANT,IPART,IPDG,IANTI,OLDERR
CHARACTER*8 CDUM
OLDERR=IERROR
C-----------------------------------------------------------------------
C Replaces all &'s in TXNAME by backslashes
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J,L
CHARACTER*1 Z
Z=CHAR(92)
C-----------------------------------------------------------------------
C PRINTS OUT DATA ON PARTON SHOWER
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I,J
IF (PRVTX) THEN
WRITE(6,10) INHAD,XFACT
C BOOST THE ENTIRE EVENT RECORD TO (IOPT=1) OR FROM (IOPT=0) ITS
C CENTRE-OF-MASS FRAME, WITH INCOMING HADRONS ON Z-AXIS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION PBOOST(5),RBOOST(3,3)
INTEGER IOPT,IHEP,BOOSTD,IHAD
SAVE BOOSTD,PBOOST,RBOOST
C I,J= 1- 6: d,u,s,c,b,t
C =11-16: e,nu_e,mu,nu_mu,tau,nu_tau
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION QSQ,CLF(7),POL1,POL2,QIF,VI,AI,VF,AF,PG,DQM,PMW,
& DEN,XRE,XIM,XSQ,VI2,AI2,VF2,AF2,PG2,PG12,DQM2,PMW2,DEN2,XRE2,
& XIM2,XSQ2,XRE12,XIM12
C QBORQQ(I) = .TRUE. if it is an antiquark or a diquark
C And stores the particle decay tables: call HWUDPR to print them
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
COMMON/HWSEED/ISEED(2)
INTEGER ISEED
INTEGER NLAST,NNEXT,NLEFT,NREST,I,J,MMWIDE,MMLONG,MMHOFF,MMVOFF
C Modified 16/01/01 by BRW to force particle on mass shell if
C p^2-m^2 < 10^-10 GeV^2 (rounding errors)
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,PMOM(4),DISP(4),PMOM2,SCALE,OFFSH
INTEGER ID
EXTERNAL HWRGEN
C modes are not the charge conjugates of the particles.
C (N.B. CP violation permits this).
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUPCM,HWUAEM,HWUALF,BRSUM,EPS,SCALE,
& BRTMP(NMXDKS),FN,X,W,Q,FAC
INTEGER HWUANT,I,IDKY,LAST,LTMP(NMXMOD),J,L,K,M,N,INDX(NMXMOD),
C If (PRNWEB) HTML to the file HW_decays/index.html
C /PART0000001.html etc.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,IUNITT,IUNTW1,IUNTW2,I,NM,J,K,
& L,M
CHARACTER*1 Z
C IF N>0 DELETE THE N ENTRIES IN IEDT FROM EVENT RECORD
C IF N<0 INSERT LINES AFTER THE -N ENTRIES IN IEDT
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER N,IEDT(*),IMAP(0:NMXHEP),IHEP,I,J,I1,I2
COMMON /HWUMAP/IMAP
C---MOVE ENTRIES AND CALCULATE MAPPING OF POINTERS
C Loads cross-section coefficients, for kinematically open channels,
C in llbar-->qqbar; lepton label IL=1-6: e,nu_e,mu,nu_mu,tau,nu_tau.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION Q2
INTEGER IL,JL,IQ
Q2=EMSCA**2
C MOVE A BLOCK OF ENTRIES IN THE EVENT RECORD
C N ENTRIES IN HEPEVT STARTING AT IFROM ARE MOVED TO AFTER ITO
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER N,IFROM,ITO,IMAP(0:NMXHEP),LFROM,LTO,I,IEDT(NMXHEP),IHEP,
$ JHEP,KHEP
COMMON /HWUMAP/IMAP
C the subdirectory HW_decays/
C ******* gives the event number 0000001 etc.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER MMWIDE,MMLONG,MMHOFF,MMVOFF,I,IST,IS,ID,MS,J,K,IUNITW,
& IUNITT
CHARACTER*1 Z
C Subroutine to handle termination of HERWIG if reaches end of event
C file
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
C--reset the number of events to the correct value
NEVHEP = NEVHEP-1
C--output information on the events
C FINALISES THE EVENT BY UNDOING THE LORENTZ BOOST IF THERE WAS ONE,
C CHECKING FOR ERRORS, AND PRINTING
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IHEP
LOGICAL CALLED
COMMON/HWDBUG/CALLED
IF (NEGWTS.AND.EVWGT.LT.ZERO) NNEGEV=NNEGEV-1
C---PRINT FIRST MAXPR EVENTS
ELSEIF (NEVHEP.LE.MAXPR) THEN
+ write(6,*) "Calling HWUEPR", NEVHEP, MAXPR
CALL HWUEPR
END IF
END
C IOPT= 2 GIVEN IWIG, RETURNS IPDG AND NWIG
C IOPT= 3 GIVEN NWIG, RETURNS IPDG AND IWIG
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IOPT,IPDG,IWIG,I
CHARACTER*8 NWIG
IF (IOPT.EQ.1) THEN
C COMPUTES CONSTANTS AND LOOKUP TABLES
C---BRW change 27/8/04: include Frixione's fix to reduce PDFSET calls
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWBVMC,HWUALF,HWUPCM,XMIN,XMAX,XPOW,QR,DQKWT,
& UQKWT,SQKWT,DIQWT,QMAX,PMAX,PTLIM,ETLIM,PGS,PTELM,X,QSCA,UPV,DNV,
& USEA,DSEA,STR,CHM,BTM,TOP,GLU,VAL(20),CLMXPW,RCLPOW,TEST,RPM(2)
WRITE (6,'(10X,A)') 'THEIR WISHES. SSph NO LONGER WORKS'
STOP
ELSE
- CALL PDFSET(PARM,VAL)
+ CALL PDFSET_HERWIG(PARM,VAL)
CALL STRUCTM(X,QSCA,UPV,DNV,USEA,DSEA,STR,CHM,BTM,TOP,GLU)
ENDIF
ENDIF
C-----------------------------------------------------------------------
C INITIALISES AN EVENT
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWRGEN,HWRGET,DUMMY
REAL TL
LOGICAL CALLED,HWRLOG
C CHOOSES MASS ACCORDING TO BREIT-WIGNER DISTRIBUTION
C--BRW fix 27/8/04: changed from mass to mass-squared BW formula
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUMBW,HWRGEN,WMX,TAU,GAM,T,TM
INTEGER ID
C--WMX IS MAX NUMBER OF WIDTHS FROM NOMINAL MASS
C Prints contents of the GUPI (Generic User Process Interface)
C common block HEPEUP
C----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER MAXNUP
PARAMETER (MAXNUP=500)
INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
C Finds iso-flavour hadrons and creates pointers for cluster decays.
C Sets CLDKWT(K) =(2J+1) spin weight normalizing largest value to 1.
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER NMXTMP
PARAMETER (NMXTMP=20)
DOUBLE PRECISION EPS,WTMX,REMMN,RWTMX,WTMP,RESTMP(91),WTMX2,
C RUNNING QUARK MASSES (MSBAR, 2-LOOP, 5 FLAVOUR, NO THRESHOLDS)
C ASSUMING RMASS(IQ) IS POLE MASS
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
DOUBLE PRECISION HWUALF,SCALE,ALFAS,P0,C1,CC,MHAT(6),RQM(6)
INTEGER IQ
LOGICAL FIRST
C-----------------------------------------------------------------------
C Subroutine to output the contents of the spin common block
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER I
C--write out the header
WRITE(6,1000)
C-----------------------------------------------------------------------
C MAKES PARTICLE TYPE 'NAME' STABLE
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER IPDG,IWIG
CHARACTER*8 NAME
CALL HWUIDT(3,IPDG,IWIG,NAME)
C 400-499 FATAL, DUMP EVENT & STOP DEAD
C 500- FATAL, STOP DEAD WITH NO DUMP
C-----------------------------------------------------------------------
- INCLUDE 'HERWIG65.INC'
+ INCLUDE 'herwig65.inc'
INTEGER ICODE
CHARACTER*6 SUBRTN
IF (ICODE.GE.0) WRITE (6,10) SUBRTN,ICODE