CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
L1RCTParameters Class Reference

#include <CondFormats/L1TObjects/interface/L1RCTParameters.h>

Public Member Functions

unsigned short calcCard (unsigned short rct_iphi, unsigned short absIeta) const
 
unsigned short calcCrate (unsigned short rct_iphi, short ieta) const
 
unsigned short calcIAbsEta (unsigned short iCrate, unsigned short iCard, unsigned short iTower) const
 
short calcIEta (unsigned short iCrate, unsigned short iCard, unsigned short iTower) const
 
unsigned short calcIPhi (unsigned short iCrate, unsigned short iCard, unsigned short iTower) const
 
unsigned short calcTower (unsigned short rct_iphi, unsigned short absIeta) const
 
double eActivityCut () const
 
const std::vector< double > & eGammaECalScaleFactors () const
 
const std::vector< double > & eGammaHCalScaleFactors () const
 
double eGammaLSB () const
 
float EGammaTPGSum (const float &ecal, const float &hcal, const unsigned &iAbsEta) const
 
unsigned eicIsolationThreshold () const
 
double eMaxForFGCut () const
 
double eMaxForHoECut () const
 
double eMinForFGCut () const
 
double eMinForHoECut () const
 
double hActivityCut () const
 
double hMinForHoECut () const
 
double hOeCut () const
 
const std::vector< double > & jetMETECalScaleFactors () const
 
const std::vector< double > & jetMETHCalScaleFactors () const
 
double jetMETLSB () const
 
float JetMETTPGSum (const float &ecal, const float &hcal, const unsigned &iAbsEta) const
 
unsigned jscQuietThresholdBarrel () const
 
unsigned jscQuietThresholdEndcap () const
 
 L1RCTParameters ()
 
 L1RCTParameters (double eGammaLSB, double jetMETLSB, double eMinForFGCut, double eMaxForFGCut, double hOeCut, double eMinForHoECut, double eMaxForHoECut, double hMinForHoECut, double eActivityCut, double hActivityCut, unsigned eicIsolationThreshold, unsigned jscQuietThresholdBarrel, unsigned jscQuietThresholdEndcap, bool noiseVetoHB, bool noiseVetoHEplus, bool noiseVetoHEminus, bool useLindsey, const std::vector< double > &eGammaECalScaleFactors, const std::vector< double > &eGammaHCalScaleFactors, const std::vector< double > &jetMETECalScaleFactors, const std::vector< double > &jetMETHCalScaleFactors, const std::vector< double > &ecal_calib, const std::vector< double > &hcal_calib, const std::vector< double > &hcal_high_calib, const std::vector< double > &cross_terms, const std::vector< double > &lowHoverE_smear, const std::vector< double > &highHoverE_smear)
 
bool noiseVetoHB () const
 
bool noiseVetoHEminus () const
 
bool noiseVetoHEplus () const
 
void print (std::ostream &s) const
 
 ~L1RCTParameters ()
 

Private Member Functions

float correctedTPGSum (const float &ecal, const float &hcal, const unsigned &index) const
 

Private Attributes

std::vector< std::vector
< double > > 		cross_terms_
 
double eActivityCut_
 
std::vector< std::vector
< double > > 		ecal_calib_
 
std::vector< double > eGammaECalScaleFactors_
 
std::vector< double > eGammaHCalScaleFactors_
 
double eGammaLSB_
 
unsigned eicIsolationThreshold_
 
double eMaxForFGCut_
 
double eMaxForHoECut_
 
double eMinForFGCut_
 
double eMinForHoECut_
 
double hActivityCut_
 
std::vector< std::vector
< double > > 		hcal_calib_
 
std::vector< std::vector
< double > > 		hcal_high_calib_
 
double hMinForHoECut_
 
double hOeCut_
 
std::vector< double > HoverE_smear_high_
 
std::vector< double > HoverE_smear_low_
 
std::vector< double > jetMETECalScaleFactors_
 
std::vector< double > jetMETHCalScaleFactors_
 
double jetMETLSB_
 
unsigned jscQuietThresholdBarrel_
 
unsigned jscQuietThresholdEndcap_
 
bool noiseVetoHB_
 
bool noiseVetoHEminus_
 
bool noiseVetoHEplus_
 
bool useCorrections_
 

Detailed Description

Description: Class to contain parameters which define RCT Lookup Tables

Usage: <usage>

Definition at line 26 of file L1RCTParameters.h.

Constructor & Destructor Documentation

L1RCTParameters::L1RCTParameters ( )
inline

Definition at line 31 of file L1RCTParameters.h.

31 {}
L1RCTParameters::L1RCTParameters ( double  eGammaLSB,
double  jetMETLSB,
double  eMinForFGCut,
double  eMaxForFGCut,
double  hOeCut,
double  eMinForHoECut,
double  eMaxForHoECut,
double  hMinForHoECut,
double  eActivityCut,
double  hActivityCut,
unsigned  eicIsolationThreshold,
unsigned  jscQuietThresholdBarrel,
unsigned  jscQuietThresholdEndcap,
bool  noiseVetoHB,
bool  noiseVetoHEplus,
bool  noiseVetoHEminus,
bool  useLindsey,
const std::vector< double > &  eGammaECalScaleFactors,
const std::vector< double > &  eGammaHCalScaleFactors,
const std::vector< double > &  jetMETECalScaleFactors,
const std::vector< double > &  jetMETHCalScaleFactors,
const std::vector< double > &  ecal_calib,
const std::vector< double > &  hcal_calib,
const std::vector< double > &  hcal_high_calib,
const std::vector< double > &  cross_terms,
const std::vector< double > &  lowHoverE_smear,
const std::vector< double > &  highHoverE_smear 
)

Definition at line 17 of file L1RCTParameters.cc.

References cross_terms_, ecal_calib_, hcal_calib_, hcal_high_calib_, and i.

