TauJetCorrector Class Reference

#include <TauJetCorrector.h>

Inheritance diagram for TauJetCorrector:

Classes
class	ParametrizationTauJet

Public Member Functions
double	correction (const LorentzVector &fJet) const override
	get correction using Jet information only More...
double	correction (const reco::Jet &) const override
	apply correction using Jet information only More...
bool	eventRequired () const override
	if correction needs event information More...
void	setParameters (std::string, int)
	TauJetCorrector (const edm::ParameterSet &fParameters)
	~TauJetCorrector () override
Public Member Functions inherited from JetCorrector
virtual double	correction (const reco::Jet &fJet, const edm::Event &fEvent, const edm::EventSetup &fSetup) const
	apply correction using all event information More...
virtual double	correction (const reco::Jet &fJet, const edm::RefToBase< reco::Jet > &fJetRef, const edm::Event &fEvent, const edm::EventSetup &fSetup) const
	apply correction using all event information More...
virtual double	correction (const reco::Jet &fJet, const edm::RefToBase< reco::Jet > &fJetRef, const edm::Event &fEvent, const edm::EventSetup &fSetup, LorentzVector &corrected) const
	Apply vectorial correction using all event information. More...
	JetCorrector ()
virtual bool	refRequired () const =0
	if correction needs the jet reference More...
virtual bool	vectorialCorrection () const
	if vectorial correction is provided More...
virtual	~JetCorrector ()

Private Types
typedef std::map< double, ParametrizationTauJet * >	ParametersMap

Private Attributes
ParametersMap	parametrization
int	type

Additional Inherited Members
Public Types inherited from JetCorrector
typedef reco::Particle::LorentzVector	LorentzVector
Static Public Member Functions inherited from JetCorrector
static const JetCorrector *	getJetCorrector (const std::string &fName, const edm::EventSetup &fSetup)
	retrieve corrector from the event setup. troughs exception if something is missing More...

Detailed Description

jet energy corrections from Taujet calibration

Definition at line 13 of file TauJetCorrector.h.

Member Typedef Documentation

ParametersMap

typedef std::map<double, ParametrizationTauJet*> TauJetCorrector::ParametersMap

private

Definition at line 42 of file TauJetCorrector.h.

Constructor & Destructor Documentation

TauJetCorrector()

TauJetCorrector::TauJetCorrector

(

const edm::ParameterSet &

fParameters

)

Definition at line 95 of file TauJetCorrector.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                               {
  type = fConfig.getParameter<int>("TauTriggerType");
  setParameters(fConfig.getParameter<std::string>("tagName"), type);
}

~TauJetCorrector()

TauJetCorrector::~TauJetCorrector ( )

override

Definition at line 100 of file TauJetCorrector.cc.

References stringResolutionProvider_cfi::parametrization.

                                  {
  for (ParametersMap::iterator ip = parametrization.begin(); ip != parametrization.end(); ip++)
    delete ip->second;
}

Member Function Documentation

correction() [1/2]

double TauJetCorrector::correction ( const LorentzVector &  fJet ) const

overridevirtual

get correction using Jet information only

Implements JetCorrector.

Definition at line 182 of file TauJetCorrector.cc.

References EgHLTOffHistBins_cfi::et, PVValHelper::eta, stringResolutionProvider_cfi::parametrization, and edm::second().

                                                                  {
  //cout<<" Start Apply Corrections "<<endl;
  if (parametrization.empty()) {
    return 1.;
  }

  double et = fJet.Et();
  double eta = fabs(fJet.Eta());

  //cout<<" Et and eta of jet "<<et<<" "<<eta<<endl;

  double etnew;
  std::map<double, ParametrizationTauJet*>::const_iterator ip = parametrization.upper_bound(eta);
  etnew = (--ip)->second->value(et, eta);

  //cout<<" The new energy found "<<etnew<<" "<<et<<endl;

  float mScale = etnew / et;

  return mScale;
}

correction() [2/2]

double TauJetCorrector::correction ( const reco::Jet &  fJet ) const

overridevirtual

apply correction using Jet information only

Implements JetCorrector.

Definition at line 204 of file TauJetCorrector.cc.

References pfMETCorrectionType0_cfi::correction, and reco::LeafCandidate::p4().

{ return correction(fJet.p4()); }

eventRequired()

bool TauJetCorrector::eventRequired ( ) const

inlineoverridevirtual

