Belle II Software development
Public Member Functions | Private Attributes | List of all members
CDCDedxHadSat Class Reference
Packages » cdc

Class to hold the hadron saturation functions. More...

#include <CDCDedxHadSat.h>

Public Member Functions

void setParameters ()
 set the parameters
 
void setParameters (std::string infile)
 set the parameters from file
 
double D2I (double cosTheta, double D) const
 hadron saturation parameterization part 1
 
double I2D (double cosTheta, double I) const
 hadron saturation parameterization part 2
 
double D2I (double cosTheta, double D, double alpha, double gamma, double delta, double power, double ratio) const
 hadron saturation parameterization part 1
 
double I2D (double cosTheta, double I, double alpha, double gamma, double delta, double power, double ratio) const
 hadron saturation parameterization part 2
 

Private Attributes

double m_alpha
 the alpha parameter for the hadron saturation correction
 
double m_gamma
 the gamma parameter for the hadron saturation correction
 
double m_delta
 the delta parameter for the hadron saturation correction
 
double m_power
 the power parameter for the hadron saturation correction
 
double m_ratio
 the ratio parameter for the hadron saturation correction
 
DBObjPtr< CDCDedxHadronCor > m_DBHadronCor
 db object for dE/dx hadron saturation parameters
 

Detailed Description

Class to hold the hadron saturation functions.

Definition at line 31 of file CDCDedxHadSat.h.

Member Function Documentation

◆ D2I() [1/2]

double D2I ( double  cosTheta,
double  D = 1 
) const

hadron saturation parameterization part 1

Definition at line 48 of file CDCDedxHadSat.cc.

49{
50
51 double absCosTheta = fabs(cosTheta);
52 double projection = pow(absCosTheta, m_power) + m_delta;
53
54 if (projection == 0) {
55 B2WARNING("Something wrong with dE/dx hadron constants!");
56 return D;
57 }
58
59 double chargeDensity = D / projection;
60 double numerator = 1 + m_alpha * chargeDensity;
61 double denominator = 1 + m_gamma * chargeDensity;
62 if (denominator == 0) {
63 B2WARNING("Something wrong with dE/dx hadron constants!");
64 return D;
65 }
66 double I = D * m_ratio * numerator / denominator;
67
68 return I;
69}
Belle2::CDCDedxHadSat::m_delta
double m_delta
the delta parameter for the hadron saturation correction
Definition: CDCDedxHadSat.h:69
Belle2::CDCDedxHadSat::m_ratio
double m_ratio
the ratio parameter for the hadron saturation correction
Definition: CDCDedxHadSat.h:71
Belle2::CDCDedxHadSat::m_gamma
double m_gamma
the gamma parameter for the hadron saturation correction
Definition: CDCDedxHadSat.h:68
Belle2::CDCDedxHadSat::m_alpha
double m_alpha
the alpha parameter for the hadron saturation correction
Definition: CDCDedxHadSat.h:67
Belle2::CDCDedxHadSat::m_power
double m_power
the power parameter for the hadron saturation correction
Definition: CDCDedxHadSat.h:70

◆ D2I() [2/2]

double D2I ( double  cosTheta,
double  D,
double  alpha,
double  gamma,
double  delta,
double  power,
double  ratio 
) const

hadron saturation parameterization part 1

Definition at line 109 of file CDCDedxHadSat.cc.

110{
111
112 double absCosTheta = fabs(cosTheta);
113 double projection = pow(absCosTheta, power) + delta;
114 if (projection == 0) {
115 B2INFO("\t HadronSaturation: D2I Something wrong with dE/dx hadron constants!");
116 return D;
117 }
118
119 double chargeDensity = D / projection;
120 double numerator = 1 + alpha * chargeDensity;
121 double denominator = 1 + gamma * chargeDensity;
122 if (denominator == 0) {
123 B2INFO("\t HadronSaturation: D2I Something wrong with dE/dx hadron constants!");
124 return D;
125 }
126
127 double I = D * ratio * numerator / denominator;
128
129 return I;
130}

◆ I2D() [1/2]

double I2D ( double  cosTheta,
double  I = 1 
) const

hadron saturation parameterization part 2

Definition at line 71 of file CDCDedxHadSat.cc.

