#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

template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

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_

Friends
class	boost::serialization::access

template<typename CondSerializationT , typename Enabled >
struct	cond::serialization::access

Detailed Description

Description: Class to contain parameters which define RCT Lookup Tables

Usage: <usage>

Definition at line 27 of file L1RCTParameters.h.

Constructor & Destructor Documentation

◆ L1RCTParameters() [1/2]

L1RCTParameters::L1RCTParameters ( )

inline

Definition at line 30 of file L1RCTParameters.h.

30 {}

◆ L1RCTParameters() [2/2]

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_, mps_fire::i, and push_back().

     : eGammaLSB_(eGammaLSB),
       jetMETLSB_(jetMETLSB),
       eMinForFGCut_(eMinForFGCut),
       eMaxForFGCut_(eMaxForFGCut),
       hOeCut_(hOeCut),
       eMinForHoECut_(eMinForHoECut),
       eMaxForHoECut_(eMaxForHoECut),
       hMinForHoECut_(hMinForHoECut),
       eActivityCut_(eActivityCut),
       hActivityCut_(hActivityCut),
       eicIsolationThreshold_(eicIsolationThreshold),
       jscQuietThresholdBarrel_(jscQuietThresholdBarrel),
       jscQuietThresholdEndcap_(jscQuietThresholdEndcap),
       noiseVetoHB_(noiseVetoHB),
       noiseVetoHEplus_(noiseVetoHEplus),
       noiseVetoHEminus_(noiseVetoHEminus),
       useCorrections_(useCorrections),
       eGammaECalScaleFactors_(eGammaECalScaleFactors),
       eGammaHCalScaleFactors_(eGammaHCalScaleFactors),
       jetMETECalScaleFactors_(jetMETECalScaleFactors),
       jetMETHCalScaleFactors_(jetMETHCalScaleFactors),
       HoverE_smear_low_(lowHoverE_smear),
       HoverE_smear_high_(highHoverE_smear) {
   ecal_calib_.resize(28);
   hcal_calib_.resize(28);
   hcal_high_calib_.resize(28);
   cross_terms_.resize(28);
 
   for (unsigned i = 0; i < ecal_calib.size(); ++i)
     ecal_calib_[i / 3].push_back(ecal_calib[i]);
   for (unsigned i = 0; i < hcal_calib.size(); ++i)
     hcal_calib_[i / 3].push_back(hcal_calib[i]);
   for (unsigned i = 0; i < hcal_high_calib.size(); ++i)
     hcal_high_calib_[i / 3].push_back(hcal_high_calib[i]);
   for (unsigned i = 0; i < cross_terms.size(); ++i)
     cross_terms_[i / 6].push_back(cross_terms[i]);
 }

◆ ~L1RCTParameters()

L1RCTParameters::~L1RCTParameters ( )

inline

Definition at line 61 of file L1RCTParameters.h.

61 { ; }

Member Function Documentation

◆ calcCard()

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

Definition at line 94 of file L1RCTParameters.cc.

Referenced by L1RCT::digiInput().

                                                                                               {
   unsigned short card = 999;
   // Note absIeta counts from 1-32 (not 0-31)
   if (absIeta <= 24) {
     card = ((absIeta - 1) / 8) * 2 + (rct_iphi % 8) / 4;
   }
   // 25 <= absIeta <= 28 (card 6)
   else if ((absIeta >= 25) && (absIeta <= 28)) {
     card = 6;
   } else {
   }
   return card;
 }

◆ calcCrate()

unsigned short L1RCTParameters::calcCrate	(	unsigned short	rct_iphi,
		short	ieta
	)		const

Definition at line 83 of file L1RCTParameters.cc.

References funct::abs(), and hcalRecHitTable_cff::ieta.

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

                                                                                    {
   unsigned short crate = rct_iphi / 8;
   if (abs(ieta) > 28)
     crate = rct_iphi / 2;
   if (ieta > 0) {
     crate = crate + 9;
   }
   return crate;
 }

◆ calcIAbsEta()

unsigned short L1RCTParameters::calcIAbsEta	(	unsigned short	iCrate,
		unsigned short	iCard,
		unsigned short	iTower
	)		const

Definition at line 171 of file L1RCTParameters.cc.

Referenced by calcIEta().

                                                                                                                     {
   unsigned short absIEta;
   if (iCard < 6)
     absIEta = (iCard / 2) * 8 + (iTower / 4) + 1;  // rm -1 JLL 7Nov07
   else if (iCard == 6) {
     // card 6, region 0
     if (iTower < 16)              // 17->16
       absIEta = 25 + iTower / 4;  // rm -1
     // card 6, region 1
     else
       absIEta = 28 - ((iTower - 16) / 4);  // 17->16
   }
   // HF regions
   else
     absIEta = 29 + iTower % 4;
   return absIEta;
 }

