Back to home page

sPhenix code displayed by LXR
 
 

    

File indexing completed on 2025-08-05 08:21:18 

0001  
0002 C*********************************************************************
0003  
0004 C...PYSUGI
0005 C...Interface to ISASUSY version 7.71.
0006 C...Warning: this interface should not be used with earlier versions
0007 C...of ISASUSY, since common block incompatibilities may then arise.
0008 C...Calls SUGRA (in ISAJET) to perform RGE evolution.
0009 C...Then converts to Gunion-Haber conventions.
0010  
0011       SUBROUTINE PYSUGI
0012       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0013  
0014       INTEGER PYK,PYCHGE,PYCOMP
0015       PARAMETER (KSUSY1=1000000,KSUSY2=2000000,KTECHN=3000000,
0016      &KEXCIT=4000000,KDIMEN=5000000)
0017  
0018 C...Date of Change
0019       CHARACTER DOC*11
0020       PARAMETER (DOC='01 May 2006')
0021  
0022 C...ISASUGRA Input:
0023       REAL MZERO,MHLF,AZERO,TANB,SGNMU,MTOP
0024 C...XISAIN contains the MSSMi inputs in natural order.
0025       COMMON /SUGXIN/ XISAIN(24),XSUGIN(7),XGMIN(14),XNRIN(4),
0026      $XAMIN(7)
0027       REAL XISAIN,XSUGIN,XGMIN,XNRIN,XAMIN
0028       SAVE /SUGXIN/
0029 C...ISASUGRA Output
0030       CHARACTER*40 ISAVER,VISAJE
0031       REAL SUPER
0032       COMMON /SSPAR/ SUPER(72)
0033       COMMON /SUGMG/ MSS(32),GSS(31),MGUTSS,GGUTSS,AGUTSS,FTGUT,
0034      $FBGUT,FTAGUT,FNGUT
0035       REAL MSS,GSS,MGUTSS,GGUTSS,AGUTSS,FTGUT,FBGUT,FTAGUT,FNGUT
0036       COMMON /SUGPAS/ XTANB,MSUSY,AMT,MGUT,MU,G2,GP,V,VP,XW,
0037      $A1MZ,A2MZ,ASMZ,FTAMZ,FBMZ,B,SIN2B,FTMT,G3MT,VEV,HIGFRZ,
0038      $FNMZ,AMNRMJ,NOGOOD,IAL3UN,ITACHY,MHPNEG,ASM3,
0039      $VUMT,VDMT,ASMTP,ASMSS,M3Q
0040       REAL XTANB,MSUSY,AMT,MGUT,MU,G2,GP,V,VP,XW,
0041      $A1MZ,A2MZ,ASMZ,FTAMZ,FBMZ,B,SIN2B,FTMT,G3MT,VEV,HIGFRZ,
0042      $FNMZ,AMNRMJ,ASM3,VUMT,VDMT,ASMTP,ASMSS,M3Q
0043       INTEGER NOGOOD,IAL3UN,ITACHY,MHPNEG
0044       INTEGER IALLOW
0045       SAVE /SUGMG/,/SSPAR/
0046 C SUPER: Filled by ISASUGRA.
0047 C SUPER(1)        = mass of ~g
0048 C SUPER(2:17)     = mass of ~u_L,~u_R,~d_L,~d_R,~s_L,~s_R,~c_L,~c_R,~b_L
0049 C                          ,~b_R,~b_1,~b_2,~t_L,~t_R,~t_1,~t_2
0050 C SUPER(18:25)    = mass of ~e_L,~e_R,~mu_L,~mu_R,~tau_L,~tau_R,~tau_1
0051 C                          ,~tau_2
0052 C SUPER(26:28)    = mass of ~nu_e,~nu_mu,~nu_tau
0053 C SUPER(29)       = Higgsino mass = - mu
0054 C SUPER(30)       = ratio v2/v1 of vev's
0055 C SUPER(31:34)    = Signed neutralino masses
0056 C SUPER(35:50)    = Neutralino mixing matrix
0057 C SUPER(51:52)    = Signed chargino masses
0058 C SUPER(53:54)    = Chargino left, right mixing angles
0059 C SUPER(55:58)    = mass of h0, H0, A0, H+
0060 C SUPER(59)       = Higgs mixing angle alpha
0061 C SUPER(60:65)    = A_t, theta_t, A_b, theta_b, A_tau, theta_tau
0062 C SUPER(66)       = Gravitino mass
0063 C SUPER(67:69)    = Top,Bottom, and Tau masses at MSUSY (not used)
0064 C SUPER(70)       = b-Yukawa at mA scale (not used)
0065 C SUPER(71:72)    = H_u, H_d vev's at MSUSY (not used)
0066 C GSS: Filled by ISASUGRA
0067 C     GSS( 1) = g_1        GSS( 2) = g_2        GSS( 3) = g_3
0068 C     GSS( 4) = y_tau      GSS( 5) = y_b        GSS( 6) = y_t
0069 C     GSS( 7) = M_1        GSS( 8) = M_2        GSS( 9) = M_3
0070 C     GSS(10) = A_tau      GSS(11) = A_b        GSS(12) = A_t
0071 C     GSS(13) = M_h12     GSS(14) = M_h22     GSS(15) = M_er2
0072 C     GSS(16) = M_el2     GSS(17) = M_dnr2    GSS(18) = M_upr2
0073 C     GSS(19) = M_upl2    GSS(20) = M_taur2   GSS(21) = M_taul2
0074 C     GSS(22) = M_btr2    GSS(23) = M_tpr2    GSS(24) = M_tpl2
0075 C     GSS(25) = mu         GSS(26) = B          GSS(27) = Y_N
0076 C     GSS(28) = M_nr       GSS(29) = A_n        GSS(30) = log(vdq)
0077 C     GSS(31) = log(vuq)
0078 C MSS: Filled by ISASUGRA
0079 C     MSS( 1) = glss     MSS( 2) = upl      MSS( 3) = upr
0080 C     MSS( 4) = dnl      MSS( 5) = dnr      MSS( 6) = stl
0081 C     MSS( 7) = str      MSS( 8) = chl      MSS( 9) = chr
0082 C     MSS(10) = b1       MSS(11) = b2       MSS(12) = t1
0083 C     MSS(13) = t2       MSS(14) = nuel     MSS(15) = numl
0084 C     MSS(16) = nutl     MSS(17) = el-      MSS(18) = er-
0085 C     MSS(19) = mul-     MSS(20) = mur-     MSS(21) = tau1
0086 C     MSS(22) = tau2     MSS(23) = z1ss     MSS(24) = z2ss
0087 C     MSS(25) = z3ss     MSS(26) = z4ss     MSS(27) = w1ss
0088 C     MSS(28) = w2ss     MSS(29) = hl0      MSS(30) = hh0
0089 C     MSS(31) = ha0      MSS(32) = h+
0090 C Unification, filled by ISASUGRA if applicable.
0091 C     MGUTSS  = M_GUT    GGUTSS  = g_GUT    AGUTSS  = alpha_GUTC
0092  
0093 C...SPYTHIA Input/Output
0094       INTEGER IMSS
0095       DOUBLE PRECISION RMSS
0096       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
0097       COMMON/PYSSMT/ZMIX(4,4),UMIX(2,2),VMIX(2,2),SMZ(4),SMW(2),
0098      &SFMIX(16,4),ZMIXI(4,4),UMIXI(2,2),VMIXI(2,2)
0099 C...SLHA Input/Output
0100       COMMON/PYLH3P/MODSEL(200),PARMIN(100),PAREXT(200),RMSOFT(0:100),
0101      &     AU(3,3),AD(3,3),AE(3,3)
0102 C...PYTHIA common blocks
0103       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0104       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0105       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0106  
0107       SAVE  /PYMSSM/,/PYSSMT/,/PYLH3P/,/PYDAT1/,/PYPARS/,/PYDAT2/
0108 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
0109       INTEGER IMODEL
0110       REAL M0,MHF,A0,MT
0111       CHARACTER*20 CHMOD(5)
0112       CHARACTER*32 FNAME
0113  
0114       COMMON /SUGNU/ XNUSUG(18)
0115       REAL XNUSUG
0116       SAVE /SUGNU/
0117  
0118       DATA CHMOD/'mSUGRA','mGMSB','non-universal SUGRA',
0119      &     'truly unified SUGRA', 'non-minimal GMSB'/
0120  
0121 C...Start by checking for incompatibilities/inconsistencies:
0122       DO 100 ICHK=2,9
0123         IF (ICHK.NE.8.AND.ICHK.NE.4.AND.IMSS(ICHK).NE.0) THEN
0124           WRITE (MSTU(11),*) '(PYSUGI:) IMSS(',ICHK,')=',IMSS(ICHK)
0125      &         ,' option not used by PYSUGI'
0126         ENDIF
0127   100 CONTINUE
0128 C...ISAJET works with REAL numbers.
0129       MZERO=REAL(RMSS(8))
0130       MHLF=REAL(RMSS(1))
0131       AZERO=REAL(RMSS(16))
0132       TANB=REAL(RMSS(5))
0133       SGNMU=REAL(RMSS(4))
0134       MTOP=REAL(PMAS(6,1))
0135       IMODEL=0
0136       IF (IMSS(1).EQ.12) THEN
0137         IMODEL=1
0138         GOTO 130
0139       ELSEIF(IMSS(1).EQ.13) THEN
0140 C...Read from isajet par file in IMSS(20)
0141         LFN=IMSS(20)
0142 C...STOP IF LFN IS ZERO (i.e. if no LFN was given).
0143         IF (LFN.EQ.0) THEN
0144           WRITE(MSTU(11),*) '(PYSUGI:) No valid unit given in IMSS(20)'
0145           GOTO 9999
0146         ENDIF
0147         WRITE(MSTU(11),*) 'READING SUSY MODEL FROM FILE...'
0148 CMrenna change to allow any susy model
0149         WRITE(MSTU(11),*) 'ENTER 1 for mSUGRA:'
0150         WRITE(MSTU(11),*) 'ENTER 2 for mGMSB:'
0151         WRITE(MSTU(11),*) 'ENTER 3 for non-universal SUGRA:'
0152         WRITE(MSTU(11),*) 'ENTER 4 for SUGRA with truly unified'//
0153      &       ' gauge couplings:'
0154         WRITE(MSTU(11),*) 'ENTER 5 for non-minimal GMSB:'
0155         READ(LFN,*) IMODEL
0156         IF (IMODEL.EQ.4) THEN
0157           IAL3UN=1
0158           IMODEL=1
0159         ENDIF
0160         IF (IMODEL.EQ.1.OR.IMODEL.EQ.3) THEN
0161           WRITE(MSTU(11),*) 'ENTER M_0, M_(1/2), A_0, tan(beta),'
0162      &         //' sgn(mu), M_t:'
0163           READ(LFN,*) M0,MHF,A0,TANB,SGNMU,MT
0164           IF (IMODEL.EQ.3) THEN
0165             IMODEL=1
0166  110        WRITE(MSTU(11),*) ' ENTER 1,...,5 for NUSUGx keyword;'
0167      &           //' 0 to continue:'
0168             WRITE(MSTU(11),*) ' NUSUG1 = GUT scale gaugino masses'
0169             WRITE(MSTU(11),*) ' NUSUG2 = GUT scale A terms'
0170             WRITE(MSTU(11),*) ' NUSUG3 = GUT scale Higgs masses'
0171             WRITE(MSTU(11),*) ' NUSUG4 = GUT scale 1st/2nd'
0172      &           //' generation masses'
0173             WRITE(MSTU(11),*)
0174      &           ' NUSUG5 = GUT scale 3rd generation masses'
0175             READ(LFN,*) INUSUG
0176             IF (INUSUG.EQ.0) THEN
0177               GOTO 120
0178             ELSEIF (INUSUG.EQ.1) THEN
0179               WRITE(MSTU(11),*) 'Enter GUT scale M_1, M_2, M_3:'
0180               READ(LFN,*) XNUSUG(1),XNUSUG(2),XNUSUG(3)
0181               IF (XNUSUG(3).LE.0.) THEN
0182                 WRITE(MSTU(11),*) ' NEGATIVE M_3 IS NOT ALLOWED'
0183                 CALL PYSTOP(109)
0184               END IF
0185             ELSEIF (INUSUG.EQ.2) THEN
0186               WRITE(MSTU(11),*) 'Enter GUT scale A_t, A_b, A_tau:'
0187               READ(LFN,*) XNUSUG(6),XNUSUG(5),XNUSUG(4)
0188             ELSEIF (INUSUG.EQ.3) THEN
0189               WRITE(MSTU(11),*) 'Enter GUT scale m_Hd, m_Hu:'
0190               READ(LFN,*) XNUSUG(7),XNUSUG(8)
0191             ELSEIF (INUSUG.EQ.4) THEN
0192               WRITE(MSTU(11),*) 'Enter GUT scale M(ul), M(dr),'
0193      &             //' M(ur), M(el), M(er):'
0194               READ(LFN,*) XNUSUG(13),XNUSUG(11),XNUSUG(12),
0195      &             XNUSUG(10),XNUSUG(9)
0196             ELSEIF (INUSUG.EQ.5) THEN
0197               WRITE(MSTU(11),*) 'Enter GUT scale M(tl), M(br), M(tr),'
0198      &              //' M(Ll), M(Lr):'
0199               READ(LFN,*) XNUSUG(18),XNUSUG(16),XNUSUG(17),
0200      &             XNUSUG(15),XNUSUG(14)
0201             ENDIF
0202             GOTO 110
0203           ENDIF
0204         ELSEIF (IMODEL.EQ.2.OR.IMODEL.EQ.5) THEN
0205           IMSS(11)=1
0206           WRITE(MSTU(11),*) 'ENTER Lambda, M_mes, N_5, tan(beta),'
0207      &         ,' sgn(mu), M_t, C_gv:'
0208           READ(LFN,*) M0,MHF,A0,TANB,SGNMU,MT,XCMGV
0209           XGMIN(7)=XCMGV
0210           XGMIN(8)=1.
0211 C...Planck scale: AMPL = 2.4 E18 GeV = {8 pi G_newton}^{1/2}
0212           AMPL=2.4D18
0213           AMGVSS=M0*MHF*XCMGV/SQRT(3D0)/AMPL
0214           IF (IMODEL.EQ.5) THEN
0215             IMODEL=2
0216             WRITE(MSTU(11),*) 'Rsl = factor multiplying gaugino'
0217      &           ,' masses at M_mes'
0218             WRITE(MSTU(11),*) 'dmH_d2, dmH_u2 = Higgs mass**2'
0219      &           ,' shifts at M_mes'
0220             WRITE(MSTU(11),*) 'd_Y = mass**2 shifts proportional to',
0221      &           ' Y at M_mes'
0222             WRITE(MSTU(11),*) 'n5_1,n5_2,n5_3 = n5 values for U(1),'
0223      &           ,'SU(2),SU(3)'
0224             WRITE(MSTU(11),*) 'ENTER Rsl, dmH_d2, dmH_u2, d_Y, n5_1,'
0225      &           ,' n5_2, n5_3'
0226             READ(LFN,*) XGMIN(8),XGMIN(9),XGMIN(10),XGMIN(11),XGMIN(12),
0227      $           XGMIN(13),XGMIN(14)
0228           ENDIF
0229         ELSE
0230           WRITE(MSTU(11),*) 'Invalid model choice.'
0231           GOTO 9999
0232         ENDIF
0233       ENDIF
0234  
0235  120  MZERO=M0
0236       MHLF=MHF
0237       AZERO=A0
0238 C     TANB=REAL(RMSS(5))
0239 C     SGNMU=REAL(RMSS(4))
0240       MTOP=MT
0241  
0242 C...Initialize MSSM parameter array
0243  130  DO 140 IPAR=1,72
0244         SUPER(IPAR)=0.0
0245  140  CONTINUE
0246 C...Call ISASUGRA
0247       CALL SUGRA(MZERO,MHLF,AZERO,TANB,SGNMU,MTOP,IMODEL)
0248 C...Check whether ISASUSY thought the model was OK.
0249       IF (NOGOOD.NE.0) THEN
0250         IF (NOGOOD.EQ.1) CALL PYERRM(26
0251      &       ,'(PYSUGI:) SUSY parameters give tachyonic particles.')
0252         IF (NOGOOD.EQ.2) CALL PYERRM(26
0253      &       ,'(PYSUGI:) SUSY parameters give no EWSB.')
0254         IF (NOGOOD.EQ.3) CALL PYERRM(26
0255      &       ,'(PYSUGI:) SUSY parameters give m(A0) < 0.')
0256         IF (NOGOOD.EQ.4) CALL PYERRM(26
0257      &       ,'(PYSUGI:) SUSY parameters give Yukawa > 100.')
0258         IF (NOGOOD.EQ.7) CALL PYERRM(26
0259      &       ,'(PYSUGI:) SUSY parameters give x_T EWSB bad.')
0260         IF (NOGOOD.EQ.8) CALL PYERRM(26
0261      &       ,'(PYSUGI:) SUSY parameters give m(h0)2 < 0.')
0262 C...Give warning, but don't stop, if LSP not ~chi_10.
0263         IF (NOGOOD.EQ.5) CALL PYERRM(16
0264      &       ,'(PYSUGI:) SUSY parameters give ~chi_10 not LSP.')
0265       ENDIF
0266 C...Warn about possible GUT scale tachyons.
0267       IF (ITACHY.NE.0) CALL PYERRM(16,
0268      &       '(PYSUGI:) Tachyonic sleptons at GUT scale.')
0269 C...Finalize spectrum (last iteration)
0270 C...(Thanks to A. Raklev for pointing this out.)
0271 C...NB: SSMSSM also calculates decays, but these are not used by Pythia.
0272       CALL SSMSSM(XISAIN(1),XISAIN(2),XISAIN(3),
0273      $ XISAIN(4),XISAIN(5),XISAIN(6),XISAIN(7),XISAIN(8),XISAIN(9),
0274      $ XISAIN(10),XISAIN(11),XISAIN(12),XISAIN(13),XISAIN(14),
0275      $ XISAIN(15),XISAIN(16),XISAIN(17),XISAIN(18),XISAIN(19),
0276      $ XISAIN(20),XISAIN(21),XISAIN(22),XISAIN(23),XISAIN(24),
0277      $ MTOP,IALLOW,1)
0278  
0279 C...M1, M2, M3.
0280       RMSS(1)=dble(GSS(7))
0281       RMSS(2)=dble(GSS(8))
0282       RMSS(3)=dble(GSS(9))
0283       RMSOFT(1)=dble(GSS(7))
0284       RMSOFT(2)=dble(GSS(8))
0285       RMSOFT(3)=dble(GSS(9))
0286 C...Mu = - Higgsino mass.
0287       RMSS(4)=-SUPER(29)
0288       RMSS(5)=TANB
0289 C...Slepton and squark masses. 2 first generations.
0290       RMSS(6)=0.5*(SUPER(18)+SUPER(20))
0291       RMSS(7)=0.5*(SUPER(19)+SUPER(21))
0292       RMSS(8)=0.25*(SUPER(2)+SUPER(4)+SUPER(6)+SUPER(8))
0293       RMSS(9)=0.25*(SUPER(3)+SUPER(5)+SUPER(7)+SUPER(9))
0294 C...Third generation.
0295       RMSS(10)=0.5*(SUPER(14)+SUPER(10))
0296       RMSS(11)=SUPER(11)
0297       RMSS(12)=SUPER(15)
0298       RMSS(13)=SUPER(22)
0299       RMSS(14)=SUPER(23)
0300 C...SLHA: store exact soft spectrum in RMSOFT
0301       RMSOFT(31)=SUPER(18)
0302       RMSOFT(32)=SUPER(20)
0303       RMSOFT(33)=SUPER(22)
0304       RMSOFT(34)=SUPER(19)
0305       RMSOFT(35)=SUPER(21)
0306       RMSOFT(36)=SUPER(23)
0307       RMSOFT(41)=0.5D0*(SUPER(2)+SUPER(4))
0308       RMSOFT(42)=0.5D0*(SUPER(6)+SUPER(8))
0309       RMSOFT(43)=0.5D0*(SUPER(10)+SUPER(14))
0310       RMSOFT(44)=SUPER(3)
0311       RMSOFT(45)=SUPER(9)
0312       RMSOFT(46)=SUPER(15)
0313       RMSOFT(47)=SUPER(5)
0314       RMSOFT(48)=SUPER(7)
0315       RMSOFT(49)=SUPER(11)
0316  
0317 C...~b, ~t, and ~tau trilinear couplings and mixing angles.
0318       RMSS(15)=SUPER(62)
0319       RMSS(16)=SUPER(60)
0320       RMSS(17)=SUPER(64)
0321       RMSS(26)=SUPER(63)
0322       RMSS(27)=SUPER(61)
0323       RMSS(28)=SUPER(65)
0324 C...SLHA trilinears
0325       DO 142 K1=1,3
0326         DO 141 K2=1,3
0327           AE(K1,K2)=0D0
0328           AU(K1,K2)=0D0
0329           AD(K1,K2)=0D0
0330  141    CONTINUE
0331  142  CONTINUE
0332       AE(3,3)=SUPER(64)
0333       AU(3,3)=SUPER(60)
0334       AD(3,3)=SUPER(62)
0335 C...Higgs mixing angle alpha (Gunion-Haber convention).
0336       RMSS(18)=-SUPER(59)
0337 C...A0 mass.
0338       RMSS(19)=SUPER(57)
0339 C...GUT scale coupling
0340       RMSS(20)=AGUTSS
0341 C...Gravitino mass (for future compatibility)
0342       RMSS(21)=MAX(RMSS(21),DBLE(SUPER(66)))
0343  
0344 C...Now we're done with RMSS. Time to fill PMAS (m > 0 required).
0345 C...Higgs sector.
0346       PMAS(PYCOMP(25),1)=ABS(SUPER(55))
0347       PMAS(PYCOMP(35),1)=ABS(SUPER(56))
0348       PMAS(PYCOMP(36),1)=ABS(SUPER(57))
0349       PMAS(PYCOMP(37),1)=ABS(SUPER(58))
0350 C...Gluino.
0351       PMAS(PYCOMP(KSUSY1+21),1)=ABS(SUPER(1))
0352 C...Squarks and Sleptons.
0353       DO 150 ILR=1,2
0354         ILRM=ILR-1
0355         PMAS(PYCOMP(ILR*KSUSY1+1),1)=ABS(SUPER(4+ILRM))
0356         PMAS(PYCOMP(ILR*KSUSY1+2),1)=ABS(SUPER(2+ILRM))
0357         PMAS(PYCOMP(ILR*KSUSY1+3),1)=ABS(SUPER(6+ILRM))
0358         PMAS(PYCOMP(ILR*KSUSY1+4),1)=ABS(SUPER(8+ILRM))
0359         PMAS(PYCOMP(ILR*KSUSY1+5),1)=ABS(SUPER(12+ILRM))
0360         PMAS(PYCOMP(ILR*KSUSY1+6),1)=ABS(SUPER(16+ILRM))
0361         PMAS(PYCOMP(ILR*KSUSY1+11),1)=ABS(SUPER(18+ILRM))
0362         PMAS(PYCOMP(ILR*KSUSY1+13),1)=ABS(SUPER(20+ILRM))
0363         PMAS(PYCOMP(ILR*KSUSY1+15),1)=ABS(SUPER(24+ILRM))
0364   150 CONTINUE
0365       PMAS(PYCOMP(KSUSY1+12),1)=ABS(SUPER(26))
0366       PMAS(PYCOMP(KSUSY1+14),1)=ABS(SUPER(27))
0367       PMAS(PYCOMP(KSUSY1+16),1)=ABS(SUPER(28))
0368 C...Neutralinos.
0369       PMAS(PYCOMP(KSUSY1+22),1)=ABS(SUPER(31))
0370       PMAS(PYCOMP(KSUSY1+23),1)=ABS(SUPER(32))
0371       PMAS(PYCOMP(KSUSY1+25),1)=ABS(SUPER(33))
0372       PMAS(PYCOMP(KSUSY1+35),1)=ABS(SUPER(34))
0373 C...Signed masses (extra minus from going to G-H convention).
0374       SMZ(1)=-SUPER(31)
0375       SMZ(2)=-SUPER(32)
0376       SMZ(3)=-SUPER(33)
0377       SMZ(4)=-SUPER(34)
0378 C...Charginos
0379       PMAS(PYCOMP(KSUSY1+24),1)=ABS(SUPER(51))
0380       PMAS(PYCOMP(KSUSY1+37),1)=ABS(SUPER(52))
0381 C...Signed masses (extra minus from going to G-H convention).
0382       SMW(1)=-SUPER(51)
0383       SMW(2)=-SUPER(52)
0384  
0385 C... Neutralino Mixing.
0386       DO 160 IN=1,4
0387         ZMIX(IN,1)= SUPER(38+4*(IN-1))
0388         ZMIX(IN,2)= SUPER(37+4*(IN-1))
0389         ZMIX(IN,3)=-SUPER(36+4*(IN-1))
0390         ZMIX(IN,4)=-SUPER(35+4*(IN-1))
0391   160 CONTINUE
0392 C...Chargino Mixing (PYTHIA same angle as HERWIG).
0393       THX=1D0
0394       THY=1D0
0395       IF (SUPER(53).GT.0) THX=-1D0
0396       IF (SUPER(54).GT.0) THY=-1D0
0397       UMIX(1,1) = -SIN(SUPER(53))
0398       UMIX(1,2) = -COS(SUPER(53))
0399       UMIX(2,1) = -THX*COS(SUPER(53))
0400       UMIX(2,2) = THX*SIN(SUPER(53))
0401       VMIX(1,1) = -SIN(SUPER(54))
0402       VMIX(1,2) = -COS(SUPER(54))
0403       VMIX(2,1) = -THY*COS(SUPER(54))
0404       VMIX(2,2) = THY*SIN(SUPER(54))
0405 C...Sfermion mixing (PYTHIA same angle as ISAJET)
0406       SFMIX(5,1)=COS(SUPER(63))
0407       SFMIX(5,2)=SIN(SUPER(63))
0408       SFMIX(5,3)=-SIN(SUPER(63))
0409       SFMIX(5,4)=COS(SUPER(63))
0410       SFMIX(6,1)=COS(SUPER(61))
0411       SFMIX(6,2)=SIN(SUPER(61))
0412       SFMIX(6,3)=-SIN(SUPER(61))
0413       SFMIX(6,4)=COS(SUPER(61))
0414       SFMIX(15,1)=COS(SUPER(65))
0415       SFMIX(15,2)=SIN(SUPER(65))
0416       SFMIX(15,3)=-SIN(SUPER(65))
0417       SFMIX(15,4)=COS(SUPER(65))
0418  
0419       IF (MSTP(122).NE.0) THEN
0420 C...Print a few lines to make the user know what's happening
0421         ISAVER=VISAJE()
0422         WRITE(MSTU(11),5000) DOC, ISAVER
0423         WRITE(MSTU(11),5100)
0424         IF (IMODEL.EQ.1) THEN
0425           WRITE(MSTU(11),5200) MZERO, MHLF, AZERO, TANB, NINT(SGNMU),
0426      &         MTOP
0427           WRITE(MSTU(11),5300)
0428         ENDIF
0429         WRITE(MSTU(11),5500) 'Pole masses'
0430         WRITE(MSTU(11),5700) (SUPER(IP),IP=2,16,2),(SUPER(IP),IP=3,17,2)
0431         WRITE(MSTU(11),5800) (SUPER(IP),IP=18,24,2),(SUPER(IP),IP=26,28)
0432      &       ,(SUPER(IP),IP=19,25,2)
0433         WRITE(MSTU(11),5900) SUPER(1),(SMZ(IP),IP=1,4), (SMW(IP)
0434      &       ,IP=1,2)
0435         WRITE(MSTU(11),5400)
0436         WRITE(MSTU(11),6000) (SUPER(IP),IP=55,58)
0437         WRITE(MSTU(11),5400)
0438         WRITE(MSTU(11),5500) 'EW scale mixing structure'
0439         WRITE(MSTU(11),6100) ((ZMIX(I,J), J=1,4),I=1,4)
0440         WRITE(MSTU(11),6200) (UMIX(1,J), J=1,2),(VMIX(1,J),J=1,2)
0441      &       ,(UMIX(2,J), J=1,2),(VMIX(2,J),J=1,2)
0442         WRITE(MSTU(11),6300) (SFMIX(5,J), J=1,2),(SFMIX(6,J),J=1,2)
0443      &       ,(SFMIX(15,J), J=1,2),(SFMIX(5,J),J=3,4),(SFMIX(6,J), J=3,4
0444      &       ),(SFMIX(15,J),J=3,4)
0445         WRITE(MSTU(11),5400)
0446         WRITE(MSTU(11),6450) RMSS(18)
0447         WRITE(MSTU(11),5400)
0448         WRITE(MSTU(11),5500) 'Couplings'
0449         WRITE(MSTU(11),6400) RMSS(15),RMSS(16),RMSS(17),RMSS(20)
0450         WRITE(MSTU(11),5400)
0451       ENDIF
0452  
0453 C...Call FeynHiggs to improve Higgs sector if requested
0454       IF (IMSS(4).EQ.3) THEN
0455         IF (MSTP(122).NE.0) WRITE(MSTU(11),'(1x,"*"/1x,"*",A)')
0456      &       ' (PYSUGI:) Now calling FeynHiggs.'
0457         CALL PYFEYN(IERR)
0458         IF (IERR.EQ.0) THEN
0459           IMSS(4)=2
0460           IF (MSTP(122).NE.0) THEN
0461             WRITE(MSTU(11),5400)
0462             WRITE(MSTU(11),5500)
0463      &           'Corrected Higgs masses and mixing'
0464             WRITE(MSTU(11),6000) PMAS(25,1),PMAS(35,1),PMAS(36,1),
0465      &           PMAS(37,1)
0466             WRITE(MSTU(11),6450) RMSS(18)
0467             WRITE(MSTU(11),5400)
0468           ENDIF
0469         ENDIF
0470       ENDIF
0471  
0472       IF (MSTP(122).NE.0) WRITE(MSTU(11),6500)
0473  
0474 C...Fix the higgs sector (in PYMSIN) using the masses and mixing angle
0475 C...output by ISASUSY.
0476       IMSS(4)=MAX(2,IMSS(4))
0477  
0478  5000 FORMAT(1x,19('*'),1x,'PYSUGI v1.52: PYTHIA/ISASUSY '
0479      &     ,'INTERFACE',1x,19('*')/1x,'*',3x,'PYSUGI: Last Change',1x,A
0480      &     ,1x,'-',1x,'P. Skands / S. Mrenna'/1x,'*',2x,A/1x,'*')
0481  5100 FORMAT(1x,'*',1x,'ISASUSY Input:'/1x,'*',1x,'----------------')
0482  5200 FORMAT(1x,'*',1x,3x,'M_0',6x,'M_1/2',5x,'A_0',3x,'Tan(beta)',
0483      &     3x,'Sgn(mu)',3x,'M_t'/1x,'*',1x,4(F8.2,1x),I8,2x,F8.2)
0484  5300 FORMAT(1x,'*'/1x,'*',1x,'ISASUSY Output:'/1x,'*',1x
0485      &     ,'----------------')
0486  5400 FORMAT(1x,'*',1x,A)
0487  5500 FORMAT(1x,'*',1x,A,':')
0488  5600 FORMAT(1x,'*',2x,2x,'M_GUT',2x,2x,'g_GUT',2x,1x,'alpha_GUT'/
0489      &       1x,'*',2x,1P,2(1x,E8.2),2x,E8.2)
0490  5700 FORMAT(1x,'*',4x,4x,'~u',2x,1x,4x,'~d',2x,1x,4x,'~s',2x,1x,
0491      &     4x,'~c',2x,1x,4x,'~b',2x,1x,2x,'~b(12)',1x,4x,'~t',2x,1x, 2x,
0492      &     '~t(12)'/1x,'*',2x,'L',1x,8(F8.2,1x)/1x,'*',2x,'R',1x,8(F8.2
0493      &     ,1x))
0494  5800 FORMAT(1x,'*'/1x,'*',4x,4x,'~e',2x,1x,3x,'~mu',2x,1x,3x,'~tau',1x
0495      &     ,1x,'~tau(12)',1x,2x,'~nu_e',1x,1x,1x,'~nu_mu',1x,1x,1x
0496      &     ,'~nu_tau'/1x,'*',2x,'L',1x,7(F8.2,1x)/1x,'*',2x,'R',1x,4(F8
0497      &     .2,1x))
0498  5900 FORMAT(1x,'*'/1x,'*',4x,4x,'~g',2x,1x,1x,'~chi_10',1x,1x,'~chi_20'
0499      &     ,1x,1x,'~chi_30',1x,1x,'~chi_40',1x,1x,'~chi_1+',1x
0500      &     ,1x,'~chi_2+'/1x,'*',3x,1x,7(F8.2,1x))
0501  6000 FORMAT(1x,'*',4x,4x,'h0',2x,1x,4x,'H0',2x,1x,4x,'A0',2x
0502      &     ,1x,4x,'H+'/1x,'*',3x,1x,5(F8.2,1x))
0503  6050 FORMAT(1x,'*'/1x,'*',4x,4x,'h0',2x,1x,4x,'H0',2x,1x,4x,'A0',2x
0504      &     ,1x,4x,'H+'/1x,'*',3x,1x,5(F8.2,1x),3x,'(Before FeynHiggs)')
0505  6100 FORMAT(1x,'*',11x,'|',3x,'~B',3x,'|',2x,'~W_3',2x,'|',2x
0506      &     ,'~H_1',2x,'|',2x,'~H_2',2x,'|'/1x,'*',3x,'~chi_10',1x,4('|'
0507      &     ,1x,F6.3,1x),'|'/1x,'*',3x,'~chi_20',1x,4('|'
0508      &     ,1x,F6.3,1x),'|'/1x,'*',3x,'~chi_30',1x,4('|'
0509      &     ,1x,F6.3,1x),'|'/1x,'*',3x,'~chi_40',1x,4('|'
0510      &     ,1x,F6.3,1x),'|')
0511  6200 FORMAT(1x,'*'/1x,'*',6x,'L',4x,'|',3x,'~W',3x,'|',3x,'~H',3x,'|'
0512      &     ,12x,'R',4x,'|',3x,'~W',3x,'|',3x,'~H',3x,'|'/1x,'*',3x
0513      &     ,'~chi_1+',1x,2('|',1x,F6.3,1x),'|',9x,'~chi_1+',1x,2('|',1x
0514      &     ,F6.3,1x),'|'/1x,'*',3x,'~chi_2+',1x,2('|',1x,F6.3,1x),'|',9x
0515      &     ,'~chi_2+',1x,2('|',1x,F6.3,1x),'|')
0516  6300 FORMAT(1x,'*'/1x,'*',8x,'|',2x,'~b_L',2x,'|',2x,'~b_R',2x,'|',8x
0517      &     ,'|',2x,'~t_L',2x,'|',2x,'~t_R',2x,'|',10x
0518      &     ,'|',1x,'~tau_L',1x,'|',1x,'~tau_R',1x,'|'/
0519      &     1x,'*',3x,'~b_1',1x,2('|',1x,F6.3,1x),'|',3x,'~t_1',1x,2('|'
0520      &     ,1x,F6.3,1x),'|',3x,'~tau_1',1x,2('|',1x,F6.3,1x),'|'/
0521      &     1x,'*',3x,'~b_2',1x,2('|',1x,F6.3,1x),'|',3x,'~t_2',1x,2('|'
0522      &     ,1x,F6.3,1x),'|',3x,'~tau_2',1x,2('|',1x,F6.3,1x),'|')
0523  6400 FORMAT(1x,'*',3x,'A_b = ',F8.2,4x,'A_t = ',F8.2,4x,'A_tau = ',F8.2
0524      &     ,4x,'Alpha_GUT = ',F8.2)
0525  6450 FORMAT(1x,'*',3x,'Alpha_Higgs = ',F8.4)
0526  6500 FORMAT(1x,32('*'),1x,'END OF PYSUGI',1x,31('*'))
0527  
0528  9999 RETURN
0529       END
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!