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::SimpleZSPJPTJetCorrector

(

)

Definition at line 20 of file SimpleZSPJPTJetCorrector.cc.

References debug_.

{ debug_ = false; }

SimpleZSPJPTJetCorrector::SimpleZSPJPTJetCorrector

(

const std::string &

fDataFile

)

Definition at line 22 of file SimpleZSPJPTJetCorrector.cc.

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

                                                                             {
  debug_ = false;
  mParameters = new JetCorrectorParameters(fDataFile, "");

  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
}

SimpleZSPJPTJetCorrector::~SimpleZSPJPTJetCorrector ( )

virtual

Definition at line 41 of file SimpleZSPJPTJetCorrector.cc.

{}

SimpleZSPJPTJetCorrector::SimpleZSPJPTJetCorrector ( const SimpleZSPJPTJetCorrector &  )

private

Member Function Documentation

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

virtual

Definition at line 48 of file SimpleZSPJPTJetCorrector.cc.

References gather_cfg::cout, debug_, PVValHelper::eta, mps_fire::i, and AlCaHLTBitMon_ParallelJobs::p.

                                                                                                          {
  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;
}

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

virtual

Definition at line 43 of file SimpleZSPJPTJetCorrector.cc.

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

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

virtual

Definition at line 97 of file SimpleZSPJPTJetCorrector.cc.

References gather_cfg::cout, debug_, PVValHelper::eta, and AlCaHLTBitMon_ParallelJobs::p.

                                                                                                            {
  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;
}

void SimpleZSPJPTJetCorrector::init

(

const std::string &

fDataFile

)

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

private

Member Data Documentation

bool SimpleZSPJPTJetCorrector::debug_

private

Definition at line 31 of file SimpleZSPJPTJetCorrector.h.

TFormula* SimpleZSPJPTJetCorrector::mFunc

private

Definition at line 30 of file SimpleZSPJPTJetCorrector.h.

JetCorrectorParameters* SimpleZSPJPTJetCorrector::mParameters

private

Definition at line 29 of file SimpleZSPJPTJetCorrector.h.