L1GctJetFinderParams Class Reference

#include <L1GctJetFinderParams.h>

Public Member Functions
uint16_t	correctedEtGct (const double correctedEt) const
	Convert the corrected Et value to a linear Et for Ht summing. More...
double	correctedEtGeV (const double et, const unsigned eta, const bool tauVeto) const
	Eta takes a value from 0-10, corresponding to jet regions running from eta=0.0 to eta=5.0. More...
unsigned	getCenForJetEtaBoundary () const
unsigned	getCenJetEtSeedGct () const
double	getCenJetEtSeedGeV () const
bool	getConvertToEnergy () const
unsigned	getCorrType () const
	Access to jet Et calibration parameters. More...
unsigned	getForJetEtSeedGct () const
double	getForJetEtSeedGeV () const
unsigned	getHtJetEtThresholdGct () const
double	getHtJetEtThresholdGeV () const
double	getHtLsbGeV () const
const std::vector< std::vector< double > > &	getJetCorrCoeffs () const
unsigned	getMHtJetEtThresholdGct () const
double	getMHtJetEtThresholdGeV () const
double	getRgnEtLsbGeV () const
const std::vector< std::vector< double > > &	getTauCorrCoeffs () const
unsigned	getTauIsoEtThresholdGct () const
double	getTauIsoEtThresholdGeV () const
unsigned	getTauJetEtSeedGct () const
double	getTauJetEtSeedGeV () const
	L1GctJetFinderParams ()
	L1GctJetFinderParams (double rgnEtLsb, double htLsb, double cJetSeed, double fJetSeed, double tJetSeed, double tauIsoEtThresh, double htJetEtThresh, double mhtJetEtThresh, unsigned etaBoundary, unsigned corrType, const std::vector< std::vector< double > > &jetCorrCoeffs, const std::vector< std::vector< double > > &tauCorrCoeffs, bool convertToEnergy, const std::vector< double > &energyConvCoeffs)
void	setHtSumParams (const double htLsb, const double htJetEtThresh, const double mhtJetEtThresh)
void	setJetEtCalibrationParams (const unsigned corrType, const std::vector< std::vector< double > > &jetCorrCoeffs, const std::vector< std::vector< double > > &tauCorrCoeffs)
void	setJetEtConvertToEnergyOff ()
void	setJetEtConvertToEnergyOn (const std::vector< double > &energyConvCoeffs)
void	setParams (const double rgnEtLsb, const double htLsb, const double cJetSeed, const double fJetSeed, const double tJetSeed, const double tauIsoEtThresh, const double htJetEtThresh, const double mhtJetEtThresh, const unsigned etaBoundary, const unsigned corrType, const std::vector< std::vector< double > > &jetCorrCoeffs, const std::vector< std::vector< double > > &tauCorrCoeffs)
void	setRegionEtLsb (const double rgnEtLsb)
void	setSlidingWindowParams (const double cJetSeed, const double fJetSeed, const double tJetSeed, const unsigned etaBoundary)
void	setTauAlgorithmParams (const double tauIsoEtThresh)
	~L1GctJetFinderParams ()

Static Public Attributes
static const unsigned	N_CENTRAL_ETA_VALUES = 7
	Number of eta bins used in correction. More...
static const unsigned	NUMBER_ETA_VALUES = 11
	Number of eta bins used in correction. More...

Private Member Functions
double	correctionFunction (const double Et, const std::vector< double > &coeffs) const
double	findCorrectedEt (const double Et, const std::vector< double > &coeffs) const
double	orcaStyleCorrect (const double Et, const std::vector< double > &coeffs) const
double	pfCorrect (const double Et, const std::vector< double > &coeffs) const
double	piecewiseCubicCorrect (const double Et, const std::vector< double > &coeffs) const
double	powerSeriesCorrect (const double Et, const std::vector< double > &coeffs) const
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)
double	simpleCorrect (const double Et, const std::vector< double > &coeffs) const