◆ calcIEta()

short L1RCTParameters::calcIEta	(	unsigned short	iCrate,
		unsigned short	iCard,
		unsigned short	iTower
	)		const

Definition at line 141 of file L1RCTParameters.cc.

References calcIAbsEta(), and l1tTowerCalibrationProducer_cfi::iEta.

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

                                                                                                         {
   unsigned short absIEta = calcIAbsEta(iCrate, iCard, iTower);
   short iEta;
   if (iCrate < 9)
     iEta = -absIEta;
   else
     iEta = absIEta;
   return iEta;
 }

◆ calcIPhi()

unsigned short L1RCTParameters::calcIPhi	(	unsigned short	iCrate,
		unsigned short	iCard,
		unsigned short	iTower
	)		const

Definition at line 152 of file L1RCTParameters.cc.

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

                                                                                                                  {
   short iPhi;
   if (iCard < 6)
     iPhi = (iCrate % 9) * 8 + (iCard % 2) * 4 + (iTower % 4);  // rm -1 7Nov07
   else if (iCard == 6) {
     // region 0
     if (iTower < 16)                           // 17->16
       iPhi = (iCrate % 9) * 8 + (iTower % 4);  // rm -1
     // region 1
     else
       iPhi = (iCrate % 9) * 8 + ((iTower - 16) % 4) + 4;  // 17 -> 16
   }
   // HF regions
   else
     iPhi = (iCrate % 9) * 2 + iTower / 4;
   return iPhi;
 }

◆ calcTower()

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

Definition at line 109 of file L1RCTParameters.cc.

References hcalRecHitTable_cff::iphi, and l1tHGCalTowerProducer_cfi::tower.

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

                                                                                                {
   unsigned short tower = 999;
   unsigned short iphi = rct_iphi;
   unsigned short regionPhi = (iphi % 8) / 4;
 
   // Note absIeta counts from 1-32 (not 0-31)
   if (absIeta <= 24) {
     // assume iphi between 0 and 71; makes towers from 0-31, mod. 7Nov07
     tower = ((absIeta - 1) % 8) * 4 + (iphi % 4);  // REMOVED +1
   }
   // 25 <= absIeta <= 28 (card 6)
   else if ((absIeta >= 25) && (absIeta <= 28)) {
     if (regionPhi == 0) {
       // towers from 0-31, modified 7Nov07 Jessica Leonard
       tower = (absIeta - 25) * 4 + (iphi % 4);  // REMOVED +1
     } else {
       tower = 28 + iphi % 4 + (25 - absIeta) * 4;  // mod. 7Nov07 JLL
     }
   }
   // absIeta >= 29 (HF regions)
   else if ((absIeta >= 29) && (absIeta <= 32)) {
     // SPECIAL DEFINITION OF REGIONPHI FOR HF SINCE HF IPHI IS 0-17
     // Sept. 19 J. Leonard
     regionPhi = iphi % 2;
     // HF MAPPING, just regions now, don't need to worry about towers
     // just calling it "tower" for convenience
     tower = (regionPhi)*4 + absIeta - 29;
   }
   return tower;
 }

◆ correctedTPGSum()

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

private

Definition at line 285 of file L1RCTParameters.cc.

References HltBtagPostValidation_cff::c, cross_terms_, MillePedeFileConverter_cfg::e, bsc_activity_cfg::ecal, ecal_calib_, h, hltEgammaHLTExtra_cfi::hcal, hcal_calib_, hcal_high_calib_, HoverE_smear_high_, HoverE_smear_low_, and conifer::pow().

Referenced by EGammaTPGSum(), and JetMETTPGSum().

                                                                                                         {
   if (index >= 28 && ecal > 120 && hcal > 120)
     return (ecal + hcal);  // return plain sum if outside of calibration range or index is too high
 
   // let's make sure we're asking for items that are there.
   if (ecal_calib_.at(index).size() != 3 || hcal_calib_.at(index).size() != 3 ||
       hcal_high_calib_.at(index).size() != 3 || cross_terms_.at(index).size() != 6 ||
       HoverE_smear_high_.size() <= index || HoverE_smear_low_.size() <= index)
     return (ecal + hcal);
 
   double e = ecal, h = hcal;
   double ec = 0.0, hc = 0.0, c = 0.0;
 
   ec = (ecal_calib_.at(index).at(0) * std::pow(e, 3.) + ecal_calib_.at(index).at(1) * std::pow(e, 2.) +
         ecal_calib_.at(index).at(2) * e);
 
   if (e + h < 23) {
     hc = (hcal_calib_.at(index).at(0) * std::pow(h, 3.) + hcal_calib_.at(index).at(1) * std::pow(h, 2.) +
           hcal_calib_.at(index).at(2) * h);
 
     c = (cross_terms_.at(index).at(0) * std::pow(e, 2.) * h + cross_terms_.at(index).at(1) * std::pow(h, 2.) * e +
          cross_terms_.at(index).at(2) * e * h + cross_terms_.at(index).at(3) * std::pow(e, 3.) * h +
          cross_terms_.at(index).at(4) * std::pow(h, 3.) * e +
          cross_terms_.at(index).at(5) * std::pow(h, 2.) * std::pow(e, 2.));
   } else {
     hc = (hcal_high_calib_.at(index).at(0) * std::pow(h, 3.) + hcal_high_calib_.at(index).at(1) * std::pow(h, 2.) +
           hcal_high_calib_.at(index).at(2) * h);
   }
 
   if (h / (e + h) >= 0.05) {
     ec *= HoverE_smear_high_.at(index);
     hc *= HoverE_smear_high_.at(index);
     c *= HoverE_smear_high_.at(index);
   } else {
     ec *= HoverE_smear_low_.at(index);
   }
   return ec + hc + c;
 }