44  :
45  eGammaLSB_(eGammaLSB),
46  jetMETLSB_(jetMETLSB),
47  eMinForFGCut_(eMinForFGCut),
48  eMaxForFGCut_(eMaxForFGCut),
49  hOeCut_(hOeCut),
50  eMinForHoECut_(eMinForHoECut),
51  eMaxForHoECut_(eMaxForHoECut),
52  hMinForHoECut_(hMinForHoECut),
53  eActivityCut_(eActivityCut),
54  hActivityCut_(hActivityCut),
55  eicIsolationThreshold_(eicIsolationThreshold),
56  jscQuietThresholdBarrel_(jscQuietThresholdBarrel),
57  jscQuietThresholdEndcap_(jscQuietThresholdEndcap),
58  noiseVetoHB_(noiseVetoHB),
59  noiseVetoHEplus_(noiseVetoHEplus),
60  noiseVetoHEminus_(noiseVetoHEminus),
61  useCorrections_(useCorrections),
62  eGammaECalScaleFactors_(eGammaECalScaleFactors),
63  eGammaHCalScaleFactors_(eGammaHCalScaleFactors),
64  jetMETECalScaleFactors_(jetMETECalScaleFactors),
65  jetMETHCalScaleFactors_(jetMETHCalScaleFactors),
66  HoverE_smear_low_(lowHoverE_smear),
67  HoverE_smear_high_(highHoverE_smear)
68 {
69  ecal_calib_.resize(28);
70  hcal_calib_.resize(28);
71  hcal_high_calib_.resize(28);
72  cross_terms_.resize(28);
73 
74  for(unsigned i = 0; i < ecal_calib.size(); ++i)
75  ecal_calib_[i/3].push_back(ecal_calib[i]);
76  for(unsigned i = 0; i < hcal_calib.size(); ++i)
77  hcal_calib_[i/3].push_back(hcal_calib[i]);
78  for(unsigned i = 0; i < hcal_high_calib.size(); ++i)
79  hcal_high_calib_[i/3].push_back(hcal_high_calib[i]);
80  for(unsigned i = 0; i < cross_terms.size(); ++i)
81  cross_terms_[i/6].push_back(cross_terms[i]);
82 }
L1RCTParameters::eGammaHCalScaleFactors
const std::vector< double > & eGammaHCalScaleFactors() const
Definition: L1RCTParameters.h:84
L1RCTParameters::jetMETHCalScaleFactors
const std::vector< double > & jetMETHCalScaleFactors() const
Definition: L1RCTParameters.h:86
i
int i
Definition: DBlmapReader.cc:9
L1RCTParameters::eMaxForFGCut
double eMaxForFGCut() const
Definition: L1RCTParameters.h:70
L1RCTParameters::jetMETECalScaleFactors_
std::vector< double > jetMETECalScaleFactors_
Definition: L1RCTParameters.h:201
L1RCTParameters::eGammaECalScaleFactors
const std::vector< double > & eGammaECalScaleFactors() const
Definition: L1RCTParameters.h:83
L1RCTParameters::eicIsolationThreshold
unsigned eicIsolationThreshold() const
Definition: L1RCTParameters.h:77
L1RCTParameters::eMinForFGCut
double eMinForFGCut() const
Definition: L1RCTParameters.h:69
L1RCTParameters::jscQuietThresholdBarrel_
unsigned jscQuietThresholdBarrel_
Definition: L1RCTParameters.h:164
L1RCTParameters::noiseVetoHEminus
bool noiseVetoHEminus() const
Definition: L1RCTParameters.h:82
L1RCTParameters::eGammaECalScaleFactors_
std::vector< double > eGammaECalScaleFactors_
Definition: L1RCTParameters.h:194
L1RCTParameters::eActivityCut
double eActivityCut() const
Definition: L1RCTParameters.h:75
L1RCTParameters::cross_terms_
std::vector< std::vector< double > > cross_terms_
Definition: L1RCTParameters.h:215
L1RCTParameters::hcal_high_calib_
std::vector< std::vector< double > > hcal_high_calib_
Definition: L1RCTParameters.h:214
L1RCTParameters::HoverE_smear_low_
std::vector< double > HoverE_smear_low_
Definition: L1RCTParameters.h:219
L1RCTParameters::jscQuietThresholdEndcap
unsigned jscQuietThresholdEndcap() const
Definition: L1RCTParameters.h:79
L1RCTParameters::noiseVetoHEplus
bool noiseVetoHEplus() const
Definition: L1RCTParameters.h:81
L1RCTParameters::jscQuietThresholdEndcap_
unsigned jscQuietThresholdEndcap_
Definition: L1RCTParameters.h:169
L1RCTParameters::jetMETECalScaleFactors
const std::vector< double > & jetMETECalScaleFactors() const
Definition: L1RCTParameters.h:85
L1RCTParameters::hActivityCut
double hActivityCut() const
Definition: L1RCTParameters.h:76
L1RCTParameters::eMinForHoECut
double eMinForHoECut() const
Definition: L1RCTParameters.h:72
L1RCTParameters::jetMETHCalScaleFactors_
std::vector< double > jetMETHCalScaleFactors_
Definition: L1RCTParameters.h:202
L1RCTParameters::jetMETLSB
double jetMETLSB() const
Definition: L1RCTParameters.h:68
L1RCTParameters::hcal_calib_
std::vector< std::vector< double > > hcal_calib_
Definition: L1RCTParameters.h:213
L1RCTParameters::noiseVetoHB
bool noiseVetoHB() const
Definition: L1RCTParameters.h:80
L1RCTParameters::eGammaLSB
double eGammaLSB() const
Definition: L1RCTParameters.h:67
L1RCTParameters::hOeCut
double hOeCut() const
Definition: L1RCTParameters.h:71
L1RCTParameters::HoverE_smear_high_
std::vector< double > HoverE_smear_high_
Definition: L1RCTParameters.h:220
L1RCTParameters::eGammaHCalScaleFactors_
std::vector< double > eGammaHCalScaleFactors_
Definition: L1RCTParameters.h:195
L1RCTParameters::ecal_calib_
std::vector< std::vector< double > > ecal_calib_
Definition: L1RCTParameters.h:212
L1RCTParameters::eicIsolationThreshold_
unsigned eicIsolationThreshold_
Definition: L1RCTParameters.h:159
L1RCTParameters::eMaxForHoECut
double eMaxForHoECut() const
Definition: L1RCTParameters.h:73
L1RCTParameters::jscQuietThresholdBarrel
unsigned jscQuietThresholdBarrel() const
Definition: L1RCTParameters.h:78
L1RCTParameters::hMinForHoECut
double hMinForHoECut() const
Definition: L1RCTParameters.h:74
L1RCTParameters::~L1RCTParameters ( )
inline

