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 28 of file L1RCTParameters.h.

Constructor & Destructor Documentation

L1RCTParameters::L1RCTParameters ( )

inline

Definition at line 33 of file L1RCTParameters.h.

References eActivityCut(), eGammaECalScaleFactors(), eGammaHCalScaleFactors(), eGammaLSB(), eicIsolationThreshold(), eMaxForFGCut(), eMaxForHoECut(), eMinForFGCut(), eMinForHoECut(), hActivityCut(), hMinForHoECut(), hOeCut(), jetMETECalScaleFactors(), jetMETHCalScaleFactors(), jetMETLSB(), jscQuietThresholdBarrel(), jscQuietThresholdEndcap(), noiseVetoHB(), noiseVetoHEminus(), and noiseVetoHEplus().

{}

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 mps_fire::i.

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

inline

Definition at line 65 of file L1RCTParameters.h.

{;}

Member Function Documentation

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

Definition at line 97 of file L1RCTParameters.cc.

Referenced by L1RCT::digiInput(), and jetMETHCalScaleFactors().

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

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

Definition at line 86 of file L1RCTParameters.cc.

References funct::abs().

Referenced by L1RCT::digiInput(), jetMETHCalScaleFactors(), 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;
}

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

Definition at line 182 of file L1RCTParameters.cc.

Referenced by calcIEta(), and jetMETHCalScaleFactors().

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

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

Definition at line 149 of file L1RCTParameters.cc.

References calcIAbsEta().

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

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

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

Definition at line 160 of file L1RCTParameters.cc.

Referenced by jetMETHCalScaleFactors(), 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;
}

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

Definition at line 114 of file L1RCTParameters.cc.

Referenced by L1RCT::digiInput(), jetMETHCalScaleFactors(), 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;
}

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

private

Definition at line 298 of file L1RCTParameters.cc.

References EnergyCorrector::c, cross_terms_, MillePedeFileConverter_cfg::e, RecoEcal_EventContent_cff::ec, digitizers_cfi::ecal, ecal_calib_, h, AnalysisDataFormats_SUSYBSMObjects::hc, digitizers_cfi::hcal, hcal_calib_, hcal_high_calib_, HoverE_smear_high_, HoverE_smear_low_, diffTreeTool::index, and funct::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;
}

double L1RCTParameters::eActivityCut ( ) const

inline

Definition at line 77 of file L1RCTParameters.h.

References eActivityCut_.

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

{return eActivityCut_;}

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

inline

Definition at line 85 of file L1RCTParameters.h.

References eGammaECalScaleFactors_.

Referenced by L1RCTParameters(), and L1RCTLookupTables::lookup().

{return eGammaECalScaleFactors_;}

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

inline

Definition at line 86 of file L1RCTParameters.h.

References eGammaHCalScaleFactors_.

Referenced by L1RCTParameters(), and L1RCTLookupTables::lookup().

{return eGammaHCalScaleFactors_;}

double L1RCTParameters::eGammaLSB ( ) const

inline

Definition at line 69 of file L1RCTParameters.h.

References eGammaLSB_.

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

{return eGammaLSB_;}

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

Definition at line 253 of file L1RCTParameters.cc.

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

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

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

unsigned L1RCTParameters::eicIsolationThreshold ( ) const

inline

Definition at line 79 of file L1RCTParameters.h.

References eicIsolationThreshold_.

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

{return eicIsolationThreshold_;}

double L1RCTParameters::eMaxForFGCut ( ) const

inline

Definition at line 72 of file L1RCTParameters.h.

References eMaxForFGCut_.

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

{return eMaxForFGCut_;}

double L1RCTParameters::eMaxForHoECut ( ) const

inline

Definition at line 75 of file L1RCTParameters.h.

References eMaxForHoECut_.

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

{return eMaxForHoECut_;}

double L1RCTParameters::eMinForFGCut ( ) const

inline

