softsusy is hosted by Hepforge, IPPP Durham
SOFTSUSY  4.1
softsusy.h
Go to the documentation of this file.
1 
15 #ifndef SOFTSUSY_H
16 #define SOFTSUSY_H
17 
18 #include <iostream>
19 #include <fstream>
20 #include <sstream>
21 #include <cstring>
22 #include <cstdlib>
23 #include <cmath>
24 #include <limits>
25 #include "def.h"
26 #include "utils.h"
27 #include "numerics.h"
28 #include "physpars.h"
29 #include "lowe.h"
30 #include "softpars.h"
31 #include "twoloophiggs.h"
32 #include "mssmUtils.h"
33 #include "higher_order.h"
34 #include "mssmjacobian.h"
35 
36 #define HR "----------------------------------------------------------"
37 
38 namespace softsusy {
39 
42  DeltaQpole = 1,
43  DeltaQmatch = 2,
44  DeltaMt = 4,
45  DeltaAlphaS = 8,
46  DeltaAlphaEm = 16,
47  DeltaAll = DeltaQpole + DeltaQmatch + DeltaMt + DeltaAlphaS + DeltaAlphaEm
48  };
49 
51  const double mxDefault = 1.9e16;
52 
54  class AltEwsbMssm {
55  private:
58  double mAcond, muCond;
59  public:
60  void setAltEwsb(double ma, double mu) { mAcond = ma; muCond = mu; };
61  void setAltEwsbMssm(const AltEwsbMssm & s) { *this = s; };
62 
63  // int displayConditionStyle() { return conditionStyle; };
64  double displayMaCond() const { return mAcond; };
65  double displayMuCond() const { return muCond; };
66  const AltEwsbMssm & displayAltEwsbMssm() const { return *this; };
67 
68  void setMaCond(double maInput) { mAcond = maInput; };
69  void setMuCond(double muInput) { muCond = muInput; };
70  };
71 
73  class MssmSoftsusy: public MssmSusy, public MssmSoftPars, public AltEwsbMssm,
74  public Approx, public RGE {
84  private:
85  sPhysical physpars;
86  drBarPars forLoops;
87  sProblem problem;
88  double msusy;
89  double minV;
90  double mw;
91  QedQcd dataSet;
92  double fracDiff;
93  bool setTbAtMX;
95  bool altEwsb;
96  bool altMt;
97  bool altAlphaS;
98  bool altAlphaEm;
99  double predMzSq;
100  double t1OV1Ms, t2OV2Ms;
101  double t1OV1Ms1loop, t2OV2Ms1loop;
106  double qewsb;
107 
108  protected:
109  void setT1OV1Ms(double t1) { t1OV1Ms = t1; }
110  void setT2OV2Ms(double t2) { t2OV2Ms = t2; }
111  void setT1OV1Ms1loop(double t1) { t1OV1Ms1loop = t1; }
112  void setT2OV2Ms1loop(double t2) { t2OV2Ms1loop = t2; }
113  double mxBC;
114 
115  void check_flags();
116  DoubleMatrix calcHiggs3L(bool is_bottom);
117  public:
118  typedef void (*TMSSMBoundaryCondition)(MssmSoftsusy&, const DoubleVector&);
119 
122 
123  // void (*boundaryCondition)(Softsusy &, const DoubleVector &);
125  MssmSoftsusy();
127  MssmSoftsusy(const MssmSusyRGE &);
129 
130  MssmSoftsusy(const MssmSoftsusy &);
135  MssmSoftsusy(const MssmSusy &, const MssmSoftPars & s,
136  const sPhysical & sp, double mu, int l,
137  int t, double hv);
138 
140  const MssmSoftsusy & operator=(const MssmSoftsusy & s);
141 
143  virtual const DoubleVector display() const;
144 
146  const MssmSoftsusy & displayMssmSoft() const;
147 
149  const sPhysical & displayPhys() const;
150 
153  TMSSMBoundaryCondition boundaryCondition,
154  double mxGuess,
155  const DoubleVector & pars, int sgnMu, double tanb,
156  const QedQcd & oneset, bool gaugeUnification,
157  bool ewsbBCscale = false,
158  int contributions = DeltaAll) const;
161 
164  const drBarPars & displayDrBarPars() const;
166  const sProblem & displayProblem() const {return problem; };
168  const QedQcd & displayDataSet() const;
170  double displaySoftA(trilinears, int, int) const;
171 
173  double displayFracDiff() const { return fracDiff; };
174  double displayMinpot() const;
175  double displayMsusy() const;
176  double displayMw() const;
178  double displayMwRun() const;
180  double displayMzRun() const;
181  double displayTadpole1Ms() const;
182  double displayTadpole2Ms() const;
183  double displayTadpole1Ms1loop() const;
184  double displayTadpole2Ms1loop() const;
185  double displayMxBC() const { return mxBC; };
187  const MssmSoftsusy & displaySoftsusy() const { return *this; }
189  double displayMz() const { return double(MZ); }
191  bool displaySetTbAtMX() const { return setTbAtMX; }
192  bool displayAltEwsb() const { return altEwsb; }
193  bool displayAltMt() const { return altMt; }
194  bool displayAltAlphaS() const { return altAlphaS; }
195  bool displayAltAlphaEm() const { return altAlphaEm; }
196  double displayPredMzSq() const { return predMzSq; }
197  double displayQewsb() const { return qewsb; }
198  double displayMatchingScale() const { return displayDataSet().displayMu(); }
199 
201  void flagMgutOutOfBounds(bool a) { problem.mgutOutOfBounds = a; };
203  void flagIrqfp(bool a) { problem.irqfp = a; };
205  void flagNonperturbative(bool a) { problem.nonperturbative = a; };
207  void flagTachyon(tachyonType a) {
210  if (close(displayMu(), MZ, 1.0e-6)) problem.tachyonWarning = a;
211  else { problem.tachyon = a;
212  if (PRINTOUT > 2) cout << tachyonNames[a] << " tachyon ";
213  }
214  };
215  void flagTachyonWarning(tachyonType a) { problem.tachyonWarning = a; }
217  void flagM3sq(bool a) { problem.m3sq = a; };
219  void flagNoConvergence(bool a) { problem.noConvergence = a; };
221  void flagNoMuConvergence(bool a) { problem.noMuConvergence = a; };
223  void flagNoRhoConvergence(bool a) { problem.noRhoConvergence = a; };
225  void flagMusqwrongsign(bool a) { problem.muSqWrongSign = a; };
228  void flagInaccurateHiggsMass(bool a) { problem.inaccurateHiggsMass = a; };
230  void flagHiggsufb(bool a) { problem.higgsUfb = a; };
232  void flagNotGlobalMin(bool a) { problem.notGlobalMin = a; };
234  void flagAllProblems(bool a, tachyonType b) { problem.irqfp = a;
235  problem.tachyon = b; problem.tachyonWarning = b;
236  problem.m3sq = a; problem.badConvergence = a;
237  problem.noConvergence = a; problem.higgsUfb = a; problem.notGlobalMin = a;
238  problem.nonperturbative = a; problem.noRhoConvergence = a;
239  problem.noMuConvergence = a; problem.muSqWrongSign = a;
240  problem.inaccurateHiggsMass = b; problem.mgutOutOfBounds = a; }
242  void flagProblemThrown(bool a) { problem.problemThrown = a; }
244  void setProblem(const sProblem a) { problem = a; }
245 
247  void setSoftsusy(const MssmSoftsusy & s) { *this = s; };
249  void setData(const QedQcd & r) { dataSet = r; };
251  void setMinpot(double);
254  void setMsusy(double);
256  void setMw(double);
258  void setPhys(const sPhysical & s) { physpars = s; };
260  void setMxBC(double mx) { mxBC = mx; };
262  void setDrBarPars(const drBarPars & s) { forLoops = s; };
264  void setSetTbAtMX(bool a) { setTbAtMX = a; }
266  void useAlternativeEwsb() { altEwsb = true; }
268  void useAlternativeMt() { altMt = true; }
270  void useAlternativeAlphaS() { altAlphaS = true; }
272  void useAlternativeAlphaEm() { altAlphaEm = true; }
274  void setPredMzSq(double a) { predMzSq = a; }
276  void set(const DoubleVector &);
278  void setQewsb(double q) { qewsb = q; };
279 
282  void setTwoLoopAlphasThresholds(bool sw);
284  void setTwoLoopMtThresholds(bool sw);
290  void setTwoLoopMbMtauThresholds(bool sw);
292  void setAllTwoLoopThresholds(bool sw);
293 
295  DoubleVector beta() const {
296  sBrevity a;
299  int i; for (i=1; i<=y.displayEnd(); i++) x(i) = y(i);
300  return x;
301  };
302 
304  void setFracDiff(double fD) { fracDiff = fD; };
306  //for use in doCalcTadpole1oneLoop
307  void H1SfSfCouplings(DoubleMatrix & lTS1Lr, DoubleMatrix & lBS1Lr,
308  DoubleMatrix & lTauS1Lr, double gmzOcthW, double mu,
309  double cosb, double v1) const;
311  //for use in doCalcTadpole1oneLoop
312  void H2SfSfCouplings(DoubleMatrix & lTS2Lr, DoubleMatrix & lBS2Lr,
313  DoubleMatrix & lTauS2Lr, double gmzOcthW, double mu,
314  double sinb) const;
316  double doCalcTad1Sfermions(DoubleMatrix lTS1Lr, DoubleMatrix lBS1Lr,
317  DoubleMatrix lTauS1Lr, double costhDRbar) const;
319  double doCalcTad2Sfermions(DoubleMatrix lTS2Lr, DoubleMatrix lBS2Lr,
320  DoubleMatrix lTauS2Lr, double costhDRbar) const;
322  //for use in doCalcTadpole1oneLoop
323  double doCalcTad1fermions(double q, double v1) const;
325  //for use in doCalcTadpole1oneLoop
326  double doCalcTad2fermions(double q) const;
328  // Follwing BPMZ Goldstone bosons are not included in this.
329  double doCalcTad1Higgs(double q, double costhDRbar2, double g,
330  double tanb) const;
331  //one loop H2 tadpole contributions from Higgs bosons in the loops
332  // Follwing BPMZ Goldstone bosons are not included in this.
333  double doCalcTad2Higgs(double q, double costhDRbar2, double g,
334  double tanb) const;
336  double doCalcTad1Neutralinos(double q, double costhDRbar, double g,
337  double cosb) const;
339  double doCalcTad2Neutralinos(double q, double costhDRbar, double g,
340  double sinb) const;
342  double doCalcTad1Charginos(double q, double g, double cosb) const;
343 
344  double doCalcTad2Charginos(double q, double g, double sinb) const;
345 
347  double doCalcTad1GaugeBosons(double q, double costhDRbar2, double g,
348  double tanb) const;
349 
350  double doCalcTad2GaugeBosons(double q, double costhDRbar2,
351  double g, double tanb) const;
353  void doTadpoles(double mt, double sinthDRbar);
355  virtual double doCalcTadpole1oneLoop(double mt, double sinthDRbar);
357  virtual double doCalcTadpole2oneLoop(double mt, double sinthDRbar);
360  virtual void calcTadpole1Ms1loop(double mt, double sinthDRbar);
363  virtual void calcTadpole2Ms1loop(double mt, double sinthDRbar);
367  DoubleMatrix addStopQCD(double p, double mt, DoubleMatrix & strong);
371  DoubleMatrix addStopStop(double p, double mt, DoubleMatrix & stop);
375  DoubleMatrix addStopSbottom(double p, double mt, DoubleMatrix & sbottom);
379  DoubleMatrix addStopHiggs(double p, double mt, DoubleMatrix & higgs);
383  DoubleMatrix addStopEweak(double p, DoubleMatrix & electroweak);
387  DoubleMatrix addStopNeutralino(double p, double mt,
388  DoubleMatrix & neutralino);
392  DoubleMatrix addStopChargino(double p, DoubleMatrix & chargino);
396  virtual void addStopCorrection(double p, DoubleMatrix & mass, double mt);
400  DoubleMatrix addSdownQCD(double p1, double p2, int family,
401  DoubleMatrix & strong);
405  DoubleMatrix addSdownHiggs(double p1, double p2, int family,
406  DoubleMatrix & higgs);
410  DoubleMatrix addSdownEweak(double p1, double p2, int family,
411  DoubleMatrix & electroweak);
415  DoubleMatrix addSdownNeutralino(double p1, double p2, int family,
416  DoubleMatrix & neutralino);
420  DoubleMatrix addSdownChargino(double p1, double p2, int family,
421  DoubleMatrix & chargino);
425  virtual void addSdownCorrection(DoubleMatrix & mass, int family);
429  DoubleMatrix addSbotQCD(double p, double mt, DoubleMatrix & strong);
433  void addSbotSfermion(double p, double mt, DoubleMatrix & stop,
434  DoubleMatrix & sbottom);
438  DoubleMatrix addSbotHiggs(double p, double mt, DoubleMatrix & higgs);
442  DoubleMatrix addSbotEweak(double p, DoubleMatrix & electroweak);
446  DoubleMatrix addSbotNeutralino(double p, double mt,
447  DoubleMatrix & neutralino);
451  DoubleMatrix addSbotChargino(double p, double mt, DoubleMatrix & chargino);
455  virtual void addSbotCorrection(double p, DoubleMatrix & mass, double mb);
459  DoubleMatrix addSlepHiggs(double p1, double p2, int family,
460  DoubleMatrix & higgs);
464  DoubleMatrix addSlepEweak(double p1, double p2, int family,
465  DoubleMatrix & electroweak);
469  void addSlepGaugino(double p1, double p2, int family,
470  DoubleMatrix & chargino, DoubleMatrix & neutralino);
474  virtual void addSlepCorrection(DoubleMatrix & mass, int family);
477  void addSlepCorrection(double p, DoubleMatrix & mass, int family);
481  void addStauSfermion(double p, double mtau, DoubleMatrix & stop,
482  DoubleMatrix & sbottom);
486  void addStauGaugino(double p, double mtau, DoubleMatrix & chargino,
487  DoubleMatrix & neutralino);
491  DoubleMatrix addStauEweak(double p, double mtau, DoubleMatrix & electroweak);
495  DoubleMatrix addStauHiggs(double p, double mtau, DoubleMatrix & higgs);
499  virtual void addStauCorrection(double p, DoubleMatrix & mass, double mtau);
503  DoubleMatrix addSupQCD(double p1, double p2, int family,
504  DoubleMatrix & strong);
508  DoubleMatrix addSupHiggs(double p1, double p2, int family,
509  DoubleMatrix & higgs);
513  DoubleMatrix addSupEweak(double p1, double p2, int family,
514  DoubleMatrix & electroweak);
518  DoubleMatrix addSupNeutralino(double p1, double p2, int family,
519  DoubleMatrix & neutralino);
523  DoubleMatrix addSupChargino(double p1, double p2, int family,
524  DoubleMatrix & chargino);
528  virtual void addSupCorrection(DoubleMatrix & mass, int family);
531  void addSnuTauSfermion(double p, double & stop, double & sbottom);
534  double addSnuTauHiggs(double p, double & higgs);
537  double addSnuTauEweak(double p, double & electroweak);
541  void addSnuTauGaugino(double p, double & chargino, double & neutralino);
545  virtual void addSnuTauCorrection(double & mass);
549  double addSnuHiggs(double p, int family, double & higgs);
553  double addSnuEweak(double p, int family, double & electroweak);
557  void addSnuGaugino(double p, int family, double & chargino,
558  double & neutralino);
562  virtual void addSnuCorrection(double & mass, int family);
572  virtual void doUpSquarks(double mt, double pizztMS, double sinthDRbarMS, int
573  accuracy);
578  virtual void doDownSquarks(double mb, double pizztMS, double sinthDRbarMS,
579  int accuracy, double mt);
584  virtual void doChargedSleptons(double mT, double pizztMS,
585  double sinthDRbarMS, int accuracy);
588  virtual void doSnu(double pizztMS, int accuracy = 0);
594  virtual void treeUpSquark(DoubleMatrix & mass, double mtrun, double pizztMS,
595  double sinthDRbarMS, int family);
601  virtual void treeDownSquark(DoubleMatrix & mass, double mbrun, double pizztMS,
602  double sinthDRbarMS, int family);
608  virtual void treeChargedSlepton(DoubleMatrix & mass, double mTrun,
609  double pizztMS, double sinthDRbarMS,
610  int family);
613  void treeSnu(double & mSq, double pizztMS, int family);
614 
616  virtual void calcDrBarPars();
619  void setNeutCurrCouplings(double sinthDRbar, double & sw2, double & guL,
620  double & gdL, double & geL, double & guR,
621  double & gdR, double & geR );
623  void calcDRTrilinears(drBarPars & eg, double vev, double beta);
624 
625  void calcDrBarHiggs(double beta, double mz2, double mw2,
626  double sinthDRbar, drBarPars & eg);
629  void treeCharginos(DoubleMatrix & mCh, double beta, double mw);
632  void treeNeutralinos(DoubleMatrix & mN, double beta, double mz,
633  double mw, double sinth);
634  // calculates the chargino and neutralino DRbar parameters.
635  //It will fill in the chargino and neutralino
637  void calcDrBarGauginos(double beta, double mw, double mz, double sinth,
638  drBarPars & eg);
642  virtual void sparticleThresholdCorrections(double tb);
643  // const DoubleVector & pars);
647  double qedSusythresh(double alphaEm, double Q) const;
651  double qcdSusythresh(double alphasMSbar, double Q);
655  double calcMs() const;
658  virtual void physical(int accuracy);
660  //to the pole mt for use in calcRunningMt
661  virtual double calcRunMtQCD() const;
663  //to the pole mt for use in calcRunningMt
664  virtual double calcRunMtStopGluino() const;
666  //to the pole mt for use in calcRunningMt
667  virtual double calcRunMtHiggs() const;
668 
670  //to the pole mt for use in calcRunningMt
671  virtual double calcRunMtNeutralinos() const;
673  //to the pole mt for use in calcRunningMt
674  virtual double calcRunMtCharginos() const;
675 
678  virtual double calcRunningMt();
680  virtual double calcRunMtauDrBarConv() const;
681  // Obtains (1 / mTAU) times 1-loop squark-chargino corrections
682  //to mtau for use in calcRunningMtau
683  virtual double calcRunMtauCharginos(double mTauSMMZ) const;
684  // Obtains (1 / mTAU) times 1-loop squark-neutralino corrections
685  //to mtau for use in calcRunningMtau
686  virtual double calcRunMtauNeutralinos(double mTauSMMZ) const;
688  // to mtau for use in calcRunningMtau
689  virtual double calcRunMtauHiggs() const;
692  virtual double calcRunningMtau() ;
693 
695  virtual double calcRunMbDrBarConv() const;
697  //to mb for use in calcRunningMb
698  virtual double calcRunMbSquarkGluino() const;
700  //to mb for use in calcRunningMb
701  virtual double calcRunMbChargino() const;
703  // to mb for use in calcRunningMb
704  virtual double calcRunMbHiggs() const;
706  // to mb for use in calcRunningMb
707  virtual double calcRunMbNeutralinos() const;
710  // rruiz: remove const
711  virtual double calcRunningMb();
712  /*
715  double calcHt(double vev);
718  double calcHb(double vev) const;
721  double calcHtau(double vev) const;
722  */
724  double calcSinthdrbar() const;
726  double getVev();
729  double getVev(double pizzt);
733  virtual void charginos(int accuracy, double piwwt);
736  virtual void addChaLoopSfermion(double p, DoubleMatrix & sigmaL, DoubleMatrix & sigmaR, DoubleMatrix & sigmaS) const;
739  virtual void addChaLoopGauge(double p, DoubleMatrix & sigmaL,
740  DoubleMatrix & sigmaR, DoubleMatrix & sigmaS,
741  DoubleMatrix b1pCha, DoubleMatrix b0pCha,
742  DoubleMatrix b1pNeut,
743  DoubleMatrix b0pNeut) const;
746  virtual void addChaLoopHiggs(double p, DoubleMatrix & sigmaL,
747  DoubleMatrix & sigmaR, DoubleMatrix & sigmaS,
748  DoubleMatrix b1pCha, DoubleMatrix b0pCha,
749  DoubleMatrix b1pNeut,
750  DoubleMatrix b0pNeut) const;
752  virtual void addCharginoLoop(double p, DoubleMatrix &);
757  virtual void neutralinos(int accuracy, double piwwt, double pizzt);
761  void addNeutLoopSfermion(double p, DoubleMatrix & sigmaL,
762  DoubleMatrix & sigmaR, DoubleMatrix & sigmaS);
764  void getNeutPassarinoVeltman(double p, double q, DoubleMatrix & b0fn,
765  DoubleMatrix & b1fn);
769  void addNeutLoopGauge(double p, DoubleMatrix & sigmaL,
770  DoubleMatrix & sigmaR, DoubleMatrix & sigmaS);
774  void addNeutLoopHiggs(double p, DoubleMatrix & sigmaL,
775  DoubleMatrix & sigmaR, DoubleMatrix & sigmaS);
777  virtual void addNeutralinoLoop(double p, DoubleMatrix &);
779  virtual void gluino(int accuracy);
782  void calcHiggsAtScale(int accuracy, double & mt, double & sinthDRbar,
783  double & piwwtMS, double & pizztMS, double q = 0.);
790  bool higgs(int accuracy, double piwwt, double pizzt);
796  // virtual void newhiggs(int accuracy, double piwwt, double pizzt);
800  virtual int rewsbMu(int sgnMu, double & mu) const;
803  virtual int rewsbM3sq(double, double &) const;
806  int rewsbM3sqTree(double, double &) const;
808  void alternativeEwsb(double mt);
816  virtual void rewsb(int sgnMu, double mt, const DoubleVector & pars,
817  double muOld = -6.66e66, double eps = 0.);
820  virtual void rewsbTreeLevel(int sgnMu);
823  double treeLevelMuSq();
833  void iterateMu(double & munew, int sgnMu, double mt,
834  int maxTries, double pizztMS, double sinthDRbar, double tol,
835  int & err);
838  double predTanb(double muSusy = -6.66e66) const;
841  double predMzsq(double & tanb, double muOld = -6.66e66, double eps = 0.);
843  double deltaEW() const;
844 
846  double ewsbCondition1TreeLevel() const;
847  double ewsbCondition2TreeLevel() const;
851  void ewsbConditions(DoubleVector & values) const;
854  void predVevs(DoubleVector & vevs, int & err);
855 
864  const DoubleVector &),
865  const DoubleVector & bcPars, double MX,
866  bool doTop = false);
870  double it1par(int numPar, const DoubleVector & bcPars);
874  double ufb3sl(double);
877  double realMinMs() const;
879  double calcBayesianNaturalness() const;
880 
882  //for p=external momentum, Q=renormalisation scale
883  virtual double piZZTHiggs(double p, double Q, double thetaWDRbar) const;
885  //for p=external momentum, Q=renormalisation scale
886  virtual double piZZTsfermions(double p, double Q) const;
889  virtual double piZZTfermions(double p, double Q, bool usePoleMt) const;
892  virtual double piZZTNeutralinos(double p, double Q, double thetaWDRbar) const;
895  virtual double piZZTCharginos(double p, double q, double thetaWDRbar) const;
898  virtual double piZZT(double p, double Q, bool usePoleMt = false) const;
901  virtual double piWWTHiggs(double p, double q, double thetaWDRbar) const;
904  virtual double piWWTfermions(double p, double Q, bool usePoleMt) const;
907  virtual double piWWTsfermions(double p, double Q) const;
910  virtual double piWWTgauginos(double p, double Q, double thetaWDRbar) const;
913  virtual double piWWT(double p, double Q, bool usePoleMt = false) const;
916  virtual void getNeutralinoCharginoHpmCoup(ComplexMatrix & apph1,
917  ComplexMatrix & apph2,
918  ComplexMatrix & bpph1,
919  ComplexMatrix & bpph2) const;
922  double piHpHmFermions(double p, double q) const;
926  double piHpHmSfermions(double p, double q, double mu) const;
929  double piHpHmGauge(double p, double q) const;
932  double piHpHmHiggs(double p, double q) const;
935  double piHpHmGauginos(double p, double q) const;
938  virtual double piHpHm(double p, double Q) const;
941  double piZGT(double p, double Q) const;
944  virtual double piAA(double p, double Q) const;
946  //self-energy: for p=external momentum, q=renormalisation scale
947  Complex pis1s1Sfermions(double p, double q, DoubleMatrix ls1tt,
948  DoubleMatrix ls1bb, DoubleMatrix ls1tautau) const;
951  Complex pis1s2Sfermions(double p, double q, DoubleMatrix ls1tt,
952  DoubleMatrix ls1bb, DoubleMatrix ls1tautau,
953  DoubleMatrix ls2tt, DoubleMatrix ls2bb,
954  DoubleMatrix ls2tautau) const;
957  Complex pis2s2Sfermions(double p, double q, DoubleMatrix ls2tt,
958  DoubleMatrix ls2bb, DoubleMatrix ls2tautau) const;
961  Complex pis1s1Fermions(double p, double q) const;
964  Complex pis2s2Fermions(double p, double q) const;
967  Complex pis1s1Higgs(double p, double q) const;
970  Complex pis1s2Higgs(double p, double q) const;
973  Complex pis2s2Higgs(double p, double q) const;
976  Complex pis1s1Neutralinos(double p, double q) const;
979  Complex pis1s2Neutralinos(double p, double q) const;
982  Complex pis2s2Neutralinos(double p, double q) const;
985  Complex pis1s1Charginos(double p, double q) const;
988  Complex pis1s2Charginos(double p, double q) const;
991  Complex pis2s2Charginos(double p, double q) const;
994  Complex pis1s1(double p, double q) const;
997  Complex pis1s2(double p, double q) const;
1000  Complex pis2s2(double p, double q) const;
1002  double sinSqThetaEff();
1004  virtual double h1s2Mix();
1013  virtual void rhohat(double & outrho, double & outsin, double alphaMZDRbar,
1014  double pizztMZ, double piwwt0, double piwwtMW,
1015  double tol, int maxTries);
1019  virtual double deltaVb(double outrho, double outsin, double alphaDRbar,
1020  double pizztMZ) const;
1024  virtual double dRho(double outrho, double outsin, double alphaDRbar,
1025  double pizztMZ, double piwwtMW);
1029  virtual double dR(double outrho, double outsin, double alphaDRbar,
1030  double pizztMZ, double piwwt0);
1032  double maxMass() const;
1037  int lsp(double & mass, int & posi, int & posj) const;
1042  int nlsp(double & mass, int & posi, int & posj) const;
1043 
1045  string printShort() const;
1047  string printLong();
1049  string printLongDrbar();
1050 
1052  virtual void printObj() {
1053  cout << this->displayMssmSusy() << this->displayMssmSoftPars();
1054  };
1055 
1057  double thet(double a, double b, double c = 0.0);
1058 
1079  void fixedPointIteration(TMSSMBoundaryCondition boundaryCondition,
1080  double mxGuess,
1081  const DoubleVector & pars, int sgnMu, double tanb,
1082  const QedQcd & oneset, bool gaugeUnification,
1083  bool ewsbBCscale = false);
1087  double lowOrg(void (*boundaryCondition)
1088  (MssmSoftsusy &, const DoubleVector &),
1089  double mxGuess,
1090  const DoubleVector & pars, int sgnMu, double tanb,
1091  const QedQcd & oneset, bool gaugeUnification,
1092  bool ewsbBCscale = false) {
1093  fixedPointIteration(boundaryCondition, mxGuess, pars, sgnMu, tanb, oneset,
1094  gaugeUnification, ewsbBCscale);
1095  return displayMxBC();
1096  };
1097 
1107 
1108  void itLowsoft(int maxTries, int sgnMu, double tol,
1109  double tanb, void (*boundaryCondition)(MssmSoftsusy &,
1110  const DoubleVector &),
1111  const DoubleVector & pars, bool gaugeUnification,
1112  bool ewsbBCscale);
1113 
1117  virtual void rpvSet(const DoubleVector & parameters);
1118 
1122  virtual void methodBoundaryCondition(const DoubleVector & pars);
1123 
1140  void isajetNumbers764
1141  (double & mtopPole, double & mGPole, double & smu, double & mA,
1142  double & tanb, double & mq1l, double & mdr, double & mur, double & meL,
1143  double & meR, double & mql3, double & mdr3, double & mur3, double & mtauL,
1144  double & mtauR, double & at, double & ab, double & atau, double & mq2l,
1145  double & msr, double & mcr, double & mmuL, double & mmuR, double & m1,
1146  double & m2) const;
1148  void isajetInterface764(const char fname[80]) const;
1151  void ssrunInterface764Inside(const char fname [80], fstream & ) const;
1155  void ssrunInterface764(const char fname [80], const char softfname [80]) const;
1159  void isawigInterface764(const char fnamein [80], const char fnameout [80],
1160  const char fnamesoft[80]) const;
1173  virtual void lesHouchesAccordOutput(ostream & out, const char model[],
1174  const DoubleVector & pars,
1175  int sgnMu, double tanb, double qMax,
1176  int numPoints,
1177  bool ewsbBCscale);
1178  /* void slha1(ostream & out, const char model[], const DoubleVector & pars,
1179  int sgnMu, double tanb, double qMax, int numPoints,
1180  bool ewsbBCscale);*/
1184  virtual void setEwsbConditions(const DoubleVector & inputs);
1187  void headerSLHA(ostream & out);
1189  virtual void spinfoSLHA(ostream & out);
1191  void modselSLHA(ostream & out, const char model[]);
1193  void sminputsSLHA(ostream & out);
1195  void minparSLHA(ostream & out, const char model [],
1196  const DoubleVector & pars, double tanb, int sgnMu,
1197  bool ewsbBCscale);
1199  virtual void extparSLHA(ostream & out, const DoubleVector & pars,
1200  bool ewsbBCscale);
1202  void massSLHA(ostream & out);
1204  virtual void higgsMSLHA(ostream & out);
1206  virtual void neutralinoCharginoMSLHA(ostream & out);
1208  virtual void sfermionsSLHA(ostream & out);
1210  void sfermionmixSLHA(ostream & out);
1212  virtual void neutralinoMixingSLHA(ostream & out);
1214  void inomixingSLHA(ostream & out);
1216  virtual void softsusySLHA(ostream & out);
1218  void alphaSLHA(ostream & out);
1220  virtual void hmixSLHA(ostream & out);
1222  void gaugeSLHA(ostream & out);
1224  virtual void yukawasSLHA(ostream & out);
1226  virtual void msoftSLHA(ostream & out);
1228  virtual void drbarSLHA(ostream & out, int numPoints, double qMax, int n);
1229 
1233  void doUfb3(double mgut);
1234 
1238  void assignHiggs(DoubleVector & higgsm, DoubleVector & higgsc,
1239  DoubleVector & dnu, DoubleVector & dnd,
1240  DoubleVector & cn, double beta) const;
1241 
1244  void assignHiggs(DoubleVector & higgsm, DoubleVector & higgsc) const;
1245 
1250  void assignHiggsSfermions(DoubleVector & higgsm, DoubleVector & higgsc,
1251  DoubleVector & dnu, DoubleVector & dnd,
1252  DoubleVector & cn, double beta) const;
1253 
1256  double smPredictionMW() const;
1257 
1260  double twoLoopGm2(double amu1Loop) const;
1261 
1264  double twoLpMt() const;
1265  double twoLpMb() const;
1266 
1271  virtual void doQuarkMixing(DoubleMatrix & mDon, DoubleMatrix & mUpq);
1272 
1277  virtual MssmSusyRGE guessAtSusyMt(double tanb, const QedQcd & oneset);
1278  };
1279 
1280  std::istream& operator>>(std::istream& left, MssmSoftsusy& s);
1281 
1283  void printShortInitialise();
1285  double rho2(double r);
1286 
1291  double minimufb3(double);
1292 
1294  double gEff(double x);
1296  double fEff(double x);
1297 
1299  double lep2Likelihood(double mh);
1302  DoubleVector mhIntegrand(double mh, const DoubleVector & y);
1305  double lnLHiggs(double mh);
1306 
1307 }
1308 
1309 #endif
1310 
1311 
1312 
matrix of complex double values dimensions (rows x cols)
Definition: linalg.h:608
drop-in replacement for the original home-grown Complex class
Definition: mycomplex.h:17
Matrix from 1..rows, 1..cols of double values.
Definition: linalg.h:214
DoubleVector is of variable length, and contains double precision.
Definition: linalg.h:35
int displayEnd() const
returns end of dimension
Definition: linalg.h:94
A different REWSB condition - given by mu and MA instead of mh1,2.
Definition: softsusy.h:54
Definition: rge.h:23
Soft SUSY breaking parameters and beta functions.
Definition: softpars.h:31
const MssmSoftPars & displayMssmSoftPars() const
Returns whole object as a const.
Definition: softpars.h:66
DoubleVector beta(const MssmSusy &) const
Returns double vector containing numerical beta functions of parameters.
Definition: softpars.cpp:2858
Contains all supersymmetric MSSM parameters, incorporating R_p MSSM.
Definition: softsusy.h:74
Complex pis1s2Sfermions(double p, double q, DoubleMatrix ls1tt, DoubleMatrix ls1bb, DoubleMatrix ls1tautau, DoubleMatrix ls2tt, DoubleMatrix ls2bb, DoubleMatrix ls2tautau) const
Definition: softsusy.cpp:11056
void isajetInterface764(const char fname[80]) const
prints a file into fname which acts as an input to isajet
Definition: softsusy.cpp:9939
double twoLoopGm2(double amu1Loop) const
Definition: softsusy.cpp:10765
Complex pis2s2(double p, double q) const
checked 28.10.02
Definition: softsusy.cpp:11783
DoubleMatrix addSbotChargino(double p, double mt, DoubleMatrix &chargino)
Definition: softsusy.cpp:5962
void calcDrBarHiggs(double beta, double mz2, double mw2, double sinthDRbar, drBarPars &eg)
Definition: softsusy.cpp:7683
double doCalcTad1Sfermions(DoubleMatrix lTS1Lr, DoubleMatrix lBS1Lr, DoubleMatrix lTauS1Lr, double costhDRbar) const
routine to calculate sfermiom contributions to H1 tadpole / v1
Definition: softsusy.cpp:392
double smPredictionMW() const
checked 22/04/06
Definition: softsusy.cpp:10737
void sminputsSLHA(ostream &out)
SMINPUTS block of SLHA.
Definition: softsusy.cpp:10371
virtual void neutralinoMixingSLHA(ostream &out)
nmix block output of SLHA
Definition: softsusy.cpp:10134
double it1par(int numPar, const DoubleVector &bcPars)
Definition: softsusy.cpp:1149
virtual double calcRunMtCharginos() const
Obtains (16 Pi^2) / mt times 1-loop Charginos corrections.
Definition: softsusy.cpp:2807
void flagNotGlobalMin(bool a)
LCT: Flags problem if not in global minimum of Higgs potential.
Definition: softsusy.h:232
virtual void addChaLoopSfermion(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS) const
Definition: softsusy.cpp:1711
virtual double dR(double outrho, double outsin, double alphaDRbar, double pizztMZ, double piwwt0)
Definition: softsusy.cpp:9518
virtual void lesHouchesAccordOutput(ostream &out, const char model[], const DoubleVector &pars, int sgnMu, double tanb, double qMax, int numPoints, bool ewsbBCscale)
SUSY Les Houches accord for interfacing to Monte-Carlos, decay programs etc.
Definition: softsusy.cpp:10532
void setProblem(const sProblem a)
Sets all problem flags.
Definition: softsusy.h:244
virtual double piZZTHiggs(double p, double Q, double thetaWDRbar) const
Calculates Higgs contribution to the transverse part of Z self-energy:
Definition: softsusy.cpp:8115
void H2SfSfCouplings(DoubleMatrix &lTS2Lr, DoubleMatrix &lBS2Lr, DoubleMatrix &lTauS2Lr, double gmzOcthW, double mu, double sinb) const
fixes trilnear H2-sfermion-sfermion couplings
Definition: softsusy.cpp:368
void isawigInterface764(const char fnamein[80], const char fnameout[80], const char fnamesoft[80]) const
Definition: softsusy.cpp:9922
DoubleVector fineTune(void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), const DoubleVector &bcPars, double MX, bool doTop=false)
Pass it an object and it'll return the fine tuning parameters.
Definition: softsusy.cpp:1229
virtual double piWWTgauginos(double p, double Q, double thetaWDRbar) const
Definition: softsusy.cpp:8477
double qcdSusythresh(double alphasMSbar, double Q)
Definition: softsusy.cpp:7031
double doCalcTad2fermions(double q) const
fixes trilnear H2-fermion-fermion couplings
Definition: softsusy.cpp:537
void addNeutLoopSfermion(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)
Definition: softsusy.cpp:2129
void gaugeSLHA(ostream &out)
sfermionic part of mixing blocks
Definition: softsusy.cpp:10230
void setSoftsusy(const MssmSoftsusy &s)
Sets whole object equal to another
Definition: softsusy.h:247
virtual void softsusySLHA(ostream &out)
SOFTSUSY comments in SLHA.
Definition: softsusy.cpp:10017
virtual double piZZTfermions(double p, double Q, bool usePoleMt) const
fermion contribution to the Transverse part of Z self-energy:
Definition: softsusy.cpp:8253
void setQewsb(double q)
Sets user-set scale qewsb.
Definition: softsusy.h:278
Complex pis1s1Charginos(double p, double q) const
Definition: softsusy.cpp:11637
virtual double deltaVb(double outrho, double outsin, double alphaDRbar, double pizztMZ) const
outrho, outsin represent the DRbar values
Definition: softsusy.cpp:9301
void flagMgutOutOfBounds(bool a)
Flags weird mgut-type problems.
Definition: softsusy.h:201
void calcDrBarGauginos(double beta, double mw, double mz, double sinth, drBarPars &eg)
masses in eg with the appropriate values on exit.
Definition: softsusy.cpp:7665
virtual double calcRunMbChargino() const
Obtains (1 / mb) times 1-loop squark-chargino corrections.
Definition: softsusy.cpp:2979
void flagNoRhoConvergence(bool a)
Flags fact that rho parameter sub iteration didn't converge.
Definition: softsusy.h:223
virtual void treeDownSquark(DoubleMatrix &mass, double mbrun, double pizztMS, double sinthDRbarMS, int family)
Definition: softsusy.cpp:6544
void massSLHA(ostream &out)
This does the job of the above method, but outputs the Mass block.
Definition: softsusy.cpp:10077
virtual double piAA(double p, double Q) const
Definition: softsusy.cpp:8819
virtual double calcRunMbDrBarConv() const
calculates factor to convert to DrBar for mb
Definition: softsusy.cpp:2943
virtual void addSdownCorrection(DoubleMatrix &mass, int family)
Definition: softsusy.cpp:5533
void ssrunInterface764(const char fname[80], const char softfname[80]) const
Definition: softsusy.cpp:9874
DoubleMatrix calcHiggs3L(bool is_bottom)
Higgs 3L corrections.
Definition: softsusy.cpp:11809
DoubleMatrix addStopNeutralino(double p, double mt, DoubleMatrix &neutralino)
Definition: softsusy.cpp:4310
virtual double piZZT(double p, double Q, bool usePoleMt=false) const
Transverse part of Z self-energy: has been checked.
Definition: softsusy.cpp:8353
DoubleMatrix addSbotQCD(double p, double mt, DoubleMatrix &strong)
Definition: softsusy.cpp:5569
void treeSnu(double &mSq, double pizztMS, int family)
Definition: softsusy.cpp:6772
DoubleMatrix addStauEweak(double p, double mtau, DoubleMatrix &electroweak)
Definition: softsusy.cpp:4969
void addSquarkCorrection(DoubleMatrix &mass)
Now these are calculated at the squark scale.
Definition: softsusy.cpp:3388
void setMxBC(double mx)
Sets the scale at which high-scale boundary conditions are applied.
Definition: softsusy.h:260
int included_thresholds
Flag allowing to choose which two-loop thresholds have to be included.
Definition: softsusy.h:121
virtual void addSnuCorrection(double &mass, int family)
Found+fixed bug 7/09/06. Thanks to J Kersten.
Definition: softsusy.cpp:3926
void addSlepCorrection(double p, DoubleMatrix &mass, int family)
virtual void addSbotCorrection(double p, DoubleMatrix &mass, double mb)
Checked 16.08.2005.
Definition: softsusy.cpp:6020
void predVevs(DoubleVector &vevs, int &err)
DH: solves the EWSB conditions for the Higgs VEVs.
Definition: softsusy.cpp:973
void inomixingSLHA(ostream &out)
This does the job of the above method, but outputs the UMIX/VMIX blocks.
Definition: softsusy.cpp:10149
double ewsbCondition1TreeLevel() const
DH: returns the values of the tree-level EWSB conditions.
Definition: softsusy.cpp:854
virtual void msoftSLHA(ostream &out)
MSOFT block of SLHA.
Definition: softsusy.cpp:10268
virtual double calcRunningMtau()
Full BPMZ expression.
Definition: softsusy.cpp:3291
const MssmSoftsusy & displaySoftsusy() const
Returns object as a const.
Definition: softsusy.h:187
double addSnuEweak(double p, int family, double &electroweak)
Definition: softsusy.cpp:3783
void flagIrqfp(bool a)
Flags Infra-red quasi fixed point breach problem.
Definition: softsusy.h:203
virtual double piZZTNeutralinos(double p, double Q, double thetaWDRbar) const
Definition: softsusy.cpp:8293
MssmSoftsusy()
Default constructor fills object with zeroes.
Definition: softsusy.cpp:67
DoubleMatrix addStopHiggs(double p, double mt, DoubleMatrix &higgs)
Definition: softsusy.cpp:4017
virtual double piWWTsfermions(double p, double Q) const
Definition: softsusy.cpp:8436
Complex pis1s2Neutralinos(double p, double q) const
Definition: softsusy.cpp:11566
double doCalcTad1GaugeBosons(double q, double costhDRbar2, double g, double tanb) const
one loop H1 tadpole contributions from Charginos in the loops
Definition: softsusy.cpp:642
virtual void printObj()
Prints whole object to standard output.
Definition: softsusy.h:1052
virtual double piWWT(double p, double Q, bool usePoleMt=false) const
Definition: softsusy.cpp:8515
void minparSLHA(ostream &out, const char model[], const DoubleVector &pars, double tanb, int sgnMu, bool ewsbBCscale)
MINPAR block of SLHA.
Definition: softsusy.cpp:10456
const drBarPars & displayDrBarPars() const
Definition: softsusy.cpp:164
virtual double piHpHm(double p, double Q) const
LCT: Charged Higgs self-energy.
Definition: softsusy.cpp:8797
virtual double calcRunningMt()
Definition: softsusy.cpp:2869
void calcHiggsAtScale(int accuracy, double &mt, double &sinthDRbar, double &piwwtMS, double &pizztMS, double q=0.)
Definition: softsusy.cpp:6782
virtual void addNeutralinoLoop(double p, DoubleMatrix &)
Adds the loop corrections on to an input tree-level neutralino mass.
Definition: softsusy.cpp:2518
void addSbotSfermion(double p, double mt, DoubleMatrix &stop, DoubleMatrix &sbottom)
Definition: softsusy.cpp:5594
double deltaEW() const
Returns fine-tuning as defined in Baer et al arXiv:1212.2655 eq 1.2.
Definition: softsusy.cpp:1210
virtual void sparticleThresholdCorrections(double tb)
Definition: softsusy.cpp:7556
virtual void getNeutralinoCharginoHpmCoup(ComplexMatrix &apph1, ComplexMatrix &apph2, ComplexMatrix &bpph1, ComplexMatrix &bpph2) const
LCT: Returns trilinear neutralino-chargino-hpm coupling in unrotated basis.
Definition: softsusy.cpp:8535
double getVev()
Calculates Higgs VEV parameter from gauge couplings and MZ.
Definition: softsusy.cpp:7543
virtual void rpvSet(const DoubleVector &parameters)
Definition: softsusy.cpp:9643
double calcMs() const
DRbar pars should be defined already for this.
Definition: softsusy.cpp:7250
double displayTadpole2Ms1loop() const
displays t_2/v_2 tadpole@1 loop
Definition: softsusy.cpp:184
virtual void doSnu(double pizztMS, int accuracy=0)
Definition: softsusy.cpp:6759
void treeCharginos(DoubleMatrix &mCh, double beta, double mw)
Definition: softsusy.cpp:7645
virtual double calcRunMtauHiggs() const
Obtains (1 / mtau) times 1-loop Higgs corrections.
Definition: softsusy.cpp:3191
void setSetTbAtMX(bool a)
Sets the setTbAtMX flag.
Definition: softsusy.h:264
void addNeutLoopHiggs(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)
Definition: softsusy.cpp:2376
double doCalcTad1Higgs(double q, double costhDRbar2, double g, double tanb) const
one loop H1 tadpole contributions from Higgs bosons in the loops
Definition: softsusy.cpp:546
Complex pis1s1(double p, double q) const
Definition: softsusy.cpp:11719
virtual double dRho(double outrho, double outsin, double alphaDRbar, double pizztMZ, double piwwtMW)
Definition: softsusy.cpp:9479
double displayMz() const
Returns value of pole MZ being used.
Definition: softsusy.h:189
DoubleMatrix addStauHiggs(double p, double mtau, DoubleMatrix &higgs)
Definition: softsusy.cpp:4846
void flagNoConvergence(bool a)
Flags fact that calculation hasn't acheived required accuracy.
Definition: softsusy.h:219
Complex pis2s2Fermions(double p, double q) const
Definition: softsusy.cpp:11252
void flagNoMuConvergence(bool a)
Flags fact that mu sub iteration didn't converge.
Definition: softsusy.h:221
sPhysical displayPhysUncertainty(TMSSMBoundaryCondition boundaryCondition, double mxGuess, const DoubleVector &pars, int sgnMu, double tanb, const QedQcd &oneset, bool gaugeUnification, bool ewsbBCscale=false, int contributions=DeltaAll) const
Displays uncertainty estimate of pole masses.
Definition: softsusy.cpp:12282
virtual double doCalcTadpole2oneLoop(double mt, double sinthDRbar)
Does the calculation of one-loop pieces of .
Definition: softsusy.cpp:729
virtual void rewsbTreeLevel(int sgnMu)
Definition: softsusy.cpp:1424
void check_flags()
check consistency of threshold flags
Definition: softsusy.cpp:7322
virtual void extparSLHA(ostream &out, const DoubleVector &pars, bool ewsbBCscale)
EXTPAR block of SLHA.
Definition: softsusy.cpp:10387
virtual double calcRunMtQCD() const
Obtains (16 Pi^2) / mt times 1-loop and 2-loop QCD corrections.
Definition: softsusy.cpp:2653
void useAlternativeAlphaS()
Use alernative alpha_s(MSSM,DR-bar) calculation.
Definition: softsusy.h:270
void H1SfSfCouplings(DoubleMatrix &lTS1Lr, DoubleMatrix &lBS1Lr, DoubleMatrix &lTauS1Lr, double gmzOcthW, double mu, double cosb, double v1) const
fixes trilnear H1-sfermion-sfermion couplings
Definition: softsusy.cpp:343
void setData(const QedQcd &r)
Sets low energy Standard Model fermion mass and gauge coupling data.
Definition: softsusy.h:249
void sfermionmixSLHA(ostream &out)
hmix block output of SLHA
Definition: softsusy.cpp:10195
double doCalcTad1fermions(double q, double v1) const
fixes trilnear H1-fermion-fermion couplings
Definition: softsusy.cpp:523
double displayMinpot() const
Returns minimum of Higgs potential.
Definition: softsusy.cpp:168
void flagNonperturbative(bool a)
Flags non-perturbative RG evolution.
Definition: softsusy.h:205
virtual void addCharginoLoop(double p, DoubleMatrix &)
Adds the loop corrections on to an input tree-level chargino mass.
Definition: softsusy.cpp:2042
virtual void doChargedSleptons(double mT, double pizztMS, double sinthDRbarMS, int accuracy)
Definition: softsusy.cpp:6671
DoubleMatrix addSupEweak(double p1, double p2, int family, DoubleMatrix &electroweak)
Definition: softsusy.cpp:6204
void itLowsoft(int maxTries, int sgnMu, double tol, double tanb, void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), const DoubleVector &pars, bool gaugeUnification, bool ewsbBCscale)
Definition: softsusy.cpp:7909
double displayMw() const
Definition: softsusy.cpp:170
void useAlternativeEwsb()
Use alernative EWSB conditions: set mu and MA(pole)
Definition: softsusy.h:266
double doCalcTad2Neutralinos(double q, double costhDRbar, double g, double sinb) const
one loop H2 tadpole contributions from Neutralinos in the loops
Definition: softsusy.cpp:602
double realMinMs() const
Definition: softsusy.cpp:7500
void setPhys(const sPhysical &s)
Sets all physical parameters.
Definition: softsusy.h:258
void setTwoLoopMtThresholds(bool sw)
Switch 2-loop threshold corrections to top mass.
Definition: softsusy.cpp:10896
double qedSusythresh(double alphaEm, double Q) const
Definition: softsusy.cpp:7103
Complex pis1s1Higgs(double p, double q) const
Definition: softsusy.cpp:11261
double doCalcTad1Charginos(double q, double g, double cosb) const
one loop H1 tadpole contributions from Charginos in the loops
Definition: softsusy.cpp:617
string printLong()
Prints a list of all sparticle/Higgs masses to standard output.
Definition: softsusy.cpp:1578
virtual void yukawasSLHA(ostream &out)
sfermionic part of mixing blocks
Definition: softsusy.cpp:10241
void flagM3sq(bool a)
Flags problem with Higgs potential minimum.
Definition: softsusy.h:217
DoubleMatrix addStopChargino(double p, DoubleMatrix &chargino)
Definition: softsusy.cpp:4371
virtual void methodBoundaryCondition(const DoubleVector &pars)
Definition: softsusy.cpp:9636
Complex pis2s2Sfermions(double p, double q, DoubleMatrix ls2tt, DoubleMatrix ls2bb, DoubleMatrix ls2tautau) const
Definition: softsusy.cpp:11129
double addSnuTauHiggs(double p, double &higgs)
Definition: softsusy.cpp:3431
virtual void addStauCorrection(double p, DoubleMatrix &mass, double mtau)
Definition: softsusy.cpp:5177
virtual double calcRunMtHiggs() const
Obtains (16 Pi^2) / mt times 1-loop Higgs corrections.
Definition: softsusy.cpp:2690
virtual void treeChargedSlepton(DoubleMatrix &mass, double mTrun, double pizztMS, double sinthDRbarMS, int family)
Definition: softsusy.cpp:6645
virtual void rhohat(double &outrho, double &outsin, double alphaMZDRbar, double pizztMZ, double piwwt0, double piwwtMW, double tol, int maxTries)
Definition: softsusy.cpp:9561
DoubleMatrix addSlepEweak(double p1, double p2, int family, DoubleMatrix &electroweak)
Definition: softsusy.cpp:4605
DoubleMatrix addStopQCD(double p, double mt, DoubleMatrix &strong)
Definition: softsusy.cpp:3955
double piHpHmFermions(double p, double q) const
LCT: Returns fermion contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8548
bool higgs(int accuracy, double piwwt, double pizzt)
Definition: softsusy.cpp:11862
virtual void calcTadpole2Ms1loop(double mt, double sinthDRbar)
Definition: softsusy.cpp:761
DoubleMatrix addSbotHiggs(double p, double mt, DoubleMatrix &higgs)
Definition: softsusy.cpp:5631
void doUfb3(double mgut)
Definition: softsusy.cpp:191
virtual void gluino(int accuracy)
Calculates pole gluino mass to 1-loop SUSY corrections.
Definition: softsusy.cpp:2602
Complex pis2s2Higgs(double p, double q) const
Definition: softsusy.cpp:11450
int lsp(double &mass, int &posi, int &posj) const
Definition: softsusy.cpp:7171
void flagMusqwrongsign(bool a)
Flags point inconsistent with electroweak symmetry breaking.
Definition: softsusy.h:225
double displayMzRun() const
Returns DRbar MZ, must be calculated by calcDrBarPars first.
Definition: softsusy.cpp:7640
sPhysical displayPhysUncertaintyScaleVariation() const
Displays uncertainty estimate of pole masses by scale variation only.
Definition: softsusy.cpp:12238
virtual void doUpSquarks(double mt, double pizztMS, double sinthDRbarMS, int accuracy)
Definition: softsusy.cpp:6470
virtual void rewsb(int sgnMu, double mt, const DoubleVector &pars, double muOld=-6.66e66, double eps=0.)
Definition: softsusy.cpp:1467
virtual void doDownSquarks(double mb, double pizztMS, double sinthDRbarMS, int accuracy, double mt)
Definition: softsusy.cpp:6572
const QedQcd & displayDataSet() const
Gives the low energy Standard Model data set used for the object.
Definition: softsusy.cpp:160
DoubleMatrix addSdownChargino(double p1, double p2, int family, DoubleMatrix &chargino)
Definition: softsusy.cpp:5487
virtual double calcRunMtauCharginos(double mTauSMMZ) const
Definition: softsusy.cpp:3149
void setPredMzSq(double a)
Set MZ^2 predicted after iteration.
Definition: softsusy.h:274
double predTanb(double muSusy=-6.66e66) const
Definition: softsusy.cpp:264
virtual void charginos(int accuracy, double piwwt)
checked
Definition: softsusy.cpp:2090
virtual double calcRunMbNeutralinos() const
Obtains (1 / mb) times 1-loop neutralino-suqrk corrections.
Definition: softsusy.cpp:12168
virtual const DoubleVector display() const
Return contents of object in a vector: for RG evolution.
Definition: softsusy.cpp:803
DoubleMatrix addStopEweak(double p, DoubleMatrix &electroweak)
Definition: softsusy.cpp:4182
virtual double piZZTCharginos(double p, double q, double thetaWDRbar) const
Definition: softsusy.cpp:8323
double ufb3sl(double)
Definition: softsusy.cpp:7003
DoubleMatrix addSupQCD(double p1, double p2, int family, DoubleMatrix &strong)
Definition: softsusy.cpp:6057
virtual void drbarSLHA(ostream &out, int numPoints, double qMax, int n)
outputs DRbar parameters at scale Q (default is at MSUSY)
Definition: softsusy.cpp:10337
Complex pis1s1Sfermions(double p, double q, DoubleMatrix ls1tt, DoubleMatrix ls1bb, DoubleMatrix ls1tautau) const
Calculates (16 Pi^2) times the sfermion contribution to Higgs.
Definition: softsusy.cpp:10945
double lowOrg(void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), double mxGuess, const DoubleVector &pars, int sgnMu, double tanb, const QedQcd &oneset, bool gaugeUnification, bool ewsbBCscale=false)
Definition: softsusy.h:1087
double displaySoftA(trilinears, int, int) const
Return a trilinear element in "SUGRA style".
Definition: softsusy.cpp:837
void treeNeutralinos(DoubleMatrix &mN, double beta, double mz, double mw, double sinth)
Definition: softsusy.cpp:7652
virtual void addSnuTauCorrection(double &mass)
Definition: softsusy.cpp:3671
virtual void addSlepCorrection(DoubleMatrix &mass, int family)
16.09.05 checked.
Definition: softsusy.cpp:4779
double treeLevelMuSq()
Definition: softsusy.cpp:1419
Complex pis1s2Higgs(double p, double q) const
Definition: softsusy.cpp:11346
virtual void spinfoSLHA(ostream &out)
This does the SPINFO block of SLHA.
Definition: softsusy.cpp:9990
virtual void calcTadpole1Ms1loop(double mt, double sinthDRbar)
checked
Definition: softsusy.cpp:713
virtual double piWWTHiggs(double p, double q, double thetaWDRbar) const
Definition: softsusy.cpp:8378
void useAlternativeAlphaEm()
Use alernative alpha_em(MSSM,DR-bar) calculation.
Definition: softsusy.h:272
void addSnuTauSfermion(double p, double &stop, double &sbottom)
Definition: softsusy.cpp:3417
void setMsusy(double)
Definition: softsusy.cpp:189
void isajetNumbers764(double &mtopPole, double &mGPole, double &smu, double &mA, double &tanb, double &mq1l, double &mdr, double &mur, double &meL, double &meR, double &mql3, double &mdr3, double &mur3, double &mtauL, double &mtauR, double &at, double &ab, double &atau, double &mq2l, double &msr, double &mcr, double &mmuL, double &mmuR, double &m1, double &m2) const
Definition: softsusy.cpp:9730
DoubleMatrix addSupNeutralino(double p1, double p2, int family, DoubleMatrix &neutralino)
Definition: softsusy.cpp:6323
double doCalcTad2Higgs(double q, double costhDRbar2, double g, double tanb) const
Definition: softsusy.cpp:566
double displayTadpole1Ms1loop() const
displays t_1/v_1 tadpole@1 loop
Definition: softsusy.cpp:180
virtual void addChaLoopHiggs(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS, DoubleMatrix b1pCha, DoubleMatrix b0pCha, DoubleMatrix b1pNeut, DoubleMatrix b0pNeut) const
Definition: softsusy.cpp:1909
double piHpHmSfermions(double p, double q, double mu) const
LCT: Returns sfermion contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8569
void addStauGaugino(double p, double mtau, DoubleMatrix &chargino, DoubleMatrix &neutralino)
Definition: softsusy.cpp:5085
void alphaSLHA(ostream &out)
sfermionic part of mixing blocks
Definition: softsusy.cpp:10127
double sinSqThetaEff()
Calculates sin^2 theta^l_eff.
Definition: softsusy.cpp:9273
double displayMsusy() const
Returns Higgs minimisation scale.
Definition: softsusy.cpp:169
void addNeutLoopGauge(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)
Definition: softsusy.cpp:2317
virtual double calcRunningMb()
Definition: softsusy.cpp:3086
double mxBC
Scale at which SUSY breaking boundary conditions set.
Definition: softsusy.h:113
virtual void neutralinos(int accuracy, double piwwt, double pizzt)
mixNeut set to diagonal = mixNeut^T mNeutralino mixNeut: checked
Definition: softsusy.cpp:2540
double calcBayesianNaturalness() const
Calculates Bayesian naturalness a la arXiv:1709.07895.
Definition: softsusy.cpp:7512
void flagAllProblems(bool a, tachyonType b)
Sets all problems equal to either true or false (contained in a)
Definition: softsusy.h:234
void addSlepGaugino(double p1, double p2, int family, DoubleMatrix &chargino, DoubleMatrix &neutralino)
Definition: softsusy.cpp:4707
void setTwoLoopAlphasThresholds(bool sw)
Definition: softsusy.cpp:10888
virtual void calcDrBarPars()
Organises tree-level calculation of all sparticle masses and mixings.
Definition: softsusy.cpp:7815
DoubleMatrix addSdownHiggs(double p1, double p2, int family, DoubleMatrix &higgs)
Definition: softsusy.cpp:5228
string printLongDrbar()
Prints a list of all running sparticle/Higgs masses to standard output.
Definition: softsusy.cpp:1643
double piZGT(double p, double Q) const
Definition: softsusy.cpp:9140
virtual double h1s2Mix()
gets h1 mixing element with Hu.
Definition: softsusy.cpp:9475
void flagHiggsufb(bool a)
Flags an inconsistent Higgs minimum.
Definition: softsusy.h:230
void setNeutCurrCouplings(double sinthDRbar, double &sw2, double &guL, double &gdL, double &geL, double &guR, double &gdR, double &geR)
Definition: softsusy.cpp:7754
double piHpHmHiggs(double p, double q) const
LCT: Returns Higgs contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8699
virtual void setEwsbConditions(const DoubleVector &inputs)
Again, another dummy - useful in alternative EWSB conditions sometimes.
Definition: softsusy.cpp:10883
const sProblem & displayProblem() const
Returns any problem flags associated with the object.
Definition: softsusy.h:166
int nlsp(double &mass, int &posi, int &posj) const
Definition: softsusy.cpp:10567
Complex pis2s2Neutralinos(double p, double q) const
Definition: softsusy.cpp:11605
double addSnuHiggs(double p, int family, double &higgs)
Definition: softsusy.cpp:3703
DoubleMatrix addSbotNeutralino(double p, double mt, DoubleMatrix &neutralino)
Definition: softsusy.cpp:5900
Complex pis1s2(double p, double q) const
Definition: softsusy.cpp:11750
void ssrunInterface764Inside(const char fname[80], fstream &) const
Definition: softsusy.cpp:9884
virtual double calcRunMtauDrBarConv() const
calculates factor to convert to DrBar for mtau
Definition: softsusy.cpp:3143
virtual void doQuarkMixing(DoubleMatrix &mDon, DoubleMatrix &mUpq)
input diagonal matrices and it'll give you back mixed ones
Definition: softsusy.cpp:10938
double calcSinthdrbar() const
Calculates DRbar sin theta_w at the current scale from gauge couplings.
Definition: softsusy.cpp:7518
double displayTadpole1Ms() const
displays t_1/v_1 tadpole
Definition: softsusy.cpp:172
void getNeutPassarinoVeltman(double p, double q, DoubleMatrix &b0fn, DoubleMatrix &b1fn)
LCT: Returns matrix of Passarino-Veltman B0 and B1 functions.
Definition: softsusy.cpp:2282
void flagInaccurateHiggsMass(bool a)
Definition: softsusy.h:228
virtual double calcRunMbSquarkGluino() const
Obtains (1 / mb) times 1-loop sbottom-Gluino corrections.
Definition: softsusy.cpp:2956
DoubleMatrix addSupHiggs(double p1, double p2, int family, DoubleMatrix &higgs)
Definition: softsusy.cpp:6075
virtual double piWWTfermions(double p, double Q, bool usePoleMt) const
Definition: softsusy.cpp:8413
virtual int rewsbM3sq(double, double &) const
returns 1 if mu < 1.0e-9
Definition: softsusy.cpp:221
void set(const DoubleVector &)
Sets total set of RGE parameters equal to elements of a vector.
Definition: softsusy.cpp:830
DoubleMatrix addSdownQCD(double p1, double p2, int family, DoubleMatrix &strong)
Definition: softsusy.cpp:5211
DoubleVector beta() const
Returns double vector containing numerical beta functions of parameters.
Definition: softsusy.h:295
void alternativeEwsb(double mt)
Fixes mH1^2 and mH2^2 in order to get muCond and mAcond correct.
Definition: softsusy.cpp:1339
void calcDRTrilinears(drBarPars &eg, double vev, double beta)
sets the Yukawas and Trilinears
Definition: softsusy.cpp:7768
DoubleMatrix addStopStop(double p, double mt, DoubleMatrix &stop)
Definition: softsusy.cpp:3981
void doTadpoles(double mt, double sinthDRbar)
Does the full 2-loop calculation of both tadpoles and sets them.
Definition: softsusy.cpp:283
void addSnuGaugino(double p, int family, double &chargino, double &neutralino)
Definition: softsusy.cpp:3863
void setTwoLoopMbMtauThresholds(bool sw)
Definition: softsusy.cpp:10904
const MssmSoftsusy & displayMssmSoft() const
Displays whole object as a const.
Definition: softsusy.cpp:158
void addSnuTauGaugino(double p, double &chargino, double &neutralino)
Definition: softsusy.cpp:3605
void flagProblemThrown(bool a)
Flags a numerical exception eg number too big/small.
Definition: softsusy.h:242
Complex pis1s1Neutralinos(double p, double q) const
Definition: softsusy.cpp:11534
DoubleMatrix addSdownNeutralino(double p1, double p2, int family, DoubleMatrix &neutralino)
Definition: softsusy.cpp:5432
double addSnuTauEweak(double p, double &electroweak)
Definition: softsusy.cpp:3522
Complex pis1s1Fermions(double p, double q) const
Definition: softsusy.cpp:11237
DoubleMatrix addStopSbottom(double p, double mt, DoubleMatrix &sbottom)
Definition: softsusy.cpp:3999
void iterateMu(double &munew, int sgnMu, double mt, int maxTries, double pizztMS, double sinthDRbar, double tol, int &err)
Definition: softsusy.cpp:1267
double twoLpMt() const
void addStauSfermion(double p, double mtau, DoubleMatrix &stop, DoubleMatrix &sbottom)
Definition: softsusy.cpp:4814
double doCalcTad2GaugeBosons(double q, double costhDRbar2, double g, double tanb) const
Definition: softsusy.cpp:656
const sPhysical & displayPhys() const
Displays physical parameters only.
Definition: softsusy.cpp:162
DoubleMatrix addSdownEweak(double p1, double p2, int family, DoubleMatrix &electroweak)
Definition: softsusy.cpp:5328
double displayFracDiff() const
Displays iteration accuracy attained.
Definition: softsusy.h:173
double maxMass() const
Returns the mass of the heaviest SUSY particle, excluding gravitino.
Definition: softsusy.cpp:7216
double thet(double a, double b, double c=0.0)
log(max(a^2, b^2, c^2) / Q^2)
Definition: softsusy.cpp:2123
int rewsbM3sqTree(double, double &) const
returns 1 if mu < 1.0e-9
Definition: softsusy.cpp:242
void setAllTwoLoopThresholds(bool sw)
Switch all two-loop thresholds to gauge-Yukawa couplings.
Definition: softsusy.cpp:10918
virtual double doCalcTadpole1oneLoop(double mt, double sinthDRbar)
Does the calculation of one-loop pieces of .
Definition: softsusy.cpp:673
virtual void physical(int accuracy)
Definition: softsusy.cpp:6850
virtual double calcRunMbHiggs() const
Obtains (1 / mb) times 1-loop Higgs corrections.
Definition: softsusy.cpp:3041
void fixedPointIteration(TMSSMBoundaryCondition boundaryCondition, double mxGuess, const DoubleVector &pars, int sgnMu, double tanb, const QedQcd &oneset, bool gaugeUnification, bool ewsbBCscale=false)
Definition: softsusy.cpp:7364
double predMzsq(double &tanb, double muOld=-6.66e66, double eps=0.)
Definition: softsusy.cpp:772
string printShort() const
Prints a list of important sparticle/Higgs masses to standard output.
Definition: softsusy.cpp:1540
virtual int rewsbMu(int sgnMu, double &mu) const
Definition: softsusy.cpp:196
virtual void sfermionsSLHA(ostream &out)
sfermions part of mass block of SLHA
Definition: softsusy.cpp:10095
DoubleMatrix addSbotEweak(double p, DoubleMatrix &electroweak)
Definition: softsusy.cpp:5786
virtual void neutralinoCharginoMSLHA(ostream &out)
neutralino and charigno part of mass block of SLHA
Definition: softsusy.cpp:10060
double displayMwRun() const
Returns DRbar MW, must be calculated by calcDrBarPars first.
Definition: softsusy.cpp:721
virtual void treeUpSquark(DoubleMatrix &mass, double mtrun, double pizztMS, double sinthDRbarMS, int family)
Definition: softsusy.cpp:3360
virtual double piZZTsfermions(double p, double Q) const
Calculates sfermion contribution to the transverse part of Z self-energy:
Definition: softsusy.cpp:8155
DoubleMatrix addSupChargino(double p1, double p2, int family, DoubleMatrix &chargino)
Definition: softsusy.cpp:6379
Complex pis1s2Charginos(double p, double q) const
Definition: softsusy.cpp:11662
DoubleMatrix addSlepHiggs(double p1, double p2, int family, DoubleMatrix &higgs)
Definition: softsusy.cpp:4518
virtual void addChaLoopGauge(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS, DoubleMatrix b1pCha, DoubleMatrix b0pCha, DoubleMatrix b1pNeut, DoubleMatrix b0pNeut) const
Definition: softsusy.cpp:1840
const MssmSoftsusy & operator=(const MssmSoftsusy &s)
Set all data in the object equal to another.
Definition: softsusy.cpp:34
virtual MssmSusyRGE guessAtSusyMt(double tanb, const QedQcd &oneset)
Definition: softsusy.cpp:7265
double piHpHmGauge(double p, double q) const
LCT: Returns pure gauge contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8683
void ewsbConditions(DoubleVector &values) const
Definition: softsusy.cpp:890
void assignHiggs(DoubleVector &higgsm, DoubleVector &higgsc, DoubleVector &dnu, DoubleVector &dnd, DoubleVector &cn, double beta) const
Definition: softsusy.cpp:4482
void setDrBarPars(const drBarPars &s)
Sets tree-level DRbar parameters.
Definition: softsusy.h:262
void setMinpot(double)
Sets potential value at minimum of Higgs potential.
Definition: softsusy.cpp:188
void flagTachyon(tachyonType a)
Flags a negative-mass squared scalar (really a CCB problem)
Definition: softsusy.h:207
virtual void addStopCorrection(double p, DoubleMatrix &mass, double mt)
As in BPMZ appendix.
Definition: softsusy.cpp:4431
virtual void addSupCorrection(DoubleMatrix &mass, int family)
As in BPMZ appendix, INCLUDING weak boson loops.
Definition: softsusy.cpp:6432
virtual void higgsMSLHA(ostream &out)
higgs part of mass block of SLHA
Definition: softsusy.cpp:10053
void headerSLHA(ostream &out)
Definition: softsusy.cpp:9980
double displayTadpole2Ms() const
displays t_2/v_2 tadpole
Definition: softsusy.cpp:176
bool displaySetTbAtMX() const
Is tan beta set at the user defined SUSY breaking scale?
Definition: softsusy.h:191
double piHpHmGauginos(double p, double q) const
Definition: softsusy.cpp:8766
virtual double calcRunMtStopGluino() const
Obtains (16 Pi^2) / mt times 1-loop stop-Gluino corrections.
Definition: softsusy.cpp:2662
virtual void hmixSLHA(ostream &out)
higgs part of mixing blocks
Definition: softsusy.cpp:10253
void useAlternativeMt()
Use alernative mt(MSSM,DR-bar) calculation.
Definition: softsusy.h:268
virtual double calcRunMtNeutralinos() const
Obtains (16 Pi^2) / mt times 1-loop Neutralino corrections.
Definition: softsusy.cpp:2740
void assignHiggsSfermions(DoubleVector &higgsm, DoubleVector &higgsc, DoubleVector &dnu, DoubleVector &dnd, DoubleVector &cn, double beta) const
some switches due to BPMZ's different conventions
Definition: softsusy.cpp:4504
void setMw(double)
sets pole MW prediction
Definition: softsusy.cpp:190
double doCalcTad1Neutralinos(double q, double costhDRbar, double g, double cosb) const
one loop H1 tadpole contributions from Neutralinos in the loops
Definition: softsusy.cpp:586
void modselSLHA(ostream &out, const char model[])
This does the MODSEL block of SLHA.
Definition: softsusy.cpp:10184
void setFracDiff(double fD)
sets fracDiff, needed for access by NmssmSoftsusy methods
Definition: softsusy.h:304
virtual double calcRunMtauNeutralinos(double mTauSMMZ) const
Definition: softsusy.cpp:3226
Complex pis2s2Charginos(double p, double q) const
Definition: softsusy.cpp:11693
double doCalcTad2Sfermions(DoubleMatrix lTS2Lr, DoubleMatrix lBS2Lr, DoubleMatrix lTauS2Lr, double costhDRbar) const
routine to calculate sfermiom contributions to (16 \pi^2) t1 / v1
Definition: softsusy.cpp:460
Contains all supersymmetric RPC-MSSM parameters.
Definition: susy.h:218
Contains all supersymmetric RPC-MSSM parameters.
Definition: susy.h:71
MssmSusy beta(sBrevity &) const
Definition: susy.cpp:300
const MssmSusy & displayMssmSusy() const
Returns whole object as a const.
Definition: susy.h:267
Quark and lepton masses and gauge couplings in QEDxQCD effective theory.
Definition: lowe.h:55
Describes a set of parameters and RGEs in general.
Definition: rge.h:49
double displayMu() const
Return renomalisation scale.
Definition: rge.h:67
double m1
decay global variable declarations
Definition: decays.cpp:14
switches (options) and parameters such as default fermion masses, required accuracy etc
Interface to higher loop corrections.
QedQcd object contains Standard Model quark and lepton masses. It integrates them using 3 loop qcd x ...
A few handy routines for the MSSM: SUSY breaking conditions etc.
Provides routines for calculating Jacobian fine-tuning.
global variable declaration
Definition: def.cpp:13
mass
used to give order of quark masses stored
Definition: lowe.h:46
DoubleVector boundaryCondition(double m0, double m12, double a0, const RpvSusyPars &r)
Definition: rpvsusypars.cpp:373
double rho2(double r)
Two-loop Standard Model corrections to rho parameter.
Definition: softsusy.cpp:9239
EPoleUncertainties
uncertainty flags
Definition: softsusy.h:41
double fEff(double x)
function used for calculating sin theta_eff
Definition: softsusy.cpp:9259
double lnLHiggs(double mh)
returns the smeared log likelihood coming from LEP2 Higgs mass bounds
Definition: softsusy.cpp:10711
trilinears
SUSY breaking trilinear parameters.
Definition: softpars.h:25
void printShortInitialise()
Prints out header line for print-short output.
Definition: softsusy.cpp:1530
double gEff(double x)
function used for calculating sin theta_eff
Definition: softsusy.cpp:9266
double MZ
global pole mass of MZ in GeV - MZCENT is defined in def.h
Definition: def.cpp:31
double sw2
Some parameters used in the computation.
Definition: def.cpp:51
int PRINTOUT
no verbose debug output
Definition: def.cpp:15
DoubleVector mhIntegrand(double mh, const DoubleVector &)
Definition: softsusy.cpp:10702
double minimufb3(double lnH2)
Definition: softsusy.cpp:6960
istream & operator>>(istream &left, flavourPhysical &s)
Formatted input of physical parameters.
Definition: flavoursoft.cpp:1220
double lep2Likelihood(double mh)
Fit to LEP2 Standard Model results.
Definition: softsusy.cpp:10691
numerical routines - differential equation solver, differentiator and function minimiser for instance
Flags for potential problems in sProblem structure, and structure for containing physical MSSM parame...
Soft SUSY breaking parameters.
DRbar values of masses and mixings in MSSM.
Definition: physpars.h:120
Contains data needed in beta function calculation to make it faster.
Definition: susy.h:31
Masses of the physical particles.
Definition: physpars.h:32
Various boolean values to flag any problems with the parameter point.
Definition: physpars.h:75
bool muSqWrongSign
mu^2 came out with wrong sign; no REWSB
Definition: physpars.h:86
tachyonType tachyonWarning
Definition: physpars.h:85
bool irqfp
Infra-red quasi fixed point breached.
Definition: physpars.h:78
bool notGlobalMin
Not in global minimum of Higgs potential.
Definition: physpars.h:89
tachyonType tachyon
Definition: physpars.h:81
bool higgsUfb
Higgs potential inconsistent with a good minimum.
Definition: physpars.h:88
bool inaccurateHiggsMass
Higgs mass is potentially inaccurate and cant be trusted.
Definition: physpars.h:93
bool noConvergence
Iteration did not converge: not always serious.
Definition: physpars.h:80
bool noRhoConvergence
Couldn't calculate electroweak rho parameter.
Definition: physpars.h:79
bool noMuConvergence
Definition: physpars.h:91
bool m3sq
m3sq came out with wrong sign; no REWSB
Definition: physpars.h:87
bool badConvergence
Nowhere near a decent solution.
Definition: physpars.h:77
bool nonperturbative
Running went non-perturbative.
Definition: physpars.h:90
const double accuracy
Approximate accuracy with which 3 body decays are calculated.
Definition: threeBodyDecays.h:38
Contains two-loop routines from Slavich et al.
bool close(double m1, double m2, double tol)
Definition: utils.cpp:68
A few handy bits and pieces - little mathematical functions and the like.