Private Attributes
unsigned	cenForJetEtaBoundary_
double	cenJetEtSeed_
bool	convertToEnergy_
unsigned	corrType_
std::vector< double >	energyConversionCoeffs_
double	forJetEtSeed_
double	htJetEtThreshold_
double	htLsb_
std::vector< std::vector< double > >	jetCorrCoeffs_
double	mhtJetEtThreshold_
double	rgnEtLsb_
std::vector< std::vector< double > >	tauCorrCoeffs_
double	tauIsoEtThreshold_
double	tauJetEtSeed_

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

Detailed Description

Definition at line 10 of file L1GctJetFinderParams.h.

Constructor & Destructor Documentation

L1GctJetFinderParams::L1GctJetFinderParams

(

)

Definition at line 17 of file L1GctJetFinderParams.cc.

                                           :
  rgnEtLsb_(0.),
  htLsb_(0.),
  cenJetEtSeed_(0.),
  forJetEtSeed_(0.),
  tauJetEtSeed_(0.),
  tauIsoEtThreshold_(0.),
  htJetEtThreshold_(0.),
  mhtJetEtThreshold_(0.),
  cenForJetEtaBoundary_(0),
  corrType_(0),
  jetCorrCoeffs_(),
  tauCorrCoeffs_(),
  convertToEnergy_(false),
  energyConversionCoeffs_()
{ }

L1GctJetFinderParams::L1GctJetFinderParams	(	double	rgnEtLsb,
		double	htLsb,
		double	cJetSeed,
		double	fJetSeed,
		double	tJetSeed,
		double	tauIsoEtThresh,
		double	htJetEtThresh,
		double	mhtJetEtThresh,
		unsigned	etaBoundary,
		unsigned	corrType,
		const std::vector< std::vector< double > > &	jetCorrCoeffs,
		const std::vector< std::vector< double > > &	tauCorrCoeffs,
		bool	convertToEnergy,
		const std::vector< double > &	energyConvCoeffs
	)

Definition at line 34 of file L1GctJetFinderParams.cc.

References corrType_, energyConversionCoeffs_, Exception, jetCorrCoeffs_, LogDebug, N_CENTRAL_ETA_VALUES, NUMBER_ETA_VALUES, and tauCorrCoeffs_.

                                                                  :
  rgnEtLsb_(rgnEtLsb),
  htLsb_(htLsb),
  cenJetEtSeed_(cJetSeed),
  forJetEtSeed_(fJetSeed),
  tauJetEtSeed_(tJetSeed),
  tauIsoEtThreshold_(tauIsoEtThresh),
  htJetEtThreshold_(htJetEtThresh),
  mhtJetEtThreshold_(mhtJetEtThresh),
  cenForJetEtaBoundary_(etaBoundary),
  corrType_(corrType),
  jetCorrCoeffs_(jetCorrCoeffs),
  tauCorrCoeffs_(tauCorrCoeffs),
  convertToEnergy_(convertToEnergy),
  energyConversionCoeffs_(energyConvCoeffs)
{ 
  // check number of eta bins
  if (jetCorrCoeffs_.size() != NUMBER_ETA_VALUES ||
      tauCorrCoeffs_.size() != N_CENTRAL_ETA_VALUES ||
      energyConversionCoeffs_.size() != NUMBER_ETA_VALUES) {

    LogDebug("L1-O2O") << "GCT jet corrections constructed with " << jetCorrCoeffs_.size() << " bins, expected " << NUMBER_ETA_VALUES << std::endl;
    LogDebug("L1-O2O") << "GCT tau corrections constructed with " << tauCorrCoeffs_.size() << " bins, expected " << N_CENTRAL_ETA_VALUES << std::endl;
    LogDebug("L1-O2O") << "GCT energy corrections constructed with " << energyConversionCoeffs_.size() << " bins, expected " << NUMBER_ETA_VALUES << std::endl;

    throw cms::Exception("InconsistentConfig") << "L1GctJetFinderParams constructed with wrong number of eta bins : " << jetCorrCoeffs_.size() << " jets, " << tauCorrCoeffs_.size() << " taus, " << energyConversionCoeffs_.size() << " energy conversion bins" << std::endl;

  }

  // check number of coefficients against expectation
  unsigned expCoeffs = 0;
  if (corrType_ == 2) expCoeffs=8;  // ORCA style
  if (corrType_ == 3) expCoeffs=4;  // Simple
  if (corrType_ == 4) expCoeffs=15;  // piecewise-cubic
  if (corrType_ == 5) expCoeffs=6;  // PF

  // correction types 1 and 4 can have any number of parameters
  if (corrType_ != 1 &&
      corrType_ != 4) {
    std::vector< std::vector<double> >::const_iterator itr;      
    for (itr=jetCorrCoeffs_.begin(); itr!=jetCorrCoeffs_.end(); ++itr) {
      if (itr->size() != expCoeffs) {
    throw cms::Exception("InconsistentConfig") << "L1GctJetFinderParams constructed with " << itr->size() << " jet correction coefficients, when " << expCoeffs << " expected" << std::endl;
      }
    }
    for (itr=tauCorrCoeffs_.begin(); itr!=tauCorrCoeffs_.end(); ++itr) {
      if (itr->size() != expCoeffs) {
    throw cms::Exception("InconsistentConfig") << "L1GctJetFinderParams constructed with " << itr->size() << " tau correction coefficients, when " << expCoeffs << " expected"<< std::endl;
      }
    }
  }
  
}

