ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Private Member Functions | Private Attributes | List of all members
OnlineLumiCalibrator Class Reference

Utility class to apply calibrations from /TDAQ/OLC/CALIBRATIONS folder. More...

#include <OnlineLumiCalibrator.h>

Collaboration diagram for OnlineLumiCalibrator:

Public Member Functions

 OnlineLumiCalibrator ()
 OnlineLumiCalibrator (const OnlineLumiCalibrator &)=default
OnlineLumiCalibrator & operator= (const OnlineLumiCalibrator &)=default
 OnlineLumiCalibrator (OnlineLumiCalibrator &&)
OnlineLumiCalibrator & operator= (OnlineLumiCalibrator &&)
bool setCalibration (const coral::AttributeList &attrList)
bool calibrateMu (float rawLumi, float &calMu) const
bool calibrateLumi (float rawLumi, float &calLumi) const
bool calibrateMu (const std::vector< float > &rawLumi, std::vector< float > &calMu) const
bool calibrateLumi (const std::vector< float > &rawLumi, std::vector< float > &calLumi) const
float getMuToLumi () const
MsgStream & dump (MsgStream &) const

Private Member Functions

float getMuVis (float rawPerBX, float sigo, float siga) const
float getMuVis2 (float rawPerBX, float sigo, float siga) const

Private Attributes

unsigned int m_nPar
std::string m_fType
float m_muToLumi
std::vector< float > m_parVec

Detailed Description

Utility class to apply calibrations from /TDAQ/OLC/CALIBRATIONS folder.

Author
E.Torrence

Definition at line 20 of file OnlineLumiCalibrator.h.

Constructor & Destructor Documentation

◆ OnlineLumiCalibrator() [1/3]

OnlineLumiCalibrator::OnlineLumiCalibrator ( )

Definition at line 21 of file OnlineLumiCalibrator.cxx.

21 :
22 m_nPar(0),
23 m_fType(""),
24 m_muToLumi(-1.) // Nonsensical intialization value
25{
26 m_parVec.clear();
27}
OnlineLumiCalibrator::m_parVec
std::vector< float > m_parVec
Definition OnlineLumiCalibrator.h:53

◆ OnlineLumiCalibrator() [2/3]

OnlineLumiCalibrator::OnlineLumiCalibrator ( const OnlineLumiCalibrator & )
default

◆ OnlineLumiCalibrator() [3/3]

OnlineLumiCalibrator::OnlineLumiCalibrator ( OnlineLumiCalibrator && other)

Definition at line 29 of file OnlineLumiCalibrator.cxx.

30 : m_nPar (other.m_nPar),
31 m_fType (std::move (other.m_fType)),
32 m_muToLumi (other.m_muToLumi),
33 m_parVec (std::move (other.m_parVec))
34{
35}

Member Function Documentation

◆ calibrateLumi() [1/2]

bool OnlineLumiCalibrator::calibrateLumi ( const std::vector< float > & rawLumi,
std::vector< float > & calLumi ) const

Definition at line 87 of file OnlineLumiCalibrator.cxx.

88{
89 calLumi.clear();
90 bool error = false;
91 float calValue;
92 for (float val : rawLumi) {
93 calValue = 0;
94 if (!calibrateLumi(val, calValue)) {
95 error = true;
96 calLumi.push_back(0.);
97 } else {
98 calLumi.push_back(calValue);
99 }
100 }
101 return error;
102}
OnlineLumiCalibrator::calibrateLumi
bool calibrateLumi(float rawLumi, float &calLumi) const
Definition OnlineLumiCalibrator.cxx:125
get_generator_info.error
error
Definition get_generator_info.py:40

◆ calibrateLumi() [2/2]

bool OnlineLumiCalibrator::calibrateLumi ( float rawLumi,
float & calLumi ) const

Definition at line 125 of file OnlineLumiCalibrator.cxx.

126{
127 calLumi = 0;
128 float calMu = 0.;
129 if (!calibrateMu(rawLumi, calMu)) return false;
130 calLumi = calMu * m_muToLumi;
131 return true;
132}
OnlineLumiCalibrator::calibrateMu
bool calibrateMu(float rawLumi, float &calMu) const
Definition OnlineLumiCalibrator.cxx:136

◆ calibrateMu() [1/2]

bool OnlineLumiCalibrator::calibrateMu ( const std::vector< float > & rawLumi,
std::vector< float > & calMu ) const

Definition at line 106 of file OnlineLumiCalibrator.cxx.

107{
108 calMu.clear();
109 bool error = false;
110 float calValue;
111 for (float val : rawLumi) {
112 calValue = 0;
113 if (!calibrateMu(val, calValue)) {
114 error = true;
115 calMu.push_back(0.);
116 } else {
117 calMu.push_back(calValue);
118 }
119 }
120 return error;
121}

◆ calibrateMu() [2/2]

bool OnlineLumiCalibrator::calibrateMu ( float rawLumi,
float & calMu ) const

Definition at line 136 of file OnlineLumiCalibrator.cxx.

