X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=HERWIG%2Fherwig6510.f;h=3b38f3402c39d74ba90d162417f81dbf6b49b8e2;hp=fd8b4743a83084a078bddc61f90f180e32893977;hb=f0667f1e9f6c8e42a086a1cf99b12276d08995f4;hpb=657679557603e320230da5df469d688018515c79 diff --git a/HERWIG/herwig6510.f b/HERWIG/herwig6510.f index fd8b4743a83..3b38f3402c3 100644 --- a/HERWIG/herwig6510.f +++ b/HERWIG/herwig6510.f @@ -152,7 +152,7 @@ C used by TAUOLA 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 @@ -256,7 +256,7 @@ C----------------------------------------------------------------------- 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 @@ -330,7 +330,7 @@ C MODIFIED 11/01/01 BY PR FOR SPIN CORRELATIONS(PROBLEM WITH ORDER 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 @@ -508,7 +508,7 @@ C 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 @@ -779,7 +779,7 @@ C FILL MISSING AREA OF DIS PHASE-SPACE WITH 2+1-JET EVENTS 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), @@ -1252,7 +1252,7 @@ C----------------------------------------------------------------------- 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, @@ -1862,7 +1862,7 @@ C----------------------------------------------------------------------- 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 @@ -1979,7 +1979,7 @@ C BRANCHING GENERATOR WITH INTERFERING GLUONS 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 @@ -2194,7 +2194,7 @@ C Makes the colour connections and performs the parton shower 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 @@ -2386,7 +2386,7 @@ C----------------------------------------------------------------------- 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, @@ -2956,7 +2956,7 @@ C 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 @@ -3070,7 +3070,7 @@ C----------------------------------------------------------------------- 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, @@ -3537,7 +3537,7 @@ C----------------------------------------------------------------------- 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, @@ -4233,7 +4233,7 @@ C----------------------------------------------------------------------- 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 @@ -4275,7 +4275,7 @@ C----------------------------------------------------------------------- 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). @@ -4494,7 +4494,7 @@ C jet started by parton no.2 interference partner 1 and spin density 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 @@ -4629,7 +4629,7 @@ C Constructs appropriate spin density/decay matrix for parton 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 @@ -4736,7 +4736,7 @@ C----------------------------------------------------------------------- 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 @@ -4900,7 +4900,7 @@ C LOGARITHMIC PART OF INTEGRAND IN EXPONENT OF SUDAKOV FORM FACTOR. 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 @@ -4982,7 +4982,7 @@ C Includes azimuthal angular correlations between branching planes 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 @@ -5136,7 +5136,7 @@ CDECK ID>, HWBTOP. 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, @@ -5297,7 +5297,7 @@ C----------------------------------------------------------------------- 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 @@ -5319,7 +5319,7 @@ C----------------------------------------------------------------------- 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), @@ -5471,7 +5471,7 @@ C----------------------------------------------------------------------- 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, @@ -5704,7 +5704,7 @@ C----------------------------------------------------------------------- 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) @@ -5769,7 +5769,7 @@ C----------------------------------------------------------------------- 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 @@ -5826,7 +5826,7 @@ C----------------------------------------------------------------------- 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 @@ -5978,7 +5978,7 @@ C----------------------------------------------------------------------- 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 @@ -6053,7 +6053,7 @@ C Modified by IGK to include BRW's colour rearrangement and 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 @@ -6394,7 +6394,7 @@ C----------------------------------------------------------------------- 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 @@ -6555,7 +6555,7 @@ C (IN SPECIAL CLUSTER CODE - SEE HWCFLA) 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, @@ -6761,7 +6761,7 @@ C----------------------------------------------------------------------- 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) @@ -6887,7 +6887,7 @@ C----------------------------------------------------------------------- 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, @@ -7115,7 +7115,7 @@ C----------------------------------------------------------------------- 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, @@ -7395,7 +7395,7 @@ C----------------------------------------------------------------------- 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) @@ -7473,7 +7473,7 @@ C----------------------------------------------------------------------- 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) @@ -7523,7 +7523,7 @@ C----------------------------------------------------------------------- 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) @@ -7572,7 +7572,7 @@ C----------------------------------------------------------------------- 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) @@ -7621,7 +7621,7 @@ C----------------------------------------------------------------------- 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) @@ -7699,7 +7699,7 @@ C----------------------------------------------------------------------- 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) @@ -7771,7 +7771,7 @@ C----------------------------------------------------------------------- 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, @@ -7814,7 +7814,7 @@ C----------------------------------------------------------------------- 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) @@ -7863,7 +7863,7 @@ C----------------------------------------------------------------------- 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) @@ -7912,7 +7912,7 @@ C----------------------------------------------------------------------- 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) @@ -7961,7 +7961,7 @@ C----------------------------------------------------------------------- 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) @@ -8010,7 +8010,7 @@ C----------------------------------------------------------------------- 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) @@ -8059,7 +8059,7 @@ C----------------------------------------------------------------------- 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) @@ -8108,7 +8108,7 @@ C----------------------------------------------------------------------- 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) @@ -8157,7 +8157,7 @@ C----------------------------------------------------------------------- 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) @@ -8206,7 +8206,7 @@ C----------------------------------------------------------------------- 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) @@ -8255,7 +8255,7 @@ C----------------------------------------------------------------------- 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) @@ -8304,7 +8304,7 @@ C----------------------------------------------------------------------- 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) @@ -8414,7 +8414,7 @@ C----------------------------------------------------------------------- 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, @@ -8644,7 +8644,7 @@ C IF BOSON CAME FROM HIGGS DECAY, GIVE BOTH THE SAME HELICITY (EPR) 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, @@ -8849,7 +8849,7 @@ C----------------------------------------------------------------------- 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, @@ -9012,7 +9012,7 @@ C IOPT=2 TO RESET COUNTERS, 1 FOR BOSON PAIR, 0 FOR ANY OTHERS 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, @@ -9189,7 +9189,7 @@ C----------------------------------------------------------------------- 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 @@ -9267,7 +9267,7 @@ C----------------------------------------------------------------------- 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 @@ -9435,7 +9435,7 @@ C----------------------------------------------------------------------- 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 @@ -9538,7 +9538,7 @@ C----------------------------------------------------------------------- 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 @@ -9589,7 +9589,7 @@ C----------------------------------------------------------------------- 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 @@ -9651,7 +9651,7 @@ C----------------------------------------------------------------------- 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 @@ -9928,7 +9928,7 @@ C----------------------------------------------------------------------- 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 @@ -9952,7 +9952,7 @@ C----------------------------------------------------------------------- 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 @@ -10017,7 +10017,7 @@ C A) FOR GAMinp=0 FIND AND DECAY HIGGS 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, @@ -10260,7 +10260,7 @@ C MODIFIED TO INCLUDE R-PARITY VIOLATING SUSY PR 9/4/99 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 @@ -10315,7 +10315,7 @@ C Subroutine to perform the first part of the heavy object decays 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 @@ -10447,7 +10447,7 @@ C Subroutine to perform the second part of the heavy object decays 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 @@ -10598,7 +10598,7 @@ C Subroutine to perform the third part of the heavy object decays 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 @@ -10793,7 +10793,7 @@ C Subroutine to perform the fourth part of the heavy object decays 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) @@ -10903,7 +10903,7 @@ C Subroutine to perform the fifth part of the heavy object decays 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 @@ -10968,7 +10968,7 @@ C----------------------------------------------------------------------- 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 @@ -11008,7 +11008,7 @@ C Performs partonic decays of hadrons containing heavy quark(s): 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 @@ -11250,7 +11250,7 @@ C----------------------------------------------------------------------- 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 @@ -11409,7 +11409,7 @@ C----------------------------------------------------------------------- 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, @@ -12141,7 +12141,7 @@ C----------------------------------------------------------------------- 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 @@ -12234,7 +12234,7 @@ C----------------------------------------------------------------------- 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 @@ -12684,7 +12684,7 @@ C IE generate the kinematics for the decay 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 @@ -12804,7 +12804,7 @@ C averages spin if not using TAUOLA 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 @@ -12870,7 +12870,7 @@ C----------------------------------------------------------------------- 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, @@ -13195,7 +13195,7 @@ C----------------------------------------------------------------------- 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, @@ -13475,7 +13475,7 @@ C Subroutine to perform the four body decays 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 @@ -13586,7 +13586,7 @@ C IOPT =-1 initialises 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) @@ -13729,7 +13729,7 @@ C----------------------------------------------------------------------- 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 @@ -13792,7 +13792,7 @@ C----------------------------------------------------------------------- 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 @@ -13820,7 +13820,7 @@ C----------------------------------------------------------------------- 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. @@ -13863,7 +13863,7 @@ C TERMINAL CALCULATIONS ON ELEMENTARY PROCESS 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 @@ -13986,7 +13986,7 @@ C----------------------------------------------------------------------- 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 @@ -14109,7 +14109,7 @@ C----------------------------------------------------------------------- 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 @@ -14263,7 +14263,7 @@ C----------------------------------------------------------------------- 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) @@ -14328,7 +14328,7 @@ C----------------------------------------------------------------------- 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, @@ -14562,7 +14562,7 @@ C----------------------------------------------------------------------- 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 @@ -14617,7 +14617,7 @@ C WHEN NEVHEP=0, CHOOSES X VALUES AND FINDS WEIGHT FOR PROCESS IPROC 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, @@ -15122,7 +15122,7 @@ C----------------------------------------------------------------------- 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 @@ -15217,7 +15217,7 @@ C...It factorises (PIFAC*ALPHA/SWEIN/RMW/RMW/SQRT(2.))**2 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) @@ -15310,7 +15310,7 @@ C SECTION ROUTINE 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)) @@ -15353,7 +15353,7 @@ C 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 @@ -15453,7 +15453,7 @@ C----------------------------------------------------------------------- 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 @@ -15497,7 +15497,7 @@ C----------------------------------------------------------------------- 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 @@ -15541,7 +15541,7 @@ C----------------------------------------------------------------------- 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 @@ -15786,7 +15786,7 @@ C----------------------------------------------------------------------- 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 @@ -15933,7 +15933,7 @@ C----------------------------------------------------------------------- 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, @@ -16329,7 +16329,7 @@ C----------------------------------------------------------------------- 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, @@ -16801,7 +16801,7 @@ C Subroutine to select virtual quark mass for HWH2PS 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 @@ -16841,7 +16841,7 @@ C Subroutine to select virtual quark mass for HWH2PS 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 @@ -17143,7 +17143,7 @@ C----------------------------------------------------------------------- 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 @@ -17212,7 +17212,7 @@ C I-L are the particles 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) @@ -17236,7 +17236,7 @@ C Z1 and Z2 are the decay products of the Z 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) @@ -17269,7 +17269,7 @@ C Z1 and Z2 are the decay products of the Z 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) @@ -17300,7 +17300,7 @@ C Z1 and Z2 are the decay products of the Z 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) @@ -17335,7 +17335,7 @@ C Z1 and Z2 are the decay products of the Z 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) @@ -17371,7 +17371,7 @@ C P1 is the polarization of the line I,J 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) @@ -17419,7 +17419,7 @@ C P1 is the polarization of the line I,J 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) @@ -17468,7 +17468,7 @@ C P1 is the polarization of the line I,J 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) @@ -17517,7 +17517,7 @@ C P1 is the polarization of the line I,J 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) @@ -17556,7 +17556,7 @@ C P1 is the polarization of the line I,J 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) @@ -17662,7 +17662,7 @@ C Matrix elements based on Ellis Ross & Terrano and Catani & Seymour 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, @@ -18276,7 +18276,7 @@ C 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 @@ -18602,7 +18602,7 @@ C gives the fourmomenta in the laboratory system for the particles 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, @@ -18804,7 +18804,7 @@ C corresponding Jacobian factor AJACOB 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, @@ -19148,7 +19148,7 @@ C---------------------------------------------------------------------- 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, @@ -19331,7 +19331,7 @@ C---------------------------------------------------------------------- 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, @@ -19619,7 +19619,7 @@ C 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 @@ -19809,7 +19809,7 @@ C----------------------------------------------------------------------- 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 @@ -19884,7 +19884,7 @@ C AFTER CUTS ON PT AND MASS OF CENTRE-OF-MASS SYSTEM 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), @@ -20155,7 +20155,7 @@ C---------------------------------------------------------------------- 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 @@ -20359,7 +20359,7 @@ C----------------------------------------------------------------------- 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, @@ -20480,7 +20480,7 @@ C If IPROC=107: --> gg, distributed as sum of light quarks. 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 @@ -20625,7 +20625,7 @@ c scheme, y_cut=1.-THMAX. 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, @@ -20799,7 +20799,7 @@ C----------------------------------------------------------------------- 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), @@ -21026,7 +21026,7 @@ C----------------------------------------------------------------------- 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, @@ -21275,7 +21275,7 @@ C----------------------------------------------------------------------- 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 @@ -21317,7 +21317,7 @@ C----------------------------------------------------------------------- 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 @@ -21478,7 +21478,7 @@ CDECK ID>, HWHEW0. 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 @@ -21645,7 +21645,7 @@ C NOTE: EXTERNAL FACTOR OF COLOR AVERAGE AND SPIN AVERAGE AND 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), @@ -21762,7 +21762,7 @@ C----------------------------------------------------------------------- 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, @@ -21955,7 +21955,7 @@ C----------------------------------------------------------------------- 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, @@ -22208,7 +22208,7 @@ C ICH returns the channel used if OPTM=.FALSE. 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), @@ -22428,7 +22428,7 @@ C ISM=2 select according to power law for IDBOZ 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 @@ -22491,7 +22491,7 @@ C----------------------------------------------------------------------- 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 @@ -22571,7 +22571,7 @@ C----------------------------------------------------------------------- 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 @@ -22684,7 +22684,7 @@ C----------------------------------------------------------------------- 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 @@ -22802,7 +22802,7 @@ C----------------------------------------------------------------------- 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 @@ -22845,7 +22845,7 @@ C viz L = - G(mu,nu) T(mu,nu) / GRVLAM 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 @@ -23208,7 +23208,7 @@ C---------------------------------------------------------------------- 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) @@ -23550,7 +23550,7 @@ C----------------------------------------------------------------------- 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, @@ -23830,7 +23830,7 @@ C----------------------------------------------------------------------- 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, @@ -24095,7 +24095,7 @@ C 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 @@ -24368,7 +24368,7 @@ C----------------------------------------------------------------------- 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 @@ -24715,7 +24715,7 @@ C IAPHIG (set in HWIGIN)=0: zero mass approximation =1: exact result 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, @@ -24801,7 +24801,7 @@ C----------------------------------------------------------------------- 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 @@ -24837,7 +24837,7 @@ C----------------------------------------------------------------------- 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 @@ -24873,7 +24873,7 @@ C----------------------------------------------------------------------- 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 @@ -24934,7 +24934,7 @@ C 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 @@ -25454,7 +25454,7 @@ C----------------------------------------------------------------------- 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), @@ -25746,7 +25746,7 @@ C QCD Higgs plus jet production; mean EVWGT = Sigma in nb*Higgs B.R. 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, @@ -25905,7 +25905,7 @@ C / (EM**2-EMH**2)**2 + (GAM(EM)*EMH**2/EM)**2 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 @@ -25996,7 +25996,7 @@ C 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 @@ -26521,7 +26521,7 @@ C----------------------------------------------------------------------- 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 @@ -26619,7 +26619,7 @@ C CALCULATE MOD SQUARED I DEFINED AS IN BARGER & PHILLIPS p433 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 @@ -26694,7 +26694,7 @@ C----------------------------------------------------------------------- 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 @@ -26980,7 +26980,7 @@ C----------------------------------------------------------------------- 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, @@ -27378,7 +27378,7 @@ C USES ALGORITHM OF BERENDS AND KLEISS: NUCL.PHYS. B260(1985)32 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, @@ -27589,7 +27589,7 @@ C----------------------------------------------------------------------- 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), @@ -27747,7 +27747,7 @@ C----------------------------------------------------------------------- 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 @@ -28141,7 +28141,7 @@ C----------------------------------------------------------------------- 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 @@ -28235,7 +28235,7 @@ C----------------------------------------------------------------------- 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 @@ -28393,7 +28393,7 @@ C----------------------------------------------------------------------- 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 @@ -28475,7 +28475,7 @@ C----------------------------------------------------------------------- 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 @@ -28545,7 +28545,7 @@ C Point-like photon/QCD direct meson production 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), @@ -28795,7 +28795,7 @@ C----------------------------------------------------------------------- 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 @@ -28881,7 +28881,7 @@ C----------------------------------------------------------------------- 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 @@ -28951,7 +28951,7 @@ C----------------------------------------------------------------------- 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, @@ -29305,7 +29305,7 @@ C----------------------------------------------------------------------- C----------------------------------------------------------------------- C IDENTIFIES HARD SUBPROCESS C----------------------------------------------------------------------- - INCLUDE 'HERWIG65.INC' + INCLUDE 'herwig65.inc' INTEGER HWRINT,ID3,ID4,IPERM,IHPR,ND3 EXTERNAL HWRINT IHPRO=IHPR @@ -29337,7 +29337,7 @@ C----------------------------------------------------------------------- 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 @@ -29424,7 +29424,7 @@ C----------------------------------------------------------------------- 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, @@ -29657,7 +29657,7 @@ C----------------------------------------------------------------------- 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), @@ -30592,7 +30592,7 @@ C----------------------------------------------------------------------- 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, @@ -30861,7 +30861,7 @@ C----------------------------------------------------------------------- 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 @@ -30964,7 +30964,7 @@ C----------------------------------------------------------------------- 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 @@ -30983,7 +30983,7 @@ C----------------------------------------------------------------------- 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, @@ -31631,7 +31631,7 @@ C Subroutine for resonant sleptons to standard model particles 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), @@ -31872,7 +31872,7 @@ C----------------------------------------------------------------------- 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, @@ -32677,7 +32677,7 @@ C----------------------------------------------------------------------- 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 @@ -32703,7 +32703,7 @@ C----------------------------------------------------------------------- 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) @@ -32803,7 +32803,7 @@ C JMUEO IS THE UNDERLYING EVENT OPTION: 1=>VETO EVENTS WITH M 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, @@ -33023,7 +33023,7 @@ C SWAP THE HARD PROCESS GENERATION PARAMETERS, 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 @@ -33066,7 +33066,7 @@ C MEAN EVWGT = SIGMA IN NB 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 @@ -33134,7 +33134,7 @@ C FOR PHOTON EXCHANGE, INTERFERENCE WITH U-CHANNEL CONTRIBUTION IS 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 @@ -33199,7 +33199,7 @@ C----------------------------------------------------------------------- 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), @@ -35882,7 +35882,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -35922,7 +35922,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -35955,7 +35955,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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, @@ -35988,7 +35988,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36029,7 +36029,7 @@ C section gluon gluon --> fermion antifermion (1st colour flow) 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 @@ -36067,7 +36067,7 @@ C section gluon gluon --> fermion antifermion (2st colour flow) 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 @@ -36105,7 +36105,7 @@ C section gluon gluon --> fermion antifermion (triple gluon piece) 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 @@ -36146,7 +36146,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36183,7 +36183,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36220,7 +36220,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36263,7 +36263,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36304,7 +36304,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36342,7 +36342,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36382,7 +36382,7 @@ C N.B. a gauge choice has been made to simplify the triple gluon vertex 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, @@ -36422,7 +36422,7 @@ C N.B. a gauge choice has been made to simplify the triple gluon vertex 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, @@ -36462,7 +36462,7 @@ C N.B. a gauge choice has been made to simplify the triple gluon vertex 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 @@ -36504,7 +36504,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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, @@ -36547,7 +36547,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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, @@ -36590,7 +36590,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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, @@ -36633,7 +36633,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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, @@ -36676,7 +36676,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36711,7 +36711,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36744,7 +36744,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36777,7 +36777,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36812,7 +36812,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36845,7 +36845,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36878,7 +36878,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36913,7 +36913,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36947,7 +36947,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -36981,7 +36981,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -37016,7 +37016,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -37050,7 +37050,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -37084,7 +37084,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -37119,7 +37119,7 @@ C Subroutine to calculate the helicity amplitudes for the 2-to-2 cross 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 @@ -37168,7 +37168,7 @@ C----------------------------------------------------------------------- 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 = @@ -37192,7 +37192,7 @@ C----------------------------------------------------------------------- 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, @@ -37884,7 +37884,7 @@ C----------------------------------------------------------------------- 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 @@ -38102,7 +38102,7 @@ C----------------------------------------------------------------------- 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, @@ -39013,7 +39013,7 @@ C----------------------------------------------------------------------- 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 @@ -39090,7 +39090,7 @@ C----------------------------------------------------------------------- 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 @@ -39111,7 +39111,7 @@ C V + 1 JET PRODUCTION, WHERE V=W (IHPRO.LT.5) OR Z (IHPRO.GE.5). 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 @@ -39308,7 +39308,7 @@ C Vector Boson production with two hard jets 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, @@ -39596,7 +39596,7 @@ C UNLESS USER SPECIFIES OTHERWISE BY MOD(IPROC,100)=1-8 RESPECTIVELY 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 @@ -39677,7 +39677,7 @@ C----------------------------------------------------------------------- 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) @@ -39818,7 +39818,7 @@ CDECK ID>, HWICHK. 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' @@ -39837,7 +39837,7 @@ C =1 PDG; =2 HERWIG numeric; =3 HERWIG character name. 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) @@ -39918,7 +39918,7 @@ C---------------------------------------------------------------------- C----------------------------------------------------------------------- C SETS INPUT PARAMETERS C---------------------------------------------------------------------- - INCLUDE 'HERWIG65.INC' + INCLUDE 'herwig65.inc' DOUBLE PRECISION FAC,ANGLE INTEGER I,J,N,L CHARACTER*28 TITLE @@ -40533,7 +40533,7 @@ C---------------------------------------------------------------------- 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 @@ -40633,7 +40633,7 @@ C----------------------------------------------------------------------- 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 @@ -40743,7 +40743,7 @@ C IOPT = 1 sets the weights for the different channels to their 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), @@ -41012,7 +41012,7 @@ C----------------------------------------------------------------------- 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 @@ -41232,7 +41232,7 @@ C----------------------------------------------------------------------- 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 @@ -41355,7 +41355,7 @@ C----------------------------------------------------------------------- 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 @@ -42569,7 +42569,7 @@ C----------------------------------------------------------------------- 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), @@ -44639,7 +44639,7 @@ C----------------------------------------------------------------------- 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), @@ -44741,7 +44741,7 @@ C----------------------------------------------------------------------- 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 @@ -44989,7 +44989,7 @@ C IPROC = 1000,... ADDS SOFT UNDERLYING EVENT 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), @@ -45274,7 +45274,7 @@ C----------------------------------------------------------------------- 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) @@ -45418,7 +45418,7 @@ C adjusts pointers, sets CMMOM (for two-body mode) and resets RSTAB 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 @@ -45574,7 +45574,7 @@ C----------------------------------------------------------------------- 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 @@ -45631,7 +45631,7 @@ C----------------------------------------------------------------------- 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 @@ -45693,7 +45693,7 @@ C----------------------------------------------------------------------- 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 @@ -45744,7 +45744,7 @@ C----------------------------------------------------------------------- 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 @@ -45893,7 +45893,7 @@ C----------------------------------------------------------------------- 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 @@ -45908,28 +45908,29 @@ CDECK ID>, 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 @@ -46014,7 +46015,7 @@ C----------------------------------------------------------------------- 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 @@ -46101,7 +46102,7 @@ C----------------------------------------------------------------------- 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, @@ -46665,7 +46666,7 @@ C 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 @@ -46874,7 +46875,7 @@ C CERN_ROOT:[DOC]PDFLIB.TXT ON VXCERN 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), @@ -47262,7 +47263,7 @@ C---FIX TO CALL SCHULER-SJOSTRAND CODE 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 @@ -47362,7 +47363,7 @@ C---USE PDFLIB NUCLEON STRUCTURE FUNCTIONS 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 @@ -47667,7 +47668,7 @@ C----------------------------------------------------------------------- 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 @@ -47695,7 +47696,7 @@ C----------------------------------------------------------------------- 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) @@ -47718,7 +47719,7 @@ C----------------------------------------------------------------------- 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 @@ -47840,7 +47841,7 @@ C----------------------------------------------------------------------- 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 @@ -47947,7 +47948,7 @@ CDECK ID>, HWSSUD. 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 @@ -48038,7 +48039,7 @@ C----------------------------------------------------------------------- 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. @@ -48069,7 +48070,7 @@ C Running electromagnetic coupling constant. 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 @@ -48152,7 +48153,7 @@ C .EQ.2 RATIO OF ABOVE TO ONE-LOOP 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 @@ -48244,7 +48245,7 @@ C----------------------------------------------------------------------- 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 @@ -48282,7 +48283,7 @@ C----------------------------------------------------------------------- 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) @@ -48301,7 +48302,7 @@ C----------------------------------------------------------------------- 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 @@ -48333,7 +48334,7 @@ C----------------------------------------------------------------------- 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 @@ -48386,7 +48387,7 @@ C labels final fermion; type given as: 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 @@ -48501,7 +48502,7 @@ C QORQQB(I) = .TRUE. if it is a quark or an antidiquark 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 @@ -59117,7 +59118,7 @@ C 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 @@ -59144,7 +59145,7 @@ C Gives warnings if: a particle has no decay modes or antiparticle's 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), @@ -59462,7 +59463,7 @@ C Designed to be printed as landscape 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 @@ -59800,7 +59801,7 @@ C EDIT THE EVENT RECORD 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 @@ -59865,7 +59866,7 @@ C----------------------------------------------------------------------- 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 @@ -59891,7 +59892,7 @@ C----------------------------------------------------------------------- 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 @@ -59940,7 +59941,7 @@ C Call HWUDPR to create particle property files in 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 @@ -60373,7 +60374,7 @@ C----------------------------------------------------------------------- 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 @@ -60389,7 +60390,7 @@ C----------------------------------------------------------------------- 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 @@ -60420,6 +60421,7 @@ C---CHECK FOR FATAL ERROR 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 @@ -60491,7 +60493,7 @@ C IOPT= 1 GIVEN IPDG, RETURNS IWIG AND NWIG 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 @@ -60549,7 +60551,7 @@ C----------------------------------------------------------------------- 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) @@ -61367,7 +61369,7 @@ C---FIX TO CALL SCHULER-SJOSTRAND CODE 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 @@ -61455,7 +61457,7 @@ C----------------------------------------------------------------------- C----------------------------------------------------------------------- C INITIALISES AN EVENT C----------------------------------------------------------------------- - INCLUDE 'HERWIG65.INC' + INCLUDE 'herwig65.inc' DOUBLE PRECISION HWRGEN,HWRGET,DUMMY REAL TL LOGICAL CALLED,HWRLOG @@ -61677,7 +61679,7 @@ C----------------------------------------------------------------------- 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 @@ -61868,7 +61870,7 @@ C---------------------------------------------------------------------- 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 @@ -61900,7 +61902,7 @@ C ( rep. enhancement factor)^2 - used in cluster decays 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, @@ -62264,7 +62266,7 @@ C----------------------------------------------------------------------- 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 @@ -62363,7 +62365,7 @@ C----------------------------------------------------------------------- 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) @@ -62394,7 +62396,7 @@ C----------------------------------------------------------------------- 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) @@ -62597,7 +62599,7 @@ C 300-399 FATAL, DUMP EVENT & TERMINATE RUN 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