if correction needs event information

Implements JetCorrector.

Definition at line 22 of file TauJetCorrector.h.

{ return false; }

setParameters()

void TauJetCorrector::setParameters	(	std::string	aCalibrationType,
		int	itype
	)

Definition at line 105 of file TauJetCorrector.cc.

References LEDCalibrationChannels::ieta, createfilelist::int, stringResolutionProvider_cfi::parametrization, muonDTDigis_cfi::pset, and edm::second().

                                                                         {
  //cout<< " Start to set parameters "<<endl;
  type = itype;
  JetCalibrationParameterSetTauJet pset(aCalibrationType);

  if ((!pset.valid()) && (aCalibrationType != "no")) {
    edm::LogError("TauJetCorrector:Jet Corrections not found ")
        << aCalibrationType
        << " not found! Cannot apply any correction ... For JetPlusTrack calibration only radii 0.5 and 0.7 are "
           "included for JetParton"
        << endl;
    return;
  }
  if (aCalibrationType == "no")
    return;

  map<int, vector<double> > pq;
  map<int, vector<double> > pg;
  map<int, vector<double> > pqcd;
  map<int, double> etaboundx;
  int iq = 0;
  int ig = 0;
  int iqcd = 0;
  int mtype = 0;
  for (int ieta = 0; ieta < pset.neta(); ieta++) {
    if (pset.type(ieta) == 1) {
      pq[iq] = pset.parameters(ieta);
      iq++;
      mtype = (int)(pset.type(ieta));
    }
    if (pset.type(ieta) == 2) {
      pg[ig] = pset.parameters(ieta);
      ig++;
      mtype = (int)(pset.type(ieta));
    }
    if (pset.type(ieta) == 3) {
      pqcd[iqcd] = pset.parameters(ieta);
      iqcd++;
      mtype = (int)(pset.type(ieta));
    }
    if (pset.type(ieta) == -1) {
      etaboundx[mtype - 1] = pset.eta(ieta);
    }
  }

  //cout<<" Number of parameters "<<iq<<" "<<ig<<" "<<iqcd<<endl;
  int mynum = 0;
  for (int ieta = 0; ieta < pset.neta(); ieta++) {
    //cout<<" New parmetrization "<<ieta<<" "<<pset.type(ieta)<<endl;

    if (pset.type(ieta) == -1)
      continue;
    if (ieta < iq + 1) {
      parametrization[pset.eta(ieta)] =
          new ParametrizationTauJet(pset.type(ieta), (*pq.find(ieta)).second, (*etaboundx.find(0)).second);
      //cout<<" ALL "<<ieta<<" "<<((*pq.find(ieta)).second)[0]<<" "<<((*pq.find(ieta)).second)[1]<<" "<<
      //((*pq.find(ieta)).second)[2]<<endl;
    }
    if (ieta > iq && ieta < iq + ig + 2) {
      mynum = ieta - iq - 1;
      parametrization[pset.eta(ieta)] =
          new ParametrizationTauJet(pset.type(ieta), (*pg.find(mynum)).second, (*etaboundx.find(1)).second);
      //cout<<" One prong "<<((*pg.find(mynum)).second)[0]<<" "<<((*pg.find(mynum)).second)[1]<<" "<<
      //((*pg.find(mynum)).second)[2]<<endl;
    }
    if (ieta > iq + ig + 1) {
      mynum = ieta - iq - ig - 2;
      //cout<<" Mynum "<<mynum<<" "<<ieta<<" "<<pset.type(ieta)<<endl;
      parametrization[pset.eta(ieta)] =
          new ParametrizationTauJet(pset.type(ieta), (*pqcd.find(mynum)).second, (*etaboundx.find(2)).second);
      //cout<<" Two prongs "<<((*pqcd.find(mynum)).second)[0]<<" "<<((*pqcd.find(mynum)).second)[1]<<" "<<
      //((*pqcd.find(mynum)).second)[2]<<endl;
    }
  }
  //cout<<" Parameters inserted into mAlgorithm "<<endl;
}

Member Data Documentation

parametrization

ParametersMap TauJetCorrector::parametrization

private

Definition at line 43 of file TauJetCorrector.h.

type

int TauJetCorrector::type

private

Definition at line 44 of file TauJetCorrector.h.

Referenced by core.autovars.NTupleVariable::makeBranch(), and core.AutoHandle.AutoHandle::ReallyLoad().