137{
138 calMu = 0.;
139
140 if (m_fType == "Polynomial") {
141
142 // Check parameter length
143 if (m_parVec.size() < 1) return false;
144
145 unsigned int nrange = m_parVec[0];
146
147 // Check parameter size again
148 if (m_parVec.size() < (8*nrange + 1)) return false;
149
150 for (unsigned int i=0; i<nrange; i++) {
151 float rmin = m_parVec[8*i+1];
152 float rmax = m_parVec[8*i+2];
153 if (rawLumi < rmax and rawLumi >= rmin) {
154 calMu += m_parVec[8*i+3] * pow(rawLumi, 0);
155 calMu += m_parVec[8*i+4] * pow(rawLumi, 1);
156 calMu += m_parVec[8*i+5] * pow(rawLumi, 2);
157 calMu += m_parVec[8*i+6] * pow(rawLumi, 3);
158 calMu += m_parVec[8*i+7] * pow(rawLumi, 4);
159 calMu += m_parVec[8*i+8] * pow(rawLumi, 5);
160 break;
161 }
162 }
163 return true;
164 }
165
166 if (m_fType == "Logarithm") {
167
168 // Check parameter length
169 if (m_parVec.size() != 1) return false;
170
171 // Check input for physically allowed range
172 if ((1.-rawLumi) <= 0.) return false;
173
174 // Convert negative luminosity to zero
175 if (rawLumi < 0.) return false;
176
177 calMu = -m_parVec[0] * log(1.-rawLumi);
178 return true;
179 }
180
181 if (m_fType == "HitLogarithm") {
182
183 // Check parameter length
184 if (m_parVec.size() != 4) return false;
185
186 // Channel count must be > 0
187 if (m_parVec[1] <= 0.) return false;
188
189 // Check input value for physically allowed range
190 if ((1.-rawLumi/m_parVec[1]) <= 0.) return false;
191
192 // Convert negative luminosity to zero
193 if (rawLumi < 0.) return false;
194
195 calMu = -m_parVec[0] * log(1.-rawLumi/m_parVec[1])/(1.-m_parVec[2]);
196 return true;
197 }
198
199 if (m_fType == "LookupTable_EventAND_Lin") {
200
201 // Check parameter size
202 if (m_parVec.size() != 6) return false;
203
204 if (rawLumi < 0.) return false;
205
206 float sigo = m_parVec[0];
207 float siga = m_parVec[1];
208 calMu = m_parVec[5] * getMuVis(rawLumi, sigo, siga);
209
210 if (calMu < 0.) { // Indicates an error, try again
211 calMu = m_parVec[5] * getMuVis2(rawLumi, sigo, siga);
212 }
213
214 if (calMu < 0.) { // Indicates an error
215 calMu = 0.;
216 return false;
217 }
218
219 return true;
220 }
221
222 if (m_fType == "LookupTable_EventAND_Log") {
223
224 // Check parameter size
225 if (m_parVec.size() != 8) return false;
226
227 if (rawLumi < 0.) return false;
228
229 float sigo = m_parVec[0];
230 float siga = m_parVec[1];
231 calMu = m_parVec[7] * getMuVis(rawLumi, sigo, siga);
232
233 if (calMu < 0.) { // Indicates an error, try again
234 calMu = m_parVec[7] * getMuVis2(rawLumi, sigo, siga);
235 }
236
237 if (calMu < 0.) { // Indicates an error
238 calMu = 0.;
239 return false;
240 }
241
242 return true;
243 }
244
245 // Unknown type
246 return false;
247}
OnlineLumiCalibrator::getMuVis2
float getMuVis2(float rawPerBX, float sigo, float siga) const
Definition OnlineLumiCalibrator.cxx:285
OnlineLumiCalibrator::getMuVis
float getMuVis(float rawPerBX, float sigo, float siga) const
Definition OnlineLumiCalibrator.cxx:251
conifer::pow
constexpr int pow(int x)
Definition conifer.h:27

◆ dump()

MsgStream & OnlineLumiCalibrator::dump ( MsgStream & stream) const

Definition at line 322 of file OnlineLumiCalibrator.cxx.

322 {
323 stream << m_fType << " Npar = " << m_nPar << " MuToLumi = " << m_muToLumi << " ParVec: " << m_parVec;
324 return stream;
325}

◆ getMuToLumi()

float OnlineLumiCalibrator::getMuToLumi ( ) const

Definition at line 80 of file OnlineLumiCalibrator.cxx.

81{
82 return m_muToLumi;
83}

◆ getMuVis()

float OnlineLumiCalibrator::getMuVis ( float rawPerBX,
float sigo,
float siga ) const
private

Definition at line 251 of file OnlineLumiCalibrator.cxx.

