English
资源说明：Spectrum and decays in the Minimal Composite Supersymmetric Standard Model
This program calculates the spectrum of the Minimal Composite Supersymmetric Standard Model.
It was hacked by Csaba Csaki and John Terning, based on the code NMSSMTools by
Ulrich Ellwanger, John F. Gunion, Cyril Hugonie, C.-C. Jean-Louis, Debottam Das, and Ana M. Teixeira
for more information on NMSSMTools see
http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html
For those familiar with NMSSMTools we have kept the same file names and structure.
HOW TO USE MCSSMTOOLS:
COMPILATION:
On Mac OS X you will need a modern fortran compiler, which can be downloaded
from http://hpc.sourceforge.net/ .
To compile, type first "make init", then "make". A first compilation
may take a while, since all subroutines of micromegas_2.2 are compiled.
The following 8 executable routines are created in the directory
"main": nmhdecay, nmhdecay_rand, nmhdecay_grid, nmspec, nmspec_rand,
nmspec_grid, nmgmsb and nmgmsb_rand.
If a subroutine in the directory "sources" was modified, one has to
type "make init" and "make" again. If a routine in the directory "main"
was modified, it suffices to type "make" again.
To delete all the already compiled codes type "make clean".
INPUT FILES:
Any name is allowed for the input file, provided it
contains the three letters "inp"; it can be of the general form
PREFIXinpSUFFIX where PREFIX and SUFFIX can contain dots etc..
The input file can be located in any directory specified by a PATH.
To run any input file PREFIXinpSUFFIX, type "run path/PREFIXinpSUFFIX",
or "./ run path/PREFIXinpSUFFIX" if the current directory is not in your $PATH
(path is optional; if absent, the input file has to be located in the
same directory as the script file "run".)
The output files are located in the directory specified by PATH.
They have the following format:
If one single point in the parameter space is evaluated:
PREFIXspectrSUFFIX, PREFIXdecaySUFFIX, PREFIXlhcsigSUFFIX and
PREFIXomegaSUFFIX (if the relic density is computed, see below)
If scans are performed:
PREFIXerrSUFFIX as well as PREFIXoutSUFFIX
However, the task to be performed by an input file must be specified in
the BLOCK MODSEL at the beginning (see the SLHA2 conventions in
B. Allanach et al., SUSY Les Houches Accord 2, arXiv:0801.0045
[hep-ph]).
The BLOCK MODSEL should contain the following four lines:
BLOCK MODSEL
9 I3 # Call micrOmegas default 0=no, 1=relic density only
13 I5 # 1: Sparticle decays via NMSDECAY
The meaning of the five integers I1, I2, I3, I4 and I5 is as follows:
I3=0: The dark matter relic density is not computed.
I3=1: The dark matter relic density is computed and checked via a call
of micromegas_2.2. This option is not possible for GMSB-like
boundary conditions. A first call of micromegas provokes the
compilation of additional subroutines, which may take a while.
In the case of a single point in parameter space (I2=0), the
relic density Omega*h^2 is given in the output files
PREFIXspectrSUFFIX as well as PREFIXomegaSUFFIX. The latter
contains in addition informations on the decomposition of the
LSP and the relevant annihilation/coannihilation processes.
The names of particles in the final states of the annihilation
and coannihilation processes are the same as in micrOMEGAS and
can be found in: G. Belanger, F. Boudjema, A. Pukhov and A. Semenov,
micrOMEGAs: A program for calculating the relic density
in the MSSM, Comput. Phys. Commun. 149 (2002) 103
[arXiv:hep-ph/0112278].
I3=2: Same as I3=1 + direct detection cross sections are computed.
In the case of a single point in parameter space (I2=0), the
BLOCK DIRECT DETECTION in PREFIXomegaSUFFIX contains:
csPsi = proton spin-independent cross section in [pb]
csNsi = neutron spin-independent cross section in [pb]
csPsd = proton spin-dependent cross section in [pb]
csNsd = neutron spin-dependent cross section in [pb]
I3=3: Same as I3=1 + the thermally averaged LSP annihilation cross section
as well as the resulting photon spectrum are computed. In the case of
a single point in parameter space (I2=0), these are written in the
BLOCK INDIRECT DETECTION of PREFIXomegaSUFFIX:
sigmaV = LSP annihilation cross section,
dN/dx = photon spectrum from LSP annihilation. N is the nb of photons
and x = log(E/M) where E is the photon energy and M the LSP mass.
I3=4: Same as I3=2+3.
Precision of the CP-even/odd/charged Higgs masses:
I4=0: 1-loop: complete contributions ~ top/bottom Yukawas
contributions ~ g1, g2, lambda and kappa to LLA
for the SM-like CP-even Higgs only
2-loop: top/bottom Yukawa contributions to LLA
I4=1: as in G. Degrassi, P. Slavich, Nucl.Phys.B825:119-150,2010,
arXiv:0907.4682 (with special thanks to P. Slavich);
corrections to the charged Higgs mass from K.H.Phan and P. Slavich:
1-loop: complete contributions ~ top/bottom Yukawas
complete contributions ~g1, g2, lambda and kappa
(except for pole masses)
2-loop: complete contributions ~ top/bottom Yukawas
I4=2: 1-loop: complete contributions ~ top/bottom Yukawas
complete contributions ~g1, g2, lambda and kappa
including pole masses (slow!)
2-loop: complete contributions ~ top/bottom Yukawas
Sparticle total widths and branching ratios:
I5=0: Not computed
I5=1: NMSDECAY is called, which computes sparticle 2-body and 3-body
branching ratios as in
SDECAY: A Fortran code for the decays of the supersymmetric
particles in the MSSM
by M. Muhlleitner (Karlsruhe, Inst. Technol.),
A. Djouadi (Orsay, LPT & CERN, Theory Division),
Y. Mambrini (Orsay, LPT),
Comput.Phys.Commun.168:46-70 (2005), hep-ph/0311167.
SDECAY should be cited whenever NMSDECAY is used.
In NMSDECAY.f in the directory sources, the flags
"flagmulti" (3-body decays)
"flagqcd" (QCD corrections to 2-body decays)
"flagloop" (loop decays)
can be switched off; otherwise a call of NMSDECAY takes about 2-3 seconds
per point in parameter space.
In MCSSMTools flagqcd is switched off.
In the versions nmhdecay.f and nmspec.f, the sparticle widths and BR's are
appended to the output file PREFIXdecaySUFFIX in SLHA2 format. If scans are
performed, the user can use the arguments of the COMMON statements in the
subroutines OUTPUT in order to define the content of the output file.
************************************************
Sample input file:
# Input file for MCSSMTools
# Based on SUSY LES HOUCHES ACCORD II
BLOCK MODSEL
9 0 # Call micrOmegas (default 0=no, 1=relic density only)
13 1 # 1: Sparticle decays via NMSDECAY
BLOCK SMINPUTS
1 127.92D0 # ALPHA_EM^-1(MZ)
2 1.16639D-5 # GF
3 .1172D0 # ALPHA_S(MZ)
4 91.187D0 # MZ
5 4.214D0 # MB(MB) (running mass)
6 171.4D0 # MTOP (pole mass)
7 1.777D0 # MTAU
BLOCK MINPAR
0 600.D0 # MSUSY (If =/= SQRT(2*MQ1+MU1+MD1)/2)
3 1.2D0 # TANB
BLOCK EXTPAR
1 800.D0 # M1 (If =/= M2/2)
2 1000.D0 # M2
3 20000.D0 # M3 (If =/= 3*M2)
12 10.D0 # AD3
13 10000.D0 # AE3
16 0.D0 # AE2 = AE1 (If =/= AE3)
33 10000.D0 # ML3
32 10000.D0 # ML2 = ML1 (If =/= ML3)
36 10000.D0 # ME3
35 10000.D0 # ME2 = ME1 (If =/= ME3)
43 300.D0 # MQ3
42 10000.D0 # MQ2 = MQ1 (If =/= MQ3)
46 300.D0 # MU3
45 20000.D0 # MU2 = MU1 (If =/= MU3)
49 10000.D0 # MD3
48 20000.D0 # MD2 = MD1 (If =/= MD3)
62 1.0D-9 # KAPPA
64 0.D0 # AKAPPA
70 1013.54D0 # A
200 -8.0D7 # Tad
201 60.D0 # f
202 400.D0 # mSing
203 300.D0 # majS
**************************************************
Content of the array PAR(I) (couplings and soft parameters at the SUSY scale):
PAR(1) = lambda, a.k.a. y, dynamical yukawa coupling y=SQRT(2) Mtop/(v sin(beta))
PAR(2) = kappa
PAR(3) = tan(beta)
PAR(4) = mu (effective mu term = lambda*s)
PAR(5) = Alambda = A/lambda
PAR(6) = Akappa
PAR(7) = mQ3**2
PAR(8) = mU3**2
PAR(9) = mD3**2
PAR(10) = mL3**2
PAR(11) = mE3**2
PAR(12) = AU3 = Alambda = A/lambda
PAR(13) = AD3
PAR(14) = AE3
PAR(15) = mQ2**2
PAR(16) = mU2**2
PAR(17) = mD2**2
PAR(18) = mL2**2
PAR(19) = mE2**2
PAR(20) = M1
PAR(21) = M2
PAR(22) = M3
PAR(23) = MA (diagonal doublet CP-odd mass matrix element)
PAR(24) = MP (diagonal singlet CP-odd mass matrix element)
PAR(25) = AE2
* Additional input parameters
* Tad = linear soft breaking term
* f = linear superpotential term for the singlet S
* mSing2 = soft breaking singlet mass
* majS = Majorana singlet mass
*
Content of the array PROB(I) (phenomenological and theoretical constraints):
PROB(I) = 0, I = 1..45: OK
PROB(1) =/= 0 chargino too light
PROB(2) =/= 0 excluded by Z -> neutralinos
PROB(3) =/= 0 charged Higgs too light
PROB(4) =/= 0 excluded by ee -> hZ
PROB(5) =/= 0 excluded by ee -> hZ, h -> bb
PROB(6) =/= 0 excluded by ee -> hZ, h -> tautau
PROB(7) =/= 0 excluded by ee -> hZ, h -> invisible
PROB(8) =/= 0 excluded by ee -> hZ, h -> 2jets
PROB(9) =/= 0 excluded by ee -> hZ, h -> 2photons
PROB(10) =/= 0 excluded by ee -> hZ, h -> AA -> 4bs
PROB(11) =/= 0 excluded by ee -> hZ, h -> AA -> 4taus
PROB(12) =/= 0 excluded by ee -> hZ, h -> AA -> 2bs 2taus
PROB(13) =/= 0 excluded by Z -> hA (Z width)
PROB(14) =/= 0 excluded by ee -> hA -> 4bs
PROB(15) =/= 0 excluded by ee -> hA -> 4taus
PROB(16) =/= 0 excluded by ee -> hA -> 2bs 2taus
PROB(17) =/= 0 excluded by ee -> hA -> AAA -> 6bs
PROB(18) =/= 0 excluded by ee -> hA -> AAA -> 6taus
PROB(19) =/= 0 excluded by ee -> Zh -> ZAA -> Z + light pairs
PROB(20) =/= 0 excluded by stop -> b l sneutrino
PROB(21) =/= 0 excluded by stop -> neutralino c
PROB(22) =/= 0 excluded by sbottom -> neutralino b
PROB(23) =/= 0 squark/gluino too light
PROB(24) =/= 0 selectron/smuon too light
PROB(25) =/= 0 stau too light
PROB(26) =/= 0 lightest neutralino is not LSP
PROB(27) =/= 0 Landau Pole in l, k, ht, hb below MGUT
PROB(28) =/= 0 unphysical global minimum
PROB(29) =/= 0 Higgs soft masses >> Msusy
PROB(30) =/= 0 excluded by WMAP (checked only if OMGFLAG=1)
PROB(31) =/= 0 eff. Higgs self-couplings in Micromegas > 1
PROB(32) =/= 0 b->s gamma more than 2 sigma away
PROB(33) =/= 0 Delta M_s more than 2 sigma away
PROB(34) =/= 0 Delta M_d more than 2 sigma away
PROB(35) =/= 0 B_s->mu+mu- more than 2 sigma away
PROB(36) =/= 0 B+-> tau+nu_tau more than 2 sigma away
PROB(37) =/= 0 (g-2)_muon more than 2 sigma away
PROB(38) =/= 0 excluded by Upsilon(1S) -> A gamma
PROB(39) =/= 0 excluded by eta_b(1S) mass difference
PROB(40) =/= 0 BR(B-->X_s mu+ mu-) more than 2 sigma away
PROB(41) =/= 0 excluded by ee -> hZ, h -> AA -> 4taus (new ALEPH analysis)
PROB(42) =/= 0 excluded by top -> b H+, H+ -> c s (CDF, D0)
PROB(43) =/= 0 excluded by top -> b H+, H+ -> tau nu_tau (D0)
PROB(44) =/= 0 excluded by top -> b H+, H+ -> W+ A1, A1 -> 2taus (CDF)
PROB(45) =/= 0 excluded by LHC: A/H -> 2taus
Output parameters:
The decay and spectrum output files can be visualized using the
spectrum program, available at http://bit.ly/mcspect.
SMASS(1-3): CP-even masses (ordered)
SCOMP(1-3,1-3): Mixing angles: if HB(I) are the bare states,
HB(I) = Re(H1), Re(H2), Re(S), and HM(I) are the mass eigenstates,
the convention is HB(I) = SUM_(J=1,3) SCOMP(J,I)*HM(J)
which is equivalent to HM(I) = SUM_(J=1,3) SCOMP(I,J)*HB(J)
PMASS(1-2): CP-odd masses (ordered)
PCOMP(1-2,1-2): Mixing angles: if AB(I) are the bare states,
AB(I) = Im(H1), Im(H2), Im(S), and AM(I) are the mass eigenstates,
the convention is
AM(I) = PCOMP(I,1)*(COSBETA*AB(1)+SINBETA*AB(2))
+ PCOMP(I,2)*AB(3)
CMASS: Charged Higgs mass
CU,CD,CV,CJ,CG(i) Reduced couplings of h1,h2,h3 (i=1,2,3) or
a1,a2 (i=4,5) to up type fermions, down type
fermions, gauge bosons, gluons and photons
Note: CV(4)=CV(5)=0
WIDTH(i) Total decay width of h1,h2,h3,a1,a2 (i=1..5)
with the following branching ratios:
BRJJ(i) h1,h2,h3,a1,a2 -> gluon gluon
BRMM(i) " -> mu mu
BRLL(i) " -> tau tau
BRSS(i) " -> ss
BRCC(i) " -> cc
BRBB(i) " -> bb
BRTT(i) " -> tt
BRWW(i) " -> WW (BRWW(4)=BRWW(5)=0)
BRZZ(i) " -> ZZ (BRZZ(4)=BRZZ(5)=0)
BRGG(i) " -> gamma gamma
BRZG(i) " -> Z gamma
BRHIGGS(i) (i=1..5) -> other Higgses, including:
BRHAA(i,j) hi -> a1a1, a1a2, a2a2 (i=1..3, j=1..3)
BRHCHC(i) hi -> h+h- (i=1..3)
BRHAZ(i,j) hi -> Zaj (i=1..3)
BRHCW(i) h1,h2,h3 -> h+W- (i=1..3), a1,a2 -> h+W- (i=4,5)
BRHHH(i) h2 -> h1h1, h3-> h1h1, h1h2, h2h2 (i=1..4)
BRAHA(i) a2 -> a1hi (i=1..3)
BRAHZ(i,j) ai -> Zhj (i=1,2, j=1..3)
BRSUSY(i) (i=1..5) -> susy particles, including:
BRNEU(i,j,k) -> neutralinos j,k (i=1..5, j,k=1..5)
BRCHA(i,j) -> charginos 11, 12, 22 (i=1..5, j=1..3)
BRHSQ(i,j) hi -> uLuL, uRuR, dLdL, dRdR, t1t1, t2t2,
t1t2, b1b1, b2b2, b1b2 (i=1..3, j=1..10)
BRASQ(i,j) ai -> t1t2, b1b2 (i=1,2, j=1,2)
BRHSL(i,j) hi -> lLlL, lRlR, nLnL, l1l1, l2l2, l1l2,
ntnt (i=1..3, j=1..7)
BRASL(i) ai -> l1l2 (i=1,2)
HCWIDTH Total decay width of the charged Higgs
with the following branching ratios:
HCBRM h+ -> mu nu_mu
HCBRL " -> tau nu_tau
HCBRSU " -> s u
HCBRBU " -> b u
HCBRSC " -> s c
HCBRBC " -> b c
HCBRBT " -> b t
HCBRWHT " -> neutral Higgs W+, including:
HCBRWH(i) " -> H1W+, H2W+, h3W+, a1W+, a2W+ (i=1..5)
HCBRSUSY " -> susy particles,including
HCBRNC(i,j)" -> neutralino i chargino j (i=1..5, j=1,2)
HCBRSQ(i) " -> uLdL, t1b1, t1b2, t2b1, t2b2 (i=1..5)
HCBRSL(i) " -> lLnL, t1nt, t2nt (i=1..3)
MNEU(i) Mass of neutralino chi_i (i=1,5, ordered in mass)
NEU(i,j) chi_i components of bino, wino, higgsino u&d, singlino
(i,j=1..5)
MCHA(i) Chargino masses
U(i,j),V(i,j) Chargino mixing matrices
Significances for Higgs detection at the LHC:
At low luminosity (30 fb^-1): in LOWSIG(X,Y), where
X=1: h1
X=2: h2
X=3: h3
X=4: a1
X=5: a2
Y=1: channel bbh/a -> bbtautau
Y=2: channel gg -> h/a -> gamma gamma
Y=3: channel gg -> h -> ZZ -> 4 leptons
Y=4: channel gg -> h -> WW -> 2 leptons 2 neutrinos
Y=5: channel WW -> h -> tautau
Y=6: channel WW -> h -> WW
Y=7: channel WW -> h -> gamma gamma
At high luminosity (300 fb^-1): in HIGSIG(X,Y), where X as above,
Y=1: channel h/a -> gamma gamma
Y=2: channel h/a -> gamma gamma lepton
Y=3: channel tth/a -> bb + X
Y=4: channel bbh/a -> bbtautau
Y=5: channel gg -> h -> ZZ -> 4 leptons
Y=6: channel gg -> h -> WW -> 2 leptons 2 neutrinos
Y=7: channel WW -> h -> tautau
Y=8: channel WW -> h -> WW
Y=9: channel WW -> h -> invisible
本源码包内暂不包含可直接显示的源代码文件，请下载源码包。
