#include <SimpleJetCorrector.h>

Public Member Functions
float	correction (const std::vector< float > &fX, const std::vector< float > &fY) const

const JetCorrectorParameters &	parameters () const

void	setInterpolation (bool fInterpolation)

	SimpleJetCorrector (const std::string &fDataFile, const std::string &fOption="")

	SimpleJetCorrector (const JetCorrectorParameters &fParameters)

Private Member Functions
float	correctionBin (unsigned fBin, const std::vector< float > &fY) const

unsigned	findInvertVar ()

float	invert (const double args, const double params) const

SimpleJetCorrector &	operator= (const SimpleJetCorrector &)

void	setFuncParameters ()

	SimpleJetCorrector (const SimpleJetCorrector &)

Private Attributes
bool	mDoInterpolation

reco::FormulaEvaluator	mFunc

unsigned	mInvertVar

JetCorrectorParameters	mParameters

Detailed Description

Definition at line 15 of file SimpleJetCorrector.h.

Constructor & Destructor Documentation

SimpleJetCorrector::SimpleJetCorrector	(	const std::string &	fDataFile,
		const std::string &	fOption = `""`
	)

Definition at line 12 of file SimpleJetCorrector.cc.

References JetCorrectorParameters::definitions(), findInvertVar(), JetCorrectorParameters::Definitions::isResponse(), mDoInterpolation, mInvertVar, and mParameters.

Referenced by parameters().

     : mParameters(fDataFile, fOption), mFunc((mParameters.definitions()).formula()) {
   mDoInterpolation = false;
   if (mParameters.definitions().isResponse())
     mInvertVar = findInvertVar();
 }

SimpleJetCorrector::SimpleJetCorrector ( const JetCorrectorParameters & fParameters )

Definition at line 22 of file SimpleJetCorrector.cc.

References JetCorrectorParameters::definitions(), findInvertVar(), JetCorrectorParameters::Definitions::isResponse(), mDoInterpolation, mInvertVar, and mParameters.

     : mParameters(fParameters), mFunc((mParameters.definitions()).formula()) {
   mDoInterpolation = false;
   if (mParameters.definitions().isResponse())
     mInvertVar = findInvertVar();
 }

SimpleJetCorrector::SimpleJetCorrector ( const SimpleJetCorrector & )

private

Member Function Documentation

float SimpleJetCorrector::correction	(	const std::vector< float > &	fX,
		const std::vector< float > &	fY
	)		const

Definition at line 32 of file SimpleJetCorrector.cc.

References newFWLiteAna::bin, JetCorrectorParameters::binIndex(), JetCorrectorParameters::binIndexN(), correctionBin(), JetCorrectorParameters::definitions(), mps_fire::i, mDoInterpolation, mParameters, JetCorrectorParameters::Definitions::nBinVar(), JetCorrectorParameters::neighbourBin(), JetCorrectorParameters::record(), mps_fire::result, createJobs::tmp, and JetCorrectorParameters::Record::xMiddle().

Referenced by setInterpolation().

                                                                                                  {
   float result = 1.;
   float tmp = 0.0;
   float cor = 0.0;
   int bin = -1;
   bin = (fX.size() <= 3 && !fX.empty()) ? mParameters.binIndexN(fX) : mParameters.binIndex(fX);
   if (bin < 0)
     return result;
   if (!mDoInterpolation)
     result = correctionBin(bin, fY);
   else {
     for (unsigned i = 0; i < mParameters.definitions().nBinVar(); i++) {
       float xMiddle[3];
       float xValue[3];
       int prevBin = mParameters.neighbourBin((unsigned)bin, i, false);
       int nextBin = mParameters.neighbourBin((unsigned)bin, i, true);
       if (prevBin >= 0 && nextBin >= 0) {
         xMiddle[0] = mParameters.record(prevBin).xMiddle(i);
         xMiddle[1] = mParameters.record(bin).xMiddle(i);
         xMiddle[2] = mParameters.record(nextBin).xMiddle(i);
         xValue[0] = correctionBin(prevBin, fY);
         xValue[1] = correctionBin(bin, fY);
         xValue[2] = correctionBin(nextBin, fY);
         cor = quadraticInterpolation(fX[i], xMiddle, xValue);
         tmp += cor;
       } else {
         cor = correctionBin(bin, fY);
         tmp += cor;
       }
     }
     result = tmp / mParameters.definitions().nBinVar();
   }
   return result;
 }