Definition at line 63 of file L1RCTParameters.h.

63 {;}

Member Function Documentation

unsigned short L1RCTParameters::calcCard ( unsigned short  rct_iphi,
unsigned short  absIeta 
) const

Definition at line 96 of file L1RCTParameters.cc.

Referenced by L1RCT::digiInput().

98 {
99  unsigned short card = 999;
100  // Note absIeta counts from 1-32 (not 0-31)
101  if (absIeta <= 24){
102  card = ((absIeta-1)/8)*2 + (rct_iphi%8)/4;
103  }
104  // 25 <= absIeta <= 28 (card 6)
105  else if ((absIeta >= 25) && (absIeta <= 28)){
106  card = 6;
107  }
108  else{}
109  return card;
110 }
unsigned short L1RCTParameters::calcCrate ( unsigned short  rct_iphi,
short  ieta 
) const

Definition at line 85 of file L1RCTParameters.cc.

References funct::abs().

Referenced by L1RCT::digiInput(), and L1RCTLutWriter::writeRcLutFile().

86 {
87  unsigned short crate = rct_iphi/8;
88  if(abs(ieta) > 28) crate = rct_iphi / 2;
89  if (ieta > 0){
90  crate = crate + 9;
91  }
92  return crate;
93 }
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
unsigned short L1RCTParameters::calcIAbsEta ( unsigned short  iCrate,
unsigned short  iCard,
unsigned short  iTower 
) const

Definition at line 181 of file L1RCTParameters.cc.

Referenced by calcIEta().

183 {
184  unsigned short absIEta;
185  if(iCard < 6)
186  absIEta = (iCard / 2) * 8 + (iTower / 4) + 1; // rm -1 JLL 7Nov07
187  else if(iCard == 6) {
188  // card 6, region 0
189  if(iTower < 16) // 17->16
190  absIEta = 25 + iTower / 4; // rm -1
191  // card 6, region 1
192  else
193  absIEta = 28 - ((iTower - 16) / 4); // 17->16
194  }
195  // HF regions
196  else
197  absIEta = 29 + iTower % 4;
198  return absIEta;
199 }
short L1RCTParameters::calcIEta ( unsigned short  iCrate,
unsigned short  iCard,
unsigned short  iTower 
) const

Definition at line 148 of file L1RCTParameters.cc.

References calcIAbsEta().

Referenced by L1RCTLookupTables::lookup(), and RctInputTextToDigi::produce().

150 {
151  unsigned short absIEta = calcIAbsEta(iCrate, iCard, iTower);
152  short iEta;
153  if(iCrate < 9) iEta = -absIEta;
154  else iEta = absIEta;
155  return iEta;
156 }
L1RCTParameters::calcIAbsEta
unsigned short calcIAbsEta(unsigned short iCrate, unsigned short iCard, unsigned short iTower) const
Definition: L1RCTParameters.cc:181
unsigned short L1RCTParameters::calcIPhi ( unsigned short  iCrate,
unsigned short  iCard,
unsigned short  iTower 
) const

Definition at line 159 of file L1RCTParameters.cc.

Referenced by L1RCTLookupTables::lookup(), and RctInputTextToDigi::produce().

162 {
163  short iPhi;
164  if(iCard < 6)
165  iPhi = (iCrate % 9) * 8 + (iCard % 2) * 4 + (iTower % 4); // rm -1 7Nov07
166  else if(iCard == 6){
167  // region 0
168  if(iTower < 16) // 17->16
169  iPhi = (iCrate % 9) * 8 + (iTower % 4); // rm -1
170  // region 1
171  else
172  iPhi = (iCrate % 9) * 8 + ((iTower - 16) % 4) + 4; // 17 -> 16
173  }
174  // HF regions
175  else
176  iPhi = (iCrate % 9) * 2 + iTower / 4;
177  return iPhi;
178 }
unsigned short L1RCTParameters::calcTower ( unsigned short  rct_iphi,
unsigned short  absIeta 
) const

Definition at line 113 of file L1RCTParameters.cc.

Referenced by L1RCT::digiInput(), and L1RCTLutWriter::writeRcLutFile().

115 {
116  unsigned short tower = 999;
117  unsigned short iphi = rct_iphi;
118  unsigned short regionPhi = (iphi % 8)/4;
119 
120  // Note absIeta counts from 1-32 (not 0-31)
121  if (absIeta <= 24){
122  // assume iphi between 0 and 71; makes towers from 0-31, mod. 7Nov07
123  tower = ((absIeta-1)%8)*4 + (iphi%4); // REMOVED +1
124  }
125  // 25 <= absIeta <= 28 (card 6)
126  else if ((absIeta >= 25) && (absIeta <= 28)){
127  if (regionPhi == 0){
128  // towers from 0-31, modified 7Nov07 Jessica Leonard
129  tower = (absIeta-25)*4 + (iphi%4); // REMOVED +1
130  }
131  else {
132  tower = 28 + iphi % 4 + (25 - absIeta) * 4; // mod. 7Nov07 JLL
133  }
134  }
135  // absIeta >= 29 (HF regions)
136  else if ((absIeta >= 29) && (absIeta <= 32)){
137  // SPECIAL DEFINITION OF REGIONPHI FOR HF SINCE HF IPHI IS 0-17
138  // Sept. 19 J. Leonard
139  regionPhi = iphi % 2;
140  // HF MAPPING, just regions now, don't need to worry about towers
141  // just calling it "tower" for convenience
142  tower = (regionPhi) * 4 + absIeta - 29;
143  }
144  return tower;
145 }
float L1RCTParameters::correctedTPGSum ( const float &  ecal,
const float &  hcal,
const unsigned &  index 
) const
private