Definition at line 71 of file L1RCTParameters.h.

References eMinForFGCut_.

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

{return eMinForFGCut_;}

double L1RCTParameters::eMinForHoECut ( ) const

inline

Definition at line 74 of file L1RCTParameters.h.

References eMinForHoECut_.

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

{return eMinForHoECut_;}

double L1RCTParameters::hActivityCut ( ) const

inline

Definition at line 78 of file L1RCTParameters.h.

References hActivityCut_.

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

{return hActivityCut_;}

double L1RCTParameters::hMinForHoECut ( ) const

inline

Definition at line 76 of file L1RCTParameters.h.

References hMinForHoECut_.

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

{return hMinForHoECut_;}

double L1RCTParameters::hOeCut ( ) const

inline

Definition at line 73 of file L1RCTParameters.h.

References hOeCut_.

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

{return hOeCut_;}

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

inline

Definition at line 87 of file L1RCTParameters.h.

References jetMETECalScaleFactors_.

Referenced by L1RCTParameters(), and L1RCTLookupTables::lookup().

{return jetMETECalScaleFactors_;}

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

inline

Definition at line 88 of file L1RCTParameters.h.

References calcCard(), calcCrate(), calcIAbsEta(), calcIEta(), calcIPhi(), calcTower(), digitizers_cfi::ecal, EGammaTPGSum(), jetMETHCalScaleFactors_, JetMETTPGSum(), print(), and alignCSCRings::s.

Referenced by L1RCTParameters(), and L1RCTLookupTables::lookup().

{return jetMETHCalScaleFactors_;}

double L1RCTParameters::jetMETLSB ( ) const

inline

Definition at line 70 of file L1RCTParameters.h.

References jetMETLSB_.

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

{return jetMETLSB_;}

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

Definition at line 202 of file L1RCTParameters.cc.

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

Referenced by L1RCTLookupTables::jetMETETCode(), and jetMETHCalScaleFactors().

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

unsigned L1RCTParameters::jscQuietThresholdBarrel ( ) const

inline

Definition at line 80 of file L1RCTParameters.h.

References jscQuietThresholdBarrel_.

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

{return jscQuietThresholdBarrel_;}

unsigned L1RCTParameters::jscQuietThresholdEndcap ( ) const

inline

Definition at line 81 of file L1RCTParameters.h.

References jscQuietThresholdEndcap_.

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

{return jscQuietThresholdEndcap_;}

bool L1RCTParameters::noiseVetoHB ( ) const

inline

Definition at line 82 of file L1RCTParameters.h.

References noiseVetoHB_.

Referenced by L1RCTParameters(), and L1RCTLookupTables::lookup().

{return noiseVetoHB_;}

bool L1RCTParameters::noiseVetoHEminus ( ) const

inline

Definition at line 84 of file L1RCTParameters.h.

References noiseVetoHEminus_.

Referenced by L1RCTParameters(), and L1RCTLookupTables::lookup().

{return noiseVetoHEminus_;}

bool L1RCTParameters::noiseVetoHEplus ( ) const

inline

Definition at line 83 of file L1RCTParameters.h.

References noiseVetoHEplus_.

Referenced by L1RCTParameters(), and L1RCTLookupTables::lookup().

{return noiseVetoHEplus_;}

void L1RCTParameters::print

(

std::ostream &

s

)

const

Definition at line 349 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_, mps_fire::i, jetMETECalScaleFactors_, jetMETHCalScaleFactors_, jetMETLSB_, jscQuietThresholdBarrel_, jscQuietThresholdEndcap_, noiseVetoHB_, noiseVetoHEminus_, noiseVetoHEplus_, alignCSCRings::s, and useCorrections_.

Referenced by L1RCTParametersTester::analyze(), L1TConfigDumper::analyze(), and jetMETHCalScaleFactors().

                                           {

    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;

}

template<class Archive >

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

private

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 226 of file L1RCTParameters.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 226 of file L1RCTParameters.h.

