C Last modification on March 10th 2001 by M.S. C ================================================================== C ================= PROGRAM HDECAY: COMMENTS ======================= C ================================================================== C C *************** C * VERSION 2.0 * C *************** C C C This program calculates the total decay widths and the branching C ratios of the C Standard Model Higgs boson (HSM) as well as those C of the neutral (HL= the light CP-even, HH= the heavy CP-even, HA= C the pseudoscalar) and the charged (HC) Higgs bosons of the Minimal C Supersymmetric extension of the Standard Model (MSSM). It includes: C C - All the decay channels which are kinematically allowed and which C have branching ratios larger than 10**(-4). C C - All QCD corrections to the fermionic and gluonic decay modes. C Most of these corrections are mapped into running masses in a C consistent way with some freedom for including high order terms. C C - Below--threshold three--body decays with off--shell top quarks C or ONE off-shell gauge boson, as well as some decays with one C off-shell Higgs boson in the MSSM. C C - Double off-shell decays: HSM,HL,HH --> W*W*,Z*Z* -->4 fermions, C which could be important for Higgs masses close to MW or MZ. C C - In the MSSM, the radiative corrections with full squark mixing and C uses the RG improved values of Higgs masses and couplings with the C main NLO corrections implemented (based on M.Carena, M. Quiros and C C.E.M. Wagner, Nucl. Phys. B461 (1996) 407, hep-ph/9508343). C C - In the MSSM, all the decays into CHARGINOS, NEUTRALINOS, SLEPTONS C and SQUARKS (with mixing in the stop and sbottom sectors). C C - Chargino, slepton and squark loops in the 2 photon decays and squark C loops in the gluonic decays (including QCD corrections). C C =================================================================== C This program has been written by A.Djouadi, J.Kalinowski and M.Spira. C For details on how to use the program see: Comp. Phys. Commun. 108 C (1998) 56, hep-ph/9704448. For any question, comment, suggestion or C complaint, please contact us at: C djouadi@lpm.univ-montp2.fr C kalino@fuw.edu.pl C Michael.Spira@cern.ch C ================ IT USES AS INPUT PARAMETERS: C C IHIGGS: =0: CALCULATE BRANCHING RATIOS OF SM HIGGS BOSON C =1: CALCULATE BRANCHING RATIOS OF MSSM h BOSON C =2: CALCULATE BRANCHING RATIOS OF MSSM H BOSON C =3: CALCULATE BRANCHING RATIOS OF MSSM A BOSON C =4: CALCULATE BRANCHING RATIOS OF MSSM H+ BOSON C =5: CALCULATE BRANCHING RATIOS OF ALL MSSM HIGGS BOSONS C C TGBET: TAN(BETA) FOR MSSM C MABEG: START VALUE OF M_A FOR MSSM AND M_H FOR SM C MAEND: END VALUE OF M_A FOR MSSM AND M_H FOR SM C NMA: NUMBER OF ITERATIONS FOR M_A C ALS(MZ): VALUE FOR ALPHA_S(M_Z) C MSBAR(1): MSBAR MASS OF STRANGE QUARK AT SCALE Q=1 GEV C MC: CHARM POLE MASS C MB: BOTTOM POLE MASS C MT: TOP POLE MASS C MTAU: TAU MASS C MMUON: MUON MASS C ALPH: INVERSE QED COUPLING C GF: FERMI CONSTANT C GAMW: W WIDTH C GAMZ: Z WIDTH C MZ: Z MASS C MW: W MASS C VUS: CKM PARAMETER V_US C VCB: CKM PARAMETER V_CB C VUB/VCB: RATIO V_UB/V_CB C 1ST AND 2ND GENERATION: C MSL1: SUSY BREAKING MASS PARAMETERS OF LEFT HANDED SLEPTONS C MER1: SUSY BREAKING MASS PARAMETERS OF RIGHT HANDED SLEPTONS C MQL1: SUSY BREAKING MASS PARAMETERS OF LEFT HANDED SUPS C MUR1: SUSY BREAKING MASS PARAMETERS OF RIGHT HANDED SUPS C MDR1: SUSY BREAKING MASS PARAMETERS OF RIGHT HANDED SDOWNS C 3RD GENERATION: C MSL: SUSY BREAKING MASS PARAMETERS OF LEFT HANDED STAUS C MER: SUSY BREAKING MASS PARAMETERS OF RIGHT HANDED STAUS C MSQ: SUSY BREAKING MASS PARAMETERS OF LEFT HANDED STOPS C MUR: SUSY BREAKING MASS PARAMETERS OF RIGHT HANDED STOPS C MDR: SUSY BREAKING MASS PARAMETERS OF RIGHT HANDED SBOTTOMS C AL: STAU TRILINEAR SOFT BREAKING TERMS C AU: STOP TRILINEAR SOFT BREAKING TERMS. C AD: SBOTTOM TRILINEAR SOFT BREAKING TERMS. C MU: SUSY HIGGS MASS PARAMETER C M2: GAUGINO MASS PARAMETER. C C NNLO (M): =0: USE O(ALPHA_S) FORMULA FOR POLE MASS --> MSBAR MASS C =1: USE O(ALPHA_S**2) FORMULA FOR POLE MASS --> MSBAR MASS C C ON-SHELL: =0: INCLUDE OFF_SHELL DECAYS H,A --> T*T*, A --> Z*H, C H --> W*H+,Z*A, H+ --> W*A, W*H, T*B C =1: EXCLUDE THE OFF-SHELL DECAYS ABOVE C C ON-SH-WZ: =0: INCLUDE DOUBLE OFF-SHELL PAIR DECAYS PHI --> W*W*,Z*Z* C =1: INCLUDE ONLY SINGLE OFF-SHELL DECAYS PHI --> W*W,Z*Z C C IPOLE: =0 COMPUTES RUNNING HIGGS MASSES (FASTER) C =1 COMPUTES POLE HIGGS MASSES C C OFF-SUSY: =0: INCLUDE DECAYS (AND LOOPS) INTO SUPERSYMMETRIC PARTICLES C =1: EXCLUDE DECAYS (AND LOOPS) INTO SUPERSYMMETRIC PARTICLES C C INIDEC: =0: PRINT OUT SUMS OF CHARGINO/NEUTRALINO/SFERMION DECAYS C =1: PRINT OUT INDIVIDUAL CHARGINO/NEUTRALINO/SFERMION DECAYS C C NF-GG: NUMBER OF LIGHT FLAVORS INCLUDED IN THE GLUONIC DECAYS C PHI --> GG* --> GQQ (3,4 OR 5) C C ======================================================================= C ============== BEGINNING OF THE MAIN PROGRAM ========================== C ======================================================================= C PROGRAM HDECAY IMPLICIT DOUBLE PRECISION (A-H,O-Z) PARAMETER(K=6,NI=87,NSA=85,NSB=86,NLA=88,NLB=89,NHA=90,NHB=91, . NHC=92,NAA=93,NAB=94,NCA=95,NCB=96,NRA=97,NRB=98, . NSUSYL=81,NSUSYA=82,NSUSYH=83,NSUSYC=84,NPAR=80, . NSUSYLA=79,NSUSYLB=78,NSUSYLC=77,NSUSYLD=76,NSUSYLE=75, . NSUSYLF=59,NSUSYHF=58, . NSUSYHA=74,NSUSYHB=73,NSUSYHC=72,NSUSYHD=71,NSUSYHE=70, . NSUSYAA=69,NSUSYAB=68,NSUSYAC=67,NSUSYAD=66,NSUSYAE=65, . NSUSYCA=64,NSUSYCB=63,NSUSYCC=62,NSUSYCD=61,NSUSYCE=60) DIMENSION GMN(4),XMN(4),GMC(2),GMST(2),GMSB(2),GMSL(2), . GMSU(2),GMSD(2),GMSE(2),GMSN(2) DIMENSION HLBRSC(2,2),HLBRSN(4,4),HHBRSC(2,2),HHBRSN(4,4), . HABRSC(2,2),HABRSN(4,4),HCBRSU(2,4), . HHBRST(2,2),HHBRSB(2,2),HCBRSTB(2,2) DIMENSION AC1(2,2),AC2(2,2),AC3(2,2), . AN1(4,4),AN2(4,4),AN3(4,4), . ACNL(2,4),ACNR(2,4) DIMENSION GLTT(2,2),GLBB(2,2),GHTT(2,2),GHBB(2,2),GCTB(2,2), . GLEE(2,2),GHEE(2,2),GCEN(2,2) COMMON/MASSES/AMS,AMC,AMB,AMT COMMON/STRANGE/AMSB COMMON/PARAM/GF,ALPH,AMTAU,AMMUON,AMZ,AMW COMMON/CKMPAR/VUS,VCB,VUB COMMON/HMASS/AMSM,AMA,AML,AMH,AMCH,AMAR COMMON/BREAK/AMEL,AMER,AMSQ,AMUR,AMDR,AL,AU,AD,AMU,AM2 COMMON/SFER1ST/AMQL1,AMUR1,AMDR1,AMEL1,AMER1 COMMON/WZWDTH/GAMC0,GAMT0,GAMT1,GAMW,GAMZ COMMON/COUP/GAT,GAB,GLT,GLB,GHT,GHB,GZAH,GZAL, . GHHH,GLLL,GHLL,GLHH,GHAA,GLAA,GLVV,GHVV, . GLPM,GHPM,B,A COMMON/ALS/XLAMBDA,AMC0,AMB0,AMT0,N0 COMMON/FLAG/IHIGGS,NNLO,IPOLE COMMON/ONSHELL/IONSH,IONWZ,IOFSUSY COMMON/OLDFASH/NFGG COMMON/WIDTHSM/SMBRB,SMBRL,SMBRM,SMBRS,SMBRC,SMBRT,SMBRG,SMBRGA, . SMBRZGA,SMBRW,SMBRZ,SMWDTH COMMON/WIDTHA/ABRB,ABRL,ABRM,ABRS,ABRC,ABRT,ABRG,ABRGA,ABRZGA, . ABRZ,AWDTH COMMON/WIDTHHL/HLBRB,HLBRL,HLBRM,HLBRS,HLBRC,HLBRT,HLBRG,HLBRGA, . HLBRZGA,HLBRW,HLBRZ,HLBRA,HLBRAZ,HLBRHW,HLWDTH COMMON/WIDTHHH/HHBRB,HHBRL,HHBRM,HHBRS,HHBRC,HHBRT,HHBRG,HHBRGA, . HHBRZGA,HHBRW,HHBRZ,HHBRH,HHBRA,HHBRAZ,HHBRHW, . HHWDTH COMMON/WIDTHHC/HCBRB,HCBRL,HCBRM,HCBRBU,HCBRS,HCBRC,HCBRT,HCBRW, . HCBRA,HCWDTH COMMON/WISUSY/HLBRSC,HLBRSN,HHBRSC,HHBRSN,HABRSC,HABRSN,HCBRSU, . HLBRCHT,HHBRCHT,HABRCHT,HLBRNET,HHBRNET,HABRNET, . HCBRCNT,HLBRSL,HHBRSL,HCBRSL,HABRSL,HABRST,HABRSB, . HHBRSQ,HHBRST,HHBRSB,HHBRSQT,HCBRSQ,HCBRSTB, . HCBRSQT,HLBRSQ,HLBRSQT COMMON/WISFER/BHLSLNL,BHLSLEL,BHLSLER,BHLSQUL,BHLSQUR,BHLSQDL, . BHLSQDR,BHLST(2,2),BHLSB(2,2),BHLSTAU(2,2), . BHHSLNL,BHHSLEL,BHHSLER,BHHSQUL,BHHSQUR,BHHSQDL, . BHHSQDR,BHHST(2,2),BHHSB(2,2),BHHSTAU(2,2), . BHASTAU,BHASB,BHAST, . BHCSL00,BHCSL11,BHCSL21,BHCSQ,BHCSTB(2,2) COMMON/SMASS/GMN,XMN,GMC,GMST,GMSB,GMSL,GMSU,GMSD,GMSE,GMSN PI = 4*DATAN(1D0) OPEN(NI,FILE='hdecay.in') OPEN(NPAR,FILE='br.input') READ(NI,101)IHIGGS READ(NI,100)TGBET READ(NI,100)AMABEG READ(NI,100)AMAEND READ(NI,101)NMA READ(NI,100)ALSMZ READ(NI,100)AMS READ(NI,100)AMC READ(NI,100)AMB READ(NI,100)AMT READ(NI,100)AMTAU READ(NI,100)AMMUON READ(NI,100)ALPH READ(NI,100)GF READ(NI,100)GAMW READ(NI,100)GAMZ READ(NI,100)AMZ READ(NI,100)AMW READ(NI,100)VUS READ(NI,100)VCB READ(NI,100)RVUB READ(NI,100)AMU READ(NI,100)AM2 READ(NI,100)AMEL1 READ(NI,100)AMER1 READ(NI,100)AMQL1 READ(NI,100)AMUR1 READ(NI,100)AMDR1 READ(NI,100)AMEL READ(NI,100)AMER READ(NI,100)AMSQ READ(NI,100)AMUR READ(NI,100)AMDR READ(NI,100)AL READ(NI,100)AU READ(NI,100)AD READ(NI,101)NNLO READ(NI,101)IONSH READ(NI,101)IONWZ READ(NI,101)IPOLE READ(NI,101)IOFSUSY READ(NI,101)INDIDEC READ(NI,101)NFGG VUB=RVUB*VCB ALPH=1.D0/ALPH AMSB = AMS C--CHECK NFGG IF(NFGG.GT.5.OR.NFGG.LT.3)THEN WRITE(6,*)'NF-GG NOT VALID. TAKING THE DEFAULT NF-GG = 3....' NFGG = 3 ENDIF 100 FORMAT(10X,G30.20) 101 FORMAT(10X,I30) IF(IHIGGS.EQ.0) THEN OPEN(NSA,FILE='br.sm1') OPEN(NSB,FILE='br.sm2') ENDIF IF(IHIGGS.EQ.1.OR.IHIGGS.EQ.5) THEN OPEN(NLA,FILE='br.l1') OPEN(NLB,FILE='br.l2') IF(IOFSUSY.EQ.0)THEN OPEN(NSUSYL,FILE='br.ls') IF(INDIDEC.NE.0)THEN OPEN(NSUSYLA,FILE='br.ls1') OPEN(NSUSYLB,FILE='br.ls2') OPEN(NSUSYLC,FILE='br.ls3') OPEN(NSUSYLD,FILE='br.ls4') OPEN(NSUSYLE,FILE='br.ls5') OPEN(NSUSYLF,FILE='br.ls6') ENDIF ENDIF ENDIF IF(IHIGGS.EQ.2.OR.IHIGGS.EQ.5) THEN OPEN(NHA,FILE='br.h1') OPEN(NHB,FILE='br.h2') OPEN(NHC,FILE='br.h3') IF(IOFSUSY.EQ.0)THEN OPEN(NSUSYH,FILE='br.hs') IF(INDIDEC.NE.0)THEN OPEN(NSUSYHA,FILE='br.hs1') OPEN(NSUSYHB,FILE='br.hs2') OPEN(NSUSYHC,FILE='br.hs3') OPEN(NSUSYHD,FILE='br.hs4') OPEN(NSUSYHE,FILE='br.hs5') OPEN(NSUSYHF,FILE='br.hs6') ENDIF ENDIF ENDIF IF(IHIGGS.EQ.3.OR.IHIGGS.EQ.5) THEN OPEN(NAA,FILE='br.a1') OPEN(NAB,FILE='br.a2') IF(IOFSUSY.EQ.0)THEN OPEN(NSUSYA,FILE='br.as') IF(INDIDEC.NE.0)THEN OPEN(NSUSYAA,FILE='br.as1') OPEN(NSUSYAB,FILE='br.as2') OPEN(NSUSYAC,FILE='br.as3') OPEN(NSUSYAD,FILE='br.as4') ENDIF ENDIF ENDIF IF(IHIGGS.EQ.4.OR.IHIGGS.EQ.5) THEN OPEN(NCA,FILE='br.c1') OPEN(NCB,FILE='br.c2') IF(IOFSUSY.EQ.0)THEN OPEN(NSUSYC,FILE='br.cs') IF(INDIDEC.NE.0)THEN OPEN(NSUSYCA,FILE='br.cs1') OPEN(NSUSYCB,FILE='br.cs2') OPEN(NSUSYCC,FILE='br.cs3') OPEN(NSUSYCD,FILE='br.cs4') ENDIF ENDIF ENDIF AMC0=AMC AMB0=AMB AMT0=AMT ACC=1.D-8 NLOOP=2 XLAMBDA=XITLA(NLOOP,ALSMZ,ACC) N0=5 CALL ALSINI(ACC) C--INITIALIZE COEFFICIENTS FOR POLYLOGARITHMS NBER = 18 CALL BERNINI(NBER) C--WRITE THE INPUT PARAMTERS TO A DATA-FILE WRITE(NPAR,8)'HIGGS = ',IHIGGS WRITE(NPAR,9)'TGBET = ',TGBET WRITE(NPAR,9)'MABEG = ',AMABEG WRITE(NPAR,9)'MAEND = ',AMAEND WRITE(NPAR,7)'NMA = ',NMA WRITE(NPAR,9)'ALS(MZ) = ',ALSMZ WRITE(NPAR,9)'MSBAR(1) = ',AMS WRITE(NPAR,9)'MC = ',AMC WRITE(NPAR,9)'MB = ',AMB WRITE(NPAR,9)'MT = ',AMT WRITE(NPAR,9)'MTAU = ',AMTAU WRITE(NPAR,9)'MMUON = ',AMMUON WRITE(NPAR,9)'ALPH = ',1.D0/ALPH WRITE(NPAR,9)'GF = ',GF WRITE(NPAR,9)'GAMW = ',GAMW WRITE(NPAR,9)'GAMZ = ',GAMZ WRITE(NPAR,9)'MZ = ',AMZ WRITE(NPAR,9)'MW = ',AMW WRITE(NPAR,9)'VUS = ',VUS WRITE(NPAR,9)'VCB = ',VCB WRITE(NPAR,9)'VUB/VCB = ',RVUB WRITE(NPAR,9)'MU = ',AMU WRITE(NPAR,9)'M2 = ',AM2 WRITE(NPAR,9)'MEL1 = ',AMEL1 WRITE(NPAR,9)'MER1 = ',AMER1 WRITE(NPAR,9)'MQL1 = ',AMQL1 WRITE(NPAR,9)'MUR1 = ',AMUR1 WRITE(NPAR,9)'MDR1 = ',AMDR1 WRITE(NPAR,9)'MEL = ',AMEL WRITE(NPAR,9)'MER = ',AMER WRITE(NPAR,9)'MSQ = ',AMSQ WRITE(NPAR,9)'MUR = ',AMUR WRITE(NPAR,9)'MDR = ',AMDR WRITE(NPAR,9)'AL = ',AL WRITE(NPAR,9)'AU = ',AU WRITE(NPAR,9)'AD = ',AD WRITE(NPAR,8)'NNLO (M) = ',NNLO WRITE(NPAR,8)'ON-SHELL = ',IONSH WRITE(NPAR,8)'ON-SH-WZ = ',IONWZ WRITE(NPAR,8)'OFF-SUSY = ',IOFSUSY WRITE(NPAR,8)'IPOLE = ',IPOLE WRITE(NPAR,8)'NF-GG = ',NFGG WRITE(NPAR,9)'LAMBDA_5 = ',XLAMBDA CLOSE(NPAR) 7 FORMAT(A11,I7) 8 FORMAT(A11,I4) 9 FORMAT(A11,G15.6) C--SETUP THE HEADS OF THE TABLES IN THE DATA-FILES IF(IHIGGS.EQ.0) THEN WRITE(NSA,70)'MHSM ','BB ','TAU TAU','MU MU ','SS ','CC ','TT ' WRITE(NSA,69) WRITE(NSA,*) WRITE(NSB,70)'MHSM ','GG ','GAM GAM','Z GAM ','WW ','ZZ ','WIDTH' WRITE(NSB,69) WRITE(NSB,*) ENDIF IF(IHIGGS.EQ.1.OR.IHIGGS.EQ.5) THEN WRITE(NLA,70)'MHL ','BB ','TAU TAU','MU MU ','SS ','CC ','TT ' WRITE(NLA,69) WRITE(NLA,*) WRITE(NLB,70)'MHL ','GG ','GAM GAM','Z GAM ','WW ','ZZ ','WIDTH' WRITE(NLB,69) WRITE(NLB,*) ENDIF IF(IHIGGS.EQ.2.OR.IHIGGS.EQ.5) THEN WRITE(NHA,70)'MHH ','BB ','TAU TAU','MU MU ','SS ','CC ','TT ' WRITE(NHA,69) WRITE(NHA,*) WRITE(NHB,72)'MHH ','GG ','GAM GAM','Z GAM ','WW ','ZZ ' WRITE(NHB,69) WRITE(NHB,*) WRITE(NHC,72)'MHH ','hh ','AA ','Z A ','W+- H-+','WIDTH ' WRITE(NHC,69) WRITE(NHC,*) ENDIF IF(IHIGGS.EQ.3.OR.IHIGGS.EQ.5) THEN WRITE(NAA,70)'MHA ','BB ','TAU TAU','MU MU ','SS ','CC ','TT ' WRITE(NAA,69) WRITE(NAA,*) WRITE(NAB,72)'MHA ','GG ','GAM GAM','Z GAM ','Z HL ','WIDTH ' WRITE(NAB,69) WRITE(NAB,*) ENDIF IF(IHIGGS.EQ.4.OR.IHIGGS.EQ.5) THEN WRITE(NCA,70)'MHC ','BC ','TAU NU ','MU NU ','SU ','CS ','TB ' WRITE(NCA,69) WRITE(NCA,*) WRITE(NCB,70)'MHC ','hW ','AW ','WIDTH ' WRITE(NCB,69) WRITE(NCB,*) ENDIF 69 FORMAT(79('_')) 70 FORMAT(A9,6(1X,A10)) 71 FORMAT(A9,4(1X,A10)) 72 FORMAT(A9,5(1X,A10)) 73 FORMAT(A9,3(1X,A10)) IWRISU = 1 DO 9999 II=1,NMA IF(NMA.NE.1)THEN AMAR = AMABEG + (AMAEND-AMABEG)/(NMA-1D0)*(II-1D0) ELSE AMAR = AMABEG ENDIF AMSM = AMAR AMA = AMAR IF(IHIGGS.NE.0)THEN C ******************************* SUSY OUTPUT CALL SUSYCP(TGBET) IF(IOFSUSY.EQ.0.AND.IWRISU.NE.0)THEN C--WRITE THE GAUGINO MASSES/ TB, MU AND M2 IN THE SUSY DATA-FILE C--WRITE THE SFERMION MASSES/ SUSY MASSES AND COUPLINGS IN SUSY DATA-FILE CALL GAUGINO(AMU,AM2,B,A,GMC,GMN,XMN,AC1,AC2,AC3, . AN1,AN2,AN3,ACNL,ACNR) CALL SFERMION(AMSQ,AMUR,AMDR,AMEL,AMER,AL,AU,AD,AMU, . GMST,GMSB,GMSL,GMSU,GMSD,GMSE,GMSN, . GLEE,GLTT,GLBB,GHEE,GHTT,GHBB, . GAEE,GATT,GABB,GCEN,GCTB) C IF(IHIGGS.EQ.1.OR.IHIGGS.EQ.5) THEN WRITE(NSUSYL,347) TGBET,AM2,AMU,AMSQ WRITE(NSUSYL,348) GMC(1),GMC(2),GMN(1),GMN(2),GMN(3),GMN(4) WRITE(NSUSYL,349) GMST(1),GMST(2),GMSU(1),GMSU(2) WRITE(NSUSYL,350) GMSB(1),GMSB(2),GMSD(1),GMSD(2) WRITE(NSUSYL,351) GMSL(1),GMSL(2),GMSE(1),GMSE(2),GMSN(1) WRITE(NSUSYL,*) WRITE(NSUSYL,*)' MHL CHARGINOS NEUTRALS '// . 'SLEPTONS SQUARKS' WRITE(NSUSYL,69) WRITE(NSUSYL,*) IF(INDIDEC.NE.0)THEN WRITE(NSUSYLA,73)'MHL ','C1 C1 ','C2 C2 ','C1 C2 ' WRITE(NSUSYLA,69) WRITE(NSUSYLA,*) WRITE(NSUSYLB,71)'MHL ','N1 N1 ','N2 N2 ','N3 N3 ','N4 N4 ' WRITE(NSUSYLB,69) WRITE(NSUSYLB,*) WRITE(NSUSYLC,70)'MHL ','N1 N2 ','N1 N3 ','N1 N4 ','N2 N3 ', . 'N2 N4 ','N3 N4 ' WRITE(NSUSYLC,69) WRITE(NSUSYLC,*) WRITE(NSUSYLD,*)' MHL SNL SNL SEL SEL '// . 'SER SER STA1 STA1 STA1 STA2 STA2 STA2' WRITE(NSUSYLD,69) WRITE(NSUSYLD,*) WRITE(NSUSYLE,*)' MHL SUL SUL SUR SUR '// . 'SDL SDL SDR SDR' WRITE(NSUSYLE,69) WRITE(NSUSYLE,*) WRITE(NSUSYLF,*)' MHL SB1 SB1 SB1 SB2 '// . 'SB2 SB2 ST1 ST1 ST1 ST2 ST2 ST2' WRITE(NSUSYLF,69) WRITE(NSUSYLF,*) ENDIF ENDIF IF(IHIGGS.EQ.2.OR.IHIGGS.EQ.5) THEN WRITE(NSUSYH,347) TGBET,AM2,AMU,AMSQ WRITE(NSUSYH,348) GMC(1),GMC(2),GMN(1),GMN(2),GMN(3),GMN(4) WRITE(NSUSYH,349) GMST(1),GMST(2),GMSU(1),GMSU(2) WRITE(NSUSYH,350) GMSB(1),GMSB(2),GMSD(1),GMSD(2) WRITE(NSUSYH,351) GMSL(1),GMSL(2),GMSE(1),GMSE(2),GMSN(1) WRITE(NSUSYH,*) WRITE(NSUSYH,*)' MHH CHARGINOS NEUTRALS '// . 'SLEPTONS SQUARKS' WRITE(NSUSYH,69) WRITE(NSUSYH,*) IF(INDIDEC.NE.0)THEN WRITE(NSUSYHA,73)'MHH ','C1 C1 ','C2 C2 ','C1 C2 ' WRITE(NSUSYHA,69) WRITE(NSUSYHA,*) WRITE(NSUSYHB,71)'MHH ','N1 N1 ','N2 N2 ','N3 N3 ','N4 N4 ' WRITE(NSUSYHB,69) WRITE(NSUSYHB,*) WRITE(NSUSYHC,70)'MHH ','N1 N2 ','N1 N3 ','N1 N4 ','N2 N3 ', . 'N2 N4 ','N3 N4 ' WRITE(NSUSYHC,69) WRITE(NSUSYHC,*) WRITE(NSUSYHD,*)' MHH SNL SNL SEL SEL '// . 'SER SER STA1 STA1 STA1 STA2 STA2 STA2' WRITE(NSUSYHD,69) WRITE(NSUSYHD,*) WRITE(NSUSYHE,*)' MHH SUL SUL SUR SUR '// . 'SDL SDL SDR SDR' WRITE(NSUSYHE,69) WRITE(NSUSYHE,*) WRITE(NSUSYHF,*)' MHH SB1 SB1 SB1 SB2 '// . 'SB2 SB2 ST1 ST1 ST1 ST2 ST2 ST2' WRITE(NSUSYHF,69) WRITE(NSUSYHF,*) ENDIF ENDIF IF(IHIGGS.EQ.3.OR.IHIGGS.EQ.5) THEN WRITE(NSUSYA,347) TGBET,AM2,AMU,AMSQ WRITE(NSUSYA,348) GMC(1),GMC(2),GMN(1),GMN(2),GMN(3),GMN(4) WRITE(NSUSYA,349) GMST(1),GMST(2),GMSU(1),GMSU(2) WRITE(NSUSYA,350) GMSB(1),GMSB(2),GMSD(1),GMSD(2) WRITE(NSUSYA,351) GMSL(1),GMSL(2),GMSE(1),GMSE(2),GMSN(1) WRITE(NSUSYA,*) WRITE(NSUSYA,*)' MHA CHARGINOS NEUTRALS '// . 'SLEPTONS SQUARKS' WRITE(NSUSYA,69) WRITE(NSUSYA,*) IF(INDIDEC.NE.0)THEN WRITE(NSUSYAA,73)'MHA ','C1 C1 ','C2 C2 ','C1 C2 ' WRITE(NSUSYAA,69) WRITE(NSUSYAA,*) WRITE(NSUSYAB,71)'MHA ','N1 N1 ','N2 N2 ','N3 N3 ','N4 N4 ' WRITE(NSUSYAB,69) WRITE(NSUSYAB,*) WRITE(NSUSYAC,70)'MHA ','N1 N2 ','N1 N3 ','N1 N4 ','N2 N3 ', . 'N2 N4 ','N3 N4 ' WRITE(NSUSYAC,69) WRITE(NSUSYAC,*) WRITE(NSUSYAD,*) WRITE(NSUSYAD,*)' MHA STA1 STA2 SB1 SB2 ST1 ST2' WRITE(NSUSYAD,69) WRITE(NSUSYAD,*) ENDIF ENDIF IF(IHIGGS.EQ.4.OR.IHIGGS.EQ.5) THEN WRITE(NSUSYC,347) TGBET,AM2,AMU,AMSQ WRITE(NSUSYC,348) GMC(1),GMC(2),GMN(1),GMN(2),GMN(3),GMN(4) WRITE(NSUSYC,349) GMST(1),GMST(2),GMSU(1),GMSU(2) WRITE(NSUSYC,350) GMSB(1),GMSB(2),GMSD(1),GMSD(2) WRITE(NSUSYC,351) GMSL(1),GMSL(2),GMSE(1),GMSE(2),GMSN(1) WRITE(NSUSYC,*) WRITE(NSUSYC,*)' MHC CHARG/NEU SLEPTONS SQUARKS' WRITE(NSUSYC,69) WRITE(NSUSYC,*) IF(INDIDEC.NE.0)THEN WRITE(NSUSYCA,70)'MHC ','C1 N1 ','C1 N2 ','C1 N3 ','C1 N4 ' WRITE(NSUSYCA,69) WRITE(NSUSYCA,*) WRITE(NSUSYCB,70)'MHC ','C2 N1 ','C2 N2 ','C2 N3 ','C2 N4 ' WRITE(NSUSYCB,69) WRITE(NSUSYCB,*) WRITE(NSUSYCC,*)' MHC SEL SNL STAU1 SNL STAU2 SNL' WRITE(NSUSYCC,69) WRITE(NSUSYCC,*) WRITE(NSUSYCD,*)' MHC SUL SDL ST1 SB1 '// . 'ST1 SB2 ST2 SB1 ST2 SB2' WRITE(NSUSYCD,69) WRITE(NSUSYCD,*) ENDIF ENDIF 347 FORMAT('TB=',G12.6,1X,'M2=',G12.6,1X,'MU=',G12.6,1X, . 'MSQ=',G12.6) 348 FORMAT('C1=',F7.3,1X,'C2=',F8.3,1X,'N1=',F7.3,1X,'N2=',F7.3,1X, . 'N3=',F8.3,1X,'N4=',F8.3) 349 FORMAT('MST1=',G12.6,1X,'MST2=',G12.6,1X, . 'MSUL=',G12.6,1X,'MSUR=',G12.6) 350 FORMAT('MSB1=',G12.6,1X,'MSB2=',G12.6,1X, . 'MSDL=',G12.6,1X,'MSDR=',G12.6) 351 FORMAT('TAU1=',F8.3,1X,'TAU2=',F8.3,1X,'EL=',F8.3,1X, . 'ER=',F8.3,1X,'NL=',F8.3) C C C ************************************************************** IWRISU = 0 ENDIF ENDIF CALL HDEC(TGBET) IF(IHIGGS.EQ.0)THEN WRITE(NSA,20)AMSM,SMBRB,SMBRL,SMBRM,SMBRS,SMBRC,SMBRT WRITE(NSB,20)AMSM,SMBRG,SMBRGA,SMBRZGA,SMBRW,SMBRZ,SMWDTH ENDIF IF(IHIGGS.EQ.1.OR.IHIGGS.EQ.5)THEN WRITE(NLA,20)AML,HLBRB,HLBRL,HLBRM,HLBRS,HLBRC,HLBRT WRITE(NLB,20)AML,HLBRG,HLBRGA,HLBRZGA,HLBRW,HLBRZ,HLWDTH IF(IOFSUSY.EQ.0)THEN WRITE(NSUSYL,21)AML,HLBRCHT,HLBRNET,HLBRSL,HLBRSQT IF(INDIDEC.NE.0)THEN WRITE(NSUSYLA,23)AML,HLBRSC(1,1),HLBRSC(2,2), . HLBRSC(1,2)+HLBRSC(2,1) WRITE(NSUSYLB,21)AML,HLBRSN(1,1),HLBRSN(2,2),HLBRSN(3,3), . HLBRSN(4,4) WRITE(NSUSYLC,20)AML,HLBRSN(1,2)+HLBRSN(2,1), . HLBRSN(1,3)+HLBRSN(3,1), . HLBRSN(1,4)+HLBRSN(4,1), . HLBRSN(2,3)+HLBRSN(3,2), . HLBRSN(2,4)+HLBRSN(4,2), . HLBRSN(3,4)+HLBRSN(4,3) WRITE(NSUSYLD,20)AML,BHLSLNL,BHLSLEL,BHLSLER,BHLSTAU(1,1), . BHLSTAU(1,2)+BHLSTAU(2,1),BHLSTAU(2,2) WRITE(NSUSYLE,21)AML,BHLSQUL,BHLSQUR,BHLSQDL,BHLSQDR WRITE(NSUSYLF,20)AML,BHLSB(1,1),BHLSB(1,2)+BHLSB(2,1),BHLSB(2,2), . BHLST(1,1),BHLST(1,2)+BHLST(2,1),BHLST(2,2) ENDIF ENDIF ENDIF IF(IHIGGS.EQ.2.OR.IHIGGS.EQ.5)THEN WRITE(NHA,20)AMH,HHBRB,HHBRL,HHBRM,HHBRS,HHBRC,HHBRT WRITE(NHB,20)AMH,HHBRG,HHBRGA,HHBRZGA,HHBRW,HHBRZ WRITE(NHC,20)AMH,HHBRH,HHBRA,HHBRAZ,HHBRHW,HHWDTH IF(IOFSUSY.EQ.0)THEN WRITE(NSUSYH,21)AMH,HHBRCHT,HHBRNET,HHBRSL,HHBRSQT IF(INDIDEC.NE.0)THEN WRITE(NSUSYHA,23)AMH,HHBRSC(1,1),HHBRSC(2,2), . HHBRSC(1,2)+HHBRSC(2,1) WRITE(NSUSYHB,21)AMH,HHBRSN(1,1),HHBRSN(2,2),HHBRSN(3,3), . HHBRSN(4,4) WRITE(NSUSYHC,20)AMH,HHBRSN(1,2)+HHBRSN(2,1), . HHBRSN(1,3)+HHBRSN(3,1), . HHBRSN(1,4)+HHBRSN(4,1), . HHBRSN(2,3)+HHBRSN(3,2), . HHBRSN(2,4)+HHBRSN(4,2), . HHBRSN(3,4)+HHBRSN(4,3) WRITE(NSUSYHD,20)AMH,BHHSLNL,BHHSLEL,BHHSLER,BHHSTAU(1,1), . BHHSTAU(1,2)+BHHSTAU(2,1),BHHSTAU(2,2) WRITE(NSUSYHE,21)AMH,BHHSQUL,BHHSQUR,BHHSQDL,BHHSQDR WRITE(NSUSYHF,20)AMH,BHHSB(1,1),BHHSB(1,2)+BHHSB(2,1),BHHSB(2,2), . BHHST(1,1),BHHST(1,2)+BHHST(2,1),BHHST(2,2) ENDIF ENDIF ENDIF IF(IHIGGS.EQ.3.OR.IHIGGS.EQ.5)THEN WRITE(NAA,20)AMA,ABRB,ABRL,ABRM,ABRS,ABRC,ABRT WRITE(NAB,22)AMA,ABRG,ABRGA,ABRZGA,ABRZ,AWDTH IF(IOFSUSY.EQ.0)THEN WRITE(NSUSYA,21)AMA,HABRCHT,HABRNET,HABRSL,HABRST+HABRSB IF(INDIDEC.NE.0)THEN WRITE(NSUSYAA,23)AMA,HABRSC(1,1),HABRSC(2,2), . HABRSC(1,2)+HABRSC(2,1) WRITE(NSUSYAB,21)AMA,HABRSN(1,1),HABRSN(2,2),HABRSN(3,3), . HABRSN(4,4) WRITE(NSUSYAC,20)AMA,HABRSN(1,2)+HABRSN(2,1), . HABRSN(1,3)+HABRSN(3,1), . HABRSN(1,4)+HABRSN(4,1), . HABRSN(2,3)+HABRSN(3,2), . HABRSN(2,4)+HABRSN(4,2), . HABRSN(3,4)+HABRSN(4,3) WRITE(NSUSYAD,23)AMA,BHASTAU,BHASB,BHAST ENDIF ENDIF ENDIF IF(IHIGGS.EQ.4.OR.IHIGGS.EQ.5)THEN WRITE(NCA,20)AMCH,HCBRB,HCBRL,HCBRM,HCBRS,HCBRC,HCBRT WRITE(NCB,23)AMCH,HCBRW,HCBRA,HCWDTH IF(IOFSUSY.EQ.0)THEN WRITE(NSUSYC,23)AMCH,HCBRCNT,HCBRSL,HCBRSQT IF(INDIDEC.NE.0)THEN WRITE(NSUSYCA,21)AMCH,HCBRSU(1,1),HCBRSU(1,2), . HCBRSU(1,3),HCBRSU(1,4) WRITE(NSUSYCB,21)AMCH,HCBRSU(2,1),HCBRSU(2,2), . HCBRSU(2,3),HCBRSU(2,4) WRITE(NSUSYCC,23)AMCH,BHCSL00,BHCSL11,BHCSL21 WRITE(NSUSYCD,22)AMCH,BHCSQ,BHCSTB(1,1),BHCSTB(1,2), . BHCSTB(2,1),BHCSTB(2,2) ENDIF ENDIF ENDIF 20 FORMAT(G12.6,6(1X,G10.4)) 21 FORMAT(G12.6,4(1X,G10.4)) 22 FORMAT(G12.6,5(1X,G10.4)) 23 FORMAT(G12.6,3(1X,G10.4)) 9999 CONTINUE CLOSE(NI) IF(IHIGGS.EQ.0) THEN CLOSE(NSA) CLOSE(NSB) ENDIF IF(IHIGGS.EQ.1.OR.IHIGGS.EQ.5) THEN CLOSE(NLA) CLOSE(NLB) CLOSE(NSUSYL) ENDIF IF(IHIGGS.EQ.2.OR.IHIGGS.EQ.5) THEN CLOSE(NHA) CLOSE(NHB) CLOSE(NHC) CLOSE(NSUSYH) ENDIF IF(IHIGGS.EQ.3.OR.IHIGGS.EQ.5) THEN CLOSE(NAA) CLOSE(NAB) CLOSE(NSUSYA) ENDIF IF(IHIGGS.EQ.4.OR.IHIGGS.EQ.5) THEN CLOSE(NCA) CLOSE(NCB) CLOSE(NSUSYC) ENDIF STOP END C ===================================================================== C =========== BEGINNING OF THE SUBROUTINE FOR THE DECAYS ============== C !!!!!!!!!!!!!! Any change below this line is at your own risk!!!!!!!! C ===================================================================== SUBROUTINE HDEC(TGBET) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DOUBLE PRECISION LAMB DIMENSION XX(4),YY(4) DIMENSION AMCHAR(2),AMNEUT(4),XMNEUT(4), . AC1(2,2),AC2(2,2),AC3(2,2), . AN1(4,4),AN2(4,4),AN3(4,4), . ACNL(2,4),ACNR(2,4), . AMST(2),AMSB(2),AMSL(2), . AMSU(2),AMSD(2),AMSE(2),AMSN(2), . GLTT(2,2),GLBB(2,2),GLEE(2,2), . GHTT(2,2),GHBB(2,2),GHEE(2,2), . GCTB(2,2),GCEN(2,2) DIMENSION GMST(2),GMSB(2),GMSL(2),GMSU(2),GMSD(2),GMSE(2), . GMSN(2) DIMENSION HLBRSC(2,2),HLBRSN(4,4),HHBRSC(2,2), . HHBRSN(4,4),HABRSC(2,2),HABRSN(4,4),HCBRSU(2,4), . HHBRST(2,2),HHBRSB(2,2),HCBRSTB(2,2) DIMENSION WHLCH(2,2),WHLNE(4,4),WHHCH(2,2),WHHNE(4,4), . WHACH(2,2),WHANE(4,4),WHCCN(2,4), . WHHST(2,2),WHHSB(2,2),WHHSTAU(2,2),WHCSTB(2,2), . WHLST(2,2),WHLSB(2,2),WHLSTAU(2,2) COMPLEX*16 CF,CG,CI1,CI2,CA,CB,CTT,CTB,CTC,CTW,CLT,CLB,CLW, . CAT,CAB,CAC,CAW,CAH,CTH,CLH,CX1,CX2,CAX1,CAX2,CTL,CAL, . CSL,CSQ,CSB1,CSB2,CST1,CST2,CSL1,CSL2, . CXL,CXQ,CXB1,CXB2,CXT1,CXT2,CXL1,CXL2 COMPLEX*16 CSEL,CSER,CSUL,CSUR,CSDL,CSDR, . CXEL,CXER,CXUL,CXUR,CXDL,CXDR COMMON/HMASS/AMSM,AMA,AML,AMH,AMCH,AMAR COMMON/CHIMASS/AMCHI COMMON/MASSES/AMS,AMC,AMB,AMT COMMON/ALS/XLAMBDA,AMC0,AMB0,AMT0,N0 COMMON/PARAM/GF,ALPH,AMTAU,AMMUON,AMZ,AMW COMMON/CKMPAR/VUS,VCB,VUB COMMON/BREAK/AMEL,AMER,AMSQ,AMUR,AMDR,AL,AU,AD,AMU,AM2 COMMON/WZWDTH/GAMC0,GAMT0,GAMT1,GAMW,GAMZ COMMON/ONSHELL/IONSH,IONWZ,IOFSUSY COMMON/OLDFASH/NFGG COMMON/FLAG/IHIGGS,NNLO,IPOLE COMMON/WIDTHSM/SMBRB,SMBRL,SMBRM,SMBRS,SMBRC,SMBRT,SMBRG,SMBRGA, . SMBRZGA,SMBRW,SMBRZ,SMWDTH COMMON/WIDTHA/ABRB,ABRL,ABRM,ABRS,ABRC,ABRT,ABRG,ABRGA,ABRZGA, . ABRZ,AWDTH COMMON/WIDTHHL/HLBRB,HLBRL,HLBRM,HLBRS,HLBRC,HLBRT,HLBRG,HLBRGA, . HLBRZGA,HLBRW,HLBRZ,HLBRA,HLBRAZ,HLBRHW,HLWDTH COMMON/WIDTHHH/HHBRB,HHBRL,HHBRM,HHBRS,HHBRC,HHBRT,HHBRG,HHBRGA, . HHBRZGA,HHBRW,HHBRZ,HHBRH,HHBRA,HHBRAZ,HHBRHW, . HHWDTH COMMON/WIDTHHC/HCBRB,HCBRL,HCBRM,HCBRBU,HCBRS,HCBRC,HCBRT,HCBRW, . HCBRA,HCWDTH COMMON/WISUSY/HLBRSC,HLBRSN,HHBRSC,HHBRSN,HABRSC,HABRSN,HCBRSU, . HLBRCHT,HHBRCHT,HABRCHT,HLBRNET,HHBRNET,HABRNET, . HCBRCNT,HLBRSL,HHBRSL,HCBRSL,HABRSL,HABRST,HABRSB, . HHBRSQ,HHBRST,HHBRSB,HHBRSQT,HCBRSQ,HCBRSTB, . HCBRSQT,HLBRSQ,HLBRSQT COMMON/WISFER/BHLSLNL,BHLSLEL,BHLSLER,BHLSQUL,BHLSQUR,BHLSQDL, . BHLSQDR,BHLST(2,2),BHLSB(2,2),BHLSTAU(2,2), . BHHSLNL,BHHSLEL,BHHSLER,BHHSQUL,BHHSQUR,BHHSQDL, . BHHSQDR,BHHST(2,2),BHHSB(2,2),BHHSTAU(2,2), . BHASTAU,BHASB,BHAST, . BHCSL00,BHCSL11,BHCSL21,BHCSQ,BHCSTB(2,2) COMMON/SMASS/AMNEUT,XMNEUT,AMCHAR,AMST,AMSB,AMSL, . AMSU,AMSD,AMSE,AMSN COMMON/COUP/GAT,GAB,GLT,GLB,GHT,GHB,GZAH,GZAL, . GHHH,GLLL,GHLL,GLHH,GHAA,GLAA,GLVV,GHVV, . GLPM,GHPM,B,A HVV(X,Y)= GF/(4.D0*PI*DSQRT(2.D0))*X**3/2.D0*BETA(Y) . *(1.D0-4.D0*Y+12.D0*Y**2) AFF(X,Y)= GF/(4*PI*DSQRT(2.D0))*X**3*Y*(BETA(Y)) HFF(X,Y)= GF/(4*PI*DSQRT(2.D0))*X**3*Y*(BETA(Y))**3 CFF(Z,TB,X,Y)= GF/(4*PI*DSQRT(2.D0))*Z**3*LAMB(X,Y) . *((1.D0-X-Y)*(X*TB**2+Y/TB**2)-4.D0*X*Y) HV(V)=3.D0*(1.D0-8.D0*V+20.D0*V**2)/DSQRT((4.D0*V-1.D0)) . *DACOS((3.D0*V-1.D0)/2.D0/DSQRT(V**3)) . -(1.D0-V)*(47.D0/2.D0*V-13.D0/2.D0+1.D0/V) . -3.D0/2.D0*(1.D0-6.D0*V+4.D0*V**2)*DLOG(V) HVH(X,Y)=0.25D0*( (1-X)*(-2+4*X-2*X**2+9*Y+9*X*Y-6*Y**2) . /(3*Y)-2*(1-X-X**2+X**3-3*Y-2*X*Y-3*X**2*Y+3*Y**2 . +3*X*Y**2-Y**3)*(-PI/2- DATAN((1-2*X+X**2-Y-X*Y)/ . ((1-X)*DSQRT(-1.D0+2*X+2*Y-(X-Y)**2))))/DSQRT(-1.D0 . +2*X-(X-Y)**2+2*Y)-(1+X**2-2*Y-2*X*Y+Y**2)*DLOG(X)) QCD0(X) = (1+X**2)*(4*SP((1-X)/(1+X)) + 2*SP((X-1)/(X+1)) . - 3*DLOG((1+X)/(1-X))*DLOG(2/(1+X)) . - 2*DLOG((1+X)/(1-X))*DLOG(X)) . - 3*X*DLOG(4/(1-X**2)) - 4*X*DLOG(X) HQCDM(X)=QCD0(X)/X+(3+34*X**2-13*X**4)/16/X**3*DLOG((1+X)/(1-X)) . + 3.D0/8/X**2*(7*X**2-1) AQCDM(X)=QCD0(X)/X + (19+2*X**2+3*X**4)/16/X*DLOG((1+X)/(1-X)) . + 3.D0/8*(7-X**2) HQCD(X)=(4.D0/3*HQCDM(BETA(X)) . +2*(4.D0/3-DLOG(X))*(1-10*X)/(1-4*X))*ASH/PI . + (29.14671D0 + RATCOUP*(1.570D0 - 2*DLOG(HIGTOP)/3 . + DLOG(X)**2/9))*(ASH/PI)**2 . + (164.14D0 - 25.77D0*5 + 0.259D0*5**2)*(ASH/PI)**3 AQCD(X)=(4.D0/3*AQCDM(BETA(X)) . +2*(4.D0/3-DLOG(X))*(1-6*X)/(1-4*X))*ASH/PI . + (29.14671D0 + RATCOUP*(23/6.D0 - DLOG(HIGTOP) . + DLOG(X)**2/6))*(ASH/PI)**2 . + (164.14D0 - 25.77D0*5 + 0.259D0*5**2)*(ASH/PI)**3 QCDH(X)=1.D0+HQCD(X) QCDA(X)=1.D0+AQCD(X) TQCDH(X)=1.D0+4.D0/3*HQCDM(BETA(X))*ASH/PI TQCDA(X)=1.D0+4.D0/3*AQCDM(BETA(X))*ASH/PI QCDC(X,Y)=1.D0+4/3.D0*ASH/PI*(9/4.D0 + (3-2*X+2*Y)/4*DLOG(X/Y) . +((1.5D0-X-Y)*LAMB(X,Y)**2+5*X*Y)/2/LAMB(X,Y) . /(1-X-Y)*DLOG(XI(X,Y)*XI(Y,X)) . + BIJ(X,Y)) . + ASH/PI*(2*(4/3.D0-DLOG(X)) . - (X*2*(4/3.D0-DLOG(X)) + Y*2*(4/3.D0-DLOG(Y)))/(1-X-Y) . - (X*2*(4/3.D0-DLOG(X))*(1-X+Y) . +Y*2*(4/3.D0-DLOG(Y))*(1+X-Y))/LAMB(X,Y)**2) QCDCI(X,Y)=1.D0+4/3.D0*ASH/PI*(3 + (Y-X)/2*DLOG(X/Y) . +(2*(1-X-Y)+LAMB(X,Y)**2)/2/LAMB(X,Y) . *DLOG(XI(X,Y)*XI(Y,X)) . + BIJ(X,Y)) . + ASH/PI*(2*(4/3.D0-DLOG(X)) + 2*(4/3.D0-DLOG(Y)) . - (X*2*(4/3.D0-DLOG(X))*(1-X+Y) . +Y*2*(4/3.D0-DLOG(Y))*(1+X-Y))/LAMB(X,Y)**2) QCDCM(X,Y)=1.D0+4/3.D0*ASH/PI*(9/4.D0 + (3-2*X+2*Y)/4*DLOG(X/Y) . +((1.5D0-X-Y)*LAMB(X,Y)**2+5*X*Y)/2/LAMB(X,Y) . /(1-X-Y)*DLOG(4*X*Y/(1-X-Y+LAMB(X,Y))**2) . + BIJ(X,Y)) QCDCMI(X,Y)=1.D0+4/3.D0*ASH/PI*(3 + (Y-X)/2*DLOG(X/Y) . +(2*(1-X-Y)*LAMB(X,Y)**2)/2/LAMB(X,Y) . *DLOG(4*X*Y/(1-X-Y+LAMB(X,Y))**2) . + BIJ(X,Y)) CQCD(Z,TB,X,Y)= GF/(4*PI*DSQRT(2.D0))*Z**3*LAMB(X,Y) . *((1.D0-X-Y)*(X*TB**2*QCDC(X,Y) . +Y/TB**2*QCDC(Y,X)) . -4.D0*X*Y*QCDCI(X,Y)) CQCDM(Z,TB,X,Y)= GF/(4*PI*DSQRT(2.D0))*Z**3*LAMB(X,Y) . *((1.D0-X-Y)*(X*TB**2*QCDCM(X,Y) . +Y/TB**2*QCDCM(Y,X)) . -4.D0*X*Y*QCDCMI(X,Y)) ELW(AMH,AMF,QF,ACF)=ALPH/PI*3.D0/2*QF**2 . *(3.D0/2-DLOG(AMH**2/AMF**2)) . +GF/8/DSQRT(2.D0)/PI**2*(ACF*AMT**2 . +AMW**2*(3*DLOG(CS)/SS-5)+AMZ**2*(0.5D0 . -3*(1-4*SS*DABS(QF))**2)) CF(CA) = -CDLOG(-(1+CDSQRT(1-CA))/(1-CDSQRT(1-CA)))**2/4 CG(CA) = CDSQRT(1-CA)/2*CDLOG(-(1+CDSQRT(1-CA))/(1-CDSQRT(1-CA))) CI1(CA,CB) = CA*CB/2/(CA-CB) . + CA**2*CB**2/2/(CA-CB)**2*(CF(CA)-CF(CB)) . + CA**2*CB/(CA-CB)**2*(CG(CA)-CG(CB)) CI2(CA,CB) = -CA*CB/2/(CA-CB)*(CF(CA)-CF(CB)) HGGQCD(ASG,NF)=1.D0+ASG/PI*(95.D0/4.D0-NF*7.D0/6.D0) SGGQCD(ASG)=ASG/PI*17.D0/6.D0 AGGQCD(ASG,NF)=1.D0+ASG/PI*(97.D0/4.D0-NF*7.D0/6.D0) HFFSELF(AMH)=1.D0+GF*AMH**2/16.D0/PI**2/DSQRT(2.D0)*2.117203D0 . -(GF*AMH**2/16.D0/PI**2/DSQRT(2.D0))**2*32.6567D0 HVVSELF(AMH)=1.D0+GF*AMH**2/16.D0/PI**2/DSQRT(2.D0)*2.800952D0 . +(GF*AMH**2/16.D0/PI**2/DSQRT(2.D0))**2*62.0308D0 PI=4D0*DATAN(1D0) SS=1.D0-(AMW/AMZ)**2 CS=1.D0-SS C--DECOUPLING THE TOP QUARK FROM ALPHAS AMT0=3.D8 C--TOP QUARK DECAY WIDTH GAMT0 = GF*AMT**3/8/DSQRT(2D0)/PI*(1-AMW**2/AMT**2)**2 . *(1+2*AMW**2/AMT**2) IF(IHIGGS.NE.0.AND.AMT.GT.AMCH+AMB)THEN GAMT1 = GF*AMT**3/8/DSQRT(2D0)/PI*(1-AMCH**2/AMT**2)**2 . *((AMB/AMT)**2*TGBET**2 + 1/TGBET**2) ELSE GAMT1 = 0 ENDIF GAMT1 = GAMT0+GAMT1 IF(IHIGGS.EQ.0)THEN C ========================================================= C SM HIGGS DECAYS C ========================================================= AMXX=AMH AMH=AMSM C ============= RUNNING MASSES RMS = RUNM(AMH,3) RMC = RUNM(AMH,4) RMB = RUNM(AMH,5) RMT = RUNM(AMH,6) RATCOUP = 1 HIGTOP = AMH**2/AMT**2 ASH=ALPHAS(AMH,2) AMC0=1.D8 AMB0=2.D8 AS3=ALPHAS(AMH,2) AMC0=AMC AS4=ALPHAS(AMH,2) AMB0=AMB C AMT0=AMT C =============== PARTIAL WIDTHS C H ---> G G C EPS=1.D-8 NFEXT = 3 ASG = AS3 CTT = 4*AMT**2/AMH**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMH**2*DCMPLX(1D0,-EPS) CAT = 2*CTT*(1+(1-CTT)*CF(CTT)) CAB = 2*CTB*(1+(1-CTB)*CF(CTB)) FQCD=HGGQCD(ASG,NFEXT) XFAC = CDABS(CAT+CAB)**2*FQCD HGG=HVV(AMH,0.D0)*(ASG/PI)**2*XFAC/8 C H ---> G G* ---> G CC TO BE ADDED TO H ---> CC NFEXT = 4 ASG = AS4 FQCD=HGGQCD(ASG,NFEXT) XFAC = CDABS(CAT+CAB)**2*FQCD DCC=HVV(AMH,0.D0)*(ASG/PI)**2*XFAC/8 - HGG C H ---> G G* ---> G BB TO BE ADDED TO H ---> BB NFEXT = 5 ASG = ASH FQCD=HGGQCD(ASG,NFEXT) XFAC = CDABS(CAT+CAB)**2*FQCD DBB=HVV(AMH,0.D0)*(ASG/PI)**2*XFAC/8 - HGG - DCC IF(NFGG.EQ.5)THEN HGG = HGG + DBB + DCC DBB = 0 DCC = 0 ELSEIF(NFGG.EQ.4)THEN HGG = HGG + DCC DCC = 0 ENDIF C H ---> MU MU IF(AMH.LE.2*AMMUON) THEN HMM = 0 ELSE HMM=HFF(AMH,(AMMUON/AMH)**2) . *(1+ELW(AMH,AMMUON,-1.D0,7.D0)) . *HFFSELF(AMH) ENDIF C H ---> TAU TAU IF(AMH.LE.2*AMTAU) THEN HLL = 0 ELSE HLL=HFF(AMH,(AMTAU/AMH)**2) . *(1+ELW(AMH,AMTAU,-1.D0,7.D0)) . *HFFSELF(AMH) ENDIF C H --> SS IF(AMH.LE.2*AMS) THEN HSS = 0 ELSE HS2=3.D0*HFF(AMH,(RMS/AMH)**2) . *QCDH(RMS**2/AMH**2) . *(1+ELW(AMH,RMS,-1.D0/3.D0,7.D0)) . *HFFSELF(AMH) IF(HS2.LT.0.D0) HS2 = 0 HS1=3.D0*HFF(AMH,(AMS/AMH)**2) . *TQCDH(AMS**2/AMH**2) . *HFFSELF(AMH) RAT = 2*AMS/AMH HSS = QQINT(RAT,HS1,HS2) ENDIF C H --> CC IF(AMH.LE.2*AMC) THEN HCC = 0 ELSE HC2=3.D0*HFF(AMH,(RMC/AMH)**2) . *QCDH(RMC**2/AMH**2) . *(1+ELW(AMH,RMC,2.D0/3.D0,7.D0)) . *HFFSELF(AMH) . + DCC IF(HC2.LT.0.D0) HC2 = 0 HC1=3.D0*HFF(AMH,(AMC/AMH)**2) . *TQCDH(AMC**2/AMH**2) . *HFFSELF(AMH) RAT = 2*AMC/AMH HCC = QQINT(RAT,HC1,HC2) ENDIF C H --> BB : IF(AMH.LE.2*AMB) THEN HBB = 0 ELSE HB2=3.D0*HFF(AMH,(RMB/AMH)**2) . *QCDH(RMB**2/AMH**2) . *(1+ELW(AMH,RMB,-1.D0/3.D0,1.D0)) . *HFFSELF(AMH) . + DBB IF(HB2.LT.0.D0) HB2 = 0 HB1=3.D0*HFF(AMH,(AMB/AMH)**2) . *TQCDH(AMB**2/AMH**2) . *HFFSELF(AMH) RAT = 2*AMB/AMH HBB = QQINT(RAT,HB1,HB2) ENDIF C H ---> TT RATCOUP = 0 IF(IONSH.EQ.0)THEN DLD=3D0 DLU=5D0 XM1 = 2D0*AMT-DLD XM2 = 2D0*AMT+DLU IF (AMH.LE.AMT+AMW+AMB) THEN HTT=0.D0 ELSEIF (AMH.LE.XM1) THEN FACTT=6.D0*GF**2*AMH**3*AMT**2/2.D0/128.D0/PI**3 CALL HTOTTS(AMH,AMT,AMB,AMW,HTTS) HTT=FACTT*HTTS ELSEIF (AMH.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 FACTT=6.D0*GF**2*XX(1)**3*AMT**2/2.D0/128.D0/PI**3 CALL HTOTTS(XX(1),AMT,AMB,AMW,HTTS) YY(1)=FACTT*HTTS FACTT=6.D0*GF**2*XX(2)**3*AMT**2/2.D0/128.D0/PI**3 CALL HTOTTS(XX(2),AMT,AMB,AMW,HTTS) YY(2)=FACTT*HTTS XMT = RUNM(XX(3),6) XY2=3.D0*HFF(XX(3),(XMT/XX(3))**2) . *QCDH(XMT**2/XX(3)**2) . *HFFSELF(XX(3)) IF(XY2.LT.0.D0) XY2 = 0 XY1=3.D0*HFF(XX(3),(AMT/XX(3))**2) . *TQCDH(AMT**2/XX(3)**2) . *HFFSELF(XX(3)) RAT = 2*AMT/XX(3) YY(3) = QQINT(RAT,XY1,XY2) XMT = RUNM(XX(4),6) XY2=3.D0*HFF(XX(4),(XMT/XX(4))**2) . *QCDH(XMT**2/XX(4)**2) . *HFFSELF(XX(4)) IF(XY2.LT.0.D0) XY2 = 0 XY1=3.D0*HFF(XX(4),(AMT/XX(4))**2) . *TQCDH(AMT**2/XX(4)**2) . *HFFSELF(XX(4)) RAT = 2*AMT/XX(4) YY(4) = QQINT(RAT,XY1,XY2) HTT = FINT(AMH,XX,YY) ELSE HT2=3.D0*HFF(AMH,(RMT/AMH)**2) . *QCDH(RMT**2/AMH**2) . *HFFSELF(AMH) IF(HT2.LT.0.D0) HT2 = 0 HT1=3.D0*HFF(AMH,(AMT/AMH)**2) . *TQCDH(AMT**2/AMH**2) . *HFFSELF(AMH) RAT = 2*AMT/AMH HTT = QQINT(RAT,HT1,HT2) ENDIF ELSE IF (AMH.LE.2.D0*AMT) THEN HTT=0.D0 ELSE HT2=3.D0*HFF(AMH,(RMT/AMH)**2) . *QCDH(RMT**2/AMH**2) . *HFFSELF(AMH) IF(HT2.LT.0.D0) HT2 = 0 HT1=3.D0*HFF(AMH,(AMT/AMH)**2) . *TQCDH(AMT**2/AMH**2) . *HFFSELF(AMH) RAT = 2*AMT/AMH HTT = QQINT(RAT,HT1,HT2) ENDIF ENDIF C H ---> GAMMA GAMMA EPS=1.D-8 CTT = 4*AMT**2/AMH**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMH**2*DCMPLX(1D0,-EPS) CTC = 4*AMC**2/AMH**2*DCMPLX(1D0,-EPS) CTL = 4*AMTAU**2/AMH**2*DCMPLX(1D0,-EPS) CTW = 4*AMW**2/AMH**2*DCMPLX(1D0,-EPS) CAW = -(2+3*CTW+3*CTW*(2-CTW)*CF(CTW)) CAT = 4/3D0 * 2*CTT*(1+(1-CTT)*CF(CTT)) CAB = 1/3D0 * 2*CTB*(1+(1-CTB)*CF(CTB)) CAC = 4/3D0 * 2*CTC*(1+(1-CTC)*CF(CTC)) CAL = 2*CTL*(1+(1-CTL)*CF(CTL)) XFAC = CDABS(CAT+CAB+CAC+CAL+CAW)**2 HGA=HVV(AMH,0.D0)*(ALPH/PI)**2/16.D0*XFAC C H ---> Z GAMMA IF(AMH.LE.AMZ)THEN HZGA=0 ELSE EPS=1.D-8 TS = SS/CS FT = -3*2D0/3*(1-4*2D0/3*SS)/DSQRT(SS*CS) FB = 3*1D0/3*(-1+4*1D0/3*SS)/DSQRT(SS*CS) CTT = 4*AMT**2/AMH**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMH**2*DCMPLX(1D0,-EPS) CTW = 4*AMW**2/AMH**2*DCMPLX(1D0,-EPS) CLT = 4*AMT**2/AMZ**2*DCMPLX(1D0,-EPS) CLB = 4*AMB**2/AMZ**2*DCMPLX(1D0,-EPS) CLW = 4*AMW**2/AMZ**2*DCMPLX(1D0,-EPS) CAT = FT*(CI1(CTT,CLT) - CI2(CTT,CLT)) CAB = FB*(CI1(CTB,CLB) - CI2(CTB,CLB)) CAW = -1/DSQRT(TS)*(4*(3-TS)*CI2(CTW,CLW) . + ((1+2/CTW)*TS - (5+2/CTW))*CI1(CTW,CLW)) XFAC = CDABS(CAT+CAB+CAW)**2 ACOUP = DSQRT(2D0)*GF*AMZ**2*SS*CS/PI**2 HZGA = GF/(4.D0*PI*DSQRT(2.D0))*AMH**3*(ALPH/PI)*ACOUP/16.D0 . *XFAC*(1-AMZ**2/AMH**2)**3 ENDIF C H ---> W W IF(IONWZ.EQ.0)THEN DLD=2D0 DLU=2D0 XM1 = 2D0*AMW-DLD XM2 = 2D0*AMW+DLU IF (AMH.LE.XM1) THEN CALL HTOVV(AMH,AMW,GAMW,HTWW) HWW = 3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/AMH**3*HTWW ELSEIF (AMH.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 CALL HTOVV(XX(1),AMW,GAMW,HTWW) YY(1)=3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/XX(1)**3*HTWW CALL HTOVV(XX(2),AMW,GAMW,HTWW) YY(2)=3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/XX(2)**3*HTWW YY(3)=HVV(XX(3),AMW**2/XX(3)**2) . *HVVSELF(XX(3)) YY(4)=HVV(XX(4),AMW**2/XX(4)**2) . *HVVSELF(XX(4)) HWW = FINT(AMH,XX,YY) ELSE HWW=HVV(AMH,AMW**2/AMH**2) . *HVVSELF(AMH) ENDIF ELSE DLD=2D0 DLU=2D0 XM1 = 2D0*AMW-DLD XM2 = 2D0*AMW+DLU IF (AMH.LE.AMW) THEN HWW=0 ELSE IF (AMH.LE.XM1) THEN CWW=3.D0*GF**2*AMW**4/16.D0/PI**3 HWW=HV(AMW**2/AMH**2)*CWW*AMH ELSE IF (AMH.LT.XM2) THEN CWW=3.D0*GF**2*AMW**4/16.D0/PI**3 XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1)=HV(AMW**2/XX(1)**2)*CWW*XX(1) YY(2)=HV(AMW**2/XX(2)**2)*CWW*XX(2) YY(3)=HVV(XX(3),AMW**2/XX(3)**2) . *HVVSELF(XX(3)) YY(4)=HVV(XX(4),AMW**2/XX(4)**2) . *HVVSELF(XX(4)) HWW = FINT(AMH,XX,YY) ELSE HWW=HVV(AMH,AMW**2/AMH**2) . *HVVSELF(AMH) ENDIF ENDIF C H ---> Z Z IF(IONWZ.EQ.0)THEN DLD=2D0 DLU=2D0 XM1 = 2D0*AMZ-DLD XM2 = 2D0*AMZ+DLU IF (AMH.LE.XM1) THEN CALL HTOVV(AMH,AMZ,GAMZ,HTZZ) HZZ = 3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/AMH**3*HTZZ ELSEIF (AMH.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 CALL HTOVV(XX(1),AMZ,GAMZ,HTZZ) YY(1)=3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/XX(1)**3*HTZZ CALL HTOVV(XX(2),AMZ,GAMZ,HTZZ) YY(2)=3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/XX(2)**3*HTZZ YY(3)=HVV(XX(3),AMZ**2/XX(3)**2)/2 . *HVVSELF(XX(3)) YY(4)=HVV(XX(4),AMZ**2/XX(4)**2)/2 . *HVVSELF(XX(4)) HZZ = FINT(AMH,XX,YY) ELSE HZZ=HVV(AMH,AMZ**2/AMH**2)/2.D0 . *HVVSELF(AMH) ENDIF ELSE DLD=2D0 DLU=2D0 XM1 = 2D0*AMZ-DLD XM2 = 2D0*AMZ+DLU IF (AMH.LE.AMZ) THEN HZZ=0 ELSE IF (AMH.LE.XM1) THEN CZZ=3.D0*GF**2*AMZ**4/192.D0/PI**3*(7-40/3.D0*SS+160/9.D0*SS**2) HZZ=HV(AMZ**2/AMH**2)*CZZ*AMH ELSE IF (AMH.LT.XM2) THEN CZZ=3.D0*GF**2*AMZ**4/192.D0/PI**3*(7-40/3.D0*SS+160/9.D0*SS**2) XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1)=HV(AMZ**2/XX(1)**2)*CZZ*XX(1) YY(2)=HV(AMZ**2/XX(2)**2)*CZZ*XX(2) YY(3)=HVV(XX(3),AMZ**2/XX(3)**2)/2 . *HVVSELF(XX(3)) YY(4)=HVV(XX(4),AMZ**2/XX(4)**2)/2 . *HVVSELF(XX(4)) HZZ = FINT(AMH,XX,YY) ELSE HZZ=HVV(AMH,AMZ**2/AMH**2)/2.D0 . *HVVSELF(AMH) ENDIF ENDIF C C ========== TOTAL WIDTH AND BRANCHING RATIOS C WTOT=HLL+HMM+HSS+HCC+HBB+HTT+HGG+HGA+HZGA+HWW+HZZ SMBRT=HTT/WTOT SMBRB=HBB/WTOT SMBRL=HLL/WTOT SMBRM=HMM/WTOT SMBRC=HCC/WTOT SMBRS=HSS/WTOT SMBRG=HGG/WTOT SMBRGA=HGA/WTOT SMBRZGA=HZGA/WTOT SMBRW=HWW/WTOT SMBRZ=HZZ/WTOT SMWDTH=WTOT AMH=AMXX endif IF(IHIGGS.GT.0)THEN C +++++++++++++++++++++++ SUSY HIGGSSES +++++++++++++++++++++++ C CALL GAUGINO(AMU,AM2,B,A,AMCHAR,AMNEUT,XMNEUT,AC1,AC2,AC3, . AN1,AN2,AN3,ACNL,ACNR) C CALL SFERMION(AMSQ,AMUR,AMDR,AMEL,AMER,AL,AU,AD,AMU, . AMST,AMSB,AMSL,AMSU,AMSD,AMSE,AMSN, . GLEE,GLTT,GLBB,GHEE,GHTT,GHBB, . GAEE,GATT,GABB,GCEN,GCTB) C ENDIF IF(IHIGGS.EQ.1.OR.IHIGGS.EQ.5)THEN C ========================================================= C LIGHT CP EVEN HIGGS DECAYS C ========================================================= C ============= RUNNING MASSES RMS = RUNM(AML,3) RMC = RUNM(AML,4) RMB = RUNM(AML,5) RMT = RUNM(AML,6) RATCOUP = GLT/GLB HIGTOP = AML**2/AMT**2 ASH=ALPHAS(AML,2) AMC0=1.D8 AMB0=2.D8 C AMT0=3.D8 AS3=ALPHAS(AML,2) AMC0=AMC AS4=ALPHAS(AML,2) AMB0=AMB C AMT0=AMT C =============== PARTIAL WIDTHS C H ---> G G EPS=1.D-8 NFEXT = 3 ASG = AS3 CTT = 4*AMT**2/AML**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AML**2*DCMPLX(1D0,-EPS) CAT = 2*CTT*(1+(1-CTT)*CF(CTT))*GLT CAB = 2*CTB*(1+(1-CTB)*CF(CTB))*GLB CTC = 4*AMC**2/AML**2*DCMPLX(1D0,-EPS) CAC = 2*CTC*(1+(1-CTC)*CF(CTC))*GLT C IF(IOFSUSY.EQ.0) THEN CSB1= 4*AMSB(1)**2/AML**2*DCMPLX(1D0,-EPS) CSB2= 4*AMSB(2)**2/AML**2*DCMPLX(1D0,-EPS) CST1= 4*AMST(1)**2/AML**2*DCMPLX(1D0,-EPS) CST2= 4*AMST(2)**2/AML**2*DCMPLX(1D0,-EPS) CXB1=-AMZ**2/AMSB(1)**2*CSB1*(1-CSB1*CF(CSB1))*GLBB(1,1) CXB2=-AMZ**2/AMSB(2)**2*CSB2*(1-CSB2*CF(CSB2))*GLBB(2,2) CXT1=-AMZ**2/AMST(1)**2*CST1*(1-CST1*CF(CST1))*GLTT(1,1) CXT2=-AMZ**2/AMST(2)**2*CST2*(1-CST2*CF(CST2))*GLTT(2,2) CSUL = 4*AMSU(1)**2/AML**2*DCMPLX(1D0,-EPS) CSUR = 4*AMSU(2)**2/AML**2*DCMPLX(1D0,-EPS) CSDL = 4*AMSD(1)**2/AML**2*DCMPLX(1D0,-EPS) CSDR = 4*AMSD(2)**2/AML**2*DCMPLX(1D0,-EPS) CXUL=2*(1.D0/2.D0-2.D0/3.D0*SS)*AMZ**2/AMSU(1)**2*DSIN(A+B) . *CSUL*(1-CSUL*CF(CSUL)) CXUR=2*(2.D0/3.D0*SS)*AMZ**2/AMSU(2)**2*DSIN(A+B) . *CSUR*(1-CSUR*CF(CSUR)) CXDL=2*(-1.D0/2.D0+1.D0/3.D0*SS)*AMZ**2/AMSD(1)**2*DSIN(A+B) . *CSDL*(1-CSDL*CF(CSDL)) CXDR=2*(-1.D0/3.D0*SS)*AMZ**2/AMSD(2)**2*DSIN(A+B) . *CSDR*(1-CSDR*CF(CSDR)) ELSE CXB1=0.D0 CXB2=0.D0 CXT1=0.D0 CXT2=0.D0 CXUL=0.D0 CXUR=0.D0 CXDL=0.D0 CXDR=0.D0 ENDIF FQCD=HGGQCD(ASG,NFEXT) SQCD=SGGQCD(ASG) XFAC = CDABS(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR)**2*FQCD . + DREAL(DCONJG(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR) . *(CXB1+CXB2+CXT1+CXT2+CXUL+CXUR+CXDL+CXDR))*SQCD HGG=HVV(AML,0.D0)*(ASG/PI)**2*XFAC/8 C H ---> G G* ---> G CC TO BE ADDED TO H ---> CC NFEXT = 4 ASG = AS4 FQCD=HGGQCD(ASG,NFEXT) SQCD=SGGQCD(ASG) XFAC = CDABS(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR)**2*FQCD . + DREAL(DCONJG(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR) . *(CXB1+CXB2+CXT1+CXT2+CXUL+CXUR+CXDL+CXDR))*SQCD DCC=HVV(AML,0.D0)*(ASG/PI)**2*XFAC/8 - HGG C H ---> G G* ---> G BB TO BE ADDED TO H ---> BB NFEXT = 5 ASG = ASH FQCD=HGGQCD(ASG,NFEXT) SQCD=SGGQCD(ASG) XFAC = CDABS(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR)**2*FQCD . + DREAL(DCONJG(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR) . *(CXB1+CXB2+CXT1+CXT2+CXUL+CXUR+CXDL+CXDR))*SQCD DBB=HVV(AML,0.D0)*(ASG/PI)**2*XFAC/8 - HGG - DCC IF(NFGG.EQ.5)THEN HGG = HGG + DBB + DCC DBB = 0 DCC = 0 ELSEIF(NFGG.EQ.4)THEN HGG = HGG + DCC DCC = 0 ENDIF C H ---> MU MU IF(AML.LE.2*AMMUON) THEN HMM = 0 ELSE HMM=HFF(AML,(AMMUON/AML)**2)*GLB**2 ENDIF C H ---> TAU TAU IF(AML.LE.2*AMTAU) THEN HLL = 0 ELSE HLL=HFF(AML,(AMTAU/AML)**2)*GLB**2 ENDIF C H --> SS IF(AML.LE.2*AMS) THEN HSS = 0 ELSE HS1=3.D0*HFF(AML,(AMS/AML)**2) . *GLB**2 . *TQCDH(AMS**2/AML**2) HS2=3.D0*HFF(AML,(RMS/AML)**2)*GLB**2 . *QCDH(RMS**2/AML**2) IF(HS2.LT.0.D0) HS2 = 0 RAT = 2*AMS/AML HSS = QQINT(RAT,HS1,HS2) ENDIF C H --> CC RATCOUP = 1 IF(AML.LE.2*AMC) THEN HCC = 0 ELSE HC1=3.D0*HFF(AML,(AMC/AML)**2) . *GLT**2 . *TQCDH(AMC**2/AML**2) HC2=3.D0*HFF(AML,(RMC/AML)**2)*GLT**2 . *QCDH(RMC**2/AML**2) . + DCC IF(HC2.LT.0.D0) HC2 = 0 RAT = 2*AMC/AML HCC = QQINT(RAT,HC1,HC2) ENDIF C H --> BB : RATCOUP = GLT/GLB IF(AML.LE.2*AMB) THEN HBB = 0 ELSE HB1=3.D0*HFF(AML,(AMB/AML)**2) . *GLB**2 . *TQCDH(AMB**2/AML**2) HB2=3.D0*HFF(AML,(RMB/AML)**2)*GLB**2 . *QCDH(RMB**2/AML**2) . + DBB IF(HB2.LT.0.D0) HB2 = 0 RAT = 2*AMB/AML HBB = QQINT(RAT,HB1,HB2) ENDIF C H ---> TT RATCOUP = 0 IF (AML.LE.2*AMT) THEN HTT=0.D0 ELSE HT1=3.D0*HFF(AML,(AMT/AML)**2)*GLT**2 . *TQCDH(AMT**2/AML**2) HT2=3.D0*HFF(AML,(RMT/AML)**2)*GLT**2 . *QCDH(RMT**2/AML**2) IF(HT2.LT.0.D0) HT2 = 0 RAT = 2*AMT/AML HTT = QQINT(RAT,HT1,HT2) ENDIF C H ---> GAMMA GAMMA EPS=1.D-8 CTT = 4*AMT**2/AML**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AML**2*DCMPLX(1D0,-EPS) CTC = 4*AMC**2/AML**2*DCMPLX(1D0,-EPS) CTL = 4*AMTAU**2/AML**2*DCMPLX(1D0,-EPS) CTW = 4*AMW**2/AML**2*DCMPLX(1D0,-EPS) CTH = 4*AMCH**2/AML**2*DCMPLX(1D0,-EPS) CAT = 4/3D0 * 2*CTT*(1+(1-CTT)*CF(CTT))*GLT CAB = 1/3D0 * 2*CTB*(1+(1-CTB)*CF(CTB))*GLB CAC = 4/3D0 * 2*CTC*(1+(1-CTC)*CF(CTC))*GLT CAL = 1.D0 * 2*CTL*(1+(1-CTL)*CF(CTL))*GLB CAW = -(2+3*CTW+3*CTW*(2-CTW)*CF(CTW))*GLVV CAH = -AMZ**2/2/AMCH**2*CTH*(1-CTH*CF(CTH))*GLPM IF(IOFSUSY.EQ.0) THEN CX1 = 4*AMCHAR(1)**2/AML**2*DCMPLX(1D0,-EPS) CX2 = 4*AMCHAR(2)**2/AML**2*DCMPLX(1D0,-EPS) CSB1= 4*AMSB(1)**2/AML**2*DCMPLX(1D0,-EPS) CSB2= 4*AMSB(2)**2/AML**2*DCMPLX(1D0,-EPS) CST1= 4*AMST(1)**2/AML**2*DCMPLX(1D0,-EPS) CST2= 4*AMST(2)**2/AML**2*DCMPLX(1D0,-EPS) CSL1= 4*AMSL(1)**2/AML**2*DCMPLX(1D0,-EPS) CSL2= 4*AMSL(2)**2/AML**2*DCMPLX(1D0,-EPS) CAX1= AMW/AMCHAR(1) * 2*CX1*(1+(1-CX1)*CF(CX1))*2*AC2(1,1) CAX2= AMW/AMCHAR(2) * 2*CX2*(1+(1-CX2)*CF(CX2))*2*AC2(2,2) CSEL = 4*AMSE(1)**2/AML**2*DCMPLX(1D0,-EPS) CSER = 4*AMSE(2)**2/AML**2*DCMPLX(1D0,-EPS) CSUL = 4*AMSU(1)**2/AML**2*DCMPLX(1D0,-EPS) CSUR = 4*AMSU(2)**2/AML**2*DCMPLX(1D0,-EPS) CSDL = 4*AMSD(1)**2/AML**2*DCMPLX(1D0,-EPS) CSDR = 4*AMSD(2)**2/AML**2*DCMPLX(1D0,-EPS) CXEL=2*(-1/2D0+SS)*AMZ**2/AMSE(1)**2*DSIN(A+B) . *CSEL*(1-CSEL*CF(CSEL)) CXER=-2*(SS)*AMZ**2/AMSE(2)**2*DSIN(A+B) . *CSER*(1-CSER*CF(CSER)) CXUL=2*4.D0/3.D0*(1.D0/2.D0-2.D0/3.D0*SS) . *AMZ**2/AMSU(1)**2*DSIN(A+B)*CSUL*(1-CSUL*CF(CSUL)) CXUR=2*4.D0/3.D0*(2.D0/3.D0*SS) . *AMZ**2/AMSU(2)**2*DSIN(A+B)*CSUR*(1-CSUR*CF(CSUR)) CXDL=2/3.D0*(-1.D0/2.D0+1.D0/3.D0*SS) . *AMZ**2/AMSD(1)**2*DSIN(A+B)*CSDL*(1-CSDL*CF(CSDL)) CXDR=2/3.D0*(-1.D0/3.D0*SS) . *AMZ**2/AMSD(2)**2*DSIN(A+B)*CSDR*(1-CSDR*CF(CSDR)) CXB1=-1/3D0*AMZ**2/AMSB(1)**2*CSB1*(1-CSB1*CF(CSB1))*GLBB(1,1) CXB2=-1/3D0*AMZ**2/AMSB(2)**2*CSB2*(1-CSB2*CF(CSB2))*GLBB(2,2) CXT1=-4/3D0*AMZ**2/AMST(1)**2*CST1*(1-CST1*CF(CST1))*GLTT(1,1) CXT2=-4/3D0*AMZ**2/AMST(2)**2*CST2*(1-CST2*CF(CST2))*GLTT(2,2) CSL1= 4*AMSL(1)**2/AML**2*DCMPLX(1D0,-EPS) CSL2= 4*AMSL(2)**2/AML**2*DCMPLX(1D0,-EPS) CXL1= -AMZ**2/AMSL(1)**2*CSL1*(1-CSL1*CF(CSL1))*GLEE(1,1) CXL2= -AMZ**2/AMSL(2)**2*CSL2*(1-CSL2*CF(CSL2))*GLEE(2,2) XFAC = CDABS(CAT+CAB+CAC+CAL+CAW+CAH+CAX1+CAX2 . + CXEL+CXER+CXUL+CXUR+CXDL+CXDR . + CXB1+CXB2+CXT1+CXT2+CXL1+CXL2)**2 ELSE XFAC = CDABS(CAT+CAB+CAC+CAL+CAW+CAH)**2 ENDIF HGA=HVV(AML,0.D0)*(ALPH/PI)**2/16.D0*XFAC C H ---> Z GAMMA IF(AML.LE.AMZ)THEN HZGA=0 ELSE TS = SS/CS FT = -3*2D0/3*(1-4*2D0/3*SS)/DSQRT(SS*CS)*GLT FB = 3*1D0/3*(-1+4*1D0/3*SS)/DSQRT(SS*CS)*GLB EPS=1.D-8 CTT = 4*AMT**2/AML**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AML**2*DCMPLX(1D0,-EPS) CTW = 4*AMW**2/AML**2*DCMPLX(1D0,-EPS) CTH = 4*AMCH**2/AML**2*DCMPLX(1D0,-EPS) CLT = 4*AMT**2/AMZ**2*DCMPLX(1D0,-EPS) CLB = 4*AMB**2/AMZ**2*DCMPLX(1D0,-EPS) CLW = 4*AMW**2/AMZ**2*DCMPLX(1D0,-EPS) CLH = 4*AMCH**2/AMZ**2*DCMPLX(1D0,-EPS) CAT = FT*(CI1(CTT,CLT) - CI2(CTT,CLT)) CAB = FB*(CI1(CTB,CLB) - CI2(CTB,CLB)) CAW = -1/DSQRT(TS)*(4*(3-TS)*CI2(CTW,CLW) . + ((1+2/CTW)*TS - (5+2/CTW))*CI1(CTW,CLW))*GLVV CAH = (1-2*SS)/DSQRT(SS*CS)*AMZ**2/2/AMCH**2*CI1(CTH,CLH)*GLPM XFAC = CDABS(CAT+CAB+CAW+CAH)**2 ACOUP = DSQRT(2D0)*GF*AMZ**2*SS*CS/PI**2 HZGA = GF/(4.D0*PI*DSQRT(2.D0))*AML**3*(ALPH/PI)*ACOUP/16.D0 . *XFAC*(1-AMZ**2/AML**2)**3 ENDIF C H ---> W W IF(IONWZ.EQ.0)THEN DLD=2D0 DLU=2D0 XM1 = 2D0*AMW-DLD XM2 = 2D0*AMW+DLU IF (AML.LE.XM1) THEN CALL HTOVV(AML,AMW,GAMW,HTWW) HWW = 3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/AML**3*HTWW*GLVV**2 ELSEIF (AML.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 CALL HTOVV(XX(1),AMW,GAMW,HTWW) YY(1)=3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/XX(1)**3*HTWW CALL HTOVV(XX(2),AMW,GAMW,HTWW) YY(2)=3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/XX(2)**3*HTWW YY(3)=HVV(XX(3),AMW**2/XX(3)**2) YY(4)=HVV(XX(4),AMW**2/XX(4)**2) HWW = FINT(AML,XX,YY)*GLVV**2 ELSE HWW=HVV(AML,AMW**2/AML**2)*GLVV**2 ENDIF ELSE DLD=2D0 DLU=2D0 XM1 = 2D0*AMW-DLD XM2 = 2D0*AMW+DLU IF (AML.LE.AMW) THEN HWW=0 ELSE IF (AML.LE.XM1) THEN CWW=3.D0*GF**2*AMW**4/16.D0/PI**3 HWW=HV(AMW**2/AML**2)*CWW*AML*GLVV**2 ELSE IF (AML.LT.XM2) THEN CWW=3.D0*GF**2*AMW**4/16.D0/PI**3 XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1)=HV(AMW**2/XX(1)**2)*CWW*XX(1) YY(2)=HV(AMW**2/XX(2)**2)*CWW*XX(2) YY(3)=HVV(XX(3),AMW**2/XX(3)**2) YY(4)=HVV(XX(4),AMW**2/XX(4)**2) HWW = FINT(AML,XX,YY)*GLVV**2 ELSE HWW=HVV(AML,AMW**2/AML**2)*GLVV**2 ENDIF ENDIF C H ---> Z Z IF(IONWZ.EQ.0)THEN DLD=2D0 DLU=2D0 XM1 = 2D0*AMZ-DLD XM2 = 2D0*AMZ+DLU IF (AML.LE.XM1) THEN CALL HTOVV(AML,AMZ,GAMZ,HTZZ) HZZ = 3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/AML**3*HTZZ*GLVV**2 ELSEIF (AML.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 CALL HTOVV(XX(1),AMZ,GAMZ,HTZZ) YY(1)=3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/XX(1)**3*HTZZ CALL HTOVV(XX(2),AMZ,GAMZ,HTZZ) YY(2)=3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/XX(2)**3*HTZZ YY(3)=HVV(XX(3),AMZ**2/XX(3)**2)/2 YY(4)=HVV(XX(4),AMZ**2/XX(4)**2)/2 HZZ = FINT(AML,XX,YY)*GLVV**2 ELSE HZZ=HVV(AML,AMZ**2/AML**2)/2.D0*GLVV**2 ENDIF ELSE DLD=2D0 DLU=2D0 XM1 = 2D0*AMZ-DLD XM2 = 2D0*AMZ+DLU IF (AML.LE.AMZ) THEN HZZ=0 ELSE IF (AML.LE.XM1) THEN CZZ=3.D0*GF**2*AMZ**4/192.D0/PI**3*(7-40/3.D0*SS+160/9.D0*SS**2) HZZ=HV(AMZ**2/AML**2)*CZZ*AML*GLVV**2 ELSE IF (AML.LT.XM2) THEN CZZ=3.D0*GF**2*AMZ**4/192.D0/PI**3*(7-40/3.D0*SS+160/9.D0*SS**2) XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1)=HV(AMZ**2/XX(1)**2)*CZZ*XX(1) YY(2)=HV(AMZ**2/XX(2)**2)*CZZ*XX(2) YY(3)=HVV(XX(3),AMZ**2/XX(3)**2)/2D0 YY(4)=HVV(XX(4),AMZ**2/XX(4)**2)/2D0 HZZ = FINT(AML,XX,YY)*GLVV**2 ELSE HZZ=HVV(AML,AMZ**2/AML**2)/2.D0*GLVV**2 ENDIF ENDIF C H ---> A A IF (AML.LE.2.D0*AMA) THEN HAA=0 ELSE HAA=GF/16.D0/DSQRT(2D0)/PI*AMZ**4/AML*BETA(AMA**2/AML**2)*GLAA**2 ENDIF C H ---> A Z IF (AML.LE.AMZ+AMA) THEN HAZ=0 ELSE CAZ=LAMB(AMA**2/AML**2,AMZ**2/AML**2) . *LAMB(AML**2/AMZ**2,AMA**2/AMZ**2)**2 HAZ=GF/8.D0/DSQRT(2D0)/PI*AMZ**4/AML*CAZ*GZAL**2 ENDIF C H ---> H+ W+ IF (AML.LE.AMW+AMCH) THEN HHW=0 ELSE CHW=LAMB(AMCH**2/AML**2,AMW**2/AML**2) . *LAMB(AML**2/AMW**2,AMCH**2/AMW**2)**2 HHW=GF/8.D0/DSQRT(2D0)/PI*AMZ**2*AMW**2/AML*CHW*GHVV**2 ENDIF C ============================ SUSY DECAYS IF(IOFSUSY.EQ.0) THEN C C HL ----> CHARGINOS C DO 711 I=1,2 DO 711 J=1,2 IF (AML.GT.AMCHAR(I)+AMCHAR(J)) THEN WHLCH(I,J)=GF*AMW**2/(2*PI*DSQRT(2.D0))/AML . *LAMB(AMCHAR(I)**2/AML**2,AMCHAR(J)**2/AML**2) . *( (AC2(I,J)**2+AC2(J,I)**2)*(AML**2-AMCHAR(I) . **2-AMCHAR(J)**2)-4.D0*AC2(I,J)*AC2(J,I)* . AMCHAR(I)*AMCHAR(J) ) ELSE WHLCH(I,J)=0.D0 ENDIF WHLCHT=WHLCH(1,1)+WHLCH(1,2)+WHLCH(2,1)+WHLCH(2,2) 711 CONTINUE C C HL ----> NEUTRALINOS C DO 712 I=1,4 DO 712 J=1,4 IF (AML.GT.AMNEUT(I)+AMNEUT(J)) THEN WHLNE(I,J)=GF*AMW**2/(2*PI*DSQRT(2.D0))/AML . *AN2(I,J)**2*(AML**2-(XMNEUT(I)+XMNEUT(J))**2) . *LAMB(AMNEUT(I)**2/AML**2,AMNEUT(J)**2/AML**2) ELSE WHLNE(I,J)=0.D0 ENDIF 712 CONTINUE WHLNET= WHLNE(1,1)+WHLNE(1,2)+WHLNE(1,3)+WHLNE(1,4) . +WHLNE(2,1)+WHLNE(2,2)+WHLNE(2,3)+WHLNE(2,4) . +WHLNE(3,1)+WHLNE(3,2)+WHLNE(3,3)+WHLNE(3,4) . +WHLNE(4,1)+WHLNE(4,2)+WHLNE(4,3)+WHLNE(4,4) CCC C HL ----> SLEPTONS C IF (AML.GT.2.D0*AMSE(1)) THEN WHLSLEL=2*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AML*DSIN(B+A)**2 . *BETA(AMSE(1)**2/AML**2)*(-0.5D0+SS)**2 ELSE WHLSLEL=0.D0 ENDIF IF (AML.GT.2.D0*AMSE(2)) THEN WHLSLER=2*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AML*DSIN(B+A)**2 . *BETA(AMSE(2)**2/AML**2)*SS**2 ELSE WHLSLER=0.D0 ENDIF IF (AML.GT.2.D0*AMSN(1)) THEN WHLSLNL=3*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AML*DSIN(B+A)**2 . *BETA(AMSN(1)**2/AML**2)*0.5D0**2 ELSE WHLSLNL=0.D0 ENDIF DO 718 I=1,2 DO 718 J=1,2 IF(AML.GT.AMSL(I)+AMSL(J)) THEN WHLSTAU(I,J)=GF*AMZ**4/2.D0/DSQRT(2.D0)/PI*GLEE(I,J)**2* . LAMB(AMSL(I)**2/AML**2,AMSL(J)**2/AML**2)/AML ELSE WHLSTAU(I,J)=0.D0 ENDIF 718 CONTINUE WHLSLT=WHLSTAU(1,1)+WHLSTAU(2,1)+WHLSTAU(1,2)+WHLSTAU(2,2) . +WHLSLEL+WHLSLER+WHLSLNL C C HL ----> SQUARKS C IF (AML.GT.2.D0*AMSU(1)) THEN WHLSQUL=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AML*DSIN(B+A)**2 . *BETA(AMSU(1)**2/AML**2)*(0.5D0-2.D0/3.D0*SS)**2 ELSE WHLSQUL=0.D0 ENDIF IF (AML.GT.2.D0*AMSU(2)) THEN WHLSQUR=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AML*DSIN(B+A)**2 . *BETA(AMSU(2)**2/AML**2)*(-2.D0/3.D0*SS)**2 ELSE WHLSQUR=0.D0 ENDIF IF (AML.GT.2.D0*AMSD(1)) THEN WHLSQDL=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AML*DSIN(B+A)**2 . *BETA(AMSD(1)**2/AML**2)*(-0.5D0+1.D0/3.D0*SS)**2 ELSE WHLSQDL=0.D0 ENDIF IF (AML.GT.2.D0*AMSD(2)) THEN WHLSQDR=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AML*DSIN(B+A)**2 . *BETA(AMSD(2)**2/AML**2)*(+1.D0/3.D0*SS)**2 ELSE WHLSQDR=0.D0 ENDIF WHLSQ=WHLSQUL+WHLSQUR+WHLSQDL+WHLSQDR C C HL ----> STOPS DO 713 I=1,2 DO 713 J=1,2 IF(AML.GT.AMST(I)+AMST(J)) THEN WHLST(I,J)=3*GF*AMZ**4/2.D0/DSQRT(2.D0)/PI*GLTT(I,J)**2* . LAMB(AMST(I)**2/AML**2,AMST(J)**2/AML**2)/AML ELSE WHLST(I,J)=0.D0 ENDIF 713 CONTINUE C C HL ----> SBOTTOMS DO 714 I=1,2 DO 714 J=1,2 IF(AML.GT.AMSB(I)+AMSB(J)) THEN WHLSB(I,J)=3*GF*AMZ**4/2.D0/DSQRT(2.D0)/PI*GLBB(I,J)**2* . LAMB(AMSB(I)**2/AML**2,AMSB(J)**2/AML**2)/AML ELSE WHLSB(I,J)=0.D0 ENDIF 714 CONTINUE C WHLSTT=WHLST(1,1)+WHLST(1,2)+WHLST(2,1)+WHLST(2,2) WHLSBB=WHLSB(1,1)+WHLSB(1,2)+WHLSB(2,1)+WHLSB(2,2) WHLSQT=WHLSTT+WHLSBB+WHLSQ ELSE WHLCHT=0.D0 WHLNET=0.D0 WHLSLT=0.D0 WHLSQT=0.D0 C--Change thanks to Elzbieta Richter-Was DO I=1,2 DO J=1,2 WHLCH(I,J)=0.D0 WHLST(I,J)=0.D0 WHLSB(I,J)=0.D0 WHLSTAU(I,J)=0.D0 ENDDO ENDDO DO I=1,4 DO J=1,4 WHLNE(I,J)=0.D0 ENDDO ENDDO ENDIF C ========== TOTAL WIDTH AND BRANCHING RATIOS WTOT=HLL+HMM+HSS+HCC+HBB+HTT+HGG+HGA+HZGA+HWW+HZZ+HAA+HAZ+HHW . +WHLCHT+WHLNET+WHLSLT+WHLSQT HLBRT=HTT/WTOT HLBRB=HBB/WTOT HLBRL=HLL/WTOT HLBRM=HMM/WTOT HLBRS=HSS/WTOT HLBRC=HCC/WTOT HLBRG=HGG/WTOT HLBRGA=HGA/WTOT HLBRZGA=HZGA/WTOT HLBRW=HWW/WTOT HLBRZ=HZZ/WTOT HLBRA=HAA/WTOT HLBRAZ=HAZ/WTOT HLBRHW=HHW/WTOT DO 811 I=1,2 DO 811 J=1,2 HLBRSC(I,J)=WHLCH(I,J)/WTOT 811 CONTINUE DO 812 I=1,4 DO 812 J=1,4 HLBRSN(I,J)=WHLNE(I,J)/WTOT 812 CONTINUE HLBRCHT=WHLCHT/WTOT HLBRNET=WHLNET/WTOT HLBRSL=WHLSLT/WTOT HLBRSQ=WHLSQ/WTOT HLBRSQT=WHLSQT/WTOT HLWDTH=WTOT BHLSLNL = WHLSLNL/WTOT BHLSLEL = WHLSLEL/WTOT BHLSLER = WHLSLER/WTOT BHLSQUL = WHLSQUL/WTOT BHLSQUR = WHLSQUR/WTOT BHLSQDL = WHLSQDL/WTOT BHLSQDR = WHLSQDR/WTOT DO I = 1,2 DO J = 1,2 BHLST(I,J) = WHLST(I,J)/WTOT BHLSB(I,J) = WHLSB(I,J)/WTOT BHLSTAU(I,J) = WHLSTAU( I,J)/WTOT ENDDO ENDDO ENDIF IF(IHIGGS.GT.1)THEN C ========================================================= C CHARGED HIGGS DECAYS C ========================================================= TB=TGBET C ============= RUNNING MASSES RMS = RUNM(AMCH,3) RMC = RUNM(AMCH,4) RMB = RUNM(AMCH,5) RMT = RUNM(AMCH,6) ASH=ALPHAS(AMCH,2) C =============== PARTIAL WIDTHS C H+ ---> MU NMU IF(AMCH.LE.AMMUON) THEN HMN = 0 ELSE HMN=CFF(AMCH,TB,(AMMUON/AMCH)**2,0.D0) ENDIF C H+ ---> TAU NTAU IF(AMCH.LE.AMTAU) THEN HLN = 0 ELSE HLN=CFF(AMCH,TB,(AMTAU/AMCH)**2,0.D0) ENDIF C H+ --> SU EPS = 1.D-12 IF(AMCH.LE.AMS+EPS) THEN HSU = 0 ELSE HSU1=3.D0*VUS**2*CQCDM(AMCH,TB,(AMS/AMCH)**2,EPS) HSU2=3.D0*VUS**2*CQCD(AMCH,TB,(RMS/AMCH)**2,EPS) IF(HSU2.LT.0.D0) HSU2 = 0 RAT = AMS/AMCH HSU = QQINT(RAT,HSU1,HSU2) ENDIF C H+ --> CS IF(AMCH.LE.AMS+AMC) THEN HSC = 0 ELSE HSC1=3.D0*CQCDM(AMCH,TB,(AMS/AMCH)**2,(AMC/AMCH)**2) HSC2=3.D0*CQCD(AMCH,TB,(RMS/AMCH)**2,(RMC/AMCH)**2) IF(HSC2.LT.0.D0) HSC2 = 0 RAT = (AMS+AMC)/AMCH HSC = QQINT(RAT,HSC1,HSC2) ENDIF C H+ --> CB IF(AMCH.LE.AMB+AMC) THEN HBC = 0 ELSE HBC1=3.D0*VCB**2*CQCDM(AMCH,TB,(AMB/AMCH)**2,(AMC/AMCH)**2) HBC2=3.D0*VCB**2*CQCD(AMCH,TB,(RMB/AMCH)**2,(RMC/AMCH)**2) IF(HBC2.LT.0.D0) HBC2 = 0 RAT = (AMB+AMC)/AMCH HBC = QQINT(RAT,HBC1,HBC2) ENDIF C H+ --> BU EPS = 1.D-12 IF(AMCH.LE.AMB+EPS) THEN HBU = 0 ELSE HBU1=3.D0*VUB**2*CQCDM(AMCH,TB,(AMB/AMCH)**2,EPS) HBU2=3.D0*VUB**2*CQCD(AMCH,TB,(RMB/AMCH)**2,EPS) IF(HBU2.LT.0.D0) HBU2 = 0 RAT = AMB/AMCH HBU = QQINT(RAT,HBU1,HBU2) ENDIF C H+ --> TB : IF(IONSH.EQ.0)THEN DLD=2D0 DLU=2D0 XM1 = AMT+AMB-DLD XM2 = AMT+AMB+DLU IF (AMCH.LE.AMW+2*AMB) THEN HBT=0.D0 ELSEIF (AMCH.LE.XM1) THEN FACTB=3.D0*GF**2*AMCH*AMT**4/32.D0/PI**3/TB**2 CALL CTOTT(AMCH,AMT,AMB,AMW,CTT0) HBT=FACTB*CTT0 ELSEIF (AMCH.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 FACTB=3.D0*GF**2*XX(1)*AMT**4/32.D0/PI**3/TB**2 CALL CTOTT(XX(1),AMT,AMB,AMW,CTT0) YY(1)=FACTB*CTT0 FACTB=3.D0*GF**2*XX(2)*AMT**4/32.D0/PI**3/TB**2 CALL CTOTT(XX(2),AMT,AMB,AMW,CTT0) YY(2)=FACTB*CTT0 XMB = RUNM(XX(3),5) XMT = RUNM(XX(3),6) XYZ2 = 3.D0*CQCD(XX(3),TB,(XMB/XX(3))**2,(XMT/XX(3))**2) IF(XYZ2.LT.0.D0) XYZ2 = 0 XYZ1 = 3.D0*CQCDM(XX(3),TB,(AMB/XX(3))**2,(AMT/XX(3))**2) RAT = (AMB+AMT)/XX(3) YY(3) = QQINT(RAT,XYZ1,XYZ2) XMB = RUNM(XX(4),5) XMT = RUNM(XX(4),6) XYZ2 = 3.D0*CQCD(XX(4),TB,(XMB/XX(4))**2,(XMT/XX(4))**2) IF(XYZ2.LT.0.D0) XYZ2 = 0 XYZ1 = 3.D0*CQCDM(XX(4),TB,(AMB/XX(4))**2,(AMT/XX(4))**2) RAT = (AMB+AMT)/XX(4) YY(4) = QQINT(RAT,XYZ1,XYZ2) HBT = FINT(AMCH,XX,YY) ELSE HBT2=3.D0*CQCD(AMCH,TB,(RMB/AMCH)**2,(RMT/AMCH)**2) IF(HBT2.LT.0.D0) HBT2 = 0 HBT1=3.D0*CQCDM(AMCH,TB,(AMB/AMCH)**2,(AMT/AMCH)**2) RAT = (AMB+AMT)/AMCH HBT = QQINT(RAT,HBT1,HBT2) ENDIF ELSE IF (AMCH.LE.AMT+AMB) THEN HBT=0.D0 ELSE HBT2=3.D0*CQCD(AMCH,TB,(RMB/AMCH)**2,(RMT/AMCH)**2) IF(HBT2.LT.0.D0) HBT2 = 0 HBT1=3.D0*CQCDM(AMCH,TB,(AMB/AMCH)**2,(AMT/AMCH)**2) RAT = (AMB+AMT)/AMCH HBT = QQINT(RAT,HBT1,HBT2) ENDIF ENDIF C H+ ---> W H IF(IONSH.EQ.0)THEN DLD=3D0 DLU=5D0 XM1 = AMW+AML-DLD XM2 = AMW+AML+DLU IF (AMCH.LT.AML) THEN HWH=0 ELSEIF (AMCH.LE.XM1) THEN IF(AMCH.LE.DABS(AMW-AML))THEN HWH=0 ELSE HWH=9.D0*GF**2/16.D0/PI**3*AMW**4*AMCH*GHVV**2 . *HVH((AML/AMCH)**2,(AMW/AMCH)**2) ENDIF ELSEIF (AMCH.LT.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1) = 9.D0*GF**2/16.D0/PI**3*AMW**4*XX(1) . *HVH((AML/XX(1))**2,(AMW/XX(1))**2) YY(2) = 9.D0*GF**2/16.D0/PI**3*AMW**4*XX(2) . *HVH((AML/XX(2))**2,(AMW/XX(2))**2) CWH=LAMB(AML**2/XX(3)**2,AMW**2/XX(3)**2) . *LAMB(XX(3)**2/AMW**2,AML**2/AMW**2)**2 YY(3)=GF/8.D0/DSQRT(2D0)/PI*AMW**4/XX(3)*CWH CWH=LAMB(AML**2/XX(4)**2,AMW**2/XX(4)**2) . *LAMB(XX(4)**2/AMW**2,AML**2/AMW**2)**2 YY(4)=GF/8.D0/DSQRT(2D0)/PI*AMW**4/XX(4)*CWH HWH = FINT(AMCH,XX,YY)*GHVV**2 ELSE CWH=LAMB(AML**2/AMCH**2,AMW**2/AMCH**2) . *LAMB(AMCH**2/AMW**2,AML**2/AMW**2)**2 HWH=GF/8.D0/DSQRT(2D0)/PI*AMW**4/AMCH*GHVV**2*CWH ENDIF ELSE IF (AMCH.LT.AMW+AML) THEN HWH=0 ELSE CWH=LAMB(AML**2/AMCH**2,AMW**2/AMCH**2) . *LAMB(AMCH**2/AMW**2,AML**2/AMW**2)**2 HWH=GF/8.D0/DSQRT(2D0)/PI*AMW**4/AMCH*GHVV**2*CWH ENDIF ENDIF C H+ ---> W A IF(IONSH.EQ.0)THEN IF (AMCH.LT.AMA) THEN HWA=0 ELSEIF (AMCH.LT.AMW+AMA) THEN IF(AMCH.LE.DABS(AMW-AMA))THEN HWA=0 ELSE HWA=9.D0*GF**2/16.D0/PI**3*AMW**4*AMCH . *HVH((AMA/AMCH)**2,(AMW/AMCH)**2) ENDIF ELSE HWA=0.D0 ENDIF ELSE IF (AMCH.LT.AMW+AMA) THEN HWA=0 ELSE HWA=0.D0 ENDIF ENDIF C ======================= SUSY DECAYS IF(IOFSUSY.EQ.0) THEN C C H+ ----> CHARGINOS+NEUTRALINOS C DO 751 I=1,2 DO 751 J=1,4 IF (AMCH.GT.AMCHAR(I)+AMNEUT(J)) THEN WHCCN(I,J)=GF*AMW**2/(2*PI*DSQRT(2.D0))/AMCH . *LAMB(AMCHAR(I)**2/AMCH**2,AMNEUT(J)**2/AMCH**2)*( . (ACNL(I,J)**2+ACNR(I,J)**2)*(AMCH**2-AMCHAR(I)**2-XMNEUT(J) . **2)-4.D0*ACNL(I,J)*ACNR(I,J)*AMCHAR(I)*XMNEUT(J) ) ELSE WHCCN(I,J)=0.D0 ENDIF 751 CONTINUE WHCCNT=WHCCN(1,1)+WHCCN(1,2)+WHCCN(1,3)+WHCCN(1,4) . +WHCCN(2,1)+WHCCN(2,2)+WHCCN(2,3)+WHCCN(2,4) C C H+ ----> SLEPTONS C IF (AMCH.GT.AMSE(1)+AMSN(1)) THEN WHCSL00=2*GF/4.D0/DSQRT(2D0)/PI*AMW**4/AMCH*DSIN(2.D0*B)**2 . *LAMB(AMSE(1)**2/AMCH**2,AMSN(1)**2/AMCH**2) ELSE WHCSL00=0.D0 ENDIF IF (AMCH.GT.AMSL(1)+AMSN(1)) THEN WHCSL11=GF/4.D0/DSQRT(2D0)/PI*AMW**4/AMCH*GCEN(1,1)**2 . *LAMB(AMSL(1)**2/AMCH**2,AMSN(1)**2/AMCH**2) ELSE WHCSL11=0.D0 ENDIF IF (AMCH.GT.AMSL(2)+AMSN(1)) THEN WHCSL21=GF/4.D0/DSQRT(2D0)/PI*AMW**4/AMCH*GCEN(1,2)**2 . *LAMB(AMSL(2)**2/AMCH**2,AMSN(1)**2/AMCH**2) ELSE WHCSL21=0.D0 ENDIF WHCSLT=WHCSL00+WHCSL11+WHCSL21 C C H+ ----> SQUARKS C IF (AMCH.GT.AMSU(1)+AMSD(1)) THEN WHCSQ=6*GF/4.D0/DSQRT(2D0)/PI*AMW**4/AMCH*DSIN(2.D0*B)**2 . *LAMB(AMSU(1)**2/AMCH**2,AMSD(1)**2/AMCH**2) ELSE WHCSQ=0.D0 ENDIF C DO 753 I=1,2 DO 753 J=1,2 IF(AMCH.GT.AMST(I)+AMSB(J)) THEN WHCSTB(I,J)=3*GF*AMW**4/4.D0/DSQRT(2.D0)/PI*GCTB(I,J)**2 . *LAMB(AMST(I)**2/AMCH**2,AMSB(J)**2/AMCH**2)/AMCH ELSE WHCSTB(I,J)=0.D0 ENDIF 753 CONTINUE C WHCSQT=WHCSQ+WHCSTB(1,1)+WHCSTB(1,2)+WHCSTB(2,1)+WHCSTB(2,2) ELSE WHCCNT=0.D0 WHCSLT=0.D0 WHCSQT=0.D0 C--Change thanks to Elzbieta Richter-Was DO I=1,2 DO J=1,2 WHCSTB(I,J)=0.D0 ENDDO ENDDO DO I=1,2 DO J=1,4 WHCCN(I,J)=0.D0 ENDDO ENDDO ENDIF C C ========== TOTAL WIDTH AND BRANCHING RATIOS C WTOT=HLN+HMN+HSU+HBU+HSC+HBC+HBT+HWH+HWA+WHCCNT+WHCSLT+WHCSQT HCBRL=HLN/WTOT HCBRM=HMN/WTOT HCBRS=HSU/WTOT HCBRBU=HBU/WTOT HCBRC=HSC/WTOT HCBRB=HBC/WTOT HCBRT=HBT/WTOT HCBRW=HWH/WTOT HCBRA=HWA/WTOT DO 851 I=1,2 DO 851 J=1,4 HCBRSU(I,J)=WHCCN(I,J)/WTOT 851 CONTINUE HCBRCNT=WHCCNT/WTOT HCBRSL=WHCSLT/WTOT HCBRSQ=WHCSQ/WTOT HCBRSQT=WHCSQT/WTOT DO 853 I=1,2 DO 853 J=1,2 HCBRSTB(I,J)=WHCSTB(I,J)/WTOT 853 CONTINUE HCWDTH=WTOT BHCSL00 = WHCSL00/WTOT BHCSL11 = WHCSL11/WTOT BHCSL21 = WHCSL21/WTOT BHCSQ = WHCSQ/WTOT DO I = 1,2 DO J = 1,2 BHCSTB(I,J) = WHCSTB(I,J)/WTOT ENDDO ENDDO GAMC0 = WTOT ENDIF IF(IHIGGS.EQ.2.OR.IHIGGS.EQ.5)THEN C ========================================================= C HEAVY CP EVEN HIGGS DECAYS C ========================================================= C ============= RUNNING MASSES RMS = RUNM(AMH,3) RMC = RUNM(AMH,4) RMB = RUNM(AMH,5) RMT = RUNM(AMH,6) RATCOUP = GHT/GHB HIGTOP = AMH**2/AMT**2 ASH=ALPHAS(AMH,2) AMC0=1.D8 AMB0=2.D8 C AMT0=3.D8 AS3=ALPHAS(AMH,2) AMC0=AMC AS4=ALPHAS(AMH,2) AMB0=AMB C AMT0=AMT C =============== PARTIAL WIDTHS C H ---> G G EPS=1.D-8 NFEXT = 3 ASG = AS3 CTT = 4*AMT**2/AMH**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMH**2*DCMPLX(1D0,-EPS) CAT = 2*CTT*(1+(1-CTT)*CF(CTT))*GHT CAB = 2*CTB*(1+(1-CTB)*CF(CTB))*GHB CTC = 4*AMC**2/AMH**2*DCMPLX(1D0,-EPS) CAC = 2*CTC*(1+(1-CTC)*CF(CTC))*GHT C IF(IOFSUSY.EQ.0) THEN CSB1= 4*AMSB(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSB2= 4*AMSB(2)**2/AMH**2*DCMPLX(1D0,-EPS) CST1= 4*AMST(1)**2/AMH**2*DCMPLX(1D0,-EPS) CST2= 4*AMST(2)**2/AMH**2*DCMPLX(1D0,-EPS) C CXB1=-AMZ**2/AMSB(1)**2*CSB1*(1-CSB1*CF(CSB1))*GHBB(1,1) CXB2=-AMZ**2/AMSB(2)**2*CSB2*(1-CSB2*CF(CSB2))*GHBB(2,2) CXT1=-AMZ**2/AMST(1)**2*CST1*(1-CST1*CF(CST1))*GHTT(1,1) CXT2=-AMZ**2/AMST(2)**2*CST2*(1-CST2*CF(CST2))*GHTT(2,2) C CSUL = 4*AMSU(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSUR = 4*AMSU(2)**2/AMH**2*DCMPLX(1D0,-EPS) CSDL = 4*AMSD(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSDR = 4*AMSD(2)**2/AMH**2*DCMPLX(1D0,-EPS) CXUL=-2*(1.D0/2.D0-2.D0/3.D0*SS)*AMZ**2/AMSU(1)**2*DCOS(A+B) . *CSUL*(1-CSUL*CF(CSUL)) CXUR=-2*(2.D0/3.D0*SS)*AMZ**2/AMSU(2)**2*DCOS(A+B) . *CSUR*(1-CSUR*CF(CSUR)) CXDL=-2*(-1.D0/2.D0+1.D0/3.D0*SS)*AMZ**2/AMSD(1)**2*DCOS(A+B) . *CSDL*(1-CSDL*CF(CSDL)) CXDR=-2*(-1.D0/3.D0*SS)*AMZ**2/AMSD(2)**2*DCOS(A+B) . *CSDR*(1-CSDR*CF(CSDR)) ELSE CXB1=0.D0 CXB2=0.D0 CXT1=0.D0 CXT2=0.D0 CXUL=0.D0 CXUR=0.D0 CXDL=0.D0 CXDR=0.D0 ENDIF FQCD=HGGQCD(ASG,NFEXT) SQCD=SGGQCD(ASG) XFAC = CDABS(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR)**2*FQCD . + DREAL(DCONJG(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR) . *(CXB1+CXB2+CXT1+CXT2+CXUL+CXUR+CXDL+CXDR))*SQCD HGG=HVV(AMH,0.D0)*(ASG/PI)**2*XFAC/8 C H ---> G G* ---> G CC TO BE ADDED TO H ---> CC NFEXT = 4 ASG = AS4 FQCD=HGGQCD(ASG,NFEXT) SQCD=SGGQCD(ASG) XFAC = CDABS(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR)**2*FQCD . + DREAL(DCONJG(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR) . *(CXB1+CXB2+CXT1+CXT2+CXUL+CXUR+CXDL+CXDR))*SQCD DCC=HVV(AMH,0.D0)*(ASG/PI)**2*XFAC/8 - HGG C H ---> G G* ---> G BB TO BE ADDED TO H ---> BB NFEXT = 5 ASG = ASH FQCD=HGGQCD(ASG,NFEXT) SQCD=SGGQCD(ASG) XFAC = CDABS(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR)**2*FQCD . + DREAL(DCONJG(CAT+CAB+CAC+CXB1+CXB2+CXT1+CXT2 . +CXUL+CXUR+CXDL+CXDR) . *(CXB1+CXB2+CXT1+CXT2+CXUL+CXUR+CXDL+CXDR))*SQCD DBB=HVV(AMH,0.D0)*(ASG/PI)**2*XFAC/8 - HGG - DCC IF(NFGG.EQ.5)THEN HGG = HGG + DBB + DCC DBB = 0 DCC = 0 ELSEIF(NFGG.EQ.4)THEN HGG = HGG + DCC DCC = 0 ENDIF C H ---> MU MU IF(AMH.LE.2*AMMUON) THEN HMM = 0 ELSE HMM=HFF(AMH,(AMMUON/AMH)**2)*GHB**2 ENDIF C H ---> LL IF(AMH.LE.2*AMTAU) THEN HLL = 0 ELSE HLL=HFF(AMH,(AMTAU/AMH)**2)*GHB**2 ENDIF C H --> SS IF(AMH.LE.2*AMS) THEN HSS = 0 ELSE HS1=3.D0*HFF(AMH,(AMS/AMH)**2) . *GHB**2 . *TQCDH(AMS**2/AMH**2) HS2=3.D0*HFF(AMH,(RMS/AMH)**2)*GHB**2 . *QCDH(RMS**2/AMH**2) IF(HS2.LT.0.D0) HS2 = 0 RAT = 2*AMS/AMH HSS = QQINT(RAT,HS1,HS2) ENDIF C H --> CC RATCOUP = 1 IF(AMH.LE.2*AMC) THEN HCC = 0 ELSE HC1=3.D0*HFF(AMH,(AMC/AMH)**2) . *GHT**2 . *TQCDH(AMC**2/AMH**2) HC2=3.D0*HFF(AMH,(RMC/AMH)**2)*GHT**2 . *QCDH(RMC**2/AMH**2) . + DCC IF(HC2.LT.0.D0) HC2 = 0 RAT = 2*AMC/AMH HCC = QQINT(RAT,HC1,HC2) ENDIF C H --> BB : RATCOUP = GHT/GHB IF(AMH.LE.2*AMB) THEN HBB = 0 ELSE HB1=3.D0*HFF(AMH,(AMB/AMH)**2) . *GHB**2 . *TQCDH(AMB**2/AMH**2) HB2=3.D0*HFF(AMH,(RMB/AMH)**2)*GHB**2 . *QCDH(RMB**2/AMH**2) . + DBB IF(HB2.LT.0.D0) HB2 = 0 RAT = 2*AMB/AMH HBB = QQINT(RAT,HB1,HB2) ENDIF C H ---> TT RATCOUP = 0 IF(IONSH.EQ.0)THEN DLD=3D0 DLU=5D0 XM1 = 2D0*AMT-DLD XM2 = 2D0*AMT+DLU IF (AMH.LE.AMT+AMW+AMB) THEN HTT=0.D0 ELSEIF (AMH.LE.XM1) THEN FACTT=6.D0*GF**2*AMH**3*AMT**2/2.D0/128.D0/PI**3 CALL HTOTT(AMH,AMT,AMB,AMW,AMCH,TB,GHT,GAT,GHVV,HTT0) HTT=FACTT*HTT0 ELSEIF (AMH.LE.XM2) THEN ZZMA=AMAR XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 CALL AMHAMA(2,XX(1),TGBET) FACTT=6.D0*GF**2*XX(1)**3*AMT**2/2.D0/128.D0/PI**3 CALL HTOTT(XX(1),AMT,AMB,AMW,AMCH,TB,GHT,GAT,GHVV,HTT0) YY(1)=FACTT*HTT0 CALL AMHAMA(2,XX(2),TGBET) FACTT=6.D0*GF**2*XX(2)**3*AMT**2/2.D0/128.D0/PI**3 CALL HTOTT(XX(2),AMT,AMB,AMW,AMCH,TB,GHT,GAT,GHVV,HTT0) YY(2)=FACTT*HTT0 CALL AMHAMA(2,XX(3),TGBET) XMT = RUNM(XX(3),6) HT1=3.D0*HFF(XX(3),(AMT/XX(3))**2)*GHT**2 . *TQCDH(AMT**2/XX(3)**2) HT2=3.D0*HFF(XX(3),(XMT/XX(3))**2)*GHT**2 . *QCDH(XMT**2/XX(3)**2) IF(HT2.LT.0.D0) HT2 = 0 RAT = 2*AMT/XX(3) YY(3) = QQINT(RAT,HT1,HT2) CALL AMHAMA(2,XX(4),TGBET) XMT = RUNM(XX(4),6) HT1=3.D0*HFF(XX(4),(AMT/XX(4))**2)*GHT**2 . *TQCDH(AMT**2/XX(4)**2) HT2=3.D0*HFF(XX(4),(XMT/XX(4))**2)*GHT**2 . *QCDH(XMT**2/XX(4)**2) IF(HT2.LT.0.D0) HT2 = 0 RAT = 2*AMT/XX(4) YY(4) = QQINT(RAT,HT1,HT2) AMA = ZZMA CALL SUSYCP(TGBET) HTT=FINT(AMH,XX,YY) ELSE HT1=3.D0*HFF(AMH,(AMT/AMH)**2)*GHT**2 . *TQCDH(AMT**2/AMH**2) HT2=3.D0*HFF(AMH,(RMT/AMH)**2)*GHT**2 . *QCDH(RMT**2/AMH**2) IF(HT2.LT.0.D0) HT2 = 0 RAT = 2*AMT/AMH HTT = QQINT(RAT,HT1,HT2) ENDIF ELSE IF (AMH.LE.2.D0*AMT) THEN HTT=0.D0 ELSE HT1=3.D0*HFF(AMH,(AMT/AMH)**2)*GHT**2 . *TQCDH(AMT**2/AMH**2) HT2=3.D0*HFF(AMH,(RMT/AMH)**2)*GHT**2 . *QCDH(RMT**2/AMH**2) IF(HT2.LT.0.D0) HT2 = 0 RAT = 2*AMT/AMH HTT = QQINT(RAT,HT1,HT2) ENDIF ENDIF C H ---> GAMMA GAMMA EPS=1.D-8 CTT = 4*AMT**2/AMH**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMH**2*DCMPLX(1D0,-EPS) CTL = 4*AMTAU**2/AMH**2*DCMPLX(1D0,-EPS) CTW = 4*AMW**2/AMH**2*DCMPLX(1D0,-EPS) CTH = 4*AMCH**2/AMH**2*DCMPLX(1D0,-EPS) CTC = 4*AMC**2/AMH**2*DCMPLX(1D0,-EPS) CAC = 4/3D0 * 2*CTC*(1+(1-CTC)*CF(CTC))*GHT CAT = 4/3D0 * 2*CTT*(1+(1-CTT)*CF(CTT))*GHT CAB = 1/3D0 * 2*CTB*(1+(1-CTB)*CF(CTB))*GHB CAL = 1.D0 * 2*CTL*(1+(1-CTL)*CF(CTL))*GHB CAW = -(2+3*CTW+3*CTW*(2-CTW)*CF(CTW))*GHVV CAH = -AMZ**2/2/AMCH**2*CTH*(1-CTH*CF(CTH))*GHPM IF(IOFSUSY.EQ.0) THEN CX1 = 4*AMCHAR(1)**2/AMH**2*DCMPLX(1D0,-EPS) CX2 = 4*AMCHAR(2)**2/AMH**2*DCMPLX(1D0,-EPS) CAX1= AMW/AMCHAR(1) * 2*CX1*(1+(1-CX1)*CF(CX1))*2*AC1(1,1) CAX2= AMW/AMCHAR(2) * 2*CX2*(1+(1-CX2)*CF(CX2))*2*AC1(2,2) CSL1= 4*AMSL(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSL2= 4*AMSL(2)**2/AMH**2*DCMPLX(1D0,-EPS) CSB1= 4*AMSB(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSB2= 4*AMSB(2)**2/AMH**2*DCMPLX(1D0,-EPS) CST1= 4*AMST(1)**2/AMH**2*DCMPLX(1D0,-EPS) CST2= 4*AMST(2)**2/AMH**2*DCMPLX(1D0,-EPS) CSEL = 4*AMSE(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSER = 4*AMSE(2)**2/AMH**2*DCMPLX(1D0,-EPS) CSUL = 4*AMSU(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSUR = 4*AMSU(2)**2/AMH**2*DCMPLX(1D0,-EPS) CSDL = 4*AMSD(1)**2/AMH**2*DCMPLX(1D0,-EPS) CSDR = 4*AMSD(2)**2/AMH**2*DCMPLX(1D0,-EPS) CXEL=-2*(-1/2D0+SS)*AMZ**2/AMSE(1)**2*DCOS(A+B) . *CSEL*(1-CSEL*CF(CSEL)) CXER=2*(SS)*AMZ**2/AMSE(2)**2*DCOS(A+B) . *CSER*(1-CSER*CF(CSER)) CXUL=-2*4.D0/3.D0*(1.D0/2.D0-2.D0/3.D0*SS) . *AMZ**2/AMSU(1)**2*DCOS(A+B)*CSUL*(1-CSUL*CF(CSUL)) CXUR=-2*4.D0/3.D0*(2.D0/3.D0*SS) . *AMZ**2/AMSU(2)**2*DCOS(A+B)*CSUR*(1-CSUR*CF(CSUR)) CXDL=-2/3.D0*(-1.D0/2.D0+1.D0/3.D0*SS) . *AMZ**2/AMSD(1)**2*DCOS(A+B)*CSDL*(1-CSDL*CF(CSDL)) CXDR=-2/3.D0*(-1.D0/3.D0*SS) . *AMZ**2/AMSD(2)**2*DCOS(A+B)*CSDR*(1-CSDR*CF(CSDR)) CXB1= -1/3D0*AMZ**2/AMSB(1)**2*CSB1*(1-CSB1*CF(CSB1))*GHBB(1,1) CXB2= -1/3D0*AMZ**2/AMSB(2)**2*CSB2*(1-CSB2*CF(CSB2))*GHBB(2,2) CXT1= -4/3D0*AMZ**2/AMST(1)**2*CST1*(1-CST1*CF(CST1))*GHTT(1,1) CXT2= -4/3D0*AMZ**2/AMST(2)**2*CST2*(1-CST2*CF(CST2))*GHTT(2,2) CXL1= -AMZ**2/AMSL(1)**2*CSL1*(1-CSL1*CF(CSL1))*GHEE(1,1) CXL2= -AMZ**2/AMSL(2)**2*CSL2*(1-CSL2*CF(CSL2))*GHEE(2,2) XFAC = CDABS(CAT+CAB+CAC+CAL+CAW+CAH+CAX1+CAX2 . + CXEL+CXER+CXUL+CXUR+CXDL+CXDR . + CXB1+CXB2+CXT1+CXT2+CXL1+CXL2)**2 ELSE XFAC = CDABS(CAT+CAB+CAC+CAL+CAW+CAH)**2 ENDIF HGA=HVV(AMH,0.D0)*(ALPH/PI)**2/16.D0*XFAC C H ---> Z GAMMA IF(AMH.LE.AMZ)THEN HZGA=0 ELSE TS = SS/CS FT = -3*2D0/3*(1-4*2D0/3*SS)/DSQRT(SS*CS)*GHT FB = 3*1D0/3*(-1+4*1D0/3*SS)/DSQRT(SS*CS)*GHB EPS=1.D-8 CTT = 4*AMT**2/AMH**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMH**2*DCMPLX(1D0,-EPS) CTW = 4*AMW**2/AMH**2*DCMPLX(1D0,-EPS) CTH = 4*AMCH**2/AMH**2*DCMPLX(1D0,-EPS) CLT = 4*AMT**2/AMZ**2*DCMPLX(1D0,-EPS) CLB = 4*AMB**2/AMZ**2*DCMPLX(1D0,-EPS) CLW = 4*AMW**2/AMZ**2*DCMPLX(1D0,-EPS) CLH = 4*AMCH**2/AMZ**2*DCMPLX(1D0,-EPS) CAT = FT*(CI1(CTT,CLT) - CI2(CTT,CLT)) CAB = FB*(CI1(CTB,CLB) - CI2(CTB,CLB)) CAW = -1/DSQRT(TS)*(4*(3-TS)*CI2(CTW,CLW) . + ((1+2/CTW)*TS - (5+2/CTW))*CI1(CTW,CLW))*GHVV CAH = (1-2*SS)/DSQRT(SS*CS)*AMZ**2/2/AMCH**2*CI1(CTH,CLH)*GHPM XFAC = CDABS(CAT+CAB+CAW+CAH)**2 ACOUP = DSQRT(2D0)*GF*AMZ**2*SS*CS/PI**2 HZGA = GF/(4.D0*PI*DSQRT(2.D0))*AMH**3*(ALPH/PI)*ACOUP/16.D0 . *XFAC*(1-AMZ**2/AMH**2)**3 ENDIF C H ---> W W IF(IONWZ.EQ.0)THEN DLD=2D0 DLU=2D0 XM1 = 2D0*AMW-DLD XM2 = 2D0*AMW+DLU IF (AMH.LE.XM1) THEN CALL HTOVV(AMH,AMW,GAMW,HTWW) HWW = 3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/AMH**3*HTWW*GHVV**2 ELSEIF (AMH.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 CALL HTOVV(XX(1),AMW,GAMW,HTWW) YY(1)=3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/XX(1)**3*HTWW CALL HTOVV(XX(2),AMW,GAMW,HTWW) YY(2)=3D0/2D0*GF*AMW**4/DSQRT(2D0)/PI/XX(2)**3*HTWW YY(3)=HVV(XX(3),AMW**2/XX(3)**2) YY(4)=HVV(XX(4),AMW**2/XX(4)**2) HWW = FINT(AMH,XX,YY)*GHVV**2 ELSE HWW=HVV(AMH,AMW**2/AMH**2)*GHVV**2 ENDIF ELSE DLD=2D0 DLU=2D0 XM1 = 2D0*AMW-DLD XM2 = 2D0*AMW+DLU IF (AMH.LE.AMW) THEN HWW=0 ELSE IF (AMH.LE.XM1) THEN CWW=3.D0*GF**2*AMW**4/16.D0/PI**3 HWW=HV(AMW**2/AMH**2)*CWW*AMH*GHVV**2 ELSE IF (AMH.LT.XM2) THEN CWW=3.D0*GF**2*AMW**4/16.D0/PI**3 XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1)=HV(AMW**2/XX(1)**2)*CWW*XX(1) YY(2)=HV(AMW**2/XX(2)**2)*CWW*XX(2) YY(3)=HVV(XX(3),AMW**2/XX(3)**2) YY(4)=HVV(XX(4),AMW**2/XX(4)**2) HWW = FINT(AMH,XX,YY)*GHVV**2 ELSE HWW=HVV(AMH,AMW**2/AMH**2)*GHVV**2 ENDIF ENDIF C H ---> Z Z IF(IONWZ.EQ.0)THEN DLD=2D0 DLU=2D0 XM1 = 2D0*AMZ-DLD XM2 = 2D0*AMZ+DLU IF (AMH.LE.XM1) THEN CALL HTOVV(AMH,AMZ,GAMZ,HTZZ) HZZ = 3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/AMH**3*HTZZ*GHVV**2 ELSEIF (AMH.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 CALL HTOVV(XX(1),AMZ,GAMZ,HTZZ) YY(1)=3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/XX(1)**3*HTZZ CALL HTOVV(XX(2),AMZ,GAMZ,HTZZ) YY(2)=3D0/4D0*GF*AMZ**4/DSQRT(2D0)/PI/XX(2)**3*HTZZ YY(3)=HVV(XX(3),AMZ**2/XX(3)**2)/2 YY(4)=HVV(XX(4),AMZ**2/XX(4)**2)/2 HZZ = FINT(AMH,XX,YY)*GHVV**2 ELSE HZZ=HVV(AMH,AMZ**2/AMH**2)/2.D0*GHVV**2 ENDIF ELSE DLD=2D0 DLU=2D0 XM1 = 2D0*AMZ-DLD XM2 = 2D0*AMZ+DLU IF (AMH.LE.AMZ) THEN HZZ=0 ELSE IF (AMH.LE.XM1) THEN CZZ=3.D0*GF**2*AMZ**4/192.D0/PI**3*(7-40/3.D0*SS+160/9.D0*SS**2) HZZ=HV(AMZ**2/AMH**2)*CZZ*AMH*GHVV**2 ELSE IF (AMH.LT.XM2) THEN CZZ=3.D0*GF**2*AMZ**4/192.D0/PI**3*(7-40/3.D0*SS+160/9.D0*SS**2) XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1)=HV(AMZ**2/XX(1)**2)*CZZ*XX(1) YY(2)=HV(AMZ**2/XX(2)**2)*CZZ*XX(2) YY(3)=HVV(XX(3),AMZ**2/XX(3)**2)/2D0 YY(4)=HVV(XX(4),AMZ**2/XX(4)**2)/2D0 HZZ = FINT(AMH,XX,YY)*GHVV**2 ELSE HZZ=HVV(AMH,AMZ**2/AMH**2)/2.D0*GHVV**2 ENDIF ENDIF C H ---> h h IF(IONSH.EQ.0)THEN ZZMA = AMAR AMREAL = AMH AMA = 20.D0 CALL SUSYCP(TGBET) AMLOW = AMH AMDEL = AMREAL - AMLOW DLD = 0.3D0*(TGBET-1.3D0) DLD = DMAX1(0.1D0,DLD) DLU=DLD AMA = ZZMA CALL SUSYCP(TGBET) XM1 = 2*AML-DLD XM2 = 2*AML+DLU IF (AMH.LE.AML) THEN HHH=0 ELSEIF (AMH.LT.XM1) THEN XH=AML**2/AMH**2 XH1=(XH-1.D0)*(2.D0-.5D0*DLOG(XH))+(1.D0-5.D0*XH) . *(DATAN((2.D0*XH-1.D0)/DSQRT(4.D0*XH-1.D0)) . -DATAN(1.D0/DSQRT(4.D0*XH-1.D0)))/DSQRT(4.D0*XH-1.D0) XH2=3*GF**2/16.D0/PI**3*AMZ**4/AMH*GHLL**2*GLB**2*AMB**2 HHH=XH1*XH2 ELSEIF (AMH.LT.XM2) THEN IFLON0 = 0 IFLON1 = 0 ZZMA=AMAR AMACRIT = AMAR AMA0 = AMAR AMA1 = AMAR 510 AMA0 = AMA0 - 1 AMA1 = AMA1 + 1 AMA = AMA0 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML) THEN IFLON0 = -1 ELSE IFLON0 = 1 ENDIF AMA = AMA1 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML) THEN IFLON1 = -1 ELSE IFLON1 = 1 ENDIF IF(IFLON0*IFLON1.NE.-1) GOTO 510 501 AMA = (AMA0+AMA1)/2 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML) THEN IF(IFLON0.EQ.-1) THEN AMA0 = AMAR ELSE AMA1 = AMAR ENDIF ELSE IF(IFLON0.EQ.-1) THEN AMA1 = AMAR ELSE AMA0 = AMAR ENDIF ENDIF AMACRIT = (AMA0+AMA1)/2 DEL = 1.D-8 AMDEL = 2*DABS(AMA1-AMA0)/(AMA1+AMA0) IF(AMDEL.GT.DEL) GOTO 501 AMA = AMACRIT CALL SUSYCP(TGBET) YM1 = AMACRIT YM2 = AMACRIT AMA0 = AMACRIT AMA1 = AMACRIT DELSTEP = 1.D0 511 AMA0 = AMA0 - DELSTEP AMA1 = AMA1 + DELSTEP AMA = AMACRIT CALL SUSYCP(TGBET) IF(AMH.LT.2*AML-DLD) THEN IFLONC = -1 ELSE IFLONC = 1 ENDIF AMA = AMA0 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML-DLD) THEN IFLON0 = -1 ELSE IFLON0 = 1 ENDIF AMA = AMA1 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML-DLD) THEN IFLON1 = -1 ELSE IFLON1 = 1 ENDIF IF(IFLON0*IFLONC.NE.-1.AND.IFLONC*IFLON1.NE.-1) GOTO 511 IF(IFLON0*IFLONC.EQ.-1) THEN AMA1 = AMACRIT IFLON1 = IFLONC ELSE AMA0 = AMACRIT IFLON0 = IFLONC ENDIF 512 AMA = (AMA0+AMA1)/2 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML-DLD) THEN IF(IFLON0.EQ.-1) THEN AMA0 = AMAR ELSE AMA1 = AMAR ENDIF ELSE IF(IFLON0.EQ.-1) THEN AMA1 = AMAR ELSE AMA0 = AMAR ENDIF ENDIF YM1 = (AMA0+AMA1)/2 DEL = 1.D-8 AMDEL = 2*DABS(AMA1-AMA0)/(AMA1+AMA0) IF(AMDEL.GT.DEL) GOTO 512 AMA = YM1 CALL SUSYCP(TGBET) AMA0 = AMACRIT AMA1 = AMACRIT DELSTEP = 1.D0 513 AMA0 = AMA0 - DELSTEP AMA1 = AMA1 + DELSTEP AMA = AMACRIT CALL SUSYCP(TGBET) IF(AMH.LT.2*AML+DLU) THEN IFLONC = -1 ELSE IFLONC = 1 ENDIF AMA = AMA0 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML+DLU) THEN IFLON0 = -1 ELSE IFLON0 = 1 ENDIF AMA = AMA1 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML+DLU) THEN IFLON1 = -1 ELSE IFLON1 = 1 ENDIF IF(IFLON0*IFLONC.NE.-1.AND.IFLONC*IFLON1.NE.-1) GOTO 513 IF(IFLON0*IFLONC.EQ.-1) THEN AMA1 = AMACRIT IFLON1 = IFLONC ELSE AMA0 = AMACRIT IFLON0 = IFLONC ENDIF 514 AMA = (AMA0+AMA1)/2 CALL SUSYCP(TGBET) IF(AMH.LT.2*AML+DLU) THEN IF(IFLON0.EQ.-1) THEN AMA0 = AMAR ELSE AMA1 = AMAR ENDIF ELSE IF(IFLON0.EQ.-1) THEN AMA1 = AMAR ELSE AMA0 = AMAR ENDIF ENDIF YM2 = (AMA0+AMA1)/2 DEL = 1.D-8 AMDEL = 2*DABS(AMA1-AMA0)/(AMA1+AMA0) IF(AMDEL.GT.DEL) GOTO 514 AMA = YM2 CALL SUSYCP(TGBET) DEL = 1.D-4 XX(1) = YM1 - DEL XX(2) = YM1 XX(3) = YM2 XX(4) = YM2 + DEL AMAR = ZZMA DO J=1,4 AMA = XX(J) CALL SUSYCP(TGBET) XX(J) = AMH IF(AMH.GE.2*AML)THEN YY(J)=GF/16D0/DSQRT(2D0)/PI*AMZ**4/XX(J) . *BETA(AML**2/XX(J)**2) ELSEIF(AMH.LE.AML)THEN YY(J) = 0 ELSE XH=AML**2/XX(J)**2 XH1=(XH-1.D0)*(2.D0-.5D0*DLOG(XH))+(1.D0-5.D0*XH) . *(DATAN((2.D0*XH-1.D0)/DSQRT(4.D0*XH-1.D0)) . -DATAN(1.D0/DSQRT(4.D0*XH-1.D0)))/DSQRT(4.D0*XH-1.D0) XH2=3*GF**2/16.D0/PI**3*AMZ**4/XX(J)*GLB**2*AMB**2 YY(J)=XH1*XH2 ENDIF ENDDO AMA = ZZMA CALL SUSYCP(TGBET) HHH = FINT(AMH,XX,YY)*GHLL**2 ELSE HHH=GF/16D0/DSQRT(2D0)/PI*AMZ**4/AMH*BETA(AML**2/AMH**2) . *GHLL**2 ENDIF ELSE IF (AMH.LE.2*AML) THEN HHH=0 ELSE HHH=GF/16D0/DSQRT(2D0)/PI*AMZ**4/AMH*BETA(AML**2/AMH**2) . *GHLL**2 ENDIF ENDIF C H ---> A A IF(IONSH.EQ.0)THEN DLD = 0.3D0*(TGBET-1.3D0) DLD = DMAX1(0.1D0,DLD) DLU=DLD ALD = DLD/2 ALU = DLU/2 XM1 = 2*AMA-DLD XM2 = 2*AMA+DLU IF (AMH.LE.AMA) THEN HAA=0 ELSEIF (AMH.LT.XM1) THEN XA=AMA**2/AMH**2 XA1=(XA-1.D0)*(2.D0-.5D0*DLOG(XA))+(1.D0-5.D0*XA) . *(DATAN((2.D0*XA-1.D0)/DSQRT(4.D0*XA-1.D0)) . -DATAN(1.D0/DSQRT(4.D0*XA-1.D0)))/DSQRT(4.D0*XA-1.D0) XA2=3*GF**2/16.D0/PI**3*AMZ**4/AMH*GHAA**2*GAB**2*AMB**2 HAA=XA1*XA2 ELSEIF (AMH.LT.XM2) THEN ZZMA=AMAR AMACRIT = AMAR AMA0 = 10.D0 AMA1 = AMAR + 50.D0 AMA = AMA0 CALL SUSYCP(TGBET) IF(AMH.LT.2*AMA) THEN IFLON0 = -1 ELSEIF(AMH.EQ.2*AMA) THEN IFLON0 = 0 AMACRIT = AMAR ELSE IFLON0 = 1 ENDIF AMA = AMA1 CALL SUSYCP(TGBET) IF(AMH.LT.2*AMA) THEN IFLON1 = -1 ELSEIF(AMH.EQ.2*AMA) THEN IFLON1 = 0 AMACRIT = AMAR ELSE IFLON1 = 1 ENDIF IF(IFLON0*IFLON1.EQ.0)THEN IFLON0 = 0 IFLON1 = 0 ENDIF IF(IFLON0.NE.IFLON1)THEN 502 AMA = (AMA0+AMA1)/2 CALL SUSYCP(TGBET) IF(AMH.LT.2*AMA) THEN IF(IFLON0.EQ.-1) THEN AMA0 = AMAR ELSE AMA1 = AMAR ENDIF ELSEIF(AMH.EQ.2*AMA) THEN IFLON0 = 0 IFLON1 = 0 AMACRIT = AMAR ELSE IF(IFLON0.EQ.-1) THEN AMA1 = AMAR ELSE AMA0 = AMAR ENDIF ENDIF IF(IFLON0.NE.0)THEN AMACRIT = (AMA0+AMA1)/2 DEL = 1.D-8 AMDEL = 2*DABS(AMA1-AMA0)/(AMA1+AMA0) IF(AMDEL.GT.DEL) GOTO 502 ENDIF ENDIF DEL = 1.D-4 XX(1) = AMACRIT - ALD - DEL XX(2) = AMACRIT - ALD XX(3) = AMACRIT + ALU XX(4) = AMACRIT + ALU + DEL DO J=1,4 AMA = XX(J) CALL SUSYCP(TGBET) XX(J) = AMH IF(AMH.GE.2*AMA)THEN YY(J)=GF/16D0/DSQRT(2D0)/PI*AMZ**4/XX(J) . *BETA(AMA**2/XX(J)**2) ELSEIF(AMH.LE.AMA)THEN YY(J) = 0 ELSE XA=AMA**2/XX(J)**2 XA1=(XA-1.D0)*(2.D0-.5D0*DLOG(XA))+(1.D0-5.D0*XA) . *(DATAN((2.D0*XA-1.D0)/DSQRT(4.D0*XA-1.D0)) . -DATAN(1.D0/DSQRT(4.D0*XA-1.D0)))/DSQRT(4.D0*XA-1.D0) XA2=3*GF**2/16.D0/PI**3*AMZ**4/XX(J)*GAB**2*AMB**2 YY(J)=XA1*XA2 ENDIF ENDDO AMA = ZZMA CALL SUSYCP(TGBET) HAA = FINT(AMH,XX,YY)*GHAA**2 ELSE HAA=GF/16D0/DSQRT(2D0)/PI*AMZ**4/AMH*BETA(AMA**2/AMH**2) . *GHAA**2 ENDIF ELSE IF (AMH.LE.2*AMA) THEN HAA=0 ELSE HAA=GF/16D0/DSQRT(2D0)/PI*AMZ**4/AMH*BETA(AMA**2/AMH**2) . *GHAA**2 ENDIF ENDIF C H ---> A Z IF(IONSH.EQ.0)THEN DLD=1D0 DLU=8D0 XM1 = AMA+AMZ-DLD XM2 = AMA+AMZ+DLU IF (AMH.LT.AMA) THEN HAZ=0 ELSEIF (AMH.LT.XM1) THEN IF(AMH.LE.DABS(AMZ-AMA))THEN HAZ=0 ELSE HAZ=9.D0*GF**2/8.D0/PI**3*AMZ**4*AMH*GZAH**2* . (7.D0/12.D0-10.D0/9.D0*SS+40.D0/27.D0*SS**2) . *HVH((AMA/AMH)**2,(AMZ/AMH)**2) ENDIF ELSEIF (AMH.LT.XM2) THEN ZZMA=AMAR 165 AMA = AMAR - 1.D0 CALL SUSYCP(TGBET) IF(AMH.LT.AMA+AMZ+DLU.AND.AMH.GT.AMA+AMZ-DLD) GOTO 165 XX(1) = AMAR-1D0 XX(2) = AMAR AMA = ZZMA CALL SUSYCP(TGBET) 166 AMA = AMAR + 1.D0 CALL SUSYCP(TGBET) IF(AMH.LT.AMA+AMZ+DLU.AND.AMH.GT.AMA+AMZ-DLD) GOTO 166 XX(3) = AMAR XX(4) = AMAR+1D0 DO IJ=1,4 AMA = XX(IJ) CALL SUSYCP(TGBET) XX(IJ) = AMH IF(AMH.LE.AMA+AMZ) THEN YY(IJ)=9.D0*GF**2/8.D0/PI**3*AMZ**4*XX(IJ)* . (7.D0/12.D0-10.D0/9.D0*SS+40.D0/27.D0*SS**2) . *HVH((AMA/XX(IJ))**2,(AMZ/XX(IJ))**2) ELSE CAZ=LAMB(AMA**2/XX(IJ)**2,AMZ**2/XX(IJ)**2) . *LAMB(XX(IJ)**2/AMZ**2,AMA**2/AMZ**2)**2 YY(IJ)=GF/8.D0/DSQRT(2D0)/PI*AMZ**4/XX(IJ)*CAZ ENDIF ENDDO AMA = ZZMA CALL SUSYCP(TGBET) HAZ = FINT(AMH,XX,YY)*GZAH**2 ELSE CAZ=LAMB(AMA**2/AMH**2,AMZ**2/AMH**2) . *LAMB(AMH**2/AMZ**2,AMA**2/AMZ**2)**2 HAZ=GF/8.D0/DSQRT(2D0)/PI*AMZ**4/AMH*CAZ*GZAH**2 ENDIF ELSE IF (AMH.LT.AMZ+AMA) THEN HAZ=0 ELSE CAZ=LAMB(AMA**2/AMH**2,AMZ**2/AMH**2) . *LAMB(AMH**2/AMZ**2,AMA**2/AMZ**2)**2 HAZ=GF/8.D0/DSQRT(2D0)/PI*AMZ**4/AMH*CAZ*GZAH**2 ENDIF ENDIF C H ---> H+ W+ IF(IONSH.EQ.0)THEN DLD=3D0 DLU=9D0 XM1 = AMCH+AMW-DLD XM2 = AMCH+AMW+DLU IF (AMH.LT.AMCH) THEN HHW=0.D0 ELSEIF (AMH.LT.XM1) THEN IF(AMH.LE.DABS(AMW-AMCH))THEN HHW=0 ELSE HHW=9.D0*GF**2/8.D0/PI**3*AMW**4*AMH*GLVV**2*2 . *HVH((AMCH/AMH)**2,(AMW/AMH)**2) ENDIF ELSEIF (AMH.LT.XM2) THEN ZZMA=AMAR 167 AMA = AMAR - 1.D0 CALL SUSYCP(TGBET) IF(AMH.LT.AMCH+AMW+DLU) GOTO 167 XX(1) = AMAR-1D0 XX(2) = AMAR AMA = ZZMA CALL SUSYCP(TGBET) 168 AMA = AMAR + 1.D0 CALL SUSYCP(TGBET) IF(AMH.GT.AMCH+AMW-DLD) GOTO 168 XX(3) = AMAR XX(4) = AMAR+1D0 AMA = XX(1) CALL SUSYCP(TGBET) XX(1) = AMH CHW=LAMB(AMCH**2/XX(1)**2,AMW**2/XX(1)**2) . *LAMB(XX(1)**2/AMW**2,AMCH**2/AMW**2)**2 YY(1)=2*GF/8.D0/DSQRT(2D0)/PI*AMW**4/XX(1)*CHW AMA = XX(2) CALL SUSYCP(TGBET) XX(2) = AMH CHW=LAMB(AMCH**2/XX(2)**2,AMW**2/XX(2)**2) . *LAMB(XX(2)**2/AMW**2,AMCH**2/AMW**2)**2 YY(2)=2*GF/8.D0/DSQRT(2D0)/PI*AMW**4/XX(2)*CHW AMA = XX(3) CALL SUSYCP(TGBET) XX(3) = AMH YY(3)=9.D0*GF**2/8.D0/PI**3*AMW**4*XX(3)*2 . *HVH((AMCH/XX(3))**2,(AMW/XX(3))**2) AMA = XX(4) CALL SUSYCP(TGBET) XX(4) = AMH YY(4)=9.D0*GF**2/8.D0/PI**3*AMW**4*XX(4)*2 . *HVH((AMCH/XX(4))**2,(AMW/XX(4))**2) AMA = ZZMA CALL SUSYCP(TGBET) HHW=FINT(AMH,XX,YY)*GLVV**2 ELSE CHW=LAMB(AMCH**2/AMH**2,AMW**2/AMH**2) . *LAMB(AMH**2/AMW**2,AMCH**2/AMW**2)**2 HHW=2*GF/8.D0/DSQRT(2D0)/PI*AMW**4/AMH*CHW*GLVV**2 ENDIF ELSE IF (AMH.LT.AMW+AMCH) THEN HHW=0.D0 ELSE CHW=LAMB(AMCH**2/AMH**2,AMW**2/AMH**2) . *LAMB(AMH**2/AMW**2,AMCH**2/AMW**2)**2 HHW=2*GF/8.D0/DSQRT(2D0)/PI*AMW**4/AMH*CHW*GLVV**2 ENDIF ENDIF C ========================== SUSY DECAYS C IF(IOFSUSY.EQ.0) THEN C HH ----> CHARGINOS DO 741 I=1,2 DO 741 J=1,2 IF (AMH.GT.AMCHAR(I)+AMCHAR(J)) THEN WHHCH(I,J)=GF*AMW**2/(2*PI*DSQRT(2.D0))/AMH . *LAMB(AMCHAR(I)**2/AMH**2,AMCHAR(J)**2/AMH**2) . *( (AC1(I,J)**2+AC1(J,I)**2)*(AMH**2-AMCHAR(I) . **2-AMCHAR(J)**2)-4.D0*AC1(I,J)*AC1(J,I)* . AMCHAR(I)*AMCHAR(J) ) ELSE WHHCH(I,J)=0.D0 ENDIF 741 CONTINUE WHHCHT=WHHCH(1,1)+WHHCH(1,2)+WHHCH(2,1)+WHHCH(2,2) C C HH ----> NEUTRALINOS DO 742 I=1,4 DO 742 J=1,4 IF (AMH.GT.AMNEUT(I)+AMNEUT(J)) THEN WHHNE(I,J)=GF*AMW**2/(2*PI*DSQRT(2.D0))/AMH . *AN1(I,J)**2*(AMH**2-(XMNEUT(I)+XMNEUT(J))**2) . *LAMB(AMNEUT(I)**2/AMH**2,AMNEUT(J)**2/AMH**2) ELSE WHHNE(I,J)=0.D0 ENDIF 742 CONTINUE WHHNET= WHHNE(1,1)+WHHNE(1,2)+WHHNE(1,3)+WHHNE(1,4) . +WHHNE(2,1)+WHHNE(2,2)+WHHNE(2,3)+WHHNE(2,4) . +WHHNE(3,1)+WHHNE(3,2)+WHHNE(3,3)+WHHNE(3,4) . +WHHNE(4,1)+WHHNE(4,2)+WHHNE(4,3)+WHHNE(4,4) C C HH ----> SLEPTONS C IF (AMH.GT.2.D0*AMSE(1)) THEN WHHSLEL=2*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AMH*DCOS(B+A)**2 . *BETA(AMSE(1)**2/AMH**2)*(-0.5D0+SS)**2 ELSE WHHSLEL=0.D0 ENDIF IF (AMH.GT.2.D0*AMSE(2)) THEN WHHSLER=2*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AMH*DCOS(B+A)**2 . *BETA(AMSE(2)**2/AMH**2)*SS**2 ELSE WHHSLER=0.D0 ENDIF IF (AMH.GT.2.D0*AMSN(1)) THEN WHHSLNL=3*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AMH*DCOS(B+A)**2 . *BETA(AMSN(1)**2/AMH**2)*0.5D0**2 ELSE WHHSLNL=0.D0 ENDIF DO 748 I=1,2 DO 748 J=1,2 IF(AMH.GT.AMSL(I)+AMSL(J)) THEN WHHSTAU(I,J)=GF*AMZ**4/2.D0/DSQRT(2.D0)/PI*GHEE(I,J)**2* . LAMB(AMSL(I)**2/AMH**2,AMSL(J)**2/AMH**2)/AMH ELSE WHHSTAU(I,J)=0.D0 ENDIF 748 CONTINUE WHHSLT=WHHSTAU(1,1)+WHHSTAU(1,2)+WHHSTAU(2,1)+WHHSTAU(2,2) . +WHHSLEL+WHHSLER+WHHSLNL C C HH ----> SQUARKS C IF (AMH.GT.2.D0*AMSU(1)) THEN WHHSQUL=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AMH*DCOS(B+A)**2 . *BETA(AMSU(1)**2/AMH**2)*(0.5D0-2.D0/3.D0*SS)**2 ELSE WHHSQUL=0.D0 ENDIF IF (AMH.GT.2.D0*AMSU(2)) THEN WHHSQUR=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AMH*DCOS(B+A)**2 . *BETA(AMSU(2)**2/AMH**2)*(-2.D0/3.D0*SS)**2 ELSE WHHSQUR=0.D0 ENDIF IF (AMH.GT.2.D0*AMSD(1)) THEN WHHSQDL=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AMH*DCOS(B+A)**2 . *BETA(AMSD(1)**2/AMH**2)*(-0.5D0+1.D0/3.D0*SS)**2 ELSE WHHSQDL=0.D0 ENDIF IF (AMH.GT.2.D0*AMSD(2)) THEN WHHSQDR=6*GF/2.D0/DSQRT(2D0)/PI*AMZ**4/AMH*DCOS(B+A)**2 . *BETA(AMSD(2)**2/AMH**2)*(+1.D0/3.D0*SS)**2 ELSE WHHSQDR=0.D0 ENDIF WHHSQ=WHHSQUL+WHHSQUR+WHHSQDL+WHHSQDR C C HH ----> STOPS DO 743 I=1,2 DO 743 J=1,2 IF(AMH.GT.AMST(I)+AMST(J)) THEN WHHST(I,J)=3*GF*AMZ**4/2.D0/DSQRT(2.D0)/PI*GHTT(I,J)**2* . LAMB(AMST(I)**2/AMH**2,AMST(J)**2/AMH**2)/AMH ELSE WHHST(I,J)=0.D0 ENDIF 743 CONTINUE C C HH ----> SBOTTOMS DO 744 I=1,2 DO 744 J=1,2 IF(AMH.GT.AMSB(I)+AMSB(J)) THEN WHHSB(I,J)=3*GF*AMZ**4/2.D0/DSQRT(2.D0)/PI*GHBB(I,J)**2* . LAMB(AMSB(I)**2/AMH**2,AMSB(J)**2/AMH**2)/AMH ELSE WHHSB(I,J)=0.D0 ENDIF 744 CONTINUE C WHHSTT=WHHST(1,1)+WHHST(1,2)+WHHST(2,1)+WHHST(2,2) WHHSBB=WHHSB(1,1)+WHHSB(1,2)+WHHSB(2,1)+WHHSB(2,2) WHHSQT=WHHSTT+WHHSBB+WHHSQ C ELSE WHHCHT=0.D0 WHHNET=0.D0 WHHSLT=0.D0 WHHSQT=0.D0 C--Change thanks to Elzbieta Richter-Was DO I=1,2 DO J=1,2 WHHCH(I,J)=0.D0 WHHST(I,J)=0.D0 WHHSB(I,J)=0.D0 WHHSTAU(I,J)=0.D0 ENDDO ENDDO DO I=1,4 DO J=1,4 WHHNE(I,J)=0.D0 ENDDO ENDDO ENDIF C C ========== TOTAL WIDTH AND BRANCHING RATIOS WTOT=HLL+HMM+HSS+HCC+HBB+HTT+HGG+HGA+HZGA+HWW+HZZ+HHH+HAA+HAZ . +HHW+WHHCHT+WHHNET+WHHSLT+WHHSQT HHBRT=HTT/WTOT HHBRB=HBB/WTOT HHBRL=HLL/WTOT HHBRM=HMM/WTOT HHBRS=HSS/WTOT HHBRC=HCC/WTOT HHBRG=HGG/WTOT HHBRGA=HGA/WTOT HHBRZGA=HZGA/WTOT HHBRW=HWW/WTOT HHBRZ=HZZ/WTOT HHBRH=HHH/WTOT HHBRA=HAA/WTOT HHBRAZ=HAZ/WTOT HHBRHW=HHW/WTOT DO 841 I=1,2 DO 841 J=1,2 HHBRSC(I,J)=WHHCH(I,J)/WTOT 841 CONTINUE DO 842 I=1,4 DO 842 J=1,4 HHBRSN(I,J)=WHHNE(I,J)/WTOT 842 CONTINUE HHBRCHT=WHHCHT/WTOT HHBRNET=WHHNET/WTOT HHBRSL=WHHSLT/WTOT HHBRSQ=WHHSQ/WTOT HHBRSQT=WHHSQT/WTOT DO 843 I=1,2 DO 843 J=1,2 HHBRST(I,J)=WHHST(I,J)/WTOT 843 CONTINUE DO 844 I=1,2 DO 844 J=1,2 HHBRSB(I,J)=WHHSB(I,J)/WTOT 844 CONTINUE HHWDTH=WTOT BHHSLNL = WHHSLNL/WTOT BHHSLEL = WHHSLEL/WTOT BHHSLER = WHHSLER/WTOT BHHSQUL = WHHSQUL/WTOT BHHSQUR = WHHSQUR/WTOT BHHSQDL = WHHSQDL/WTOT BHHSQDR = WHHSQDR/WTOT DO I = 1,2 DO J = 1,2 BHHST(I,J) = WHHST(I,J)/WTOT BHHSB(I,J) = WHHSB(I,J)/WTOT BHHSTAU(I,J) = WHHSTAU( I,J)/WTOT ENDDO ENDDO ENDIF IF(IHIGGS.EQ.3.OR.IHIGGS.EQ.5)THEN C C ========================================================= C CP ODD HIGGS DECAYS C ========================================================= C ============= RUNNING MASSES RMS = RUNM(AMA,3) RMC = RUNM(AMA,4) RMB = RUNM(AMA,5) RMT = RUNM(AMA,6) RATCOUP = GAT/GAB HIGTOP = AMA**2/AMT**2 ASH=ALPHAS(AMA,2) AMC0=1.D8 AMB0=2.D8 C AMT0=3.D8 AS3=ALPHAS(AMA,2) AMC0=AMC AS4=ALPHAS(AMA,2) AMB0=AMB C AMT0=AMT C =============== PARTIAL WIDTHS C A ---> G G EPS=1.D-8 NFEXT = 3 ASG = AS3 CTT = 4*AMT**2/AMA**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMA**2*DCMPLX(1D0,-EPS) CAT = CTT*CF(CTT)*GAT CAB = CTB*CF(CTB)*GAB CTC = 4*AMC**2/AMA**2*DCMPLX(1D0,-EPS) CAC = CTC*CF(CTC)*GAT FQCD=AGGQCD(ASG,NFEXT) XFAC = CDABS(CAT+CAB+CAC)**2*FQCD HGG=GF/(16.D0*PI*DSQRT(2.D0))*AMA**3*(ASG/PI)**2*XFAC C A ---> G G* ---> G CC TO BE ADDED TO A ---> CC NFEXT = 4 ASG = AS4 FQCD=AGGQCD(ASG,NFEXT) XFAC = CDABS(CAT+CAB+CAC)**2*FQCD DCC=GF/(16.D0*PI*DSQRT(2.D0))*AMA**3*(ASG/PI)**2*XFAC . - HGG C A ---> G G* ---> G BB TO BE ADDED TO A ---> BB NFEXT = 5 ASG = ASH FQCD=AGGQCD(ASG,NFEXT) XFAC = CDABS(CAT+CAB+CAC)**2*FQCD DBB=GF/(16.D0*PI*DSQRT(2.D0))*AMA**3*(ASG/PI)**2*XFAC . - HGG - DCC IF(NFGG.EQ.5)THEN HGG = HGG + DBB + DCC DBB = 0 DCC = 0 ELSEIF(NFGG.EQ.4)THEN HGG = HGG + DCC DCC = 0 ENDIF C A ---> MU MU IF(AMA.LE.2*AMMUON) THEN HMM = 0 ELSE HMM=AFF(AMA,(AMMUON/AMA)**2)*GAB**2 ENDIF C A ---> LL IF(AMA.LE.2*AMTAU) THEN HLL = 0 ELSE HLL=AFF(AMA,(AMTAU/AMA)**2)*GAB**2 ENDIF C A --> SS IF(AMA.LE.2*AMS) THEN HSS = 0 ELSE HS1=3.D0*HFF(AMA,(AMS/AMA)**2) . *GAB**2 . *TQCDA(AMS**2/AMA**2) HS2=3.D0*AFF(AMA,(RMS/AMA)**2) . *GAB**2 . *QCDA(RMS**2/AMA**2) IF(HS2.LT.0.D0) HS2 = 0 RAT = 2*AMS/AMA HSS = QQINT(RAT,HS1,HS2) ENDIF C A --> CC RATCOUP = 1 IF(AMA.LE.2*AMC) THEN HCC = 0 ELSE HC1=3.D0*HFF(AMA,(AMC/AMA)**2) . *GAT**2 . *TQCDA(AMC**2/AMA**2) HC2=3.D0*HFF(AMA,(RMC/AMA)**2) . *GAT**2 . *QCDA(RMC**2/AMA**2) . + DCC IF(HC2.LT.0.D0) HC2 = 0 RAT = 2*AMC/AMA HCC = QQINT(RAT,HC1,HC2) ENDIF C A --> BB : RATCOUP = GAT/GAB IF(AMA.LE.2*AMB) THEN HBB = 0 ELSE HB1=3.D0*HFF(AMA,(AMB/AMA)**2) . *GAB**2 . *TQCDA(AMB**2/AMA**2) HB2=3.D0*HFF(AMA,(RMB/AMA)**2) . *GAB**2 . *QCDA(RMB**2/AMA**2) . + DBB IF(HB2.LT.0.D0) HB2 = 0 RAT = 2*AMB/AMA HBB = QQINT(RAT,HB1,HB2) ENDIF C A --> TT : RATCOUP = 0 IF(IONSH.EQ.0)THEN DLD=3D0 DLU=4D0 XM1 = 2D0*AMT-DLD XM2 = 2D0*AMT+DLU IF (AMA.LE.AMT+AMW+AMB) THEN HTT=0.D0 ELSEIF (AMA.LE.XM1) THEN FACTT=6.D0*GF**2*AMA**3*AMT**2/2.D0/128.D0/PI**3*GAT**2 CALL ATOTT(AMA,AMT,AMB,AMW,AMCH,ATT0) HTT=FACTT*ATT0 ELSEIF (AMA.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 FACTT=6.D0*GF**2*XX(1)**3*AMT**2/2.D0/128.D0/PI**3 CALL ATOTT(XX(1),AMT,AMB,AMW,AMCH,ATT0) YY(1)=FACTT*ATT0 FACTT=6.D0*GF**2*XX(2)**3*AMT**2/2.D0/128.D0/PI**3 CALL ATOTT(XX(2),AMT,AMB,AMW,AMCH,ATT0) YY(2)=FACTT*ATT0 XMT = RUNM(XX(3),6) XYZ1 =3.D0*AFF(XX(3),(AMT/XX(3))**2) . *TQCDA(AMT**2/XX(3)**2) XYZ2 =3.D0*AFF(XX(3),(XMT/XX(3))**2) . *QCDA(XMT**2/XX(3)**2) IF(XYZ2.LT.0.D0) XYZ2 = 0 RAT = 2*AMT/XX(3) YY(3) = QQINT(RAT,XYZ1,XYZ2) XMT = RUNM(XX(4),6) XYZ1 =3.D0*AFF(XX(4),(AMT/XX(4))**2) . *TQCDA(AMT**2/XX(4)**2) XYZ2 =3.D0*AFF(XX(4),(XMT/XX(4))**2) . *QCDA(XMT**2/XX(4)**2) IF(XYZ2.LT.0.D0) XYZ2 = 0 RAT = 2*AMT/XX(4) YY(4) = QQINT(RAT,XYZ1,XYZ2) HTT = FINT(AMA,XX,YY)*GAT**2 ELSE HT1=3.D0*AFF(AMA,(AMT/AMA)**2)*GAT**2 . *TQCDA(AMT**2/AMA**2) HT2=3.D0*AFF(AMA,(RMT/AMA)**2)*GAT**2 . *QCDA(RMT**2/AMA**2) IF(HT2.LT.0.D0) HT2 = 0 RAT = 2*AMT/AMA HTT = QQINT(RAT,HT1,HT2) ENDIF ELSE IF (AMA.LE.2.D0*AMT) THEN HTT=0.D0 ELSE HT1=3.D0*AFF(AMA,(AMT/AMA)**2)*GAT**2 . *TQCDA(AMT**2/AMA**2) HT2=3.D0*AFF(AMA,(RMT/AMA)**2)*GAT**2 . *QCDA(RMT**2/AMA**2) IF(HT2.LT.0.D0) HT2 = 0 RAT = 2*AMT/AMA HTT = QQINT(RAT,HT1,HT2) ENDIF ENDIF C A ---> GAMMA GAMMA EPS=1.D-8 CTT = 4*AMT**2/AMA**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMA**2*DCMPLX(1D0,-EPS) CAT = 4/3D0 * CTT*CF(CTT)*GAT CAB = 1/3D0 * CTB*CF(CTB)*GAB CTC = 4*AMC**2/AMA**2*DCMPLX(1D0,-EPS) CAC = 4/3D0 * CTC*CF(CTC)*GAT CTL = 4*AMTAU**2/AMA**2*DCMPLX(1D0,-EPS) CAL = 1.D0 * CTL*CF(CTL)*GAB IF(IOFSUSY.EQ.0) THEN CX1 = 4*AMCHAR(1)**2/AMA**2*DCMPLX(1D0,-EPS) CX2 = 4*AMCHAR(2)**2/AMA**2*DCMPLX(1D0,-EPS) CAX1= AMW/AMCHAR(1) * CX1*CF(CX1) * 2*AC3(1,1) CAX2= AMW/AMCHAR(2) * CX2*CF(CX2) * 2*AC3(2,2) XFAC = CDABS(CAT+CAB+CAC+CAL+CAX1+CAX2)**2 ELSE XFAC = CDABS(CAT+CAB+CAC+CAL)**2 ENDIF HGA=GF/(32.D0*PI*DSQRT(2.D0))*AMA**3*(ALPH/PI)**2*XFAC C A ---> Z GAMMA IF(AMA.LE.AMZ)THEN HZGA=0 ELSE TS = SS/CS FT = -3*2D0/3*(1-4*2D0/3*SS)/DSQRT(SS*CS)*GAT FB = 3*1D0/3*(-1+4*1D0/3*SS)/DSQRT(SS*CS)*GAB EPS=1.D-8 CTT = 4*AMT**2/AMA**2*DCMPLX(1D0,-EPS) CTB = 4*AMB**2/AMA**2*DCMPLX(1D0,-EPS) CLT = 4*AMT**2/AMZ**2*DCMPLX(1D0,-EPS) CLB = 4*AMB**2/AMZ**2*DCMPLX(1D0,-EPS) CAT = FT*(- CI2(CTT,CLT)) CAB = FB*(- CI2(CTB,CLB)) XFAC = CDABS(CAT+CAB)**2 ACOUP = DSQRT(2D0)*GF*AMZ**2*SS*CS/PI**2 HZGA = GF/(4.D0*PI*DSQRT(2.D0))*AMA**3*(ALPH/PI)*ACOUP/16.D0 . *XFAC*(1-AMZ**2/AMA**2)**3 ENDIF C A ---> H Z* ---> HFF IF(IONSH.EQ.0)THEN DLD=3D0 DLU=5D0 XM1 = AML+AMZ-DLD XM2 = AML+AMZ+DLU IF (AMA.LE.AML) THEN HAZ=0 ELSEIF (AMA.LE.XM1) THEN IF (AMA.LE.DABS(AMZ-AML)) THEN HAZ=0 ELSE HAZ=9.D0*GF**2/8.D0/PI**3*AMZ**4*AMA*GZAL**2* . (7.D0/12.D0-10.D0/9.D0*SS+40.D0/27.D0*SS**2) . *HVH((AML/AMA)**2,(AMZ/AMA)**2) ENDIF ELSEIF (AMA.LE.XM2) THEN XX(1) = XM1-1D0 XX(2) = XM1 XX(3) = XM2 XX(4) = XM2+1D0 YY(1)=9.D0*GF**2/8.D0/PI**3*AMZ**4*XX(1)* . (7.D0/12.D0-10.D0/9.D0*SS+40.D0/27.D0*SS**2) . *HVH((AML/XX(1))**2,(AMZ/XX(1))**2) YY(2)=9.D0*GF**2/8.D0/PI**3*AMZ**4*XX(2)* . (7.D0/12.D0-10.D0/9.D0*SS+40.D0/27.D0*SS**2) . *HVH((AML/XX(2))**2,(AMZ/XX(2))**2) CAZ=LAMB(AML**2/XX(3)**2,AMZ**2/XX(3)**2) . *LAMB(XX(3)**2/AMZ**2,AML**2/AMZ**2)**2 YY(3)=GF/8D0/DSQRT(2D0)/PI*AMZ**4/XX(3)*CAZ CAZ=LAMB(AML**2/XX(4)**2,AMZ**2/XX(4)**2) . *LAMB(XX(4)**2/AMZ**2,AML**2/AMZ**2)**2 YY(4)=GF/8D0/DSQRT(2D0)/PI*AMZ**4/XX(4)*CAZ HAZ = FINT(AMA,XX,YY)*GZAL**2 ELSE CAZ=LAMB(AML**2/AMA**2,AMZ**2/AMA**2) . *LAMB(AMA**2/AMZ**2,AML**2/AMZ**2)**2 HAZ=GF/8D0/DSQRT(2D0)/PI*AMZ**4/AMA*GZAL**2*CAZ ENDIF ELSE IF (AMA.LE.AMZ+AML) THEN HAZ=0 ELSE CAZ=LAMB(AML**2/AMA**2,AMZ**2/AMA**2) . *LAMB(AMA**2/AMZ**2,AML**2/AMZ**2)**2 HAZ=GF/8D0/DSQRT(2D0)/PI*AMZ**4/AMA*GZAL**2*CAZ ENDIF ENDIF C C ========================== SUSY DECAYS C IF(IOFSUSY.EQ.0) THEN C A ----> CHARGINOS DO 731 I=1,2 DO 731 J=1,2 IF (AMA.GT.AMCHAR(I)+AMCHAR(J)) THEN WHACH(I,J)=GF*AMW**2/(2*PI*DSQRT(2.D0))/AMA . *LAMB(AMCHAR(I)**2/AMA**2,AMCHAR(J)**2/AMA**2) . *( (AC3(I,J)**2+AC3(J,I)**2)*(AMA**2-AMCHAR(I) . **2-AMCHAR(J)**2)+4.D0*AC3(I,J)*AC3(J,I)* . AMCHAR(I)*AMCHAR(J) ) ELSE WHACH(I,J)=0.D0 ENDIF 731 CONTINUE WHACHT=WHACH(1,1)+WHACH(1,2)+WHACH(2,1)+WHACH(2,2) C A ----> NEUTRALINOS DO 732 I=1,4 DO 732 J=1,4 IF (AMA.GT.AMNEUT(I)+AMNEUT(J)) THEN WHANE(I,J)=GF*AMW**2/(2*PI*DSQRT(2.D0))/AMA . *AN3(I,J)**2*(AMA**2-(XMNEUT(I)-XMNEUT(J))**2) . *LAMB(AMNEUT(I)**2/AMA**2,AMNEUT(J)**2/AMA**2) ELSE WHANE(I,J)=0.D0 ENDIF 732 CONTINUE WHANET= WHANE(1,1)+WHANE(1,2)+WHANE(1,3)+WHANE(1,4) . +WHANE(2,1)+WHANE(2,2)+WHANE(2,3)+WHANE(2,4) . +WHANE(3,1)+WHANE(3,2)+WHANE(3,3)+WHANE(3,4) . +WHANE(4,1)+WHANE(4,2)+WHANE(4,3)+WHANE(4,4) C A ----> STAU'S C IF(AMA.GT.AMSL(1)+AMSL(2)) THEN WHASL=GF*AMZ**4/DSQRT(2.D0)/PI*GAEE**2* . LAMB(AMSL(1)**2/AMA**2,AMSL(2)**2/AMA**2)/AMA ELSE WHASL=0.D0 ENDIF C C A ----> STOPS C IF(AMA.GT.AMST(1)+AMST(2)) THEN WHAST=3*GF*AMZ**4/DSQRT(2.D0)/PI*GATT**2* . LAMB(AMST(1)**2/AMA**2,AMST(2)**2/AMA**2)/AMA ELSE WHAST=0.D0 ENDIF C C A ----> SBOTTOMS C IF(AMA.GT.AMSB(1)+AMSB(2)) THEN WHASB=3*GF*AMZ**4/DSQRT(2.D0)/PI*GABB**2* . LAMB(AMSB(1)**2/AMA**2,AMSB(2)**2/AMA**2)/AMA ELSE WHASB=0.D0 ENDIF C ELSE WHACHT=0.D0 WHANET=0.D0 WHASL=0.D0 WHAST=0.D0 WHASB=0.D0 C--Change thanks to Elzbieta Richter-Was DO I=1,2 DO J=1,2 WHACH(I,J)=0.D0 ENDDO ENDDO DO I=1,4 DO J=1,4 WHANE(I,J)=0.D0 ENDDO ENDDO ENDIF C C ========== TOTAL WIDTH AND BRANCHING RATIOS WTOT=HLL+HMM+HSS+HCC+HBB+HGG+HGA+HZGA+HAZ+HTT . +WHACHT+WHANET+WHASL+WHAST+WHASB ABRT=HTT/WTOT ABRB=HBB/WTOT ABRL=HLL/WTOT ABRM=HMM/WTOT ABRS=HSS/WTOT ABRC=HCC/WTOT ABRG=HGG/WTOT ABRGA=HGA/WTOT ABRZGA=HZGA/WTOT ABRZ=HAZ/WTOT DO 831 I=1,2 DO 831 J=1,2 HABRSC(I,J)=WHACH(I,J)/WTOT 831 CONTINUE DO 832 I=1,4 DO 832 J=1,4 HABRSN(I,J)=WHANE(I,J)/WTOT 832 CONTINUE HABRCHT=WHACHT/WTOT HABRNET=WHANET/WTOT HABRSL=WHASL/WTOT HABRST=WHAST/WTOT HABRSB=WHASB/WTOT C AWDTH=WTOT BHASTAU = WHASL/WTOT BHASB = WHASB/WTOT BHAST = WHAST/WTOT C ============================================================== ENDIF RETURN END DOUBLE PRECISION FUNCTION BIJ(X,Y) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DOUBLE PRECISION LAMB BIJ = (1-X-Y)/LAMB(X,Y)*( . 4*SP(XI(X,Y)*XI(Y,X)) . - 2*SP(-XI(X,Y)) - 2*SP(-XI(Y,X)) . + 2*DLOG(XI(X,Y)*XI(Y,X))*DLOG(1-XI(X,Y)*XI(Y,X)) . - DLOG(XI(X,Y))*DLOG(1+XI(X,Y)) . - DLOG(XI(Y,X))*DLOG(1+XI(Y,X)) . ) . -4*(DLOG(1-XI(X,Y)*XI(Y,X)) . +XI(X,Y)*XI(Y,X)/(1-XI(X,Y)*XI(Y,X))*DLOG(XI(X,Y)*XI(Y,X))) . + (LAMB(X,Y)+X-Y)/LAMB(X,Y)*(DLOG(1+XI(X,Y)) . - XI(X,Y)/(1+XI(X,Y))*DLOG(XI(X,Y))) . + (LAMB(X,Y)-X+Y)/LAMB(X,Y)*(DLOG(1+XI(Y,X)) . - XI(Y,X)/(1+XI(Y,X))*DLOG(XI(Y,X))) RETURN END DOUBLE PRECISION FUNCTION BETA(X) IMPLICIT DOUBLE PRECISION (A-H,O-Z) BETA=DSQRT(1.D0-4.D0*X) RETURN END DOUBLE PRECISION FUNCTION LAMB(X,Y) IMPLICIT DOUBLE PRECISION (A-H,O-Z) LAMB=DSQRT((1.D0-X-Y)**2-4.D0*X*Y) RETURN END DOUBLE PRECISION FUNCTION XI(X,Y) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DOUBLE PRECISION LAMB XI = 2*X/(1-X-Y+LAMB(X,Y)) RETURN END C ***************************************************************** C ************* SUBROUTINE FOR THE SUSY COUPLINGS ***************** C ***************************************************************** SUBROUTINE SUSYCP(TGBET) IMPLICIT DOUBLE PRECISION (A-H,M,O-Z) DOUBLE PRECISION LA1,LA2,LA3,LA4,LA5,LA6,LA7,LA3T COMPLEX*16 F0 DIMENSION MST(2),GLTT(2,2),GHTT(2,2), . MSB(2),GLBB(2,2),GHBB(2,2) COMMON/FLAG/IHIGGS,NNLO,IPOLE COMMON/MASSES/AMS,AMC,AMB,AMT COMMON/HMASS/AMSM,AMA,AML,AMH,AMCH,AMAR COMMON/CHIMASS/AMCHI COMMON/HSELF/LA1,LA2,LA3,LA4,LA5,LA6,LA7 COMMON/BREAK/AMEL,AMER,AMSQ,AMUR,AMDR,AL,AU,AD,AMU,AM2 COMMON/PARAM/GF,ALPH,AMTAU,AMMUON,AMZ,AMW COMMON/COUP/GAT,GAB,GLT,GLB,GHT,GHB,GZAH,GZAL, . GHHH,GLLL,GHLL,GLHH,GHAA,GLAA,GLVV,GHVV, . GLPM,GHPM,B,A PI=4*DATAN(1D0) V=1.D0/DSQRT(DSQRT(2.D0)*GF) BET=DATAN(TGBET) SB = DSIN(BET) CB = DCOS(BET) AMAR = AMA C ============ HEAVIEST CHARGINO MASS NEEDED FOR SUBH ========== AMCHI2=AM2**2+AMU**2+2.D0*AMW**2+DSQRT((AM2**2-AMU**2)**2 . +4.D0*AMW**4*DCOS(2.D0*BET)**2+4.D0*AMW**2* . (AM2**2+AMU**2+2.D0*AMU*AM2*DSIN(2.D0*BET) ) ) AMCHI=DSQRT(0.5D0*AMCHI2) C =============================================================== C ========== RUNNING MASSES CALL SUBH(AMA,TGBET,AMSQ,AMUR,AMDR,AMT,AU,AD,AMU,AMCHI, . AMLR,AMHR,AMCH,SA,CA,TANBA) LA3T=LA3+LA4+LA5 AMA2=AMAR**2 AML2=AMLR**2 AMH2=AMHR**2 AMP2=AMCH**2 C ========== HIGGS COUPLINGS SBMA = SB*CA-CB*SA CBMA = CB*CA+SB*SA SBPA = SB*CA+CB*SA CBPA = CB*CA-SB*SA S2A = 2*SA*CA C2A = CA**2-SA**2 S2B = 2*SB*CB C2B = CB**2-SB**2 GLZZ = 1/V/2*AML2*SBMA GHZZ = 1/V/2*AMH2*CBMA GLWW = 2*GLZZ GHWW = 2*GHZZ GLAZ = 1/V*(AML2-AMA2)*CBMA GHAZ = -1/V*(AMH2-AMA2)*SBMA GLPW = -1/V*(AMP2-AML2)*CBMA GLMW = GLPW GHPW = 1/V*(AMP2-AMH2)*SBMA GHMW = GHPW GAPW = 1/V*(AMP2-AMA2) GAMW = -GAPW GHHH = V/2*(LA1*CA**3*CB + LA2*SA**3*SB + LA3T*SA*CA*SBPA . + LA6*CA**2*(3*SA*CB+CA*SB) + LA7*SA**2*(3*CA*SB+SA*CB)) GLLL = -V/2*(LA1*SA**3*CB - LA2*CA**3*SB + LA3T*SA*CA*CBPA . - LA6*SA**2*(3*CA*CB-SA*SB) + LA7*CA**2*(3*SA*SB-CA*CB)) GLHH = -3*V/2*(LA1*CA**2*CB*SA - LA2*SA**2*SB*CA . + LA3T*(SA**3*CB-CA**3*SB+2*SBMA/3) . - LA6*CA*(CB*C2A-SA*SBPA) - LA7*SA*(C2A*SB+CA*SBPA)) GHLL = 3*V/2*(LA1*SA**2*CB*CA + LA2*CA**2*SB*SA . + LA3T*(SA**3*SB+CA**3*CB-2*CBMA/3) . - LA6*SA*(CB*C2A+CA*CBPA) + LA7*CA*(C2A*SB+SA*CBPA)) GLAA = -V/2*(LA1*SB**2*CB*SA - LA2*CB**2*SB*CA . - LA3T*(SB**3*CA-CB**3*SA) + 2*LA5*SBMA . - LA6*SB*(CB*SBPA+SA*C2B) - LA7*CB*(C2B*CA-SB*SBPA)) GHAA = V/2*(LA1*SB**2*CB*CA + LA2*CB**2*SB*SA . + LA3T*(SB**3*SA+CB**3*CA) - 2*LA5*CBMA . - LA6*SB*(CB*CBPA+CA*C2B) + LA7*CB*(SB*CBPA+SA*C2B)) GLPM = 2*GLAA + V*(LA5 - LA4)*SBMA GHPM = 2*GHAA + V*(LA5 - LA4)*CBMA GLZZ = 2*GLZZ GHZZ = 2*GHZZ GLLL = 6*GLLL GHHH = 6*GHHH GLHH = 2*GLHH GHLL = 2*GHLL GLAA = 2*GLAA GHAA = 2*GHAA XNORM = AMZ**2/V GLLL = GLLL/XNORM GHLL = GHLL/XNORM GLHH = GLHH/XNORM GHHH = GHHH/XNORM GHAA = GHAA/XNORM GLAA = GLAA/XNORM GLPM = GLPM/XNORM GHPM = GHPM/XNORM GAT=1.D0/TGBET GAB=TGBET GLT=CA/SB GLB=-SA/CB GHT=SA/SB GHB=CA/CB GZAL=-CBMA GZAH=SBMA GLVV=SBMA GHVV=CBMA B=BET A=DATAN(SA/CA) IF(CA.LT.0D0)THEN IF(SA.LT.0D0)THEN A = A-PI ELSE A = A+PI ENDIF ENDIF C =============================================================== C ========== POLE MASSES IF(IPOLE.EQ.1) THEN MT=RUNM(AMT,6) MB=RUNM(AMT,5) SW2=1.D0-AMW**2/AMZ**2 C===== STOP MASSES MSTL2=AMSQ**2+(0.5D0-2.D0/3.D0*SW2)*AMZ**2*DCOS(2.D0*B) MSTR2=AMUR**2+2.D0/3.D0*SW2*AMZ**2*DCOS(2.D0*B) MLRT=AU-AMU/TGBET DELT=(MSTL2-MSTR2)**2+4*MT**2*MLRT**2 MST12=MT**2+0.5D0*(MSTL2+MSTR2-DSQRT(DELT)) MST22=MT**2+0.5D0*(MSTL2+MSTR2+DSQRT(DELT)) IF(MST12.LT.0.D0)GOTO 111 MST(1)=DSQRT(MST12) MST(2)=DSQRT(MST22) IF(MSTL2.EQ.MSTR2) THEN THET = PI/4 ELSE THET=0.5D0*DATAN(2.D0*MT*MLRT / (MSTL2-MSTR2) ) IF(MSTL2.GT.MSTR2) THET = THET + PI/2 ENDIF CST= DCOS(THET) SST= DSIN(THET) C===== SBOTTOM MASSES MSBL2=AMSQ**2+(-0.5D0+1.D0/3.D0*SW2)*AMZ**2*DCOS(2.D0*B) MSBR2=AMDR**2-1.D0/3.D0*SW2*AMZ**2*DCOS(2.D0*B) MLRB=AD-AMU*TGBET DELB=(MSBL2-MSBR2)**2+4*MB**2*MLRB**2 MSB12=MB**2+0.5D0*(MSBL2+MSBR2-DSQRT(DELB)) MSB22=MB**2+0.5D0*(MSBL2+MSBR2+DSQRT(DELB)) IF(MSB12.LT.0.D0)GOTO 111 MSB(1)=DSQRT(MSB12) MSB(2)=DSQRT(MSB22) IF(MSBL2.EQ.MSBR2) THEN THEB = PI/4 ELSE THEB=0.5D0*DATAN(2.D0*MB*MLRB / (MSBL2-MSBR2) ) IF(MSBL2.GT.MSBR2) THEB = THEB + PI/2 ENDIF CSB= DCOS(THEB) SSB= DSIN(THEB) C===== LIGHT HIGGS COUPLINGS GLTT(1,1)=-SBPA*(0.5D0*CST**2-2.D0/3.D0*SW2*DCOS(2*THET) ) . +MT**2/AMZ**2*GLT + MT*SST*CST/AMZ**2*(AU*GLT+AMU*GHT) GLTT(2,2)=-SBPA*(0.5D0*SST**2+2.D0/3.D0*SW2*DCOS(2*THET) ) . +MT**2/AMZ**2*GLT - MT*SST*CST/AMZ**2*(AU*GLT+AMU*GHT) GLTT(1,2)=-2*SBPA*SST*CST*(2.D0/3.D0*SW2-0.25D0) . + MT*DCOS(2*THET)/2.D0/AMZ**2*(AU*GLT+AMU*GHT) GLTT(2,1)=-2*SBPA*SST*CST*(2.D0/3.D0*SW2-0.25D0) . + MT*DCOS(2*THET)/2.D0/AMZ**2*(AU*GLT+AMU*GHT) GLBB(1,1)=-SBPA*(-0.5D0*CSB**2+1.D0/3.D0*SW2*DCOS(2*THEB)) . +MB**2/AMZ**2*GLB + MB*SSB*CSB/AMZ**2*(AD*GLB-AMU*GHB) GLBB(2,2)=-SBPA*(-0.5D0*SSB**2-1.D0/3.D0*SW2*DCOS(2*THEB)) . +MB**2/AMZ**2*GLB - MB*SSB*CSB/AMZ**2*(AD*GLB-AMU*GHB) GLBB(1,2)=-2*SBPA*SSB*CSB*(-1.D0/3.D0*SW2+0.25D0) . + MB*DCOS(2*THEB)/2.D0/AMZ**2*(AD*GLB-AMU*GHB) GLBB(2,1)=-2*SBPA*SSB*CSB*(-1.D0/3.D0*SW2+0.25D0) . + MB*DCOS(2*THEB)/2.D0/AMZ**2*(AD*GLB-AMU*GHB) C===== HEAVY HIGGS COUPLINGS GHTT(1,1)=CBPA*(0.5D0*CST**2-2.D0/3.D0*SW2*DCOS(2*THET)) . +MT**2/AMZ**2*GHT + MT*SST*CST/AMZ**2*(AU*GHT-AMU*GLT) GHTT(2,2)=CBPA*(0.5D0*SST**2+2.D0/3.D0*SW2*DCOS(2*THET)) . +MT**2/AMZ**2*GHT - MT*SST*CST/AMZ**2*(AU*GHT-AMU*GLT) GHTT(1,2)=2*CBPA*SST*CST*(2.D0/3.D0*SW2-0.25D0) . +MT*DCOS(2*THET)/2.D0/AMZ**2*(AU*GHT-AMU*GLT) GHTT(2,1)=2*CBPA*SST*CST*(2.D0/3.D0*SW2-0.25D0) . + MT*DCOS(2*THET)/2.D0/AMZ**2*(AU*GHT-AMU*GLT) GHBB(1,1)=CBPA*(-0.5D0*CSB**2+1.D0/3.D0*SW2*DCOS(2*THEB)) . +MB**2/AMZ**2*GHB + MB*SSB*CSB/AMZ**2*(AD*GHB+AMU*GLB) GHBB(2,2)=CBPA*(-0.5D0*SSB**2-1.D0/3.D0*SW2*DCOS(2*THEB)) . + MB**2/AMZ**2*GHB - MB*SSB*CSB/AMZ**2*(AD*GHB+AMU*GLB) GHBB(1,2)=2*CBPA*SSB*CSB*(-1.D0/3.D0*SW2+0.25D0) . + MB*DCOS(2*THEB)/2.D0/AMZ**2*(AD*GHB+AMU*GLB) GHBB(2,1)=2*CBPA*SSB*CSB*(-1.D0/3.D0*SW2+0.25D0) . + MB*DCOS(2*THEB)/2.D0/AMZ**2*(AD*GHB+AMU*GLB) C===== PSEUDOSCALAR HIGGS COUPLINGS GATT=-MT/2.D0/AMZ**2*(AMU+AU*GAT) GABB=-MB/2.D0/AMZ**2*(AMU+AD*GAB) C======= LOOP CORRECTIONS XDLT=GF/(2.D0*DSQRT(2.D0)*PI**2)*GLT**2*(-2.D0*MT**2+0.5D0*AML2) . *DREAL(F0(MT,MT,AML2)) . *3*MT**2 XDLB=GF/(2.D0*DSQRT(2.D0)*PI**2)*GLB**2*(-2.D0*MB**2+0.5D0*AML2) . *DREAL(F0(MB,MB,AML2)) . *3*MB**2 C--BUG IN CARENA ET AL. FIXED . +GF/(2.D0*DSQRT(2.D0)*PI**2)*GLB**2*(0.5D0*AML2) . *DLOG(MB**2/MT**2) . *3*MB**2 XDHT=GF/(2.D0*DSQRT(2.D0)*PI**2)*GHT**2*(-2.D0*MT**2+0.5D0*AMH2) . *DREAL(F0(MT,MT,AMH2)) . *3*MT**2 XDHB=GF/(2.D0*DSQRT(2.D0)*PI**2)*GHB**2*(-2.D0*MB**2+0.5D0*AMH2) . *DREAL(F0(MB,MB,AMH2)) . *3*MB**2 C--BUG IN CARENA ET AL. FIXED . +GF/(2.D0*DSQRT(2.D0)*PI**2)*GHB**2*(0.5D0*AMH2) . *DLOG(MB**2/MT**2) . *3*MB**2 XDAT=GF/(2.D0*DSQRT(2.D0)*PI**2)*GAT**2*(-0.5D0*AMA2) . *DREAL(F0(MT,MT,AMA2)) . *3*MT**2 XDAB=GF/(2.D0*DSQRT(2.D0)*PI**2)*GAB**2*(-0.5D0*AMA2) . *DREAL(F0(MB,MB,AMA2)) . *3*MB**2 C--BUG IN CARENA ET AL. FIXED . +GF/(2.D0*DSQRT(2.D0)*PI**2)*GAB**2*(-0.5D0*AMA2) . *DLOG(MB**2/MT**2) . *3*MB**2 XDLST=0.D0 XDLSB=0.D0 XDHST=0.D0 XDHSB=0.D0 DO 311 I=1,2 DO 311 J=1,2 XDLST=XDLST+GF/(2.D0*DSQRT(2.D0)*PI**2)*GLTT(I,J)**2* . DREAL(F0(MST(I),MST(J),AML2)) . *3*AMZ**4 XDLSB=XDLSB+GF/(2.D0*DSQRT(2.D0)*PI**2)*GLBB(I,J)**2* . DREAL(F0(MSB(I),MSB(J),AML2)) . *3*AMZ**4 XDHST=XDHST+GF/(2.D0*DSQRT(2.D0)*PI**2)*GHTT(I,J)**2* . DREAL(F0(MST(I),MST(J),AMH2)) . *3*AMZ**4 XDHSB=XDHSB+GF/(2.D0*DSQRT(2.D0)*PI**2)*GHBB(I,J)**2* . DREAL(F0(MSB(I),MSB(J),AMH2)) . *3*AMZ**4 311 CONTINUE XDAST=GF/(1.D0*DSQRT(2.D0)*PI**2)*GATT**2* . DREAL(F0(MST(1),MST(2),AMA2)) . *3*AMZ**4 XDASB=GF/(1.D0*DSQRT(2.D0)*PI**2)*GABB**2* . DREAL(F0(MSB(1),MSB(2),AMA2)) . *3*AMZ**4 AML=DSQRT(AML2+XDLT+XDLB+XDLST+XDLSB) AMH=DSQRT(AMH2+XDHT+XDHB+XDHST+XDHSB) AMA=DSQRT(AMA2+XDAT+XDAB+XDAST+XDASB) ELSE AML=AMLR AMH=AMHR AMA=AMAR ENDIF RETURN 111 STOP END C ===================== THE FUNCTION F0 =============== COMPLEX FUNCTION F0*16(M1,M2,QSQ) IMPLICIT REAL*8 (A-H,M,O-Z) COMPLEX*16 CD,CR,CQ2,IEPS,CBET,CXX M1SQ = M1*M1 M2SQ = M2*M2 AQSQ = DABS(QSQ) IEPS = DCMPLX(1.D0,1.D-12) CQ2 = QSQ*IEPS CD = (M1SQ-M2SQ)/CQ2 CR = CDSQRT((1+CD)**2 - 4*M1SQ/CQ2) IF(QSQ.EQ.0.D0) THEN F0 = 0.D0 ELSE IF(M1.EQ.M2) THEN F0 = -2.D0 + CR*CDLOG(-(1+CR)/(1-CR)) ELSE CBET = CDSQRT(1-4*M1*M2/(CQ2 - (M1-M2)**2)) CXX = (CBET-1)/(CBET+1) F0 = -1 + ((QSQ+M2SQ-M1SQ)/2/QSQ - M2SQ/(M2SQ-M1SQ)) . *DLOG(M2SQ/M1SQ) . - (QSQ-(M1-M2)**2)/QSQ*CBET*CDLOG(CXX) ENDIF ENDIF RETURN END C ************************************************************ C SUBROUTINE FOR HSM ---> V*V* ---> 4F C ************************************************************ SUBROUTINE HTOVV(AMH,AMV,GAMV,HTVV) IMPLICIT DOUBLE PRECISION (A-H,O-Z) COMMON/VVOFF/AMH1,AMV1,GAMV1 COMMON/PREC/IP EXTERNAL FTOVV1 IP=20 AMH1=AMH AMV1=AMV GAMV1=GAMV DLT=1D0/IP SUM=0D0 DO 1 I=1,IP UU=DLT*I DD=UU-DLT CALL QGAUS1(FTOVV1,DD,UU,RES) SUM=SUM+RES 1 CONTINUE HTVV=SUM RETURN END DOUBLE PRECISION FUNCTION FTOVV1(XX) IMPLICIT DOUBLE PRECISION (A-H,O-Z) COMMON/FIRST/X1 COMMON/PREC/IP EXTERNAL FTOVV2 X1=XX DLT=1D0/IP SUM=0D0 DO 1 I=1,IP UU=DLT*I DD=UU-DLT CALL QGAUS2(FTOVV2,DD,UU,RES) SUM=SUM+RES 1 CONTINUE FTOVV1=SUM RETURN END DOUBLE PRECISION FUNCTION FTOVV2(XX) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION YY(2) COMMON/FIRST/X1 YY(1)=X1 YY(2)=XX FTOVV2=FTOVV(YY) RETURN END DOUBLE PRECISION FUNCTION FTOVV(XX) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DOUBLE PRECISION LAMB DIMENSION XX(2) COMMON/VVOFF/AMH,AMW,GAMW LAMB(X,Y)=DSQRT((1.D0-X-Y)**2-4.D0*X*Y) PI=4D0*DATAN(1D0) ICASE = 1 IF(ICASE.EQ.0)THEN YY = AMH**2 Y1 = DATAN((YY-AMW**2)/AMW/GAMW) Y2 = -DATAN((AMW**2)/AMW/GAMW) DJAC = Y1-Y2 T1 = TAN(Y1*XX(1)+Y2*(1.D0-XX(1))) SP = AMW**2 + AMW*GAMW*T1 YY = (AMH-DSQRT(SP))**2 Y1 = DATAN((YY-AMW**2)/AMW/GAMW) Y2 = -DATAN((AMW**2)/AMW/GAMW) DJAC = DJAC*(Y1-Y2) T2 = TAN(Y1*XX(2)+Y2*(1.D0-XX(2))) SM = AMW**2 + AMW*GAMW*T2 AM2=AMH**2 GAM = AM2*LAMB(SP/AM2,SM/AM2)*(1+LAMB(SP/AM2,SM/AM2)**2*AMH**4 . /SP/SM/12) PRO1 = SP/AMW**2 PRO2 = SM/AMW**2 FTOVV = PRO1*PRO2*GAM*DJAC/PI**2 ELSE SP = AMH**2*XX(1) SM = (AMH-DSQRT(SP))**2*XX(2) DJAC = AMH**2*(AMH-DSQRT(SP))**2/PI**2 AM2=AMH**2 GAM = AM2*LAMB(SP/AM2,SM/AM2)*(1+LAMB(SP/AM2,SM/AM2)**2*AMH**4 . /SP/SM/12) PRO1 = SP*GAMW/AMW/((SP-AMW**2)**2+AMW**2*GAMW**2) PRO2 = SM*GAMW/AMW/((SM-AMW**2)**2+AMW**2*GAMW**2) FTOVV = PRO1*PRO2*GAM*DJAC ENDIF RETURN END C ************************************************************ C SUBROUTINE FOR HSM ---> TT* ---> TBW C ************************************************************ SUBROUTINE HTOTTS(AMH,AMT,AMB,AMW,HTTS) IMPLICIT REAL*8(A-Z) INTEGER IP,K COMMON/PREC1/IP EXTERNAL FUNSTT1 COMMON/IKSY0/X1,X2,M1,M2,M3,ECM,S COMMON/TOP0/AMH0,AMT0,AMB0,AMW0 AMH0=AMH AMT0=AMT AMB0=AMB AMW0=AMW IP=5 M1=AMB M2=AMT M3=AMW C FIRST INTEGRATE OVER X2, i.e. (1+3) SYSTEM C CHECK WHETHER ENOUGH PHASE SPACE MASTOT=M1+M2+M3 IF(MASTOT.GE.AMH) GOTO 12 ECM=AMH S=ECM**2 U1=(ECM-M2)**2 D1=(M1+M3)**2 U=(S-D1+M2**2)/s D=(S-U1+M2**2)/s DEL=(U-D)/IP U=D+DEL XSEC=0.D0 DO K=1,IP CALL QGAUS1(FUNSTT1,D,U,SS) D=U U=D+DEL XSEC=XSEC+SS ENDDO HTTS=XSEC 12 CONTI