Definition at line 239 of file L1RCTParameters.cc.

References trackerHits::c, cross_terms_, alignCSCRings::e, patCandidatesForDimuonsSequences_cff::ecal, ecal_calib_, h, AnalysisDataFormats_SUSYBSMObjects::hc, patCandidatesForDimuonsSequences_cff::hcal, hcal_calib_, hcal_high_calib_, HoverE_smear_high_, HoverE_smear_low_, getHLTprescales::index, and funct::pow().

Referenced by EGammaTPGSum(), and JetMETTPGSum().

240 {
241  if(index >= 28 && ecal > 120 && hcal > 120) return (ecal + hcal); // return plain sum if outside of calibration range or index is too high
242 
243  // let's make sure we're asking for items that are there.
244  if(ecal_calib_.at(index).size() != 3 || hcal_calib_.at(index).size() != 3 ||
245  hcal_high_calib_.at(index).size() != 3 || cross_terms_.at(index).size() != 6 ||
246  HoverE_smear_high_.size() <= index || HoverE_smear_low_.size() <= index)
247  return (ecal+hcal);
248 
249  double e = ecal, h = hcal;
250  double ec = 0.0, hc = 0.0, c = 0.0;
251 
252  ec = (ecal_calib_.at(index).at(0)*std::pow(e,3.) +
253  ecal_calib_.at(index).at(1)*std::pow(e,2.) +
254  ecal_calib_.at(index).at(2)*e);
255 
256  if(e + h < 23)
257  {
258  hc = (hcal_calib_.at(index).at(0)*std::pow(h,3.) +
259  hcal_calib_.at(index).at(1)*std::pow(h,2.) +
260  hcal_calib_.at(index).at(2)*h);
261 
262  c = (cross_terms_.at(index).at(0)*std::pow(e, 2.)*h +
263  cross_terms_.at(index).at(1)*std::pow(h, 2.)*e +
264  cross_terms_.at(index).at(2)*e*h +
265  cross_terms_.at(index).at(3)*std::pow(e, 3.)*h +
266  cross_terms_.at(index).at(4)*std::pow(h, 3.)*e +
267  cross_terms_.at(index).at(5)*std::pow(h, 2.)*std::pow(e, 2.));
268  }
269  else
270  {
271  hc = (hcal_high_calib_.at(index).at(0)*std::pow(h,3.) +
272  hcal_high_calib_.at(index).at(1)*std::pow(h,2.) +
273  hcal_high_calib_.at(index).at(2)*h);
274  }
275 
276  if(h/(e+h) >= 0.05)
277  {
278  ec *= HoverE_smear_high_.at(index);
279  hc *= HoverE_smear_high_.at(index);
280  c *= HoverE_smear_high_.at(index);
281  }
282  else
283  {
284  ec *= HoverE_smear_low_.at(index);
285  }
286  return ec + hc + c;
287 }
L1RCTParameters::cross_terms_
std::vector< std::vector< double > > cross_terms_
Definition: L1RCTParameters.h:215
L1RCTParameters::hcal_high_calib_
std::vector< std::vector< double > > hcal_high_calib_
Definition: L1RCTParameters.h:214
L1RCTParameters::HoverE_smear_low_
std::vector< double > HoverE_smear_low_
Definition: L1RCTParameters.h:219
L1RCTParameters::hcal_calib_
std::vector< std::vector< double > > hcal_calib_
Definition: L1RCTParameters.h:213
L1RCTParameters::HoverE_smear_high_
std::vector< double > HoverE_smear_high_
Definition: L1RCTParameters.h:220
patCandidatesForDimuonsSequences_cff.ecal
tuple ecal
Definition: patCandidatesForDimuonsSequences_cff.py:70
h
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
L1RCTParameters::ecal_calib_
std::vector< std::vector< double > > ecal_calib_
Definition: L1RCTParameters.h:212
AnalysisDataFormats_SUSYBSMObjects::hc
susybsm::HSCParticleCollection hc
Definition: classes.h:25
patCandidatesForDimuonsSequences_cff.hcal
tuple hcal
Definition: patCandidatesForDimuonsSequences_cff.py:71
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
double L1RCTParameters::eActivityCut ( ) const
inline

Definition at line 75 of file L1RCTParameters.h.

References eActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

75 {return eActivityCut_;}
const std::vector<double>& L1RCTParameters::eGammaECalScaleFactors ( ) const
inline

Definition at line 83 of file L1RCTParameters.h.

References eGammaECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

83 {return eGammaECalScaleFactors_;}
L1RCTParameters::eGammaECalScaleFactors_
std::vector< double > eGammaECalScaleFactors_
Definition: L1RCTParameters.h:194
const std::vector<double>& L1RCTParameters::eGammaHCalScaleFactors ( ) const
inline

Definition at line 84 of file L1RCTParameters.h.

References eGammaHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

84 {return eGammaHCalScaleFactors_;}
L1RCTParameters::eGammaHCalScaleFactors_
std::vector< double > eGammaHCalScaleFactors_
Definition: L1RCTParameters.h:195
double L1RCTParameters::eGammaLSB ( ) const
inline

Definition at line 67 of file L1RCTParameters.h.

References eGammaLSB_.

Referenced by L1RCTLookupTables::convertEcal(), and L1RCTLookupTables::eGammaETCode().