float SimpleJetCorrector::correctionBin	(	unsigned	fBin,
		const std::vector< float > &	fY
	)		const

private

Definition at line 69 of file SimpleJetCorrector.cc.

References JetCorrectorParameters::definitions(), reco::FormulaEvaluator::evaluate(), mps_fire::i, invert(), JetCorrectorParameters::Definitions::isResponse(), mFunc, mParameters, N, JetCorrectorParameters::Record::parameters(), CalibrationSummaryClient_cfi::params, JetCorrectorParameters::record(), JetCorrectorParameters::size(), and x.

Referenced by correction(), and parameters().

                                                                                        {
   if (fBin >= mParameters.size()) {
     std::stringstream sserr;
     sserr << "wrong bin: " << fBin << ": only " << mParameters.size() << " available!";
     handleError("SimpleJetCorrector", sserr.str());
   }
   unsigned N = fY.size();
   if (N > 4) {
     std::stringstream sserr;
     sserr << "two many variables: " << N << " maximum is 4";
     handleError("SimpleJetCorrector", sserr.str());
   }
   const std::vector<float>& par = mParameters.record(fBin).parameters();
 
   double params[par.size() - 2 * N];
   for (unsigned int i = 2 * N; i < par.size(); i++) {
     params[i - 2 * N] = par[i];
   }
   double x[4] = {};
   for (unsigned i = 0; i < N; i++) {
     x[i] = (fY[i] < par[2 * i]) ? par[2 * i] : (fY[i] > par[2 * i + 1]) ? par[2 * i + 1] : fY[i];
   }
   if (mParameters.definitions().isResponse()) {
     return invert(x, params);
   }
   return mFunc.evaluate(reco::formula::ArrayAdaptor(x, N), reco::formula::ArrayAdaptor(params, par.size() - 2 * N));
 }

unsigned SimpleJetCorrector::findInvertVar ( )

private

Definition at line 99 of file SimpleJetCorrector.cc.

References JetCorrectorParameters::definitions(), mps_fire::i, mParameters, JetCorrectorParameters::Definitions::parVar(), and mps_fire::result.

Referenced by parameters(), and SimpleJetCorrector().

                                            {
   unsigned result = 9999;
   std::vector<std::string> vv = mParameters.definitions().parVar();
   for (unsigned i = 0; i < vv.size(); i++)
     if (vv[i] == "JetPt") {
       result = i;
       break;
     }
   if (result >= vv.size())
     handleError("SimpleJetCorrector", "Response inversion is required but JetPt is not specified as parameter");
   return result;
 }

float SimpleJetCorrector::invert	(	const double *	args,
		const double *	params
	)		const

private

Definition at line 114 of file SimpleJetCorrector.cc.

References MillePedeFileConverter_cfg::e, reco::FormulaEvaluator::evaluate(), mFunc, mInvertVar, metDiagnosticParameterSet_cfi::nMax, reco::FormulaEvaluator::numberOfParameters(), common_cff::precision, createJobs::tmp, and x.

Referenced by correctionBin(), and parameters().

                                                                                {
   unsigned nMax = 50;
   float precision = 0.0001;
   float rsp = 1.0;
   float e = 1.0;
   double x[4];
   unsigned nLoop = 0;
 
   // 4 dimensions (x, y, z, t) supported in TFormula
   memcpy(&x, args, sizeof(double) * 4);
 
   while (e > precision && nLoop < nMax) {
     rsp = mFunc.evaluate(reco::formula::ArrayAdaptor(x, 4),
                          reco::formula::ArrayAdaptor(params, mFunc.numberOfParameters()));
     float tmp = x[mInvertVar] * rsp;
     e = fabs(tmp - args[mInvertVar]) / args[mInvertVar];
     x[mInvertVar] = args[mInvertVar] / rsp;
     nLoop++;
   }
   return 1. / rsp;
 }