CKM.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef CKM_H
#define CKM_H
 
#include <math.h>
#include "gslpp.h"
 
class CKM {
public:
    CKM();
 
    void computeCKMwithWolfenstein(double Lambda_v, double A_v, double Rho_v, double Eta_v);
    
    void computeCKM(double Vus_v, double Vcb_v, double Vub_v, double gamma_v, bool useVud = false);
 
    void computeCKMfromAngles(double s12_in, double s23_in, double s13_in, double delta_in);
 
    const gslpp::matrix<gslpp::complex> getCKM() const 
    {
        return V;
    }
 
    // Wolfenstein parameters
    
    const double getRhoBar() const 
    {
        return Rho;
    }
    
    const double getEtaBar() const
    {
        return Eta;
    }
    
    const double getLambda() const
    {
        return Lambda;
    }
    
    const double getA() const
    {
        return A;
    }
 
    const double getRho() const
    {
        return (s13 * cos(delta) / s12 / s23);
    }
 
    const double getEta() const
    {
        return (s13 * sin(delta) / s12 / s23);
    }
 
    // Gilman parameterization
 
    const double gets12() const
    {
        return s12;
    }
 
    const double gets13() const
    {
        return s13;
    }
 
    const double gets23() const
    {
        return s23;
    }
 
    const double getc12() const
    {
        return c12;
    }
 
    const double getc23() const
    {
        return c23;
    }
 
    const double getc13() const
    {
        return c13;
    }
 
    const double getdelta() const
    {
        return delta;
    }
 
    // J_CP
 
    const double getJcp() const
    {
        return Eta * pow(A * pow(Lambda, 3), 2);
    }
 
    //Complex values of CKM elements
 
    const gslpp::complex getV_ud() const
    {
        return V(0, 0);
    }
    
    const gslpp::complex getV_us() const
    {
        return V(0, 1);
    }
 
    const gslpp::complex getV_ub() const
    {
        return V(0, 2);
    }
 
    const gslpp::complex getV_cd() const
    {
        return V(1, 0);
    }
 
    const gslpp::complex getV_cs() const
    {
        return V(1, 1);
    }
 
    const gslpp::complex getV_cb() const
    {
        return V(1, 2);
    }
 
    const gslpp::complex getV_td() const
    {
        return V(2, 0);
    }
 
    const gslpp::complex getV_ts() const
    {
        return V(2, 1);
    }
 
    const gslpp::complex getV_tb() const
    {
        return V(2, 2);
    }
    
    const double computeBeta() const;
 
    const double computeGamma() const;
 
    const double computeAlpha() const;
 
    const double computeBetas() const;
 
    const gslpp::complex computelamt() const;
 
    const gslpp::complex computelamc() const;
 
    const gslpp::complex computelamu() const;
 
    const gslpp::complex computelamt_d() const;
 
    const gslpp::complex computelamc_d() const;
 
    const gslpp::complex computelamu_d() const;
 
    const gslpp::complex computelamt_s() const;
 
    const gslpp::complex computelamc_s() const;
 
    const gslpp::complex computelamu_s() const;
 
    const double computeRt() const;
 
    const double computeRts() const;
 
    const double computeRb() const;
    
private:
    void computeCKMfromAngles();
 
    double Rho, Eta, Lambda, A; 
    double s12, s13, s23; 
    double c12, c23, c13; 
    double delta; 
 
    gslpp::matrix<gslpp::complex> V; 
 
};
 
#endif  /* CKM_H */