a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
BXqllObservables.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2014 HEPfit Collaboration
3  *
4  *
5  * For the licensing terms see doc/COPYING.
6  */
7 
8 #include "BXqllObservables.h"
9 #include "StandardModel.h"
10 #include <TMath.h>
11 
12 /*******************************************************************************
13  * Observables *
14  * ****************************************************************************/
15 
16 
17 R_BXqll::R_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i)
18 : ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i)
19 {
20 // 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;
21 // lep = lep_i;
22 // quark = quark_i;
23 //
24  setParametersForObservable(myBXqll.initializeBXqllParameters());
25 }
26 
27 double R_BXqll::computeThValue()
28 {
29 // double q_min = getBinMin();
30 // double q_max = getBinMax();
31 // return (myBXqll.integrate_Rquark(q_min, q_max, LOWQ2));
32  return myBXqll.getR_LOWQ2(0.15);
33 }
34 
35 
36 HT_BXqll::HT_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i)
37 : ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i)
38 {
39 // 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;
40 // lep = lep_i;
41 // quark = quark_i;
42 //
43  setParametersForObservable(myBXqll.initializeBXqllParameters());
44 }
45 
46 double HT_BXqll::computeThValue()
47 {
48  double q_min = getBinMin();
49  double q_max = getBinMax();
50  return (myBXqll.integrateH("T", q_min, q_max));
51 }
52 
53 
54 HL_BXqll::HL_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i)
55 : ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i)
56 {
57 // 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;
58 // lep = lep_i;
59 // quark = quark_i;
60 //
61  setParametersForObservable(myBXqll.initializeBXqllParameters());
62 }
63 
64 double HL_BXqll::computeThValue()
65 {
66  double q_min = getBinMin();
67  double q_max = getBinMax();
68  return (myBXqll.integrateH("L", q_min, q_max));
69 }
70 
71 
72 HA_BXqll::HA_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i)
73 : ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i)
74 {
75 // 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;
76 // lep = lep_i;
77 // quark = quark_i;
78 //
79  setParametersForObservable(myBXqll.initializeBXqllParameters());
80 }
81 
82 double HA_BXqll::computeThValue()
83 {
84  double q_min = getBinMin();
85  double q_max = getBinMax();
86  return (myBXqll.integrateH("A", q_min, q_max));
87 }
88 
89 
90 BR_BXqll::BR_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i)
91 : ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i)
92 {
93 // 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;
94 // lep = lep_i;
95 // quark = quark_i;
96 //
97  setParametersForObservable(myBXqll.initializeBXqllParameters());
98 }
99 
100 double BR_BXqll::computeThValue()
101 {
102  double q_min = getBinMin();
103  double q_max = getBinMax();
104  return (myBXqll.integrateH("TL", q_min, q_max));
105 }
106 
107 
108 AFB_BXqll::AFB_BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i)
109 : ThObservable(SM_i), myBXqll(SM_i, quark_i, lep_i)
110 {
111 // 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;
112 // lep = lep_i;
113 // quark = quark_i;
114 //
115  setParametersForObservable(myBXqll.initializeBXqllParameters());
116 }
117 
118 double AFB_BXqll::computeThValue()
119 {
120  double q_min = getBinMin();
121  double q_max = getBinMax();
122  return (3. * myBXqll.integrateH("A", q_min, q_max) / 4.);
123 }
HT_BXqll::computeThValue
double computeThValue()
The binned observable in .
Definition: BXqllObservables.cpp:45
ThObservable::setParametersForObservable
void setParametersForObservable(std::vector< std::string > parametersForObservable_i)
A set method to get the parameters for the specific observable.
Definition: ThObservable.h:109
R_BXqll::R_BXqll
R_BXqll(const StandardModel &SM_i, QCD::quark quark_i, QCD::lepton lep_i)
Constructor.
Definition: BXqllObservables.cpp:16
BXqll::getR_LOWQ2
double getR_LOWQ2(double sh)
dGamma/ds for in the low dilepton invariant mass region.
Definition: BXqll.cpp:116
ThObservable::getBinMax
double getBinMax()
A get method to get the maximum value of the bin.
Definition: ThObservable.h:91
HL_BXqll::computeThValue
double computeThValue()
The binned observable in .
Definition: BXqllObservables.cpp:63
HT_BXqll::myBXqll
BXqll myBXqll
Definition: BXqllObservables.h:83
StandardModel.h
HT_BXqll::HT_BXqll
HT_BXqll(const StandardModel &SM_i, QCD::quark quark_i, QCD::lepton lep_i)
Constructor.
Definition: BXqllObservables.cpp:35
AFB_BXqll::myBXqll
BXqll myBXqll
Definition: BXqllObservables.h:230
R_BXqll::computeThValue
double computeThValue()
The binned observable in .
Definition: BXqllObservables.cpp:26
StandardModel
A model class for the Standard Model.
Definition: StandardModel.h:477
HA_BXqll::HA_BXqll
HA_BXqll(const StandardModel &SM_i, QCD::quark quark_i, QCD::lepton lep_i)
Constructor.
Definition: BXqllObservables.cpp:71
BXqll::integrateH
double integrateH(std::string obs, double q_min, double q_max)
The integral of each observable as defined in .
Definition: BXqll.cpp:241
HA_BXqll::myBXqll
BXqll myBXqll
Definition: BXqllObservables.h:157
BR_BXqll::myBXqll
BXqll myBXqll
Definition: BXqllObservables.h:193
BR_BXqll::computeThValue
double computeThValue()
The binned observable .
Definition: BXqllObservables.cpp:99
HL_BXqll::HL_BXqll
HL_BXqll(const StandardModel &SM_i, QCD::quark quark_i, QCD::lepton lep_i)
Constructor.
Definition: BXqllObservables.cpp:53
BXqll::initializeBXqllParameters
std::vector< std::string > initializeBXqllParameters()
A method for initializing the parameters necessary for BXqll.
Definition: BXqll.cpp:39
BXqllObservables.h
QCD::quark
quark
An enum type for quarks.
Definition: QCD.h:323
BR_BXqll::BR_BXqll
BR_BXqll(const StandardModel &SM_i, QCD::quark quark_i, QCD::lepton lep_i)
Constructor.
Definition: BXqllObservables.cpp:89
R_BXqll::myBXqll
BXqll myBXqll
Definition: BXqllObservables.h:46
ThObservable
A class for a model prediction of an observable.
Definition: ThObservable.h:25
ThObservable::getBinMin
double getBinMin()
A get method to get the minimum value of the bin.
Definition: ThObservable.h:82
AFB_BXqll::computeThValue
double computeThValue()
The binned observable .
Definition: BXqllObservables.cpp:117
HA_BXqll::computeThValue
double computeThValue()
The binned observable in .
Definition: BXqllObservables.cpp:81
HL_BXqll::myBXqll
BXqll myBXqll
Definition: BXqllObservables.h:120
AFB_BXqll::AFB_BXqll
AFB_BXqll(const StandardModel &SM_i, QCD::quark quark_i, QCD::lepton lep_i)
Constructor.
Definition: BXqllObservables.cpp:107
QCD::lepton
lepton
An enum type for leptons.
Definition: QCD.h:310