67 {return eGammaLSB_;}
float L1RCTParameters::EGammaTPGSum ( const float &  ecal,
const float &  hcal,
const unsigned &  iAbsEta 
) const

Definition at line 221 of file L1RCTParameters.cc.

References correctedTPGSum(), eGammaECalScaleFactors_, eGammaHCalScaleFactors_, query::result, and useCorrections_.

Referenced by L1RCTLookupTables::eGammaETCode().

222 {
223  float ecal_c = ecal*eGammaECalScaleFactors_.at(iAbsEta-1);
224  float hcal_c = hcal*eGammaHCalScaleFactors_.at(iAbsEta-1);
225  float result = ecal_c + hcal_c;
226 
227  if(useCorrections_)
228  {
229  if(eGammaHCalScaleFactors_.at(iAbsEta-1) != 0)
230  hcal_c = hcal;
231 
232  result = correctedTPGSum(ecal_c,hcal_c,iAbsEta-1);
233  }
234 
235  return result;
236 }
L1RCTParameters::eGammaECalScaleFactors_
std::vector< double > eGammaECalScaleFactors_
Definition: L1RCTParameters.h:194
L1RCTParameters::correctedTPGSum
float correctedTPGSum(const float &ecal, const float &hcal, const unsigned &index) const
Definition: L1RCTParameters.cc:239
query.result
tuple result
Definition: query.py:137
patCandidatesForDimuonsSequences_cff.ecal
tuple ecal
Definition: patCandidatesForDimuonsSequences_cff.py:70
L1RCTParameters::eGammaHCalScaleFactors_
std::vector< double > eGammaHCalScaleFactors_
Definition: L1RCTParameters.h:195
patCandidatesForDimuonsSequences_cff.hcal
tuple hcal
Definition: patCandidatesForDimuonsSequences_cff.py:71
unsigned L1RCTParameters::eicIsolationThreshold ( ) const
inline

Definition at line 77 of file L1RCTParameters.h.

References eicIsolationThreshold_.

Referenced by L1RCTLutWriter::analyze(), and L1RCTElectronIsolationCard::calcElectronCandidates().

77 {return eicIsolationThreshold_;}
L1RCTParameters::eicIsolationThreshold_
unsigned eicIsolationThreshold_
Definition: L1RCTParameters.h:159
double L1RCTParameters::eMaxForFGCut ( ) const
inline

Definition at line 70 of file L1RCTParameters.h.

References eMaxForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

70 {return eMaxForFGCut_;}
double L1RCTParameters::eMaxForHoECut ( ) const
inline

Definition at line 73 of file L1RCTParameters.h.

References eMaxForHoECut_.

Referenced by L1RCTLookupTables::activityBit(), and L1RCTLookupTables::hOeFGVetoBit().

73 {return eMaxForHoECut_;}
double L1RCTParameters::eMinForFGCut ( ) const
inline

Definition at line 69 of file L1RCTParameters.h.

References eMinForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

69 {return eMinForFGCut_;}
double L1RCTParameters::eMinForHoECut ( ) const
inline

Definition at line 72 of file L1RCTParameters.h.

References eMinForHoECut_.

Referenced by L1RCTLookupTables::activityBit(), and L1RCTLookupTables::hOeFGVetoBit().

72 {return eMinForHoECut_;}
double L1RCTParameters::hActivityCut ( ) const
inline

Definition at line 76 of file L1RCTParameters.h.

References hActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

76 {return hActivityCut_;}
double L1RCTParameters::hMinForHoECut ( ) const
inline

Definition at line 74 of file L1RCTParameters.h.

References hMinForHoECut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

74 {return hMinForHoECut_;}
double L1RCTParameters::hOeCut ( ) const
inline

Definition at line 71 of file L1RCTParameters.h.

References hOeCut_.

Referenced by L1RCTLookupTables::activityBit(), and L1RCTLookupTables::hOeFGVetoBit().

71 {return hOeCut_;}
const std::vector<double>& L1RCTParameters::jetMETECalScaleFactors ( ) const
inline

Definition at line 85 of file L1RCTParameters.h.

References jetMETECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

85 {return jetMETECalScaleFactors_;}
L1RCTParameters::jetMETECalScaleFactors_
std::vector< double > jetMETECalScaleFactors_
Definition: L1RCTParameters.h:201
const std::vector<double>& L1RCTParameters::jetMETHCalScaleFactors ( ) const
inline

Definition at line 86 of file L1RCTParameters.h.

References jetMETHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

86 {return jetMETHCalScaleFactors_;}
L1RCTParameters::jetMETHCalScaleFactors_
std::vector< double > jetMETHCalScaleFactors_
Definition: L1RCTParameters.h:202
double L1RCTParameters::jetMETLSB ( ) const
inline

Definition at line 68 of file L1RCTParameters.h.

References jetMETLSB_.

Referenced by L1RCTLookupTables::convertHcal(), L1RCTLookupTables::jetMETETCode(), and L1RCTLookupTables::lookup().

68 {return jetMETLSB_;}
float L1RCTParameters::JetMETTPGSum ( const float &  ecal,
const float &  hcal,
const unsigned &  iAbsEta 
) const

Definition at line 201 of file L1RCTParameters.cc.

References correctedTPGSum(), eGammaECalScaleFactors_, jetMETECalScaleFactors_, jetMETHCalScaleFactors_, query::result, and useCorrections_.

Referenced by L1RCTLookupTables::jetMETETCode().

