#include <EtaDependentPileup.h>

Inheritance diagram for fftjetcms::EtaDependentPileup:

Public Member Functions
	EtaDependentPileup (const fftjet::LinearInterpolator2d &i, double inputRhoFactor, double outputRhoFactor)

bool	isPhiDependent () const override

double	operator() (double eta, double phi, const reco::FFTJetPileupSummary &summary) const override

	~EtaDependentPileup () override

Public Member Functions inherited from fftjetcms::AbsPileupCalculator
virtual	~AbsPileupCalculator ()

Private Attributes
double	etaMax_

double	etaMin_

double	inputRhoFactor_

fftjet::LinearInterpolator2d	interp_

double	outputRhoFactor_

double	rhoMax_

double	rhoMin_

double	rhoStep_

Detailed Description

Definition at line 17 of file EtaDependentPileup.h.

Constructor & Destructor Documentation

fftjetcms::EtaDependentPileup::EtaDependentPileup	(	const fftjet::LinearInterpolator2d &	i,
		double	inputRhoFactor,
		double	outputRhoFactor
	)

Definition at line 6 of file EtaDependentPileup.cc.

References cms::cuda::assert(), etaMax_, etaMin_, rhoMax_, rhoMin_, and rhoStep_.

       : interp_(interp),
         inputRhoFactor_(inputRhoFactor),
         outputRhoFactor_(outputRhoFactor),
         etaMin_(interp.xMin()),
         etaMax_(interp.xMax()),
         rhoMin_(interp.yMin()),
         rhoMax_(interp.yMax()),
         rhoStep_((rhoMax_ - rhoMin_) / interp.ny()) {
     const double etaStep = (etaMax_ - etaMin_) / interp.nx();
     etaMin_ += etaStep * 0.5;
     etaMax_ -= etaStep * 0.5;
     rhoMin_ += rhoStep_ * 0.5;
     rhoMax_ -= rhoStep_ * 0.5;
     assert(etaMin_ < etaMax_);
     assert(rhoMin_ < rhoMax_);
   }

fftjetcms::EtaDependentPileup::~EtaDependentPileup ( )

inlineoverride

Definition at line 21 of file EtaDependentPileup.h.

21 {}

Member Function Documentation

bool fftjetcms::EtaDependentPileup::isPhiDependent ( ) const

inlineoverridevirtual

Implements fftjetcms::AbsPileupCalculator.

Definition at line 25 of file EtaDependentPileup.h.

25 { return false; }

double fftjetcms::EtaDependentPileup::operator()	(	double	eta,
		double	phi,
		const reco::FFTJetPileupSummary &	summary
	)		const

overridevirtual

Implements fftjetcms::AbsPileupCalculator.

Definition at line 26 of file EtaDependentPileup.cc.

References etaMax_, etaMin_, inputRhoFactor_, interp_, outputRhoFactor_, reco::FFTJetPileupSummary::pileupRho(), rhoMax_, rhoMin_, rhoStep_, and relativeConstraints::value.

                                                                                                           {
     double value = 0.0;
     const double rho = summary.pileupRho() * inputRhoFactor_;
     if (eta < etaMin_)
       eta = etaMin_;
     if (eta > etaMax_)
       eta = etaMax_;
     if (rho >= rhoMin_ && rho <= rhoMax_)
       value = interp_(eta, rho);
     else {
       double x0, x1;
       if (rho < rhoMin_) {
         x0 = rhoMin_;
         x1 = rhoMin_ + rhoStep_ * 0.5;
       } else {
         x0 = rhoMax_;
         x1 = rhoMax_ - rhoStep_ * 0.5;
       }
       const double z0 = interp_(eta, x0);
       const double z1 = interp_(eta, x1);
       value = z0 + (z1 - z0) * ((rho - x0) / (x1 - x0));
     }
     return (value > 0.0 ? value : 0.0) * outputRhoFactor_;
   }