nmssmsoftpars.h

Go to the documentation of this file.

 
#ifndef NMSSMSOFTPARS_H
#define NMSSMSOFTPARS_H
 
#include <cmath>
#include "linalg.h"
#include "utils.h"
#include "numerics.h"
#include "nmssmsusy.h"
#include "softpars.h"
 
namespace softsusy {
  
  const static int numSoftParsNmssm = 83 + numNMssmPars;
 
  class SoftParsNmssm; class NmssmSoftsusy;
  void addBetaSoftParsNmssm
  (nmsBrevity & a, const NmssmSusyPars & n, const MssmSusy & m, 
   const MssmSoftPars & msoft, const SoftParsNmssm & nsoft, DoubleVector & dmG, 
   double & dmH1sq, double & dmH2sq, double & dm3sq, double & dmSsq, 
   double & dmSpsq, double & dxiSsq, DoubleMatrix & dmq, DoubleMatrix & dmu, 
   DoubleMatrix & dmd, DoubleMatrix & dme, DoubleMatrix dml, 
   DoubleMatrix & dhu, DoubleMatrix & dhd, DoubleMatrix & dhe, double & dhlam, 
   double & dhkap);
  class SoftParsNmssm {
  private:
    double alambda, akappa;
    double mSsq, mSpsq, xiS;
    bool Z3; 
  public:
    SoftParsNmssm()
      : alambda(0.0), akappa(0.0), mSsq(0.0), mSpsq(0.0), xiS(0.0), Z3(true) {};
    SoftParsNmssm(const SoftParsNmssm & s)
    : alambda(s.displayTrialambda()), akappa(s.displayTriakappa()), 
      mSsq(s.displayMsSquared()), mSpsq(s.displayMspSquared()), 
      xiS(s.displayXiS()), Z3(s.Z3) {};
    SoftParsNmssm(double al, double ak, double ms, double msp, double x)
      : alambda(al), akappa(ak), mSsq(ms), mSpsq(msp), xiS(x), Z3(true) {}
    const SoftParsNmssm & operator=(const SoftParsNmssm & s);
 
    
    const SoftParsNmssm & displaySoftParsNmssm() const { return *this; }
    //PA: Return trilinear soft mass $a_\lambda S H_u H_d$
    double displayTrialambda() const { return alambda; };
    //PA: Return trillinear soft mass $a_\kappa S S S$
    double displayTriakappa() const { return akappa; };
    double displayMsSquared() const { return mSsq; }; 
    double displayMspSquared() const { return mSpsq; }; 
    double displayXiS() const { return xiS; }; 
    double displaySoftAlambda(double lam) const;
    double displaySoftAkappa(double ka) const;
    //    const DoubleVector display(const NmssmSusy & n, const MssmSoftPars & s)       const;
    void display(DoubleVector & y, int & k) const;
    bool displayZ3() const { return Z3; }
 
    void setZ3(bool z) { Z3 = z; }
    void setTrialambda(double al) { alambda = al; };
    void setTriakappa(double ak) { akappa = ak; };
    void setMsSquared(double f) { mSsq = f; }; 
    void setMspSquared(double f) { mSpsq = f; }; 
    void setXiS(double f) { xiS = f; }; 
    void setSoftParsNmssm(SoftParsNmssm const & s) { *this = s; }
    //    void set(const DoubleVector &);    
    void set(const DoubleVector &, int & k);    
 
    void anomalousDeriv(DoubleMatrix & gEE, DoubleMatrix & gLL,
                        DoubleMatrix & gQQ, DoubleMatrix & gUU,
                        DoubleMatrix & gDD,
                        double & gH1H1, double & gH2H2, double & gSS) const;
    
    void addAmsb(double maux, const NmssmSusy & Nms, MssmSoftPars & m);
    void standardSugra(double m0, double m12, double a0, 
                       const NmssmSusyPars & n, const MssmSusy & m,
                       const MssmSoftPars & msoft);
    void standardsemiSugra(double m0,  double m12, double a0, double Al, 
                           double Ak, const NmssmSusy & n, 
                           const MssmSoftPars & m, double mS=0.0);
    void universalScalars(double m0);
    void semiuniversalScalars(double m0);
    void universalTrilinears(double a0, const NmssmSusy & n);
    void semiuniversalTrilinears(double a0, double al, double ak, 
                                 const NmssmSusy & n);
    
    void inputSoftParsOnly();
  };
  
  ostream & operator <<(ostream &left, const SoftParsNmssm &s);
} // namespace softsusy
 
#endif