JetPartonCorrector Class Reference

#include <JetPartonCorrector.h>

Inheritance diagram for JetPartonCorrector:

Public Member Functions
double	correction (const LorentzVector &fJet) const override
	get correction using Jet information only More...
bool	eventRequired () const override
	if correction needs event information More...
	JetPartonCorrector (const edm::ParameterSet &fConfig)
void	setParameters (std::string aCalibrationType, double aJetFinderRadius, int aPartonMixture)
	~JetPartonCorrector () override
Public Member Functions inherited from JetCorrector
virtual double	correction (const reco::Jet &fJet) const =0
	apply correction using Jet information only More...
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, JetPartonNamespace::ParametrizationJetParton * >	ParametersMap

Private Attributes
ParametersMap	parametrization
double	theJetFinderRadius
int	thePartonMixture

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

Definition at line 21 of file JetPartonCorrector.h.

Member Typedef Documentation

typedef std::map<double, JetPartonNamespace::ParametrizationJetParton*> JetPartonCorrector::ParametersMap

private

Definition at line 34 of file JetPartonCorrector.h.

Constructor & Destructor Documentation

JetPartonCorrector::JetPartonCorrector

(

const edm::ParameterSet &

fConfig

)

Definition at line 143 of file JetPartonCorrector.cc.

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

                                                                     {
  thePartonMixture = fConfig.getParameter<int>("MixtureType");
  theJetFinderRadius = fConfig.getParameter<double>("Radius");
  setParameters(fConfig.getParameter<std::string>("tagName"), theJetFinderRadius, thePartonMixture);
}

JetPartonCorrector::~JetPartonCorrector ( )

override

Definition at line 149 of file JetPartonCorrector.cc.

References stringResolutionProvider_cfi::parametrization.

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

Member Function Documentation

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

overridevirtual

get correction using Jet information only

Implements JetCorrector.

Definition at line 197 of file JetPartonCorrector.cc.

References gather_cfg::cout, EgHLTOffHistBins_cfi::et, PVValHelper::eta, HLT_2018_cff::eta1, HLT_2018_cff::eta2, stringResolutionProvider_cfi::parametrization, and edm::second().

                                                                     {
  if (parametrization.empty()) {
    return 1.;
  }

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

  //if(eta<10) { eta=abs(fJet.getY()); }

  double etnew;
  std::map<double, JetPartonNamespace::ParametrizationJetParton*>::const_iterator ip = parametrization.upper_bound(eta);
  if (ip == parametrization.begin()) {
    etnew = ip->second->value(et, eta);
  } else if (ip == parametrization.end()) {
    etnew = (--ip)->second->value(et, eta);
  } else {
    double eta2 = ip->first;
    double et2 = ip->second->value(et, eta);
    ip--;
    double eta1 = ip->first;
    double et1 = ip->second->value(et, eta);

    etnew = (eta2 * et1 - eta1 * et2 + eta * et2 - eta * et1) / (eta2 - eta1);
  }
  cout << " JetParton::The new energy found " << etnew << " " << et << endl;
  float mScale = 1000.;

  if (et > 0.001)
    mScale = etnew / et;

  return mScale;
}

bool JetPartonCorrector::eventRequired ( ) const

inlineoverridevirtual

if correction needs event information

Implements JetCorrector.

Definition at line 31 of file JetPartonCorrector.h.

{ return false; }

void JetPartonCorrector::setParameters	(	std::string	aCalibrationType,
		double	aJetFinderRadius,
		int	aPartonMixture
	)

Definition at line 154 of file JetPartonCorrector.cc.

References JetPartonNamespace::JetPartonCalibrationParameterSet::eta(), LEDCalibrationChannels::ieta, JetPartonNamespace::JetPartonCalibrationParameterSet::neta(), JetPartonNamespace::JetPartonCalibrationParameterSet::parameters(), stringResolutionProvider_cfi::parametrization, muonDTDigis_cfi::pset, edm::second(), JetPartonNamespace::JetPartonCalibrationParameterSet::type(), and JetPartonNamespace::JetPartonCalibrationParameterSet::valid().

                                                                                                              {
  theJetFinderRadius = aJetFinderRadius;
  thePartonMixture = aPartonMixture;

  JetPartonNamespace::JetPartonCalibrationParameterSet pset(aCalibrationType);

  if ((!pset.valid()) && (aCalibrationType != "no")) {
    edm::LogError("JetPartonCorrector: 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;
  int iq = 0;
  int ig = 0;
  int iqcd = 0;
  for (int ieta = 0; ieta < pset.neta(); ieta++) {
    if (pset.type(ieta) == 1) {
      pq[iq] = pset.parameters(ieta);
      iq++;
    };
    if (pset.type(ieta) == 2) {
      pg[ig] = pset.parameters(ieta);
      ig++;
    };
    if (pset.type(ieta) == 3) {
      pqcd[iqcd] = pset.parameters(ieta);
      iqcd++;
    };
  }

  for (int ieta = 0; ieta < iq; ieta++) {
    parametrization[pset.eta(ieta)] = new JetPartonNamespace::ParametrizationJetParton(
        thePartonMixture, (*pq.find(ieta)).second, (*pg.find(ieta)).second, (*pqcd.find(ieta)).second);
  }
}

Member Data Documentation

ParametersMap JetPartonCorrector::parametrization

private

Definition at line 35 of file JetPartonCorrector.h.

double JetPartonCorrector::theJetFinderRadius

private

Definition at line 37 of file JetPartonCorrector.h.

int JetPartonCorrector::thePartonMixture

private

Definition at line 36 of file JetPartonCorrector.h.