202 {
203  float ecal_c = ecal*jetMETECalScaleFactors_.at(iAbsEta-1);
204  float hcal_c = hcal*jetMETHCalScaleFactors_.at(iAbsEta-1);
205  float result = ecal_c + hcal_c;
206 
207  if(useCorrections_)
208  {
209  if(jetMETHCalScaleFactors_.at(iAbsEta-1) != 0)
210  hcal_c = hcal;
211 
212  if(jetMETECalScaleFactors_.at(iAbsEta-1) != 0)
213  ecal_c = ecal*eGammaECalScaleFactors_.at(iAbsEta-1); // Use eGamma Corrections
214 
215  result = correctedTPGSum(ecal_c,hcal_c,iAbsEta-1);
216  }
217 
218  return result;
219 }
L1RCTParameters::jetMETECalScaleFactors_
std::vector< double > jetMETECalScaleFactors_
Definition: L1RCTParameters.h:201
L1RCTParameters::eGammaECalScaleFactors_
std::vector< double > eGammaECalScaleFactors_
Definition: L1RCTParameters.h:194
L1RCTParameters::correctedTPGSum
float correctedTPGSum(const float &ecal, const float &hcal, const unsigned &index) const
Definition: L1RCTParameters.cc:239
L1RCTParameters::jetMETHCalScaleFactors_
std::vector< double > jetMETHCalScaleFactors_
Definition: L1RCTParameters.h:202
query.result
tuple result
Definition: query.py:137
patCandidatesForDimuonsSequences_cff.ecal
tuple ecal
Definition: patCandidatesForDimuonsSequences_cff.py:70
patCandidatesForDimuonsSequences_cff.hcal
tuple hcal
Definition: patCandidatesForDimuonsSequences_cff.py:71
unsigned L1RCTParameters::jscQuietThresholdBarrel ( ) const
inline

Definition at line 78 of file L1RCTParameters.h.

References jscQuietThresholdBarrel_.

Referenced by L1RCTLutWriter::analyze(), and L1RCTJetSummaryCard::fillQuietBits().

78 {return jscQuietThresholdBarrel_;}
L1RCTParameters::jscQuietThresholdBarrel_
unsigned jscQuietThresholdBarrel_
Definition: L1RCTParameters.h:164
unsigned L1RCTParameters::jscQuietThresholdEndcap ( ) const
inline

Definition at line 79 of file L1RCTParameters.h.

References jscQuietThresholdEndcap_.

Referenced by L1RCTLutWriter::analyze(), and L1RCTJetSummaryCard::fillQuietBits().

79 {return jscQuietThresholdEndcap_;}
L1RCTParameters::jscQuietThresholdEndcap_
unsigned jscQuietThresholdEndcap_
Definition: L1RCTParameters.h:169
bool L1RCTParameters::noiseVetoHB ( ) const
inline

Definition at line 80 of file L1RCTParameters.h.

References noiseVetoHB_.

Referenced by L1RCTLookupTables::lookup().

80 {return noiseVetoHB_;}
bool L1RCTParameters::noiseVetoHEminus ( ) const
inline

Definition at line 82 of file L1RCTParameters.h.

References noiseVetoHEminus_.

Referenced by L1RCTLookupTables::lookup().

82 {return noiseVetoHEminus_;}
bool L1RCTParameters::noiseVetoHEplus ( ) const
inline

Definition at line 81 of file L1RCTParameters.h.

References noiseVetoHEplus_.

Referenced by L1RCTLookupTables::lookup().

81 {return noiseVetoHEplus_;}
void L1RCTParameters::print ( std::ostream &  s) const

Definition at line 290 of file L1RCTParameters.cc.

References cross_terms_, eActivityCut_, ecal_calib_, eGammaECalScaleFactors_, eGammaHCalScaleFactors_, eGammaLSB_, eicIsolationThreshold_, eMaxForFGCut_, eMaxForHoECut_, eMinForFGCut_, eMinForHoECut_, end, hActivityCut_, hcal_calib_, hcal_high_calib_, hMinForHoECut_, hOeCut_, HoverE_smear_high_, HoverE_smear_low_, i, j, jetMETECalScaleFactors_, jetMETHCalScaleFactors_, jetMETLSB_, jscQuietThresholdBarrel_, jscQuietThresholdEndcap_, noiseVetoHB_, noiseVetoHEminus_, noiseVetoHEplus_, and useCorrections_.