L1GctJetFinderParams::~L1GctJetFinderParams

(

)

Definition at line 102 of file L1GctJetFinderParams.cc.

{}

Member Function Documentation

uint16_t L1GctJetFinderParams::correctedEtGct

(

const double

correctedEt

)

const

Convert the corrected Et value to a linear Et for Ht summing.

Convert the corrected Et value to an integer Et for Ht summing.

Definition at line 209 of file L1GctJetFinderParams.cc.

References htLsb_, and L1GctStaticParameters::jetCalibratedEtMax.

Referenced by getMHtJetEtThresholdGct().

{
  double scaledEt = correctedEt / htLsb_;

  uint16_t jetEtOut = static_cast<uint16_t>(scaledEt);
  
  if(jetEtOut > L1GctStaticParameters::jetCalibratedEtMax) {
    return L1GctStaticParameters::jetCalibratedEtMax;
  } else {
    return jetEtOut;
  }
}

double L1GctJetFinderParams::correctedEtGeV	(	const double	et,
		const unsigned	eta,
		const bool	tauVeto
	)		const

Eta takes a value from 0-10, corresponding to jet regions running from eta=0.0 to eta=5.0.

Definition at line 185 of file L1GctJetFinderParams.cc.

References cenForJetEtaBoundary_, convertToEnergy_, correctionFunction(), energyConversionCoeffs_, jetCorrCoeffs_, NUMBER_ETA_VALUES, mps_fire::result, and tauCorrCoeffs_.

Referenced by getMHtJetEtThresholdGct(), and L1GctJetEtCalibrationLut::value().

                                               {

  if (eta>=NUMBER_ETA_VALUES) {
    return 0;
  } else {
    double result=0;
    if ((eta>=cenForJetEtaBoundary_) || tauVeto) {
      // Use jetCorrCoeffs for central and forward jets.
      // In forward eta bins we ignore the tau flag (as in the firmware)
      result=correctionFunction(et, jetCorrCoeffs_.at(eta));
    } else {
      // Use tauCorrCoeffs for tau jets (in central eta bins)
      result=correctionFunction(et, tauCorrCoeffs_.at(eta));
    }
    if (convertToEnergy_)  { result *= energyConversionCoeffs_.at(eta); }
    return result;
  }

}

double L1GctJetFinderParams::correctionFunction ( const double  Et, const std::vector< double > &  coeffs  ) const

private

Definition at line 226 of file L1GctJetFinderParams.cc.

References corrType_, orcaStyleCorrect(), pfCorrect(), piecewiseCubicCorrect(), powerSeriesCorrect(), mps_fire::result, and simpleCorrect().

Referenced by correctedEtGeV(), and getTauCorrCoeffs().

{
  double result=0;
  switch (corrType_)
  {
  case 0:  // no correction
    result = Et;
    break;
  case 1:   // power series correction
    result = powerSeriesCorrect(Et, coeffs);
    break;
  case 2:  // ORCA style correction
    result = orcaStyleCorrect(Et, coeffs);
    break;
  case 3:  // simple correction (JetMEt style)
    result = simpleCorrect(Et, coeffs);
    break;
  case 4:  // piecwise cubic correction a la Greg Landsberg et al
    result = piecewiseCubicCorrect(Et, coeffs);
    break;
  case 5:  // PF correction
    result = pfCorrect(Et, coeffs);
    break;
  default:
    result = Et;      
  }
  return result;
}