◆ eActivityCut()

double L1RCTParameters::eActivityCut ( ) const

inline

Definition at line 73 of file L1RCTParameters.h.

References eActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

73 { return eActivityCut_; }

L1RCTParameters::eActivityCut_

double eActivityCut_

Definition: L1RCTParameters.h:143

◆ eGammaECalScaleFactors()

const std::vector<double>& L1RCTParameters::eGammaECalScaleFactors ( ) const

inline

Definition at line 81 of file L1RCTParameters.h.

References eGammaECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

81 { return eGammaECalScaleFactors_; }

L1RCTParameters::eGammaECalScaleFactors_

std::vector< double > eGammaECalScaleFactors_

Definition: L1RCTParameters.h:191

◆ eGammaHCalScaleFactors()

const std::vector<double>& L1RCTParameters::eGammaHCalScaleFactors ( ) const

inline

Definition at line 82 of file L1RCTParameters.h.

References eGammaHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

82 { return eGammaHCalScaleFactors_; }

L1RCTParameters::eGammaHCalScaleFactors_

std::vector< double > eGammaHCalScaleFactors_

Definition: L1RCTParameters.h:192

◆ eGammaLSB()

double L1RCTParameters::eGammaLSB ( ) const

inline

Definition at line 65 of file L1RCTParameters.h.

References eGammaLSB_.

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

65 { return eGammaLSB_; }

L1RCTParameters::eGammaLSB_

double eGammaLSB_

Definition: L1RCTParameters.h:109

◆ EGammaTPGSum()

float L1RCTParameters::EGammaTPGSum	(	const float &	ecal,
		const float &	hcal,
		const unsigned &	iAbsEta
	)		const

Definition at line 240 of file L1RCTParameters.cc.

References correctedTPGSum(), eGammaECalScaleFactors_, eGammaHCalScaleFactors_, Exception, createfilelist::int, mps_fire::result, and useCorrections_.

Referenced by L1RCTLookupTables::eGammaETCode().

                                                                                                        {
   // We never deal with HF in this function (EG objects don't use hcal at all, really)
   if (iAbsEta < 1 || iAbsEta > 28)
     throw cms::Exception("L1RCTParameters invalid function call") << "Eta out of range in MET TPGSum: " << iAbsEta;
   float ecal_c = ecal * eGammaECalScaleFactors_.at(iAbsEta - 1);
   float hcal_c = hcal * eGammaHCalScaleFactors_.at(iAbsEta - 1);
 
   // scale factors will either be length 28 for legacy, or 28*(# et bins + 1) where the first set of 28 is an average over the et bins
   // The first set of 28 provides a legacy fallthrough option
   // Currently, # et bins is 9
   if (eGammaECalScaleFactors_.size() == 28 * 10) {
     int et_bin = ((int)floor(ecal) / 5);
     // lowest bin (1) is 0-10GeV
     if (et_bin < 1)
       et_bin = 1;
     // highest bin (9) is 45GeV and up
     if (et_bin > 9)
       et_bin = 9;
     ecal_c = ecal * eGammaECalScaleFactors_.at(et_bin * 28 + iAbsEta - 1);
   }
   if (eGammaHCalScaleFactors_.size() == 28 * 10) {
     int ht_bin = ((int)floor(hcal) / 5);
     // lowest bin (1) is 0-10GeV
     if (ht_bin < 1)
       ht_bin = 1;
     // highest bin (9) is 45GeV and up
     if (ht_bin > 9)
       ht_bin = 9;
     hcal_c = hcal * eGammaHCalScaleFactors_.at(ht_bin * 28 + iAbsEta - 1);
   }
 
   float result = ecal_c + hcal_c;
 
   // defunct section (polynomial-parameterized corrections)
   if (useCorrections_) {
     if (eGammaHCalScaleFactors_.at(iAbsEta - 1) != 0)
       hcal_c = hcal;
 
     result = correctedTPGSum(ecal_c, hcal_c, iAbsEta - 1);
   }
 
   return result;
 }

