|
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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
#include <algorithm>
#include <iostream>
#include <set>
#include "CondFormats/HcalObjects/interface/HcalSiPMCharacteristics.h"
#include "CondFormats/HcalObjects/interface/HcalObjectAddons.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
HcalSiPMCharacteristics::HcalSiPMCharacteristics(const HcalSiPMCharacteristicsAddons::Helper& helper)
: mPItems(helper.mPItems.begin(), helper.mPItems.end()) {
initialize();
}
HcalSiPMCharacteristics::~HcalSiPMCharacteristics() {}
// copy-ctor
HcalSiPMCharacteristics::HcalSiPMCharacteristics(const HcalSiPMCharacteristics& src)
: mPItems(src.mPItems), mPItemsByType(src.mPItemsByType) {}
// copy assignment operator
HcalSiPMCharacteristics& HcalSiPMCharacteristics::operator=(const HcalSiPMCharacteristics& rhs) {
HcalSiPMCharacteristics temp(rhs);
temp.swap(*this);
return *this;
}
// public swap function
void HcalSiPMCharacteristics::swap(HcalSiPMCharacteristics& other) {
std::swap(mPItems, other.mPItems);
std::swap(mPItemsByType, other.mPItemsByType);
}
// move constructor
HcalSiPMCharacteristics::HcalSiPMCharacteristics(HcalSiPMCharacteristics&& other) : HcalSiPMCharacteristics() {
other.swap(*this);
}
const HcalSiPMCharacteristics::PrecisionItem* HcalSiPMCharacteristics::findByType(int type) const {
const HcalSiPMCharacteristics::PrecisionItem* retItem = nullptr;
for (unsigned int i = 0; i < getTypes(); i++) {
auto iter = &mPItems.at(i);
if (type == iter->type_)
retItem = iter;
}
return retItem;
//NOT WORKING:
//PrecisionItem target(type, 0, 0, 0, 0, 0, 0, 0);
//return HcalObjectAddons::findByT<PrecisionItem,HcalSiPMCharacteristicsAddons::LessByType>(&target,mPItemsByType);
}
HcalSiPMCharacteristicsAddons::Helper::Helper() {}
bool HcalSiPMCharacteristicsAddons::Helper::loadObject(
int type, int pixels, float parLin1, float parLin2, float parLin3, float crossTalk, int auxi1, float auxi2) {
HcalSiPMCharacteristics::PrecisionItem target(type, pixels, parLin1, parLin2, parLin3, crossTalk, auxi1, auxi2);
auto iter = mPItems.find(target);
if (iter != mPItems.end()) {
edm::LogWarning("HCAL") << "HcalSiPMCharacteristics::loadObject type " << type << " already exists with pixels "
<< iter->pixels_ << " NoLinearity parameters " << iter->parLin1_ << ":" << iter->parLin2_
<< ":" << iter->parLin3_ << " CrossTalk parameter " << iter->crossTalk_ << " new values "
<< pixels << ", " << parLin1 << ", " << parLin2 << ", " << parLin3 << ", " << crossTalk
<< ", " << auxi1 << " and " << auxi2 << " are ignored";
return false;
} else {
mPItems.insert(target);
return true;
}
}
int HcalSiPMCharacteristics::getPixels(int type) const {
const HcalSiPMCharacteristics::PrecisionItem* item = findByType(type);
return (item ? item->pixels_ : 0);
}
std::vector<float> HcalSiPMCharacteristics::getNonLinearities(int type) const {
const HcalSiPMCharacteristics::PrecisionItem* item = findByType(type);
std::vector<float> pars;
pars.reserve(3);
if (item) {
pars.push_back(item->parLin1_);
pars.push_back(item->parLin2_);
pars.push_back(item->parLin3_);
}
return pars;
}
float HcalSiPMCharacteristics::getCrossTalk(int type) const {
const PrecisionItem* item = findByType(type);
return (item ? item->crossTalk_ : 0);
}
int HcalSiPMCharacteristics::getAuxi1(int type) const {
const HcalSiPMCharacteristics::PrecisionItem* item = findByType(type);
return (item ? item->auxi1_ : 0);
}
float HcalSiPMCharacteristics::getAuxi2(int type) const {
const HcalSiPMCharacteristics::PrecisionItem* item = findByType(type);
return (item ? item->auxi2_ : 0);
}
void HcalSiPMCharacteristics::sortByType() {
HcalObjectAddons::sortByT<PrecisionItem, HcalSiPMCharacteristicsAddons::LessByType>(mPItems, mPItemsByType);
}
void HcalSiPMCharacteristics::initialize() { HcalSiPMCharacteristics::sortByType(); }
|