softpars.h

Go to the documentation of this file.

 
#ifndef SOFTPARS_H
#define SOFTPARS_H
 
#include <cmath>
#include "susy.h"
#include "linalg.h"
#include "utils.h"
#include "numerics.h"
 
namespace softsusy {
  
  typedef enum {mQl=1, mUr, mDr, mLl, mEr} softMasses;
  typedef enum {UA=1, DA, EA} trilinears;
  
  const static int numSoftParsMssm = 78 + numSusyPars;
  
  class MssmSoftPars {
  private:
    DoubleVector mGaugino; 
    DoubleMatrix ua, da, ea; 
    DoubleMatrix mQLsq, mURsq, mDRsq, mLLsq, mSEsq; 
    double m3sq, mH1sq, mH2sq; 
    double m32;         
  public:
    MssmSoftPars();
    //  MssmSoftPars(const MssmSusy &);
    MssmSoftPars(const MssmSoftPars &);
    const MssmSoftPars & operator=(const MssmSoftPars & s);
    MssmSoftPars(const DoubleVector & mG, const
                 DoubleMatrix & aU, const DoubleMatrix & aD, const DoubleMatrix
                 & aE, const DoubleMatrix & mQl, const DoubleMatrix & mUr, const
                 DoubleMatrix & mDr, const DoubleMatrix & mLl, const
                 DoubleMatrix & mEr, double m3sq, double mH1sq, double mH2sq,
                 double mGravitino);
    
    const MssmSoftPars & displayMssmSoftPars() const { return *this; }
    const DoubleMatrix & displayTrilinear(trilinears) const;
    double displayTrilinear(trilinears, int i, int j) const;
    const DoubleMatrix & displaySoftMassSquared(softMasses) const;
    double displaySoftMassSquared(softMasses, int i, int j) const;
    
    double displayGravitino() const; 
    double displayM3Squared() const;     
    double displayMh1Squared() const;    
    double displayMh2Squared() const;    
    DoubleVector displayGaugino() const; 
    double displayGaugino(int i) const;  
    void display(DoubleVector & v, int & k) const;
    const DoubleVector display2() const;
    
    void setM32(double);
    void setSoftPars(MssmSoftPars const &);
    void setSoftMassElement(softMasses, int, int, double);
    void setSoftMassMatrix(softMasses, const DoubleMatrix &);
    void setTrilinearMatrix(trilinears, const DoubleMatrix &);
    void setTrilinearElement(trilinears k, int i, int j, double a);
    void setGauginoMass(int, double);
    void setAllGauginos(const DoubleVector &);
    void setM3Squared(double);  
    void setMh1Squared(double); 
    void setMh2Squared(double); 
    void set(const DoubleVector & y);
 
    DoubleVector beta(const MssmSusy &) const; 
    MssmSoftPars beta2(const MssmSusy & xx) const;
    MssmSoftPars beta2(const MssmSusy & xx, sBrevity& a) const;
    void anomalousDeriv(const MssmSusy & xx, DoubleMatrix & gEE, 
                        DoubleMatrix & gLL,
                        DoubleMatrix & gQQ, DoubleMatrix & gUU,
                        DoubleMatrix & gDD, 
                        double & gH1H1, double & gH2H2) const;
    void yTildes(const MssmSusy & xx, 
                 DoubleMatrix & yu, DoubleMatrix & yd, DoubleMatrix &ye) const;
    void u1R_PQflip(MssmSusy & s);
    void universal(MssmSusy & s, double m0,  double m12,  double a0,  double mu,
                   double m3sq);
    void addAmsb(const MssmSusy &, double m32);
    void universalTrilinears(const MssmSusy & a, double a0);
    void minimalGmsb(const MssmSusy & xx, int n5, double LAMBDA, 
                     double mMess, double cgrav);  
    void standardSugra(const MssmSusy & xx, double m0,  double m12, double a0);
    void universalScalars(double m0);
    void universalGauginos(double m12);
    
    void inputSoftParsOnly();
  };
  
  ostream & operator <<(ostream &left, const MssmSoftPars &s);
  istream & operator >>(istream &left, MssmSoftPars &s);
  
  
} // namespace softsusy
 
#endif