BXqllObservables.cpp

Go to the documentation of this file.

/* 
 * Copyright (C) 2014 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#include "BXqllObservables.h"
#include "StandardModel.h"
#include <TMath.h>
 
/*******************************************************************************
 * Observables                                                                 *
 * ****************************************************************************/
 
 
R_BXqll::R_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i) 
: ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i) 
{  
//    if (SM.getModelName().compare("StandardModel") != 0) std::cout << "\nWARNING: B to Xq l+ l-: R_BXqll not implemented in: " + SM.getModelName() + " model, returning Standard Model value.\n" << std::endl;
//    lep = lep_i;
//    quark = quark_i;
//    
    setParametersForObservable(myBXqll.initializeBXqllParameters());
}
 
double R_BXqll::computeThValue()
{ 
//    double q_min = getBinMin();
//    double q_max = getBinMax();
//    return (myBXqll.integrate_Rquark(q_min, q_max, LOWQ2));
    return myBXqll.getR_LOWQ2(0.15);
}
 
 
HT_BXqll::HT_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i) 
: ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i) 
{  
//    if (SM.getModelName().compare("StandardModel") != 0) std::cout << "\nWARNING: B to Xq l+ l-: Rlow_BXqll not implemented in: " + SM.getModelName() + " model, returning Standard Model value.\n" << std::endl;
//    lep = lep_i;
//    quark = quark_i;
//    
    setParametersForObservable(myBXqll.initializeBXqllParameters());
}
 
double HT_BXqll::computeThValue()
{
    double q_min = getBinMin();
    double q_max = getBinMax();
    return (myBXqll.integrateH("T", q_min, q_max));
}
 
 
HL_BXqll::HL_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i) 
: ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i) 
{  
//    if (SM.getModelName().compare("StandardModel") != 0) std::cout << "\nWARNING: B to Xq l+ l-: Rlow_BXqll not implemented in: " + SM.getModelName() + " model, returning Standard Model value.\n" << std::endl;
//    lep = lep_i;
//    quark = quark_i;
//    
    setParametersForObservable(myBXqll.initializeBXqllParameters());
}
 
double HL_BXqll::computeThValue()
{
    double q_min = getBinMin();
    double q_max = getBinMax();
    return (myBXqll.integrateH("L", q_min, q_max));
}
 
 
HA_BXqll::HA_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i) 
: ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i) 
{  
//    if (SM.getModelName().compare("StandardModel") != 0) std::cout << "\nWARNING: B to Xq l+ l-: Rlow_BXqll not implemented in: " + SM.getModelName() + " model, returning Standard Model value.\n" << std::endl;
//    lep = lep_i;
//    quark = quark_i;
//    
    setParametersForObservable(myBXqll.initializeBXqllParameters());
}
 
double HA_BXqll::computeThValue()
{
    double q_min = getBinMin();
    double q_max = getBinMax();
    return (myBXqll.integrateH("A", q_min, q_max));
}
 
 
BR_BXqll::BR_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i) 
: ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i) 
{  
//    if (SM.getModelName().compare("StandardModel") != 0) std::cout << "\nWARNING: B to Xq l+ l-: Rlow_BXqll not implemented in: " + SM.getModelName() + " model, returning Standard Model value.\n" << std::endl;
//    lep = lep_i;
//    quark = quark_i;
//    
    setParametersForObservable(myBXqll.initializeBXqllParameters());
}
 
double BR_BXqll::computeThValue()
{
    double q_min = getBinMin();
    double q_max = getBinMax();
    return (myBXqll.integrateH("TL", q_min, q_max));
}
 
 
AFB_BXqll::AFB_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i) 
: ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i) 
{  
//    if (SM.getModelName().compare("StandardModel") != 0) std::cout << "\nWARNING: B to Xq l+ l-: Rlow_BXqll not implemented in: " + SM.getModelName() + " model, returning Standard Model value.\n" << std::endl;
//    lep = lep_i;
//    quark = quark_i;
//    
    setParametersForObservable(myBXqll.initializeBXqllParameters());
}
 
double AFB_BXqll::computeThValue()
{
    double q_min = getBinMin();
    double q_max = getBinMax();
    return (3. * myBXqll.integrateH("A", q_min, q_max) / 4.);
}