C-------------------------------------------------------------------- SUBROUTINE ISALHA(IMODEL,IMODIN,MT) C-------------------------------------------------------------------- C C Print SUGRA parameters in 'Les Houches accord 3' (LHA3) format C C. Balazs, Jan. 5 2005, v1.0 C (History: C. Balazs, July 24 2003, v0.1) C C IMODEL = model type for SUGRA C IMODIN = input model type to control formatting C MT = top mass CsB ISAJET common blocks from SUGPRT ... +CDE,SSLUN +CDE,SUGXIN +CDE,SUGMG +CDE,SUGPAS +CDE,SSPAR +CDE,SUGNU CsB Additional ISAJET common blocks ... +CDE,SSSM +CDE,CONST CsB ... explicitly included from v7.71 ! COMMON/SSLUN/LOUT ! INTEGER LOUT ! SAVE /SSLUN/ !C XSUGIN contains the inputs to SUGRA: !C XSUGIN(1) = M_0 XSUGIN(2) = M_(1/2) XSUGIN(3) = A_0 !C XSUGIN(4) = tan(beta) XSUGIN(5) = sgn(mu) XSUGIN(6) = M_t !C XSUGIN(7) = SUG BC scale !C XGMIN(1) = LAM XGMIN(2) = M_MES XGMIN(3) = XN5 !C XGMIN(4) = tan(beta) XGMIN(5) = sgn(mu) XGMIN(6) = M_t !C XGMIN(7) = CGRAV XGMIN(8) =RSL XGMIN(9) = DEL_HD !C XGMIN(10) = DEL_HU XGMIN(11) = DY XGMIN(12) = N5_1 !C XGMIN(13) = N5_2 XGMIN(14) = N5_3 !C XNRIN(1) = M_N3 XNRIN(2) = M_MAJ XNRIN(3) = ANSS !C XNRIN(4) = M_N3SS !C XISAIN contains the MSSMi inputs in natural order. ! COMMON /SUGXIN/ XISAIN(24),XSUGIN(7),XGMIN(14),XNRIN(4), ! $XAMIN(7) ! REAL XISAIN,XSUGIN,XGMIN,XNRIN,XAMIN ! SAVE /SUGXIN/ !C Frozen couplings from RG equations: !C GSS( 1) = g_1 GSS( 2) = g_2 GSS( 3) = g_3 !C GSS( 4) = y_tau GSS( 5) = y_b GSS( 6) = y_t !C GSS( 7) = M_1 GSS( 8) = M_2 GSS( 9) = M_3 !C GSS(10) = A_tau GSS(11) = A_b GSS(12) = A_t !C GSS(13) = M_hd^2 GSS(14) = M_hu^2 GSS(15) = M_er^2 !C GSS(16) = M_el^2 GSS(17) = M_dnr^2 GSS(18) = M_upr^2 !C GSS(19) = M_upl^2 GSS(20) = M_taur^2 GSS(21) = M_taul^2 !C GSS(22) = M_btr^2 GSS(23) = M_tpr^2 GSS(24) = M_tpl^2 !C GSS(25) = mu GSS(26) = B GSS(27) = Y_N !C GSS(28) = M_nr GSS(29) = A_n GSS(30) = vdq !C GSS(31) = vuq !C Masses: !C MSS( 1) = glss MSS( 2) = upl MSS( 3) = upr !C MSS( 4) = dnl MSS( 5) = dnr MSS( 6) = stl !C MSS( 7) = str MSS( 8) = chl MSS( 9) = chr !C MSS(10) = b1 MSS(11) = b2 MSS(12) = t1 !C MSS(13) = t2 MSS(14) = nuel MSS(15) = numl !C MSS(16) = nutl MSS(17) = el- MSS(18) = er- !C MSS(19) = mul- MSS(20) = mur- MSS(21) = tau1 !C MSS(22) = tau2 MSS(23) = z1ss MSS(24) = z2ss !C MSS(25) = z3ss MSS(26) = z4ss MSS(27) = w1ss !C MSS(28) = w2ss MSS(29) = hl0 MSS(30) = hh0 !C MSS(31) = ha0 MSS(32) = h+ !C Unification: !C MGUTSS = M_GUT GGUTSS = g_GUT AGUTSS = alpha_GUT ! COMMON /SUGMG/ MSS(32),GSS(31),MGUTSS,GGUTSS,AGUTSS,FTGUT, ! $FBGUT,FTAGUT,FNGUT ! REAL MSS,GSS,MGUTSS,GGUTSS,AGUTSS,FTGUT,FBGUT,FTAGUT,FNGUT ! SAVE /SUGMG/ ! COMMON /SUGPAS/ XTANB,MSUSY,AMT,MGUT,MU,G2,GP,V,VP,XW, ! $A1MZ,A2MZ,ASMZ,FTAMZ,FBMZ,B,SIN2B,FTMT,G3MT,VEV,HIGFRZ, ! $FNMZ,AMNRMJ,NOGOOD,IAL3UN,ITACHY,MHPNEG,ASM3, ! $VUMT,VDMT,ASMTP,ASMSS,M3Q ! REAL XTANB,MSUSY,AMT,MGUT,MU,G2,GP,V,VP,XW, ! $A1MZ,A2MZ,ASMZ,FTAMZ,FBMZ,B,SIN2B,FTMT,G3MT,VEV,HIGFRZ, ! $FNMZ,AMNRMJ,ASM3,VUMT,VDMT,ASMTP,ASMSS,M3Q ! INTEGER NOGOOD,IAL3UN,ITACHY,MHPNEG ! SAVE /SUGPAS/ !C SUSY parameters !C AMGLSS = gluino mass !C AMULSS = up-left squark mass !C AMELSS = left-selectron mass !C AMERSS = right-slepton mass !C AMNiSS = sneutrino mass for generation i !C TWOM1 = Higgsino mass = - mu !C RV2V1 = ratio v2/v1 of vev's !C AMTLSS,AMTRSS = left,right stop masses !C AMT1SS,AMT2SS = light,heavy stop masses !C AMBLSS,AMBRSS = left,right sbottom masses !C AMB1SS,AMB2SS = light,heavy sbottom masses !C AMLLSS,AMLRSS = left,right stau masses !C AML1SS,AML2SS = light,heavy stau masses !C AMZiSS = signed mass of Zi !C ZMIXSS = Zi mixing matrix !C AMWiSS = signed Wi mass !C GAMMAL,GAMMAR = Wi left, right mixing angles !C AMHL,AMHH,AMHA = neutral Higgs h0, H0, A0 masses !C AMHC = charged Higgs H+ mass !C ALFAH = Higgs mixing angle !C AAT = stop trilinear term !C THETAT = stop mixing angle !C AAB = sbottom trilinear term !C THETAB = sbottom mixing angle !C AAL = stau trilinear term !C THETAL = stau mixing angle !C AMGVSS = gravitino mass !C MTQ = top mass at weak scale !C MBQ = bottom mass at weak scale !C MLQ = tau mass at weak scale !C FBMA = b-Yukawa at mA scale !C VUQ = Hu vev at MSUSY !C VDQ = Hd vev at MSUSY ! COMMON/SSPAR/AMGLSS,AMULSS,AMURSS,AMDLSS,AMDRSS,AMSLSS ! $,AMSRSS,AMCLSS,AMCRSS,AMBLSS,AMBRSS,AMB1SS,AMB2SS ! $,AMTLSS,AMTRSS,AMT1SS,AMT2SS,AMELSS,AMERSS,AMMLSS,AMMRSS ! $,AMLLSS,AMLRSS,AML1SS,AML2SS,AMN1SS,AMN2SS,AMN3SS ! $,TWOM1,RV2V1,AMZ1SS,AMZ2SS,AMZ3SS,AMZ4SS,ZMIXSS(4,4) ! $,AMW1SS,AMW2SS ! $,GAMMAL,GAMMAR,AMHL,AMHH,AMHA,AMHC,ALFAH,AAT,THETAT ! $,AAB,THETAB,AAL,THETAL,AMGVSS,MTQ,MBQ,MLQ,FBMA, ! $VUQ,VDQ ! REAL AMGLSS,AMULSS,AMURSS,AMDLSS,AMDRSS,AMSLSS ! $,AMSRSS,AMCLSS,AMCRSS,AMBLSS,AMBRSS,AMB1SS,AMB2SS ! $,AMTLSS,AMTRSS,AMT1SS,AMT2SS,AMELSS,AMERSS,AMMLSS,AMMRSS ! $,AMLLSS,AMLRSS,AML1SS,AML2SS,AMN1SS,AMN2SS,AMN3SS ! $,TWOM1,RV2V1,AMZ1SS,AMZ2SS,AMZ3SS,AMZ4SS,ZMIXSS ! $,AMW1SS,AMW2SS ! $,GAMMAL,GAMMAR,AMHL,AMHH,AMHA,AMHC,ALFAH,AAT,THETAT ! $,AAB,THETAB,AAL,THETAL,AMGVSS,MTQ,MBQ,MLQ,FBMA,VUQ,VDQ ! REAL AMZISS(4) ! EQUIVALENCE (AMZISS(1),AMZ1SS) ! SAVE /SSPAR/ !C XNUSUG contains non-universal GUT scale soft terms for SUGRA: !C XNUSUG(1)=M1 XNUSUG(2)=M2 XNUSUG(3)=M3 !C XNUSUG(4)=A_tau XNUSUG(5)=A_b XNUSUG(6)=A_t !C XNUSUG(7)=m_Hd XNUSUG(8)=m_Hu XNUSUG(9)=m_eR XNUSUG(10)=m_eL !C XNUSUG(11)=m_dR XNUSUG(12)=m_uR XNUSUG(13)=m_uL XNUSUG(14)=m_lR !C XNUSUG(15)=m_lL XNUSUG(16)=m_bR XNUSUG(17)=m_tR XNUSUG(18)=m_tL !C ! COMMON /SUGNU/ XNUSUG(18) ! REAL XNUSUG ! SAVE /SUGNU/ !C Standard model parameters !C AMUP,...,AMTP = quark masses !C AME,AMMU,AMTAU = lepton masses !C AMW,AMZ = W,Z masses !C GAMW,GAMZ = W,Z widths !C ALFAEM,SN2THW,ALFA3 = SM couplings !C ALQCD4 = 4 flavor lambda ! COMMON/SSSM/AMUP,AMDN,AMST,AMCH,AMBT,AMTP,AME,AMMU,AMTAU ! $,AMW,AMZ,GAMW,GAMZ,ALFAEM,SN2THW,ALFA2,ALFA3,ALQCD4 ! REAL AMUP,AMDN,AMST,AMCH,AMBT,AMTP,AME,AMMU,AMTAU ! $,AMW,AMZ,GAMW,GAMZ,ALFAEM,SN2THW,ALFA2,ALFA3,ALQCD4 ! SAVE /SSSM/ ! REAL PI,SQRT2,ALFA,GF,UNITS ! COMMON/CONST/PI,SQRT2,ALFA,GF,UNITS ! SAVE /CONST/ CsB End of ISAJET common blocks REAL GPX,SIN2W,ALEMI,AS,ASMB,MBMB INTEGER IMODEL,J,K,IMODIN CsB Local ISAJET related variables CHARACTER*40 VERSN,VISAJE CsB Local LHA3 related variables Integer iPDG(33),ISA2LHA3(33),I2L3GSS(50) DIMENSION CHAF(33),ModelDescr(7),SoftParaLHA(50) CHARACTER CHAF*16, ModelDescr*50,SoftParaLHA*16 CsB These are the particle names the masses for which LHA3 wants C (in this order) DATA CHAF / $' top',' h^0',' H^0',' A^0',' H^+', $' dnl',' upl',' stl',' chl',' b1',' t1', $' el-',' nuel',' mul-',' numl',' tau1',' nutl', $' glss',' z1ss',' z2ss',' w1ss',' z3ss',' z4ss',' w2ss', $' dnr',' upr',' str',' chr',' b2',' t2', $' er-',' mur-',' tau2' / CsB These are the PDG codes of the above DATA iPDG / & 6, 25, 35, 36, 37, &1000001,1000002,1000003,1000004,1000005,1000006, &1000011,1000012,1000013,1000014,1000015,1000016, &1000021,1000022,1000023,1000024,1000025,1000035,1000037, &2000001,2000002,2000003,2000004,2000005,2000006, &2000011,2000013,2000015 / CsB These are the MSS indices of the above Data ISA2LHA3 / $ 0,29,30,31,32, $ 4, 2, 6, 8,10,12, $17,14,19,15,21,16, $ 1,23,24,27,25,26,28, $ 5, 3, 7, 9,11,13, $18,20,22/ CsB These are the soft parameters of which LHA3 wants (in this order) DATA SoftParaLHA / ,'M_1(Q)','M_2(Q)','M_3(Q)',' ',' ', ,' ',' ',' ',' ',' ', ,' ',' ',' ',' ',' ', ,' ',' ',' ',' ',' ', ,' ',' ',' ',' ',' ', ,' ',' ',' ',' ',' ', ,'MeL(Q)','MmuL(Q)','MtauL(Q)','MeR(Q)','MmuR(Q)', ,'MtauR(Q)',' ',' ',' ',' ', ,'MqL1(Q)','MqL2(Q)','MqL3(Q)','MuR(Q)','McR(Q)', ,'MtR(Q)','MdR(Q)','MsR(Q)','MbR(Q)',' '/ CsB These are GSS indices of the above Data I2L3GSS / $ 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, $ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, $ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, $ 16,16,21,15,15,20, 0, 0, 0, 0, $ 19,19,24,18,18,23,17,17,22, 0/ Data ModelDescr / $'Minimal supergravity (mSUGRA,CMSSM) model', $'Minimal gauge mediated (GMSB) model', $'Non-universal supergravity model', $'Supergravity model with truly unified couplings', $'Non-minimal gauge mediated (GMSB) model', $'Supergravity model with right-handed neutrinos', $'Anomaly-mediated SUSY breaking model' / Logical Testing Testing = .False. If (iModIn.NE.1) then Print 'This ISALHA version only writes an mSUGRA spectrum' Return End If C C Entry C PI=4.*ATAN(1.) GPX=SQRT(.6)*GSS(1) SIN2W=GPX**2/(GSS(2)**2+GPX**2) ALEMI=4*PI/GSS(2)**2/SIN2W AS=GSS(3)**2/4./PI CsB Open LHA3 output file LOUT = 91 Open(LOUT,FILE='ISALHA3.out',FORM='FORMATTED') CsB Write LHA3 header WRITE(LOUT,7000) . ' ISAJET SUSY parameters in Les Houches accord format' WRITE(LOUT,7000) . ' Created by ISALHA3. Last revision by C. Balazs on 2005 Jan. 5' VERSN=VISAJE() VERSN=VERSN(14:) WRITE(LOUT,7001) 'SPINFO', , 'Program information' WRITE(LOUT,7012) 1, 'ISASUGRA from ISAJET ', , 'Spectrum Calculator' WRITE(LOUT,7012) 2, VERSN, , 'Version number' WRITE(LOUT,7001) 'MODSEL', 'Model selection' WRITE(LOUT,7010) 1, IMODIN, ModelDescr(IMODIN) Call SETCON ASMB=SUALFS(AMBT**2,.36,AMTP,3) MBMB=AMBT*(1.-4*ASMB/3./PI) WRITE(LOUT,7001) 'SMINPUTS', 'Standard Model inputs' WRITE(LOUT,7011) 1, ALEMI, 'alpha_em^(-1)' ! (MZ) SM MSbar' WRITE(LOUT,7011) 2, GF, 'G_Fermi' WRITE(LOUT,7011) 3, AS, 'alpha_s(M_Z)' WRITE(LOUT,7011) 4, AMZ, 'm_{Z}(pole)' ! WRITE(LOUT,7011) 5, AMBT, 'm_{b}(pole)' WRITE(LOUT,7011) 5, MBMB, 'm_{b}(m_{b})' WRITE(LOUT,7011) 6, AMTP, 'm_{top}(pole)' WRITE(LOUT,7011) 7, AMTAU, 'm_{tau}(pole)' WRITE(LOUT,7001) 'MINPAR', 'SUSY breaking input parameters' C C Print inputs and GUT couplings for SUGRA/AMSB models C IF(IMODEL.EQ.1.OR.IMODEL.EQ.7) THEN IF(IMODEL.EQ.1) THEN WRITE(LOUT,7011) 1, XSUGIN(1), 'm_0' WRITE(LOUT,7011) 2, XSUGIN(2), 'm_{1/2}' WRITE(LOUT,7011) 3, XSUGIN(4), 'tan(beta)' WRITE(LOUT,7011) 4, XSUGIN(5), 'sign(mu)' WRITE(LOUT,7011) 5, XSUGIN(3), 'A_0' ! WRITE(LOUT,1000) XSUGIN(1),XSUGIN(2),XSUGIN(3),XSUGIN(4), ! $ XSUGIN(5),XSUGIN(6) !1000 FORMAT( ! $ ' M_0, M_(1/2), A_0, tan(beta), sgn(mu), M_t =' ! $ /4F10.3,2X,F6.1,F10.3) ELSE IF (IMODEL.EQ.7) THEN WRITE(LOUT,1018) XSUGIN(1),XSUGIN(2),XSUGIN(4),XSUGIN(5), $ XSUGIN(6) 1018 FORMAT( $ ' M_0, M_(3/2), tan(beta), sgn(mu), M_t =' $ /3F10.3,2X,F6.1,2F10.3) END IF CsB Goto 1234 C C Write out non-universal GUT scale parameters IF(XNUSUG(1).LT.1.E19.OR.XNUSUG(2).LT.1.E19.OR.XNUSUG(3) $ .LT.1.E19) THEN WRITE(LOUT,1010) XNUSUG(1),XNUSUG(2),XNUSUG(3) 1010 FORMAT(/' M_1(GUT)= ',F8.2,' M_2(GUT)= ',F8.2, $ ' M_3(GUT)= ',F8.2) END IF IF(XNUSUG(4).LT.1.E19.OR.XNUSUG(5).LT.1.E19.OR.XNUSUG(6) $ .LT.1.E19) THEN WRITE(LOUT,1011) XNUSUG(4),XNUSUG(5),XNUSUG(6) 1011 FORMAT(/' A_tau(GUT)= ',F8.2,' A_b(GUT)= ',F8.2, $ ' A_t(GUT)= ',F8.2) END IF IF(XNUSUG(7).LT.1.E19.OR.XNUSUG(8).LT.1.E19) THEN WRITE(LOUT,1012) XNUSUG(7),XNUSUG(8) 1012 FORMAT(/' M_Hd(GUT)= ',F8.2,' M_Hu(GUT)= ',F8.2) END IF IF (XNUSUG(9).LT.1.E19.OR.XNUSUG(10).LT.1.E19) THEN WRITE(LOUT,1013) XNUSUG(9),XNUSUG(10) 1013 FORMAT(/' M_eR(GUT)= ',F8.2,' M_eL(GUT)= ',F8.2) END IF IF(XNUSUG(11).LT.1.E19.OR.XNUSUG(12).LT.1.E19.OR.XNUSUG(13) $ .LT.1.E19) THEN WRITE(LOUT,1014) XNUSUG(11),XNUSUG(12),XNUSUG(13) 1014 FORMAT(' M_dR(GUT)= ',F8.2,' M_uR(GUT)= ',F8.2, $ ' M_uL(GUT)=',F8.2) END IF IF(XNUSUG(14).LT.1.E19.OR.XNUSUG(15).LT.1.E19) THEN WRITE(LOUT,1015) XNUSUG(14),XNUSUG(15) 1015 FORMAT(/' M_tauR(GUT)= ',F8.2,' M_tauL(GUT)= ',F8.2) END IF IF(XNUSUG(16).LT.1.E19.OR.XNUSUG(17).LT.1.E19.OR.XNUSUG(18) $ .LT.1.E19) THEN WRITE(LOUT,1016) XNUSUG(16),XNUSUG(17),XNUSUG(18) 1016 FORMAT(' M_bR(GUT)= ',F8.2,' M_tR(GUT)= ',F8.2, $ ' M_tL(GUT)=',F8.2) END IF IF(XSUGIN(7).NE.0) THEN WRITE(LOUT,1026) XSUGIN(7) 1026 FORMAT(' Q_max= ',E12.4) ENDIF C C Right-handed neutrino parameters IF (XNRIN(2).LT.1.E19) THEN WRITE(LOUT,1017) XNRIN(1),XNRIN(2),XNRIN(3),XNRIN(4), $ FNMZ,FNGUT 1017 FORMAT(' Right-handed neutrino parameters:'/ $ ' M(nu_tau)=',E10.3,' M(N_R) =',E10.3, $ ' A_N=',F8.2,' M(NRSS)=',F8.2/ $ ' FN(M_Z) =',F8.4, ' FN(M_GUT) =',F8.4) END IF C C Unification results WRITE(LOUT,1001) MGUTSS,GGUTSS,AGUTSS 1001 FORMAT(/' ISASUGRA unification:'/' M_GUT =',E10.3, $ ' g_GUT =',F5.3,3X,' alpha_GUT =',F5.3) WRITE(LOUT,999) FTGUT,FBGUT,FTAGUT 999 FORMAT(' FT_GUT =',F6.3, $ ' FB_GUT =',F6.3,3X,' FL_GUT =',F6.3) C C Print inputs for GMSB models C ELSE IF (IMODEL.EQ.2) THEN WRITE(LOUT,1002) (XGMIN(J),J=1,7) 1002 FORMAT( $ ' Lambda, M_mes, N_5, tan(beta), sgn(mu), M_t, C_grav=' $ /2E10.3,2F10.3,2X,F6.1,F10.3,1X,E10.3) WRITE(LOUT,1020) (XGMIN(J),J=8,14) 1020 FORMAT(/' GMSB2 model input:'/ $ ' Rsl, dmH_d^2, dmH_u^2, d_Y, N5_1, N5_2, N5_3=' $ /F7.3,1X,E10.3,1X,E10.3,1X,E10.3,2X,3F7.3) WRITE(LOUT,1003) AMGVSS 1003 FORMAT(/' M(gravitino)=',E10.3) END IF C C Weak scale couplings C WRITE(LOUT,1004) ALEMI,SIN2W,AS 1004 FORMAT(/' 1/alpha_em =',F8.2,2X, $' sin**2(thetaw) =',F6.4,2X,' alpha_s =',F5.3) WRITE(LOUT,1005) GSS(7),GSS(8),GSS(9) 1005 FORMAT(' M_1 =',F8.2,2X, $' M_2 =',F8.2,' M_3 =',F8.2) WRITE(LOUT,1006) MU,B,HIGFRZ 1006 FORMAT(' mu(Q) =',F8.2,2X, $' B(Q) =',F8.2,' Q =',F8.2) WRITE(LOUT,1007) GSS(13),GSS(14) 1007 FORMAT(' M_H1^2 =',E10.3,' M_H2^2 =',E10.3) 1234 Continue C C Print mass spectrum from ISASUGRA C WRITE(LOUT,7000) ' ' WRITE(LOUT,6999) ' M_{GUT} =', MGUTSS WRITE(LOUT,7001) 'MASS', 'Scalar and gaugino mass spectrum' WRITE(LOUT,7000) ' PDG code mass particle' If (Testing) then WRITE(LOUT,2000) MSS(1),MSS(2),MSS(3),MSS(4),MSS(5),MSS(10), $ MSS(11),MSS(12),MSS(13),MSS(14),MSS(17),MSS(18),MSS(16), $ MSS(21),MSS(22),MSS(23),MSS(24),MSS(25),MSS(26),MSS(27), $ MSS(28),MSS(29),MSS(30),MSS(31),MSS(32) 2000 FORMAT(/' ISAJET masses (with signs):'/ $ ' M(GL) =',F9.2/ $ ' M(UL) =',F9.2,' M(UR) =',F9.2,' M(DL) =',F9.2, $ ' M(DR) =',F9.2/ $ ' M(B1) =',F9.2,' M(B2) =',F9.2,' M(T1) =',F9.2, $ ' M(T2) =',F9.2/ $ ' M(SN) =',F9.2,' M(EL) =',F9.2,' M(ER) =',F9.2/ $ ' M(NTAU)=',F9.2,' M(TAU1)=',F9.2,' M(TAU2)=',F9.2/ $ ' M(Z1) =',F9.2,' M(Z2) =',F9.2,' M(Z3) =',F9.2, $ ' M(Z4) =',F9.2/ $ ' M(W1) =',F9.2,' M(W2) =',F9.2/ $ ' M(HL) =',F9.2,' M(HH) =',F9.2,' M(HA) =',F9.2, $ ' M(H+) =',F9.2) EndIf ! WRITE(LOUT,7013) iPDG(1), MT, CHAF(1) WRITE(LOUT,7013) 24, AMW, ' W^+' DO 370 I=2,33 sg = 1. CsB The signs of the (EW) gaugino masses are flipped according to ISAWIG If (iPDG(I).Eq.1000022 .or. iPDG(I).Eq.1000023 .or. . iPDG(I).Eq.1000024 .or. iPDG(I).Eq.1000025 .or. . iPDG(I).Eq.1000035 .or. iPDG(I).Eq.1000037) sg = -1. WRITE(LOUT,7013) iPDG(I), sg*MSS(ISA2LHA3(I)), CHAF(I) 370 CONTINUE C SUSY scale RMSUSY = HIGFRZ !!! check this WRITE(LOUT,7000) ' Higgs mixing' WRITE(LOUT,7001) 'ALPHA','Effective Higgs mixing parameter' WRITE(LOUT,7016) -ALFAH, 'alpha' ! Sign flips for LHA3 If (Testing) then WRITE(LOUT,2001) THETAT,THETAB,THETAL,ALFAH 2001 FORMAT(/,' theta_t=',F9.4,' theta_b=',F9.4, $ ' theta_l=',F9.4,' alpha_h=',F9.4) EndIf C C Write out chargino /neutralino masses/eigenvectors C If (Testing) then WRITE(LOUT,3100) AMZ1SS,AMZ2SS,AMZ3SS,AMZ4SS 3100 FORMAT(/' NEUTRALINO MASSES (SIGNED) =',4F10.3) DO 100 J=1,4 WRITE(LOUT,3200) J,(ZMIXSS(K,J),K=1,4) 3200 FORMAT(' EIGENVECTOR ',I1,' =',4F10.5) 100 CONTINUE WRITE(LOUT,3300) AMW1SS,AMW2SS 3300 FORMAT(/' CHARGINO MASSES (SIGNED) =',2F10.3) WRITE(LOUT,3400) GAMMAL,GAMMAR 3400 FORMAT(' GAMMAL, GAMMAR =',2F10.5/) EndIf CsB For the mixing matrices I follow ISAWIG1200 to the letter WRITE(LOUT,7001) 'STOPMIX','stop mixing matrix' WRITE(LOUT,7021) 1, 1, COS(THETAT), 'O_{11}' WRITE(LOUT,7021) 1, 2, -SIN(THETAT), 'O_{12}' WRITE(LOUT,7021) 2, 1, SIN(THETAT), 'O_{21}' WRITE(LOUT,7021) 2, 2, COS(THETAT), 'O_{22}' WRITE(LOUT,7001) 'SBOTMIX','sbottom mixing matrix' WRITE(LOUT,7021) 1, 1, COS(THETAB), 'O_{11}' WRITE(LOUT,7021) 1, 2, -SIN(THETAB), 'O_{12}' WRITE(LOUT,7021) 2, 1, SIN(THETAB), 'O_{21}' WRITE(LOUT,7021) 2, 2, COS(THETAB), 'O_{22}' WRITE(LOUT,7001) 'STAUMIX','stau mixing matrix' WRITE(LOUT,7021) 1, 1, COS(THETAL), 'O_{11}' WRITE(LOUT,7021) 1, 2, -SIN(THETAL), 'O_{12}' WRITE(LOUT,7021) 2, 1, SIN(THETAL), 'O_{21}' WRITE(LOUT,7021) 2, 2, COS(THETAL), 'O_{22}' WRITE(LOUT,7001) 'NMIX','neutralino mixing matrix' CsB ... in ascending mass order (rows) and in the order C (bino, w3ino, higgs1, higgs2) (columns) DO I1=1,4 DO I2=1,4 sg = 1. If (I2.GT.2) sg = -1. J1 = 5 - I2 J2 = I1 WRITE(LOUT,7021) I1, I2, sg*ZMIXSS(J1,J2) EndDo EndDo THX=SIGN(1.,1./TAN(GAMMAL)) THY=SIGN(1.,1./TAN(GAMMAR)) WRITE(LOUT,7001) 'UMIX','chargino U mixing matrix' WRITE(LOUT,7021) 1, 1, -1.0*SIN(GAMMAL), 'U_{11}' WRITE(LOUT,7021) 1, 2, -1.0*COS(GAMMAL), 'U_{12}' WRITE(LOUT,7021) 2, 1, -THX*COS(GAMMAL), 'U_{21}' WRITE(LOUT,7021) 2, 2, THX*SIN(GAMMAL), 'U_{22}' WRITE(LOUT,7001) 'VMIX','chargino V mixing matrix' WRITE(LOUT,7021) 1, 1, -1.0*SIN(GAMMAR), 'V_{11}' WRITE(LOUT,7021) 1, 2, -1.0*COS(GAMMAR), 'V_{12}' WRITE(LOUT,7021) 2, 1, -THY*COS(GAMMAR), 'V_{21}' WRITE(LOUT,7021) 2, 2, THY*SIN(GAMMAR), 'V_{22}' WRITE(LOUT,7002) 'GAUGE',RMSUSY !!! check: are these at Q=RMSUSY? WRITE(LOUT,7011) 1, SQRT(.6)*GSS(1), 'g`' WRITE(LOUT,7011) 2, GSS(2), 'g_2' WRITE(LOUT,7011) 3, GSS(3), 'g_3' WRITE(LOUT,7002) 'YU',RMSUSY WRITE(LOUT,7021) 3, 3, GSS( 6), 'y_t' !!! check: are these at Q=RMSUSY? WRITE(LOUT,7002) 'YD',RMSUSY WRITE(LOUT,7021) 3, 3, GSS( 5), 'y_b' WRITE(LOUT,7002) 'YE',RMSUSY WRITE(LOUT,7021) 3, 3, GSS( 4), 'y_tau' WRITE(LOUT,7002) 'HMIX',RMSUSY,'Higgs mixing parameters' !!! check: are these at Q=RMSUSY? WRITE(LOUT,7011) 1, MU, 'mu(Q)' WRITE(LOUT,7011) 2, XtanB, 'tan(beta)(M_{GUT})' WRITE(LOUT,7011) 3, Sqrt(2.)*VEV, 'Higgs vev at Q' WRITE(LOUT,7011) 4, MSS(31)**2, 'm_A^2(Q)' WRITE(LOUT,7002) 'MSOFT',RMSUSY, , 'DRbar SUSY breaking parameters' !!! check: are these at Q=RMSUSY? Do I=1,3 If (I2L3GSS(I).NE.0) , WRITE(LOUT,7011) I, GSS(I2L3GSS(I)), SoftParaLHA(I) End Do Do I=4,50 If (I2L3GSS(I).NE.0) !!! Fix sign - if necessary , WRITE(LOUT,7011) I, Sqrt(Abs(GSS(I2L3GSS(I)))), SoftParaLHA(I) End Do WRITE(LOUT,7002) 'AU',RMSUSY WRITE(LOUT,7021) 1, 1, GSS(12), 'A_u' WRITE(LOUT,7021) 2, 2, GSS(12), 'A_c' WRITE(LOUT,7021) 3, 3, GSS(12), 'A_t' WRITE(LOUT,7002) 'AD',RMSUSY WRITE(LOUT,7021) 1, 1, GSS(11), 'A_d' WRITE(LOUT,7021) 2, 2, GSS(11), 'A_s' WRITE(LOUT,7021) 3, 3, GSS(11), 'A_b' WRITE(LOUT,7002) 'AE',RMSUSY WRITE(LOUT,7021) 1, 1, GSS(10), 'A_e' WRITE(LOUT,7021) 2, 2, GSS(10), 'A_mu' WRITE(LOUT,7021) 3, 3, GSS(10), 'A_tau' C C Print ISAJET MSSMi equivalent input C If (Testing) then WRITE(LOUT,3000) 3000 FORMAT(/' ISAJET equivalent input:') WRITE(LOUT,3001) MSS(1),MU,MSS(31),XSUGIN(4) 3001 FORMAT(' MSSMA: ',4F8.2) WRITE(LOUT,3002) SQRT(GSS(19)),SQRT(GSS(17)),SQRT(GSS(18)), $ SQRT(GSS(16)),SQRT(GSS(15)) 3002 FORMAT(' MSSMB: ',5F8.2) WRITE(LOUT,3003) SIGN(1.,GSS(24))*SQRT(ABS(GSS(24))), $ SQRT(GSS(22)),SIGN(1.,GSS(23))*SQRT(ABS(GSS(23))), $ SQRT(GSS(21)),SQRT(GSS(20)),GSS(12),GSS(11),GSS(10) 3003 FORMAT(' MSSMC: ',8F8.2) WRITE(LOUT,3004) 3004 FORMAT(' MSSMD: SAME AS MSSMB (DEFAULT)') WRITE(LOUT,3005) GSS(7),GSS(8) 3005 FORMAT(' MSSME: ',2F8.2) EndIf Close(91) CsB LHA3 format statements C Formats for user information printout. 5000 FORMAT(1x,17('*'),1x,'ISALHA3 v0.1: SUSY SPECTRUM ' & ,'INTERFACE',1x,17('*')/1x,'*',3x & ,'ISALHA3: Last Change',1x,A,1x,'-',1x,'C. Balazs') 5001 FORMAT(1x,'*',3x,'Writing spectrum file on unit: ',I3) 5002 FORMAT(1x,'*',3x,'Reading spectrum file on unit: ',I3) 5003 FORMAT(1x,'*',3x,'Spectrum Calculator was: ',A,' version ',A) 5100 FORMAT(1x,'*',1x,'Model parameters:'/1x,'*',1x,'----------------') 5200 FORMAT(1x,'*',1x,3x,'m_0',6x,'m_{1/2}',5x,'A_0',3x,'tan(beta)', & 3x,'sgn(mu)',3x,'m_t'/1x,'*',1x,4(F8.2,1x),I8,2x,F8.2) 5300 FORMAT(1x,'*'/1x,'*',1x,'Model spectrum :'/1x,'*',1x & ,'----------------') 5400 FORMAT(1x,'*',1x,A) 5500 FORMAT(1x,'*',1x,A,':') 5600 FORMAT(1x,'*',2x,2x,'M_GUT',2x,2x,'g_GUT',2x,1x,'alpha_GUT'/ & 1x,'*',2x,1P,2(1x,E8.2),2x,E8.2) 5700 FORMAT(1x,'*',4x,4x,'~d',2x,1x,4x,'~u',2x,1x,4x,'~s',2x,1x, & 4x,'~c',2x,1x,1x,'~b(12)',1x,1x,1x,'~t(12)'/1x,'*',2x,'L',1x & ,6(F8.2,1x)/1x,'*',2x,'R',1x,6(F8.2,1x)) 5800 FORMAT(1x,'*'/1x,'*',4x,4x,'~e',2x,1x,2x,'~nu_e',2x,1x,3x,'~mu',2x & ,1x,1x,'~nu_mu',1x,1x,'~tau(12)',1x,1x,'~nu_tau'/1x,'*',2x & ,'L',1x,6(F8.2,1x)/1x,'*',2x,'R',1x,6(F8.2,1x)) 5900 FORMAT(1x,'*'/1x,'*',4x,4x,'~g',2x,1x,1x,'~chi_10',1x,1x,'~chi_20' & ,1x,1x,'~chi_30',1x,1x,'~chi_40',1x,1x,'~chi_1+',1x & ,1x,'~chi_2+'/1x,'*',3x,1x,7(F8.2,1x)) 6000 FORMAT(1x,'*'/1x,'*',4x,4x,'h0',2x,1x,4x,'H0',2x,1x,4x,'A0',2x & ,1x,4x,'H+'/1x,'*',3x,1x,5(F8.2,1x)) 6100 FORMAT(1x,'*',11x,'|',3x,'~B',3x,'|',2x,'~W_3',2x,'|',2x & ,'~H_1',2x,'|',2x,'~H_2',2x,'|'/1x,'*',3x,'~chi_10',1x,4('|' & ,1x,F6.3,1x),'|'/1x,'*',3x,'~chi_20',1x,4('|' & ,1x,F6.3,1x),'|'/1x,'*',3x,'~chi_30',1x,4('|' & ,1x,F6.3,1x),'|'/1x,'*',3x,'~chi_40',1x,4('|' & ,1x,F6.3,1x),'|') 6200 FORMAT(1x,'*'/1x,'*',6x,'L',4x,'|',3x,'~W',3x,'|',3x,'~H',3x,'|' & ,12x,'R',4x,'|',3x,'~W',3x,'|',3x,'~H',3x,'|'/1x,'*',3x & ,'~chi_1+',1x,2('|',1x,F6.3,1x),'|',9x,'~chi_1+',1x,2('|',1x & ,F6.3,1x),'|'/1x,'*',3x,'~chi_2+',1x,2('|',1x,F6.3,1x),'|',9x & ,'~chi_2+',1x,2('|',1x,F6.3,1x),'|') 6300 FORMAT(1x,'*'/1x,'*',8x,'|',2x,'~b_L',2x,'|',2x,'~b_R',2x,'|',8x & ,'|',2x,'~t_L',2x,'|',2x,'~t_R',2x,'|',10x & ,'|',1x,'~tau_L',1x,'|',1x,'~tau_R',1x,'|'/ & 1x,'*',3x,'~b_1',1x,2('|',1x,F6.3,1x),'|',3x,'~t_1',1x,2('|' & ,1x,F6.3,1x),'|',3x,'~tau_1',1x,2('|',1x,F6.3,1x),'|'/ & 1x,'*',3x,'~b_2',1x,2('|',1x,F6.3,1x),'|',3x,'~t_2',1x,2('|' & ,1x,F6.3,1x),'|',3x,'~tau_2',1x,2('|',1x,F6.3,1x),'|') 6400 FORMAT(1x,'*',3x,' A_b = ',F8.2,4x,' A_t = ',F8.2,4x & ,'A_tau = ',F8.2) 6450 FORMAT(1x,'*',3x,'alpha = ',F8.2,4x,'tan(beta) = ',F8.2,4x & ,' mu = ',F8.2) 6500 FORMAT(1x,32('*'),1x,'END OF ISALHA3',1x,31('*')) C Format to use for comments 6999 FORMAT('# ',A,1x,E16.8) 7000 FORMAT('# ',A) C Format to use for block statements 7001 FORMAT('Block',1x,A,3x,'#',1x,A) 7002 FORMAT('Block',1x,A,1x,'Q=',1P,E16.8,0P,3x,'#',1x,A) C Indexed Int 7010 FORMAT(1x,I5,1x,I5,3x,'#',1x,A) C Indexed Double 7011 FORMAT(1x,I5,3x,1P,E16.8,0P,3x,'#',1x,A) C Indexed Char(12) 7012 FORMAT(1x,I5,3x,A27,3x,'#',1x,A) C Long Indexed Double 7013 FORMAT(1x,I9,3x,1P,E16.8,0P,3x,'#',1x,A) C Indexed Double without leading integer 7016 FORMAT(8x,1P,E16.8,0P,3x,'#',1x,A) C Double Matrix 7022 FORMAT(1x,I2,1x,I2,3x,1P,E16.8,3x,E16.8,0P,3x,'#',1x,A) C Single matrix 7021 FORMAT(1x,I2,1x,I2,3x,1P,E16.8,0P,3x,'#',1x,A) C Write Decay Table 7500 FORMAT('Decay',1x,I9,1x,'WIDTH=',1P,E16.8,0P,3x,'#',1x,A) 7501 FORMAT(4x,1P,E16.8,0P,3x,I2,3x,'IDA=',1x,5(1x,I9),3x,'#',1x,A) RETURN END