A class for calculating the Higgs and sparticle spectra at tree level. More...

#include <SUSYSpectrum.h>

Detailed Description

A class for calculating the Higgs and sparticle spectra at tree level.

Author: HEPfit Collaboration

Copyright: GNU General Public License

The SUSYSpectrum class calculates the all the sparticle masses and their mixing matrices at tree-level.

Definition at line 23 of file SUSYSpectrum.h.

Public Member Functions
bool	CalcChargino (gslpp::matrix< gslpp::complex > &U_i, gslpp::matrix< gslpp::complex > &V_i, gslpp::vector< double > &mch_i)
	Computes the chargino spectrum at tree level. More...

bool	CalcHiggs (double mh[4], gslpp::complex &saeff_i)
	Computes the Higgs spectrum at tree level. More...

bool	CalcNeutralino (gslpp::matrix< gslpp::complex > &N_i, gslpp::vector< double > &mneu_i)
	Computes the neutralino spectrum at tree level. More...

bool	CalcSdown (gslpp::matrix< gslpp::complex > &Rd_i, gslpp::vector< double > &m_sd2_i)
	Computes the down-type squark spectrum at tree level. More...

bool	CalcSelectron (gslpp::matrix< gslpp::complex > &Rl_i, gslpp::vector< double > &m_se2_i)
	Computes the charged-slepton spectrum at tree level. More...

bool	CalcSneutrino (gslpp::matrix< gslpp::complex > &Rn_i, gslpp::vector< double > &m_sn2_i)
	Computes the sneutrino spectrum at tree level. More...

bool	CalcSpectrum ()
	Computes the SUSY spectrum without the Higgs part at tree level. More...

bool	CalcSup (gslpp::matrix< gslpp::complex > &Ru_i, gslpp::vector< double > &m_su2_i)
	Computes the up-type squark spectrum at tree level. More...

gslpp::vector< double >	getMch () const
	Gets the Chargino spectrum at tree-level. More...

gslpp::matrix< gslpp::complex >	getMchargino () const
	Gets the Chargino spectrum at tree-level. More...

double	getMHa () const
	Gets the pseudo-scalar Higgs mass. More...

double	getMHh () const
	Gets the heavy Higgs mass. More...

double	getMHl () const
	Gets the light Higgs mass. More...

double	getMHp () const
	Gets the charged Higgs mass. More...

gslpp::vector< double >	getMneu () const
	Gets the Neutralino spectrum at tree-level. More...

gslpp::matrix< gslpp::complex >	getMneutralino () const
	Gets the Neutralino spectrum at tree-level. More...

gslpp::vector< double >	getMsd2 () const
	Gets the Down-squark spectrum at tree-level. More...

gslpp::matrix< gslpp::complex >	getMsdown2 () const
	Gets the Down-squark mass matrix at tree-level. More...

gslpp::vector< double >	getMse2 () const
	Gets the Slepton spectrum at tree-level. More...

gslpp::matrix< gslpp::complex >	getMselectron2 () const
	Gets the Slepton mass matrix at tree-level. More...

gslpp::vector< double >	getMsn2 () const
	Gets the Sneutrino spectrum at tree-level. More...

gslpp::matrix< gslpp::complex >	getMsneutrino2 () const
	Gets the Sneutrino mass matrix at tree-level. More...

gslpp::vector< double >	getMsu2 () const
	Gets the Up-squark spectrum at tree-level. More...

gslpp::matrix< gslpp::complex >	getMsup2 () const
	Gets the Up-squark mass matrix at tree-level. More...

gslpp::matrix< gslpp::complex >	getN () const
	Gets the Neutralino mixing matrix. More...

gslpp::matrix< gslpp::complex >	getRd () const
	Gets the Down-squark mixing matrix. More...

gslpp::matrix< gslpp::complex >	getRl () const
	Gets the Slepton mixing matrix. More...

gslpp::matrix< gslpp::complex >	getRn () const
	Gets the Sneutrino mixing matrix. More...

gslpp::matrix< gslpp::complex >	getRu () const
	Gets the Up-squark mixing matrix. More...

gslpp::matrix< gslpp::complex >	getU () const
	Gets the Chargino mixing matrix U. More...

gslpp::matrix< gslpp::complex >	getV () const
	Gets the Chargino mixing matrix V. More...