double L1GctJetFinderParams::findCorrectedEt ( const double  Et, const std::vector< double > &  coeffs  ) const

private

Referenced by getTauCorrCoeffs().

unsigned L1GctJetFinderParams::getCenForJetEtaBoundary ( ) const

inline

Definition at line 45 of file L1GctJetFinderParams.h.

References cenForJetEtaBoundary_.

Referenced by operator<<(), L1GlobalCaloTrigger::setJetFinderParams(), and L1GctJetFinderBase::setJetFinderParams().

{ return cenForJetEtaBoundary_; }

unsigned L1GctJetFinderParams::getCenJetEtSeedGct ( ) const

inline

Definition at line 49 of file L1GctJetFinderParams.h.

References cenJetEtSeed_, and rgnEtLsb_.

Referenced by L1GlobalCaloTrigger::setJetFinderParams(), and L1GctJetFinderBase::setJetFinderParams().

{ return static_cast<unsigned>(cenJetEtSeed_/rgnEtLsb_); }

double L1GctJetFinderParams::getCenJetEtSeedGeV ( ) const

inline

Definition at line 39 of file L1GctJetFinderParams.h.

References cenJetEtSeed_.

Referenced by operator<<().

{ return cenJetEtSeed_; }

bool L1GctJetFinderParams::getConvertToEnergy ( ) const

inline

Definition at line 46 of file L1GctJetFinderParams.h.

References convertToEnergy_.

{ return convertToEnergy_; }

unsigned L1GctJetFinderParams::getCorrType ( ) const

inline

Access to jet Et calibration parameters.

Definition at line 92 of file L1GctJetFinderParams.h.

References corrType_.

Referenced by operator<<().

{ return corrType_; }

unsigned L1GctJetFinderParams::getForJetEtSeedGct ( ) const

inline

Definition at line 50 of file L1GctJetFinderParams.h.

References forJetEtSeed_, and rgnEtLsb_.

Referenced by L1GctJetFinderBase::setJetFinderParams().

{ return static_cast<unsigned>(forJetEtSeed_/rgnEtLsb_); }

double L1GctJetFinderParams::getForJetEtSeedGeV ( ) const

inline

Definition at line 40 of file L1GctJetFinderParams.h.

References forJetEtSeed_.

Referenced by operator<<().

{ return forJetEtSeed_; }

unsigned L1GctJetFinderParams::getHtJetEtThresholdGct ( ) const

inline

Definition at line 53 of file L1GctJetFinderParams.h.

References htJetEtThreshold_, and htLsb_.

Referenced by L1GctJetFinderBase::setJetFinderParams().

{ return static_cast<unsigned>(htJetEtThreshold_/htLsb_); }

double L1GctJetFinderParams::getHtJetEtThresholdGeV ( ) const

inline

Definition at line 43 of file L1GctJetFinderParams.h.

References htJetEtThreshold_.

Referenced by L1GctValidation::analyze(), and operator<<().

{ return htJetEtThreshold_; }

double L1GctJetFinderParams::getHtLsbGeV ( ) const

inline

Definition at line 38 of file L1GctJetFinderParams.h.

References htLsb_.

Referenced by L1GctValidation::analyze(), L1GetHistLimits::getHistLimits(), operator<<(), L1ExtraParticlesProd::produce(), and L1GctGlobalEnergyAlgos::setJetFinderParams().

{ return htLsb_; }

const std::vector< std::vector<double> >& L1GctJetFinderParams::getJetCorrCoeffs ( ) const

inline

Definition at line 93 of file L1GctJetFinderParams.h.

References jetCorrCoeffs_.

Referenced by operator<<().

{ return jetCorrCoeffs_; }

unsigned L1GctJetFinderParams::getMHtJetEtThresholdGct ( ) const

inline

Definition at line 54 of file L1GctJetFinderParams.h.