290  {
291 
292  s << "\nPrinting record L1RCTParametersRcd" <<endl;
293  s << "\n\"Parameter description\" \n \"Parameter name\" \"Value\" "
294  << endl;
295  s << "\ne/gamma least significant bit energy transmitted from receiver cards to EIC cards. \n "
296  << "eGammaLSB = " << eGammaLSB_ << endl ;
297  s << "\nLSB of region Et scale from RCT to GCT (GeV) \n "
298  << "jetMETLSB = " << jetMETLSB_ << endl ;
299  s << "\nminimum ECAL Et for which fine-grain veto is applied (GeV) \n "
300  << " eMinForFGCut = " << eMinForFGCut_ << endl ;
301  s << "\nmaximum ECAL Et for which fine-grain veto is applied (GeV) \n "
302  << "eMaxForFGCut = " << eMaxForFGCut_ << endl ;
303  s << "\nmaximum value of (HCAL Et / ECAL Et) \n "
304  << "hOeCut = " << hOeCut_ << endl ;
305  s << "\nminimum ECAL Et for which H/E veto is applied (GeV) \n "
306  << "eMinForHoECut = " << eMinForHoECut_ << endl ;
307  s << "\nmaximum ECAL Et for which H/E veto is applied (GeV) \n "
308  << "eMaxForHoECut = " << eMaxForHoECut_ << endl ;
309  s << "\nminimum HCAL Et for which H/E veto is applied (GeV) \n "
310  << "hMinForHoECut = " << hMinForHoECut_ << endl ;
311  s << "\nECAL Et threshold above which tau activity bit is set (GeV) \n "
312  << "eActivityCut = " << eActivityCut_ << endl ;
313  s << "\nHCAL Et threshold above which tau activity bit is set (GeV) \n "
314  << "hActivityCut = " << hActivityCut_ << endl ;
315  s << "\nNeighboring trigger tower energy minimum threshold that marks candidate as non-isolated. (LSB bits) \n "
316  << "eicIsolationThreshold = " << eicIsolationThreshold_ << endl ;
317  s << "\nIf jetMet energy in RCT Barrel Region is below this value, a quiet bit is set. (LSB bits)\n "
318  << "jscQuietThreshBarrel = " << jscQuietThresholdBarrel_ << endl ;
319  s << "\nIf jetMet energy in RCT Endcap Region is below this value, a quiet bit is set. (LSB bits) \n "
320  << "jscQuietThreshEndcap = " << jscQuietThresholdEndcap_ << endl ;
321  s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT Barrel \n "
322  << "noiseVetoHB = " << noiseVetoHB_ << endl ;
323  s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT Encap+ \n "
324  << "noiseVetoHEplus = " << noiseVetoHEplus_ << endl ;
325  s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT Endcap- \n "
326  << "noiseVetoHEminus = " << noiseVetoHEminus_ << endl ;
327 
328  s << "\n\neta-dependent multiplicative factors for ECAL Et before summation \n "
329  << "eGammaECal Scale Factors " << endl;
330  s << "ieta ScaleFactor" <<endl;
331  for(int i = 0 ; i<28; i++)
332  s << setw(4) << i+1 << " " << eGammaECalScaleFactors_.at(i) <<endl;
333 
334  s << "\n\neta-dependent multiplicative factors for HCAL Et before summation \n "
335  <<"eGammaHCal Scale Factors "<<endl;
336  s << "ieta ScaleFactor" <<endl;
337  for(int i = 0 ; i<28; i++)
338  s << setw(4) << i+1 << " " << eGammaHCalScaleFactors_.at(i) <<endl;
339 
340  s << "\n\neta-dependent multiplicative factors for ECAL Et before summation \n "
341  <<"jetMETECal Scale Factors "<<endl;
342  s << "ieta ScaleFactor" <<endl;
343  for(int i = 0 ; i<28; i++)
344  s<< setw(4) << i+1 << " " << jetMETECalScaleFactors_.at(i) <<endl;
345 
346  s << "\n\neta-dependent multiplicative factors for HCAL Et before summation \n"
347  <<"jetMETHCal Scale Factors "<<endl;
348  s << "ieta ScaleFactor" <<endl;
349  for(int i = 0 ; i<28; i++)
350  s << setw(4) <<i+1 << " " << jetMETHCalScaleFactors_.at(i) <<endl;
351 
352 
353  if(useCorrections_) {
354  s<< "\n\nUSING calibration variables " <<endl;
355 
356 
357  s << "\n\nH over E smear low Correction Factors "<<endl;
358  s << "ieta Correction Factor" <<endl;
359  for(int i = 0 ; i<28; i++)
360  s << setw(4) <<i+1 << " " << HoverE_smear_low_.at(i) <<endl;
361 
362 
363  s << "\n\nH over E smear high Correction Factors "<<endl;
364  s <<"ieta Correction Factor" <<endl;
365  for(int i = 0 ; i<28; i++)
366  s << setw(4) <<i+1 << " " << HoverE_smear_high_.at(i) <<endl;
367 
368  s << "\n\necal calibrations "<<endl;
369  s << "ieta CorrFactor1 CorrFactor2 CorrFactor3" <<endl;
370  int end = ecal_calib_[0].size();
371  for(int i = 0 ; i<28; i++) {
372  s << setw(4) << i;
373  for(int j = 0; j< end ; j++)
374  s << setw(11) << setprecision(8) << ecal_calib_[i][j] ;
375 
376  s << endl;
377 
378  }
379 
380  s <<"\n\nhcal calibrations "<<endl;
381  s <<"ieta CorrFactor1 CorrFactor2 CorrFactor3" <<endl;
382  end = hcal_calib_[0].size();
383  for(int i = 0 ; i<28; i++) {
384  s << setw(4) << i;
385  for(int j = 0; j< end ; j++)
386  s << setw(11) << setprecision(8) << hcal_calib_[i][j] ;
387 
388  s << endl;
389 
390  }
391  s <<"\n\nhcal_high calibrations "<<endl;
392  s <<"ieta CorrFactor1 CorrFactor2 CorrFactor3" <<endl;
393  end = hcal_high_calib_[0].size();
394  for(int i = 0 ; i<28; i++) {
395  s << setw(4) << i;
396  for(int j = 0; j< end ; j++)
397  s << setw(11) << setprecision(8) << hcal_high_calib_[i][j] ;
398 
399  s << endl;
400 
401  }
402  end = cross_terms_[0].size();
403  s <<"\n\ncross terms calibrations "<<endl;
404  s <<"ieta CorrFactor1 CorrFactor2 CorrFactor3 CorrFactor4 CorrFactor5 CorrFactor6" <<endl;
405  for(int i = 0 ; i<28; i++) {
406  s << setw(4) << i;
407  for(int j = 0; j< end ; j++)
408  s << setw(11) << setprecision(8) << cross_terms_[i][j] ;
409 
410  s << endl;
411 
412  }
413 
414  } else
415  s<< "\n\nNOT USING calibration variables " <<endl;
416 
417  s << "\n\n" <<endl;
418 
419 }
i
int i
Definition: DBlmapReader.cc:9
L1RCTParameters::jetMETECalScaleFactors_
std::vector< double > jetMETECalScaleFactors_
Definition: L1RCTParameters.h:201
L1RCTParameters::jscQuietThresholdBarrel_
unsigned jscQuietThresholdBarrel_
Definition: L1RCTParameters.h:164
L1RCTParameters::eGammaECalScaleFactors_
std::vector< double > eGammaECalScaleFactors_
Definition: L1RCTParameters.h:194
L1RCTParameters::cross_terms_
std::vector< std::vector< double > > cross_terms_
Definition: L1RCTParameters.h:215
L1RCTParameters::hcal_high_calib_
std::vector< std::vector< double > > hcal_high_calib_
Definition: L1RCTParameters.h:214
L1RCTParameters::HoverE_smear_low_
std::vector< double > HoverE_smear_low_
Definition: L1RCTParameters.h:219
L1RCTParameters::jscQuietThresholdEndcap_
unsigned jscQuietThresholdEndcap_
Definition: L1RCTParameters.h:169
L1RCTParameters::jetMETHCalScaleFactors_
std::vector< double > jetMETHCalScaleFactors_
Definition: L1RCTParameters.h:202
j
int j
Definition: DBlmapReader.cc:9
L1RCTParameters::hcal_calib_
std::vector< std::vector< double > > hcal_calib_
Definition: L1RCTParameters.h:213
end
#define end
Definition: vmac.h:37
L1RCTParameters::HoverE_smear_high_
std::vector< double > HoverE_smear_high_
Definition: L1RCTParameters.h:220
L1RCTParameters::eGammaHCalScaleFactors_
std::vector< double > eGammaHCalScaleFactors_
Definition: L1RCTParameters.h:195
L1RCTParameters::ecal_calib_
std::vector< std::vector< double > > ecal_calib_
Definition: L1RCTParameters.h:212
L1RCTParameters::eicIsolationThreshold_
unsigned eicIsolationThreshold_
Definition: L1RCTParameters.h:159