void	SortSfermionMasses (gslpp::vector< double > &m_sf2, gslpp::matrix< gslpp::complex > &Rf) const

	SUSYSpectrum (const SUSY &SUSY_in)
	A SUSYSpectrum constructor. More...

Private Attributes
gslpp::vector< double >	m_sd2

gslpp::vector< double >	m_se2

gslpp::vector< double >	m_sn2

gslpp::vector< double >	m_su2

gslpp::vector< double >	mch
	Stores the tree-level Up-squark, Down-squark, Sneutrino, and Slepton mass-squared eigenvalues. More...

gslpp::matrix< gslpp::complex >	Mchargino
	Stores the tree-level Chargino and Neutralino mass matrix. More...

double	mh [4]
	Stores the tree-level Higgs spectrum. More...

gslpp::vector< double >	mneu

gslpp::matrix< gslpp::complex >	Mneutralino

gslpp::matrix< gslpp::complex >	Msdown2

gslpp::matrix< gslpp::complex >	Mselectron2

gslpp::matrix< gslpp::complex >	Msneutrino2

gslpp::matrix< gslpp::complex >	Msup2
	Stores the tree-level Up-squark, Down-squark, Sneutrino, and Slepton mass matrix. More...

const SUSY &	mySUSY

gslpp::matrix< gslpp::complex >	N

gslpp::matrix< gslpp::complex >	Rd

gslpp::matrix< gslpp::complex >	Rl

gslpp::matrix< gslpp::complex >	Rn

gslpp::matrix< gslpp::complex >	Ru
	Stores the tree-level Up-squark, Down-squark, Sneutrino, and Slepton mixing matrices. More...

gslpp::complex	saeff
	Stores the Sine of tree-level CP-even mixing angle. More...

gslpp::matrix< gslpp::complex >	U
	Stores the tree-level Chargino and Neutralino mixing matrices. More...

gslpp::matrix< gslpp::complex >	V

Constructor & Destructor Documentation

◆ SUSYSpectrum()

SUSYSpectrum::SUSYSpectrum ( const SUSY & SUSY_in )

A SUSYSpectrum constructor.

Parameters

[in] SUSY_in An object of SUSY class.

Definition at line 15 of file SUSYSpectrum.cpp.

: mySUSY(SUSY_in), Mchargino(2,2,0.), Mneutralino(4,4,0.),
        U(2,2,0.), V(2,2,0.), N(4,4,0.),
        Msup2(6,0.), Msdown2(6,0.), Msneutrino2(6,0.), Mselectron2(6,0.),
        mch(2,0.), mneu(4,0.), m_su2(6,0.), m_sd2(6,0.), m_sn2(6,0.), m_se2(6,0.),
        Ru(6,6,0.), Rd(6,6,0.), Rn(6,6,0.), Rl(6,6,0.)
{
}

Member Function Documentation

◆ CalcChargino()

bool SUSYSpectrum::CalcChargino	(	gslpp::matrix< gslpp::complex > &	U_i,
		gslpp::matrix< gslpp::complex > &	V_i,
		gslpp::vector< double > &	mch_i
	)

Computes the chargino spectrum at tree level.

SLHA: diag(mChi) = U^* M V^+.

Definition at line 52 of file SUSYSpectrum.cpp.

{
    double Mw = mySUSY.Mw_tree();
 
    Mchargino.assign(0, 0, mySUSY.getM2());
    Mchargino.assign(0, 1, sqrt(2) * Mw * mySUSY.getSinb());
    //
    Mchargino.assign(1, 0, sqrt(2) * Mw * mySUSY.getCosb());
    Mchargino.assign(1, 1, mySUSY.getMuH());
 
    /*
     * M.singularvalue(&U, &V, &S) decompose M into
     * M = U diag(S) V^+.
     */
    gslpp::matrix<gslpp::complex> Utmp(2, 2, 0.), Vtmp(2, 2, 0.);
    Mchargino.singularvalue(Utmp, Vtmp, mch_i);
 
    U_i = Utmp.transpose();
    V_i = Vtmp.hconjugate();
 
    return true;
}

◆ CalcHiggs()

bool SUSYSpectrum::CalcHiggs	(	double	mh[4],
		gslpp::complex &	saeff_i
	)

Computes the Higgs spectrum at tree level.

Definition at line 24 of file SUSYSpectrum.cpp.

