|
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
|
/*
* See header file for a description of this class.
*
* \author N. Amapane, G. Cerminara
*/
#include "CondFormats/DTObjects/interface/DTRecoConditions.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <iostream>
#include <algorithm>
#include <TFormula.h>
using std::cout;
using std::endl;
using std::map;
using std::string;
using std::vector;
DTRecoConditions::DTRecoConditions() : formula(nullptr), formulaType(0), expression("[0]"), theVersion(0) {}
DTRecoConditions::DTRecoConditions(const DTRecoConditions& iOther)
: formula(nullptr),
formulaType(0),
expression(iOther.expression),
payload(iOther.payload),
theVersion(iOther.theVersion) {}
const DTRecoConditions& DTRecoConditions::operator=(const DTRecoConditions& iOther) {
delete formula.load();
formula = nullptr;
formulaType = 0;
expression = iOther.expression;
payload = iOther.payload;
theVersion = iOther.theVersion;
return *this;
}
DTRecoConditions::~DTRecoConditions() { delete formula.load(); }
float DTRecoConditions::get(const DTWireId& wireid, double* x) const {
map<uint32_t, vector<double> >::const_iterator slIt = payload.find(wireid.superlayerId().rawId());
if (slIt == payload.end()) {
throw cms::Exception("InvalidInput") << "The SLId: " << wireid.superlayerId() << " is not in the paylaod map";
}
const vector<double>& par = slIt->second;
// Initialize if this is the first call
if (formulaType == 0) {
if (expression == "[0]") {
formulaType = 1;
} else if (expression == "par[step]") {
formulaType = 2;
} else {
std::unique_ptr<TFormula> temp{new TFormula("DTExpr", expression.c_str())};
TFormula* expected = nullptr;
if (formula.compare_exchange_strong(expected, temp.get())) {
//This thread set the value
temp.release();
}
formulaType = 99;
}
}
if (formulaType == 1 || par.size() == 1) {
// Return value is simply a constant. Assume this is the case also if only one parameter exists.
return par[0];
} else if (formulaType == 2) {
// Special case: par[i] represents the value for step i
return par[unsigned(x[0])];
} else {
// Use the formula corresponding to expression.
return (*formula).EvalPar(x, par.data());
}
}
void DTRecoConditions::set(const DTWireId& wireid, const std::vector<double>& values) {
payload[wireid.superlayerId()] = values;
}
DTRecoConditions::const_iterator DTRecoConditions::begin() const { return payload.begin(); }
DTRecoConditions::const_iterator DTRecoConditions::end() const { return payload.end(); }
|