Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
pysugi.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pysugi.f
1 
2 C*********************************************************************
3 
4 C...PYSUGI
5 C...Interface to ISASUSY version 7.71.
6 C...Warning: this interface should not be used with earlier versions
7 C...of ISASUSY, since common block incompatibilities may then arise.
8 C...Calls SUGRA (in ISAJET) to perform RGE evolution.
9 C...Then converts to Gunion-Haber conventions.
10 
11  SUBROUTINE pysugi
12  IMPLICIT DOUBLE PRECISION(a-h, o-z)
13 
14  INTEGER pyk,pychge,pycomp
15  parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
16  &kexcit=4000000,kdimen=5000000)
17 
18 C...Date of Change
19  CHARACTER doc*11
20  parameter(doc='01 May 2006')
21 
22 C...ISASUGRA Input:
23  REAL mzero,mhlf,azero,tanb,sgnmu,mtop
24 C...XISAIN contains the MSSMi inputs in natural order.
25  COMMON /sugxin/ xisain(24),xsugin(7),xgmin(14),xnrin(4),
26  $xamin(7)
27  REAL xisain,xsugin,xgmin,xnrin,xamin
28  SAVE /sugxin/
29 C...ISASUGRA Output
30  CHARACTER*40 isaver,visaje
31  REAL super
32  COMMON /sspar/ super(72)
33  COMMON /sugmg/ mss(32),gss(31),mgutss,ggutss,agutss,ftgut,
34  $fbgut,ftagut,fngut
35  REAL mss,gss,mgutss,ggutss,agutss,ftgut,fbgut,ftagut,fngut
36  COMMON /sugpas/ xtanb,msusy,amt,mgut,mu,g2,gp,v,vp,xw,
37  $a1mz,a2mz,asmz,ftamz,fbmz,b,sin2b,ftmt,g3mt,vev,higfrz,
38  $fnmz,amnrmj,nogood,ial3un,itachy,mhpneg,asm3,
39  $vumt,vdmt,asmtp,asmss,m3q
40  REAL xtanb,msusy,amt,mgut,mu,g2,gp,v,vp,xw,
41  $a1mz,a2mz,asmz,ftamz,fbmz,b,sin2b,ftmt,g3mt,vev,higfrz,
42  $fnmz,amnrmj,asm3,vumt,vdmt,asmtp,asmss,m3q
43  INTEGER nogood,ial3un,itachy,mhpneg
44  INTEGER iallow
45  SAVE /sugmg/,/sspar/
46 C SUPER: Filled by ISASUGRA.
47 C SUPER(1) = mass of ~g
48 C SUPER(2:17) = mass of ~u_L,~u_R,~d_L,~d_R,~s_L,~s_R,~c_L,~c_R,~b_L
49 C ,~b_R,~b_1,~b_2,~t_L,~t_R,~t_1,~t_2
50 C SUPER(18:25) = mass of ~e_L,~e_R,~mu_L,~mu_R,~tau_L,~tau_R,~tau_1
51 C ,~tau_2
52 C SUPER(26:28) = mass of ~nu_e,~nu_mu,~nu_tau
53 C SUPER(29) = Higgsino mass = - mu
54 C SUPER(30) = ratio v2/v1 of vev's
55 C SUPER(31:34) = Signed neutralino masses
56 C SUPER(35:50) = Neutralino mixing matrix
57 C SUPER(51:52) = Signed chargino masses
58 C SUPER(53:54) = Chargino left, right mixing angles
59 C SUPER(55:58) = mass of h0, H0, A0, H+
60 C SUPER(59) = Higgs mixing angle alpha
61 C SUPER(60:65) = A_t, theta_t, A_b, theta_b, A_tau, theta_tau
62 C SUPER(66) = Gravitino mass
63 C SUPER(67:69) = Top,Bottom, and Tau masses at MSUSY (not used)
64 C SUPER(70) = b-Yukawa at mA scale (not used)
65 C SUPER(71:72) = H_u, H_d vev's at MSUSY (not used)
66 C GSS: Filled by ISASUGRA
67 C GSS( 1) = g_1 GSS( 2) = g_2 GSS( 3) = g_3
68 C GSS( 4) = y_tau GSS( 5) = y_b GSS( 6) = y_t
69 C GSS( 7) = M_1 GSS( 8) = M_2 GSS( 9) = M_3
70 C GSS(10) = A_tau GSS(11) = A_b GSS(12) = A_t
71 C GSS(13) = M_h12 GSS(14) = M_h22 GSS(15) = M_er2
72 C GSS(16) = M_el2 GSS(17) = M_dnr2 GSS(18) = M_upr2
73 C GSS(19) = M_upl2 GSS(20) = M_taur2 GSS(21) = M_taul2
74 C GSS(22) = M_btr2 GSS(23) = M_tpr2 GSS(24) = M_tpl2
75 C GSS(25) = mu GSS(26) = B GSS(27) = Y_N
76 C GSS(28) = M_nr GSS(29) = A_n GSS(30) = log(vdq)
77 C GSS(31) = log(vuq)
78 C MSS: Filled by ISASUGRA
79 C MSS( 1) = glss MSS( 2) = upl MSS( 3) = upr
80 C MSS( 4) = dnl MSS( 5) = dnr MSS( 6) = stl
81 C MSS( 7) = str MSS( 8) = chl MSS( 9) = chr
82 C MSS(10) = b1 MSS(11) = b2 MSS(12) = t1
83 C MSS(13) = t2 MSS(14) = nuel MSS(15) = numl
84 C MSS(16) = nutl MSS(17) = el- MSS(18) = er-
85 C MSS(19) = mul- MSS(20) = mur- MSS(21) = tau1
86 C MSS(22) = tau2 MSS(23) = z1ss MSS(24) = z2ss
87 C MSS(25) = z3ss MSS(26) = z4ss MSS(27) = w1ss
88 C MSS(28) = w2ss MSS(29) = hl0 MSS(30) = hh0
89 C MSS(31) = ha0 MSS(32) = h+
90 C Unification, filled by ISASUGRA if applicable.
91 C MGUTSS = M_GUT GGUTSS = g_GUT AGUTSS = alpha_GUTC
92 
93 C...SPYTHIA Input/Output
94  INTEGER imss
95  DOUBLE PRECISION rmss
96  common/pymssm/imss(0:99),rmss(0:99)
97  common/pyssmt/zmix(4,4),umix(2,2),vmix(2,2),smz(4),smw(2),
98  &sfmix(16,4),zmixi(4,4),umixi(2,2),vmixi(2,2)
99 C...SLHA Input/Output
100  common/pylh3p/modsel(200),parmin(100),parext(200),rmsoft(0:100),
101  & au(3,3),ad(3,3),ae(3,3)
102 C...PYTHIA common blocks
103  common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
104  common/pypars/mstp(200),parp(200),msti(200),pari(200)
105  common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
106 
107  SAVE /pymssm/,/pyssmt/,/pylh3p/,/pydat1/,/pypars/,/pydat2/
108 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
109  INTEGER imodel
110  REAL m0,mhf,a0,mt
111  CHARACTER*20 chmod(5)
112  CHARACTER*32 fname
113 
114  COMMON /sugnu/ xnusug(18)
115  REAL xnusug
116  SAVE /sugnu/
117 
118  DATA chmod/'mSUGRA','mGMSB','non-universal SUGRA',
119  & 'truly unified SUGRA', 'non-minimal GMSB'/
120 
121 C...Start by checking for incompatibilities/inconsistencies:
122  DO 100 ichk=2,9
123  IF (ichk.NE.8.AND.ichk.NE.4.AND.imss(ichk).NE.0) THEN
124  WRITE (mstu(11),*) '(PYSUGI:) IMSS(',ichk,')=',imss(ichk)
125  & ,' option not used by PYSUGI'
126  ENDIF
127  100 CONTINUE
128 C...ISAJET works with REAL numbers.
129  mzero=REAL(rmss(8))
130  mhlf=REAL(rmss(1))
131  azero=REAL(rmss(16))
132  tanb=REAL(rmss(5))
133  sgnmu=REAL(rmss(4))
134  mtop=REAL(pmas(6,1))
135  imodel=0
136  IF (imss(1).EQ.12) THEN
137  imodel=1
138  goto 130
139  ELSEIF(imss(1).EQ.13) THEN
140 C...Read from isajet par file in IMSS(20)
141  lfn=imss(20)
142 C...STOP IF LFN IS ZERO (i.e. if no LFN was given).
143  IF (lfn.EQ.0) THEN
144  WRITE(mstu(11),*) '(PYSUGI:) No valid unit given in IMSS(20)'
145  goto 9999
146  ENDIF
147  WRITE(mstu(11),*) 'READING SUSY MODEL FROM FILE...'
148 CMrenna change to allow any susy model
149  WRITE(mstu(11),*) 'ENTER 1 for mSUGRA:'
150  WRITE(mstu(11),*) 'ENTER 2 for mGMSB:'
151  WRITE(mstu(11),*) 'ENTER 3 for non-universal SUGRA:'
152  WRITE(mstu(11),*) 'ENTER 4 for SUGRA with truly unified'//
153  & ' gauge couplings:'
154  WRITE(mstu(11),*) 'ENTER 5 for non-minimal GMSB:'
155  READ(lfn,*) imodel
156  IF (imodel.EQ.4) THEN
157  ial3un=1
158  imodel=1
159  ENDIF
160  IF (imodel.EQ.1.OR.imodel.EQ.3) THEN
161  WRITE(mstu(11),*) 'ENTER M_0, M_(1/2), A_0, tan(beta),'
162  & //' sgn(mu), M_t:'
163  READ(lfn,*) m0,mhf,a0,tanb,sgnmu,mt
164  IF (imodel.EQ.3) THEN
165  imodel=1
166  110 WRITE(mstu(11),*) ' ENTER 1,...,5 for NUSUGx keyword;'
167  & //' 0 to continue:'
168  WRITE(mstu(11),*) ' NUSUG1 = GUT scale gaugino masses'
169  WRITE(mstu(11),*) ' NUSUG2 = GUT scale A terms'
170  WRITE(mstu(11),*) ' NUSUG3 = GUT scale Higgs masses'
171  WRITE(mstu(11),*) ' NUSUG4 = GUT scale 1st/2nd'
172  & //' generation masses'
173  WRITE(mstu(11),*)
174  & ' NUSUG5 = GUT scale 3rd generation masses'
175  READ(lfn,*) inusug
176  IF (inusug.EQ.0) THEN
177  goto 120
178  ELSEIF (inusug.EQ.1) THEN
179  WRITE(mstu(11),*) 'Enter GUT scale M_1, M_2, M_3:'
180  READ(lfn,*) xnusug(1),xnusug(2),xnusug(3)
181  IF (xnusug(3).LE.0.) THEN
182  WRITE(mstu(11),*) ' NEGATIVE M_3 IS NOT ALLOWED'
183  CALL pystop(109)
184  END IF
185  ELSEIF (inusug.EQ.2) THEN
186  WRITE(mstu(11),*) 'Enter GUT scale A_t, A_b, A_tau:'
187  READ(lfn,*) xnusug(6),xnusug(5),xnusug(4)
188  ELSEIF (inusug.EQ.3) THEN
189  WRITE(mstu(11),*) 'Enter GUT scale m_Hd, m_Hu:'
190  READ(lfn,*) xnusug(7),xnusug(8)
191  ELSEIF (inusug.EQ.4) THEN
192  WRITE(mstu(11),*) 'Enter GUT scale M(ul), M(dr),'
193  & //' M(ur), M(el), M(er):'
194  READ(lfn,*) xnusug(13),xnusug(11),xnusug(12),
195  & xnusug(10),xnusug(9)
196  ELSEIF (inusug.EQ.5) THEN
197  WRITE(mstu(11),*) 'Enter GUT scale M(tl), M(br), M(tr),'
198  & //' M(Ll), M(Lr):'
199  READ(lfn,*) xnusug(18),xnusug(16),xnusug(17),
200  & xnusug(15),xnusug(14)
201  ENDIF
202  goto 110
203  ENDIF
204  ELSEIF (imodel.EQ.2.OR.imodel.EQ.5) THEN
205  imss(11)=1
206  WRITE(mstu(11),*) 'ENTER Lambda, M_mes, N_5, tan(beta),'
207  & ,' sgn(mu), M_t, C_gv:'
208  READ(lfn,*) m0,mhf,a0,tanb,sgnmu,mt,xcmgv
209  xgmin(7)=xcmgv
210  xgmin(8)=1.
211 C...Planck scale: AMPL = 2.4 E18 GeV = {8 pi G_newton}^{1/2}
212  ampl=2.4d18
213  amgvss=m0*mhf*xcmgv/sqrt(3d0)/ampl
214  IF (imodel.EQ.5) THEN
215  imodel=2
216  WRITE(mstu(11),*) 'Rsl = factor multiplying gaugino'
217  & ,' masses at M_mes'
218  WRITE(mstu(11),*) 'dmH_d2, dmH_u2 = Higgs mass**2'
219  & ,' shifts at M_mes'
220  WRITE(mstu(11),*) 'd_Y = mass**2 shifts proportional to',
221  & ' Y at M_mes'
222  WRITE(mstu(11),*) 'n5_1,n5_2,n5_3 = n5 values for U(1),'
223  & ,'SU(2),SU(3)'
224  WRITE(mstu(11),*) 'ENTER Rsl, dmH_d2, dmH_u2, d_Y, n5_1,'
225  & ,' n5_2, n5_3'
226  READ(lfn,*) xgmin(8),xgmin(9),xgmin(10),xgmin(11),xgmin(12),
227  $ xgmin(13),xgmin(14)
228  ENDIF
229  ELSE
230  WRITE(mstu(11),*) 'Invalid model choice.'
231  goto 9999
232  ENDIF
233  ENDIF
234 
235  120 mzero=m0
236  mhlf=mhf
237  azero=a0
238 C TANB=REAL(RMSS(5))
239 C SGNMU=REAL(RMSS(4))
240  mtop=mt
241 
242 C...Initialize MSSM parameter array
243  130 DO 140 ipar=1,72
244  super(ipar)=0.0
245  140 CONTINUE
246 C...Call ISASUGRA
247  CALL sugra(mzero,mhlf,azero,tanb,sgnmu,mtop,imodel)
248 C...Check whether ISASUSY thought the model was OK.
249  IF (nogood.NE.0) THEN
250  IF (nogood.EQ.1) CALL pyerrm(26
251  & ,'(PYSUGI:) SUSY parameters give tachyonic particles.')
252  IF (nogood.EQ.2) CALL pyerrm(26
253  & ,'(PYSUGI:) SUSY parameters give no EWSB.')
254  IF (nogood.EQ.3) CALL pyerrm(26
255  & ,'(PYSUGI:) SUSY parameters give m(A0) < 0.')
256  IF (nogood.EQ.4) CALL pyerrm(26
257  & ,'(PYSUGI:) SUSY parameters give Yukawa > 100.')
258  IF (nogood.EQ.7) CALL pyerrm(26
259  & ,'(PYSUGI:) SUSY parameters give x_T EWSB bad.')
260  IF (nogood.EQ.8) CALL pyerrm(26
261  & ,'(PYSUGI:) SUSY parameters give m(h0)2 < 0.')
262 C...Give warning, but don't stop, if LSP not ~chi_10.
263  IF (nogood.EQ.5) CALL pyerrm(16
264  & ,'(PYSUGI:) SUSY parameters give ~chi_10 not LSP.')
265  ENDIF
266 C...Warn about possible GUT scale tachyons.
267  IF (itachy.NE.0) CALL pyerrm(16,
268  & '(PYSUGI:) Tachyonic sleptons at GUT scale.')
269 C...Finalize spectrum (last iteration)
270 C...(Thanks to A. Raklev for pointing this out.)
271 C...NB: SSMSSM also calculates decays, but these are not used by Pythia.
272  CALL ssmssm(xisain(1),xisain(2),xisain(3),
273  $ xisain(4),xisain(5),xisain(6),xisain(7),xisain(8),xisain(9),
274  $ xisain(10),xisain(11),xisain(12),xisain(13),xisain(14),
275  $ xisain(15),xisain(16),xisain(17),xisain(18),xisain(19),
276  $ xisain(20),xisain(21),xisain(22),xisain(23),xisain(24),
277  $ mtop,iallow,1)
278 
279 C...M1, M2, M3.
280  rmss(1)=dble(gss(7))
281  rmss(2)=dble(gss(8))
282  rmss(3)=dble(gss(9))
283  rmsoft(1)=dble(gss(7))
284  rmsoft(2)=dble(gss(8))
285  rmsoft(3)=dble(gss(9))
286 C...Mu = - Higgsino mass.
287  rmss(4)=-super(29)
288  rmss(5)=tanb
289 C...Slepton and squark masses. 2 first generations.
290  rmss(6)=0.5*(super(18)+super(20))
291  rmss(7)=0.5*(super(19)+super(21))
292  rmss(8)=0.25*(super(2)+super(4)+super(6)+super(8))
293  rmss(9)=0.25*(super(3)+super(5)+super(7)+super(9))
294 C...Third generation.
295  rmss(10)=0.5*(super(14)+super(10))
296  rmss(11)=super(11)
297  rmss(12)=super(15)
298  rmss(13)=super(22)
299  rmss(14)=super(23)
300 C...SLHA: store exact soft spectrum in RMSOFT
301  rmsoft(31)=super(18)
302  rmsoft(32)=super(20)
303  rmsoft(33)=super(22)
304  rmsoft(34)=super(19)
305  rmsoft(35)=super(21)
306  rmsoft(36)=super(23)
307  rmsoft(41)=0.5d0*(super(2)+super(4))
308  rmsoft(42)=0.5d0*(super(6)+super(8))
309  rmsoft(43)=0.5d0*(super(10)+super(14))
310  rmsoft(44)=super(3)
311  rmsoft(45)=super(9)
312  rmsoft(46)=super(15)
313  rmsoft(47)=super(5)
314  rmsoft(48)=super(7)
315  rmsoft(49)=super(11)
316 
317 C...~b, ~t, and ~tau trilinear couplings and mixing angles.
318  rmss(15)=super(62)
319  rmss(16)=super(60)
320  rmss(17)=super(64)
321  rmss(26)=super(63)
322  rmss(27)=super(61)
323  rmss(28)=super(65)
324 C...SLHA trilinears
325  DO 142 k1=1,3
326  DO 141 k2=1,3
327  ae(k1,k2)=0d0
328  au(k1,k2)=0d0
329  ad(k1,k2)=0d0
330  141 CONTINUE
331  142 CONTINUE
332  ae(3,3)=super(64)
333  au(3,3)=super(60)
334  ad(3,3)=super(62)
335 C...Higgs mixing angle alpha (Gunion-Haber convention).
336  rmss(18)=-super(59)
337 C...A0 mass.
338  rmss(19)=super(57)
339 C...GUT scale coupling
340  rmss(20)=agutss
341 C...Gravitino mass (for future compatibility)
342  rmss(21)=max(rmss(21),dble(super(66)))
343 
344 C...Now we're done with RMSS. Time to fill PMAS (m > 0 required).
345 C...Higgs sector.
346  pmas(pycomp(25),1)=abs(super(55))
347  pmas(pycomp(35),1)=abs(super(56))
348  pmas(pycomp(36),1)=abs(super(57))
349  pmas(pycomp(37),1)=abs(super(58))
350 C...Gluino.
351  pmas(pycomp(ksusy1+21),1)=abs(super(1))
352 C...Squarks and Sleptons.
353  DO 150 ilr=1,2
354  ilrm=ilr-1
355  pmas(pycomp(ilr*ksusy1+1),1)=abs(super(4+ilrm))
356  pmas(pycomp(ilr*ksusy1+2),1)=abs(super(2+ilrm))
357  pmas(pycomp(ilr*ksusy1+3),1)=abs(super(6+ilrm))
358  pmas(pycomp(ilr*ksusy1+4),1)=abs(super(8+ilrm))
359  pmas(pycomp(ilr*ksusy1+5),1)=abs(super(12+ilrm))
360  pmas(pycomp(ilr*ksusy1+6),1)=abs(super(16+ilrm))
361  pmas(pycomp(ilr*ksusy1+11),1)=abs(super(18+ilrm))
362  pmas(pycomp(ilr*ksusy1+13),1)=abs(super(20+ilrm))
363  pmas(pycomp(ilr*ksusy1+15),1)=abs(super(24+ilrm))
364  150 CONTINUE
365  pmas(pycomp(ksusy1+12),1)=abs(super(26))
366  pmas(pycomp(ksusy1+14),1)=abs(super(27))
367  pmas(pycomp(ksusy1+16),1)=abs(super(28))
368 C...Neutralinos.
369  pmas(pycomp(ksusy1+22),1)=abs(super(31))
370  pmas(pycomp(ksusy1+23),1)=abs(super(32))
371  pmas(pycomp(ksusy1+25),1)=abs(super(33))
372  pmas(pycomp(ksusy1+35),1)=abs(super(34))
373 C...Signed masses (extra minus from going to G-H convention).
374  smz(1)=-super(31)
375  smz(2)=-super(32)
376  smz(3)=-super(33)
377  smz(4)=-super(34)
378 C...Charginos
379  pmas(pycomp(ksusy1+24),1)=abs(super(51))
380  pmas(pycomp(ksusy1+37),1)=abs(super(52))
381 C...Signed masses (extra minus from going to G-H convention).
382  smw(1)=-super(51)
383  smw(2)=-super(52)
384 
385 C... Neutralino Mixing.
386  DO 160 in=1,4
387  zmix(in,1)= super(38+4*(in-1))
388  zmix(in,2)= super(37+4*(in-1))
389  zmix(in,3)=-super(36+4*(in-1))
390  zmix(in,4)=-super(35+4*(in-1))
391  160 CONTINUE
392 C...Chargino Mixing (PYTHIA same angle as HERWIG).
393  thx=1d0
394  thy=1d0
395  IF (super(53).GT.0) thx=-1d0
396  IF (super(54).GT.0) thy=-1d0
397  umix(1,1) = -sin(super(53))
398  umix(1,2) = -cos(super(53))
399  umix(2,1) = -thx*cos(super(53))
400  umix(2,2) = thx*sin(super(53))
401  vmix(1,1) = -sin(super(54))
402  vmix(1,2) = -cos(super(54))
403  vmix(2,1) = -thy*cos(super(54))
404  vmix(2,2) = thy*sin(super(54))
405 C...Sfermion mixing (PYTHIA same angle as ISAJET)
406  sfmix(5,1)=cos(super(63))
407  sfmix(5,2)=sin(super(63))
408  sfmix(5,3)=-sin(super(63))
409  sfmix(5,4)=cos(super(63))
410  sfmix(6,1)=cos(super(61))
411  sfmix(6,2)=sin(super(61))
412  sfmix(6,3)=-sin(super(61))
413  sfmix(6,4)=cos(super(61))
414  sfmix(15,1)=cos(super(65))
415  sfmix(15,2)=sin(super(65))
416  sfmix(15,3)=-sin(super(65))
417  sfmix(15,4)=cos(super(65))
418 
419  IF (mstp(122).NE.0) THEN
420 C...Print a few lines to make the user know what's happening
421  isaver=visaje()
422  WRITE(mstu(11),5000) doc, isaver
423  WRITE(mstu(11),5100)
424  IF (imodel.EQ.1) THEN
425  WRITE(mstu(11),5200) mzero, mhlf, azero, tanb, nint(sgnmu),
426  & mtop
427  WRITE(mstu(11),5300)
428  ENDIF
429  WRITE(mstu(11),5500) 'Pole masses'
430  WRITE(mstu(11),5700) (super(ip),ip=2,16,2),(super(ip),ip=3,17,2)
431  WRITE(mstu(11),5800) (super(ip),ip=18,24,2),(super(ip),ip=26,28)
432  & ,(super(ip),ip=19,25,2)
433  WRITE(mstu(11),5900) super(1),(smz(ip),ip=1,4), (smw(ip)
434  & ,ip=1,2)
435  WRITE(mstu(11),5400)
436  WRITE(mstu(11),6000) (super(ip),ip=55,58)
437  WRITE(mstu(11),5400)
438  WRITE(mstu(11),5500) 'EW scale mixing structure'
439  WRITE(mstu(11),6100) ((zmix(i,j), j=1,4),i=1,4)
440  WRITE(mstu(11),6200) (umix(1,j), j=1,2),(vmix(1,j),j=1,2)
441  & ,(umix(2,j), j=1,2),(vmix(2,j),j=1,2)
442  WRITE(mstu(11),6300) (sfmix(5,j), j=1,2),(sfmix(6,j),j=1,2)
443  & ,(sfmix(15,j), j=1,2),(sfmix(5,j),j=3,4),(sfmix(6,j), j=3,4
444  & ),(sfmix(15,j),j=3,4)
445  WRITE(mstu(11),5400)
446  WRITE(mstu(11),6450) rmss(18)
447  WRITE(mstu(11),5400)
448  WRITE(mstu(11),5500) 'Couplings'
449  WRITE(mstu(11),6400) rmss(15),rmss(16),rmss(17),rmss(20)
450  WRITE(mstu(11),5400)
451  ENDIF
452 
453 C...Call FeynHiggs to improve Higgs sector if requested
454  IF (imss(4).EQ.3) THEN
455  IF (mstp(122).NE.0) WRITE(mstu(11),'(1x,"*"/1x,"*",A)')
456  & ' (PYSUGI:) Now calling FeynHiggs.'
457  CALL pyfeyn(ierr)
458  IF (ierr.EQ.0) THEN
459  imss(4)=2
460  IF (mstp(122).NE.0) THEN
461  WRITE(mstu(11),5400)
462  WRITE(mstu(11),5500)
463  & 'Corrected Higgs masses and mixing'
464  WRITE(mstu(11),6000) pmas(25,1),pmas(35,1),pmas(36,1),
465  & pmas(37,1)
466  WRITE(mstu(11),6450) rmss(18)
467  WRITE(mstu(11),5400)
468  ENDIF
469  ENDIF
470  ENDIF
471 
472  IF (mstp(122).NE.0) WRITE(mstu(11),6500)
473 
474 C...Fix the higgs sector (in PYMSIN) using the masses and mixing angle
475 C...output by ISASUSY.
476  imss(4)=max(2,imss(4))
477 
478  5000 FORMAT(1x,19('*'),1x,'PYSUGI v1.52: PYTHIA/ISASUSY '
479  & ,'INTERFACE',1x,19('*')/1x,'*',3x,'PYSUGI: Last Change',1x,a
480  & ,1x,'-',1x,'P. Skands / S. Mrenna'/1x,'*',2x,a/1x,'*')
481  5100 FORMAT(1x,'*',1x,'ISASUSY Input:'/1x,'*',1x,'----------------')
482  5200 FORMAT(1x,'*',1x,3x,'M_0',6x,'M_1/2',5x,'A_0',3x,'Tan(beta)',
483  & 3x,'Sgn(mu)',3x,'M_t'/1x,'*',1x,4(f8.2,1x),i8,2x,f8.2)
484  5300 FORMAT(1x,'*'/1x,'*',1x,'ISASUSY Output:'/1x,'*',1x
485  & ,'----------------')
486  5400 FORMAT(1x,'*',1x,a)
487  5500 FORMAT(1x,'*',1x,a,':')
488  5600 FORMAT(1x,'*',2x,2x,'M_GUT',2x,2x,'g_GUT',2x,1x,'alpha_GUT'/
489  & 1x,'*',2x,1p,2(1x,e8.2),2x,e8.2)
490  5700 FORMAT(1x,'*',4x,4x,'~u',2x,1x,4x,'~d',2x,1x,4x,'~s',2x,1x,
491  & 4x,'~c',2x,1x,4x,'~b',2x,1x,2x,'~b(12)',1x,4x,'~t',2x,1x, 2x,
492  & '~t(12)'/1x,'*',2x,'L',1x,8(f8.2,1x)/1x,'*',2x,'R',1x,8(f8.2
493  & ,1x))
494  5800 FORMAT(1x,'*'/1x,'*',4x,4x,'~e',2x,1x,3x,'~mu',2x,1x,3x,'~tau',1x
495  & ,1x,'~tau(12)',1x,2x,'~nu_e',1x,1x,1x,'~nu_mu',1x,1x,1x
496  & ,'~nu_tau'/1x,'*',2x,'L',1x,7(f8.2,1x)/1x,'*',2x,'R',1x,4(f8
497  & .2,1x))
498  5900 FORMAT(1x,'*'/1x,'*',4x,4x,'~g',2x,1x,1x,'~chi_10',1x,1x,'~chi_20'
499  & ,1x,1x,'~chi_30',1x,1x,'~chi_40',1x,1x,'~chi_1+',1x
500  & ,1x,'~chi_2+'/1x,'*',3x,1x,7(f8.2,1x))
501  6000 FORMAT(1x,'*',4x,4x,'h0',2x,1x,4x,'H0',2x,1x,4x,'A0',2x
502  & ,1x,4x,'H+'/1x,'*',3x,1x,5(f8.2,1x))
503  6050 FORMAT(1x,'*'/1x,'*',4x,4x,'h0',2x,1x,4x,'H0',2x,1x,4x,'A0',2x
504  & ,1x,4x,'H+'/1x,'*',3x,1x,5(f8.2,1x),3x,'(Before FeynHiggs)')
505  6100 FORMAT(1x,'*',11x,'|',3x,'~B',3x,'|',2x,'~W_3',2x,'|',2x
506  & ,'~H_1',2x,'|',2x,'~H_2',2x,'|'/1x,'*',3x,'~chi_10',1x,4('|'
507  & ,1x,f6.3,1x),'|'/1x,'*',3x,'~chi_20',1x,4('|'
508  & ,1x,f6.3,1x),'|'/1x,'*',3x,'~chi_30',1x,4('|'
509  & ,1x,f6.3,1x),'|'/1x,'*',3x,'~chi_40',1x,4('|'
510  & ,1x,f6.3,1x),'|')
511  6200 FORMAT(1x,'*'/1x,'*',6x,'L',4x,'|',3x,'~W',3x,'|',3x,'~H',3x,'|'
512  & ,12x,'R',4x,'|',3x,'~W',3x,'|',3x,'~H',3x,'|'/1x,'*',3x
513  & ,'~chi_1+',1x,2('|',1x,f6.3,1x),'|',9x,'~chi_1+',1x,2('|',1x
514  & ,f6.3,1x),'|'/1x,'*',3x,'~chi_2+',1x,2('|',1x,f6.3,1x),'|',9x
515  & ,'~chi_2+',1x,2('|',1x,f6.3,1x),'|')
516  6300 FORMAT(1x,'*'/1x,'*',8x,'|',2x,'~b_L',2x,'|',2x,'~b_R',2x,'|',8x
517  & ,'|',2x,'~t_L',2x,'|',2x,'~t_R',2x,'|',10x
518  & ,'|',1x,'~tau_L',1x,'|',1x,'~tau_R',1x,'|'/
519  & 1x,'*',3x,'~b_1',1x,2('|',1x,f6.3,1x),'|',3x,'~t_1',1x,2('|'
520  & ,1x,f6.3,1x),'|',3x,'~tau_1',1x,2('|',1x,f6.3,1x),'|'/
521  & 1x,'*',3x,'~b_2',1x,2('|',1x,f6.3,1x),'|',3x,'~t_2',1x,2('|'
522  & ,1x,f6.3,1x),'|',3x,'~tau_2',1x,2('|',1x,f6.3,1x),'|')
523  6400 FORMAT(1x,'*',3x,'A_b = ',f8.2,4x,'A_t = ',f8.2,4x,'A_tau = ',f8.2
524  & ,4x,'Alpha_GUT = ',f8.2)
525  6450 FORMAT(1x,'*',3x,'Alpha_Higgs = ',f8.4)
526  6500 FORMAT(1x,32('*'),1x,'END OF PYSUGI',1x,31('*'))
527 
528  9999 RETURN
529  END