◆ eicIsolationThreshold()

unsigned L1RCTParameters::eicIsolationThreshold ( ) const

inline

Definition at line 75 of file L1RCTParameters.h.

References eicIsolationThreshold_.

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

75 { return eicIsolationThreshold_; }

L1RCTParameters::eicIsolationThreshold_

unsigned eicIsolationThreshold_

Definition: L1RCTParameters.h:156

◆ eMaxForFGCut()

double L1RCTParameters::eMaxForFGCut ( ) const

inline

Definition at line 68 of file L1RCTParameters.h.

References eMaxForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

68 { return eMaxForFGCut_; }

L1RCTParameters::eMaxForFGCut_

double eMaxForFGCut_

Definition: L1RCTParameters.h:121

◆ eMaxForHoECut()

double L1RCTParameters::eMaxForHoECut ( ) const

inline

Definition at line 71 of file L1RCTParameters.h.

References eMaxForHoECut_.

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

71 { return eMaxForHoECut_; }

L1RCTParameters::eMaxForHoECut_

double eMaxForHoECut_

Definition: L1RCTParameters.h:133

◆ eMinForFGCut()

double L1RCTParameters::eMinForFGCut ( ) const

inline

Definition at line 67 of file L1RCTParameters.h.

References eMinForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

67 { return eMinForFGCut_; }

L1RCTParameters::eMinForFGCut_

double eMinForFGCut_

Definition: L1RCTParameters.h:117

◆ eMinForHoECut()

double L1RCTParameters::eMinForHoECut ( ) const

inline

Definition at line 70 of file L1RCTParameters.h.

References eMinForHoECut_.

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

70 { return eMinForHoECut_; }

L1RCTParameters::eMinForHoECut_

double eMinForHoECut_

Definition: L1RCTParameters.h:129

◆ hActivityCut()

double L1RCTParameters::hActivityCut ( ) const

inline

Definition at line 74 of file L1RCTParameters.h.

References hActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

74 { return hActivityCut_; }

L1RCTParameters::hActivityCut_

double hActivityCut_

Definition: L1RCTParameters.h:149

◆ hMinForHoECut()

double L1RCTParameters::hMinForHoECut ( ) const

inline

Definition at line 72 of file L1RCTParameters.h.

References hMinForHoECut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

72 { return hMinForHoECut_; }

L1RCTParameters::hMinForHoECut_

double hMinForHoECut_

Definition: L1RCTParameters.h:137

◆ hOeCut()

double L1RCTParameters::hOeCut ( ) const

inline

Definition at line 69 of file L1RCTParameters.h.

References hOeCut_.

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

69 { return hOeCut_; }

L1RCTParameters::hOeCut_

double hOeCut_

Definition: L1RCTParameters.h:125

◆ jetMETECalScaleFactors()

const std::vector<double>& L1RCTParameters::jetMETECalScaleFactors ( ) const

inline

Definition at line 83 of file L1RCTParameters.h.

References jetMETECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

83 { return jetMETECalScaleFactors_; }

L1RCTParameters::jetMETECalScaleFactors_

std::vector< double > jetMETECalScaleFactors_

Definition: L1RCTParameters.h:198

◆ jetMETHCalScaleFactors()

const std::vector<double>& L1RCTParameters::jetMETHCalScaleFactors ( ) const

inline

Definition at line 84 of file L1RCTParameters.h.

References jetMETHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

84 { return jetMETHCalScaleFactors_; }

L1RCTParameters::jetMETHCalScaleFactors_

std::vector< double > jetMETHCalScaleFactors_

Definition: L1RCTParameters.h:199

◆ jetMETLSB()

double L1RCTParameters::jetMETLSB ( ) const

inline

Definition at line 66 of file L1RCTParameters.h.

References jetMETLSB_.

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

66 { return jetMETLSB_; }

L1RCTParameters::jetMETLSB_

double jetMETLSB_

Definition: L1RCTParameters.h:113

◆ JetMETTPGSum()

float L1RCTParameters::JetMETTPGSum	(	const float &	ecal,
		const float &	hcal,
		const unsigned &	iAbsEta
	)		const

Definition at line 189 of file L1RCTParameters.cc.

References correctedTPGSum(), eGammaECalScaleFactors_, Exception, createfilelist::int, jetMETECalScaleFactors_, jetMETHCalScaleFactors_, mps_fire::result, and useCorrections_.

