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 16 of file SimpleZSPJPTJetCorrector.h.

Constructor & Destructor Documentation

SimpleZSPJPTJetCorrector::SimpleZSPJPTJetCorrector

(

)

Definition at line 21 of file SimpleZSPJPTJetCorrector.cc.

{ debug_=false;}

SimpleZSPJPTJetCorrector::SimpleZSPJPTJetCorrector

(

const std::string &

fDataFile

)

Definition at line 24 of file SimpleZSPJPTJetCorrector.cc.

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

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

SimpleZSPJPTJetCorrector::~SimpleZSPJPTJetCorrector ( )

virtual

Definition at line 45 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 53 of file SimpleZSPJPTJetCorrector.cc.

References gather_cfg::cout, eta, V0MonitoringClient_cfi::factor, 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 48 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 99 of file SimpleZSPJPTJetCorrector.cc.

References gather_cfg::cout, eta, V0MonitoringClient_cfi::factor, 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 32 of file SimpleZSPJPTJetCorrector.h.

TFormula* SimpleZSPJPTJetCorrector::mFunc

private

Definition at line 31 of file SimpleZSPJPTJetCorrector.h.

JetCorrectorParameters* SimpleZSPJPTJetCorrector::mParameters

private

Definition at line 30 of file SimpleZSPJPTJetCorrector.h.