{
    double Mw = mySUSY.Mw_tree();
    double Mz = mySUSY.getMz();
    double cos2b = 2.0 * mySUSY.getCosb() * mySUSY.getCosb() - 1.0;
 
    /* charged Higgs mass */
    mh[3] = mySUSY.mHptree;
 
    /* pseudo-scalar Higgs mass */
    mh[2] = sqrt(mh[3] * mh[3] - Mw * Mw);
 
    double temp = mh[2] * mh[2] + Mz * Mz;
    double temp1 = 2.0 * mh[2] * Mz * cos2b;
    double temp2 = sqrt(temp * temp - temp1 * temp1);
 
    /* light Higgs mass */
    mh[0] = sqrt((temp - temp2)/2.0);
    
    /* heavy Higgs mass */
    mh[1] = sqrt((temp + temp2)/2.0);
    
    /* CP-even Higgs mixing angle*/
    saeff_i = sin(atan((mh[2]*mh[2]+Mz*Mz)*sqrt(1-cos2b*cos2b)/((mh[2]*mh[2]-Mz*Mz)*cos2b))/2.0);
 
    return true;
}

◆ CalcNeutralino()

bool SUSYSpectrum::CalcNeutralino	(	gslpp::matrix< gslpp::complex > &	N_i,
		gslpp::vector< double > &	mneu_i
	)

Computes the neutralino spectrum at tree level.

Definition at line 78 of file SUSYSpectrum.cpp.

{
    double Mw = mySUSY.Mw_tree();
    double Mz = mySUSY.getMz();
    double cW2 = Mw*Mw/Mz/Mz;
    double cW = sqrt(cW2);
    double sW = sqrt(1.0 - cW2);
    double sb = mySUSY.getSinb();
    double cb = mySUSY.getCosb();
 
    Mneutralino.assign(0, 0, mySUSY.getM1());
    Mneutralino.assign(0, 2, - Mz * sW * cb);
    Mneutralino.assign(0, 3, Mz * sW * sb);
    //
    Mneutralino.assign(1, 1, mySUSY.getM2());
    Mneutralino.assign(1, 2, Mz * cW * cb);
    Mneutralino.assign(1, 3, - Mz * cW * sb);
    //
    Mneutralino.assign(2, 0, Mneutralino(0,2));
    Mneutralino.assign(2, 1, Mneutralino(1,2));
    Mneutralino.assign(2, 3, - mySUSY.getMuH());
    //
    Mneutralino.assign(3, 0, Mneutralino(0,3));
    Mneutralino.assign(3, 1, Mneutralino(1,3));
    Mneutralino.assign(3, 2, Mneutralino(2,3));
 
    /*
     * M.singularvalue(&U, &V, &S) decompose M into
     * M = U diag(S) V^+.
     */
    gslpp::matrix<gslpp::complex> Ntemp1(4, 4, 0.), Ntemp2(4, 4, 0.);
    Mneutralino.singularvalue(Ntemp1, Ntemp2, mneu_i);
 
    //gslpp::matrix<gslpp::complex> Nleft = Ntemp1.transpose();
    gslpp::matrix<gslpp::complex> Nright = Ntemp2.hconjugate();
    //std::cout << "Nleft = " << Nleft << std::endl;
    //std::cout << "Nright = " << Nright << std::endl;
 
    /* adopt N=Nright as N^* M N^+ = diag, not as N M N^+.
     * As a result, a phase rotation is required for N. */
    gslpp::matrix<gslpp::complex> Mdiag_tmp(4, 4, 0.);
    Mdiag_tmp = Nright.hconjugate().transpose() * Mneutralino * Nright.hconjugate();
    //std::cout << "Mdiag_tmp = " << Mdiag_tmp << std::endl;
    gslpp::vector<gslpp::complex> v1(4, 0.);
    for(int i = 0; i < 4; i++)
        v1.assign(i, gslpp::complex(1., Mdiag_tmp(i,i).arg()/2.0, true));
    Nright = gslpp::matrix<gslpp::complex>(v1) * Nright;
 
    N_i = Nright;
 
    return true;
}

◆ CalcSdown()

bool SUSYSpectrum::CalcSdown	(	gslpp::matrix< gslpp::complex > &	Rd_i,
		gslpp::vector< double > &	m_sd2_i
	)

Computes the down-type squark spectrum at tree level.

Definition at line 172 of file SUSYSpectrum.cpp.

