RunScans = {   }; 
 
(*----------------------------------------------*) 
(*   SETUP                                      *) 
(*----------------------------------------------*) 

(* Define location, in- and output of SPheno *) 
 
LoadSettings="DefaultSettings.m.NMSSM";  
 
 
DEFINITION[a_][IncludeMicrOmegas]=False; 
 
(* ------------------------- CALCHEP -------------------------*) 
 
DEFINITION[a_][IncludeCalcHep]=True; 


DEFINITION[a_][CalcHepRuns]={{"~/Desktop/Temp/toolbox-1.1.7/build/CH_NMSSM/eeZH","./n_calchep -blind \"[[[[[[[[[{[[[{}}\"","prt_1"} }
 
DEFINITION[a_][IncludeHiggsBounds]=True;

DEFINITION[a_][IncludeMicrOmegas]=True;
DEFINITION[a_][DarkMatterCandidate]=1000022;

 
RunScans = {NMSSMdemo}; 
 
(*----------------------------------------------*)
(*    PARAMETER RANGES                          *)
(*----------------------------------------------*)

DEFINITION[a_][Blocks]={MODSEL,SMINPUTS,SPhenoInput,MINPAR,EXTPAR}; 
 
DEFINITION[a_][MODSEL]={ 
{{1},{Value->1}}, 
{{2},{Value->1}}, 
{{6},{Value->1}} 
}; 
DEFINITION[a_][SMINPUTS]={ 
{{2},{Value->1.166390*10^-5}}, 
{{3},{Value->0.1172}}, 
{{4},{Value->91.18760}}, 
{{5},{Value->4.2}}, 
{{6},{Value->172.9}}, 
{{7},{Value->1.777}} 
}; 
DEFINITION[a_][MINPAR]={ 
{{1},{Min->400, Max->2500, Steps->7, Distribution->LINEAR}} (*m0*), 
{{2},{Min->400, Max->2500, Steps->7, Distribution->LINEAR}} (*m12*), 
{{3},{Value->10}} (*TanBeta*), 
{{5},{Value->-2500}} (*Azero*)}; 
 
DEFINITION[a_][EXTPAR]={ 
{{61},{Value -> 0.09}} (*LambdaInput*), 
{{62},{Value -> 0.11}} (*KappaInput*), 
{{63},{Value->-1500}} (*ALambdaInput*), 
{{64},{Value->-36}} (*AKappaInput*), 
{{65},{Value->13689}} (*vSInput*)}; 

DEFINITION[a_][SPhenoInput]={ 
{{1},{Value->-1}}, (* error level *)
{{2},{Value->1}},  (* SPA conventions *) 
{{11},{Value->1}}, (* Calculate widhts and BRs *)
{{12},{Value->0.0001}}, (* minimal BR to write out *) 
{{13},{Value->0}}, (* Enable 3-body decays *) 
{{34},{Value->0.0001}}, (* precision of masses *) 
{{38},{Value->2}}, (* 1/2 - Loop RGEs *) 
{{50},{Value->0}}, (* Majorana convention *) 
{{51},{Value->0}}, (* Switch to CKM matrix *) 
{{55},{Value->1}}, (* 1 - Loop masses *)
{{57},{Value->1}}, (* low energy constraints *)
{{60},{Value->1}}, (* Include possible, kinetic mixing *)
{{65},{Value->1}}, (* Solution tadpole equation *)
{{75},{Value->1}}, (* Write WHIZARD files *) 
{{76},{Value->1}},  (* Write HiggsBounds files *) 
{{86},{Value->0.}}  (* Maximal width to be counted as invisible in Higgs decays; -1: only LSP *) 
}; 
(*----------------------------------------------*)
(*    PLOTS                                     *)
(*----------------------------------------------*) 

DEFINITION[a_][Plots]={ 
{P3D, {MINPAR[1],MINPAR[2],MASS[25]}, {Frame->True, Axes->False,FrameLabel->{UseLaTeX["$m_0$"],UseLaTeX["$M_{1/2}$"]},ContourLabels->True},"m0_m12_h1.eps"},
{P3D, {MINPAR[1],MINPAR[2],Log[10,NMHMIX[1,3]^2]}, {Frame->True, Axes->False,FrameLabel->{UseLaTeX["$m_0$"],UseLaTeX["$M_{1/2}$"]},ContourLabels->True},"m0_m12_SingletFraction.eps"},
{P2D, {MINPAR[1],{CHEP[1]}},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$m_0$"],UseLaTeX["$e^+ e^- \\to Z h_1$"]}}},"m0_eeZH.eps"},
{P2D, {MINPAR[2],{CHEP[1]}},{{Frame->True, Axes->False, FrameLabel->{seLaTeX["$M_{1/2}$"],UseLaTeX["$e^+ e^- \\to Z h_1$"]}}},"m12_eeZH.eps"},
{P2D, {MINPAR[1],{DAKRMATTER[1]}},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$m_0$"],UseLaTeX["$\\Omega h^2$"]}}},"m0_DM.eps"},
{P2D, {MINPAR[2],MASS[1000022]},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$M_{1/2}$"],UseLaTeX["$M_{\\tilde{\\chi}^0_{1}}$"]}}},"m12_mN1.eps"},
{P2D, {MASS[1000022],{Log[10,DARKMATTER[1]]}},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$M_{\\tilde{\\chi}^0_1}$"],UseLaTeX["$\\Omega h^2$"]}}},"mN1_DM.eps"},
{P2D, {Abs[MASS[1000022]-MASS[1000011]],{Log[10,DARKMATTER[1]]}},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$M_{\\tilde{\\chi}^0_1}-m_{\\tau_1}$"],UseLaTeX["$\\Omega h^2$"]}}},"DeltaM_DM.eps"},
{P2D, {MASS[25],{HIGGSBOUNDS[7]}},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$m_{h_1}$"],"allowed"}}},"m25_HB-bounds.eps"},
{P2D, {MASS[25],{HIGGSBOUNDS[9]}},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$m_{h_1}$"],"obsratio"}}},"m25_HB-saturation.eps"},
{P2D, {MASS[25],{25[FinalParticles[22, 22]]*10^3}},{{Frame->True, Axes->False, FrameLabel->{UseLaTeX["$m_{h_1}$"],UseLaTeX["$10^{3} BR(h_1\\to \\gamma\\gamma)$"]}}},"mh1_BR_hGammaGamma.eps"}
}; 
