html/softsusy_8h_source.html

 #ifndef SOFTSUSY_H

 #define SOFTSUSY_H


 #include <iostream>

 #include <fstream>

 #include <sstream>

 #include <cstring>

 #include <cstdlib>

 #include <cmath>

 #include <limits>

 #include "def.h"

 #include "utils.h"

 #include "numerics.h"

 #include "physpars.h"

 #include "lowe.h"

 #include "softpars.h"

 #include "twoloophiggs.h"

 #include "mssmUtils.h"

 #include "higher_order.h"

 #include "mssmjacobian.h"


 #define HR "----------------------------------------------------------"


 namespace softsusy {


   enum EPoleUncertainties {

      DeltaQpole = 1,

      DeltaQmatch = 2,

      DeltaMt = 4,

      DeltaAlphaS = 8,

      DeltaAlphaEm = 16,

      DeltaAll = DeltaQpole + DeltaQmatch + DeltaMt + DeltaAlphaS + DeltaAlphaEm

   };


   const double mxDefault = 1.9e16;


   class AltEwsbMssm {

   private:

     double mAcond, muCond;

   public:

      void setAltEwsb(double ma, double mu) { mAcond = ma; muCond = mu; };

      void setAltEwsbMssm(const AltEwsbMssm & s) { *this = s; };


     //   int displayConditionStyle() { return conditionStyle; };

      double displayMaCond() const { return mAcond; };

      double displayMuCond() const { return muCond; };

      const AltEwsbMssm & displayAltEwsbMssm() const { return *this; };


      void setMaCond(double maInput) { mAcond = maInput;  };

      void setMuCond(double muInput) { muCond = muInput;  };

   };


   class MssmSoftsusy: public MssmSusy, public MssmSoftPars, public AltEwsbMssm,

                       public Approx, public RGE {

   private:

     sPhysical physpars;

     drBarPars forLoops;

     sProblem problem;

     double msusy;

     double minV;

     double mw;

     QedQcd dataSet;

     double fracDiff;

     bool setTbAtMX;

     bool altEwsb;

     bool altMt;

     bool altAlphaS;

     bool altAlphaEm;

     double predMzSq;

     double t1OV1Ms, t2OV2Ms;

     double t1OV1Ms1loop, t2OV2Ms1loop;

     double qewsb;


   protected:

     void setT1OV1Ms(double t1) { t1OV1Ms = t1; }

     void setT2OV2Ms(double t2) { t2OV2Ms = t2; }

     void setT1OV1Ms1loop(double t1) { t1OV1Ms1loop = t1; }

     void setT2OV2Ms1loop(double t2) { t2OV2Ms1loop = t2; }

     double mxBC;


     void check_flags();

     DoubleMatrix calcHiggs3L(bool is_bottom);

   public:

     typedef void (*TMSSMBoundaryCondition)(MssmSoftsusy&, const DoubleVector&);


     int included_thresholds;


     //  void (*boundaryCondition)(Softsusy &, const DoubleVector &);

     MssmSoftsusy();

     MssmSoftsusy(const MssmSusyRGE &);


     MssmSoftsusy(const MssmSoftsusy &);

     MssmSoftsusy(const MssmSusy &, const MssmSoftPars & s,

                  const sPhysical & sp, double mu, int l,

                  int t, double hv);


     const MssmSoftsusy & operator=(const MssmSoftsusy & s);


     virtual const DoubleVector display() const;


      const MssmSoftsusy & displayMssmSoft() const;


      const sPhysical & displayPhys() const;


     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;

     sPhysical displayPhysUncertaintyScaleVariation() const;


      const drBarPars & displayDrBarPars() const;

     const sProblem & displayProblem() const {return problem; };

      const QedQcd & displayDataSet() const;

     double displaySoftA(trilinears, int, int) const;


      double displayFracDiff() const { return fracDiff; };

     double displayMinpot() const;

     double displayMsusy() const;

     double displayMw() const;

     double displayMwRun() const;

     double displayMzRun() const;

     double displayTadpole1Ms() const;

     double displayTadpole2Ms() const;

     double displayTadpole1Ms1loop() const;

     double displayTadpole2Ms1loop() const;

     double displayMxBC() const { return mxBC; };

     const MssmSoftsusy & displaySoftsusy() const { return *this; }

     double displayMz() const { return double(MZ); }

     bool displaySetTbAtMX() const { return setTbAtMX; }

     bool displayAltEwsb() const { return altEwsb; }

     bool displayAltMt() const { return altMt; }

     bool displayAltAlphaS() const { return altAlphaS; }

     bool displayAltAlphaEm() const { return altAlphaEm; }

     double displayPredMzSq() const { return predMzSq; }

     double displayQewsb() const { return qewsb; }

     double displayMatchingScale() const { return displayDataSet().displayMu(); }


     void flagMgutOutOfBounds(bool a) { problem.mgutOutOfBounds = a; };

     void flagIrqfp(bool a) { problem.irqfp = a; };

     void flagNonperturbative(bool a) { problem.nonperturbative = a; };

     void flagTachyon(tachyonType a) {

       if (close(displayMu(), MZ, 1.0e-6)) problem.tachyonWarning = a;

       else { problem.tachyon = a;

         if (PRINTOUT > 2) cout << tachyonNames[a] << " tachyon ";

       }

     };

     void flagTachyonWarning(tachyonType a) { problem.tachyonWarning = a; }

     void flagM3sq(bool a) { problem.m3sq = a; };

     void flagNoConvergence(bool a) { problem.noConvergence = a; };

     void flagNoMuConvergence(bool a) { problem.noMuConvergence = a; };

     void flagNoRhoConvergence(bool a) { problem.noRhoConvergence = a; };

     void flagMusqwrongsign(bool a) { problem.muSqWrongSign = a; };

     void flagInaccurateHiggsMass(bool a) { problem.inaccurateHiggsMass = a; };

     void flagHiggsufb(bool a) { problem.higgsUfb = a; };

     void flagNotGlobalMin(bool a) { problem.notGlobalMin = a; };

     void flagAllProblems(bool a, tachyonType b) { problem.irqfp = a;

       problem.tachyon = b; problem.tachyonWarning = b;

       problem.m3sq = a; problem.badConvergence = a;

       problem.noConvergence = a; problem.higgsUfb = a; problem.notGlobalMin = a;

       problem.nonperturbative = a; problem.noRhoConvergence = a;

       problem.noMuConvergence = a; problem.muSqWrongSign = a;

       problem.inaccurateHiggsMass = b; problem.mgutOutOfBounds = a; }

     void flagProblemThrown(bool a) { problem.problemThrown = a; }

     void setProblem(const sProblem a) { problem = a; }


     void setSoftsusy(const MssmSoftsusy & s) { *this = s; };

     void setData(const QedQcd & r) { dataSet = r; };

     void setMinpot(double);

     void setMsusy(double);

     void setMw(double);

     void setPhys(const sPhysical & s) { physpars = s; };

     void setMxBC(double mx) { mxBC = mx; };

     void setDrBarPars(const drBarPars & s) { forLoops = s; };

     void setSetTbAtMX(bool a) { setTbAtMX = a; }

     void useAlternativeEwsb() { altEwsb = true; }

     void useAlternativeMt() { altMt = true; }

     void useAlternativeAlphaS() { altAlphaS = true; }

     void useAlternativeAlphaEm() { altAlphaEm = true; }

     void setPredMzSq(double a) { predMzSq = a; }

     void set(const DoubleVector &);

     void setQewsb(double q) { qewsb = q; };


     void setTwoLoopAlphasThresholds(bool sw);

     void setTwoLoopMtThresholds(bool sw);

     void setTwoLoopMbMtauThresholds(bool sw);

     void setAllTwoLoopThresholds(bool sw);


     DoubleVector beta() const {

       sBrevity a;

       DoubleVector y(MssmSusy::beta(a).display());

       DoubleVector x(MssmSoftPars::beta(displayMssmSusy()));

       int i; for (i=1; i<=y.displayEnd(); i++) x(i) = y(i);

       return x;

     };


     void setFracDiff(double fD) { fracDiff = fD; };

     //for use in doCalcTadpole1oneLoop

     void H1SfSfCouplings(DoubleMatrix & lTS1Lr, DoubleMatrix & lBS1Lr,

                          DoubleMatrix & lTauS1Lr, double gmzOcthW, double mu,

                          double cosb, double v1) const;

     //for use in doCalcTadpole1oneLoop

     void H2SfSfCouplings(DoubleMatrix & lTS2Lr, DoubleMatrix & lBS2Lr,

                          DoubleMatrix & lTauS2Lr, double gmzOcthW, double mu,

                          double sinb) const;

     double doCalcTad1Sfermions(DoubleMatrix lTS1Lr,  DoubleMatrix lBS1Lr,

                                DoubleMatrix lTauS1Lr, double costhDRbar) const;

     double doCalcTad2Sfermions(DoubleMatrix lTS2Lr, DoubleMatrix lBS2Lr,

                                DoubleMatrix lTauS2Lr, double costhDRbar) const;

     //for use in doCalcTadpole1oneLoop

     double doCalcTad1fermions(double q, double v1) const;

     //for use in doCalcTadpole1oneLoop

     double doCalcTad2fermions(double q) const;

     // Follwing BPMZ Goldstone bosons are not included in this.

     double doCalcTad1Higgs(double q, double costhDRbar2, double g,

                            double tanb) const;

     //one loop H2 tadpole contributions from Higgs bosons in the loops

     // Follwing BPMZ Goldstone bosons are not included in this.

     double doCalcTad2Higgs(double q, double costhDRbar2, double g,

                            double tanb) const;

     double doCalcTad1Neutralinos(double q, double costhDRbar, double g,

                                  double cosb) const;

     double doCalcTad2Neutralinos(double q, double costhDRbar, double g,

                                  double sinb) const;

     double doCalcTad1Charginos(double q, double g, double cosb) const;


     double doCalcTad2Charginos(double q, double g, double sinb) const;


     double doCalcTad1GaugeBosons(double q, double costhDRbar2, double g,

                                  double tanb) const;


     double doCalcTad2GaugeBosons(double q, double costhDRbar2,

                                  double g, double tanb) const;

     void doTadpoles(double mt, double sinthDRbar);

     virtual double doCalcTadpole1oneLoop(double mt, double sinthDRbar);

     virtual double doCalcTadpole2oneLoop(double mt, double sinthDRbar);

     virtual void calcTadpole1Ms1loop(double mt, double sinthDRbar);

     virtual void calcTadpole2Ms1loop(double mt, double sinthDRbar);

     DoubleMatrix addStopQCD(double p, double mt, DoubleMatrix & strong);

     DoubleMatrix addStopStop(double p, double mt, DoubleMatrix & stop);

     DoubleMatrix addStopSbottom(double p, double mt, DoubleMatrix & sbottom);

     DoubleMatrix addStopHiggs(double p, double mt, DoubleMatrix & higgs);

     DoubleMatrix addStopEweak(double p, DoubleMatrix & electroweak);

     DoubleMatrix addStopNeutralino(double p, double mt,

                                    DoubleMatrix & neutralino);

     DoubleMatrix addStopChargino(double p, DoubleMatrix & chargino);

     virtual void addStopCorrection(double p, DoubleMatrix & mass, double mt);

     DoubleMatrix addSdownQCD(double p1, double p2, int family,

                              DoubleMatrix & strong);

     DoubleMatrix addSdownHiggs(double p1, double p2, int family,

                                DoubleMatrix & higgs);

     DoubleMatrix addSdownEweak(double p1, double p2, int family,

                                DoubleMatrix & electroweak);

     DoubleMatrix addSdownNeutralino(double p1, double p2, int family,

                                     DoubleMatrix & neutralino);

     DoubleMatrix addSdownChargino(double p1, double p2, int family,

                                   DoubleMatrix & chargino);

     virtual void addSdownCorrection(DoubleMatrix & mass, int family);

     DoubleMatrix addSbotQCD(double p, double mt, DoubleMatrix & strong);

     void addSbotSfermion(double p, double mt, DoubleMatrix & stop,

                          DoubleMatrix & sbottom);

     DoubleMatrix addSbotHiggs(double p, double mt, DoubleMatrix & higgs);

     DoubleMatrix addSbotEweak(double p, DoubleMatrix & electroweak);

     DoubleMatrix addSbotNeutralino(double p, double mt,

                                    DoubleMatrix & neutralino);

     DoubleMatrix addSbotChargino(double p, double mt, DoubleMatrix & chargino);

     virtual void addSbotCorrection(double p, DoubleMatrix & mass, double mb);

     DoubleMatrix addSlepHiggs(double p1, double p2, int family,

                               DoubleMatrix & higgs);

     DoubleMatrix addSlepEweak(double p1, double p2, int family,

                               DoubleMatrix & electroweak);

     void addSlepGaugino(double p1, double p2, int family,

                         DoubleMatrix & chargino, DoubleMatrix & neutralino);

     virtual void addSlepCorrection(DoubleMatrix & mass, int family);

     void addSlepCorrection(double p, DoubleMatrix & mass, int family);

     void addStauSfermion(double p, double mtau, DoubleMatrix & stop,

                          DoubleMatrix & sbottom);

     void addStauGaugino(double p, double mtau, DoubleMatrix & chargino,

                         DoubleMatrix & neutralino);

     DoubleMatrix addStauEweak(double p, double mtau, DoubleMatrix & electroweak);

     DoubleMatrix addStauHiggs(double p, double mtau, DoubleMatrix & higgs);

     virtual void addStauCorrection(double p, DoubleMatrix & mass, double mtau);

     DoubleMatrix addSupQCD(double p1, double p2, int family,

                            DoubleMatrix & strong);

     DoubleMatrix addSupHiggs(double p1, double p2, int family,

                              DoubleMatrix & higgs);

     DoubleMatrix addSupEweak(double p1, double p2, int family,

                              DoubleMatrix & electroweak);

     DoubleMatrix addSupNeutralino(double p1, double p2, int family,

                                   DoubleMatrix & neutralino);

     DoubleMatrix addSupChargino(double p1, double p2, int family,

                                 DoubleMatrix & chargino);

     virtual void addSupCorrection(DoubleMatrix & mass, int family);

     void addSnuTauSfermion(double p, double & stop, double & sbottom);

     double addSnuTauHiggs(double p, double & higgs);

     double addSnuTauEweak(double p, double & electroweak);

     void addSnuTauGaugino(double p, double & chargino, double & neutralino);

     virtual void addSnuTauCorrection(double & mass);

     double addSnuHiggs(double p, int family, double & higgs);

     double addSnuEweak(double p, int family, double & electroweak);

     void addSnuGaugino(double p, int family, double & chargino,

                        double & neutralino);

     virtual void addSnuCorrection(double & mass, int family);

     void addSquarkCorrection(DoubleMatrix & mass);

     virtual void doUpSquarks(double mt, double pizztMS, double sinthDRbarMS, int

                              accuracy);

     virtual void doDownSquarks(double mb, double pizztMS, double sinthDRbarMS,

                                int accuracy, double mt);

     virtual void doChargedSleptons(double mT, double pizztMS,

                                    double sinthDRbarMS, int accuracy);

     virtual void doSnu(double pizztMS, int accuracy = 0);

     virtual void treeUpSquark(DoubleMatrix & mass, double mtrun, double pizztMS,

                               double sinthDRbarMS, int family);

     virtual void treeDownSquark(DoubleMatrix & mass, double mbrun, double pizztMS,

                                 double sinthDRbarMS, int family);

     virtual void treeChargedSlepton(DoubleMatrix & mass, double mTrun,

                                     double pizztMS, double sinthDRbarMS,

                                     int family);

     void treeSnu(double & mSq, double pizztMS, int family);


     virtual void calcDrBarPars();

     void setNeutCurrCouplings(double sinthDRbar, double & sw2, double & guL,

                               double & gdL, double & geL, double & guR,

                               double & gdR, double & geR );

     void calcDRTrilinears(drBarPars & eg, double vev, double beta);


     void calcDrBarHiggs(double beta, double mz2, double mw2,

                         double sinthDRbar, drBarPars & eg);

     void treeCharginos(DoubleMatrix & mCh, double beta, double mw);

     void treeNeutralinos(DoubleMatrix & mN, double beta, double mz,

                          double mw, double sinth);

     // calculates the chargino and neutralino DRbar parameters.

     //It will fill in the chargino and neutralino

     void calcDrBarGauginos(double beta, double mw, double mz, double sinth,

                            drBarPars & eg);

     virtual void sparticleThresholdCorrections(double tb);

     //                                       const DoubleVector & pars);

     double qedSusythresh(double alphaEm, double Q) const;

     double qcdSusythresh(double alphasMSbar, double Q);

     double calcMs() const;

     virtual void physical(int accuracy);

     //to the pole mt for use in calcRunningMt

     virtual double calcRunMtQCD() const;

     //to the pole mt for use in calcRunningMt

     virtual double calcRunMtStopGluino() const;

     //to the pole mt for use in calcRunningMt

     virtual double calcRunMtHiggs() const;


     //to the pole mt for use in calcRunningMt

     virtual double calcRunMtNeutralinos() const;

     //to the pole mt for use in calcRunningMt

     virtual double calcRunMtCharginos() const;


     virtual double calcRunningMt();

     virtual double calcRunMtauDrBarConv() const;

     // Obtains (1 / mTAU)  times 1-loop squark-chargino corrections

     //to mtau for use in calcRunningMtau

     virtual double calcRunMtauCharginos(double mTauSMMZ) const;

     // Obtains (1 / mTAU)  times 1-loop squark-neutralino corrections

     //to mtau for use in calcRunningMtau

     virtual double calcRunMtauNeutralinos(double mTauSMMZ) const;

     // to mtau for use in calcRunningMtau

     virtual double calcRunMtauHiggs() const;

     virtual double calcRunningMtau() ;


     virtual double calcRunMbDrBarConv() const;

     //to mb for use in calcRunningMb

     virtual double calcRunMbSquarkGluino() const;

     //to mb for use in calcRunningMb

     virtual double calcRunMbChargino() const;

     // to mb for use in calcRunningMb

     virtual double calcRunMbHiggs() const;

     // to mb for use in calcRunningMb

     virtual double calcRunMbNeutralinos() const;

     // rruiz: remove const

     virtual double calcRunningMb();

     /*

    double calcHt(double vev);

    double calcHb(double vev) const;

    double calcHtau(double vev) const;

     */

     double calcSinthdrbar() const;

     double getVev();

     double getVev(double pizzt);

     virtual void charginos(int accuracy, double piwwt);

     virtual void addChaLoopSfermion(double p, DoubleMatrix & sigmaL, DoubleMatrix & sigmaR, DoubleMatrix & sigmaS) const;

     virtual void addChaLoopGauge(double p, DoubleMatrix & sigmaL,

                                  DoubleMatrix & sigmaR, DoubleMatrix & sigmaS,

                                  DoubleMatrix b1pCha, DoubleMatrix b0pCha,

                                  DoubleMatrix b1pNeut,

                                  DoubleMatrix b0pNeut) const;

     virtual void addChaLoopHiggs(double p, DoubleMatrix & sigmaL,

                                  DoubleMatrix & sigmaR, DoubleMatrix & sigmaS,

                                  DoubleMatrix b1pCha, DoubleMatrix b0pCha,

                                  DoubleMatrix b1pNeut,

                                  DoubleMatrix b0pNeut) const;

     virtual void addCharginoLoop(double p, DoubleMatrix &);

     virtual void neutralinos(int accuracy, double piwwt, double pizzt);

     void addNeutLoopSfermion(double p, DoubleMatrix & sigmaL,

                              DoubleMatrix & sigmaR, DoubleMatrix & sigmaS);

     void getNeutPassarinoVeltman(double p, double q, DoubleMatrix & b0fn,

                                  DoubleMatrix & b1fn);

     void addNeutLoopGauge(double p, DoubleMatrix & sigmaL,

                           DoubleMatrix & sigmaR, DoubleMatrix & sigmaS);

     void addNeutLoopHiggs(double p, DoubleMatrix & sigmaL,

                           DoubleMatrix & sigmaR, DoubleMatrix & sigmaS);

     virtual void addNeutralinoLoop(double p, DoubleMatrix &);

     virtual void gluino(int accuracy);

     void calcHiggsAtScale(int accuracy, double & mt, double & sinthDRbar,

                           double & piwwtMS, double & pizztMS, double q = 0.);

     bool higgs(int accuracy, double piwwt, double pizzt);

     //  virtual void newhiggs(int accuracy, double piwwt, double pizzt);

     virtual int rewsbMu(int sgnMu, double & mu) const;

     virtual int rewsbM3sq(double, double &) const;

     int rewsbM3sqTree(double, double &) const;

     void alternativeEwsb(double mt);

     virtual void rewsb(int sgnMu, double mt, const DoubleVector & pars,

                        double muOld = -6.66e66, double eps = 0.);

     virtual void rewsbTreeLevel(int sgnMu);

     double treeLevelMuSq();

     void iterateMu(double & munew, int sgnMu, double mt,

                    int maxTries, double pizztMS, double sinthDRbar, double tol,

                    int  & err);

     double predTanb(double muSusy = -6.66e66) const;

     double predMzsq(double & tanb, double muOld = -6.66e66, double eps = 0.);

     double deltaEW() const;


     double ewsbCondition1TreeLevel() const;

     double ewsbCondition2TreeLevel() const;

     void ewsbConditions(DoubleVector & values) const;

     void predVevs(DoubleVector & vevs, int & err);


     DoubleVector fineTune(void (*boundaryCondition)(MssmSoftsusy &,

                                                     const DoubleVector &),

                           const DoubleVector & bcPars, double MX,

                           bool doTop = false);

     double it1par(int numPar, const DoubleVector & bcPars);

     double ufb3sl(double);

     double realMinMs() const;

     double calcBayesianNaturalness() const;


     //for p=external momentum, Q=renormalisation scale

     virtual double piZZTHiggs(double p, double Q, double thetaWDRbar) const;

     //for p=external momentum, Q=renormalisation scale

     virtual double piZZTsfermions(double p, double Q) const;

     virtual double piZZTfermions(double p, double Q, bool usePoleMt) const;

     virtual double piZZTNeutralinos(double p, double Q, double thetaWDRbar) const;

     virtual double piZZTCharginos(double p, double q, double thetaWDRbar) const;

     virtual double piZZT(double p, double Q, bool usePoleMt = false) const;

     virtual double piWWTHiggs(double p, double q, double thetaWDRbar) const;

     virtual double piWWTfermions(double p, double Q, bool usePoleMt) const;

     virtual double piWWTsfermions(double p, double Q) const;

     virtual double piWWTgauginos(double p, double Q, double thetaWDRbar) const;

     virtual double piWWT(double p, double Q, bool usePoleMt = false) const;

     virtual void getNeutralinoCharginoHpmCoup(ComplexMatrix & apph1,

                                               ComplexMatrix & apph2,

                                               ComplexMatrix & bpph1,

                                               ComplexMatrix & bpph2) const;

     double piHpHmFermions(double p, double q) const;

     double piHpHmSfermions(double p, double q, double mu) const;

     double piHpHmGauge(double p, double q) const;

     double piHpHmHiggs(double p, double q) const;

     double piHpHmGauginos(double p, double q) const;

     virtual double piHpHm(double p, double Q) const;

     double piZGT(double p, double Q) const;

     virtual double piAA(double p, double Q) const;

     //self-energy: for p=external momentum, q=renormalisation scale

     Complex pis1s1Sfermions(double p, double q,  DoubleMatrix ls1tt,

                             DoubleMatrix ls1bb,  DoubleMatrix ls1tautau) const;

     Complex pis1s2Sfermions(double p, double q,  DoubleMatrix ls1tt,

                             DoubleMatrix ls1bb,  DoubleMatrix ls1tautau,

                             DoubleMatrix ls2tt,  DoubleMatrix ls2bb,

                             DoubleMatrix ls2tautau) const;

     Complex pis2s2Sfermions(double p, double q, DoubleMatrix ls2tt,

                             DoubleMatrix ls2bb,  DoubleMatrix ls2tautau) const;

     Complex pis1s1Fermions(double p, double q) const;

     Complex pis2s2Fermions(double p, double q) const;

     Complex pis1s1Higgs(double p, double q) const;

     Complex pis1s2Higgs(double p, double q) const;

     Complex pis2s2Higgs(double p, double q) const;

     Complex pis1s1Neutralinos(double p, double q) const;

     Complex pis1s2Neutralinos(double p, double q) const;

     Complex pis2s2Neutralinos(double p, double q) const;

     Complex pis1s1Charginos(double p, double q) const;

     Complex pis1s2Charginos(double p, double q) const;

     Complex pis2s2Charginos(double p, double q) const;

     Complex pis1s1(double p, double q) const;

     Complex pis1s2(double p, double q) const;

     Complex pis2s2(double p, double q) const;

     double sinSqThetaEff();

     virtual double h1s2Mix();

     virtual void rhohat(double & outrho, double & outsin, double alphaMZDRbar,

                         double pizztMZ, double piwwt0, double piwwtMW,

                         double tol, int maxTries);

     virtual double deltaVb(double outrho, double outsin, double alphaDRbar,

                            double pizztMZ) const;

     virtual double dRho(double outrho, double outsin, double alphaDRbar,

                         double pizztMZ, double piwwtMW);

     virtual double dR(double outrho, double outsin, double alphaDRbar,

                       double pizztMZ, double piwwt0);

     double maxMass() const;

     int lsp(double & mass, int & posi, int & posj) const;

     int nlsp(double & mass, int & posi, int & posj) const;


     string printShort() const;

     string printLong();

     string printLongDrbar();


     virtual void printObj() {

       cout << this->displayMssmSusy() << this->displayMssmSoftPars();

     };


     double thet(double a, double b, double c = 0.0);


     void fixedPointIteration(TMSSMBoundaryCondition boundaryCondition,

                              double mxGuess,

                              const DoubleVector & pars, int sgnMu, double tanb,

                              const QedQcd & oneset, bool gaugeUnification,

                              bool ewsbBCscale =  false);

     double lowOrg(void (*boundaryCondition)

                   (MssmSoftsusy &, const DoubleVector &),

                   double mxGuess,

                   const DoubleVector & pars, int sgnMu, double tanb,

                   const QedQcd & oneset, bool gaugeUnification,

                   bool ewsbBCscale =  false) {

       fixedPointIteration(boundaryCondition, mxGuess, pars, sgnMu, tanb, oneset,

                           gaugeUnification, ewsbBCscale);

       return displayMxBC();

     };


     void itLowsoft(int maxTries, int sgnMu, double tol,

                    double tanb, void (*boundaryCondition)(MssmSoftsusy &,

                                                           const DoubleVector &),

                    const DoubleVector & pars, bool gaugeUnification,

                    bool ewsbBCscale);


     virtual void rpvSet(const DoubleVector & parameters);


     virtual void methodBoundaryCondition(const DoubleVector & pars);


     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;

     void isajetInterface764(const char fname[80]) const;

     void ssrunInterface764Inside(const char fname [80], fstream & ) const;

     void ssrunInterface764(const char fname [80], const char softfname [80]) const;

     void isawigInterface764(const char fnamein [80], const char fnameout [80],

                             const char fnamesoft[80]) const;

     virtual void lesHouchesAccordOutput(ostream & out, const char model[],

                                         const DoubleVector & pars,

                                         int sgnMu, double tanb, double qMax,

                                         int numPoints,

                                         bool ewsbBCscale);

     /*    void slha1(ostream & out, const char model[], const DoubleVector & pars,

                int sgnMu, double tanb, double qMax, int numPoints,

                bool ewsbBCscale);*/

     virtual void setEwsbConditions(const DoubleVector & inputs);

     void headerSLHA(ostream & out);

     virtual void spinfoSLHA(ostream & out);

     void modselSLHA(ostream & out, const char model[]);

     void sminputsSLHA(ostream & out);

     void minparSLHA(ostream & out, const char model [],

                     const DoubleVector & pars, double tanb, int sgnMu,

                     bool ewsbBCscale);

     virtual void extparSLHA(ostream & out, const DoubleVector & pars,

                             bool ewsbBCscale);

     void massSLHA(ostream & out);

     virtual void higgsMSLHA(ostream & out);

     virtual void neutralinoCharginoMSLHA(ostream & out);

     virtual void sfermionsSLHA(ostream & out);

     void sfermionmixSLHA(ostream & out);

     virtual void neutralinoMixingSLHA(ostream & out);

     void inomixingSLHA(ostream & out);

     virtual void softsusySLHA(ostream & out);

     void alphaSLHA(ostream & out);

     virtual void hmixSLHA(ostream & out);

     void gaugeSLHA(ostream & out);

     virtual void yukawasSLHA(ostream & out);

     virtual void msoftSLHA(ostream & out);

     virtual void drbarSLHA(ostream & out, int numPoints, double qMax, int n);


     void doUfb3(double mgut);


     void assignHiggs(DoubleVector & higgsm, DoubleVector & higgsc,

                      DoubleVector & dnu, DoubleVector & dnd,

                      DoubleVector & cn, double beta) const;


     void assignHiggs(DoubleVector & higgsm, DoubleVector & higgsc) const;


     void assignHiggsSfermions(DoubleVector & higgsm, DoubleVector & higgsc,

                               DoubleVector & dnu, DoubleVector & dnd,

                               DoubleVector & cn, double beta) const;


     double smPredictionMW() const;


     double twoLoopGm2(double amu1Loop) const;


     double twoLpMt() const;

     double twoLpMb() const;


     virtual void doQuarkMixing(DoubleMatrix & mDon, DoubleMatrix & mUpq);


     virtual MssmSusyRGE guessAtSusyMt(double tanb, const QedQcd & oneset);

   };


   std::istream& operator>>(std::istream& left, MssmSoftsusy& s);


   void printShortInitialise();

   double rho2(double r);


   double minimufb3(double);


   double gEff(double x);

   double fEff(double x);


   double lep2Likelihood(double mh);

   DoubleVector mhIntegrand(double mh, const DoubleVector & y);

   double lnLHiggs(double mh);


 }


 #endif