252{
253
254 //std::cout << "getMuVis("<<rawPerBX<<","<<sigo<<","<<siga<<")"<<std::endl;
255
256 float mu, y, dy;
257 float munew = 0.01;
258 float a = (sigo/siga + 1) / 2.;
259 float b = sigo/siga;
260
261 // Set a fixed number of cycles, but break if we converge faster
262 for (int i=0; i<30; i++) {
263 mu = munew;
264 y = rawPerBX - 1. - exp(-b * mu) + 2. * exp(-a * mu);
265 dy = b * exp(-b * mu) - 2. * a * exp(-a * mu);
266 munew = mu-y/dy;
267
268 //std::cout << i <<" - munew: " << munew << " mu:" << mu << " diff:"
269 // << fabs(munew-mu) << std::endl;
270
271 // Protect against unphysical values
272 if (munew <= 0.) return -1.;
273 if (fabs(munew-mu)/munew < 1.E-5) break;
274 }
275
276 return munew;
277}
a
static Double_t a
Definition LArPhysWaveHECTool.cxx:38
y
#define y
CaloNoise_fillDB.mu
int mu
Definition CaloNoise_fillDB.py:173

◆ getMuVis2()

float OnlineLumiCalibrator::getMuVis2 ( float rawPerBX,
float sigo,
float siga ) const
private

Definition at line 285 of file OnlineLumiCalibrator.cxx.

286{
287 float muvl=1e-10;
288 float muvu=100.;
289 float muvm=10.;
290 float sr=sigo/siga;
291 float rbxl=rpbx(sr,muvl);
292 float rbxu=rpbx(sr,muvu);
293 float rbxm=rpbx(sr,muvm);
294
295 // Check that starting value is in the valid range
296 if (rawPerBX < rbxl || rawPerBX > rbxu) return -1.;
297
298 int i=0;
299 while (++i <= 50) {
300 if (rbxl<rawPerBX && rbxm>rawPerBX) {
301 rbxu=rbxm;
302 muvu=muvm;
303 muvm=0.5*(muvu+muvl);
304 } else {
305 rbxl=rbxm;
306 muvl=muvm;
307 muvm=0.5*(muvu+muvl);
308 }
309
310 rbxm = rpbx(sr, muvm);
311 //std::cout << i <<" - munew: " << muvu << " mu:" << muvl << " diff:" << fabs(muvu-muvl) << std::endl;
312
313 if ((muvu-muvl)/muvl < 1e-5) return muvm;
314 }
315
316 // Didn't converge
317 return -1.;
318}
rpbx
float rpbx(float sr, float mu)
Definition OnlineLumiCalibrator.cxx:280

◆ operator=() [1/2]

OnlineLumiCalibrator & OnlineLumiCalibrator::operator= ( const OnlineLumiCalibrator & )
default

◆ operator=() [2/2]

OnlineLumiCalibrator & OnlineLumiCalibrator::operator= ( OnlineLumiCalibrator && other)

Definition at line 39 of file OnlineLumiCalibrator.cxx.

40{
41 if (this != &other) {
42 m_nPar = other.m_nPar;
43 m_fType = std::move (other.m_fType);
44 m_muToLumi = other.m_muToLumi;
45 m_parVec = std::move (other.m_parVec);
46 }
47 return *this;
48}

◆ setCalibration()

bool OnlineLumiCalibrator::setCalibration ( const coral::AttributeList & attrList)

Definition at line 53 of file OnlineLumiCalibrator.cxx.

54{
55 if (attrList["NumOfParameters"].isNull()) return false;
56 m_nPar = attrList["NumOfParameters"].data<uint32_t>();
57
58 if (attrList["Function"].isNull()) return false;
59 m_fType = attrList["Function"].data<std::string>();
60
61 if (attrList["MuToLumi"].isNull()) return false;
62 m_muToLumi = attrList["MuToLumi"].data<float>();
63
64 if (attrList["Parameters"].isNull()) return false;
65 const coral::Blob& blob = attrList["Parameters"].data<coral::Blob>();
66
67
68 // Verify length
69 if ( static_cast<uint32_t>( blob.size() ) != 4*m_nPar) return false;
70
71 // Read calibration parameters
72 const float* p = static_cast<const float*>(blob.startingAddress());
73 for (unsigned int i=0; i<m_nPar; i++, p++)
74 m_parVec.push_back(*p);
75
76 return true;
77}
CaloCondBlobAlgs_fillNoiseFromASCII.blob
blob
Definition CaloCondBlobAlgs_fillNoiseFromASCII.py:95
xAOD::uint32_t
setEventNumber uint32_t
Definition EventInfo_v1.cxx:127

Member Data Documentation

◆ m_fType

std::string OnlineLumiCalibrator::m_fType
private

Definition at line 51 of file OnlineLumiCalibrator.h.

◆ m_muToLumi

float OnlineLumiCalibrator::m_muToLumi
private

Definition at line 52 of file OnlineLumiCalibrator.h.

◆ m_nPar

unsigned int OnlineLumiCalibrator::m_nPar
private

Definition at line 50 of file OnlineLumiCalibrator.h.

◆ m_parVec

std::vector<float> OnlineLumiCalibrator::m_parVec
private

Definition at line 53 of file OnlineLumiCalibrator.h.

The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.17.0