|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
#ifndef CondEx_Efficiency_H
#define CondEx_Efficiency_H
/* example of polymorphic condition
* LUT, function, mixed....
* this is just a prototype: classes do not need to be defined and declared in the same file
* at the moment though all derived classes better sit in the same package together with the base one
*/
#include "CondFormats/Serialization/interface/Serializable.h"
#include <cmath>
#include <iostream>
namespace condex {
/* very simple base class
* trivial inheritance, no template tricks
*/
class Efficiency {
public:
Efficiency() {}
virtual ~Efficiency() {}
virtual void initialize() { std::cout << "initializing base class Efficiency" << std::endl; }
float operator()(float pt, float eta) const { return value(pt, eta); }
virtual float value(float pt, float eta) const = 0;
COND_SERIALIZABLE;
};
class ParametricEfficiencyInPt : public Efficiency {
public:
ParametricEfficiencyInPt() : cutLow(0), cutHigh(0), low(0), high(0) {}
ParametricEfficiencyInPt(float cm, float ch, float el, float eh) : cutLow(cm), cutHigh(ch), low(el), high(eh) {}
private:
float value(float pt, float) const override {
if (pt < low)
return cutLow;
if (pt > high)
return cutHigh;
return cutLow + (pt - low) / (high - low) * (cutHigh - cutLow);
}
float cutLow, cutHigh;
float low, high;
COND_SERIALIZABLE;
};
class ParametricEfficiencyInEta : public Efficiency {
public:
ParametricEfficiencyInEta() : cutLow(0), cutHigh(0), low(0), high(0) {}
ParametricEfficiencyInEta(float cmin, float cmax, float el, float eh)
: cutLow(cmin), cutHigh(cmax), low(el), high(eh) {}
private:
float value(float, float eta) const override {
eta = std::abs(eta);
if (eta < low)
return cutLow;
if (eta > high)
return cutHigh;
return cutLow + (eta - low) / (high - low) * (cutHigh - cutLow);
}
float cutLow, cutHigh;
float low, high;
COND_SERIALIZABLE;
};
} // namespace condex
#endif
|