|
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
|
// -*- C++ -*-
//
// Package: L1TObjects
// Class : L1CaloEtScale
//
// Implementation:
// <Notes on implementation>
//
// Author:
// Created: Wed Sep 27 17:18:27 CEST 2006
// $Id:
#include "CondFormats/L1TObjects/interface/L1CaloEtScale.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <stdexcept>
using std::endl;
using std::ostream;
using std::vector;
// default constructor (testing only!)
L1CaloEtScale::L1CaloEtScale() : m_linScaleMax(0x3ff), m_rankScaleMax(0x3f), m_linearLsb(1.0) {
for (unsigned i = 0; i < m_rankScaleMax; i++) {
m_thresholds.push_back(m_linearLsb * i);
}
}
// ctor that provides backwards compatibility with fixed max scale values
// OK to use this with e/gamma and rank scales
L1CaloEtScale::L1CaloEtScale(const double linearLsbInGeV, const vector<double>& thresholdsInGeV)
: m_linScaleMax(0x3ff), m_rankScaleMax(0x3f), m_linearLsb(linearLsbInGeV), m_thresholds(thresholdsInGeV) {
// protect against too many thresholds!
// while ( m_threshold.size() > (L1GctJetScale::maxRank+1) ) {
// m_thresholds.pop_back();
// }
}
// ctor that sets scale max values
L1CaloEtScale::L1CaloEtScale(const unsigned linScaleMax,
const unsigned rankScaleMax,
const double linearLsbInGeV,
const vector<double>& thresholdsInGeV)
: m_linScaleMax(linScaleMax),
m_rankScaleMax(rankScaleMax),
m_linearLsb(linearLsbInGeV),
m_thresholds(thresholdsInGeV) {}
L1CaloEtScale::~L1CaloEtScale() {}
// convert from linear Et to rank
uint16_t L1CaloEtScale::rank(const uint16_t linear) const { return rank((linear & m_linScaleMax) * m_linearLsb); }
/// convert from Et in GeV to rank
uint16_t L1CaloEtScale::rank(const double EtInGeV) const {
uint16_t out = 0;
for (unsigned i = 0; i < m_thresholds.size() && i < (unsigned)(m_rankScaleMax + 1); i++) {
if (EtInGeV >= m_thresholds.at(i)) {
out = i;
}
}
return out & m_rankScaleMax;
}
// convert from rank to Et/GeV
double L1CaloEtScale::et(const uint16_t rank) const {
// return bin centre, except for highest bin
// if (rank < m_thresholds.size()-1) {
// return (m_thresholds[rank+1]+m_thresholds[rank]) / 2;
// }
// else {
// return m_thresholds.back();
// }
// return bin lower edge
try {
return m_thresholds.at(rank);
} catch (std::out_of_range const&) {
throw cms::Exception("OutOfRange") << "Index out of range in L1CaloEtScale::et(rank)" << std::endl;
}
}
void L1CaloEtScale::print(ostream& s) const {
s << "L1CaloEtScale :" << endl;
s << " Input scale max = " << m_linScaleMax << endl;
s << " Input LSB = " << m_linearLsb << " GeV" << endl;
s << " Rank scale max = " << m_rankScaleMax << endl;
for (unsigned i = 0; i < m_thresholds.size(); i++) {
s << " Threshold " << i << " = " << m_thresholds[i] << " GeV" << endl;
}
}
std::ostream& operator<<(std::ostream& os, const L1CaloEtScale obj) {
os << "L1CaloEtScale :" << endl;
os << " Input scale max = " << obj.linScaleMax() << endl;
os << " Input LSB = " << obj.linearLsb() << " GeV" << endl;
os << " Rank scale max = " << obj.rankScaleMax() << endl;
for (unsigned i = 0; i < obj.getThresholds().size(); i++) {
os << " Threshold " << i << " = " << obj.getThresholds().at(i) << " GeV" << endl;
}
return os;
}
|