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
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.

{}

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.

    : 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]);
}

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

~L1RCTParameters()

L1RCTParameters::~L1RCTParameters ( )

inline

Definition at line 61 of file L1RCTParameters.h.

{ ; }

Member Function Documentation

calcCard()

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

Definition at line 94 of file L1RCTParameters.cc.

                                                                                              {
  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;
}

Referenced by L1RCT::digiInput().

calcCrate()

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

Definition at line 83 of file L1RCTParameters.cc.

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

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

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

calcIAbsEta()

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

Definition at line 171 of file L1RCTParameters.cc.

                                                                                                                    {
  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;
}

Referenced by calcIEta().

calcIEta()

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

Definition at line 141 of file L1RCTParameters.cc.

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

References calcIAbsEta(), and L1TowerCalibrationProducer_cfi::iEta.

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

calcIPhi()

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

Definition at line 152 of file L1RCTParameters.cc.

                                                                                                                 {
  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;
}

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

calcTower()

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

Definition at line 109 of file L1RCTParameters.cc.

                                                                                               {
  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;
}

References LEDCalibrationChannels::iphi, and hgcalTowerProducer_cfi::tower.

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

correctedTPGSum()

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

private

Definition at line 285 of file L1RCTParameters.cc.

                                                                                                        {
  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;
}

References 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 funct::pow().

Referenced by EGammaTPGSum(), and JetMETTPGSum().

eActivityCut()

double L1RCTParameters::eActivityCut ( ) const

inline

Definition at line 73 of file L1RCTParameters.h.

{ return eActivityCut_; }

References eActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

eGammaECalScaleFactors()

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

inline

Definition at line 81 of file L1RCTParameters.h.

{ return eGammaECalScaleFactors_; }

References eGammaECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

eGammaHCalScaleFactors()

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

inline

Definition at line 82 of file L1RCTParameters.h.

{ return eGammaHCalScaleFactors_; }

References eGammaHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

eGammaLSB()

double L1RCTParameters::eGammaLSB ( ) const

inline

Definition at line 65 of file L1RCTParameters.h.

{ return eGammaLSB_; }

References eGammaLSB_.

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

EGammaTPGSum()

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

Definition at line 240 of file L1RCTParameters.cc.

                                                                                                       {
  // 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;
}

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

Referenced by L1RCTLookupTables::eGammaETCode().

eicIsolationThreshold()

unsigned L1RCTParameters::eicIsolationThreshold ( ) const

inline

Definition at line 75 of file L1RCTParameters.h.

{ return eicIsolationThreshold_; }

References eicIsolationThreshold_.

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

eMaxForFGCut()

double L1RCTParameters::eMaxForFGCut ( ) const

inline

Definition at line 68 of file L1RCTParameters.h.

{ return eMaxForFGCut_; }

References eMaxForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

eMaxForHoECut()

double L1RCTParameters::eMaxForHoECut ( ) const

inline

Definition at line 71 of file L1RCTParameters.h.

{ return eMaxForHoECut_; }

References eMaxForHoECut_.

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

eMinForFGCut()

double L1RCTParameters::eMinForFGCut ( ) const

inline

Definition at line 67 of file L1RCTParameters.h.

{ return eMinForFGCut_; }

References eMinForFGCut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

eMinForHoECut()

double L1RCTParameters::eMinForHoECut ( ) const

inline

Definition at line 70 of file L1RCTParameters.h.

{ return eMinForHoECut_; }

References eMinForHoECut_.

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

hActivityCut()

double L1RCTParameters::hActivityCut ( ) const

inline

Definition at line 74 of file L1RCTParameters.h.

{ return hActivityCut_; }

References hActivityCut_.

Referenced by L1RCTLookupTables::activityBit().

hMinForHoECut()

double L1RCTParameters::hMinForHoECut ( ) const

inline

Definition at line 72 of file L1RCTParameters.h.

{ return hMinForHoECut_; }

References hMinForHoECut_.

Referenced by L1RCTLookupTables::hOeFGVetoBit().

hOeCut()

double L1RCTParameters::hOeCut ( ) const

inline

Definition at line 69 of file L1RCTParameters.h.

{ return hOeCut_; }

References hOeCut_.

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

jetMETECalScaleFactors()

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

inline

Definition at line 83 of file L1RCTParameters.h.

{ return jetMETECalScaleFactors_; }

References jetMETECalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

jetMETHCalScaleFactors()

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

inline

Definition at line 84 of file L1RCTParameters.h.

{ return jetMETHCalScaleFactors_; }

References jetMETHCalScaleFactors_.

Referenced by L1RCTLookupTables::lookup().

jetMETLSB()

double L1RCTParameters::jetMETLSB ( ) const

inline

Definition at line 66 of file L1RCTParameters.h.

{ return jetMETLSB_; }

References jetMETLSB_.

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

JetMETTPGSum()

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

Definition at line 189 of file L1RCTParameters.cc.

                                                                                                       {
  // 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;
}

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

Referenced by L1RCTLookupTables::jetMETETCode().

jscQuietThresholdBarrel()

unsigned L1RCTParameters::jscQuietThresholdBarrel ( ) const

inline

Definition at line 76 of file L1RCTParameters.h.

{ return jscQuietThresholdBarrel_; }

References jscQuietThresholdBarrel_.

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

jscQuietThresholdEndcap()

unsigned L1RCTParameters::jscQuietThresholdEndcap ( ) const

inline

Definition at line 77 of file L1RCTParameters.h.

{ return jscQuietThresholdEndcap_; }

References jscQuietThresholdEndcap_.

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

noiseVetoHB()

bool L1RCTParameters::noiseVetoHB ( ) const

inline

Definition at line 78 of file L1RCTParameters.h.

{ return noiseVetoHB_; }

References noiseVetoHB_.

Referenced by L1RCTLookupTables::lookup().

noiseVetoHEminus()

bool L1RCTParameters::noiseVetoHEminus ( ) const

inline

Definition at line 80 of file L1RCTParameters.h.

{ return noiseVetoHEminus_; }

References noiseVetoHEminus_.

Referenced by L1RCTLookupTables::lookup().

noiseVetoHEplus()

bool L1RCTParameters::noiseVetoHEplus ( ) const

inline

Definition at line 79 of file L1RCTParameters.h.

{ return noiseVetoHEplus_; }

References noiseVetoHEplus_.

Referenced by L1RCTLookupTables::lookup().

print()

void L1RCTParameters::print

(

std::ostream &

s

)

const

Definition at line 324 of file L1RCTParameters.cc.

                                               {
  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;
}

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_.

Referenced by L1RCTParametersTester::analyze().

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().