ComplexMatrix
matrix of complex double values dimensions (rows x cols)
Definition: linalg.h:608

Complex
drop-in replacement for the original home-grown Complex class
Definition: mycomplex.h:17

DoubleMatrix
Matrix from 1..rows, 1..cols of double values.
Definition: linalg.h:214

DoubleVector
DoubleVector is of variable length, and contains double precision.
Definition: linalg.h:35

DoubleVector::displayEnd
int displayEnd() const
returns end of dimension
Definition: linalg.h:94

softsusy::AltEwsbMssm
A different REWSB condition - given by mu and MA instead of mh1,2.
Definition: softsusy.h:54

softsusy::Approx
Definition: rge.h:23

softsusy::MssmSoftPars
Soft SUSY breaking parameters and beta functions.
Definition: softpars.h:31

softsusy::MssmSoftPars::displayMssmSoftPars
const MssmSoftPars & displayMssmSoftPars() const
Returns whole object as a const.
Definition: softpars.h:66

softsusy::MssmSoftPars::beta
DoubleVector beta(const MssmSusy &) const
Returns double vector containing numerical beta functions of parameters.
Definition: softpars.cpp:2858

softsusy::MssmSoftsusy
Contains all supersymmetric MSSM parameters, incorporating R_p MSSM.
Definition: softsusy.h:74

