SUSYSpectrum.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef SUSYSPECTRUM_H
#define SUSYSPECTRUM_H
 
#include <gslpp.h>
 
class SUSY;
 
class SUSYSpectrum {
public:
 
    SUSYSpectrum(const SUSY& SUSY_in);
 
    bool CalcHiggs(double mh[4], gslpp::complex& saeff_i);
 
    bool CalcChargino(gslpp::matrix<gslpp::complex>& U_i, gslpp::matrix<gslpp::complex>& V_i, gslpp::vector<double>& mch_i);
 
    bool CalcNeutralino(gslpp::matrix<gslpp::complex>& N_i, gslpp::vector<double>& mneu_i);
 
    bool CalcSup(gslpp::matrix<gslpp::complex>& Ru_i, gslpp::vector<double>& m_su2_i);
 
    bool CalcSdown(gslpp::matrix<gslpp::complex>& Rd_i, gslpp::vector<double>& m_sd2_i);
 
    bool CalcSneutrino(gslpp::matrix<gslpp::complex>& Rn_i, gslpp::vector<double>& m_sn2_i);
 
    bool CalcSelectron(gslpp::matrix<gslpp::complex>& Rl_i, gslpp::vector<double>& m_se2_i);
 
    bool CalcSpectrum();
    
    void SortSfermionMasses(gslpp::vector<double>& m_sf2, gslpp::matrix<gslpp::complex>& Rf) const;
 
 
    double getMHl() const
    {
        return mh[0];
    }
 
    double getMHh() const
    {
        return mh[1];
    }
 
    double getMHa() const
    {
        return mh[2];
    }
 
    double getMHp() const
    {
        return mh[3];
    }
 
 
    gslpp::matrix<gslpp::complex> getMchargino() const
    {
        return Mchargino;
    }
 
    gslpp::vector<double> getMch() const
    {
        return mch;
    }
 
    gslpp::matrix<gslpp::complex> getU() const
    {
        return U;
    }
 
    gslpp::matrix<gslpp::complex> getV() const
    {
        return V;
    }
 
 
    gslpp::matrix<gslpp::complex> getMneutralino() const
    {
        return Mneutralino;
    }
 
    gslpp::vector<double> getMneu() const
    {
        return mneu;
    }
 
    gslpp::matrix<gslpp::complex> getN() const
    {
        return N;
    }
 
 
    gslpp::matrix<gslpp::complex> getMsup2() const
    {
        return Msup2;
    }
 
    gslpp::matrix<gslpp::complex> getMsdown2() const
    {
        return Msdown2;
    }
    
    gslpp::vector<double> getMsu2() const
    {
        return m_su2;
    }
 
    gslpp::vector<double> getMsd2() const
    {
        return m_sd2;
    }
 
    gslpp::matrix<gslpp::complex> getRu() const
    {
        return Ru;
    }
 
    gslpp::matrix<gslpp::complex> getRd() const
    {
        return Rd;
    }
 
    
 
    gslpp::matrix<gslpp::complex> getMsneutrino2() const
    {
        return Msneutrino2;
    }
 
    gslpp::matrix<gslpp::complex> getMselectron2() const
    {
        return Mselectron2;
    }
 
    gslpp::vector<double> getMsn2() const
    {
        return m_sn2;
    }
 
    gslpp::vector<double> getMse2() const
    {
        return m_se2;
    }
 
    gslpp::matrix<gslpp::complex> getRn() const
    {
        return Rn;
    }
 
    gslpp::matrix<gslpp::complex> getRl() const
    {
        return Rl;
    }
 
 
private:
    const SUSY& mySUSY;
 
    double mh[4];
 
    gslpp::complex saeff;
 
    gslpp::matrix<gslpp::complex> Mchargino, Mneutralino;
    
    gslpp::matrix<gslpp::complex> U, V, N;
 
    gslpp::matrix<gslpp::complex> Msup2, Msdown2, Msneutrino2, Mselectron2;
 
    gslpp::vector<double> mch, mneu, m_su2, m_sd2, m_sn2, m_se2;
 
    gslpp::matrix<gslpp::complex> Ru, Rd, Rn, Rl;
 
};
 
#endif  /* SUSYSPECTRUM_H */