{
    double Mw = mySUSY.Mw_tree();
    double Mz2 = mySUSY.getMz()*mySUSY.getMz();
    double sW2 = 1.0 - Mw*Mw/Mz2;
    double cos2b = 2.0 * mySUSY.getCosb() * mySUSY.getCosb() - 1.0;
    gslpp::matrix<gslpp::complex> Id3 = gslpp::matrix<gslpp::complex>::Id(3);
 
    /* DRbar down-type quark masses at scale Q */
    gslpp::matrix<double> Md(3, 3, 0.);
    Md(0, 0) = mySUSY.Mq_Q(mySUSY.DOWN);
    Md(1, 1) = mySUSY.Mq_Q(mySUSY.STRANGE);
    Md(2, 2) = mySUSY.Mq_Q(mySUSY.BOTTOM);
 
    gslpp::matrix<gslpp::complex> dLL( mySUSY.msQhat2 + Md * Md
                         + cos2b * Mz2 * (- 1.0/2.0 + 1.0/3.0 * sW2) * Id3 );
    gslpp::matrix<gslpp::complex> dLR( mySUSY.v1()/sqrt(2.0) * mySUSY.getTDhat().hconjugate()
                         - mySUSY.getMuH() * Md * mySUSY.getTanb() );
    gslpp::matrix<gslpp::complex> dRR( mySUSY.msDhat2 + Md * Md - cos2b * Mz2 /3.0 * sW2 * Id3 );
    for(int i = 0; i < 3; i++)
        for(int j = 0; j < 3; j++) {
            Msdown2.assign(i,   j,   dLL(i,j));
            Msdown2.assign(i,   j+3, dLR(i,j));
            Msdown2.assign(i+3, j,   dLR(j,i).conjugate());
            Msdown2.assign(i+3, j+3, dRR(i,j));
        }
 
    /*
     * M.singularvalue(&U, &V, &S) decompose M into
     * M = U diag(S) V^+.
     */
    gslpp::matrix<gslpp::complex> RdTmp(6, 6, 0.);
    Msdown2.eigensystem(RdTmp, m_sd2_i);
 
    Rd_i = RdTmp.hconjugate();
 
    return true;
 
}

◆ CalcSelectron()

bool SUSYSpectrum::CalcSelectron	(	gslpp::matrix< gslpp::complex > &	Rl_i,
		gslpp::vector< double > &	m_se2_i
	)

Computes the charged-slepton spectrum at tree level.

Definition at line 245 of file SUSYSpectrum.cpp.

{
    double Mw = mySUSY.Mw_tree();
    double Mz2 = mySUSY.getMz()*mySUSY.getMz();
    double sW2 = 1.0 - Mw*Mw/Mz2;
    double cos2b = 2.0 * mySUSY.getCosb() * mySUSY.getCosb() - 1.0;
    gslpp::matrix<gslpp::complex> Id3 = gslpp::matrix<gslpp::complex>::Id(3);
 
    gslpp::matrix<double> Me(3, 3, 0.);
    Me(0, 0) = mySUSY.Ml_Q(mySUSY.ELECTRON);
    Me(1, 1) = mySUSY.Ml_Q(mySUSY.MU);
    Me(2, 2) = mySUSY.Ml_Q(mySUSY.TAU);
 
    gslpp::matrix<gslpp::complex> eLL( mySUSY.msLhat2 + Me * Me
                         + cos2b * Mz2 * (- 1.0/2.0 + sW2) * Id3 );
    gslpp::matrix<gslpp::complex> eLR( mySUSY.v1()/sqrt(2.0) * mySUSY.getTEhat().hconjugate()
                         - mySUSY.getMuH() * Me * mySUSY.getTanb() );
    gslpp::matrix<gslpp::complex> eRR( mySUSY.msEhat2 + Me * Me - cos2b * Mz2 * sW2 * Id3 );
 
        for(int i = 0; i < 3; i++)
        {
            for(int j = 0; j < 3; j++)
            {
                Mselectron2.assign(i,   j,   eLL(i,j));
                Mselectron2.assign(i,   j+3, eLR(i,j));
                Mselectron2.assign(i+3, j,   eLR(j,i).conjugate());
                Mselectron2.assign(i+3, j+3, eRR(i,j));
            }
        }
 
    /*
     * M.singularvalue(&U, &V, &S) decompose M into
     * M = U diag(S) V^+.
     */
    gslpp::matrix<gslpp::complex> RlTmp(6, 6, 0.);
    Mselectron2.eigensystem(RlTmp, m_se2_i);
    
    Rl_i = RlTmp.hconjugate();
 
    return true;
 
}

