definition of the couplings for the different BGL models More...

#include <Formulas.h>

Collaboration diagram for BGLmodels::Mixes:

Public Member Functions
	Mixes (ex tanb0, ex cp0, int genL=2, int genQ=2, int lup=0, int qup=0, int mssm=0)

ex	mass (const Fermion &f) const

double	massnum (const Fermion &f) const

void	appendtolst (lst &reps) const

Public Attributes
lst	reps

vector< Matrix >	V

multivector< Matrix, 2 >	M

multivector< Matrix, 2 >	N

multivector< Matrix, 2 >	VN

multivector< Matrix, 2 >	N_

multivector< Matrix, 2 >	VN_

lst	replacements

ex	cp

ex	tanb

Detailed Description

definition of the couplings for the different BGL models

Definition at line 105 of file Formulas.h.

Constructor & Destructor Documentation

BGLmodels::Mixes::Mixes	(	ex	tanb0,
		ex	cp0,
		int	genL = `2`,
		int	genQ = `2`,
		int	lup = `0`,
		int	qup = `0`,
		int	mssm = `0`
	)

inline

Definition at line 108 of file Formulas.h.

References std::Matrix::conjugate(), BGLmodels::iDown, BGLmodels::iUp, BGLmodels::tLepton, and BGLmodels::tQuark.

                                                                              :M(Matrix(),2,2),N(Matrix(),2,2),VN(Matrix(),2,2),N_(Matrix(),2,2),VN_(Matrix(),2,2),cp(cp0)
   {
           tanb=tanb0;
           M[tLepton][iDown]=Matrix(possymbol("m_e"),possymbol("m_\\mu"),possymbol("m_\\tau"));
           M[tQuark][iUp]=Matrix(possymbol("m_u"),possymbol("m_c"),possymbol("m_t"));
           M[tQuark][iDown]=Matrix(possymbol("m_d"),possymbol("m_s"),possymbol("m_b"));
           const char * ln[3]={"1","2","3"};
           const char * ll[3]={"e","\\mu","\\tau"};
           const char * lu[3]={"u","c","t"};
           const char * ld[3]={"d","s","b"};
 
         V.push_back(Matrix("U",ln,ll));
         V.push_back(Matrix("V",lu,ld));
         
         
         int up[2];
         up[0]=lup;
         up[1]=qup;
         
         vector< Matrix > delta;
         
           vector<int> gL(3,0);
           gL[genL]=1;
           //Leptons
           delta.push_back(Matrix(gL[0],gL[1],gL[2]));
           //Quarks
           vector<int> gQ(3,0);
           gQ[genQ]=1;
           delta.push_back(Matrix(gQ[0],gQ[1],gQ[2]));
           
       for(uint i=0;i<2;i++){
                 if(mssm){
                 //Nu
             N_[i][0]=0;
         //Nd
         N_[i][1]=0;
         //VNd
         VN_[i][1]=Matrix(tanb)*V[i];
         //Nu*V
         VN_[i][0]=Matrix(1/tanb)*V[i];
         }
         else if(up[i]){
                 //Nu
                 N_[i][0]=Matrix(tanb)+Matrix(-tanb-1/tanb)*V[i]*delta[i]*V[i].conjugate();
         //Nd
         N_[i][1]=Matrix(tanb)+Matrix(-tanb-1/tanb)*delta[i];
         //VNd
         VN_[i][1]=Matrix(tanb)*V[i]+Matrix(-tanb-1/tanb)*V[i]*delta[i];
         //Nu*V
         VN_[i][0]=Matrix(-1)*(Matrix(tanb)*V[i]+Matrix(-tanb-1/tanb)*V[i]*delta[i]);
                 }else{
                 //Nu
                 N_[i][0]=Matrix(tanb)+Matrix(-tanb-1/tanb)*delta[i];
         //Nd
         N_[i][1]=Matrix(tanb)+Matrix(-tanb-1/tanb)*V[i].conjugate()*delta[i]*V[i];
         //VNd
         VN_[i][1]=Matrix(tanb)*V[i]+Matrix(-tanb-1/tanb)*delta[i]*V[i];
         //Nu*V
         VN_[i][0]=Matrix(-1)*(Matrix(tanb)*V[i]+Matrix(-tanb-1/tanb)*delta[i]*V[i]);
                 }
                 
                 N[i][0]=N_[i][0]*M[i][0];
                 N[i][1]=N_[i][1]*M[i][1];
                 VN[i][1]=VN_[i][1]*M[i][1];
                 VN[i][0]=M[i][0]*VN_[i][0];
       }
       appendtolst(replacements);  
       
   
   }

Here is the call graph for this function:

Member Function Documentation

void BGLmodels::Mixes::appendtolst ( lst & reps ) const

inline

Definition at line 181 of file Formulas.h.

                                     {//,vector<ex>& var_errors
           reps.append(M[0][1][0][0]==0.510998910e-3);
           reps.append(M[0][1][1][1]==105.6583715e-3);
           reps.append(M[0][1][2][2]==1776.82e-3);
           
           reps.append(M[1][0][0][0]==2.4e-3);
           reps.append(M[1][0][1][1]==1.29);
           reps.append(M[1][0][2][2]==172.9);
           
           reps.append(M[1][1][0][0]==5.3e-3);
           reps.append(M[1][1][1][1]==95e-3);
           reps.append(M[1][1][2][2]==4.2);
           
           vector< Matrix > Vn;
         Vn.push_back(Matrix(33.6*M_PI/180,9.11*M_PI/180,40.4*M_PI/180,M_PI/4).conjugate());
         Vn.push_back(Matrix(13.04*M_PI/180,0.201*M_PI/180,2.38*M_PI/180,1.2));
         for(uint i=0; i<2;i++)
                 for(uint j=0; j<3;j++)
                         for(uint k=0; k<3;k++)
                                 reps.append(V[i][j][k]==Vn[i][j][k]);
           }