softsusy::MssmSoftsusy::pis1s2Sfermions
Complex pis1s2Sfermions(double p, double q, DoubleMatrix ls1tt, DoubleMatrix ls1bb, DoubleMatrix ls1tautau, DoubleMatrix ls2tt, DoubleMatrix ls2bb, DoubleMatrix ls2tautau) const
Definition: softsusy.cpp:11056

softsusy::MssmSoftsusy::isajetInterface764
void isajetInterface764(const char fname[80]) const
prints a file into fname which acts as an input to isajet
Definition: softsusy.cpp:9939

softsusy::MssmSoftsusy::twoLoopGm2
double twoLoopGm2(double amu1Loop) const
Definition: softsusy.cpp:10765

softsusy::MssmSoftsusy::pis2s2
Complex pis2s2(double p, double q) const
checked 28.10.02
Definition: softsusy.cpp:11783

softsusy::MssmSoftsusy::addSbotChargino
DoubleMatrix addSbotChargino(double p, double mt, DoubleMatrix &chargino)
Definition: softsusy.cpp:5962

softsusy::MssmSoftsusy::calcDrBarHiggs
void calcDrBarHiggs(double beta, double mz2, double mw2, double sinthDRbar, drBarPars &eg)
Definition: softsusy.cpp:7683

softsusy::MssmSoftsusy::doCalcTad1Sfermions
double doCalcTad1Sfermions(DoubleMatrix lTS1Lr, DoubleMatrix lBS1Lr, DoubleMatrix lTauS1Lr, double costhDRbar) const
routine to calculate sfermiom contributions to H1 tadpole / v1
Definition: softsusy.cpp:392

softsusy::MssmSoftsusy::smPredictionMW
double smPredictionMW() const
checked 22/04/06
Definition: softsusy.cpp:10737