Member Data Documentation

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

private

Definition at line 217 of file L1RCTParameters.h.

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

double L1RCTParameters::eActivityCut_

private

Definition at line 148 of file L1RCTParameters.h.

Referenced by eActivityCut(), and print().

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

private

Definition at line 214 of file L1RCTParameters.h.

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

std::vector<double> L1RCTParameters::eGammaECalScaleFactors_

private

Definition at line 196 of file L1RCTParameters.h.

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

std::vector<double> L1RCTParameters::eGammaHCalScaleFactors_

private

Definition at line 197 of file L1RCTParameters.h.

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

double L1RCTParameters::eGammaLSB_

private

Definition at line 114 of file L1RCTParameters.h.

Referenced by eGammaLSB(), and print().

unsigned L1RCTParameters::eicIsolationThreshold_

private

Definition at line 161 of file L1RCTParameters.h.

Referenced by eicIsolationThreshold(), and print().

double L1RCTParameters::eMaxForFGCut_

private

Definition at line 126 of file L1RCTParameters.h.

Referenced by eMaxForFGCut(), and print().

double L1RCTParameters::eMaxForHoECut_

private

Definition at line 138 of file L1RCTParameters.h.

Referenced by eMaxForHoECut(), and print().

double L1RCTParameters::eMinForFGCut_

private

Definition at line 122 of file L1RCTParameters.h.

Referenced by eMinForFGCut(), and print().

double L1RCTParameters::eMinForHoECut_

private

Definition at line 134 of file L1RCTParameters.h.

Referenced by eMinForHoECut(), and print().

double L1RCTParameters::hActivityCut_

private

Definition at line 154 of file L1RCTParameters.h.

Referenced by hActivityCut(), and print().

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

private

Definition at line 215 of file L1RCTParameters.h.

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

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

private

Definition at line 216 of file L1RCTParameters.h.

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

double L1RCTParameters::hMinForHoECut_

private

Definition at line 142 of file L1RCTParameters.h.

Referenced by hMinForHoECut(), and print().

double L1RCTParameters::hOeCut_

private

Definition at line 130 of file L1RCTParameters.h.

Referenced by hOeCut(), and print().

std::vector<double> L1RCTParameters::HoverE_smear_high_

private

Definition at line 222 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), and print().

std::vector<double> L1RCTParameters::HoverE_smear_low_

private

Definition at line 221 of file L1RCTParameters.h.

Referenced by correctedTPGSum(), and print().

std::vector<double> L1RCTParameters::jetMETECalScaleFactors_

private

Definition at line 203 of file L1RCTParameters.h.

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

std::vector<double> L1RCTParameters::jetMETHCalScaleFactors_

private

Definition at line 204 of file L1RCTParameters.h.

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

double L1RCTParameters::jetMETLSB_

private

Definition at line 118 of file L1RCTParameters.h.

Referenced by jetMETLSB(), and print().

unsigned L1RCTParameters::jscQuietThresholdBarrel_

private

Definition at line 166 of file L1RCTParameters.h.

Referenced by jscQuietThresholdBarrel(), and print().

unsigned L1RCTParameters::jscQuietThresholdEndcap_

private

Definition at line 171 of file L1RCTParameters.h.

Referenced by jscQuietThresholdEndcap(), and print().

bool L1RCTParameters::noiseVetoHB_

private

Definition at line 177 of file L1RCTParameters.h.

Referenced by noiseVetoHB(), and print().

bool L1RCTParameters::noiseVetoHEminus_

private

Definition at line 187 of file L1RCTParameters.h.

Referenced by noiseVetoHEminus(), and print().

bool L1RCTParameters::noiseVetoHEplus_

private

Definition at line 182 of file L1RCTParameters.h.

Referenced by noiseVetoHEplus(), and print().

bool L1RCTParameters::useCorrections_

private

Definition at line 190 of file L1RCTParameters.h.

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