References correctedEtGct(), correctedEtGeV(), stringResolutionProvider_cfi::et, stringResolutionProvider_cfi::eta, htLsb_, mhtJetEtThreshold_, setHtSumParams(), setJetEtCalibrationParams(), setJetEtConvertToEnergyOff(), setJetEtConvertToEnergyOn(), setParams(), setRegionEtLsb(), setSlidingWindowParams(), and setTauAlgorithmParams().

Referenced by L1GctJetFinderBase::setJetFinderParams().

{ return static_cast<unsigned>(mhtJetEtThreshold_/htLsb_); }

double L1GctJetFinderParams::getMHtJetEtThresholdGeV ( ) const

inline

Definition at line 44 of file L1GctJetFinderParams.h.

References mhtJetEtThreshold_.

Referenced by L1GctValidation::analyze(), and operator<<().

{ return mhtJetEtThreshold_; }

double L1GctJetFinderParams::getRgnEtLsbGeV ( ) const

inline

Definition at line 37 of file L1GctJetFinderParams.h.

References rgnEtLsb_.

Referenced by L1GctValidation::analyze(), and operator<<().

{ return rgnEtLsb_; }

const std::vector< std::vector<double> >& L1GctJetFinderParams::getTauCorrCoeffs ( ) const

inline

Definition at line 94 of file L1GctJetFinderParams.h.

References correctionFunction(), findCorrectedEt(), orcaStyleCorrect(), pfCorrect(), piecewiseCubicCorrect(), powerSeriesCorrect(), simpleCorrect(), and tauCorrCoeffs_.

Referenced by operator<<().

{ return tauCorrCoeffs_; }

unsigned L1GctJetFinderParams::getTauIsoEtThresholdGct ( ) const

inline

Definition at line 52 of file L1GctJetFinderParams.h.

References rgnEtLsb_, and tauIsoEtThreshold_.

Referenced by L1GctJetFinderBase::setJetFinderParams().

{ return static_cast<unsigned>(tauIsoEtThreshold_/rgnEtLsb_); }

double L1GctJetFinderParams::getTauIsoEtThresholdGeV ( ) const

inline

Definition at line 42 of file L1GctJetFinderParams.h.

References tauIsoEtThreshold_.

Referenced by operator<<().

{ return tauIsoEtThreshold_; }

unsigned L1GctJetFinderParams::getTauJetEtSeedGct ( ) const

inline

Definition at line 51 of file L1GctJetFinderParams.h.

References rgnEtLsb_, and tauJetEtSeed_.

Referenced by L1GlobalCaloTrigger::setJetFinderParams(), and L1GctJetFinderBase::setJetFinderParams().

{ return static_cast<unsigned>(tauJetEtSeed_/rgnEtLsb_); }

double L1GctJetFinderParams::getTauJetEtSeedGeV ( ) const

inline

Definition at line 41 of file L1GctJetFinderParams.h.

References tauJetEtSeed_.

Referenced by operator<<().

{ return tauJetEtSeed_; }

double L1GctJetFinderParams::orcaStyleCorrect ( const double  Et, const std::vector< double > &  coeffs  ) const

private

Definition at line 264 of file L1GctJetFinderParams.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, patCaloMETCorrections_cff::C, mathSSE::sqrt(), and ctppsDiamondLocalTracks_cfi::threshold.

Referenced by correctionFunction(), and getTauCorrCoeffs().

{

  // The coefficients are arranged in groups of four
  // The first in each group is a threshold value of Et.

  std::vector<double>::const_iterator next_coeff=coeffs.begin();
  while (next_coeff != coeffs.end()) {
    double threshold = *next_coeff++;
    double A = *next_coeff++;
    double B = *next_coeff++;
    double C = *next_coeff++;
    if (Et>threshold) {
      // This function is an inverse quadratic:
      //   (input Et) = A + B*(output Et) + C*(output Et)^2
      return 2*(Et-A)/(B+sqrt(B*B-4*A*C+4*Et*C));
    }
    // If we are below all specified thresholds (or the vector is empty), return output=input.
  }
  return Et;
}

double L1GctJetFinderParams::pfCorrect ( const double  Et, const std::vector< double > &  coeffs  ) const

private

Definition at line 329 of file L1GctJetFinderParams.cc.