softsusy::MssmSoftsusy::sminputsSLHA
void sminputsSLHA(ostream &out)
SMINPUTS block of SLHA.
Definition: softsusy.cpp:10371

softsusy::MssmSoftsusy::neutralinoMixingSLHA
virtual void neutralinoMixingSLHA(ostream &out)
nmix block output of SLHA
Definition: softsusy.cpp:10134

softsusy::MssmSoftsusy::it1par
double it1par(int numPar, const DoubleVector &bcPars)
Definition: softsusy.cpp:1149

softsusy::MssmSoftsusy::calcRunMtCharginos
virtual double calcRunMtCharginos() const
Obtains (16 Pi^2) / mt times 1-loop Charginos corrections.
Definition: softsusy.cpp:2807

softsusy::MssmSoftsusy::flagNotGlobalMin
void flagNotGlobalMin(bool a)
LCT: Flags problem if not in global minimum of Higgs potential.
Definition: softsusy.h:232

softsusy::MssmSoftsusy::addChaLoopSfermion
virtual void addChaLoopSfermion(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS) const
Definition: softsusy.cpp:1711

softsusy::MssmSoftsusy::dR
virtual double dR(double outrho, double outsin, double alphaDRbar, double pizztMZ, double piwwt0)
Definition: softsusy.cpp:9518

softsusy::MssmSoftsusy::lesHouchesAccordOutput
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

softsusy::MssmSoftsusy::setProblem
void setProblem(const sProblem a)
Sets all problem flags.
Definition: softsusy.h:244

softsusy::MssmSoftsusy::piZZTHiggs
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

softsusy::MssmSoftsusy::H2SfSfCouplings
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

softsusy::MssmSoftsusy::isawigInterface764
void isawigInterface764(const char fnamein[80], const char fnameout[80], const char fnamesoft[80]) const
Definition: softsusy.cpp:9922

softsusy::MssmSoftsusy::fineTune
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

softsusy::MssmSoftsusy::piWWTgauginos
virtual double piWWTgauginos(double p, double Q, double thetaWDRbar) const
Definition: softsusy.cpp:8477

softsusy::MssmSoftsusy::qcdSusythresh
double qcdSusythresh(double alphasMSbar, double Q)
Definition: softsusy.cpp:7031

softsusy::MssmSoftsusy::doCalcTad2fermions
double doCalcTad2fermions(double q) const
fixes trilnear H2-fermion-fermion couplings
Definition: softsusy.cpp:537

softsusy::MssmSoftsusy::addNeutLoopSfermion
void addNeutLoopSfermion(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)
Definition: softsusy.cpp:2129

softsusy::MssmSoftsusy::gaugeSLHA
void gaugeSLHA(ostream &out)
sfermionic part of mixing blocks
Definition: softsusy.cpp:10230

softsusy::MssmSoftsusy::setSoftsusy
void setSoftsusy(const MssmSoftsusy &s)
Sets whole object equal to another
Definition: softsusy.h:247

softsusy::MssmSoftsusy::softsusySLHA
virtual void softsusySLHA(ostream &out)
SOFTSUSY comments in SLHA.
Definition: softsusy.cpp:10017

softsusy::MssmSoftsusy::piZZTfermions
virtual double piZZTfermions(double p, double Q, bool usePoleMt) const
fermion contribution to the Transverse part of Z self-energy:
Definition: softsusy.cpp:8253

softsusy::MssmSoftsusy::setQewsb
void setQewsb(double q)
Sets user-set scale qewsb.
Definition: softsusy.h:278

softsusy::MssmSoftsusy::pis1s1Charginos
Complex pis1s1Charginos(double p, double q) const
Definition: softsusy.cpp:11637

softsusy::MssmSoftsusy::deltaVb
virtual double deltaVb(double outrho, double outsin, double alphaDRbar, double pizztMZ) const
outrho, outsin represent the DRbar values
Definition: softsusy.cpp:9301

softsusy::MssmSoftsusy::flagMgutOutOfBounds
void flagMgutOutOfBounds(bool a)
Flags weird mgut-type problems.
Definition: softsusy.h:201

softsusy::MssmSoftsusy::calcDrBarGauginos
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

softsusy::MssmSoftsusy::calcRunMbChargino
virtual double calcRunMbChargino() const
Obtains (1 / mb) times 1-loop squark-chargino corrections.
Definition: softsusy.cpp:2979

softsusy::MssmSoftsusy::flagNoRhoConvergence
void flagNoRhoConvergence(bool a)
Flags fact that rho parameter sub iteration didn't converge.
Definition: softsusy.h:223

softsusy::MssmSoftsusy::treeDownSquark
virtual void treeDownSquark(DoubleMatrix &mass, double mbrun, double pizztMS, double sinthDRbarMS, int family)
Definition: softsusy.cpp:6544

softsusy::MssmSoftsusy::massSLHA
void massSLHA(ostream &out)
This does the job of the above method, but outputs the Mass block.
Definition: softsusy.cpp:10077

softsusy::MssmSoftsusy::piAA
virtual double piAA(double p, double Q) const
Definition: softsusy.cpp:8819

softsusy::MssmSoftsusy::calcRunMbDrBarConv
virtual double calcRunMbDrBarConv() const
calculates factor to convert to DrBar for mb
Definition: softsusy.cpp:2943

softsusy::MssmSoftsusy::addSdownCorrection
virtual void addSdownCorrection(DoubleMatrix &mass, int family)
Definition: softsusy.cpp:5533

softsusy::MssmSoftsusy::ssrunInterface764
void ssrunInterface764(const char fname[80], const char softfname[80]) const
Definition: softsusy.cpp:9874

softsusy::MssmSoftsusy::calcHiggs3L
DoubleMatrix calcHiggs3L(bool is_bottom)
Higgs 3L corrections.
Definition: softsusy.cpp:11809

softsusy::MssmSoftsusy::addStopNeutralino
DoubleMatrix addStopNeutralino(double p, double mt, DoubleMatrix &neutralino)
Definition: softsusy.cpp:4310

softsusy::MssmSoftsusy::piZZT
virtual double piZZT(double p, double Q, bool usePoleMt=false) const
Transverse part of Z self-energy: has been checked.
Definition: softsusy.cpp:8353

softsusy::MssmSoftsusy::addSbotQCD
DoubleMatrix addSbotQCD(double p, double mt, DoubleMatrix &strong)
Definition: softsusy.cpp:5569

softsusy::MssmSoftsusy::treeSnu
void treeSnu(double &mSq, double pizztMS, int family)
Definition: softsusy.cpp:6772

softsusy::MssmSoftsusy::addStauEweak
DoubleMatrix addStauEweak(double p, double mtau, DoubleMatrix &electroweak)
Definition: softsusy.cpp:4969

softsusy::MssmSoftsusy::addSquarkCorrection
void addSquarkCorrection(DoubleMatrix &mass)
Now these are calculated at the squark scale.
Definition: softsusy.cpp:3388

softsusy::MssmSoftsusy::setMxBC
void setMxBC(double mx)
Sets the scale at which high-scale boundary conditions are applied.
Definition: softsusy.h:260

softsusy::MssmSoftsusy::included_thresholds
int included_thresholds
Flag allowing to choose which two-loop thresholds have to be included.
Definition: softsusy.h:121

softsusy::MssmSoftsusy::addSnuCorrection
virtual void addSnuCorrection(double &mass, int family)
Found+fixed bug 7/09/06. Thanks to J Kersten.
Definition: softsusy.cpp:3926

softsusy::MssmSoftsusy::addSlepCorrection
void addSlepCorrection(double p, DoubleMatrix &mass, int family)

softsusy::MssmSoftsusy::addSbotCorrection
virtual void addSbotCorrection(double p, DoubleMatrix &mass, double mb)
Checked 16.08.2005.
Definition: softsusy.cpp:6020

softsusy::MssmSoftsusy::predVevs
void predVevs(DoubleVector &vevs, int &err)
DH: solves the EWSB conditions for the Higgs VEVs.
Definition: softsusy.cpp:973

softsusy::MssmSoftsusy::inomixingSLHA
void inomixingSLHA(ostream &out)
This does the job of the above method, but outputs the UMIX/VMIX blocks.
Definition: softsusy.cpp:10149

softsusy::MssmSoftsusy::ewsbCondition1TreeLevel
double ewsbCondition1TreeLevel() const
DH: returns the values of the tree-level EWSB conditions.
Definition: softsusy.cpp:854

softsusy::MssmSoftsusy::msoftSLHA
virtual void msoftSLHA(ostream &out)
MSOFT block of SLHA.
Definition: softsusy.cpp:10268

softsusy::MssmSoftsusy::calcRunningMtau
virtual double calcRunningMtau()
Full BPMZ expression.
Definition: softsusy.cpp:3291

softsusy::MssmSoftsusy::displaySoftsusy
const MssmSoftsusy & displaySoftsusy() const
Returns object as a const.
Definition: softsusy.h:187

softsusy::MssmSoftsusy::addSnuEweak
double addSnuEweak(double p, int family, double &electroweak)
Definition: softsusy.cpp:3783

softsusy::MssmSoftsusy::flagIrqfp
void flagIrqfp(bool a)
Flags Infra-red quasi fixed point breach problem.
Definition: softsusy.h:203

softsusy::MssmSoftsusy::piZZTNeutralinos
virtual double piZZTNeutralinos(double p, double Q, double thetaWDRbar) const
Definition: softsusy.cpp:8293

softsusy::MssmSoftsusy::MssmSoftsusy
MssmSoftsusy()
Default constructor fills object with zeroes.
Definition: softsusy.cpp:67

softsusy::MssmSoftsusy::addStopHiggs
DoubleMatrix addStopHiggs(double p, double mt, DoubleMatrix &higgs)
Definition: softsusy.cpp:4017

softsusy::MssmSoftsusy::piWWTsfermions
virtual double piWWTsfermions(double p, double Q) const
Definition: softsusy.cpp:8436

softsusy::MssmSoftsusy::pis1s2Neutralinos
Complex pis1s2Neutralinos(double p, double q) const
Definition: softsusy.cpp:11566

softsusy::MssmSoftsusy::doCalcTad1GaugeBosons
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

softsusy::MssmSoftsusy::printObj
virtual void printObj()
Prints whole object to standard output.
Definition: softsusy.h:1052

softsusy::MssmSoftsusy::piWWT
virtual double piWWT(double p, double Q, bool usePoleMt=false) const
Definition: softsusy.cpp:8515

softsusy::MssmSoftsusy::minparSLHA
void minparSLHA(ostream &out, const char model[], const DoubleVector &pars, double tanb, int sgnMu, bool ewsbBCscale)
MINPAR block of SLHA.
Definition: softsusy.cpp:10456

softsusy::MssmSoftsusy::displayDrBarPars
const drBarPars & displayDrBarPars() const
Definition: softsusy.cpp:164

softsusy::MssmSoftsusy::piHpHm
virtual double piHpHm(double p, double Q) const
LCT: Charged Higgs self-energy.
Definition: softsusy.cpp:8797

softsusy::MssmSoftsusy::calcRunningMt
virtual double calcRunningMt()
Definition: softsusy.cpp:2869

