SUSYSpectrum.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * 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();


    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;
    }
    
    void SortSfermionMasses(gslpp::vector<double>& m_sf2, gslpp::matrix<gslpp::complex>& Rf) const;


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 */