SimpleZSPJPTJetCorrector Class Reference

classes declaration More...

#include <SimpleZSPJPTJetCorrector.h>

Public Member Functions
virtual double	correctionEtEtaPhiP (double fEt, double fEta, double fPhi, double fP) const
virtual double	correctionPtEtaPhiE (double fPt, double fEta, double fPhi, double fE) const
virtual double	correctionPUEtEtaPhiP (double fEt, double fEta, double fPhi, double fP) const
void	init (const std::string &fDataFile)
	SimpleZSPJPTJetCorrector ()
	SimpleZSPJPTJetCorrector (const std::string &fDataFile)
virtual	~SimpleZSPJPTJetCorrector ()

Private Member Functions
SimpleZSPJPTJetCorrector &	operator= (const SimpleZSPJPTJetCorrector &)
	SimpleZSPJPTJetCorrector (const SimpleZSPJPTJetCorrector &)

Private Attributes
bool	debug_
TFormula *	mFunc
JetCorrectorParameters *	mParameters

Detailed Description

classes declaration

Definition at line 15 of file SimpleZSPJPTJetCorrector.h.

Constructor & Destructor Documentation

SimpleZSPJPTJetCorrector() [1/3]

SimpleZSPJPTJetCorrector::SimpleZSPJPTJetCorrector

(

)

Definition at line 20 of file SimpleZSPJPTJetCorrector.cc.

{ debug_ = false; }

SimpleZSPJPTJetCorrector() [2/3]

SimpleZSPJPTJetCorrector::SimpleZSPJPTJetCorrector

(

const std::string &

fDataFile

)

Definition at line 22 of file SimpleZSPJPTJetCorrector.cc.

                                                                             {
  debug_ = false;
  mParameters = new JetCorrectorParameters(fDataFile, "");
 
  //std::cout<<" Formula "<<((mParameters->definitions()).formula()).c_str()<<std::endl;
 
  mFunc = new TFormula("function", ((mParameters->definitions()).formula()).c_str());
 
  // Read parameters
  if (debug_) {
    std::cout << " Size of parameters as read by SimpleJetCorrectorParameters " << mParameters->size() << std::endl;
    for (unsigned int i = 0; i < mParameters->size(); i++) {
      const std::vector<float> p = mParameters->record(i).parameters();
      for (std::vector<float>::const_iterator j = p.begin(); j < p.end(); j++) {
        std::cout << " Parameter number " << mParameters->record(i).xMin(0) << " " << mParameters->record(i).xMax(0)
                  << " " << (*j) << std::endl;
      }
    }
  }  // debug
}

References gather_cfg::cout, pfMETCorrectionType0_cfi::formula, mps_fire::i, dqmiolumiharvest::j, and AlCaHLTBitMon_ParallelJobs::p.

~SimpleZSPJPTJetCorrector()

SimpleZSPJPTJetCorrector::~SimpleZSPJPTJetCorrector ( )

virtual

Definition at line 43 of file SimpleZSPJPTJetCorrector.cc.

{}

SimpleZSPJPTJetCorrector() [3/3]

SimpleZSPJPTJetCorrector::SimpleZSPJPTJetCorrector ( const SimpleZSPJPTJetCorrector &  )

private

Member Function Documentation

correctionEtEtaPhiP()

double SimpleZSPJPTJetCorrector::correctionEtEtaPhiP ( double  fEt, double  fEta, double  fPhi, double  fP  ) const

virtual

Definition at line 50 of file SimpleZSPJPTJetCorrector.cc.

                                                                                                          {
  double et = fEt;
  double eta = fabs(fEta);
  double factor = 1.;
 
  // Define Eta bin for parametrization
  std::vector<float> xx;
  xx.push_back(eta);
  int band = mParameters->binIndex(xx);
 
  if (band < 0)
    return factor;
 
  const std::vector<float> p = mParameters->record(band).parameters();
 
  // Set parameters
  for (unsigned int i = 0; i < p.size(); i++) {
    mFunc->SetParameter(i, p[i]);
  }
 
  if (debug_) {
    cout << " Et and eta of jet " << et << " " << eta << " bin " << band << " " << p[1] << " " << p[2] << " " << p[3]
         << " " << p[4] << " " << p[5] << endl;
  }
 
  //  double koef = 1. - p[2]*exp(p[3]*et)-p[4]*exp(p[5]*et);
 
  double koef = 1. - mFunc->Eval(et);
 
  // If et calojet less then some value - use correction on the boundary
 
  if (et < p[0])
    koef = 1. - mFunc->Eval(p[0]);
 
  //
  if (koef <= 0.000001) {
    if (debug_)
      std::cout << "SimpleZSPJPTJetCorrector::Problem with ZSP corrections " << koef << std::endl;
    koef = 1.;
  }
 
  double etnew = et / koef;
 
  if (debug_)
    cout << " The new energy found " << etnew << " " << et << " " << koef << endl;
 
  return etnew / et;
}

References gather_cfg::cout, EgHLTOffHistBins_cfi::et, PVValHelper::eta, DQMScaleToClient_cfi::factor, mps_fire::i, AlCaHLTBitMon_ParallelJobs::p, and geometryCSVtoXML::xx.

correctionPtEtaPhiE()

double SimpleZSPJPTJetCorrector::correctionPtEtaPhiE ( double  fPt, double  fEta, double  fPhi, double  fE  ) const

virtual

Definition at line 45 of file SimpleZSPJPTJetCorrector.cc.

                                                                                                          {
  double costhetainv = cosh(fEta);
  return correctionEtEtaPhiP(fE / costhetainv, fEta, fPhi, fPt * costhetainv);
}

correctionPUEtEtaPhiP()

double SimpleZSPJPTJetCorrector::correctionPUEtEtaPhiP ( double  fEt, double  fEta, double  fPhi, double  fP  ) const

virtual

Definition at line 99 of file SimpleZSPJPTJetCorrector.cc.

                                                                                                            {
  double et = fEt;
  double eta = fabs(fEta);
  double factor = 1.;
 
  // Define Eta bin for parametrization
  std::vector<float> xx;
  xx.push_back(eta);
  int band = mParameters->binIndex(xx);
 
  if (band < 0)
    return factor;
 
  const std::vector<float> p = mParameters->record(band).parameters();
 
  if (debug_) {
    cout << " Et and eta of jet " << et << " " << eta << " bin " << band << std::endl;
  }
 
  double koef = (et - p[2]) / et;
  double etnew = et / koef;
 
  if (debug_)
    cout << " The new energy found " << etnew << " " << et << endl;
 
  return etnew / et;
}

References gather_cfg::cout, EgHLTOffHistBins_cfi::et, PVValHelper::eta, DQMScaleToClient_cfi::factor, AlCaHLTBitMon_ParallelJobs::p, and geometryCSVtoXML::xx.

init()

void SimpleZSPJPTJetCorrector::init

(

const std::string &

fDataFile

)

operator=()

SimpleZSPJPTJetCorrector& SimpleZSPJPTJetCorrector::operator= ( const SimpleZSPJPTJetCorrector &  )

private

Member Data Documentation

debug_

bool SimpleZSPJPTJetCorrector::debug_

private

Definition at line 31 of file SimpleZSPJPTJetCorrector.h.

mFunc

TFormula* SimpleZSPJPTJetCorrector::mFunc

private

Definition at line 30 of file SimpleZSPJPTJetCorrector.h.

mParameters

JetCorrectorParameters* SimpleZSPJPTJetCorrector::mParameters

private

Definition at line 29 of file SimpleZSPJPTJetCorrector.h.