softsusy::MssmSoftsusy::calcHiggsAtScale
void calcHiggsAtScale(int accuracy, double &mt, double &sinthDRbar, double &piwwtMS, double &pizztMS, double q=0.)
Definition: softsusy.cpp:6782

softsusy::MssmSoftsusy::addNeutralinoLoop
virtual void addNeutralinoLoop(double p, DoubleMatrix &)
Adds the loop corrections on to an input tree-level neutralino mass.
Definition: softsusy.cpp:2518

softsusy::MssmSoftsusy::addSbotSfermion
void addSbotSfermion(double p, double mt, DoubleMatrix &stop, DoubleMatrix &sbottom)
Definition: softsusy.cpp:5594

softsusy::MssmSoftsusy::deltaEW
double deltaEW() const
Returns fine-tuning as defined in Baer et al arXiv:1212.2655 eq 1.2.
Definition: softsusy.cpp:1210

softsusy::MssmSoftsusy::sparticleThresholdCorrections
virtual void sparticleThresholdCorrections(double tb)
Definition: softsusy.cpp:7556

softsusy::MssmSoftsusy::getNeutralinoCharginoHpmCoup
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

softsusy::MssmSoftsusy::getVev
double getVev()
Calculates Higgs VEV parameter from gauge couplings and MZ.
Definition: softsusy.cpp:7543

softsusy::MssmSoftsusy::rpvSet
virtual void rpvSet(const DoubleVector &parameters)
Definition: softsusy.cpp:9643

softsusy::MssmSoftsusy::calcMs
double calcMs() const
DRbar pars should be defined already for this.
Definition: softsusy.cpp:7250

softsusy::MssmSoftsusy::displayTadpole2Ms1loop
double displayTadpole2Ms1loop() const
displays t_2/v_2 tadpole@1 loop
Definition: softsusy.cpp:184

softsusy::MssmSoftsusy::doSnu
virtual void doSnu(double pizztMS, int accuracy=0)
Definition: softsusy.cpp:6759

softsusy::MssmSoftsusy::treeCharginos
void treeCharginos(DoubleMatrix &mCh, double beta, double mw)
Definition: softsusy.cpp:7645

softsusy::MssmSoftsusy::calcRunMtauHiggs
virtual double calcRunMtauHiggs() const
Obtains (1 / mtau) times 1-loop Higgs corrections.
Definition: softsusy.cpp:3191

softsusy::MssmSoftsusy::setSetTbAtMX
void setSetTbAtMX(bool a)
Sets the setTbAtMX flag.
Definition: softsusy.h:264

softsusy::MssmSoftsusy::addNeutLoopHiggs
void addNeutLoopHiggs(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)
Definition: softsusy.cpp:2376

softsusy::MssmSoftsusy::doCalcTad1Higgs
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

softsusy::MssmSoftsusy::pis1s1
Complex pis1s1(double p, double q) const
Definition: softsusy.cpp:11719

softsusy::MssmSoftsusy::dRho
virtual double dRho(double outrho, double outsin, double alphaDRbar, double pizztMZ, double piwwtMW)
Definition: softsusy.cpp:9479

softsusy::MssmSoftsusy::displayMz
double displayMz() const
Returns value of pole MZ being used.
Definition: softsusy.h:189

softsusy::MssmSoftsusy::addStauHiggs
DoubleMatrix addStauHiggs(double p, double mtau, DoubleMatrix &higgs)
Definition: softsusy.cpp:4846

softsusy::MssmSoftsusy::flagNoConvergence
void flagNoConvergence(bool a)
Flags fact that calculation hasn't acheived required accuracy.
Definition: softsusy.h:219

softsusy::MssmSoftsusy::pis2s2Fermions
Complex pis2s2Fermions(double p, double q) const
Definition: softsusy.cpp:11252

softsusy::MssmSoftsusy::flagNoMuConvergence
void flagNoMuConvergence(bool a)
Flags fact that mu sub iteration didn't converge.
Definition: softsusy.h:221

softsusy::MssmSoftsusy::displayPhysUncertainty
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

softsusy::MssmSoftsusy::doCalcTadpole2oneLoop
virtual double doCalcTadpole2oneLoop(double mt, double sinthDRbar)
Does the calculation of one-loop pieces of .
Definition: softsusy.cpp:729

softsusy::MssmSoftsusy::rewsbTreeLevel
virtual void rewsbTreeLevel(int sgnMu)
Definition: softsusy.cpp:1424

softsusy::MssmSoftsusy::check_flags
void check_flags()
check consistency of threshold flags
Definition: softsusy.cpp:7322

softsusy::MssmSoftsusy::extparSLHA
virtual void extparSLHA(ostream &out, const DoubleVector &pars, bool ewsbBCscale)
EXTPAR block of SLHA.
Definition: softsusy.cpp:10387

softsusy::MssmSoftsusy::calcRunMtQCD
virtual double calcRunMtQCD() const
Obtains (16 Pi^2) / mt times 1-loop and 2-loop QCD corrections.
Definition: softsusy.cpp:2653

softsusy::MssmSoftsusy::useAlternativeAlphaS
void useAlternativeAlphaS()
Use alernative alpha_s(MSSM,DR-bar) calculation.
Definition: softsusy.h:270

softsusy::MssmSoftsusy::H1SfSfCouplings
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

softsusy::MssmSoftsusy::setData
void setData(const QedQcd &r)
Sets low energy Standard Model fermion mass and gauge coupling data.
Definition: softsusy.h:249

softsusy::MssmSoftsusy::sfermionmixSLHA
void sfermionmixSLHA(ostream &out)
hmix block output of SLHA
Definition: softsusy.cpp:10195

softsusy::MssmSoftsusy::doCalcTad1fermions
double doCalcTad1fermions(double q, double v1) const
fixes trilnear H1-fermion-fermion couplings
Definition: softsusy.cpp:523

softsusy::MssmSoftsusy::displayMinpot
double displayMinpot() const
Returns minimum of Higgs potential.
Definition: softsusy.cpp:168

softsusy::MssmSoftsusy::flagNonperturbative
void flagNonperturbative(bool a)
Flags non-perturbative RG evolution.
Definition: softsusy.h:205

softsusy::MssmSoftsusy::addCharginoLoop
virtual void addCharginoLoop(double p, DoubleMatrix &)
Adds the loop corrections on to an input tree-level chargino mass.
Definition: softsusy.cpp:2042

softsusy::MssmSoftsusy::doChargedSleptons
virtual void doChargedSleptons(double mT, double pizztMS, double sinthDRbarMS, int accuracy)
Definition: softsusy.cpp:6671

softsusy::MssmSoftsusy::addSupEweak
DoubleMatrix addSupEweak(double p1, double p2, int family, DoubleMatrix &electroweak)
Definition: softsusy.cpp:6204

softsusy::MssmSoftsusy::itLowsoft
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

softsusy::MssmSoftsusy::displayMw
double displayMw() const
Definition: softsusy.cpp:170

softsusy::MssmSoftsusy::useAlternativeEwsb
void useAlternativeEwsb()
Use alernative EWSB conditions: set mu and MA(pole)
Definition: softsusy.h:266

softsusy::MssmSoftsusy::doCalcTad2Neutralinos
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

softsusy::MssmSoftsusy::realMinMs
double realMinMs() const
Definition: softsusy.cpp:7500

softsusy::MssmSoftsusy::setPhys
void setPhys(const sPhysical &s)
Sets all physical parameters.
Definition: softsusy.h:258

softsusy::MssmSoftsusy::setTwoLoopMtThresholds
void setTwoLoopMtThresholds(bool sw)
Switch 2-loop threshold  corrections to top mass.
Definition: softsusy.cpp:10896

softsusy::MssmSoftsusy::qedSusythresh
double qedSusythresh(double alphaEm, double Q) const
Definition: softsusy.cpp:7103

softsusy::MssmSoftsusy::pis1s1Higgs
Complex pis1s1Higgs(double p, double q) const
Definition: softsusy.cpp:11261

softsusy::MssmSoftsusy::doCalcTad1Charginos
double doCalcTad1Charginos(double q, double g, double cosb) const
one loop H1 tadpole contributions from Charginos in the loops
Definition: softsusy.cpp:617

softsusy::MssmSoftsusy::printLong
string printLong()
Prints a list of all sparticle/Higgs masses to standard output.
Definition: softsusy.cpp:1578

softsusy::MssmSoftsusy::yukawasSLHA
virtual void yukawasSLHA(ostream &out)
sfermionic part of mixing blocks
Definition: softsusy.cpp:10241

softsusy::MssmSoftsusy::flagM3sq
void flagM3sq(bool a)
Flags problem with Higgs potential minimum.
Definition: softsusy.h:217

softsusy::MssmSoftsusy::addStopChargino
DoubleMatrix addStopChargino(double p, DoubleMatrix &chargino)
Definition: softsusy.cpp:4371

softsusy::MssmSoftsusy::methodBoundaryCondition
virtual void methodBoundaryCondition(const DoubleVector &pars)
Definition: softsusy.cpp:9636

softsusy::MssmSoftsusy::pis2s2Sfermions
Complex pis2s2Sfermions(double p, double q, DoubleMatrix ls2tt, DoubleMatrix ls2bb, DoubleMatrix ls2tautau) const
Definition: softsusy.cpp:11129

softsusy::MssmSoftsusy::addSnuTauHiggs
double addSnuTauHiggs(double p, double &higgs)
Definition: softsusy.cpp:3431

softsusy::MssmSoftsusy::addStauCorrection
virtual void addStauCorrection(double p, DoubleMatrix &mass, double mtau)
Definition: softsusy.cpp:5177

softsusy::MssmSoftsusy::calcRunMtHiggs
virtual double calcRunMtHiggs() const
Obtains (16 Pi^2) / mt times 1-loop Higgs corrections.
Definition: softsusy.cpp:2690

softsusy::MssmSoftsusy::treeChargedSlepton
virtual void treeChargedSlepton(DoubleMatrix &mass, double mTrun, double pizztMS, double sinthDRbarMS, int family)
Definition: softsusy.cpp:6645

softsusy::MssmSoftsusy::rhohat
virtual void rhohat(double &outrho, double &outsin, double alphaMZDRbar, double pizztMZ, double piwwt0, double piwwtMW, double tol, int maxTries)
Definition: softsusy.cpp:9561

softsusy::MssmSoftsusy::addSlepEweak
DoubleMatrix addSlepEweak(double p1, double p2, int family, DoubleMatrix &electroweak)
Definition: softsusy.cpp:4605

softsusy::MssmSoftsusy::addStopQCD
DoubleMatrix addStopQCD(double p, double mt, DoubleMatrix &strong)
Definition: softsusy.cpp:3955

softsusy::MssmSoftsusy::piHpHmFermions
double piHpHmFermions(double p, double q) const
LCT: Returns fermion contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8548

softsusy::MssmSoftsusy::higgs
bool higgs(int accuracy, double piwwt, double pizzt)
Definition: softsusy.cpp:11862

softsusy::MssmSoftsusy::calcTadpole2Ms1loop
virtual void calcTadpole2Ms1loop(double mt, double sinthDRbar)
Definition: softsusy.cpp:761

