|
|
|
#include <L1GetHistLimits.h>
Classes | |
| struct | L1HistLimits |
Public Member Functions | |
| L1GetHistLimits (const edm::EventSetup &evSetup) | |
| const std::vector< float > & | l1HistBinThresholds (const L1GtObject &l1GtObject, const std::string &quantity) |
| const L1HistLimits & | l1HistLimits (const L1GtObject &l1GtObject, const std::string &quantity, const double histMinValue, const double histMaxValue) |
| const L1HistLimits & | l1HistLimits (const L1GtObject &l1GtObject, const std::string &quantity) |
| const double | l1HistLowerBinValue (const L1GtObject &l1GtObject, const std::string &quantity) |
| const int | l1HistNrBins (const L1GtObject &l1GtObject, const std::string &quantity) |
| const double | l1HistUpperBinValue (const L1GtObject &l1GtObject, const std::string &quantity) |
| virtual | ~L1GetHistLimits () |
Private Member Functions | |
| void | getHistLimits (const L1GtObject &l1GtObject, const std::string &quantity) |
Private Attributes | |
| const edm::EventSetup & | m_evSetup |
| L1HistLimits | m_l1HistLimits |
| all limits for a histogram | |
Description: use L1 scales to define histogram limits for L1 trigger objects.
Implementation: <TODO: enter implementation details>
$Date$ $Revision$
Definition at line 54 of file L1GetHistLimits.h.
| L1GetHistLimits::L1GetHistLimits | ( | const edm::EventSetup & | evSetup | ) | [explicit] |
Definition at line 38 of file L1GetHistLimits.cc.
:
m_evSetup(evSetup) {
//
}
| L1GetHistLimits::~L1GetHistLimits | ( | ) | [virtual] |
Definition at line 46 of file L1GetHistLimits.cc.
{
// empty
}
| void L1GetHistLimits::getHistLimits | ( | const L1GtObject & | l1GtObject, |
| const std::string & | quantity | ||
| ) | [private] |
for a L1 trigger object and a given quantity, compute the number of bins, the lower limit of the first bin, the upper limit of the last bin and the vector of bin thresholds
Definition at line 52 of file L1GetHistLimits.cc.
References L1GetHistLimits::L1HistLimits::binThresholds, BPTX, CenJet, L1CaloGeometry::emJetPhiBinHighEdge(), L1CaloGeometry::emJetPhiBinLowEdge(), ETM, L1CaloGeometry::etSumPhiBinHighEdge(), L1CaloGeometry::etSumPhiBinLowEdge(), ETT, ForJet, edm::EventSetup::get(), L1CaloGeometry::globalEtaBinHighEdge(), L1CaloGeometry::globalEtaBinLowEdge(), HfBitCounts, HfRingEtSums, HTM, L1CaloGeometry::htSumPhiBinHighEdge(), L1CaloGeometry::htSumPhiBinLowEdge(), HTT, IsoEG, j, JetCounts, L1GctEtHad::kEtHadMaxValue, L1GctEtMiss::kEtMissMaxValue, L1GctEtTotal::kEtTotalMaxValue, LogDebug, L1GetHistLimits::L1HistLimits::lowerBinValue, m_evSetup, m_l1HistLimits, Mu, NoIsoEG, L1GetHistLimits::L1HistLimits::nrBins, L1CaloGeometry::numberGctCentralEtaBinsPerHalf(), L1CaloGeometry::numberGctEmJetPhiBins(), L1CaloGeometry::numberGctEtSumPhiBins(), L1CaloGeometry::numberGctForwardEtaBinsPerHalf(), L1CaloGeometry::numberGctHtSumPhiBins(), edm::ESHandle< T >::product(), R3BINS, R3MAX, R3MIN, rad2deg(), TauJet, TechTrig, and L1GetHistLimits::L1HistLimits::upperBinValue.
Referenced by l1HistBinThresholds(), l1HistLimits(), l1HistLowerBinValue(), l1HistNrBins(), and l1HistUpperBinValue().
{
m_l1HistLimits.nrBins = 0;
m_l1HistLimits.lowerBinValue = 0;
m_l1HistLimits.upperBinValue = 0;
m_l1HistLimits.binThresholds.clear();
// number of objects is independent of the object type
// and hardcoded in the actual version
if (quantity == "NrObjects") {
m_l1HistLimits.nrBins = 15;
m_l1HistLimits.lowerBinValue = -0.5;
m_l1HistLimits.upperBinValue = 14.5;
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = m_l1HistLimits.lowerBinValue
+ iBin * (m_l1HistLimits.upperBinValue
- m_l1HistLimits.lowerBinValue)
/ m_l1HistLimits.nrBins;
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
return;
}
switch (l1GtObject) {
case Mu: {
if (quantity == "PT") {
edm::ESHandle<L1MuTriggerPtScale> muPtScale;
m_evSetup.get<L1MuTriggerPtScaleRcd>().get(muPtScale);
m_l1HistLimits.nrBins = muPtScale->getPtScale()->getNBins();
m_l1HistLimits.lowerBinValue
= muPtScale->getPtScale()->getScaleMin();
m_l1HistLimits.upperBinValue
= muPtScale->getPtScale()->getScaleMax();
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= muPtScale->getPtScale()->getValue(iBin);
}
// last limit for muon is set too high (10^6) - resize the last bin
float lastBinSize = m_l1HistLimits.upperBinValue
- m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins - 1];
// limit the size of the last bin to maximum 200 GeV
float maxLaxtBinsize = 200;
if (lastBinSize >= maxLaxtBinsize) {
m_l1HistLimits.upperBinValue
= m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins - 1]
+ 2.* (m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins - 1]
- m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins - 2]);
LogDebug("L1GetHistLimits")
<< "\n L1ExtraMuon: PT histogram"
<< "\nm_l1HistLimits.upperBinValue truncated to = "
<< m_l1HistLimits.upperBinValue << std::endl;
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
edm::ESHandle<L1MuTriggerScales> muScales;
m_evSetup.get<L1MuTriggerScalesRcd>().get(muScales);
if (quantity == "eta") {
// eta scale defined for positive values - need to be symmetrized
int histNrBinsHalf = muScales->getGMTEtaScale()->getNBins();
m_l1HistLimits.lowerBinValue
= muScales->getGMTEtaScale()->getScaleMin();
m_l1HistLimits.upperBinValue
= muScales->getGMTEtaScale()->getScaleMax();
m_l1HistLimits.nrBins = 2 * histNrBinsHalf;
m_l1HistLimits.lowerBinValue
= -m_l1HistLimits.upperBinValue;
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
int iBin = 0;
for (int j = histNrBinsHalf; j > 0; j--, iBin++) {
m_l1HistLimits.binThresholds[iBin] = (-1)
* muScales->getGMTEtaScale()->getValue(j);
}
for (int j = 0; j <= histNrBinsHalf; j++, iBin++) {
m_l1HistLimits.binThresholds[iBin]
= muScales->getGMTEtaScale()->getValue(j);
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else {
m_l1HistLimits.nrBins = muScales->getPhiScale()->getNBins();
m_l1HistLimits.lowerBinValue = rad2deg(
muScales->getPhiScale()->getScaleMin());
m_l1HistLimits.upperBinValue = rad2deg(
muScales->getPhiScale()->getScaleMax());
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins
+ 1);
for (int iBin = 0; iBin <= m_l1HistLimits.nrBins; iBin++) {
m_l1HistLimits.binThresholds[iBin] = rad2deg(
muScales->getPhiScale()->getValue(iBin));
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
}
}
}
break;
case NoIsoEG:
case IsoEG: {
// common scales for NoIsoEG and IsoEG
if (quantity == "ET") {
edm::ESHandle<L1CaloEtScale> emScale;
m_evSetup.get<L1EmEtScaleRcd>().get(emScale);
std::vector<double> emThresholds = emScale->getThresholds();
m_l1HistLimits.nrBins = emThresholds.size();
m_l1HistLimits.lowerBinValue = emThresholds.at(0);
// FIXME high edge retrieval in the scale definition
// now, last bin has the same width like the last but one
m_l1HistLimits.upperBinValue
= emThresholds[m_l1HistLimits.nrBins - 1]
+ (emThresholds[m_l1HistLimits.nrBins - 1]
- emThresholds[m_l1HistLimits.nrBins - 2]);
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= static_cast<float>(emThresholds[iBin]);
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
edm::ESHandle<L1CaloGeometry> caloGeomESH;
m_evSetup.get<L1CaloGeometryRecord>().get(caloGeomESH);
const L1CaloGeometry* caloGeomScales = caloGeomESH.product();
if (quantity == "eta") {
m_l1HistLimits.nrBins
= 2 * (caloGeomScales->numberGctCentralEtaBinsPerHalf()
+ caloGeomScales->numberGctForwardEtaBinsPerHalf());
m_l1HistLimits.lowerBinValue
= caloGeomScales->globalEtaBinLowEdge(0);
m_l1HistLimits.upperBinValue
= caloGeomScales->globalEtaBinHighEdge(
m_l1HistLimits.nrBins - 1);
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= caloGeomScales->globalEtaBinLowEdge(iBin);
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else {
m_l1HistLimits.nrBins
= caloGeomScales->numberGctEmJetPhiBins();
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins
+ 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = rad2deg(
caloGeomScales->emJetPhiBinLowEdge(iBin));
// treat correctly the 10 deg anti-clockwise rotation
if (rad2deg(caloGeomScales->emJetPhiBinHighEdge(iBin))
< m_l1HistLimits.binThresholds[iBin]) {
m_l1HistLimits.binThresholds[iBin]
= m_l1HistLimits.binThresholds[iBin] - 360.;
}
}
m_l1HistLimits.lowerBinValue
= m_l1HistLimits.binThresholds[0];
// last bin upper limit
m_l1HistLimits.upperBinValue = rad2deg(
caloGeomScales->emJetPhiBinHighEdge(
m_l1HistLimits.nrBins - 1));
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
}
}
}
break;
case CenJet:
case ForJet:
case TauJet: {
// common scales for all jets
if (quantity == "ET") {
edm::ESHandle<L1CaloEtScale> jetScale;
m_evSetup.get<L1JetEtScaleRcd>().get(jetScale);
std::vector<double> jetThresholds = jetScale->getThresholds();
m_l1HistLimits.nrBins = jetThresholds.size();
m_l1HistLimits.lowerBinValue = jetThresholds.at(0);
// FIXME high edge retrieval in the scale definition
// now, last bin has the same width like the last but one
m_l1HistLimits.upperBinValue
= jetThresholds[m_l1HistLimits.nrBins - 1]
+ (jetThresholds[m_l1HistLimits.nrBins - 1]
- jetThresholds[m_l1HistLimits.nrBins - 2]);
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= static_cast<float>(jetThresholds[iBin]);
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
edm::ESHandle<L1CaloGeometry> caloGeomESH;
m_evSetup.get<L1CaloGeometryRecord>().get(caloGeomESH);
const L1CaloGeometry* caloGeomScales = caloGeomESH.product();
if (quantity == "eta") {
m_l1HistLimits.nrBins
= 2 * (caloGeomScales->numberGctCentralEtaBinsPerHalf()
+ caloGeomScales->numberGctForwardEtaBinsPerHalf());
m_l1HistLimits.lowerBinValue
= caloGeomScales->globalEtaBinLowEdge(0);
m_l1HistLimits.upperBinValue
= caloGeomScales->globalEtaBinHighEdge(
m_l1HistLimits.nrBins - 1);
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= caloGeomScales->globalEtaBinLowEdge(iBin);
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else {
m_l1HistLimits.nrBins
= caloGeomScales->numberGctEmJetPhiBins();
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins
+ 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = rad2deg(
caloGeomScales->emJetPhiBinLowEdge(iBin));
// treat correctly the 10 deg anti-clockwise rotation
if (rad2deg(caloGeomScales->emJetPhiBinHighEdge(iBin))
< m_l1HistLimits.binThresholds[iBin]) {
m_l1HistLimits.binThresholds[iBin]
= m_l1HistLimits.binThresholds[iBin] - 360.;
}
}
m_l1HistLimits.lowerBinValue
= m_l1HistLimits.binThresholds[0];
// last bin upper limit
m_l1HistLimits.upperBinValue = rad2deg(
caloGeomScales->emJetPhiBinHighEdge(
m_l1HistLimits.nrBins - 1));
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
}
}
}
break;
case ETM: {
if (quantity == "ET") {
edm::ESHandle<L1CaloEtScale> etMissScale;
m_evSetup.get<L1JetEtScaleRcd>().get(etMissScale);
const double etSumLSB = etMissScale->linearLsb() ;
m_l1HistLimits.nrBins = L1GctEtMiss::kEtMissMaxValue;
m_l1HistLimits.lowerBinValue = 0;
m_l1HistLimits.upperBinValue = (m_l1HistLimits.nrBins + 1)
* etSumLSB;
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = iBin * etSumLSB;
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
edm::ESHandle<L1CaloGeometry> caloGeomESH;
m_evSetup.get<L1CaloGeometryRecord>().get(caloGeomESH);
const L1CaloGeometry* caloGeomScales = caloGeomESH.product();
if (quantity == "eta") {
// do nothing, eta is not defined for ETM
} else {
m_l1HistLimits.nrBins
= caloGeomScales->numberGctEtSumPhiBins();
m_l1HistLimits.lowerBinValue = rad2deg(
caloGeomScales->etSumPhiBinLowEdge(0));
m_l1HistLimits.upperBinValue = rad2deg(
caloGeomScales->etSumPhiBinHighEdge(
m_l1HistLimits.nrBins - 1));
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins
+ 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = rad2deg(
caloGeomScales->etSumPhiBinLowEdge(iBin));
}
// last bin upper limit
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
}
}
}
break;
case ETT: {
if (quantity == "ET") {
edm::ESHandle<L1CaloEtScale> etMissScale;
m_evSetup.get<L1JetEtScaleRcd>().get(etMissScale);
const double etSumLSB = etMissScale->linearLsb() ;
m_l1HistLimits.nrBins = L1GctEtTotal::kEtTotalMaxValue;
m_l1HistLimits.lowerBinValue = 0;
m_l1HistLimits.upperBinValue = (m_l1HistLimits.nrBins + 1)
* etSumLSB;
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = iBin * etSumLSB;
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
// do nothing, eta and phi are not defined for ETT
}
}
break;
case HTT: {
if (quantity == "ET") {
edm::ESHandle< L1GctJetFinderParams > jetFinderParams ;
m_evSetup.get< L1GctJetFinderParamsRcd >().get( jetFinderParams ) ;
double htSumLSB = jetFinderParams->getHtLsbGeV();
m_l1HistLimits.nrBins = L1GctEtHad::kEtHadMaxValue;
m_l1HistLimits.lowerBinValue = 0;
m_l1HistLimits.upperBinValue = (m_l1HistLimits.nrBins + 1)
* htSumLSB;
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = iBin * htSumLSB;
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
// do nothing, eta and phi are not defined for HTT
}
}
break;
case HTM: {
if (quantity == "ET") {
edm::ESHandle<L1CaloEtScale> htMissScale;
m_evSetup.get<L1HtMissScaleRcd>().get(htMissScale);
const std::vector<double>& htThresholds =
htMissScale->getThresholds();
m_l1HistLimits.nrBins = htThresholds.size();
m_l1HistLimits.lowerBinValue = htThresholds[0];
// FIXME high edge retrieval in the scale definition
// now, last bin has the same width like the last but one
m_l1HistLimits.upperBinValue
= htThresholds[m_l1HistLimits.nrBins - 1]
+ (htThresholds[m_l1HistLimits.nrBins - 1]
- htThresholds[m_l1HistLimits.nrBins - 2]);
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= static_cast<float>(htThresholds[iBin]);
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
edm::ESHandle<L1CaloGeometry> caloGeomESH;
m_evSetup.get<L1CaloGeometryRecord>().get(caloGeomESH);
const L1CaloGeometry* caloGeomScales = caloGeomESH.product();
if (quantity == "eta") {
// do nothing, eta is not defined for HTM
} else {
m_l1HistLimits.nrBins
= caloGeomScales->numberGctHtSumPhiBins();
m_l1HistLimits.lowerBinValue = rad2deg(
caloGeomScales->htSumPhiBinLowEdge(0));
m_l1HistLimits.upperBinValue = rad2deg(
caloGeomScales->htSumPhiBinHighEdge(
m_l1HistLimits.nrBins - 1));
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins
+ 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin] = rad2deg(
caloGeomScales->htSumPhiBinLowEdge(iBin));
}
// last bin upper limit
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
}
}
}
break;
case JetCounts: {
}
break;
case HfBitCounts: {
// there are no scales for HfBitCounts, so one implements a fixed scale
// use same values as GCT for 3 bits
const unsigned int R3BINS = 8;
const float R3MIN = -0.5;
const float R3MAX = 7.5;
m_l1HistLimits.nrBins
= R3BINS;
m_l1HistLimits.lowerBinValue
= R3MIN;
m_l1HistLimits.upperBinValue
= R3MAX;
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins
+ 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= R3MIN + iBin*(R3MAX - R3MIN)/R3BINS;
}
// last bin upper limit
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
}
break;
case HfRingEtSums: {
if (quantity == "ET") {
edm::ESHandle<L1CaloEtScale> hfRingEtScale;
m_evSetup.get<L1HfRingEtScaleRcd>().get(hfRingEtScale);
const std::vector<double>& hfRingEtThresholds =
hfRingEtScale->getThresholds();
m_l1HistLimits.nrBins = hfRingEtThresholds.size();
m_l1HistLimits.lowerBinValue = hfRingEtThresholds[0];
// FIXME high edge retrieval in the scale definition
// now, last bin has the same width like the last but one
m_l1HistLimits.upperBinValue
= hfRingEtThresholds[m_l1HistLimits.nrBins - 1]
+ (hfRingEtThresholds[m_l1HistLimits.nrBins - 1]
- hfRingEtThresholds[m_l1HistLimits.nrBins - 2]);
m_l1HistLimits.binThresholds.resize(m_l1HistLimits.nrBins + 1);
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
m_l1HistLimits.binThresholds[iBin]
= static_cast<float>(hfRingEtThresholds[iBin]);
}
// set last bin upper edge
m_l1HistLimits.binThresholds[m_l1HistLimits.nrBins]
= m_l1HistLimits.upperBinValue;
} else if (quantity == "eta" || quantity == "phi") {
// do nothing, eta and phi are not defined for HfRingEtSums
}
}
break;
case TechTrig:
case Castor:
case BPTX:
default: {
// do nothing, for these cases ET/PT, eta and phi are not defined
}
break;
}
}
| const std::vector< float > & L1GetHistLimits::l1HistBinThresholds | ( | const L1GtObject & | l1GtObject, |
| const std::string & | quantity | ||
| ) |
for a L1 trigger object and a given quantity, return the vector of bin thresholds
Definition at line 816 of file L1GetHistLimits.cc.
References L1GetHistLimits::L1HistLimits::binThresholds, getHistLimits(), and m_l1HistLimits.
{
getHistLimits(l1GtObject, quantity);
return m_l1HistLimits.binThresholds;
}
| const L1GetHistLimits::L1HistLimits & L1GetHistLimits::l1HistLimits | ( | const L1GtObject & | l1GtObject, |
| const std::string & | quantity, | ||
| const double | histMinValue, | ||
| const double | histMaxValue | ||
| ) |
for a L1 trigger object and a given quantity, return the real limits for a histogram given an arbitrary range
Definition at line 691 of file L1GetHistLimits.cc.
References L1GetHistLimits::L1HistLimits::binThresholds, getHistLimits(), edm::isDebugEnabled(), LogDebug, LogTrace, L1GetHistLimits::L1HistLimits::lowerBinValue, m_l1HistLimits, L1GetHistLimits::L1HistLimits::nrBins, and L1GetHistLimits::L1HistLimits::upperBinValue.
{
getHistLimits(l1GtObject, quantity);
bool foundLowerBinValue = false;
bool foundUpperBinValue = false;
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
if (m_l1HistLimits.binThresholds[iBin] <= histMinValue) {
m_l1HistLimits.lowerBinValue = m_l1HistLimits.binThresholds[iBin];
foundLowerBinValue = true;
break;
}
}
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
if (m_l1HistLimits.binThresholds[iBin] > histMaxValue) {
m_l1HistLimits.upperBinValue = m_l1HistLimits.binThresholds[iBin];
foundUpperBinValue = true;
break;
}
}
if (foundLowerBinValue && foundUpperBinValue) {
int countBins = -1;
std::vector<float> binThresh;
binThresh.reserve(m_l1HistLimits.binThresholds.size());
for (int iBin = 0; iBin < m_l1HistLimits.nrBins; ++iBin) {
if ((m_l1HistLimits.binThresholds[iBin] >= histMinValue)
&& m_l1HistLimits.binThresholds[iBin] < histMaxValue) {
m_l1HistLimits.upperBinValue
= m_l1HistLimits.binThresholds[iBin];
countBins++;
binThresh.push_back(m_l1HistLimits.binThresholds[iBin]);
}
}
m_l1HistLimits.nrBins = countBins;
m_l1HistLimits.binThresholds.clear();
m_l1HistLimits.binThresholds = binThresh;
// FIXME last bin untested.
} else {
m_l1HistLimits.nrBins = 0;
m_l1HistLimits.lowerBinValue = 0;
m_l1HistLimits.upperBinValue = 0;
m_l1HistLimits.binThresholds.clear();
LogDebug("L1GetHistLimits") << "\n Histogram limits for L1GtObject"
<< l1GtObject << " and quantity " << quantity
<< " within the required range [" << histMinValue << ", "
<< histMaxValue << "] not found."
<< "\n The range is not included in the original histogram range."
<< std::endl;
return m_l1HistLimits;
}
if (edm::isDebugEnabled()) {
LogDebug("L1GetHistLimits") << "\n Histogram limits for L1GtObject"
<< l1GtObject << " and quantity " << quantity
<< "\n Number of bins: " << m_l1HistLimits.nrBins
<< "\n Lower limit of first bin: "
<< m_l1HistLimits.lowerBinValue
<< "\n Upper limit of last bin: "
<< m_l1HistLimits.upperBinValue << std::endl;
int binThreshSize = static_cast<int>(m_l1HistLimits.binThresholds.size());
if (binThreshSize != (m_l1HistLimits.nrBins + 1)) {
LogTrace("L1GetHistLimits")
<< "\n Warning: inconsistent nrBins and binThresholds size"
<< "\n Number of bins nrBins = " << m_l1HistLimits.nrBins
<< "\n binThresholds size = " << binThreshSize
<< "\n Please fix the L1GetLimits class.\n\n" << std::endl;
}
for (int iBin = 0; iBin < binThreshSize; ++iBin) {
LogTrace("L1GetHistLimits") << " Bin " << std::right
<< std::setw(5) << iBin << ": "
<< m_l1HistLimits.binThresholds[iBin] << std::endl;
}
}
return m_l1HistLimits;
}
| const L1GetHistLimits::L1HistLimits & L1GetHistLimits::l1HistLimits | ( | const L1GtObject & | l1GtObject, |
| const std::string & | quantity | ||
| ) |
for a L1 trigger object and a given quantity, return all limits for a histogram
Definition at line 652 of file L1GetHistLimits.cc.
References L1GetHistLimits::L1HistLimits::binThresholds, getHistLimits(), edm::isDebugEnabled(), LogDebug, LogTrace, L1GetHistLimits::L1HistLimits::lowerBinValue, m_l1HistLimits, L1GetHistLimits::L1HistLimits::nrBins, and L1GetHistLimits::L1HistLimits::upperBinValue.
Referenced by L1ExtraDQM::L1ExtraMonElement< CollectionType >::bookHistograms().
{
getHistLimits(l1GtObject, quantity);
if (edm::isDebugEnabled()) {
LogDebug("L1GetHistLimits") << "\n Histogram limits for L1GtObject "
<< l1GtObject << " and quantity " << quantity
<< "\n Number of bins: " << m_l1HistLimits.nrBins
<< "\n Lower limit of first bin: "
<< m_l1HistLimits.lowerBinValue
<< "\n Upper limit of last bin: "
<< m_l1HistLimits.upperBinValue << std::endl;
int binThreshSize = static_cast<int>(m_l1HistLimits.binThresholds.size());
if (binThreshSize != (m_l1HistLimits.nrBins + 1)) {
LogTrace("L1GetHistLimits")
<< "\n Warning: inconsistent nrBins and binThresholds size"
<< "\n Number of bins nrBins = " << m_l1HistLimits.nrBins
<< "\n binThresholds size = " << binThreshSize
<< "\n Please fix the L1GetLimits class.\n\n" << std::endl;
}
for (int iBin = 0; iBin < binThreshSize; ++iBin) {
LogTrace("L1GetHistLimits") << " Bin " << std::right
<< std::setw(5) << iBin << ": "
<< m_l1HistLimits.binThresholds[iBin] << std::endl;
}
}
return m_l1HistLimits;
}
| const double L1GetHistLimits::l1HistLowerBinValue | ( | const L1GtObject & | l1GtObject, |
| const std::string & | quantity | ||
| ) |
for a L1 trigger object and a given quantity, return the lower limit of the first bin in the histogram
Definition at line 800 of file L1GetHistLimits.cc.
References getHistLimits(), L1GetHistLimits::L1HistLimits::lowerBinValue, and m_l1HistLimits.
{
getHistLimits(l1GtObject, quantity);
return m_l1HistLimits.lowerBinValue;
}
| const int L1GetHistLimits::l1HistNrBins | ( | const L1GtObject & | l1GtObject, |
| const std::string & | quantity | ||
| ) |
for a L1 trigger object and a given quantity, return the numbers of bins for a given histogram
Definition at line 792 of file L1GetHistLimits.cc.
References getHistLimits(), m_l1HistLimits, and L1GetHistLimits::L1HistLimits::nrBins.
{
getHistLimits(l1GtObject, quantity);
return m_l1HistLimits.nrBins;
}
| const double L1GetHistLimits::l1HistUpperBinValue | ( | const L1GtObject & | l1GtObject, |
| const std::string & | quantity | ||
| ) |
for a L1 trigger object and a given quantity, return the upper limit of the last bin in the histogram
Definition at line 808 of file L1GetHistLimits.cc.
References getHistLimits(), m_l1HistLimits, and L1GetHistLimits::L1HistLimits::upperBinValue.
{
getHistLimits(l1GtObject, quantity);
return m_l1HistLimits.upperBinValue;
}
const edm::EventSetup& L1GetHistLimits::m_evSetup [private] |
Definition at line 121 of file L1GetHistLimits.h.
Referenced by getHistLimits().
L1HistLimits L1GetHistLimits::m_l1HistLimits [private] |
all limits for a histogram
Definition at line 124 of file L1GetHistLimits.h.
Referenced by getHistLimits(), l1HistBinThresholds(), l1HistLimits(), l1HistLowerBinValue(), l1HistNrBins(), and l1HistUpperBinValue().
1.7.3