◆ CalcSneutrino()

bool SUSYSpectrum::CalcSneutrino	(	gslpp::matrix< gslpp::complex > &	Rn_i,
		gslpp::vector< double > &	m_sn2_i
	)

Computes the sneutrino spectrum at tree level.

Definition at line 212 of file SUSYSpectrum.cpp.

{
    double Mz2 = mySUSY.getMz()*mySUSY.getMz();
    double cos2b = 2.0 * mySUSY.getCosb() * mySUSY.getCosb() - 1.0;
    gslpp::matrix<gslpp::complex> Id3 = gslpp::matrix<gslpp::complex>::Id(3);
 
    gslpp::matrix<gslpp::complex> nLL( mySUSY.msLhat2 + cos2b * Mz2 /2.0 * Id3 );
    gslpp::matrix<gslpp::complex> nLR( mySUSY.v2()/sqrt(2.0) * mySUSY.getTNhat().hconjugate() );
    gslpp::matrix<gslpp::complex> nRR( mySUSY.msNhat2 );
    for(int i = 0; i < 3; i++) {
        for(int j = 0; j < 3; j++) {
            Msneutrino2.assign(i,   j,   nLL(i,j));
            Msneutrino2.assign(i,   j+3, nLR(i,j));
            Msneutrino2.assign(i+3, j,   nLR(j,i).conjugate());
            Msneutrino2.assign(i+3, j+3, nRR(i,j));
        }
    }
    /*
     * M.singularvalue(&U, &V, &S) decompose M into
     * M = U diag(S) V^+.
     */
    gslpp::matrix<gslpp::complex> RnTmp(6, 6, 0.);
    Msneutrino2.eigensystem(RnTmp, m_sn2_i);
 
    m_sn2_i(0)=std::numeric_limits<double>::max();
    m_sn2_i(1)=std::numeric_limits<double>::max();
    m_sn2_i(2)=std::numeric_limits<double>::max();
    Rn_i = RnTmp.hconjugate();
 
    return true;
 
}

◆ CalcSpectrum()

bool SUSYSpectrum::CalcSpectrum ( )

Computes the SUSY spectrum without the Higgs part at tree level.

◆ CalcSup()

bool SUSYSpectrum::CalcSup	(	gslpp::matrix< gslpp::complex > &	Ru_i,
		gslpp::vector< double > &	m_su2_i
	)

Computes the up-type squark spectrum at tree level.

Definition at line 131 of file SUSYSpectrum.cpp.

{
    double Mw = mySUSY.Mw_tree();
    double Mz2 = mySUSY.getMz()*mySUSY.getMz();
    double sW2 = 1.0 - Mw*Mw/Mz2;
    double cos2b = 2.0 * mySUSY.getCosb() * mySUSY.getCosb() - 1.0;
    gslpp::matrix<gslpp::complex> CKM(mySUSY.getVCKM());
    gslpp::matrix<gslpp::complex> Id3 = gslpp::matrix<gslpp::complex>::Id(3);
    
    /* DRbar up-type quark masses at scale Q */
    gslpp::matrix<double> Mu(3, 3, 0.);
    Mu(0, 0) = mySUSY.Mq_Q(mySUSY.UP);
    Mu(1, 1) = mySUSY.Mq_Q(mySUSY.CHARM);
    Mu(2, 2) = mySUSY.Mq_Q(mySUSY.TOP);
 
    gslpp::matrix<gslpp::complex> uLL( CKM * mySUSY.msQhat2 * CKM.hconjugate() + Mu * Mu
                         + cos2b * Mz2 * (1.0/2.0 - 2.0/3.0 * sW2) * Id3 );
    gslpp::matrix<gslpp::complex> uLR( mySUSY.v2()/sqrt(2.0) * mySUSY.getTUhat().hconjugate()
                         - mySUSY.getMuH() * Mu / mySUSY.getTanb() );
    gslpp::matrix<gslpp::complex> uRR( mySUSY.msUhat2 + Mu * Mu + cos2b * Mz2 * 2.0/3.0 * sW2 * Id3 );
    for(int i = 0; i < 3; i++)
        for(int j = 0; j < 3; j++) {
            Msup2.assign(i,   j,   uLL(i,j));
            Msup2.assign(i,   j+3, uLR(i,j));
            Msup2.assign(i+3, j,   uLR(j,i).conjugate());
            Msup2.assign(i+3, j+3, uRR(i,j));
        }
 
    /*
     * M.singularvalue(&U, &V, &S) decompose M into
     * M = U diag(S) V^+.
     */
    gslpp::matrix<gslpp::complex> RuTmp(6, 6, 0.);
    Msup2.eigensystem(RuTmp, m_su2_i);
 
    Ru_i = RuTmp.hconjugate();
 
    return true;
 
}