softsusy::MssmSoftsusy::addSbotHiggs
DoubleMatrix addSbotHiggs(double p, double mt, DoubleMatrix &higgs)
Definition: softsusy.cpp:5631

softsusy::MssmSoftsusy::doUfb3
void doUfb3(double mgut)
Definition: softsusy.cpp:191

softsusy::MssmSoftsusy::gluino
virtual void gluino(int accuracy)
Calculates pole gluino mass to 1-loop SUSY corrections.
Definition: softsusy.cpp:2602

softsusy::MssmSoftsusy::pis2s2Higgs
Complex pis2s2Higgs(double p, double q) const
Definition: softsusy.cpp:11450

softsusy::MssmSoftsusy::lsp
int lsp(double &mass, int &posi, int &posj) const
Definition: softsusy.cpp:7171

softsusy::MssmSoftsusy::flagMusqwrongsign
void flagMusqwrongsign(bool a)
Flags point inconsistent with electroweak symmetry breaking.
Definition: softsusy.h:225

softsusy::MssmSoftsusy::displayMzRun
double displayMzRun() const
Returns DRbar MZ, must be calculated by calcDrBarPars first.
Definition: softsusy.cpp:7640

softsusy::MssmSoftsusy::displayPhysUncertaintyScaleVariation
sPhysical displayPhysUncertaintyScaleVariation() const
Displays uncertainty estimate of pole masses by scale variation only.
Definition: softsusy.cpp:12238

softsusy::MssmSoftsusy::doUpSquarks
virtual void doUpSquarks(double mt, double pizztMS, double sinthDRbarMS, int accuracy)
Definition: softsusy.cpp:6470

softsusy::MssmSoftsusy::rewsb
virtual void rewsb(int sgnMu, double mt, const DoubleVector &pars, double muOld=-6.66e66, double eps=0.)
Definition: softsusy.cpp:1467

softsusy::MssmSoftsusy::doDownSquarks
virtual void doDownSquarks(double mb, double pizztMS, double sinthDRbarMS, int accuracy, double mt)
Definition: softsusy.cpp:6572

softsusy::MssmSoftsusy::displayDataSet
const QedQcd & displayDataSet() const
Gives the low energy Standard Model data set used for the object.
Definition: softsusy.cpp:160

softsusy::MssmSoftsusy::addSdownChargino
DoubleMatrix addSdownChargino(double p1, double p2, int family, DoubleMatrix &chargino)
Definition: softsusy.cpp:5487

softsusy::MssmSoftsusy::calcRunMtauCharginos
virtual double calcRunMtauCharginos(double mTauSMMZ) const
Definition: softsusy.cpp:3149

softsusy::MssmSoftsusy::setPredMzSq
void setPredMzSq(double a)
Set MZ^2 predicted after iteration.
Definition: softsusy.h:274

softsusy::MssmSoftsusy::predTanb
double predTanb(double muSusy=-6.66e66) const
Definition: softsusy.cpp:264

softsusy::MssmSoftsusy::charginos
virtual void charginos(int accuracy, double piwwt)
checked
Definition: softsusy.cpp:2090

softsusy::MssmSoftsusy::calcRunMbNeutralinos
virtual double calcRunMbNeutralinos() const
Obtains (1 / mb) times 1-loop neutralino-suqrk corrections.
Definition: softsusy.cpp:12168

softsusy::MssmSoftsusy::display
virtual const DoubleVector display() const
Return contents of object in a vector: for RG evolution.
Definition: softsusy.cpp:803

softsusy::MssmSoftsusy::addStopEweak
DoubleMatrix addStopEweak(double p, DoubleMatrix &electroweak)
Definition: softsusy.cpp:4182

softsusy::MssmSoftsusy::piZZTCharginos
virtual double piZZTCharginos(double p, double q, double thetaWDRbar) const
Definition: softsusy.cpp:8323

softsusy::MssmSoftsusy::ufb3sl
double ufb3sl(double)
Definition: softsusy.cpp:7003

softsusy::MssmSoftsusy::addSupQCD
DoubleMatrix addSupQCD(double p1, double p2, int family, DoubleMatrix &strong)
Definition: softsusy.cpp:6057

softsusy::MssmSoftsusy::drbarSLHA
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

softsusy::MssmSoftsusy::pis1s1Sfermions
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

softsusy::MssmSoftsusy::lowOrg
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

softsusy::MssmSoftsusy::displaySoftA
double displaySoftA(trilinears, int, int) const
Return a trilinear element in "SUGRA style".
Definition: softsusy.cpp:837

softsusy::MssmSoftsusy::treeNeutralinos
void treeNeutralinos(DoubleMatrix &mN, double beta, double mz, double mw, double sinth)
Definition: softsusy.cpp:7652

softsusy::MssmSoftsusy::addSnuTauCorrection
virtual void addSnuTauCorrection(double &mass)
Definition: softsusy.cpp:3671

softsusy::MssmSoftsusy::addSlepCorrection
virtual void addSlepCorrection(DoubleMatrix &mass, int family)
16.09.05 checked.
Definition: softsusy.cpp:4779

softsusy::MssmSoftsusy::treeLevelMuSq
double treeLevelMuSq()
Definition: softsusy.cpp:1419

softsusy::MssmSoftsusy::pis1s2Higgs
Complex pis1s2Higgs(double p, double q) const
Definition: softsusy.cpp:11346

softsusy::MssmSoftsusy::spinfoSLHA
virtual void spinfoSLHA(ostream &out)
This does the SPINFO block of SLHA.
Definition: softsusy.cpp:9990

softsusy::MssmSoftsusy::calcTadpole1Ms1loop
virtual void calcTadpole1Ms1loop(double mt, double sinthDRbar)
checked
Definition: softsusy.cpp:713

softsusy::MssmSoftsusy::piWWTHiggs
virtual double piWWTHiggs(double p, double q, double thetaWDRbar) const
Definition: softsusy.cpp:8378

softsusy::MssmSoftsusy::useAlternativeAlphaEm
void useAlternativeAlphaEm()
Use alernative alpha_em(MSSM,DR-bar) calculation.
Definition: softsusy.h:272

softsusy::MssmSoftsusy::addSnuTauSfermion
void addSnuTauSfermion(double p, double &stop, double &sbottom)
Definition: softsusy.cpp:3417

softsusy::MssmSoftsusy::setMsusy
void setMsusy(double)
Definition: softsusy.cpp:189

softsusy::MssmSoftsusy::isajetNumbers764
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

softsusy::MssmSoftsusy::addSupNeutralino
DoubleMatrix addSupNeutralino(double p1, double p2, int family, DoubleMatrix &neutralino)
Definition: softsusy.cpp:6323

softsusy::MssmSoftsusy::doCalcTad2Higgs
double doCalcTad2Higgs(double q, double costhDRbar2, double g, double tanb) const
Definition: softsusy.cpp:566

softsusy::MssmSoftsusy::displayTadpole1Ms1loop
double displayTadpole1Ms1loop() const
displays t_1/v_1 tadpole@1 loop
Definition: softsusy.cpp:180

softsusy::MssmSoftsusy::addChaLoopHiggs
virtual void addChaLoopHiggs(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS, DoubleMatrix b1pCha, DoubleMatrix b0pCha, DoubleMatrix b1pNeut, DoubleMatrix b0pNeut) const
Definition: softsusy.cpp:1909

softsusy::MssmSoftsusy::piHpHmSfermions
double piHpHmSfermions(double p, double q, double mu) const
LCT: Returns sfermion contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8569

softsusy::MssmSoftsusy::addStauGaugino
void addStauGaugino(double p, double mtau, DoubleMatrix &chargino, DoubleMatrix &neutralino)
Definition: softsusy.cpp:5085

softsusy::MssmSoftsusy::alphaSLHA
void alphaSLHA(ostream &out)
sfermionic part of mixing blocks
Definition: softsusy.cpp:10127

softsusy::MssmSoftsusy::sinSqThetaEff
double sinSqThetaEff()
Calculates sin^2 theta^l_eff.
Definition: softsusy.cpp:9273

softsusy::MssmSoftsusy::displayMsusy
double displayMsusy() const
Returns Higgs minimisation scale.
Definition: softsusy.cpp:169

softsusy::MssmSoftsusy::addNeutLoopGauge
void addNeutLoopGauge(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)
Definition: softsusy.cpp:2317

softsusy::MssmSoftsusy::calcRunningMb
virtual double calcRunningMb()
Definition: softsusy.cpp:3086

softsusy::MssmSoftsusy::mxBC
double mxBC
Scale at which SUSY breaking boundary conditions set.
Definition: softsusy.h:113

softsusy::MssmSoftsusy::neutralinos
virtual void neutralinos(int accuracy, double piwwt, double pizzt)
mixNeut set to diagonal = mixNeut^T mNeutralino mixNeut: checked
Definition: softsusy.cpp:2540

softsusy::MssmSoftsusy::calcBayesianNaturalness
double calcBayesianNaturalness() const
Calculates Bayesian naturalness a la arXiv:1709.07895.
Definition: softsusy.cpp:7512

softsusy::MssmSoftsusy::flagAllProblems
void flagAllProblems(bool a, tachyonType b)
Sets all problems equal to either true or false (contained in a)
Definition: softsusy.h:234

softsusy::MssmSoftsusy::addSlepGaugino
void addSlepGaugino(double p1, double p2, int family, DoubleMatrix &chargino, DoubleMatrix &neutralino)
Definition: softsusy.cpp:4707

softsusy::MssmSoftsusy::setTwoLoopAlphasThresholds
void setTwoLoopAlphasThresholds(bool sw)
Definition: softsusy.cpp:10888

softsusy::MssmSoftsusy::calcDrBarPars
virtual void calcDrBarPars()
Organises tree-level calculation of all sparticle masses and mixings.
Definition: softsusy.cpp:7815

softsusy::MssmSoftsusy::addSdownHiggs
DoubleMatrix addSdownHiggs(double p1, double p2, int family, DoubleMatrix &higgs)
Definition: softsusy.cpp:5228

softsusy::MssmSoftsusy::printLongDrbar
string printLongDrbar()
Prints a list of all running sparticle/Higgs masses to standard output.
Definition: softsusy.cpp:1643

softsusy::MssmSoftsusy::piZGT
double piZGT(double p, double Q) const
Definition: softsusy.cpp:9140

softsusy::MssmSoftsusy::h1s2Mix
virtual double h1s2Mix()
gets h1 mixing element with Hu.
Definition: softsusy.cpp:9475

softsusy::MssmSoftsusy::flagHiggsufb
void flagHiggsufb(bool a)
Flags an inconsistent Higgs minimum.
Definition: softsusy.h:230

softsusy::MssmSoftsusy::setNeutCurrCouplings
void setNeutCurrCouplings(double sinthDRbar, double &sw2, double &guL, double &gdL, double &geL, double &guR, double &gdR, double &geR)
Definition: softsusy.cpp:7754

softsusy::MssmSoftsusy::piHpHmHiggs
double piHpHmHiggs(double p, double q) const
LCT: Returns Higgs contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8699

