#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 24 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 36 of file SimpleJetCorrector.cc.

References stringResolutionProvider_cfi::bin, JetCorrectorParameters::binIndex(), JetCorrectorParameters::binIndexN(), correctionBin(), JetCorrectorParameters::definitions(), mps_fire::i, mDoInterpolation, mParameters, JetCorrectorParameters::Definitions::nBinVar(), JetCorrectorParameters::neighbourBin(), JetCorrectorParameters::record(), mps_fire::result, 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.size()>0) ? 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 79 of file SimpleJetCorrector.cc.

References JetCorrectorParameters::definitions(), reco::FormulaEvaluator::evaluate(), mps_fire::i, invert(), JetCorrectorParameters::Definitions::isResponse(), mFunc, mParameters, N, JetCorrectorParameters::Record::parameters(), 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 114 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 131 of file SimpleJetCorrector.cc.

References MillePedeFileConverter_cfg::e, reco::FormulaEvaluator::evaluate(), mFunc, mInvertVar, reco::FormulaEvaluator::numberOfParameters(), 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;
 }