SimpleJetCorrectionUncertainty Class Reference

#include <SimpleJetCorrectionUncertainty.h>

Public Member Functions
const JetCorrectorParameters &	parameters () const
	SimpleJetCorrectionUncertainty ()
	SimpleJetCorrectionUncertainty (const std::string &fDataFile)
	SimpleJetCorrectionUncertainty (const JetCorrectorParameters &fParameters)
float	uncertainty (const std::vector< float > &fX, float fY, bool fDirection) const
	~SimpleJetCorrectionUncertainty ()

Private Member Functions
int	findBin (const std::vector< float > &v, float x) const
float	linearInterpolation (float fZ, const float fX[2], const float fY[2]) const
SimpleJetCorrectionUncertainty &	operator= (const SimpleJetCorrectionUncertainty &)
	SimpleJetCorrectionUncertainty (const SimpleJetCorrectionUncertainty &)
float	uncertaintyBin (unsigned fBin, float fY, bool fDirection) const

Private Attributes
JetCorrectorParameters *	mParameters

Detailed Description

Definition at line 10 of file SimpleJetCorrectionUncertainty.h.

Constructor & Destructor Documentation

SimpleJetCorrectionUncertainty::SimpleJetCorrectionUncertainty

(

)

Definition at line 8 of file SimpleJetCorrectionUncertainty.cc.

References mParameters.

{
  mParameters = new JetCorrectorParameters();
}

SimpleJetCorrectionUncertainty::SimpleJetCorrectionUncertainty

(

const std::string &

fDataFile

)

Definition at line 13 of file SimpleJetCorrectionUncertainty.cc.

References mParameters.

{
  mParameters = new JetCorrectorParameters(fDataFile);
}

SimpleJetCorrectionUncertainty::SimpleJetCorrectionUncertainty

(

const JetCorrectorParameters &

fParameters

)

Definition at line 18 of file SimpleJetCorrectionUncertainty.cc.

References mParameters.

{
  mParameters = new JetCorrectorParameters(fParameters);
}

SimpleJetCorrectionUncertainty::~SimpleJetCorrectionUncertainty

(

)

Definition at line 23 of file SimpleJetCorrectionUncertainty.cc.

References mParameters.

{
  delete mParameters;
}

SimpleJetCorrectionUncertainty::SimpleJetCorrectionUncertainty ( const SimpleJetCorrectionUncertainty &  )

private

Member Function Documentation

int SimpleJetCorrectionUncertainty::findBin ( const std::vector< float > &  v, float  x  ) const

private

Definition at line 102 of file SimpleJetCorrectionUncertainty.cc.

References i, and n.

Referenced by uncertaintyBin().

{
  int i;
  int n = v.size()-1;
  if (n<=0) return -1;
  if (x<v[0] || x>=v[n])
    return -1;
  for(i=0;i<n;i++)
   {
     if (x>=v[i] && x<v[i+1])
       return i;
   }
  return 0; 
}

float SimpleJetCorrectionUncertainty::linearInterpolation ( float  fZ, const float  fX[2], const float  fY[2]  ) const

private

Definition at line 79 of file SimpleJetCorrectionUncertainty.cc.

References a, b, and alignCSCRings::r.

Referenced by uncertaintyBin().

{
  // Linear interpolation through the points (x[i],y[i]). First find the line that
  // is defined by the points and then calculate the y(z).
  float r = 0;
  if (fX[0] == fX[1])
    {
      if (fY[0] == fY[1])
        r = fY[0];
      else {
    edm::LogError("SimpleJetCorrectionUncertainty")<<" interpolation error";
    return -999.0;
      }
    } 
  else   
    {
      float a = (fY[1]-fY[0])/(fX[1]-fX[0]);
      float b = (fY[0]*fX[1]-fY[1]*fX[0])/(fX[1]-fX[0]);
      r = a*fZ+b;
    }
  return r;
}

SimpleJetCorrectionUncertainty& SimpleJetCorrectionUncertainty::operator= ( const SimpleJetCorrectionUncertainty &  )

private

const JetCorrectorParameters& SimpleJetCorrectionUncertainty::parameters ( void  ) const

inline

Definition at line 17 of file SimpleJetCorrectionUncertainty.h.

References mParameters.

Referenced by JetCorrectionUncertainty::getUncertainty(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

{return *mParameters;}

float SimpleJetCorrectionUncertainty::uncertainty	(	const std::vector< float > &	fX,
		float	fY,
		bool	fDirection
	)		const

Definition at line 28 of file SimpleJetCorrectionUncertainty.cc.

References newFWLiteAna::bin, JetCorrectorParameters::binIndex(), mParameters, query::result, and uncertaintyBin().

Referenced by JetCorrectionUncertainty::getUncertainty().

{
  float result = 1.;
  int bin = mParameters->binIndex(fX);
  if (bin<0) {
    edm::LogError("SimpleJetCorrectionUncertainty")<<" bin variables out of range";
    result = -999.0;
  } else 
    result = uncertaintyBin((unsigned)bin,fY,fDirection);
  return result;
}

float SimpleJetCorrectionUncertainty::uncertaintyBin ( unsigned  fBin, float  fY, bool  fDirection  ) const

private

Definition at line 40 of file SimpleJetCorrectionUncertainty.cc.

References newFWLiteAna::bin, edm::hlt::Exception, findBin(), i, linearInterpolation(), mParameters, N, AlCaHLTBitMon_ParallelJobs::p, JetCorrectorParameters::Record::parameters(), JetCorrectorParameters::record(), query::result, JetCorrectorParameters::size(), and relativeConstraints::value.

Referenced by uncertainty().

{
  if (fBin >= mParameters->size()) { 
    edm::LogError("SimpleJetCorrectionUncertainty")<<" wrong bin: "<<fBin<<": only "<<mParameters->size()<<" are available";
    return -999.0;
  }
  const std::vector<float>& p = mParameters->record(fBin).parameters();
  if ((p.size() % 3) != 0)
    throw cms::Exception ("SimpleJetCorrectionUncertainty")<<"wrong # of parameters: multiple of 3 expected, "<<p.size()<< " got";
  std::vector<float> yGrid,value;
  unsigned int N = p.size()/3;
  float result = -1.0;
  for(unsigned i=0;i<N;i++)
    {
      unsigned ind = 3*i;
      yGrid.push_back(p[ind]);
      if (fDirection)// true = UP
        value.push_back(p[ind+1]);
      else // false = DOWN
        value.push_back(p[ind+2]); 
    }
  if (fY <= yGrid[0])
    result = value[0];  
  else if (fY >= yGrid[N-1])
    result = value[N-1]; 
  else
    {
      int bin = findBin(yGrid,fY); 
      float vx[2],vy[2];
      for(int i=0;i<2;i++)
        {
          vx[i] = yGrid[bin+i]; 
          vy[i] = value[bin+i];
        } 
      result = linearInterpolation(fY,vx,vy);
    }
  return result;
}

Member Data Documentation

JetCorrectorParameters* SimpleJetCorrectionUncertainty::mParameters

private

Definition at line 26 of file SimpleJetCorrectionUncertainty.h.

Referenced by parameters(), SimpleJetCorrectionUncertainty(), uncertainty(), uncertaintyBin(), and ~SimpleJetCorrectionUncertainty().