HOW TO USE NMSSMTOOLS: After the download of NMSSMTools_x.tgz (where x is the version number), type "tar zxvf NMSSMTools_x.tgz". This will create a directory named NMSSMTools_x. Go into this directory. There you find the directories -- "sources", "main" and "microRun", which contain source files; -- "EXPCON", which contains data files corresponding to experimental constraints; -- "SAMPLES", which contains sample input and output files; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% COMPILATION: To compile, type first "make init", then "make". A first compilation may take a while, since many subroutines of micromegas are compiled. The following executable routines are created in the directory "main": nmhdecay, nmhdecay_rand, nmhdecay_grid, nmhdecay_mcmc, nmspec, nmspec_rand, nmspec_grid, nmspec_mcmc, nmgmsb, nmgmsb_rand, nmspec_grid, nmspec_mcmc, nmhdecayCPV, nmhdecayCPV_rand, nmhdecayCPV_grid, nmhdecayCPV_mcmc. 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: From the version 2.0.0 onwards, the task to be performed by an input file is independent from its name. Any name is allowed, 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". (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 (includes the previous PREFIXdecaySUFFIX) 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 at least the first four among the following lines: BLOCK MODSEL 3 1 # NMSSM PARTICLE CONTENT 1 I0 # IMOD (0=general NMSSM, 1=SUGRA, 2=GMSB # 3=Z3 inv NMSSM, 4=Z3 inv SUGRA, 5=Z3 inv GMSB) 5 I1 # CPV (0=no, 2=yes) 10 I2 # ISCAN (0=no scan, 1=grid scan, 2=random scan, 3=MCMC) 9 I3 # |OMGFLAG|=0: no (default), =1: relic density only, # =2: dir. det. rate, =3: indir. det. rate, =4: both, # OMGFLAG>0: 0.107 VV,VV* (default); 1: H->VV,VV*,V*V* 15 I9 # Precision for micromegas (defalt=0): # +0/1: fast computation on/off # +0/2: Beps=1d-3, 1d-6 # +0/4: virtual Ws off/on 16 I10 # 0: Output contains privately defined BLOCK's (default) # 1: Restricts the output to BLOCK's defined by SLHA1/2 # (useful as param_card.dat for MadGraph) 17 I11 # Theoretical constraints (1=yes, 0=no, default=1) 18 I12 # LEP/Tevatron constraints(1=yes, 0=no, default=1) 19 I13 # LHC Higgs constraints (1=yes, 0=no, default=1) 20 I14 # Flavor constraints (1=yes, 0=no, default=1) 21 I15 # EDM constraints (1=yes, 0=no, default=1) 21 I16 # CMS charg(neutral)ino constraints (1=yes, 0=no, default=1) The meaning of the integers I1..I10 is as follows: I0=0: general NMSSM with parameters specified at the SUSY breaking scale (an average of the squark masses, unless specified by the user). As in the previous versions of NMSSMTools, the input parameters have to be given in the BLOCK EXTPAR following the SLHA2 conventions. Z_3 violating terms are allowed. See the example in SAMPLES/inp.dat. I0=1: mSUGRA-like boundary conditions at the GUT scale with universal scalar masses m0, gaugino masses M12 and trilinear couplings A0 (to be specified in the BLOCK MINPAR). Z_3 violating terms are allowed. See the example in SAMPLES/inpsp.dat. I0=2: GMSB-like boundary conditions as in the paper U. Ellwanger et al., arXiv:0803.2962 [hep-ph]. Z_3 violating terms are allowed See the example in SAMPLES/inpgm.dat. I0=3: Z_3 invariant NMSSM with parameters specified at the SUSY breaking scale (an average of the squark masses, unless specified by the user). As in the previous versions of NMSSMTools, the input parameters have to be given in the BLOCK EXTPAR following the SLHA2 conventions. Z_3 violating terms as input generate an error message. See the example in SAMPLES/inpZ3.dat. I0=4: Z_3 invariant mSUGRA-like boundary conditions at the GUT scale with universalscalar masses m0, gaugino masses M12 and trilinear couplings A0 (to be specified in the BLOCK MINPAR). Z_3 violating terms as input generate an error message. See the example in SAMPLES/inpsp.dat. I0=5: Z_3 invariant GMSB-like boundary conditions. Z_3 violating terms as input generate an error message. I1=0: CP conserving parameters I1=2: CP violating parameters I2=1: A scan over a grid in parameter space is performed. The boundaries in parameter space as well as the corresponding numbers of steps have to be specified. For possible scans in the general NMSSM see the example SAMPLES/gridinp.dat. For possible scans with mSUGRA-like boundary conditions see the example SAMPLES/gridinpsp.dat. For possible scans with GMSB-like boundary conditions see the example SAMPLES/gridinpgm.dat The output file PREFIXerrSUFFIX gives the number of points which have passed all constraints, and the corresponding range of input parameters. The output file PREFIXoutSUFFIX contains details of points which have passed constraints; the latter can be modified by the user by editing the corresponding routines in the directory "main": The output file PREFIXoutSUFFIX is created in the subroutine OUTPUT near the end of the file nmhdecay_grid.f (general NMSSM), or nmspec_grid.f (mSUGRA) or nmgmsb_grid.f (GMSB). If the first line in the subroutine OUTPUT reads "IF(IFAIL.EQ.0)THEN", properties of points with phenomenological problems are not written into PREFIXoutSUFFIX. If this line reads "IF(IFAIL.EQ.0 .OR. IFAIL.EQ.10)THEN", points with phenomenological problems are written as well. The properties which are listed in PREFIXoutSUFFIX depend on the content of the array RES and should be specified by the user. The meaning of the various arrays containing Higgs and sparticle masses and mixing angles, Higgs branching ratios is given at the beginning of the file nmhdecay_grid.f (general NMSSM), nmspec_grid.f (mSUGRA) or nmspec_grid.f (GMSB). For convenience, we list the content of the array PAR(I) (the couplings and soft terms at the SUSY scale) as well as the content of the array PROB(I) (phenomenological and some theoretical constraints) below. I2=2: A random scan in parameter space is performed. The boundaries in parameter space as well as the total number of steps have to be specified. For possible scans in the general NMSSM see the example SAMPLES/randinp.dat. For possible scans with mSUGRA-like boundary conditions see the example SAMPLES/randinpsp.dat. For possible scans with GMSB-like boundary conditions see the example SAMPLES/randinpgm.dat For the content of the output files PREFIXerrSUFFIX and PREFIXoutSUFFIX (and the treatment of the latter) see the description of the case I2=1 above; the corresponding MAIN routines for random scans in the directory "main" are denoted by nmhdecay_rand.f, nmspec_rand.f and nmgmsb_rand.f. I2=3: A Markov Chain Monte Carlo scan in parameter space is performed. The relative and the minimum size of each step for all parameters and the total number of steps have to be specified. For possible scans in the general NMSSM see the example SAMPLES/mcmcinp.dat. For possible scans with mSUGRA-like boundary conditions see the example SAMPLES/mcmcinpsp.dat. For possible scans with GMSB-like boundary conditions see the examples SAMPLES/mcmcinpgm.dat. For the content of the output files PREFIXerrSUFFIX and PREFIXoutSUFFIX (and the treatment of the latter) see the description of the case I2=1 above; the corresponding MAIN routine for a MCMC scans in the directory "main" are denoted by nmhdecay_mcmc.f, nmspec_mcmc.f and nmgmsb_mcmc.f. 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. A first call of micromegas provokes the compilation of additional subroutines, which may take a while. For GMSB models, only the non-thermal relic density of the gravitino (coming from NLSP to gravitino decay) is computed, i.e. micromegas computes the NLSP relic density and multiplies the result by the mass ratio of gravitino over NLSP. 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]. If I3>0 constraints on relic density from Planck +/- 10% are checked: 0.107 < Omega h^2 < 0.131 If I3<0 only the upper bound is imposed: Omega h^2 < 0.131 |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 NDMCROSSSECT 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] Constraints from LUX (arXiv:1310.8214) are checked |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 ANNIHILATION of PREFIXomegaSUFFIX: sigmaV = LSP annihilation cross section. |I3|=4: Same as |I3|=2+3. I4=0: Precision of the CP-even/odd/charged Higgs masses: 1-loop: complete contributions ~ top/bottom Yukawas contributions ~ g1, g2, lambda and kappa to LLA 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*alpha_S 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*alpha_S I4=1,2: for Z_3 invariant NMSSM only (I0>2). I4=3: same as I4=0 + computation of Higgs mass uncertainty by varying the RG scale by *2 and /2. I4=4: same as I4=1 + computation of Higgs mass uncertainty by varying the RG scale by *2 and /2. I4=5: same as I4=2 + computation of Higgs mass uncertainty by varying the RG scale by *2 and /2. I4=6: same as I4=3 + check that uncertainty is < 3 GeV. I4=7: same as I4=4 + check that uncertainty is < 3 GeV. I4=8: same as I4=5 + check that uncertainty is < 3 GeV. I4>2: for general NMSSM at the SUSY scale only (I0=1,3) I5=1: Constraints on (g-2)_muon are imposed (default value) I5=0: They are not imposed (for scanning versions only, I2=/=0) I6=0: no particular GMSB model assumed above the messenger scale MS and Alambda can be given as free input parameters at MMESS I6=1,3: Delgado/Giudice/Slavich model (from arXiv:0706.3873[hep-ph]) with unified singlet-messenger coupling at the GUT scale xi_U The value of xi_U is given in BLOCK EXTPAR, switch 73 I6=2,4: Delgado/Giudice/Slavich model without universality with a triplet coupling LTT and a doublet coupling LPP given at the messenger scale in BLOCK EXTPAR, switches 74, 75 I6=1,2: Evans/Ibe/Yanagida model (from arXiv:1107.3006[hep-ph]) without messenger-Higgs mixing. The parameters LU, LD, LT LB, LL are given at MMESS in BLOCK EXTPAR, switches 76 to 80. I6=1,2: Evans/Ibe/Yanagida model with messenger-Higgs mixing. LT, LB are given at MMESS in BLOCK EXTPAR, switches 78, 79. The others are given by LD=L*LT/HT, LU=L*LB/HB, LL=HL*LB/HB If I6=/=0: ALAMBDA cannot be given at MMESS as it is related to the extra messenger couplings above MMESS. XIF=XIS=MU'=MS'^2=0. KAPPA and MUEFF and MS are computed from the minimisation eqs. On the other hand MS is also related to the extra couplings. DMIN, the maximal relative deviation allowed between the computed and expected value for MS, has to be given in BLOCK EXTPAR (switch 0). If the deviation is larger than DMIN then: IFAIL=21 (if previously = 0) IFAIL=22 (if previously = 10). I7=0: Sparticle total widths and branching ratios not computed I7=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 the versions nmhdecay.f and nmspec.f, the sparticle widths and BR's are appended to the output file PREFIXspectrSUFFIX 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. I8=0: No double-offshell H->V*V* decays (faster) I8=1: Double-offshell H->V*V* decays (slower) I9=0+0/1: micromegas fast computation of resonances on/off +0/2: minimal Boltzman suppression for coannihilation = 1d-3/1d-6 +0/4: virtual Ws off/on I10=1: The name of the output file PREFIXspectrSUFFIX remains unchanged, but it can be renamed as param_card.dat, useful for MadGraph. At least for the version MG5_aMC_v2_2_3 the following procedure seems to work: 1) Place a proc_card.dat in the directory MG5_aMC_v2_2_3. The proc_card.dat should contain a line "import model nmssm", possibly adding "-modelname" (which changes the syntax for Higgs bosons into h01, h02, h03, a01, a02). Generate the process typing "./bin/mg5aMC proc_card.dat". This generates a directory "PROC_nmssm_x" where "x" is a number. 2) Go into this directory. Replace the (default) param_card.dat in the subdirectory "Cards" by your desired param_card.dat. Edit the run_card.dat if desired. 3) Generate events typing "./bin/generate_events" and reply to the subsequent questions as desired. Inside the directory "PROC_nmssm_x", this procedure can be repeated with different cards. Good luck! I11=1: Landau poles and false minima constraints are imposed (default value) I11=0: They are not imposed (for scanning versions only, I2=/=0) I12=1: LEP/Tevatron Higgs+sparticle constraints are imposed (default value) I12=0: They are not imposed (for scanning versions only, I2=/=0) I13=1: LHC Higgs constraints are imposed (default value) I13=0: They are not imposed (for scanning versions only, I2=/=0) I14=1: Upsilon, B and K decay constraints are imposed (default value) I14=0: They are not imposed (for scanning versions only, I2=/=0) I15=1: EDM constraints are imposed (default value) I15=0: They are not imposed (for scanning versions only, I2=/=0, only for CPV models, I1=2) I16=1: CMS charg(neutral)ino constraints from arXiv:1801.03957[hep-ex] are imposed (default value) I16=0: They are not imposed (for scanning versions only, I2=/=0, not for GMSB models, I0=/=2) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% INPUT/OUTPUT parameters *************************************************************************************** CP conserving case *************************************************************************************** Content of the array PAR (parameters at the SUSY scale): PAR(1) = lambda PAR(2) = kappa PAR(3) = tan(beta) PAR(4) = mu (effective mu term = lambda*s) PAR(5) = Alambda (if MA is not an input) 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 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 Output parameters: 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) AMASS(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 CB(I) Reduced couplings of h1,h2,h3 (i=1,2,3) or a1,a2 (i=4,5) to b-quarks including DELMB corrections 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 *************************************************************************************** CP violating case *************************************************************************************** Content of the arrays REALP, IMAGP, PAR, MYPHASES (parameters at the SUSY scale): REALP(1) = Re(lambda) REALP(2) = Re(kappa) REALP(3) = Re(M1) REALP(4) = Re(M2) REALP(5) = Re(M3) REALP(6) = Re(AU3) REALP(7) = Re(AD3) REALP(8) = Re(AE3) REALP(9) = Re(XIF) REALP(10) = Re(XIS) REALP(11) = Re(MUP) REALP(12) = Re(MSP) REALP(13) = Re(M3H) REALP(14) = Re(MUEFF) IMAGP(1) = Im(lambda) IMAGP(2) = Im(kappa) IMAGP(3) = Im(M1) IMAGP(4) = Im(M2) IMAGP(5) = Im(M3) IMAGP(6) = Im(AU3) IMAGP(7) = Im(AD3) IMAGP(8) = Im(AE3) IMAGP(9) = Im(XIF) IMAGP(10) = Im(XIS) IMAGP(11) = Im(MUP) IMAGP(12) = Im(MSP) IMAGP(13) = Im(M3H) [ IMAGP(14) = Im(MUEFF) -> aligned with lambda = REALP(14)*IMAGP(1)/REALP(1) ] [ PAR(1) = |lambda| -> Sqrt( REALP(1)**2 + IMAGP(1)**2 ) ] [ PAR(2) = sgn(Re(kappa)).|kappa| -> sgn(REALP(2)).Sqrt( REALP(2)**2 + IMAGP(2)**2 ) ] PAR(3) = tan(beta) [ PAR(4) = sgn(Re(mueff)).|mueff| -> aligned with lambda = REALP(14)*PAR(1)/REALP(1) ] PAR(5) = Re(Alambda) (if (MA,XIF) is not an input) PAR(6) = Re(Akappa) (if (MP,XIS) is not an input) PAR(7) = mQ3**2 PAR(8) = mU3**2 PAR(9) = mD3**2 PAR(10) = mL3**2 PAR(11) = mE3**2 PAR(12) = AU3 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) = sgn(Re(M1)).|M1| -> sgn(REALP(3)).Sqrt( REALP(3)**2 + IMAGP(3)**2 ) ] [ PAR(21) = sgn(Re(M2)).|M2| -> sgn(REALP(4)).Sqrt( REALP(4)**2 + IMAGP(4)**2 ) ] [ PAR(22) = sgn(Re(M3)).|M3| -> sgn(REALP(5)).Sqrt( REALP(5)**2 + IMAGP(5)**2 ) ] PAR(23) = MA (diagonal doublet CP-odd mass matrix element at tree-level) PAR(24) = MP (diagonal singlet CP-odd mass matrix element at tree-level) PAR(25) = AE2 Phases in [-Pi/2,Pi/2] [ MYPHASES(1) = Phi_lambda -> arctan( IMAGP(1) / REALP(1) ) ] [ MYPHASES(2) = Phi_kappa -> arctan( IMAGP(2) / REALP(2) ) ] [ MYPHASES(3) = Phi_M1 -> arctan( IMAGP(3) / REALP(3) ) ] [ MYPHASES(4) = Phi_M2 -> arctan( IMAGP(4) / REALP(4) ) ] [ MYPHASES(5) = Phi_M3 -> arctan( IMAGP(5) / REALP(5) ) ] [ MYPHASES(6) = Phi_AU3 -> arctan( IMAGP(6) / REALP(6) ) ] [ MYPHASES(7) = Phi_AD3 -> arctan( IMAGP(7) / REALP(7) ) ] [ MYPHASES(8) = Phi_AE3 -> arctan( IMAGP(8) / REALP(8) ) ] [ MYPHASES(9) = Phi_AU2 (unused) -> 0 ] [ MYPHASES(10) = Phi_AD2 (unused) -> 0 ] [ MYPHASES(11) = Phi_AE2 (unused) -> 0 ] [ MYPHASES(12) = Phi_XIS -> arctan( IMAGP(10) / REALP(10) ) ] [ MYPHASES(13) = Phi_MSP -> arctan( IMAGP(12) / REALP(12) ) ] [ MYPHASES(14) = Phi_XIF -> arctan( IMAGP(9) / REALP(9) ) ] [ MYPHASES(15) = Phi_MUP -> arctan( IMAGP(11) / REALP(11) ) ] [ MYPHASES(16) = Phi_M3H -> arctan( IMAGP(13) / REALP(13) ) ] Output parameters: MH0(1-5): Higgs masses (squared, ordered) XH0(1-5,1-5): Mixing matrix in the basis H1_R=Hu_R, H2_R=Hd_R, S_R, A, S_I where the Goldstone mode has been rotated away MHC: Charged Higgs mass (squared) CU,CD,CV,CJ,CG(i),CZG(i) Reduced scalar couplings of hi (i=1..5) to up type fermions, down type fermions, gauge bosons, gluons, photons and (Z+photon) CUP,CDP,CJP,CGP(i),CZGP(i) idem for pseudoscalar component CB(i) Reduced couplings of hi (i=1..5) to b-quarks including DELMB corrections CBP(i) idem for pseudoscalar component WIDTH(i) Total decay width of hi (i=1..5) with the following branching ratios: BRJJ(i) hi (i=1..5) -> gluon gluon BREE(i) " -> e+ e- BRMM(i) " -> mu mu BRLL(i) " -> tau tau BRSS(i) " -> ss BRCC(i) " -> cc BRBB(i) " -> bb BRTT(i) " -> tt BRWW(i) " -> WW BRZZ(i) " -> ZZ BRGG(i) " -> gamma gamma BRZG(i) " -> Z gamma BRHIGGS(i) " -> other Higgses, including: BRHCHC(i) " -> h+h- BRHAZ(i,j) " -> Zhj (j=1..4) BRHCW(i) " -> h+W- BRHHH(i,j) " -> h1h1 (1), h1h2 (2), h2h2 (3), h1h3 (4), h2h3 (5), h3h3 (6), h1h4 (7), h2h4 (8), h3h4 (9), h4h4 (10) BRSUSY(i) " -> susy particles, including: BRNEU(i,j,k) -> neutralinos j,k (j,k=1..5) BRCHA(i,j) -> charginos 11, 12, 22 (j=1..3) BRHSQ(i,j) " -> uLuL, uRuR, dLdL, dRdR, t1t1, t2t2, t1t2, b1b1, b2b2, b1b2 (j=1..10) BRHSL(i,j) " -> lLlL, lRlR, nLnL, l1l1, l2l2, l1l2, ntnt (i=1..3, j=1..7) 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) " -> hiW+ (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, squared, ordered in mass) NEU(i,j,k) chi_i components of bino, wino, higgsino u&d, singlino (i,j=1..5), real part (k=1) and imaginary part (k=2) MCH2(i) Chargino masses (squared) U(i,j,k),V(i,j,k) Chargino mixing matrices, real part (k=1) and imaginary part (k=2) MGl Gluino mass Sfermion masses: MSU2(i) S-up/charm, masses squared (i=1->L,2->R), DRbar MSU2P(i) idem, including QCD+Yuk. corrections MSD2(i) S-down/strange, masses squared (i=1->L,2->R), DRbar MSD2P(i) idem, including QCD+Yuk. corrections MSE2(i) S-electron, masses squared (i=1->L,2->R), DRbar MSMU2(i) S-muon, masses squared ordered (i=1,2), DRbar UMU(i,j,k) S-muon rotation matrix (i->mass, j=1->L, 2->R)* MSNE2 Sneutrino 1st gen., mass squared MST2(i) S-top, masses squared ordered (i=1->L,2->R), DRbar MST2P(i) idem, including QCD+Yuk. corrections UT(i,j,k) S-top rotation matrix (i->mass, j=1->L, 2->R)* MSB2(i) S-bottom, masses squared ordered (i=1->L,2->R), DRbar MSB2P(i) idem, including QCD+Yuk. corrections UB(i,j,k) S-bottom rotation matrix (i->mass, j=1->L, 2->R)* MSL2(i) S-tau, masses squared ordered (i=1,2), DRbar UL(i,j,k) S-tau rotation matrix (i->mass, j=1->L, 2->R)* MSNT2 Sneutrino 3rd gen., mass squared * [real part (k=1) and imaginary part (k=2)] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Content of the array PROB (phenomenological and theoretical constraints): PROB(I) = 0, I = 1..62: 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 or < 511 keV 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 DM relic density (checked only if OMGFLAG=/=0) PROB(31) =/= 0 excluded by DM SI WIMP-nucleon xs (checked if |OMGFLAG|=2 or 4) 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 measurement PROB(40) =/= 0 BR(B-->X_s mu+ mu-) more than 2 sigma away PROB(41) =/= 0 excluded by ee -> hZ, h -> AA -> 4taus (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 t -> bH+ (LHC) PROB(46) =/= 0 No Higgs in the MHmin-MHmax GeV range PROB(47) =/= 0 chi2gam > chi2max PROB(48) =/= 0 chi2bb > chi2max PROB(49) =/= 0 chi2zz > chi2max PROB(51) =/= 0: excluded by H/A->tautau PROB(52) =/= 0: Excluded by H->AA->4leptons/2lept.+2b (LHC) PROB(53) =/= 0: excluded by ggF->H/A->gamgam (65GeV < M < 122GeV, ATLAS) PROB(55) =/= 0: b -> d gamma more than 2 sigma away PROB(56) =/= 0: B_d -> mu+ mu- more than 2 sigma away PROB(57) =/= 0: b -> s nu nubar more than 2 sigma away PROB(58) =/= 0: b -> c tau nu more than 2 sigma away (as SM) PROB(59) =/= 0: K -> pi nu nubar more than 2 sigma away PROB(60) =/= 0: DMK / epsK more than 2 sigma away PROB(61) =/= 0 excluded by DM SD WIMP-neutron xs (checked if |OMGFLAG|=2 or 4) PROB(62) =/= 0 excluded by DM SD WIMP-proton xs (checked if |OMGFLAG|=2 or 4) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% INPUT OPTIONS (see the examples in the directory "SAMPLES"): The definitions of the parameters XIF, XIS, MU', MS'^2, M3H^2 in the general NMSSM (GNMSSM) can be found in Phys.Rept. 496 (2010) 1 (arXiv:0910.1785). *************************************************************************************** INPUT OPTIONS FOR NMHDECAY: *************************************************************************************** Sparticles: M2, AU3, AD3, AE3, ML3, ME3, MQ3, MU3, MD3 have to be given. Unless specified otherwise, the code assumes M1=M2/2, M3=3*M2, AE1=AE2=AE3, ML1=ML2=ML3, ME1=ME2=ME3, MQ1=MQ2=MQ3, MU1=MU2=MU3, MD1=MD2=MD3 NMSSM-specific parameters in the Higgs sector: The soft Higgs mass terms MHU^2, MHD^2 and MS^2 are always computed in terms of MZ, TANB and MUEFF. Always to specify: TANB, LAMBDA and MUEFF Note: values of commented parameters (including KAPPA) are assumed = 0, unless they are computed as output. Options for the parameter set "ALAMBDA, MA, XIF" where MA=diagonal doublet-doublet entry in the CP-odd mass matrix: Two out of these three parameters have to be specified, the remaining one is computed. (If all are commented, MA is computed assuming ALAMBDA=XIF=0, if only MA is input, ALAMBDA is computed assuming XIF=0, if only ALAMBDA is input, MA is computed assuming XIF=0, if only XIF is input, MA is computed assuming ALMBDA=0) Options for the parameter set "AKAPPA, MP, XIS" where MP=diagonal singlet-singlet entry in the CP-odd mass matrix: Two out of these three parameters have to be specified, the remaining one is computed. (If all are commented, MP is computed assuming AKAPPA=XIS=0, if only MP is input, AKAPPA is computed assuming XIS=0, if only AKAPPA is input, MP is computed assuming XIS=0, if only XIS is input, MP is computed assuming AKAPPA=0) Note: If KAPPA=0, AKAPPA must be =0 In the output file, the values of the computed parameters are shown in the BLOCK EXTPAR, but commented. PARTICULAR CASES: Z_3-invariant NMSSM: Always to specify: TANB, LAMBDA, MUEFF, KAPPA, and in addition - either ALAMBDA or MA, and - either AKAPPA or MP. If too many of these latter parameters are specified, XIF and/or XIS are computed and indicated in the output. nMSSM (kappa=Akappa=0, but XIF, XIS non-vanishing): Always to specify: TANB, LAMBDA, MUEFF, and in addition - two parameters of the set (XIF, ALAMBDA, MA), the remaining parameter is computed, and - MP or XIS, the remaining parameter is computed. *************************************************************************************** INPUT OPTIONS FOR NMHDECAY_RAND: *************************************************************************************** Scans can be performed over all independent parameters of the Higgs sector (see above), and the gaugino masses M1, M2 and M3. If M1 and/or M3 are not specified, M1=M2/2 and M3=3*M2. For each parameter "P", P_max=P_min if P_max is not specified. P_min=0 by default, except for TANB, LAMBDA and MUEFF which have to be specified, and except for the parameters which are computed. Unless P_min=0 by default, P_max cannot be given if P_min is not known. (If KAPPA=0, AKAPPA is set = 0.) The user has to specify the total number of points. *************************************************************************************** INPUT OPTIONS FOR NMHDECAY_GRID: *************************************************************************************** The choice of independent parameters is as above. The user has to specify the number of points "NP" for each parameter to be scanned. By default NP=1. If P_max > P_min and NP=1, P_min is used. If P_max=P_min and NP>1, the user is alerted by a warning. *************************************************************************************** INPUT OPTIONS FOR NMHDECAY_MCMC: *************************************************************************************** For MCMC scans, the user has to provide a starting point with parameters as for NMHDECAY. For each parameter P, the user has to give the relative and the minimum step size PDEV and PMIN (resp.). Then, at each step, the new value of P is computed as: PNEW=P+MAX(DABS(P),PMIN)*PDEV*GAU where GAU is a random number with a gaussian distribution (mean=0, width=1). Whether the new point is taken as the new starting point for next step of the random walk or not depends on the value of a penalty function which is a sum of all the discrepancies of the chosen point with respect to experimental bounds. The total number of scanned points is given in BLOCK STEPS as NTOT and the initial random seed as ISEED. *************************************************************************************** INPUT OPTIONS FOR NMSPEC: *************************************************************************************** The user has to specify TANB, LAMBDA, M0, M12 (unless all M1, M2 and M3 are specified separately) and A0. ALAMBDA and AKAPPA (both at MGUT) can be specified independently from A0, if desired. MHD^2 and MHU^2 (both at MGUT) can be specified independently from m0, if MUEFF is NOT specified. Default values for all other parameters (unless computed in terms of others) are =0. By default, KAPPA is computed in terms of the other parameters. If KAPPA is input: XIF is computed and hence =/= 0 in general. If neither KAPPA nor XIF are given as inputs then KAPPA is computed assuming XIF=0. Both parameters cannot be simultaneously inputs. By default, the singlet soft mass MS^2 at MGUT is computed in terms of the other parameters and hence =/= m0 in general. If MS^2 is input, XIS is computed instead and hence =/= 0 in general. If neither MS^2 nor XIS are given as inputs then MS^2 is computed assuming XIS=0. Both parameters cannot be simultaneously inputs. By default, MUEFF is computed in terms of the other parameters, but then its sign SIGMU (FLOAT, not INTEGER) must be given. If MUEFF is an input parameter, KAPPA, XIF and XIS are also interpreted as input parameter (with default value =0 unless specified otherwise). Then the Higgs soft masses MHD^2, MHU^2 and MS^2 at MGUT are computed in terms of the other parameters and hence =/= m0 in general. By default, XIF=XIS=MU'=MS'^2=M3H^2=0, but can be specified otherwise (at MGUT). Exceptions: If MUEFF is NOT an input parameter and - if KAPPA is input: XIF is computed and hence =/= 0 in general; - if MS^2 is input: XIS is computed and hence =/= 0 in general. PARTICULAR CASES: cNMSSM (fully constrained Z_3-invariant NMSSM): Not possible as input; at least MS^2 at MGUT is an output sNMSSM (semi-constrained Z_3-invariant NMSSM, with non-universal singlet sector): Inputs TANB, LAMBDA, M0, M12, A0, ALAMBDA, AKAPPA -> KAPPA, MUEFF and MS^2 are computed. NUH-NMSSM (Z_3-invariant NMSSM with non-universal Higgs masses MS^2, MHD^2 and MHU^2 at MGUT): Inputs TANB, LAMBDA, KAPPA, M0, M12, A0 (possibly also ALAMBDA, AKAPPA), MUEFF -> MS^2, MHD^2 and MHU^2 are computed. nMSSM (no self coupling for the singlet): Inputs TANB, LAMBDA, KAPPA=AKAPPA=0, M0, M12, A0 (possibly also ALAMBDA), XIS -> XIF, MUEFF, MS^2 are computed; or Inputs TANB, LAMBDA, KAPPA=AKAPPA=0, M0, M12, A0 (possibly also ALAMBDA), MS^2 -> XIF, MUEFF, XIS are computed. Inputs TANB, LAMBDA, KAPPA=AKAPPA=0, M0, M12, A0 (possibly also ALAMBDA), MUEFF, XIF, XIS -> MHU^2, MHD^2, MS^2 are computed. *************************************************************************************** INPUT OPTIONS FOR NMSPEC_RAND: *************************************************************************************** Scans can be performed over all independent parameters (as above). For each parameter "P", P_max=P_min if P_max is not specified, P_max cannot be given if P_min is not known. Same default value assumed for P_min as for P. The user has to specify the total number of points. *************************************************************************************** INPUT OPTIONS FOR NMSPEC_GRID: *************************************************************************************** The choice of independent parameters is as above. The user has to specify the number of points "NP" for each parameter to be scanned. By default NP=1. If P_max > P_min and NP=1, P_min is used. If P_max=P_min and NP>1, the user is alerted by a warning. *************************************************************************************** INPUT OPTIONS FOR NMSPEC_MCMC: *************************************************************************************** For MCMC scans, the user has to provide a starting point with parameters as for NMSPEC. For each parameter P, the user has to give the relative and the minimum step size PDEV and PMIN (resp.). Then, at each step, the new value of P is computed as: PNEW=P+MAX(DABS(P),PMIN)*PDEV*GAU where GAU is a random number with a gaussian distribution (mean=0, width=1). Whether the new point is taken as the new starting point for next step of the random walk or not depends on the value of a penalty function which is a sum of all the discrepancies of the chosen point with respect to experimental bounds. The total number of scanned points is given in BLOCK STEPS as NTOT and the initial random seed as ISEED. *************************************************************************************** INPUT OPTIONS FOR NMGMSB: *************************************************************************************** If NMGMSB = 0: The definitions of the parameters specific to GMSB models can be found in JHEP 0805 (2008) 044 (arXiv:0803.2962). The user has to specify TANB, LAMBDA, MSUSY (SUSY scale), MMESS (messenger scale) N5 (number of 5-plets) and SIGMU (sign of MUEFF which is always computed). ALAMBDA (at MMESS) can be specified if =/=0. AKAPPA is always =3*ALAMBDA. Default values for all other parameters (unless computed in terms of others) are 0. By default, KAPPA is computed in terms of the other parameters. If KAPPA is input: XIF is computed and hence =/= 0 in general. If neither KAPPA nor XIF are given as inputs then KAPPA is computed assuming XIF=0. Both parameters cannot be simultaneously inputs. By default, the singlet soft mass MS^2 at MGUT is computed in terms of the other parameters. If MS^2 is input, XIS is computed instead and hence =/= 0 in general. If neither MS^2 nor XIS are given as inputs then MS^2 is computed assuming XIS=0. Both parameters cannot be simultaneously inputs. By default, XIF=XIS=MU'=MS'^2=0, but can be specified otherwise (at MMESS). - if KAPPA is input: XIF is computed and hence =/= 0 in general; - if MS^2 is input: XIS is computed and hence =/= 0 in general. If NMGMSB =/= 0: The definitions of the parameters specific to GMSB models above MMESS can be found in arXiv:0706.3873[hep-ph] and arXiv:1107.3006[hep-ph]. The user has to specify TANB, LAMBDA, MSUSY (SUSY scale), MMESS (messenger scale) N5 (number of 5-plets) and SIGMU (sign of MUEFF which is always computed from minimisation eqs.). The extra parameters XIU (at MGUT), LTT, LTT, LU, LD, LT, LB, LL (at MMESS) can be given in BLOCK EXTPAR (switches 73-80). Soft parameters at MMESS are related to these extra couplings. MUEFF, KAPPA and MS are computed from minimisation eqs. DMIN, the maximal relative deviation allowed between the computed and expected value for MS, has to be given in BLOCK EXTPAR (switch 0). *************************************************************************************** INPUT OPTIONS FOR NMGMSB_RAND: *************************************************************************************** Scans can be performed over all independent parameters (as above). For each parameter "P", P_max=P_min if P_max is not specified, P_max cannot be given if P_min is not known. Same default value assumed for P_min as for P. The user has to specify the total number of points. *************************************************************************************** INPUT OPTIONS FOR NMGMSB_GRID: *************************************************************************************** The choice of independent parameters is as above. The user has to specify the number of points "NP" for each parameter to be scanned. By default NP=1. If P_max > P_min and NP=1, P_min is used. If P_max=P_min and NP>1, the user is alerted by a warning. *************************************************************************************** INPUT OPTIONS FOR NMGMSB_MCMC: *************************************************************************************** For MCMC scans, the user has to provide a starting point with parameters as for NMGMSB. For each parameter P, the user has to give the relative and the minimum step size PDEV and PMIN (resp.). Then, at each step, the new value of P is computed as: PNEW=P+MAX(DABS(P),PMIN)*PDEV*GAU where GAU is a random number with a gaussian distribution (mean=0, width=1). Whether the new point is taken as the new starting point for next step of the random walk or not depends on the value of a penalty function which is a sum of all the discrepancies of the chosen point with respect to experimental bounds. The total number of scanned points is given in BLOCK STEPS as NTOT and the initial random seed as ISEED. *************************************************************************************** INPUT OPTIONS FOR NMHDECAYCPV: *************************************************************************************** Same input parameters as NMHDECAY for the real part of the parameters at SUSY scale. Imaginary parts of the parameters are given in a separate block except Im(MUEFF): aligned with LAMBDA Im(ALAMBDA), Im(AKAPPA): computed from minimisation equations If KAPPA=0 (nMSSM) then AKAPPA must be zero and Im(XIS) is computed from the second minimisation equation (hence in general =/=0). For the scanning versions (RAND, GRID and MCMC) the input parameters are as above.