|
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
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
|
#include "DataFormats/Luminosity/interface/LumiSummary.h"
#include <iomanip>
#include <ostream>
#include <iostream>
float LumiSummary::avgInsDelLumi() const {
size_t iIndex = lumiversion_.rfind('v');
//i.e. not "-1" and not "DIP", "-1" and "DIP" lumi are already corrected and unit conversion included in the raw data.
if (iIndex != std::string::npos) {
return avginsdellumi_ * 1000.0;
}
return avginsdellumi_;
}
float LumiSummary::intgDelLumi() const { return this->avgInsDelLumi() * float(this->lumiSectionLength()); }
float LumiSummary::avgInsDelLumiErr() const { return avginsdellumierr_; }
float LumiSummary::intgRecLumi() const { return this->avgInsRecLumi() * float(this->lumiSectionLength()); }
short LumiSummary::lumiSecQual() const { return lumisecqual_; }
unsigned long long LumiSummary::deadcount() const { return deadcount_; }
unsigned long long LumiSummary::bitzerocount() const { return bitzerocount_; }
float LumiSummary::deadFrac() const {
//definition: deadcount/bitzerocount
//if no trigger data, return deadfraction 1.0,mask out this LS
//if bitzerocount=0, return -1.0 meaning no beam
if (lumiversion_ == "DIP") {
return float(deadcount_) / float(bitzerocount_);
}
if (l1data_.empty())
return 1.0;
if (bitzerocount_ == 0)
return -1.0;
return float(deadcount_) / float(bitzerocount_);
}
float LumiSummary::liveFrac() const {
//1-deadfraction
//else if deadfraction<0 meaning no beam, live fraction=0
//
if (deadFrac() < 0)
return 0;
return 1 - deadFrac();
}
float LumiSummary::lumiSectionLength() const {
//numorbits*3564*24.95e-09
return numorbit_ * 3564.0 * 24.95e-9;
}
unsigned int LumiSummary::lsNumber() const { return lsnumber_; }
unsigned int LumiSummary::startOrbit() const { return startorbit_; }
unsigned int LumiSummary::numOrbit() const { return numorbit_; }
bool LumiSummary::isValid() const { return avginsdellumi_ > 0; }
LumiSummary::L1 LumiSummary::l1info(unsigned int idx) const { return l1data_.at(idx); }
LumiSummary::HLT LumiSummary::hltinfo(unsigned int idx) const { return hltdata_.at(idx); }
size_t LumiSummary::nTriggerLine() const { return l1data_.size(); }
size_t LumiSummary::nHLTPath() const { return hltdata_.size(); }
float LumiSummary::avgInsRecLumi() const { return this->avgInsDelLumi() * liveFrac(); }
float LumiSummary::avgInsRecLumiErr() const { return avginsdellumierr_ * liveFrac(); }
bool LumiSummary::isProductEqual(LumiSummary const& next) const {
return (avginsdellumi_ == next.avginsdellumi_ && avginsdellumierr_ == next.avginsdellumierr_ &&
lumisecqual_ == next.lumisecqual_ && deadcount_ == next.deadcount_ && lsnumber_ == next.lsnumber_ &&
startorbit_ == next.startorbit_ && numorbit_ == next.numorbit_ && l1data_.size() == next.l1data_.size() &&
hltdata_.size() == next.hltdata_.size() && lumiversion_ == next.lumiversion_);
}
std::string LumiSummary::lumiVersion() const { return lumiversion_; }
void LumiSummary::setLumiVersion(const std::string& lumiversion) { lumiversion_ = lumiversion; }
void LumiSummary::setLumiData(float instlumi, float instlumierr, short lumiquality) {
avginsdellumi_ = instlumi;
avginsdellumierr_ = instlumierr;
lumisecqual_ = lumiquality;
}
void LumiSummary::setDeadCount(unsigned long long deadcount) { deadcount_ = deadcount; }
void LumiSummary::setBitZeroCount(unsigned long long bitzerocount) { bitzerocount_ = bitzerocount; }
void LumiSummary::setlsnumber(unsigned int lsnumber) { lsnumber_ = lsnumber; }
void LumiSummary::setOrbitData(unsigned int startorbit, unsigned int numorbit) {
startorbit_ = startorbit;
numorbit_ = numorbit;
}
void LumiSummary::swapL1Data(std::vector<L1>& l1data) { l1data_.swap(l1data); }
void LumiSummary::swapHLTData(std::vector<HLT>& hltdata) { hltdata_.swap(hltdata); }
void LumiSummary::copyL1Data(const std::vector<L1>& l1data) { l1data_.assign(l1data.begin(), l1data.end()); }
void LumiSummary::copyHLTData(const std::vector<HLT>& hltdata) { hltdata_.assign(hltdata.begin(), hltdata.end()); }
std::ostream& operator<<(std::ostream& s, const LumiSummary& lumiSummary) {
s << "\nDumping LumiSummary\n\n";
if (!lumiSummary.isValid()) {
s << " === Invalid Lumi values === \n";
}
s << " lumiVersion = " << lumiSummary.lumiVersion() << "\n";
s << " avgInsDelLumi = " << lumiSummary.avgInsDelLumi() << "\n";
s << " avgIntgDelLumi = " << lumiSummary.intgDelLumi() << "\n";
s << " avgInsDelLumiErr = " << lumiSummary.avgInsDelLumiErr() << "\n";
s << " lumiSecQual = " << lumiSummary.lumiSecQual() << "\n";
s << " deadCount = " << lumiSummary.deadcount() << "\n";
s << " bitZeroCount = " << lumiSummary.bitzerocount() << "\n";
s << " deadFrac = " << (float)lumiSummary.deadFrac() << "\n";
s << " liveFrac = " << (float)lumiSummary.liveFrac() << "\n";
s << " lsNumber = " << lumiSummary.lsNumber() << "\n";
s << " startOrbit = " << lumiSummary.startOrbit() << "\n";
s << " numOrbit = " << lumiSummary.numOrbit() << "\n";
s << " avgInsRecLumi = " << lumiSummary.avgInsRecLumi() << "\n";
s << " avgInsRecLumiErr = " << lumiSummary.avgInsRecLumiErr() << "\n\n";
s << std::setw(15) << "l1nameidx";
s << std::setw(15) << "l1prescale";
s << "\n";
size_t nTriggers = lumiSummary.nTriggerLine();
size_t nHLTPath = lumiSummary.nHLTPath();
for (unsigned int i = 0; i < nTriggers; ++i) {
s << std::setw(15);
s << lumiSummary.l1info(i).triggernameidx;
s << std::setw(15);
s << lumiSummary.l1info(i).prescale;
s << "\n";
}
s << std::setw(15) << "hltpathidx";
s << std::setw(15) << "hltprescale";
s << "\n";
for (unsigned int i = 0; i < nHLTPath; ++i) {
s << std::setw(15);
s << lumiSummary.hltinfo(i).pathnameidx;
s << std::setw(15);
s << lumiSummary.hltinfo(i).prescale;
s << "\n";
}
return s << "\n";
}
|