◆ getMch()

gslpp::vector< double > SUSYSpectrum::getMch ( ) const

inline

Gets the Chargino spectrum at tree-level.

Returns: the Chargino mass eigenvalues

Definition at line 132 of file SUSYSpectrum.h.

    {
        return mch;
    }

◆ getMchargino()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getMchargino ( ) const

inline

Gets the Chargino spectrum at tree-level.

Returns: the Chargino mass eigenvalues

Definition at line 123 of file SUSYSpectrum.h.

    {
        return Mchargino;
    }

◆ getMHa()

double SUSYSpectrum::getMHa ( ) const

inline

Gets the pseudo-scalar Higgs mass.

Returns: The pseudo-scalar Higgs mass.

Definition at line 103 of file SUSYSpectrum.h.

    {
        return mh[2];
    }

◆ getMHh()

double SUSYSpectrum::getMHh ( ) const

inline

Gets the heavy Higgs mass.

Returns: The heavy Higgs mass.

Definition at line 94 of file SUSYSpectrum.h.

    {
        return mh[1];
    }

◆ getMHl()

double SUSYSpectrum::getMHl ( ) const

inline

Gets the light Higgs mass.

Returns: The light Higgs mass.

Definition at line 85 of file SUSYSpectrum.h.

    {
        return mh[0];
    }

◆ getMHp()

double SUSYSpectrum::getMHp ( ) const

inline

Gets the charged Higgs mass.

Returns: The charged Higgs mass.

Definition at line 112 of file SUSYSpectrum.h.

    {
        return mh[3];
    }

◆ getMneu()

gslpp::vector< double > SUSYSpectrum::getMneu ( ) const

inline

Gets the Neutralino spectrum at tree-level.

Returns: the Neutralino mass eigenvalues

Definition at line 170 of file SUSYSpectrum.h.

    {
        return mneu;
    }

◆ getMneutralino()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getMneutralino ( ) const

inline

Gets the Neutralino spectrum at tree-level.

Returns: the Neutralino mass eigenvalues

Definition at line 161 of file SUSYSpectrum.h.

    {
        return Mneutralino;
    }

◆ getMsd2()

gslpp::vector< double > SUSYSpectrum::getMsd2 ( ) const

inline

Gets the Down-squark spectrum at tree-level.

Returns: the Down-squark mass-squared eigenvalues

Definition at line 217 of file SUSYSpectrum.h.

    {
        return m_sd2;
    }

◆ getMsdown2()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getMsdown2 ( ) const

inline

Gets the Down-squark mass matrix at tree-level.

Returns: the Down-squark mass matrix

Definition at line 199 of file SUSYSpectrum.h.

    {
        return Msdown2;
    }

◆ getMse2()

gslpp::vector< double > SUSYSpectrum::getMse2 ( ) const

inline

Gets the Slepton spectrum at tree-level.

Returns: the Slepton mass-squared eigenvalues

Definition at line 274 of file SUSYSpectrum.h.

    {
        return m_se2;
    }

◆ getMselectron2()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getMselectron2 ( ) const

inline

Gets the Slepton mass matrix at tree-level.

Returns: the Slepton mass matrix

Definition at line 256 of file SUSYSpectrum.h.

    {
        return Mselectron2;
    }

◆ getMsn2()

gslpp::vector< double > SUSYSpectrum::getMsn2 ( ) const

inline

Gets the Sneutrino spectrum at tree-level.

Returns: the Sneutrino mass-squared eigenvalues

Definition at line 265 of file SUSYSpectrum.h.

    {
        return m_sn2;
    }

◆ getMsneutrino2()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getMsneutrino2 ( ) const

inline

Gets the Sneutrino mass matrix at tree-level.

Returns: the Sneutrino mass matrix

Definition at line 247 of file SUSYSpectrum.h.

    {
        return Msneutrino2;
    }