References JetChargeProducer_cfi::exp, and funct::pow().

Referenced by correctionFunction(), and getTauCorrCoeffs().

{  
  //
  // corr_factor = [0]+[1]/(pow(log10(x),2)+[2])+[3]*exp(-[4]*(log10(x)-[5])*(log10(x)-[5])) 
  // Et_out = Et_in * corr_factor
  //

  return et * (coeffs.at(0)+coeffs.at(1)/(pow(log10(et),2)+coeffs.at(2))+coeffs.at(3)*exp(-coeffs.at(4)*(log10(et)-coeffs.at(5))*(log10(et)-coeffs.at(5))));

}

double L1GctJetFinderParams::piecewiseCubicCorrect ( const double  Et, const std::vector< double > &  coeffs  ) const

private

Definition at line 298 of file L1GctJetFinderParams.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, patCaloMETCorrections_cff::C, and ctppsDiamondLocalTracks_cfi::threshold.

Referenced by correctionFunction(), and getTauCorrCoeffs().

{

  // The correction fuction is a set of 3rd order polynomials
  //    Et_out = Et_in + (p0 + p1*Et_in + p2*Et_in^2 + p3*Et_in^3)
  // with different coefficients for different energy ranges.
  // The parameters are arranged in groups of five.
  // The first in each group is a threshold value of input Et,
  // followed by the four coefficients for the cubic function.
  double etOut = Et;
  std::vector<double>::const_iterator next_coeff=coeffs.begin();
  while (next_coeff != coeffs.end()) {

    // Read the coefficients from the vector
    double threshold = *next_coeff++;
    double A = *next_coeff++; //p0
    double B = *next_coeff++; //p1
    double C = *next_coeff++; //p2
    double D = *next_coeff++; //p3

    // Check we are in the right energy range and make correction
    if (Et>threshold) {
      etOut += (A + etOut*(B + etOut*(C + etOut*D))) ;
      break;
    }

  }
  return etOut;

}

double L1GctJetFinderParams::powerSeriesCorrect ( const double  Et, const std::vector< double > &  coeffs  ) const

private

Definition at line 255 of file L1GctJetFinderParams.cc.

References mps_fire::i, and funct::pow().

Referenced by correctionFunction(), and getTauCorrCoeffs().

{
  double corrEt = Et;
  for (unsigned i=0; i<coeffs.size();i++) {
    corrEt += coeffs.at(i)*pow(Et,(int)i); 
  }
  return corrEt;
}

template<class Archive >

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

private

void L1GctJetFinderParams::setHtSumParams	(	const double	htLsb,
		const double	htJetEtThresh,
		const double	mhtJetEtThresh
	)

Definition at line 146 of file L1GctJetFinderParams.cc.

References htJetEtThreshold_, htLsb_, and mhtJetEtThreshold_.

Referenced by getMHtJetEtThresholdGct(), and setParams().

{
  htLsb_ = htLsb;
  htJetEtThreshold_ = htJetEtThresh;
  mhtJetEtThreshold_ = mhtJetEtThresh;
}

void L1GctJetFinderParams::setJetEtCalibrationParams	(	const unsigned	corrType,
		const std::vector< std::vector< double > > &	jetCorrCoeffs,
		const std::vector< std::vector< double > > &	tauCorrCoeffs
	)

Definition at line 125 of file L1GctJetFinderParams.cc.

References corrType_, jetCorrCoeffs_, and tauCorrCoeffs_.

Referenced by getMHtJetEtThresholdGct(), and setParams().

{
  corrType_ = corrType;
  jetCorrCoeffs_ = jetCorrCoeffs;
  tauCorrCoeffs_ = tauCorrCoeffs;
}

void L1GctJetFinderParams::setJetEtConvertToEnergyOff

(

)

Definition at line 140 of file L1GctJetFinderParams.cc.

References convertToEnergy_, and energyConversionCoeffs_.

Referenced by getMHtJetEtThresholdGct().

{
  convertToEnergy_ = false;
  energyConversionCoeffs_.clear();
}

void L1GctJetFinderParams::setJetEtConvertToEnergyOn

(

const std::vector< double > &

energyConvCoeffs

)

Definition at line 134 of file L1GctJetFinderParams.cc.

