mssmjacobian.h

Go to the documentation of this file.

 
#ifndef MSSMJACOBIAN_H
#define MSSMJACOBIAN_H
 
#include "linalg.h"
 
#include <utility>
#include <limits>
 
namespace softsusy {
 
  class MssmSoftsusy;
 
  class MssmJacobian {
  public:
 
    explicit MssmJacobian(bool doTop = false);
    ~MssmJacobian();
 
    bool displayIncludeTopFlag() const { return includeTop; }
 
    void setIncludeTopFlag(bool flag) { includeTop = flag; }
 
    bool displayUseRunningMassesFlag() const { return useRunningMasses; }
 
    void setUseRunningMassesFlag(bool flag) { useRunningMasses = flag; }
 
    bool displayUseSugraTrilinearsFlag() const { return useSugraTrilinears; }
 
    void setUseSugraTrilinearsFlag(bool flag) { useSugraTrilinears = flag; }
 
    void resetError() { hasError = false; }
 
    bool displayError() const { return hasError; }
 
    DoubleMatrix displayInverseRGFlowJacobian() const { return invJacRGFlow; }
 
    DoubleMatrix displayInverseRGFlowJacobianErrors() const {
      return invJacRGFlowErrors;
    }
 
    DoubleMatrix displayInverseEWSBJacobian() const { return invJacEWSB; }
 
    DoubleMatrix displayInverseEWSBJacobianErrors() const {
      return invJacEWSBErrors;
    }
 
    double calcFTInverseJacobian(const MssmSoftsusy& m);
 
    double calcFTInverseJacobian(const MssmSoftsusy& m, double mx);
 
    double calcDeltaJ(const MssmSoftsusy& m);
 
    double calcDeltaJ(const MssmSoftsusy& m, double mx);
 
  private:
    enum Parameters { Mzsq, Tanb, Mtsq, SMu, M3Sq, Mh1Sq, Mh2Sq, Yt };
 
    DoubleMatrix invJacRGFlow;
    DoubleMatrix invJacEWSB;
 
    DoubleMatrix invJacRGFlowErrors;
    DoubleMatrix invJacEWSBErrors;
 
    bool includeTop;
    bool useRunningMasses;
    bool useSugraTrilinears;
    bool hasError;
 
    struct EWSBPars {
      MssmSoftsusy* model;
      Parameters independent;
      Parameters dependent;
      DoubleVector outputs;
      bool useRunningMasses;
      bool useSugraTrilinears;
 
      EWSBPars()
        : model(0), independent(Mzsq), dependent(SMu),
          outputs(2), useRunningMasses(false),
          useSugraTrilinears(false)
        {}
    };
 
    struct RGFlowPars {
      MssmSoftsusy* model;
      Parameters independent;
      Parameters dependent;
      double toScale;
      bool useSugraTrilinears;
    };
 
    static double calcMz(MssmSoftsusy& model, bool getRunningMass = false);
    static double calcMt(MssmSoftsusy& model, bool getRunningMass = false);
 
    static double calcRunningParameter(double x, void* parameters);
    std::pair<double,double> calcRGDerivative(MssmSoftsusy& model,
                                              Parameters dep, Parameters indep,
                                              double toScale);
    double calcRGFlowJacobian(MssmSoftsusy& model, double startScale,
                              double endScale);
    static double calcEWSBOutput(double x, void* parameters);
    std::pair<double,double> calcEWSBOutputDerivative(MssmSoftsusy& model,
                                                      Parameters dep,
                                                      Parameters indep);
    double calcInverseEWSBJacobian(MssmSoftsusy& model);
  };
 
} 
 
#endif