FFTJetPileupEstimator Class Reference

#include <RecoJets/FFTJetProducers/plugins/FFTJetPileupEstimator.cc>

Inheritance diagram for FFTJetPileupEstimator:

Public Member Functions
	FFTJetPileupEstimator (const edm::ParameterSet &)
	~FFTJetPileupEstimator ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Protected Member Functions
void	beginJob ()
void	endJob ()
void	produce (edm::Event &, const edm::EventSetup &)
Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Private Member Functions
template<class Ptr >
void	checkConfig (const Ptr &ptr, const char *message)
	FFTJetPileupEstimator ()
	FFTJetPileupEstimator (const FFTJetPileupEstimator &)
FFTJetPileupEstimator &	operator= (const FFTJetPileupEstimator &)

Private Attributes
std::auto_ptr < fftjet::Functor1< double, double > >	calibrationCurve
double	cdfvalue
unsigned	filterNumber
edm::InputTag	inputLabel
std::string	outputLabel
double	ptToDensityFactor
std::auto_ptr < fftjet::Functor1< double, double > >	uncertaintyCurve
std::vector< double >	uncertaintyZones

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
typedef WorkerT< EDProducer >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Description: estimates the actual pileup

Implementation: [Notes on implementation]

Definition at line 43 of file FFTJetPileupEstimator.cc.

Constructor & Destructor Documentation

FFTJetPileupEstimator::FFTJetPileupEstimator ( const edm::ParameterSet &  ps )

explicit

Definition at line 80 of file FFTJetPileupEstimator.cc.

References calibrationCurve, checkConfig(), fftjetcms::fftjet_Function_parser(), edm::ParameterSet::getParameter(), outputLabel, and uncertaintyCurve.

    : init_param(edm::InputTag, inputLabel),
      init_param(std::string, outputLabel),
      init_param(double, cdfvalue),
      init_param(double, ptToDensityFactor),
      init_param(unsigned, filterNumber),
      init_param(std::vector<double>, uncertaintyZones)
{
    calibrationCurve = fftjet_Function_parser(
        ps.getParameter<edm::ParameterSet>("calibrationCurve"));
    checkConfig(calibrationCurve, "bad calibration curve definition");

    uncertaintyCurve = fftjet_Function_parser(
        ps.getParameter<edm::ParameterSet>("uncertaintyCurve"));
    checkConfig(uncertaintyCurve, "bad uncertainty curve definition");

    produces<reco::FFTJetPileupSummary>(outputLabel);
}

FFTJetPileupEstimator::~FFTJetPileupEstimator

(

)

Definition at line 100 of file FFTJetPileupEstimator.cc.

{
}

FFTJetPileupEstimator::FFTJetPileupEstimator ( )

private

FFTJetPileupEstimator::FFTJetPileupEstimator ( const FFTJetPileupEstimator &  )

private

Member Function Documentation

void FFTJetPileupEstimator::beginJob ( void  )

protectedvirtual

Reimplemented from edm::EDProducer.

Definition at line 171 of file FFTJetPileupEstimator.cc.

{
}

template<class Ptr >

void FFTJetPileupEstimator::checkConfig ( const Ptr &  ptr, const char *  message  )

inlineprivate

Definition at line 61 of file FFTJetPileupEstimator.cc.

References NULL.

Referenced by FFTJetPileupEstimator().

    {
        if (ptr.get() == NULL)
            throw cms::Exception("FFTJetBadConfig") << message << std::endl;
    }

void FFTJetPileupEstimator::endJob ( void  )

protectedvirtual

Reimplemented from edm::EDProducer.

Definition at line 176 of file FFTJetPileupEstimator.cc.

{
}

FFTJetPileupEstimator& FFTJetPileupEstimator::operator= ( const FFTJetPileupEstimator &  )

private

void FFTJetPileupEstimator::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

protectedvirtual

Implements edm::EDProducer.

Definition at line 109 of file FFTJetPileupEstimator.cc.