References convertToEnergy_, and energyConversionCoeffs_.

Referenced by getMHtJetEtThresholdGct().

{
  convertToEnergy_ = true;
  energyConversionCoeffs_ = energyConvCoeffs;
}

void L1GctJetFinderParams::setParams	(	const double	rgnEtLsb,
		const double	htLsb,
		const double	cJetSeed,
		const double	fJetSeed,
		const double	tJetSeed,
		const double	tauIsoEtThresh,
		const double	htJetEtThresh,
		const double	mhtJetEtThresh,
		const unsigned	etaBoundary,
		const unsigned	corrType,
		const std::vector< std::vector< double > > &	jetCorrCoeffs,
		const std::vector< std::vector< double > > &	tauCorrCoeffs
	)

Definition at line 160 of file L1GctJetFinderParams.cc.

References setHtSumParams(), setJetEtCalibrationParams(), setRegionEtLsb(), setSlidingWindowParams(), and setTauAlgorithmParams().

Referenced by getMHtJetEtThresholdGct().

{
  setRegionEtLsb (rgnEtLsb);
  setSlidingWindowParams(cJetSeed, fJetSeed, tJetSeed, etaBoundary);
  setJetEtCalibrationParams(corrType, jetCorrCoeffs, tauCorrCoeffs);
  setHtSumParams(htLsb, htJetEtThresh, mhtJetEtThresh);
  setTauAlgorithmParams(tauIsoEtThresh);
}

void L1GctJetFinderParams::setRegionEtLsb

(

const double

rgnEtLsb

)

Definition at line 109 of file L1GctJetFinderParams.cc.

References rgnEtLsb_.

Referenced by getMHtJetEtThresholdGct(), and setParams().

{
  rgnEtLsb_ = rgnEtLsb;
}

void L1GctJetFinderParams::setSlidingWindowParams	(	const double	cJetSeed,
		const double	fJetSeed,
		const double	tJetSeed,
		const unsigned	etaBoundary
	)

Definition at line 114 of file L1GctJetFinderParams.cc.

References cenForJetEtaBoundary_, cenJetEtSeed_, forJetEtSeed_, and tauJetEtSeed_.

Referenced by getMHtJetEtThresholdGct(), and setParams().

{
  cenJetEtSeed_ = cJetSeed;
  forJetEtSeed_ = fJetSeed;
  tauJetEtSeed_ = tJetSeed;
  cenForJetEtaBoundary_ = etaBoundary;
}

void L1GctJetFinderParams::setTauAlgorithmParams

(

const double

tauIsoEtThresh

)

Definition at line 155 of file L1GctJetFinderParams.cc.

References tauIsoEtThreshold_.

Referenced by getMHtJetEtThresholdGct(), and setParams().

{
  tauIsoEtThreshold_ = tauIsoEtThresh;
}

double L1GctJetFinderParams::simpleCorrect ( const double  Et, const std::vector< double > &  coeffs  ) const

private

Definition at line 286 of file L1GctJetFinderParams.cc.

References funct::pow().

Referenced by correctionFunction(), and getTauCorrCoeffs().

{
  // function is :
  //    Et_out = A + B/[ (log10(Et_in))^C + D ] + E/Et_in
  // 
  //    fitcor_as_str = "[0]+[1]/(pow(log10(x),[2])+[3]) + [4]/x"
  // 

  return coeffs.at(0) + coeffs.at(1)/(pow(log10(Et),coeffs.at(2))+coeffs.at(3)) + (coeffs.at(4)/Et);

}

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 134 of file L1GctJetFinderParams.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 134 of file L1GctJetFinderParams.h.

Member Data Documentation

unsigned L1GctJetFinderParams::cenForJetEtaBoundary_

private

Definition at line 122 of file L1GctJetFinderParams.h.

Referenced by correctedEtGeV(), getCenForJetEtaBoundary(), and setSlidingWindowParams().

double L1GctJetFinderParams::cenJetEtSeed_

private

Definition at line 116 of file L1GctJetFinderParams.h.

Referenced by getCenJetEtSeedGct(), getCenJetEtSeedGeV(), and setSlidingWindowParams().

bool L1GctJetFinderParams::convertToEnergy_

