#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.
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.
unsigned	getCenForJetEtaBoundary () const
unsigned	getCenJetEtSeedGct () const
double	getCenJetEtSeedGeV () const
bool	getConvertToEnergy () const
unsigned	getCorrType () const
	Access to jet Et calibration parameters.
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, std::vector< std::vector< double > > jetCorrCoeffs, std::vector< std::vector< double > > tauCorrCoeffs, bool convertToEnergy, 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.
static const unsigned	NUMBER_ETA_VALUES = 11
	Number of eta bins used in correction.
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
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_

Detailed Description

Definition at line 8 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,
		std::vector< std::vector< double > >	jetCorrCoeffs,
		std::vector< std::vector< double > >	tauCorrCoeffs,
		bool	convertToEnergy,
		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.

{
  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, query::result, and tauCorrCoeffs_.

Referenced by 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(), query::result, and simpleCorrect().

Referenced by correctedEtGeV().

{
  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]`

unsigned L1GctJetFinderParams::getCenForJetEtaBoundary ( ) const [inline]

Definition at line 43 of file L1GctJetFinderParams.h.

References cenForJetEtaBoundary_.

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

{ return cenForJetEtaBoundary_; }

unsigned L1GctJetFinderParams::getCenJetEtSeedGct ( ) const [inline]

Definition at line 47 of file L1GctJetFinderParams.h.

References cenJetEtSeed_, and rgnEtLsb_.

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

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

double L1GctJetFinderParams::getCenJetEtSeedGeV ( ) const [inline]

Definition at line 37 of file L1GctJetFinderParams.h.

References cenJetEtSeed_.

Referenced by operator<<().

{ return cenJetEtSeed_; }

bool L1GctJetFinderParams::getConvertToEnergy ( ) const [inline]

Definition at line 44 of file L1GctJetFinderParams.h.

References convertToEnergy_.

{ return convertToEnergy_; }

unsigned L1GctJetFinderParams::getCorrType ( ) const [inline]

Access to jet Et calibration parameters.

Definition at line 90 of file L1GctJetFinderParams.h.

References corrType_.

Referenced by operator<<().

{ return corrType_; }

unsigned L1GctJetFinderParams::getForJetEtSeedGct ( ) const [inline]

Definition at line 48 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 38 of file L1GctJetFinderParams.h.

References forJetEtSeed_.

Referenced by operator<<().

{ return forJetEtSeed_; }

unsigned L1GctJetFinderParams::getHtJetEtThresholdGct ( ) const [inline]

Definition at line 51 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 41 of file L1GctJetFinderParams.h.

References htJetEtThreshold_.

Referenced by operator<<().

{ return htJetEtThreshold_; }

double L1GctJetFinderParams::getHtLsbGeV ( ) const [inline]

Definition at line 36 of file L1GctJetFinderParams.h.

References htLsb_.

Referenced by operator<<(), and L1GctGlobalEnergyAlgos::setJetFinderParams().

{ return htLsb_; }

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

Definition at line 91 of file L1GctJetFinderParams.h.

References jetCorrCoeffs_.

Referenced by operator<<().

{ return jetCorrCoeffs_; }

unsigned L1GctJetFinderParams::getMHtJetEtThresholdGct ( ) const [inline]

Definition at line 52 of file L1GctJetFinderParams.h.

References htLsb_, and mhtJetEtThreshold_.

Referenced by L1GctJetFinderBase::setJetFinderParams().

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

double L1GctJetFinderParams::getMHtJetEtThresholdGeV ( ) const [inline]

Definition at line 42 of file L1GctJetFinderParams.h.

References mhtJetEtThreshold_.

Referenced by operator<<().

{ return mhtJetEtThreshold_; }

double L1GctJetFinderParams::getRgnEtLsbGeV ( ) const [inline]

Definition at line 35 of file L1GctJetFinderParams.h.

References rgnEtLsb_.

Referenced by operator<<().

{ return rgnEtLsb_; }

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

Definition at line 92 of file L1GctJetFinderParams.h.

References tauCorrCoeffs_.

Referenced by operator<<().

{ return tauCorrCoeffs_; }

unsigned L1GctJetFinderParams::getTauIsoEtThresholdGct ( ) const [inline]

Definition at line 50 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 40 of file L1GctJetFinderParams.h.

References tauIsoEtThreshold_.

Referenced by operator<<().

{ return tauIsoEtThreshold_; }

unsigned L1GctJetFinderParams::getTauJetEtSeedGct ( ) const [inline]

Definition at line 49 of file L1GctJetFinderParams.h.

References rgnEtLsb_, and tauJetEtSeed_.

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

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

double L1GctJetFinderParams::getTauJetEtSeedGeV ( ) const [inline]

Definition at line 39 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 funct::A, funct::C, mathSSE::sqrt(), and dtDQMClient_cfg::threshold.

Referenced by correctionFunction().

{

  // 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 funct::exp(), and funct::pow().

Referenced by correctionFunction().

{  
  //
  // 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 funct::A, funct::C, funct::D, and dtDQMClient_cfg::threshold.

Referenced by correctionFunction().

{

  // 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 i, and funct::pow().

Referenced by correctionFunction().

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

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

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

void L1GctJetFinderParams::setJetEtConvertToEnergyOff ( )

Definition at line 140 of file L1GctJetFinderParams.cc.

References convertToEnergy_, and energyConversionCoeffs_.

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

void L1GctJetFinderParams::setJetEtConvertToEnergyOn ( const std::vector< double > & energyConvCoeffs )

Definition at line 134 of file L1GctJetFinderParams.cc.

References convertToEnergy_, and energyConversionCoeffs_.

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

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

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

}