Member Data Documentation

std::vector<std::vector<double> > L1RCTParameters::cross_terms_
private

Definition at line 215 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), L1RCTParameters(), and print().

double L1RCTParameters::eActivityCut_
private

Definition at line 146 of file L1RCTParameters.h.

Referenced by eActivityCut(), and print().

std::vector<std::vector<double> > L1RCTParameters::ecal_calib_
private

Definition at line 212 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), L1RCTParameters(), and print().

std::vector<double> L1RCTParameters::eGammaECalScaleFactors_
private

Definition at line 194 of file L1RCTParameters.h.

Referenced by eGammaECalScaleFactors(), EGammaTPGSum(), JetMETTPGSum(), and print().

std::vector<double> L1RCTParameters::eGammaHCalScaleFactors_
private

Definition at line 195 of file L1RCTParameters.h.

Referenced by eGammaHCalScaleFactors(), EGammaTPGSum(), and print().

double L1RCTParameters::eGammaLSB_
private

Definition at line 112 of file L1RCTParameters.h.

Referenced by eGammaLSB(), and print().

unsigned L1RCTParameters::eicIsolationThreshold_
private

Definition at line 159 of file L1RCTParameters.h.

Referenced by eicIsolationThreshold(), and print().

double L1RCTParameters::eMaxForFGCut_
private

Definition at line 124 of file L1RCTParameters.h.

Referenced by eMaxForFGCut(), and print().

double L1RCTParameters::eMaxForHoECut_
private

Definition at line 136 of file L1RCTParameters.h.

Referenced by eMaxForHoECut(), and print().

double L1RCTParameters::eMinForFGCut_
private

Definition at line 120 of file L1RCTParameters.h.

Referenced by eMinForFGCut(), and print().

double L1RCTParameters::eMinForHoECut_
private

Definition at line 132 of file L1RCTParameters.h.

Referenced by eMinForHoECut(), and print().

double L1RCTParameters::hActivityCut_
private

Definition at line 152 of file L1RCTParameters.h.

Referenced by hActivityCut(), and print().

std::vector<std::vector<double> > L1RCTParameters::hcal_calib_
private

Definition at line 213 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), L1RCTParameters(), and print().

std::vector<std::vector<double> > L1RCTParameters::hcal_high_calib_
private

Definition at line 214 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), L1RCTParameters(), and print().

double L1RCTParameters::hMinForHoECut_
private

Definition at line 140 of file L1RCTParameters.h.

Referenced by hMinForHoECut(), and print().

double L1RCTParameters::hOeCut_
private

Definition at line 128 of file L1RCTParameters.h.

Referenced by hOeCut(), and print().

std::vector<double> L1RCTParameters::HoverE_smear_high_
private

Definition at line 220 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), and print().

std::vector<double> L1RCTParameters::HoverE_smear_low_
private

Definition at line 219 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), and print().

std::vector<double> L1RCTParameters::jetMETECalScaleFactors_
private

Definition at line 201 of file L1RCTParameters.h.

Referenced by jetMETECalScaleFactors(), JetMETTPGSum(), and print().

std::vector<double> L1RCTParameters::jetMETHCalScaleFactors_
private

Definition at line 202 of file L1RCTParameters.h.

Referenced by jetMETHCalScaleFactors(), JetMETTPGSum(), and print().

double L1RCTParameters::jetMETLSB_
private

Definition at line 116 of file L1RCTParameters.h.

Referenced by jetMETLSB(), and print().

unsigned L1RCTParameters::jscQuietThresholdBarrel_
private

Definition at line 164 of file L1RCTParameters.h.

Referenced by jscQuietThresholdBarrel(), and print().

unsigned L1RCTParameters::jscQuietThresholdEndcap_
private

Definition at line 169 of file L1RCTParameters.h.

Referenced by jscQuietThresholdEndcap(), and print().

bool L1RCTParameters::noiseVetoHB_
private

Definition at line 175 of file L1RCTParameters.h.

Referenced by noiseVetoHB(), and print().

bool L1RCTParameters::noiseVetoHEminus_
private

Definition at line 185 of file L1RCTParameters.h.

Referenced by noiseVetoHEminus(), and print().

bool L1RCTParameters::noiseVetoHEplus_
private

Definition at line 180 of file L1RCTParameters.h.

Referenced by noiseVetoHEplus(), and print().

bool L1RCTParameters::useCorrections_
private

Definition at line 188 of file L1RCTParameters.h.

Referenced by EGammaTPGSum(), JetMETTPGSum(), and print().

Generated for CMSSW Reference Manual by doxygen 1.8.5