private

Definition at line 130 of file L1GctJetFinderParams.h.

Referenced by correctedEtGeV(), getConvertToEnergy(), setJetEtConvertToEnergyOff(), and setJetEtConvertToEnergyOn().

unsigned L1GctJetFinderParams::corrType_

private

Definition at line 125 of file L1GctJetFinderParams.h.

Referenced by correctionFunction(), getCorrType(), L1GctJetFinderParams(), and setJetEtCalibrationParams().

std::vector<double> L1GctJetFinderParams::energyConversionCoeffs_

private

Definition at line 131 of file L1GctJetFinderParams.h.

Referenced by correctedEtGeV(), L1GctJetFinderParams(), setJetEtConvertToEnergyOff(), and setJetEtConvertToEnergyOn().

double L1GctJetFinderParams::forJetEtSeed_

private

Definition at line 117 of file L1GctJetFinderParams.h.

Referenced by getForJetEtSeedGct(), getForJetEtSeedGeV(), and setSlidingWindowParams().

double L1GctJetFinderParams::htJetEtThreshold_

private

Definition at line 120 of file L1GctJetFinderParams.h.

Referenced by getHtJetEtThresholdGct(), getHtJetEtThresholdGeV(), and setHtSumParams().

double L1GctJetFinderParams::htLsb_

private

Definition at line 113 of file L1GctJetFinderParams.h.

Referenced by correctedEtGct(), getHtJetEtThresholdGct(), getHtLsbGeV(), getMHtJetEtThresholdGct(), and setHtSumParams().

std::vector< std::vector<double> > L1GctJetFinderParams::jetCorrCoeffs_

private

Definition at line 126 of file L1GctJetFinderParams.h.

Referenced by correctedEtGeV(), getJetCorrCoeffs(), L1GctJetFinderParams(), and setJetEtCalibrationParams().

double L1GctJetFinderParams::mhtJetEtThreshold_

private

Definition at line 121 of file L1GctJetFinderParams.h.

Referenced by getMHtJetEtThresholdGct(), getMHtJetEtThresholdGeV(), and setHtSumParams().

const unsigned L1GctJetFinderParams::N_CENTRAL_ETA_VALUES = 7

static

Number of eta bins used in correction.

Definition at line 15 of file L1GctJetFinderParams.h.

Referenced by L1GctConfigProducers::L1GctConfigProducers(), L1GctJetFinderParams(), and L1GctJetFinderParamsOnlineProd::newObject().

const unsigned L1GctJetFinderParams::NUMBER_ETA_VALUES = 11

static

Number of eta bins used in correction.

Definition at line 14 of file L1GctJetFinderParams.h.

Referenced by correctedEtGeV(), L1GctConfigProducers::L1GctConfigProducers(), L1GctEmulator::L1GctEmulator(), L1GctJetFinderParams(), L1GctJetFinderParamsOnlineProd::newObject(), L1GctJetFinderBase::setEnergySumMasks(), and L1GctJetFinderBase::setJetEtCalibrationLuts().

double L1GctJetFinderParams::rgnEtLsb_

private

Definition at line 112 of file L1GctJetFinderParams.h.

Referenced by getCenJetEtSeedGct(), getForJetEtSeedGct(), getRgnEtLsbGeV(), getTauIsoEtThresholdGct(), getTauJetEtSeedGct(), and setRegionEtLsb().

std::vector< std::vector<double> > L1GctJetFinderParams::tauCorrCoeffs_

private

Definition at line 127 of file L1GctJetFinderParams.h.

Referenced by correctedEtGeV(), getTauCorrCoeffs(), L1GctJetFinderParams(), and setJetEtCalibrationParams().

double L1GctJetFinderParams::tauIsoEtThreshold_

private

Definition at line 119 of file L1GctJetFinderParams.h.

Referenced by getTauIsoEtThresholdGct(), getTauIsoEtThresholdGeV(), and setTauAlgorithmParams().

double L1GctJetFinderParams::tauJetEtSeed_

private

Definition at line 118 of file L1GctJetFinderParams.h.

Referenced by getTauJetEtSeedGct(), getTauJetEtSeedGeV(), and setSlidingWindowParams().