softsusy::MssmSoftsusy::setEwsbConditions
virtual void setEwsbConditions(const DoubleVector &inputs)
Again, another dummy - useful in alternative EWSB conditions sometimes.
Definition: softsusy.cpp:10883

softsusy::MssmSoftsusy::displayProblem
const sProblem & displayProblem() const
Returns any problem flags associated with the object.
Definition: softsusy.h:166

softsusy::MssmSoftsusy::nlsp
int nlsp(double &mass, int &posi, int &posj) const
Definition: softsusy.cpp:10567

softsusy::MssmSoftsusy::pis2s2Neutralinos
Complex pis2s2Neutralinos(double p, double q) const
Definition: softsusy.cpp:11605

softsusy::MssmSoftsusy::addSnuHiggs
double addSnuHiggs(double p, int family, double &higgs)
Definition: softsusy.cpp:3703

softsusy::MssmSoftsusy::addSbotNeutralino
DoubleMatrix addSbotNeutralino(double p, double mt, DoubleMatrix &neutralino)
Definition: softsusy.cpp:5900

softsusy::MssmSoftsusy::pis1s2
Complex pis1s2(double p, double q) const
Definition: softsusy.cpp:11750

softsusy::MssmSoftsusy::ssrunInterface764Inside
void ssrunInterface764Inside(const char fname[80], fstream &) const
Definition: softsusy.cpp:9884

softsusy::MssmSoftsusy::calcRunMtauDrBarConv
virtual double calcRunMtauDrBarConv() const
calculates factor to convert to DrBar for mtau
Definition: softsusy.cpp:3143

softsusy::MssmSoftsusy::doQuarkMixing
virtual void doQuarkMixing(DoubleMatrix &mDon, DoubleMatrix &mUpq)
input diagonal matrices and it'll give you back mixed ones
Definition: softsusy.cpp:10938

softsusy::MssmSoftsusy::calcSinthdrbar
double calcSinthdrbar() const
Calculates DRbar sin theta_w at the current scale from gauge couplings.
Definition: softsusy.cpp:7518

softsusy::MssmSoftsusy::displayTadpole1Ms
double displayTadpole1Ms() const
displays t_1/v_1 tadpole
Definition: softsusy.cpp:172

softsusy::MssmSoftsusy::getNeutPassarinoVeltman
void getNeutPassarinoVeltman(double p, double q, DoubleMatrix &b0fn, DoubleMatrix &b1fn)
LCT: Returns matrix of Passarino-Veltman B0 and B1 functions.
Definition: softsusy.cpp:2282

softsusy::MssmSoftsusy::flagInaccurateHiggsMass
void flagInaccurateHiggsMass(bool a)
Definition: softsusy.h:228

softsusy::MssmSoftsusy::calcRunMbSquarkGluino
virtual double calcRunMbSquarkGluino() const
Obtains (1 / mb) times 1-loop sbottom-Gluino corrections.
Definition: softsusy.cpp:2956

softsusy::MssmSoftsusy::addSupHiggs
DoubleMatrix addSupHiggs(double p1, double p2, int family, DoubleMatrix &higgs)
Definition: softsusy.cpp:6075

softsusy::MssmSoftsusy::piWWTfermions
virtual double piWWTfermions(double p, double Q, bool usePoleMt) const
Definition: softsusy.cpp:8413

softsusy::MssmSoftsusy::rewsbM3sq
virtual int rewsbM3sq(double, double &) const
returns 1 if mu < 1.0e-9
Definition: softsusy.cpp:221

softsusy::MssmSoftsusy::set
void set(const DoubleVector &)
Sets total set of RGE parameters equal to elements of a vector.
Definition: softsusy.cpp:830

softsusy::MssmSoftsusy::addSdownQCD
DoubleMatrix addSdownQCD(double p1, double p2, int family, DoubleMatrix &strong)
Definition: softsusy.cpp:5211

softsusy::MssmSoftsusy::beta
DoubleVector beta() const
Returns double vector containing numerical beta functions of parameters.
Definition: softsusy.h:295

softsusy::MssmSoftsusy::alternativeEwsb
void alternativeEwsb(double mt)
Fixes mH1^2 and mH2^2 in order to get muCond and mAcond correct.
Definition: softsusy.cpp:1339

softsusy::MssmSoftsusy::calcDRTrilinears
void calcDRTrilinears(drBarPars &eg, double vev, double beta)
sets the Yukawas and Trilinears
Definition: softsusy.cpp:7768

softsusy::MssmSoftsusy::addStopStop
DoubleMatrix addStopStop(double p, double mt, DoubleMatrix &stop)
Definition: softsusy.cpp:3981

softsusy::MssmSoftsusy::doTadpoles
void doTadpoles(double mt, double sinthDRbar)
Does the full 2-loop calculation of both tadpoles and sets them.
Definition: softsusy.cpp:283

softsusy::MssmSoftsusy::addSnuGaugino
void addSnuGaugino(double p, int family, double &chargino, double &neutralino)
Definition: softsusy.cpp:3863

softsusy::MssmSoftsusy::setTwoLoopMbMtauThresholds
void setTwoLoopMbMtauThresholds(bool sw)
Definition: softsusy.cpp:10904

softsusy::MssmSoftsusy::displayMssmSoft
const MssmSoftsusy & displayMssmSoft() const
Displays whole object as a const.
Definition: softsusy.cpp:158

softsusy::MssmSoftsusy::addSnuTauGaugino
void addSnuTauGaugino(double p, double &chargino, double &neutralino)
Definition: softsusy.cpp:3605

softsusy::MssmSoftsusy::flagProblemThrown
void flagProblemThrown(bool a)
Flags a numerical exception eg number too big/small.
Definition: softsusy.h:242

softsusy::MssmSoftsusy::pis1s1Neutralinos
Complex pis1s1Neutralinos(double p, double q) const
Definition: softsusy.cpp:11534

softsusy::MssmSoftsusy::addSdownNeutralino
DoubleMatrix addSdownNeutralino(double p1, double p2, int family, DoubleMatrix &neutralino)
Definition: softsusy.cpp:5432

softsusy::MssmSoftsusy::addSnuTauEweak
double addSnuTauEweak(double p, double &electroweak)
Definition: softsusy.cpp:3522

softsusy::MssmSoftsusy::pis1s1Fermions
Complex pis1s1Fermions(double p, double q) const
Definition: softsusy.cpp:11237

softsusy::MssmSoftsusy::addStopSbottom
DoubleMatrix addStopSbottom(double p, double mt, DoubleMatrix &sbottom)
Definition: softsusy.cpp:3999

softsusy::MssmSoftsusy::iterateMu
void iterateMu(double &munew, int sgnMu, double mt, int maxTries, double pizztMS, double sinthDRbar, double tol, int &err)
Definition: softsusy.cpp:1267

softsusy::MssmSoftsusy::twoLpMt
double twoLpMt() const

softsusy::MssmSoftsusy::addStauSfermion
void addStauSfermion(double p, double mtau, DoubleMatrix &stop, DoubleMatrix &sbottom)
Definition: softsusy.cpp:4814

softsusy::MssmSoftsusy::doCalcTad2GaugeBosons
double doCalcTad2GaugeBosons(double q, double costhDRbar2, double g, double tanb) const
Definition: softsusy.cpp:656

softsusy::MssmSoftsusy::displayPhys
const sPhysical & displayPhys() const
Displays physical parameters only.
Definition: softsusy.cpp:162

softsusy::MssmSoftsusy::addSdownEweak
DoubleMatrix addSdownEweak(double p1, double p2, int family, DoubleMatrix &electroweak)
Definition: softsusy.cpp:5328

softsusy::MssmSoftsusy::displayFracDiff
double displayFracDiff() const
Displays iteration accuracy attained.
Definition: softsusy.h:173

softsusy::MssmSoftsusy::maxMass
double maxMass() const
Returns the mass of the heaviest SUSY particle, excluding gravitino.
Definition: softsusy.cpp:7216

softsusy::MssmSoftsusy::thet
double thet(double a, double b, double c=0.0)
log(max(a^2, b^2, c^2) / Q^2)
Definition: softsusy.cpp:2123

softsusy::MssmSoftsusy::rewsbM3sqTree
int rewsbM3sqTree(double, double &) const
returns 1 if mu < 1.0e-9
Definition: softsusy.cpp:242

softsusy::MssmSoftsusy::setAllTwoLoopThresholds
void setAllTwoLoopThresholds(bool sw)
Switch all two-loop thresholds to gauge-Yukawa couplings.
Definition: softsusy.cpp:10918

softsusy::MssmSoftsusy::doCalcTadpole1oneLoop
virtual double doCalcTadpole1oneLoop(double mt, double sinthDRbar)
Does the calculation of one-loop pieces of .
Definition: softsusy.cpp:673

softsusy::MssmSoftsusy::physical
virtual void physical(int accuracy)
Definition: softsusy.cpp:6850

softsusy::MssmSoftsusy::calcRunMbHiggs
virtual double calcRunMbHiggs() const
Obtains (1 / mb) times 1-loop Higgs corrections.
Definition: softsusy.cpp:3041

softsusy::MssmSoftsusy::fixedPointIteration
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

softsusy::MssmSoftsusy::predMzsq
double predMzsq(double &tanb, double muOld=-6.66e66, double eps=0.)
Definition: softsusy.cpp:772

softsusy::MssmSoftsusy::printShort
string printShort() const
Prints a list of important sparticle/Higgs masses to standard output.
Definition: softsusy.cpp:1540

softsusy::MssmSoftsusy::rewsbMu
virtual int rewsbMu(int sgnMu, double &mu) const
Definition: softsusy.cpp:196

softsusy::MssmSoftsusy::sfermionsSLHA
virtual void sfermionsSLHA(ostream &out)
sfermions part of mass block of SLHA
Definition: softsusy.cpp:10095

softsusy::MssmSoftsusy::addSbotEweak
DoubleMatrix addSbotEweak(double p, DoubleMatrix &electroweak)
Definition: softsusy.cpp:5786

softsusy::MssmSoftsusy::neutralinoCharginoMSLHA
virtual void neutralinoCharginoMSLHA(ostream &out)
neutralino and charigno part of mass block of SLHA
Definition: softsusy.cpp:10060

softsusy::MssmSoftsusy::displayMwRun
double displayMwRun() const
Returns DRbar MW, must be calculated by calcDrBarPars first.
Definition: softsusy.cpp:721

softsusy::MssmSoftsusy::treeUpSquark
virtual void treeUpSquark(DoubleMatrix &mass, double mtrun, double pizztMS, double sinthDRbarMS, int family)
Definition: softsusy.cpp:3360

softsusy::MssmSoftsusy::piZZTsfermions
virtual double piZZTsfermions(double p, double Q) const
Calculates sfermion contribution to the transverse part of Z self-energy:
Definition: softsusy.cpp:8155

softsusy::MssmSoftsusy::addSupChargino
DoubleMatrix addSupChargino(double p1, double p2, int family, DoubleMatrix &chargino)
Definition: softsusy.cpp:6379

softsusy::MssmSoftsusy::pis1s2Charginos
Complex pis1s2Charginos(double p, double q) const
Definition: softsusy.cpp:11662