72{
73
74 double absCosTheta = fabs(cosTheta);
75 double projection = pow(absCosTheta, m_power) + m_delta;
76
77 if (projection == 0 || m_ratio == 0) {
78 B2WARNING("Something wrong with dE/dx hadron constants!");
79 return I;
80 }
81 double a = m_alpha / projection;
82 double b = 1 - m_gamma / projection * (I / m_ratio);
83 double c = -1.0 * I / m_ratio;
84
85 if (b == 0 && a == 0) {
86 B2WARNING("both a and b coefficients for hadron correction are 0");
87 return I;
88 }
89
90 double discr = b * b - 4.0 * a * c;
91 if (discr < 0) {
92 B2WARNING("negative discriminant; return uncorrected value");
93 return I;
94 }
95
96 double D = (a != 0) ? (-b + sqrt(discr)) / (2.0 * a) : -c / b;
97 if (D < 0) {
98 B2WARNING("D is less 0! will try another solution");
99 D = (a != 0) ? (-b - sqrt(discr)) / (2.0 * a) : -c / b;
100 if (D < 0) {
101 B2WARNING("D is still less 0! just return uncorrected value");
102 return I;
103 }
104 }
105 return D;
106}
Belle2::sqrt
double sqrt(double a)
sqrt for double
Definition: beamHelpers.h:28

◆ I2D() [2/2]

double I2D ( double  cosTheta,
double  I,
double  alpha,
double  gamma,
double  delta,
double  power,
double  ratio 
) const

hadron saturation parameterization part 2

Definition at line 133 of file CDCDedxHadSat.cc.

134{
135
136 double absCosTheta = fabs(cosTheta);
137 double projection = pow(absCosTheta, power) + delta;
138
139 if (projection == 0 || ratio == 0) {
140 B2INFO("\t HadronSaturation: I2D Something wrong with dE/dx hadron constants!");
141 return I;
142 }
143
144 double a = alpha / projection;
145 double b = 1 - (gamma / projection) * (I / ratio);
146 double c = -1.0 * I / ratio;
147
148 if (b == 0 && a == 0) {
149 B2INFO("\t HadronSaturation: both a and b coefficients for hadron correction are 0");
150 return I;
151 }
152
153 double discr = b * b - 4.0 * a * c;
154 if (discr < 0) {
155 B2INFO("negative discriminant; return uncorrected value");
156 return I;
157 }
158
159 double D = (a != 0) ? (-b + sqrt(discr)) / (2.0 * a) : -c / b;
160
161 if (D < 0) {
162 B2INFO("D is less 0! will try another solution");
163 D = (a != 0) ? (-b - sqrt(discr)) / (2.0 * a) : -c / b;
164 if (D < 0) {
165 B2INFO("D is still less 0! just return uncorrected value");
166 return I;
167 }
168 }
169
170 return D;
171}

◆ setParameters() [1/2]

void setParameters ( )

set the parameters

Definition at line 32 of file CDCDedxHadSat.cc.

33{
34
35 if (!m_DBHadronCor || m_DBHadronCor->getSize() == 0) {
36 B2FATAL("No hadron correction parameters!");
37 } else {
38 m_alpha = m_DBHadronCor->getHadronPar(0);
39 m_gamma = m_DBHadronCor->getHadronPar(1);
40 m_delta = m_DBHadronCor->getHadronPar(2);
41 m_power = m_DBHadronCor->getHadronPar(3);
42 m_ratio = m_DBHadronCor->getHadronPar(4);
43 }
44}
Belle2::CDCDedxHadSat::m_DBHadronCor
DBObjPtr< CDCDedxHadronCor > m_DBHadronCor
db object for dE/dx hadron saturation parameters
Definition: CDCDedxHadSat.h:73

◆ setParameters() [2/2]

void setParameters ( std::string  infile)

set the parameters from file

Definition at line 13 of file CDCDedxHadSat.cc.

14{
15
16 B2INFO("\n\t Hadron Saturation: Using parameters from file --> " << infile);
17
18 std::ifstream fin;
19 fin.open(infile.c_str());
20
21 if (!fin.good()) B2FATAL("\tWARNING: CANNOT FIND parameters.txt!");
22
23 fin >> m_alpha ;
24 fin >> m_gamma ;
25 fin >> m_delta ;
26 fin >> m_power ;
27 fin >> m_ratio ;
28
29 fin.close();
30}

Member Data Documentation

◆ m_alpha

double m_alpha
private

the alpha parameter for the hadron saturation correction

Definition at line 67 of file CDCDedxHadSat.h.

◆ m_DBHadronCor

DBObjPtr<CDCDedxHadronCor> m_DBHadronCor
private

db object for dE/dx hadron saturation parameters

Definition at line 73 of file CDCDedxHadSat.h.

◆ m_delta

double m_delta
private

the delta parameter for the hadron saturation correction

Definition at line 69 of file CDCDedxHadSat.h.

◆ m_gamma

double m_gamma
private

the gamma parameter for the hadron saturation correction

Definition at line 68 of file CDCDedxHadSat.h.

◆ m_power

double m_power
private

the power parameter for the hadron saturation correction

Definition at line 70 of file CDCDedxHadSat.h.

◆ m_ratio

double m_ratio
private

the ratio parameter for the hadron saturation correction

Definition at line 71 of file CDCDedxHadSat.h.

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