+++ /dev/null
-#include "isajet/pilot.h"
- SUBROUTINE SSMHN(MHLNEG)
-C-----------------------------------------------------------------------
-C
-C Calculate HL, HH masses and ALFAH
-C (scalar Higgs mixing angle) using radiative
-C corrections calculated by M. Bisset
-C and save results in /SSPAR/.
-C
-C Both top and bottom couplings are now
-C included. Non-degenerate mixed squark
-C masses and A-terms are also included.
-C The D-terms from the squark mass matrix
-C (terms prop. to g**2 * Yukawa coupling)
-C are included as an option:
-C INRAD = 1 ==> D-TERMS ON
-C INRAD = 2 ==> D-TERMS OFF .
-C
-C 10/18/93 D-terms are now turned on.
-C INRAD = 1
-C
-C There is an arbitrary mass scale that must
-C chosen to avoid dimensionful logarithms.
-C The choice does not matter if D-terms are
-C not included, but it does matter if D-terms
-C are included.
-C
-C Arbitrary mass scale updated to
-C QQQ = HIGFRZ = SQRT(AMTLSS*AMTRSS)
-C with running masses to include dominant 2-loop
-C effects. 12/10/96 H. Baer
-C
-C It is assumed that the A-terms are real.
-C Complex A-terms are allowed
-C (unless RTT=0 or RBB=0 --see below) in
-C this subroutine, but the imaginary parts
-C are now set to zero.
-C
-C-----------------------------------------------------------------------
-#if defined(CERNLIB_IMPNONE)
- IMPLICIT NONE
-#endif
-#include "isajet/sslun.inc"
-#include "isajet/sssm.inc"
-#include "isajet/sspar.inc"
-C
- REAL PI,PI2,SR2,G2,GP2,GGP,GG1,GG2
- REAL TANB,COTB,COSB,SINB,BE
- REAL SINB2,COSB2,COS2B,SIN2B
- REAL V2,VP2,V,VP,VVP,VPVM,VVPP
- REAL MT2,MB2,FT2,FB2,MW2,ZAP,QQQ2
- REAL EP,EP2,RR,MHP2
- REAL ATI,ABI,ATR,ABR,AT2,AB2
- REAL TLRM,BLRM,TLRP,BLRP
- REAL MST1SQ,MST2SQ,MSB1SQ,MSB2SQ
-
- REAL RTT,TTT1,TEMPT,TM1BT
- REAL TEMPS,T1RD,T2RD,T1RPD,T2RPD
- REAL CT1,A1,A2,T1RR,T2RR
- REAL CT5,A5,A6,T1RPRP,T2RPRP
- REAL A9,T1RRP,T2RRP
- REAL TEMPSQ,DT1,DT2,VRRT,VRPRPT,VRRPT
- REAL ALPHAT,LAT
-C
- REAL RBB,BBB1,TEMPB,TM1BB
- REAL B1RD,B2RD,B1RPD,B2RPD
- REAL CB3,A3,A4,B1RR,B2RR
- REAL CB7,A7,A8,B1RPRP,B2RPRP
- REAL A10,B1RRP,B2RRP
- REAL DB1,DB2,VRRB,VRPRPB,VRRPB
- REAL ALPHAB,LAB
-C
- REAL DVRR,DVRPRP,DVRRP,TEMPH
- REAL MHL2,MHH2,TRACEM,TPAL,TANAH
- REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS,HIGFRZ
- DOUBLE PRECISION SSMQCD
- INTEGER INRAD,MHLNEG
-C
- MHLNEG=0
- PI=4.*ATAN(1.)
- PI2 = PI**2
- SR2=SQRT(2.)
- G2=4.*PI*ALFAEM/SN2THW
- GP2=G2*SN2THW/(1.-SN2THW)
- HIGFRZ=SQRT(AMTLSS*AMTRSS)
- ASMB=SUALFS(AMBT**2,.36,AMTP,3)
- MBMB=AMBT*(1.-4*ASMB/3./PI)
- MBQ=SSMQCD(DBLE(MBMB),DBLE(HIGFRZ))
- ASMT=SUALFS(AMTP**2,.36,AMTP,3)
- MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
- $(ASMT/PI)**2)
- MTQ=SSMQCD(DBLE(MTMT),DBLE(HIGFRZ))
- MT2=MTQ**2
- MB2=MBQ**2
- MW2=AMW**2
- EP=TWOM1
- EP2=EP**2
- RR=RV2V1
- MHP2=AMHA**2
- TANB=1.0/RR
- COTB=RR
- BE=ATAN(1./RV2V1)
- SINB=SIN(BE)
- COSB=COS(BE)
- SINB2=SINB**2
- COSB2=COSB**2
- SIN2B=SIN(2.*BE)
- COS2B=COS(2.*BE)
- V2=2.0*MW2*SINB2/G2
- VP2=2.0*MW2*COSB2/G2
- V=SQRT(V2)
- VP=SQRT(VP2)
- VVP=SQRT(V2*VP2)
- VPVM=VP2-V2
- GGP=G2+GP2
- GG1=G2-5.0*GP2/3.0
- GG2=G2-GP2/3.0
- VVPP=2.0*AMZ**2/GGP
- FT2=MT2/V2
- FB2=MB2/VP2
-C
- TLRM=AMTLSS**2-AMTRSS**2
- BLRM=AMBLSS**2-AMBRSS**2
- TLRP=AMTLSS**2+AMTRSS**2
- BLRP=AMBLSS**2+AMBRSS**2
-C
-C Higgs mass matrix
-C
-C (AAT and AAB are also assumed to be real)
-C
- ATR=AAT
- ABR=AAB
- ATI=0.0
- ABI=0.0
- AT2=ATR**2+ATI**2
- AB2=ABR**2+ABI**2
-C
- MST1SQ=AMT1SS**2
- MST2SQ=AMT2SS**2
- MSB1SQ=AMB1SS**2
- MSB2SQ=AMB2SS**2
- INRAD=1
- QQQ2=HIGFRZ**2
-C
- ZAP = 1.0
-C
-C STOP TERMS
-C
- RTT=(TLRM+VPVM*ZAP*GG1/4.0)**2
- $ +4.0*MT2*(EP*COTB+ATR)**2+4.0*MT2*ATI**2
- RTT=SQRT(RTT)
-C
-C calculate 2M1*B term
-C
- TTT1=0.5*TLRP+MT2+VPVM*ZAP*GGP/8.0
- IF(RTT.NE.0.0) THEN
- TEMPT=4.0*EP*FT2*VVP*ATI**2/(RTT**2)
- TM1BT=-2.0*FT2*(TEMPT+ATR)*TTT1
- $ *LOG(MST2SQ/MST1SQ)/RTT
- TM1BT=TM1BT-FT2*ATR
- $ *LOG(MST1SQ*MST2SQ/QQQ2/QQQ2)
- TM1BT=TM1BT+FT2*(2.0*TEMPT-ATR)
- TM1BT=3.0*EP*TM1BT/32.0/PI2
-C
-C calculate first derivatives w.r.t H_R
-C divided by sqrt(2) * v
-C
- TEMPS=-ZAP*GG1*(TLRM+ZAP*GG1*VPVM/4.0)/2.0
- TEMPS=TEMPS+4.0*FT2*(AT2+EP*COTB*ATR)
- TEMPS=TEMPS/RTT/4.0
- T1RD=FT2-ZAP*GGP/8.0-TEMPS
- T2RD=FT2-ZAP*GGP/8.0+TEMPS
-C
-C calculate first derivatives w.r.t H_R'
-C divided by sqrt(2) * v'
-C
- TEMPS=ZAP*GG1*(TLRM+ZAP*GG1*VPVM/4.0)/2.0
- TEMPS=TEMPS+4.0*FT2*EP*(EP+TANB*ATR)
- TEMPS=TEMPS/RTT/4.0
- T1RPD=ZAP*GGP/8.0-TEMPS
- T2RPD=ZAP*GGP/8.0+TEMPS
-C
-C calculate second derivatives w.r.t. H_R
-C
- CT1=-V*ZAP*GG1*(TLRM+ZAP*GG1*VPVM/4.0)/SR2
- CT1=CT1+4.0*SR2*FT2*V*(EP*COTB*ATR+AT2)
- A1=-CT1**2/(RTT**3)/8.0
- A2=-ZAP*GG1*(TLRM+ZAP*GG1*VPVM/4.0)/2.0
- A2=A2+V2*ZAP*GG1**2/4.0+4.0*FT2*AT2
- A2=A2/RTT/4.0
- T1RR=FT2-ZAP*GGP/8.0-A1-A2
- T2RR=FT2-ZAP*GGP/8.0+A1+A2
-C
-C calculate second derivatives w.r.t. H_R'
-C
- CT5=VP*ZAP*GG1*(TLRM+ZAP*GG1*VPVM/4.0)/SR2
- CT5=CT5+4.0*SR2*FT2*VP*EP*(EP+TANB*ATR)
- A5=-CT5**2/(RTT**3)/8.0
- A6=ZAP*GG1*(TLRM+ZAP*GG1*VPVM/4.0)/2.0
- A6=A6+VP2*ZAP*GG1**2/4.0+4.0*FT2*EP2
- A6=A6/RTT/4.0
- T1RPRP=ZAP*GGP/8.0-A5-A6
- T2RPRP=ZAP*GGP/8.0+A5+A6
-C
-C calculate second derivatives w.r.t. H_R and H_R'
-C
- A9=-VVP*ZAP*(GG1**2)/4.0+4.0*FT2*EP*ATR
- A9=A9/RTT/4.0
- T1RRP=CT1*CT5/(RTT**3)/8.0-A9
- T2RRP=-CT1*CT5/(RTT**3)/8.0+A9
-C
-C calculate D^2 V / D^2 H_R
-C
- TEMPSQ=MST1SQ*(T1RR-T1RD)
- DT1=2.0*(2.0*V2*T1RD**2+TEMPSQ)*LOG(MST1SQ/QQQ2)
- DT1=DT1+6.0*V2*T1RD**2+TEMPSQ
- TEMPSQ=MST2SQ*(T2RR-T2RD)
- DT2=2.0*(2.0*V2*T2RD**2+TEMPSQ)*LOG(MST2SQ/QQQ2)
- DT2=DT2+6.0*V2*T2RD**2+TEMPSQ
- VRRT=DT1+DT2-8.0*FT2*MT2*LOG(MT2/QQQ2)-12.0*FT2*MT2
- VRRT=-TM1BT*COTB+3.0*VRRT/32.0/PI2
-C
-C calculate D^2 V / D^2 H'_R
-C
- TEMPSQ=MST1SQ*(T1RPRP-T1RPD)
- DT1=2.0*(2.0*VP2*T1RPD**2+TEMPSQ)*LOG(MST1SQ/QQQ2)
- DT1=DT1+6.0*VP2*T1RPD**2+TEMPSQ
- TEMPSQ=MST2SQ*(T2RPRP-T2RPD)
- DT2=2.0*(2.0*VP2*T2RPD**2+TEMPSQ)*LOG(MST2SQ/QQQ2)
- DT2=DT2+6.0*VP2*T2RPD**2+TEMPSQ
- VRPRPT=-TM1BT*TANB+3.0*(DT1+DT2)/32.0/PI2
-C
-C calculate D^2 V / D^H_R D^H_R'
-C
- DT1=2.0*VVP*T1RD*T1RPD+MST1SQ*T1RRP
- DT1=2.0*DT1*LOG(MST1SQ/QQQ2)
- DT1=DT1+6.0*VVP*T1RD*T1RPD+MST1SQ*T1RRP
- DT2=2.0*VVP*T2RD*T2RPD+MST2SQ*T2RRP
- DT2=2.0*DT2*LOG(MST2SQ/QQQ2)
- DT2=DT2+6.0*VVP*T2RD*T2RPD+MST2SQ*T2RRP
- VRRPT=TM1BT+3.0*(DT1+DT2)/32.0/PI2
-C
- ELSE IF(RTT.EQ.0.0) THEN
-C
- ALPHAT=TLRP/2.0+MT2+ZAP*GGP*VPVM/8.0
- LAT=2.0*LOG(ALPHAT/QQQ2)+3.0
-C
-C calculate D^2 V / D^2 H_R
-C
- VRRT=V2*(GGP**2+GG1**2)/16.0-MT2*GGP
- VRRT=ZAP*VRRT*LAT+8.0*FT2*MT2*LOG(ALPHAT/MT2)
- VRRT=3.0*VRRT/32.0/PI2
-C
-C calculate D^2 V / D^2 H_R'
-C
- VRPRPT=ZAP*VP2*(GGP**2+GG1**2)/16.0
- VRPRPT=3.0*(VRPRPT*LAT)/32.0/PI2
-C
-C calculate D^2 V / D^H_R D^H_R'
-C
- VRRPT=FT2*GGP-(GGP**2+GG1**2)/8.0
- VRRPT=ZAP*VVP*VRRPT*LAT/2.0
- VRRPT=3.0*VRRPT/32.0/PI2
-C
-C
- ENDIF
-C
-C SBOTTOM TERMS
-C
- RBB=(BLRM-VPVM*ZAP*GG2/4.0)**2
- $ +4.0*MB2*(EP*TANB+ABR)**2+4.0*MB2*ABI**2
- RBB=SQRT(RBB)
-C IF(RBB.EQ.0.0.AND.ABI.NE.0.0) THEN
-C WRITE(6,*) 'RBB=0, ABI NOT 0'
-C WRITE(6,*) 'ERROR: THIS CASE NOT COVERED YET'
-C GO TO 1000
-C ENDIF
-C
- IF(RBB.NE.0.0) THEN
-C
-C calculate 2M1*B term
-C
- BBB1=0.5*BLRP+MB2-VPVM*ZAP*GGP/8.0
- TEMPB=4.0*EP*FB2*VVP*ABI**2/(RBB**2)
- TM1BB=-2.0*FB2*(TEMPB+ABR)*BBB1
- $ *LOG(MSB2SQ/MSB1SQ)/RBB
- TM1BB=TM1BB-FB2*ABR
- $ *LOG(MSB1SQ*MSB2SQ/QQQ2/QQQ2)
- TM1BB=TM1BB+FB2*(2.0*TEMPB-ABR)
- TM1BB=3.0*EP*TM1BB/32.0/PI2
-C
-C calculate first derivatives w.r.t H_R
-C divided by sqrt(2) * v
-C
- TEMPS=ZAP*GG2*(BLRM-ZAP*GG2*VPVM/4.0)/2.0
- TEMPS=TEMPS+4.0*FB2*EP*(EP+COTB*ABR)
- TEMPS=TEMPS/RBB/4.0
- B1RD=ZAP*GGP/8.0-TEMPS
- B2RD=ZAP*GGP/8.0+TEMPS
-
-C calculate first derivatives w.r.t H_R'
-C divided by sqrt(2) * v'
-C
- TEMPS=-ZAP*GG2*(BLRM-ZAP*GG2*VPVM/4.0)/2.0
- TEMPS=TEMPS+4.0*FB2*(AB2+EP*TANB*ABR)
- TEMPS=TEMPS/RBB/4.0
- B1RPD=FB2-ZAP*GGP/8.0-TEMPS
- B2RPD=FB2-ZAP*GGP/8.0+TEMPS
-C
-C calculate second derivatives w.r.t. H_R
-C
- CB3=V*ZAP*GG2*(BLRM-ZAP*GG2*VPVM/4.0)/SR2
- CB3=CB3+4.0*SR2*FB2*V*EP*(EP+COTB*ABR)
- A3=-CB3**2/(RBB**3)/8.0
- A4=ZAP*GG2*(BLRM-ZAP*GG2*VPVM/4.0)/2.0
- A4=A4+V2*ZAP*GG2**2/4.0+4.0*FB2*EP2
- A4=A4/RBB/4.0
- B1RR=ZAP*GGP/8.0-A3-A4
- B2RR=ZAP*GGP/8.0+A3+A4
-C
-C calculate second derivatives w.r.t. H_R'
-C
- CB7=-VP*ZAP*GG2*(BLRM-ZAP*GG2*VPVM/4.0)/SR2
- CB7=CB7+4.0*SR2*FB2*VP*(AB2+EP*TANB*ABR)
- A7=-CB7**2/(RBB**3)/8.0
- A8=-ZAP*GG2*(BLRM-ZAP*GG2*VPVM/4.0)/2.0
- A8=A8+VP2*ZAP*GG2**2/4.0+4.0*FB2*AB2
- A8=A8/RBB/4.0
- B1RPRP=FB2-ZAP*GGP/8.0-A7-A8
- B2RPRP=FB2-ZAP*GGP/8.0+A7+A8
-C
-C calculate second derivatives w.r.t. H_R and H_R'
-C
- A10=-VVP*ZAP*(GG2**2)/4.0+4.0*FB2*EP*ABR
- A10=A10/RBB/4.0
- B1RRP=CB3*CB7/(RBB**3)/8.0-A10
- B2RRP=-CB3*CB7/(RBB**3)/8.0+A10
-C
-C calculate D^2 V / D^2 H_R
-C
- TEMPSQ=MSB1SQ*(B1RR-B1RD)
- DB1=2.0*(2.0*V2*B1RD**2+TEMPSQ)*LOG(MSB1SQ/QQQ2)
- DB1=DB1+6.0*V2*B1RD**2+TEMPSQ
- TEMPSQ=MSB2SQ*(B2RR-B2RD)
- DB2=2.0*(2.0*V2*B2RD**2+TEMPSQ)*LOG(MSB2SQ/QQQ2)
- DB2=DB2+6.0*V2*B2RD**2+TEMPSQ
- VRRB=-TM1BB*COTB+3.0*(DB1+DB2)/32.0/PI2
-C
-C calculate D^2 V / D^2 H'_R
-C
- TEMPSQ=MSB1SQ*(B1RPRP-B1RPD)
- DB1=2.0*(2.0*VP2*B1RPD**2+TEMPSQ)*LOG(MSB1SQ/QQQ2)
- DB1=DB1+6.0*VP2*B1RPD**2+TEMPSQ
- TEMPSQ=MSB2SQ*(B2RPRP-B2RPD)
- DB2=2.0*(2.0*VP2*B2RPD**2+TEMPSQ)*LOG(MSB2SQ/QQQ2)
- DB2=DB2+6.0*VP2*B2RPD**2+TEMPSQ
- VRPRPB=DB1+DB2
- VRPRPB=DB1+DB2-8.0*FB2*MB2*LOG(MB2/QQQ2)-12.0*FB2*MB2
- VRPRPB=-TM1BB*TANB+3.0*VRPRPB/32.0/PI2
-C
-C calculate D^2 V / D H_R D H'_R
-C
- DB1=2.0*VVP*B1RD*B1RPD+MSB1SQ*B1RRP
- DB1=2.0*DB1*LOG(MSB1SQ/QQQ2)
- DB1=DB1+6.0*VVP*B1RD*B1RPD+MSB1SQ*B1RRP
- DB2=2.0*VVP*B2RD*B2RPD+MSB2SQ*B2RRP
- DB2=2.0*DB2*LOG(MSB2SQ/QQQ2)
- DB2=DB2+6.0*VVP*B2RD*B2RPD+MSB2SQ*B2RRP
- VRRPB=TM1BB+3.0*(DB1+DB2)/32.0/PI2
-
- ELSE IF(RBB.EQ.0.0) THEN
-C
- ALPHAB=BLRP/2.0+MB2-ZAP*GGP*VPVM/8.0
- LAB=2.0*LOG(ALPHAB/QQQ2)+3.0
-C
-C calculate D^2 V / D^2 H_R
-C
- VRRB=ZAP*V2*(GGP**2 + GG2**2)/16.0
- VRRB=3.0*(VRRB*LAB)/32.0/PI2
-C
-C calculate D^2 V / D^2 H_R'
-C
- VRPRPB=VP2*(GGP**2+GG2**2)/16.0-MB2*GGP
- VRPRPB=ZAP*VRPRPB*LAB+8.0*FB2*MB2*LOG(ALPHAB/MB2)
- VRPRPB=3.0*VRPRPB/32.0/PI2
-C
-C calculate D^2 V / D^H_R D^H_R'
-C
- VRRPB=FB2*GGP-(GGP**2+GG2**2)/8.0
- VRRPB=ZAP*VVP*VRRPB*LAB/2.0
- VRRPB=3.0*VRRPB/32.0/PI2
-C
- ENDIF
-C
- DVRR=VRRT+VRRB+VP2*MHP2/VVPP + V2*GGP/2.0
- DVRPRP=VRPRPT+VRPRPB+V2*MHP2/VVPP + VP2*GGP/2.0
- DVRRP=VRRPT+VRRPB-VVP*MHP2/VVPP - VVP*GGP/2.0
-C TEMPH is always non-negative:
- TEMPH=(DVRR-DVRPRP)**2+4*DVRRP**2
- TEMPH=0.5*SQRT(TEMPH)
- MHL2=0.5*(DVRR+DVRPRP)-TEMPH
- MHH2=0.5*(DVRR+DVRPRP)+TEMPH
- IF(MHL2.LT.0.0) THEN
- MHLNEG=1
-C WRITE(LOUT,*) 'SSMHN: ERROR: MHL**2 < 0.0 FOR PARAMETERS:'
-C WRITE(LOUT,*) 'MHP =', AMHA, 'TANB =', 1.0/RR
-C WRITE(LOUT,*) 'MSTL=', AMTLSS, 'MSBL=', AMBLSS
-C WRITE(LOUT,*) 'MSTR=', AMTRSS, 'MSBR=', AMBRSS
-C WRITE(LOUT,*) 'AT=', AAT, 'AB=', AAB
-C WRITE(LOUT,*) 'MU=-2M1=', -EP
-C WRITE(LOUT,*) 'MT=', AMTP, 'MB=', AMBT
-C WRITE(LOUT,*) 'D-TERMS? 1=YES 2=NO :', INRAD
-C WRITE(LOUT,*) 'MASS SCALE (QQQ)=', SQRT(QQQ2)
- AMHH=SQRT(MHH2)
- AMHL=SQRT(ABS(MHL2))
- GO TO 1000
- ENDIF
- AMHL=SQRT(MHL2)
- AMHH=SQRT(MHH2)
-
-C
-C Now calculate mixing angle ALFAH
-C
- TRACEM=DVRR-DVRPRP
- TPAL=TRACEM**2 + 4.0*DVRRP**2
- TANAH=TRACEM+SQRT(TPAL)
- IF(DVRRP.EQ.0.0) THEN
- WRITE(LOUT,*) 'SSMHN: OFF-DIAGONAL TERM OF SCALAR HIGGS',
- $ ' MASS MATRIX IS ZERO '
- IF(TANAH.NE.0.0) THEN
- WRITE(LOUT,*) 'SSMHN: WARNING: TAN(ALFAH) FORMULA',
- $ ' YIELDS INFINITY'
- ELSE IF(TANAH.EQ.0.0) THEN
- WRITE(LOUT,*) 'SSMHN: WARNING: TAN(ALFAH) FORMULA',
- $ ' YIELDS 0/0 '
- ENDIF
- IF(DVRR.GT.DVRPRP) THEN
- WRITE(LOUT,*) 'SSMHN: DVRR > DVRPRP ==> SET ALFAH=PI/2'
- ALFAH = PI/2.0
- ELSE IF (DVRR .LT. DVRPRP) THEN
- WRITE(LOUT,*) 'SSMHN: DVRR < DVRPRP ==> SET ALFAH=0'
- ALFAH = 0.0
- ELSE IF (DVRR .EQ. DVRPRP) THEN
- WRITE(LOUT,*) 'SSMHN: DVRR = DVRPRP ==> ALFAH INDETERMINANT'
- WRITE(LOUT,*) 'SETTING SCALAR MIXING ANGLE ALPHA=PI/4'
- ALFAH=PI/4.0
- ENDIF
- GO TO 1000
- ENDIF
- TANAH = -0.5*TANAH/DVRRP
- ALFAH = ATAN(TANAH)
-C
-1000 RETURN
- END