softsusy::MssmSoftsusy::addSlepHiggs
DoubleMatrix addSlepHiggs(double p1, double p2, int family, DoubleMatrix &higgs)
Definition: softsusy.cpp:4518

softsusy::MssmSoftsusy::addChaLoopGauge
virtual void addChaLoopGauge(double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS, DoubleMatrix b1pCha, DoubleMatrix b0pCha, DoubleMatrix b1pNeut, DoubleMatrix b0pNeut) const
Definition: softsusy.cpp:1840

softsusy::MssmSoftsusy::operator=
const MssmSoftsusy & operator=(const MssmSoftsusy &s)
Set all data in the object equal to another.
Definition: softsusy.cpp:34

softsusy::MssmSoftsusy::guessAtSusyMt
virtual MssmSusyRGE guessAtSusyMt(double tanb, const QedQcd &oneset)
Definition: softsusy.cpp:7265

softsusy::MssmSoftsusy::piHpHmGauge
double piHpHmGauge(double p, double q) const
LCT: Returns pure gauge contributions to charged Higgs self-energy.
Definition: softsusy.cpp:8683

softsusy::MssmSoftsusy::ewsbConditions
void ewsbConditions(DoubleVector &values) const
Definition: softsusy.cpp:890

softsusy::MssmSoftsusy::assignHiggs
void assignHiggs(DoubleVector &higgsm, DoubleVector &higgsc, DoubleVector &dnu, DoubleVector &dnd, DoubleVector &cn, double beta) const
Definition: softsusy.cpp:4482

softsusy::MssmSoftsusy::setDrBarPars
void setDrBarPars(const drBarPars &s)
Sets tree-level DRbar parameters.
Definition: softsusy.h:262

softsusy::MssmSoftsusy::setMinpot
void setMinpot(double)
Sets potential value at minimum of Higgs potential.
Definition: softsusy.cpp:188

softsusy::MssmSoftsusy::flagTachyon
void flagTachyon(tachyonType a)
Flags a negative-mass squared scalar (really a CCB problem)
Definition: softsusy.h:207

softsusy::MssmSoftsusy::addStopCorrection
virtual void addStopCorrection(double p, DoubleMatrix &mass, double mt)
As in BPMZ appendix.
Definition: softsusy.cpp:4431

softsusy::MssmSoftsusy::addSupCorrection
virtual void addSupCorrection(DoubleMatrix &mass, int family)
As in BPMZ appendix, INCLUDING weak boson loops.
Definition: softsusy.cpp:6432

softsusy::MssmSoftsusy::higgsMSLHA
virtual void higgsMSLHA(ostream &out)
higgs part of mass block of SLHA
Definition: softsusy.cpp:10053

softsusy::MssmSoftsusy::headerSLHA
void headerSLHA(ostream &out)
Definition: softsusy.cpp:9980

softsusy::MssmSoftsusy::displayTadpole2Ms
double displayTadpole2Ms() const
displays t_2/v_2 tadpole
Definition: softsusy.cpp:176

softsusy::MssmSoftsusy::displaySetTbAtMX
bool displaySetTbAtMX() const
Is tan beta set at the user defined SUSY breaking scale?
Definition: softsusy.h:191

softsusy::MssmSoftsusy::piHpHmGauginos
double piHpHmGauginos(double p, double q) const
Definition: softsusy.cpp:8766

softsusy::MssmSoftsusy::calcRunMtStopGluino
virtual double calcRunMtStopGluino() const
Obtains (16 Pi^2) / mt times 1-loop stop-Gluino corrections.
Definition: softsusy.cpp:2662

softsusy::MssmSoftsusy::hmixSLHA
virtual void hmixSLHA(ostream &out)
higgs part of mixing blocks
Definition: softsusy.cpp:10253

softsusy::MssmSoftsusy::useAlternativeMt
void useAlternativeMt()
Use alernative mt(MSSM,DR-bar) calculation.
Definition: softsusy.h:268

softsusy::MssmSoftsusy::calcRunMtNeutralinos
virtual double calcRunMtNeutralinos() const
Obtains (16 Pi^2) / mt times 1-loop Neutralino corrections.
Definition: softsusy.cpp:2740

softsusy::MssmSoftsusy::assignHiggsSfermions
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

softsusy::MssmSoftsusy::setMw
void setMw(double)
sets pole MW prediction
Definition: softsusy.cpp:190

softsusy::MssmSoftsusy::doCalcTad1Neutralinos
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

softsusy::MssmSoftsusy::modselSLHA
void modselSLHA(ostream &out, const char model[])
This does the MODSEL block of SLHA.
Definition: softsusy.cpp:10184

softsusy::MssmSoftsusy::setFracDiff
void setFracDiff(double fD)
sets fracDiff, needed for access by NmssmSoftsusy methods
Definition: softsusy.h:304

softsusy::MssmSoftsusy::calcRunMtauNeutralinos
virtual double calcRunMtauNeutralinos(double mTauSMMZ) const
Definition: softsusy.cpp:3226

softsusy::MssmSoftsusy::pis2s2Charginos
Complex pis2s2Charginos(double p, double q) const
Definition: softsusy.cpp:11693

softsusy::MssmSoftsusy::doCalcTad2Sfermions
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

softsusy::MssmSusyRGE
Contains all supersymmetric RPC-MSSM parameters.
Definition: susy.h:218

softsusy::MssmSusy
Contains all supersymmetric RPC-MSSM parameters.
Definition: susy.h:71

softsusy::MssmSusy::beta
MssmSusy beta(sBrevity &) const
Definition: susy.cpp:300

softsusy::MssmSusy::displayMssmSusy
const MssmSusy & displayMssmSusy() const
Returns whole object as a const.
Definition: susy.h:267

softsusy::QedQcd
Quark and lepton masses and gauge couplings in QEDxQCD effective theory.
Definition: lowe.h:55

softsusy::RGE
Describes a set of parameters and RGEs in general.
Definition: rge.h:49

softsusy::RGE::displayMu
double displayMu() const
Return renomalisation scale.
Definition: rge.h:67

m1
double m1
decay global variable declarations
Definition: decays.cpp:14

def.h
switches (options) and parameters such as default fermion masses, required accuracy etc

higher_order.h
Interface to higher loop corrections.

lowe.h
QedQcd object contains Standard Model quark and lepton masses. It integrates them using 3 loop qcd x ...

mssmUtils.h
A few handy routines for the MSSM: SUSY breaking conditions etc.

mssmjacobian.h
Provides routines for calculating Jacobian fine-tuning.

softsusy
global variable declaration
Definition: def.cpp:13

softsusy::mass
mass
used to give order of quark masses stored
Definition: lowe.h:46

softsusy::boundaryCondition
DoubleVector boundaryCondition(double m0, double m12, double a0, const RpvSusyPars &r)
Definition: rpvsusypars.cpp:373

softsusy::rho2
double rho2(double r)
Two-loop Standard Model corrections to rho parameter.
Definition: softsusy.cpp:9239

softsusy::EPoleUncertainties
EPoleUncertainties
uncertainty flags
Definition: softsusy.h:41

softsusy::fEff
double fEff(double x)
function used for calculating sin theta_eff
Definition: softsusy.cpp:9259

softsusy::lnLHiggs
double lnLHiggs(double mh)
returns the smeared log likelihood coming from LEP2 Higgs mass bounds
Definition: softsusy.cpp:10711

softsusy::trilinears
trilinears
SUSY breaking trilinear parameters.
Definition: softpars.h:25

softsusy::printShortInitialise
void printShortInitialise()
Prints out header line for print-short output.
Definition: softsusy.cpp:1530

softsusy::gEff
double gEff(double x)
function used for calculating sin theta_eff
Definition: softsusy.cpp:9266

softsusy::MZ
double MZ
global pole mass of MZ in GeV - MZCENT is defined in def.h
Definition: def.cpp:31

softsusy::sw2
double sw2
Some parameters used in the computation.
Definition: def.cpp:51

softsusy::PRINTOUT
int PRINTOUT
no verbose debug output
Definition: def.cpp:15

softsusy::mhIntegrand
DoubleVector mhIntegrand(double mh, const DoubleVector &)
Definition: softsusy.cpp:10702

softsusy::minimufb3
double minimufb3(double lnH2)
Definition: softsusy.cpp:6960

softsusy::operator>>
istream & operator>>(istream &left, flavourPhysical &s)
Formatted input of physical parameters.
Definition: flavoursoft.cpp:1220

softsusy::lep2Likelihood
double lep2Likelihood(double mh)
Fit to LEP2 Standard Model results.
Definition: softsusy.cpp:10691

numerics.h
numerical routines - differential equation solver, differentiator and function minimiser for instance

physpars.h
Flags for potential problems in sProblem structure, and structure for containing physical MSSM parame...

softpars.h
Soft SUSY breaking parameters.

softsusy::drBarPars
DRbar values of masses and mixings in MSSM.
Definition: physpars.h:120

softsusy::sBrevity
Contains data needed in beta function calculation to make it faster.
Definition: susy.h:31

softsusy::sPhysical
Masses of the physical particles.
Definition: physpars.h:32

softsusy::sProblem
Various boolean values to flag any problems with the parameter point.
Definition: physpars.h:75

softsusy::sProblem::muSqWrongSign
bool muSqWrongSign
mu^2 came out with wrong sign; no REWSB
Definition: physpars.h:86

softsusy::sProblem::tachyonWarning
tachyonType tachyonWarning
Definition: physpars.h:85

softsusy::sProblem::irqfp
bool irqfp
Infra-red quasi fixed point breached.
Definition: physpars.h:78

softsusy::sProblem::notGlobalMin
bool notGlobalMin
Not in global minimum of Higgs potential.
Definition: physpars.h:89

softsusy::sProblem::tachyon
tachyonType tachyon
Definition: physpars.h:81

softsusy::sProblem::higgsUfb
bool higgsUfb
Higgs potential inconsistent with a good minimum.
Definition: physpars.h:88

softsusy::sProblem::inaccurateHiggsMass
bool inaccurateHiggsMass
Higgs mass is potentially inaccurate and cant be trusted.
Definition: physpars.h:93

softsusy::sProblem::noConvergence
bool noConvergence
Iteration did not converge: not always serious.
Definition: physpars.h:80

softsusy::sProblem::noRhoConvergence
bool noRhoConvergence
Couldn't calculate electroweak rho parameter.
Definition: physpars.h:79

softsusy::sProblem::noMuConvergence
bool noMuConvergence
Definition: physpars.h:91

softsusy::sProblem::m3sq
bool m3sq
m3sq came out with wrong sign; no REWSB
Definition: physpars.h:87

softsusy::sProblem::badConvergence
bool badConvergence
Nowhere near a decent solution.
Definition: physpars.h:77

softsusy::sProblem::nonperturbative
bool nonperturbative
Running went non-perturbative.
Definition: physpars.h:90

accuracy
const double accuracy
Approximate accuracy with which 3 body decays are calculated.
Definition: threeBodyDecays.h:38

twoloophiggs.h
Contains two-loop routines from Slavich et al.

close
bool close(double m1, double m2, double tol)
Definition: utils.cpp:68

utils.h
A few handy bits and pieces - little mathematical functions and the like.