12 IMPLICIT DOUBLE PRECISION(a-
h, o-
z)
15 parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
16 &kexcit=4000000,kdimen=5000000)
20 parameter(doc=
'01 May 2006')
23 REAL mzero,mhlf,azero,tanb,sgnmu,mtop
25 COMMON /sugxin/ xisain(24),xsugin(7),xgmin(14),xnrin(4),
27 REAL xisain,xsugin,xgmin,xnrin,xamin
32 COMMON /sspar/ super(72)
33 COMMON /sugmg/ mss(32),gss(31),mgutss,ggutss,agutss,ftgut,
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
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)
100 common/pylh3p/modsel(200),parmin(100),parext(200),rmsoft(0:100),
101 & au(3,3),ad(3,3),ae(3,3)
103 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
105 common/
pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
111 CHARACTER*20 chmod(5)
114 COMMON /sugnu/ xnusug(18)
118 DATA chmod/
'mSUGRA',
'mGMSB',
'non-universal SUGRA',
119 &
'truly unified SUGRA',
'non-minimal GMSB'/
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'
136 IF (imss(1).EQ.12)
THEN
139 ELSEIF(imss(1).EQ.13)
THEN
144 WRITE(mstu(11),*)
'(PYSUGI:) No valid unit given in IMSS(20)'
147 WRITE(mstu(11),*)
'READING SUSY MODEL FROM FILE...'
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:'
156 IF (imodel.EQ.4)
THEN
160 IF (imodel.EQ.1.OR.imodel.EQ.3)
THEN
161 WRITE(mstu(11),*)
'ENTER M_0, M_(1/2), A_0, tan(beta),'
163 READ(lfn,*) m0,mhf,a0,tanb,sgnmu,
mt
164 IF (imodel.EQ.3)
THEN
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'
174 &
' NUSUG5 = GUT scale 3rd generation masses'
176 IF (inusug.EQ.0)
THEN
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'
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),'
199 READ(lfn,*) xnusug(18),xnusug(16),xnusug(17),
200 & xnusug(15),xnusug(14)
204 ELSEIF (imodel.EQ.2.OR.imodel.EQ.5)
THEN
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
213 amgvss=m0*mhf*xcmgv/sqrt(3d0)/ampl
214 IF (imodel.EQ.5)
THEN
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',
222 WRITE(mstu(11),*)
'n5_1,n5_2,n5_3 = n5 values for U(1),'
224 WRITE(mstu(11),*)
'ENTER Rsl, dmH_d2, dmH_u2, d_Y, n5_1,'
226 READ(lfn,*) xgmin(8),xgmin(9),xgmin(10),xgmin(11),xgmin(12),
227 $ xgmin(13),xgmin(14)
230 WRITE(mstu(11),*)
'Invalid model choice.'
247 CALL
sugra(mzero,mhlf,azero,tanb,sgnmu,mtop,imodel)
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.')
263 IF (nogood.EQ.5) CALL
pyerrm(16
264 & ,
'(PYSUGI:) SUSY parameters give ~chi_10 not LSP.')
267 IF (itachy.NE.0) CALL
pyerrm(16,
268 &
'(PYSUGI:) Tachyonic sleptons at GUT scale.')
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),
283 rmsoft(1)=dble(gss(7))
284 rmsoft(2)=dble(gss(8))
285 rmsoft(3)=dble(gss(9))
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))
295 rmss(10)=0.5*(super(14)+super(10))
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))
342 rmss(21)=
max(rmss(21),dble(super(66)))
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))
351 pmas(
pycomp(ksusy1+21),1)=abs(super(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))
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))
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))
379 pmas(
pycomp(ksusy1+24),1)=abs(super(51))
380 pmas(
pycomp(ksusy1+37),1)=abs(super(52))
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))
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))
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))
419 IF (
mstp(122).NE.0)
THEN
422 WRITE(mstu(11),5000) doc, isaver
424 IF (imodel.EQ.1)
THEN
425 WRITE(mstu(11),5200) mzero, mhlf, azero, tanb, nint(sgnmu),
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)
436 WRITE(mstu(11),6000) (super(ip),ip=55,58)
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)
446 WRITE(mstu(11),6450) rmss(18)
448 WRITE(mstu(11),5500)
'Couplings'
449 WRITE(mstu(11),6400) rmss(15),rmss(16),rmss(17),rmss(20)
454 IF (imss(4).EQ.3)
THEN
455 IF (
mstp(122).NE.0)
WRITE(mstu(11),
'(1x,"*"/1x,"*",A)')
456 &
' (PYSUGI:) Now calling FeynHiggs.'
460 IF (
mstp(122).NE.0)
THEN
463 &
'Corrected Higgs masses and mixing'
464 WRITE(mstu(11),6000) pmas(25,1),pmas(35,1),pmas(36,1),
466 WRITE(mstu(11),6450) rmss(18)
472 IF (
mstp(122).NE.0)
WRITE(mstu(11),6500)
476 imss(4)=
max(2,imss(4))
478 5000
FORMAT(1
x,19(
'*'),1
x,
'PYSUGI v1.52: PYTHIA/ISASUSY '
479 & ,
'INTERFACE',1
x,19(
'*')/1
x,
'*',3
x,
'PYSUGI: Last Change',1
x,a
480 & ,1
x,
'-',1
x,
'P. Skands / S. Mrenna'/1
x,
'*',2
x,a/1
x,
'*')
481 5100
FORMAT(1
x,
'*',1
x,
'ISASUSY Input:'/1
x,
'*',1
x,
'----------------')
482 5200
FORMAT(1
x,
'*',1
x,3
x,
'M_0',6
x,
'M_1/2',5
x,
'A_0',3
x,
'Tan(beta)',
483 & 3
x,
'Sgn(mu)',3
x,
'M_t'/1
x,
'*',1
x,4(f8.2,1
x),i8,2
x,f8.2)
484 5300
FORMAT(1
x,
'*'/1
x,
'*',1
x,
'ISASUSY Output:'/1
x,
'*',1
x
485 & ,
'----------------')
486 5400
FORMAT(1
x,
'*',1
x,a)
487 5500
FORMAT(1
x,
'*',1
x,a,
':')
488 5600
FORMAT(1
x,
'*',2
x,2
x,
'M_GUT',2
x,2
x,
'g_GUT',2
x,1
x,
'alpha_GUT'/
489 & 1
x,
'*',2
x,1
p,2(1
x,e8.2),2
x,e8.2)
490 5700
FORMAT(1
x,
'*',4
x,4
x,
'~u',2
x,1
x,4
x,
'~d',2
x,1
x,4
x,
'~s',2
x,1
x,
491 & 4
x,
'~c',2
x,1
x,4
x,
'~b',2
x,1
x,2
x,
'~b(12)',1
x,4
x,
'~t',2
x,1
x, 2
x,
492 &
'~t(12)'/1
x,
'*',2
x,
'L',1
x,8(f8.2,1
x)/1
x,
'*',2
x,
'R',1
x,8(f8.2
494 5800
FORMAT(1
x,
'*'/1
x,
'*',4
x,4
x,
'~e',2
x,1
x,3
x,
'~mu',2
x,1
x,3
x,
'~tau',1
x
495 & ,1
x,
'~tau(12)',1
x,2
x,
'~nu_e',1
x,1
x,1
x,
'~nu_mu',1
x,1
x,1
x
496 & ,
'~nu_tau'/1
x,
'*',2
x,
'L',1
x,7(f8.2,1
x)/1
x,
'*',2
x,
'R',1
x,4(f8
498 5900
FORMAT(1
x,
'*'/1
x,
'*',4
x,4
x,
'~g',2
x,1
x,1
x,
'~chi_10',1
x,1
x,
'~chi_20'
499 & ,1
x,1
x,
'~chi_30',1
x,1
x,
'~chi_40',1
x,1
x,
'~chi_1+',1
x
500 & ,1
x,
'~chi_2+'/1
x,
'*',3
x,1
x,7(f8.2,1
x))
501 6000
FORMAT(1
x,
'*',4
x,4
x,
'h0',2
x,1
x,4
x,
'H0',2
x,1
x,4
x,
'A0',2
x
502 & ,1
x,4
x,
'H+'/1
x,
'*',3
x,1
x,5(f8.2,1
x))
503 6050
FORMAT(1
x,
'*'/1
x,
'*',4
x,4
x,
'h0',2
x,1
x,4
x,
'H0',2
x,1
x,4
x,
'A0',2
x
504 & ,1
x,4
x,
'H+'/1
x,
'*',3
x,1
x,5(f8.2,1
x),3
x,
'(Before FeynHiggs)')
505 6100
FORMAT(1
x,
'*',11
x,
'|',3
x,
'~B',3
x,
'|',2
x,
'~W_3',2
x,
'|',2
x
506 & ,
'~H_1',2
x,
'|',2
x,
'~H_2',2
x,
'|'/1
x,
'*',3
x,
'~chi_10',1
x,4(
'|'
507 & ,1
x,f6.3,1
x),
'|'/1
x,
'*',3
x,
'~chi_20',1
x,4(
'|'
508 & ,1
x,f6.3,1
x),
'|'/1
x,
'*',3
x,
'~chi_30',1
x,4(
'|'
509 & ,1
x,f6.3,1
x),
'|'/1
x,
'*',3
x,
'~chi_40',1
x,4(
'|'
511 6200
FORMAT(1
x,
'*'/1
x,
'*',6
x,
'L',4
x,
'|',3
x,
'~W',3
x,
'|',3
x,
'~H',3
x,
'|'
512 & ,12
x,
'R',4
x,
'|',3
x,
'~W',3
x,
'|',3
x,
'~H',3
x,
'|'/1
x,
'*',3
x
513 & ,
'~chi_1+',1
x,2(
'|',1
x,f6.3,1
x),
'|',9
x,
'~chi_1+',1
x,2(
'|',1
x
514 & ,f6.3,1
x),
'|'/1
x,
'*',3
x,
'~chi_2+',1
x,2(
'|',1
x,f6.3,1
x),
'|',9
x
515 & ,
'~chi_2+',1
x,2(
'|',1
x,f6.3,1
x),
'|')
516 6300
FORMAT(1
x,
'*'/1
x,
'*',8
x,
'|',2
x,
'~b_L',2
x,
'|',2
x,
'~b_R',2
x,
'|',8
x
517 & ,
'|',2
x,
'~t_L',2
x,
'|',2
x,
'~t_R',2
x,
'|',10
x
518 & ,
'|',1
x,
'~tau_L',1
x,
'|',1
x,
'~tau_R',1
x,
'|'/
519 & 1
x,
'*',3
x,
'~b_1',1
x,2(
'|',1
x,f6.3,1
x),
'|',3
x,
'~t_1',1
x,2(
'|'
520 & ,1
x,f6.3,1
x),
'|',3
x,
'~tau_1',1
x,2(
'|',1
x,f6.3,1
x),
'|'/
521 & 1
x,
'*',3
x,
'~b_2',1
x,2(
'|',1
x,f6.3,1
x),
'|',3
x,
'~t_2',1
x,2(
'|'
522 & ,1
x,f6.3,1
x),
'|',3
x,
'~tau_2',1
x,2(
'|',1
x,f6.3,1
x),
'|')
523 6400
FORMAT(1
x,
'*',3
x,
'A_b = ',f8.2,4
x,
'A_t = ',f8.2,4
x,
'A_tau = ',f8.2
524 & ,4
x,
'Alpha_GUT = ',f8.2)
525 6450
FORMAT(1
x,
'*',3
x,
'Alpha_Higgs = ',f8.4)
526 6500
FORMAT(1
x,32(
'*'),1
x,
'END OF PYSUGI',1
x,31(
'*'))