◆ getMsu2()

gslpp::vector< double > SUSYSpectrum::getMsu2 ( ) const

inline

Gets the Up-squark spectrum at tree-level.

Returns: the Up-squark mass-squared eigenvalues

Definition at line 208 of file SUSYSpectrum.h.

    {
        return m_su2;
    }

◆ getMsup2()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getMsup2 ( ) const

inline

Gets the Up-squark mass matrix at tree-level.

Returns: the Up-squark mass matrix

Definition at line 190 of file SUSYSpectrum.h.

    {
        return Msup2;
    }

◆ getN()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getN ( ) const

inline

Gets the Neutralino mixing matrix.

Returns: the Neutralino mixing matrix

Definition at line 179 of file SUSYSpectrum.h.

    {
        return N;
    }

◆ getRd()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getRd ( ) const

inline

Gets the Down-squark mixing matrix.

Returns: the Down-squark mixing matrix

Definition at line 235 of file SUSYSpectrum.h.

    {
        return Rd;
    }

◆ getRl()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getRl ( ) const

inline

Gets the Slepton mixing matrix.

Returns: the Slepton mixing matrix

Definition at line 292 of file SUSYSpectrum.h.

    {
        return Rl;
    }

◆ getRn()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getRn ( ) const

inline

Gets the Sneutrino mixing matrix.

Returns: the Sneutrino mixing matrix

Definition at line 283 of file SUSYSpectrum.h.

    {
        return Rn;
    }

◆ getRu()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getRu ( ) const

inline

Gets the Up-squark mixing matrix.

Returns: the Up-squark mixing matrix

Definition at line 226 of file SUSYSpectrum.h.

    {
        return Ru;
    }

◆ getU()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getU ( ) const

inline

Gets the Chargino mixing matrix U.

Returns: the Chargino mixing matrix U

Definition at line 141 of file SUSYSpectrum.h.

    {
        return U;
    }

◆ getV()

gslpp::matrix< gslpp::complex > SUSYSpectrum::getV ( ) const

inline

Gets the Chargino mixing matrix V.

Returns: the Chargino mixing matrix V

Definition at line 150 of file SUSYSpectrum.h.

    {
        return V;
    }

◆ SortSfermionMasses()

void SUSYSpectrum::SortSfermionMasses	(	gslpp::vector< double > &	m_sf2,
		gslpp::matrix< gslpp::complex > &	Rf
	)		const

Sort sfermion masses in increasing order.

Parameters

[in,out]	m_sf2	A vector of sfermion mass squared.
[in,out]	Rf	The corresponding rotation matrix.

Definition at line 288 of file SUSYSpectrum.cpp.

{
    int newIndex[6];
    for (int i = 0; i < 6; i++)
        newIndex[i] = i;
 
    /* sort sfermion masses in increasing order */
    for (int i = 0; i < 5; i++)
        for (int k = i + 1; k < 6; k++)
            if (m_sf2(i) > m_sf2(k)) {
                std::swap(m_sf2(i), m_sf2(k));
                std::swap(newIndex[i], newIndex[k]);
            }
 
    /* sort the corresponding rotation matrix, where the first(second) index
     * denotes mass(gauge) eigenstates. */
    gslpp::matrix<gslpp::complex> myRf(6, 6, 0.);
    for (int i = 0; i < 6; i++)
        for (int k = 0; k < 6; k++)
            myRf.assign(k, i, Rf(newIndex[k], i));
    Rf = myRf;
}

Member Data Documentation

◆ m_sd2

gslpp::vector<double> SUSYSpectrum::m_sd2

private

Definition at line 330 of file SUSYSpectrum.h.

◆ m_se2

gslpp::vector<double> SUSYSpectrum::m_se2

private

Definition at line 330 of file SUSYSpectrum.h.

◆ m_sn2

gslpp::vector<double> SUSYSpectrum::m_sn2

private

Definition at line 330 of file SUSYSpectrum.h.

◆ m_su2

gslpp::vector<double> SUSYSpectrum::m_su2

private

Definition at line 330 of file SUSYSpectrum.h.

◆ mch

gslpp::vector<double> SUSYSpectrum::mch

private

Stores the tree-level Up-squark, Down-squark, Sneutrino, and Slepton mass-squared eigenvalues.

Definition at line 330 of file SUSYSpectrum.h.

◆ Mchargino