Referenced by L1RCTLookupTables::jetMETETCode().

                                                                                                        {
   // We never deal with HF in this function (see note below)
   if (iAbsEta < 1 || iAbsEta > 28)
     throw cms::Exception("L1RCTParameters invalid function call") << "Eta out of range in MET TPGSum: " << iAbsEta;
   float ecal_c = ecal * jetMETECalScaleFactors_.at(iAbsEta - 1);
   float hcal_c = hcal * jetMETHCalScaleFactors_.at(iAbsEta - 1);
 
   // scale factors will either be length 28 for legacy, or 28*(# et bins + 1) where the first set is an average over the et bins
   // The first set provides a legacy fallthrough option
   // Currently, # et bins is 9
   if (jetMETECalScaleFactors_.size() == 28 * 10) {
     int et_bin = ((int)floor(ecal) / 5);
     // lowest bin (1) is 0-10GeV
     if (et_bin < 1)
       et_bin = 1;
     // highest bin (9) is 45GeV and up
     if (et_bin > 9)
       et_bin = 9;
     ecal_c = ecal * jetMETECalScaleFactors_.at(et_bin * 28 + iAbsEta - 1);
   }
 
   // We may be interested in HF jets, in which case, there are four more scale factors
   // for the 4 HF regions.  HOWEVER, we will never expect to see them accessed from
   // this function.  HF scaling is done in L1Trigger/RegionalCaloTrigger/src/L1RCTLookupTables
   if (jetMETHCalScaleFactors_.size() == 32 * 10) {
     int ht_bin = ((int)floor(hcal) / 5);
     // lowest bin (1) is 0-10GeV
     if (ht_bin < 1)
       ht_bin = 1;
     // highest bin (9) is 45GeV and up
     if (ht_bin > 9)
       ht_bin = 9;
     hcal_c = hcal * jetMETHCalScaleFactors_.at(ht_bin * 32 + iAbsEta - 1);
   }
 
   float result = ecal_c + hcal_c;
 
   // defunct section (polynomial-parameterized corrections)
   if (useCorrections_) {
     if (jetMETHCalScaleFactors_.at(iAbsEta - 1) != 0)
       hcal_c = hcal;
 
     if (jetMETECalScaleFactors_.at(iAbsEta - 1) != 0)
       ecal_c = ecal * eGammaECalScaleFactors_.at(iAbsEta - 1);  // Use eGamma Corrections
 
     result = correctedTPGSum(ecal_c, hcal_c, iAbsEta - 1);
   }
 
   return result;
 }

◆ jscQuietThresholdBarrel()

unsigned L1RCTParameters::jscQuietThresholdBarrel ( ) const

inline

Definition at line 76 of file L1RCTParameters.h.

References jscQuietThresholdBarrel_.

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

76 { return jscQuietThresholdBarrel_; }

L1RCTParameters::jscQuietThresholdBarrel_

unsigned jscQuietThresholdBarrel_

Definition: L1RCTParameters.h:161

◆ jscQuietThresholdEndcap()

unsigned L1RCTParameters::jscQuietThresholdEndcap ( ) const

inline

Definition at line 77 of file L1RCTParameters.h.

References jscQuietThresholdEndcap_.

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

77 { return jscQuietThresholdEndcap_; }

L1RCTParameters::jscQuietThresholdEndcap_

unsigned jscQuietThresholdEndcap_

Definition: L1RCTParameters.h:166

◆ noiseVetoHB()

bool L1RCTParameters::noiseVetoHB ( ) const

inline

Definition at line 78 of file L1RCTParameters.h.

References noiseVetoHB_.

Referenced by L1RCTLookupTables::lookup().

78 { return noiseVetoHB_; }

L1RCTParameters::noiseVetoHB_

bool noiseVetoHB_

Definition: L1RCTParameters.h:172

◆ noiseVetoHEminus()

bool L1RCTParameters::noiseVetoHEminus ( ) const

inline

Definition at line 80 of file L1RCTParameters.h.

References noiseVetoHEminus_.

Referenced by L1RCTLookupTables::lookup().

80 { return noiseVetoHEminus_; }

L1RCTParameters::noiseVetoHEminus_

bool noiseVetoHEminus_

Definition: L1RCTParameters.h:182

◆ noiseVetoHEplus()

bool L1RCTParameters::noiseVetoHEplus ( ) const

inline

Definition at line 79 of file L1RCTParameters.h.

References noiseVetoHEplus_.

Referenced by L1RCTLookupTables::lookup().

79 { return noiseVetoHEplus_; }

L1RCTParameters::noiseVetoHEplus_

bool noiseVetoHEplus_

Definition: L1RCTParameters.h:177

◆ print()

void L1RCTParameters::print ( std::ostream & s ) const