References cdfvalue, edm::hlt::Exception, filterNumber, edm::Event::getByLabel(), h, i, LaserDQM_cfg::input, inputLabel, outputLabel, ptToDensityFactor, edm::Event::put(), edmLumisInFiles::summary, and uncertaintyZones.

{
    edm::Handle<TH2D> input;
    iEvent.getByLabel(inputLabel, input);

    const TH2D& h(*input);
    const unsigned nScales = h.GetXaxis()->GetNbins();
    const unsigned nCdfvalues = h.GetYaxis()->GetNbins();

    const unsigned fixedCdfvalueBin = static_cast<unsigned>(
        std::floor(cdfvalue*nCdfvalues));
    if (fixedCdfvalueBin >= nCdfvalues)
    {
        throw cms::Exception("FFTJetBadConfig") 
            << "Bad cdf value" << std::endl;
    }
    if (filterNumber >= nScales)
    {
        throw cms::Exception("FFTJetBadConfig") 
            << "Bad filter number" << std::endl;
    }

    // Simple fixed-point pile-up estimate
    const double curve = h.GetBinContent(filterNumber+1U,
                                         fixedCdfvalueBin+1U);
    const double pileupRho = ptToDensityFactor*(*calibrationCurve)(curve);
    const double rhoUncert = ptToDensityFactor*(*uncertaintyCurve)(curve);

    // Determine the uncertainty zone of the estimate. The "curve"
    // has to be above or equal to uncertaintyZones[i]  but below
    // uncertaintyZones[i + 1] (the second condition is also satisfied
    // by i == uncertaintyZones.size() - 1). Of course, it is assumed
    // that the vector of zones is configured appropriately -- the zone
    // boundaries must be presented in the increasing order.
    int uncertaintyCode = -1;
    if (!uncertaintyZones.empty())
    {
        const unsigned nZones = uncertaintyZones.size();
        for (unsigned i = 0; i < nZones; ++i)
            if (curve >= uncertaintyZones[i])
            {
                if (i == nZones - 1U)
                {
                    uncertaintyCode = i;
                    break;
                }
                else if (curve < uncertaintyZones[i + 1])
                {
                    uncertaintyCode = i;
                    break;
                }
            }
    }

    std::auto_ptr<reco::FFTJetPileupSummary> summary(
        new reco::FFTJetPileupSummary(curve, pileupRho,
                                      rhoUncert, uncertaintyCode));
    iEvent.put(summary, outputLabel);
}

Member Data Documentation

std::auto_ptr<fftjet::Functor1<double,double> > FFTJetPileupEstimator::calibrationCurve

private

Definition at line 73 of file FFTJetPileupEstimator.cc.

Referenced by FFTJetPileupEstimator().

double FFTJetPileupEstimator::cdfvalue

private

Definition at line 69 of file FFTJetPileupEstimator.cc.

Referenced by produce().

unsigned FFTJetPileupEstimator::filterNumber

private

Definition at line 71 of file FFTJetPileupEstimator.cc.

Referenced by produce().

edm::InputTag FFTJetPileupEstimator::inputLabel

private

Definition at line 67 of file FFTJetPileupEstimator.cc.

Referenced by produce().

std::string FFTJetPileupEstimator::outputLabel

private

Definition at line 68 of file FFTJetPileupEstimator.cc.

Referenced by FFTJetPileupEstimator(), and produce().

double FFTJetPileupEstimator::ptToDensityFactor

private

Definition at line 70 of file FFTJetPileupEstimator.cc.

Referenced by produce().

std::auto_ptr<fftjet::Functor1<double,double> > FFTJetPileupEstimator::uncertaintyCurve

private

Definition at line 74 of file FFTJetPileupEstimator.cc.

Referenced by FFTJetPileupEstimator().

std::vector<double> FFTJetPileupEstimator::uncertaintyZones

private

Definition at line 72 of file FFTJetPileupEstimator.cc.

Referenced by produce().