gslpp::matrix<gslpp::complex> SUSYSpectrum::Mchargino

private

Stores the tree-level Chargino and Neutralino mass matrix.

Definition at line 315 of file SUSYSpectrum.h.

◆ mh

double SUSYSpectrum::mh[4]

private

Stores the tree-level Higgs spectrum.

Definition at line 305 of file SUSYSpectrum.h.

◆ mneu

gslpp::vector<double> SUSYSpectrum::mneu

private

Definition at line 330 of file SUSYSpectrum.h.

◆ Mneutralino

gslpp::matrix<gslpp::complex> SUSYSpectrum::Mneutralino

private

Definition at line 315 of file SUSYSpectrum.h.

◆ Msdown2

gslpp::matrix<gslpp::complex> SUSYSpectrum::Msdown2

private

Definition at line 325 of file SUSYSpectrum.h.

◆ Mselectron2

gslpp::matrix<gslpp::complex> SUSYSpectrum::Mselectron2

private

Definition at line 325 of file SUSYSpectrum.h.

◆ Msneutrino2

gslpp::matrix<gslpp::complex> SUSYSpectrum::Msneutrino2

private

Definition at line 325 of file SUSYSpectrum.h.

◆ Msup2

gslpp::matrix<gslpp::complex> SUSYSpectrum::Msup2

private

Stores the tree-level Up-squark, Down-squark, Sneutrino, and Slepton mass matrix.

Definition at line 325 of file SUSYSpectrum.h.

◆ mySUSY

const SUSY& SUSYSpectrum::mySUSY

private

Definition at line 300 of file SUSYSpectrum.h.

◆ N

gslpp::matrix<gslpp::complex> SUSYSpectrum::N

private

Definition at line 320 of file SUSYSpectrum.h.

◆ Rd

gslpp::matrix<gslpp::complex> SUSYSpectrum::Rd

private

Definition at line 335 of file SUSYSpectrum.h.

◆ Rl

gslpp::matrix<gslpp::complex> SUSYSpectrum::Rl

private

Definition at line 335 of file SUSYSpectrum.h.

◆ Rn

gslpp::matrix<gslpp::complex> SUSYSpectrum::Rn

private

Definition at line 335 of file SUSYSpectrum.h.

◆ Ru

gslpp::matrix<gslpp::complex> SUSYSpectrum::Ru

private

Stores the tree-level Up-squark, Down-squark, Sneutrino, and Slepton mixing matrices.

Definition at line 335 of file SUSYSpectrum.h.

◆ saeff

gslpp::complex SUSYSpectrum::saeff

private

Stores the Sine of tree-level CP-even mixing angle.

Definition at line 310 of file SUSYSpectrum.h.

◆ U

gslpp::matrix<gslpp::complex> SUSYSpectrum::U

private

Stores the tree-level Chargino and Neutralino mixing matrices.

Definition at line 320 of file SUSYSpectrum.h.

◆ V

gslpp::matrix<gslpp::complex> SUSYSpectrum::V

private

Definition at line 320 of file SUSYSpectrum.h.

The documentation for this class was generated from the following files:

Detailed Description

Public Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ SUSYSpectrum()

Member Function Documentation

◆ CalcChargino()

◆ CalcHiggs()

◆ CalcNeutralino()

◆ CalcSdown()

◆ CalcSelectron()

◆ CalcSneutrino()

◆ CalcSpectrum()

◆ CalcSup()

◆ getMch()

◆ getMchargino()

◆ getMHa()

◆ getMHh()

◆ getMHl()

◆ getMHp()

◆ getMneu()

◆ getMneutralino()

◆ getMsd2()

◆ getMsdown2()

◆ getMse2()

◆ getMselectron2()

◆ getMsn2()

◆ getMsneutrino2()

◆ getMsu2()

◆ getMsup2()

◆ getN()

◆ getRd()

◆ getRl()

◆ getRn()

◆ getRu()

◆ getU()

◆ getV()

◆ SortSfermionMasses()

Member Data Documentation

◆ m_sd2

◆ m_se2

◆ m_sn2

◆ m_su2

◆ mch

◆ Mchargino

◆ mh

◆ mneu

◆ Mneutralino

◆ Msdown2

◆ Mselectron2

◆ Msneutrino2

◆ Msup2

◆ mySUSY

◆ N

◆ Rd

◆ Rl

◆ Rn

◆ Ru

◆ saeff

◆ U

◆ V