Definition at line 324 of file L1RCTParameters.cc.

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

                                                {
   s << "\nPrinting record L1RCTParametersRcd" << endl;
   s << "\n\"Parameter description\" \n  \"Parameter name\"  \"Value\" " << endl;
   s << "\ne/gamma least significant bit energy transmitted from receiver cards to EIC cards. \n  "
     << "eGammaLSB = " << eGammaLSB_ << endl;
   s << "\nLSB of region Et scale from RCT to GCT (GeV) \n  "
     << "jetMETLSB = " << jetMETLSB_ << endl;
   s << "\nminimum ECAL Et for which fine-grain veto is applied (GeV) \n  "
     << " eMinForFGCut = " << eMinForFGCut_ << endl;
   s << "\nmaximum ECAL Et for which fine-grain veto is applied (GeV) \n  "
     << "eMaxForFGCut = " << eMaxForFGCut_ << endl;
   s << "\nmaximum value of (HCAL Et / ECAL Et) \n  "
     << "hOeCut = " << hOeCut_ << endl;
   s << "\nminimum ECAL Et for which H/E veto is applied (GeV) \n  "
     << "eMinForHoECut = " << eMinForHoECut_ << endl;
   s << "\nmaximum ECAL Et for which H/E veto is applied (GeV) \n  "
     << "eMaxForHoECut = " << eMaxForHoECut_ << endl;
   s << "\nminimum HCAL Et for which H/E veto is applied (GeV)  \n  "
     << "hMinForHoECut = " << hMinForHoECut_ << endl;
   s << "\nECAL Et threshold above which tau activity bit is set (GeV)  \n  "
     << "eActivityCut = " << eActivityCut_ << endl;
   s << "\nHCAL Et threshold above which tau activity bit is set (GeV)  \n  "
     << "hActivityCut = " << hActivityCut_ << endl;
   s << "\nNeighboring trigger tower energy minimum threshold that marks candidate as non-isolated. (LSB bits) \n  "
     << "eicIsolationThreshold = " << eicIsolationThreshold_ << endl;
   s << "\nIf jetMet energy in RCT Barrel Region is below this value, a quiet bit is set. (LSB bits)\n  "
     << "jscQuietThreshBarrel = " << jscQuietThresholdBarrel_ << endl;
   s << "\nIf jetMet energy in RCT Endcap Region is below this value, a quiet bit is set. (LSB bits) \n  "
     << "jscQuietThreshEndcap = " << jscQuietThresholdEndcap_ << endl;
   s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT "
        "Barrel \n  "
     << "noiseVetoHB = " << noiseVetoHB_ << endl;
   s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT "
        "Encap+ \n  "
     << "noiseVetoHEplus = " << noiseVetoHEplus_ << endl;
   s << "\nWhen set to TRUE, HCAL energy is ignored if no ECAL energy is present in corresponding trigger tower for RCT "
        "Endcap- \n  "
     << "noiseVetoHEminus = " << noiseVetoHEminus_ << endl;
 
   auto printScalefactors = [&s](const std::vector<double>& sf) {
     if (sf.size() == 10 * 28) {
       s << "et bin  ieta  ScaleFactor" << endl;
       for (unsigned i = 0; i < sf.size(); i++)
         s << setw(6) << i / 28 << "  " << setw(4) << i % 28 + 1 << "  " << sf.at(i) << endl;
     } else if (sf.size() == 10 * 32)  // jet HCAL (HF regions are 29-32)
     {
       s << "et bin  ieta  ScaleFactor" << endl;
       for (unsigned i = 0; i < sf.size(); i++)
         s << setw(6) << i / 32 << "  " << setw(4) << i % 32 + 1 << "  " << sf.at(i) << endl;
     } else {
       s << "ieta  ScaleFactor" << endl;
       for (unsigned i = 0; i < sf.size(); i++)
         s << setw(4) << i + 1 << "  " << sf.at(i) << endl;
     }
   };
 
   s << "\n\neta-dependent multiplicative factors for ECAL Et before summation \n  "
     << "eGammaECal Scale Factors " << endl;
   printScalefactors(eGammaECalScaleFactors_);
 
   s << "\n\neta-dependent multiplicative factors for HCAL Et before summation \n  "
     << "eGammaHCal Scale Factors " << endl;
   printScalefactors(eGammaHCalScaleFactors_);
 
   s << "\n\neta-dependent multiplicative factors for ECAL Et before summation \n  "
     << "jetMETECal Scale Factors " << endl;
   printScalefactors(jetMETECalScaleFactors_);
 
   s << "\n\neta-dependent multiplicative factors for HCAL Et before summation \n"
     << "jetMETHCal Scale Factors " << endl;
   printScalefactors(jetMETHCalScaleFactors_);
 
   if (useCorrections_) {
     s << "\n\nUSING calibration variables " << endl;
 
     s << "\n\nH over E smear low Correction Factors " << endl;
     s << "ieta  Correction Factor" << endl;
     for (int i = 0; i < 28; i++)
       s << setw(4) << i + 1 << "  " << HoverE_smear_low_.at(i) << endl;
 
     s << "\n\nH over E smear high Correction Factors " << endl;
     s << "ieta  Correction Factor" << endl;
     for (int i = 0; i < 28; i++)
       s << setw(4) << i + 1 << "  " << HoverE_smear_high_.at(i) << endl;
 
     s << "\n\necal calibrations " << endl;
     s << "ieta  CorrFactor1  CorrFactor2  CorrFactor3" << endl;
     int end = ecal_calib_[0].size();
     for (int i = 0; i < 28; i++) {
       s << setw(4) << i;
       for (int j = 0; j < end; j++)
         s << setw(11) << setprecision(8) << ecal_calib_[i][j];
 
       s << endl;
     }
 
     s << "\n\nhcal calibrations " << endl;
     s << "ieta  CorrFactor1  CorrFactor2  CorrFactor3" << endl;
     end = hcal_calib_[0].size();
     for (int i = 0; i < 28; i++) {
       s << setw(4) << i;
       for (int j = 0; j < end; j++)
         s << setw(11) << setprecision(8) << hcal_calib_[i][j];
 
       s << endl;
     }
     s << "\n\nhcal_high calibrations " << endl;
     s << "ieta  CorrFactor1  CorrFactor2  CorrFactor3" << endl;
     end = hcal_high_calib_[0].size();
     for (int i = 0; i < 28; i++) {
       s << setw(4) << i;
       for (int j = 0; j < end; j++)
         s << setw(11) << setprecision(8) << hcal_high_calib_[i][j];
 
       s << endl;
     }
     end = cross_terms_[0].size();
     s << "\n\ncross terms calibrations " << endl;
     s << "ieta  CorrFactor1  CorrFactor2  CorrFactor3  CorrFactor4  CorrFactor5  CorrFactor6" << endl;
     for (int i = 0; i < 28; i++) {
       s << setw(4) << i;
       for (int j = 0; j < end; j++)
         s << setw(11) << setprecision(8) << cross_terms_[i][j];
 
       s << endl;
     }
 
   } else
     s << "\n\nNOT USING calibration variables " << endl;
 
   s << "\n\n" << endl;
 }

◆ serialize()

template<class Archive >

void L1RCTParameters::serialize	(	Archive &	ar,
		const unsigned int	version
	)

private

Friends And Related Function Documentation

◆ boost::serialization::access

friend class boost::serialization::access

friend

Definition at line 218 of file L1RCTParameters.h.

◆ cond::serialization::access

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 218 of file L1RCTParameters.h.

Member Data Documentation

◆ cross_terms_

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

private

Definition at line 211 of file L1RCTParameters.h.

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

◆ eActivityCut_

double L1RCTParameters::eActivityCut_

private

Definition at line 143 of file L1RCTParameters.h.

Referenced by eActivityCut(), and print().

◆ ecal_calib_

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

private

Definition at line 208 of file L1RCTParameters.h.

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

◆ eGammaECalScaleFactors_

std::vector<double> L1RCTParameters::eGammaECalScaleFactors_

private

Definition at line 191 of file L1RCTParameters.h.

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

◆ eGammaHCalScaleFactors_

std::vector<double> L1RCTParameters::eGammaHCalScaleFactors_

private

Definition at line 192 of file L1RCTParameters.h.

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

◆ eGammaLSB_

double L1RCTParameters::eGammaLSB_

private

Definition at line 109 of file L1RCTParameters.h.

Referenced by eGammaLSB(), and print().

◆ eicIsolationThreshold_

unsigned L1RCTParameters::eicIsolationThreshold_

private

Definition at line 156 of file L1RCTParameters.h.

Referenced by eicIsolationThreshold(), and print().

◆ eMaxForFGCut_

double L1RCTParameters::eMaxForFGCut_

private

Definition at line 121 of file L1RCTParameters.h.

Referenced by eMaxForFGCut(), and print().

◆ eMaxForHoECut_

double L1RCTParameters::eMaxForHoECut_

private

Definition at line 133 of file L1RCTParameters.h.

Referenced by eMaxForHoECut(), and print().

◆ eMinForFGCut_

double L1RCTParameters::eMinForFGCut_

private

Definition at line 117 of file L1RCTParameters.h.

Referenced by eMinForFGCut(), and print().

◆ eMinForHoECut_

double L1RCTParameters::eMinForHoECut_

private

Definition at line 129 of file L1RCTParameters.h.

Referenced by eMinForHoECut(), and print().

◆ hActivityCut_

double L1RCTParameters::hActivityCut_

private

Definition at line 149 of file L1RCTParameters.h.

Referenced by hActivityCut(), and print().

◆ hcal_calib_

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

private

Definition at line 209 of file L1RCTParameters.h.

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

◆ hcal_high_calib_

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

private

Definition at line 210 of file L1RCTParameters.h.

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

◆ hMinForHoECut_

double L1RCTParameters::hMinForHoECut_

private

Definition at line 137 of file L1RCTParameters.h.

Referenced by hMinForHoECut(), and print().

◆ hOeCut_

double L1RCTParameters::hOeCut_

private

Definition at line 125 of file L1RCTParameters.h.

Referenced by hOeCut(), and print().

◆ HoverE_smear_high_

std::vector<double> L1RCTParameters::HoverE_smear_high_

private

Definition at line 216 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), and print().

◆ HoverE_smear_low_

std::vector<double> L1RCTParameters::HoverE_smear_low_

private

Definition at line 215 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), and print().

◆ jetMETECalScaleFactors_

std::vector<double> L1RCTParameters::jetMETECalScaleFactors_

private

Definition at line 198 of file L1RCTParameters.h.

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

◆ jetMETHCalScaleFactors_

std::vector<double> L1RCTParameters::jetMETHCalScaleFactors_

private

Definition at line 199 of file L1RCTParameters.h.

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

◆ jetMETLSB_

double L1RCTParameters::jetMETLSB_

private

Definition at line 113 of file L1RCTParameters.h.

Referenced by jetMETLSB(), and print().

◆ jscQuietThresholdBarrel_

unsigned L1RCTParameters::jscQuietThresholdBarrel_

private

Definition at line 161 of file L1RCTParameters.h.

Referenced by jscQuietThresholdBarrel(), and print().

◆ jscQuietThresholdEndcap_

unsigned L1RCTParameters::jscQuietThresholdEndcap_

private

Definition at line 166 of file L1RCTParameters.h.

Referenced by jscQuietThresholdEndcap(), and print().

◆ noiseVetoHB_

bool L1RCTParameters::noiseVetoHB_

private

Definition at line 172 of file L1RCTParameters.h.

Referenced by noiseVetoHB(), and print().

◆ noiseVetoHEminus_

bool L1RCTParameters::noiseVetoHEminus_

private

Definition at line 182 of file L1RCTParameters.h.

Referenced by noiseVetoHEminus(), and print().

◆ noiseVetoHEplus_

bool L1RCTParameters::noiseVetoHEplus_

private

Definition at line 177 of file L1RCTParameters.h.

Referenced by noiseVetoHEplus(), and print().

◆ useCorrections_

bool L1RCTParameters::useCorrections_

private

Definition at line 185 of file L1RCTParameters.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ L1RCTParameters() [1/2]

◆ L1RCTParameters() [2/2]

◆ ~L1RCTParameters()

Member Function Documentation

◆ calcCard()

◆ calcCrate()

◆ calcIAbsEta()

◆ calcIEta()

◆ calcIPhi()

◆ calcTower()

◆ correctedTPGSum()

◆ eActivityCut()

◆ eGammaECalScaleFactors()

◆ eGammaHCalScaleFactors()

◆ eGammaLSB()

◆ EGammaTPGSum()

◆ eicIsolationThreshold()

◆ eMaxForFGCut()

◆ eMaxForHoECut()

◆ eMinForFGCut()

◆ eMinForHoECut()

◆ hActivityCut()

◆ hMinForHoECut()

◆ hOeCut()

◆ jetMETECalScaleFactors()

◆ jetMETHCalScaleFactors()

◆ jetMETLSB()

◆ JetMETTPGSum()

◆ jscQuietThresholdBarrel()

◆ jscQuietThresholdEndcap()

◆ noiseVetoHB()

◆ noiseVetoHEminus()

◆ noiseVetoHEplus()

◆ print()

◆ serialize()

Friends And Related Function Documentation

◆ boost::serialization::access

◆ cond::serialization::access

Member Data Documentation

◆ cross_terms_

◆ eActivityCut_

◆ ecal_calib_

◆ eGammaECalScaleFactors_

◆ eGammaHCalScaleFactors_

◆ eGammaLSB_

◆ eicIsolationThreshold_

◆ eMaxForFGCut_

◆ eMaxForHoECut_

◆ eMinForFGCut_

◆ eMinForHoECut_

◆ hActivityCut_

◆ hcal_calib_

◆ hcal_high_calib_

◆ hMinForHoECut_

◆ hOeCut_

◆ HoverE_smear_high_

◆ HoverE_smear_low_

◆ jetMETECalScaleFactors_

◆ jetMETHCalScaleFactors_

◆ jetMETLSB_

◆ jscQuietThresholdBarrel_

◆ jscQuietThresholdEndcap_

◆ noiseVetoHB_

◆ noiseVetoHEminus_

◆ noiseVetoHEplus_

◆ useCorrections_