fftjetcms::FFTJetInterface Class Reference

#include <FFTJetInterface.h>

Inheritance diagram for fftjetcms::FFTJetInterface:

Public Member Functions
	FFTJetInterface ()=delete
	FFTJetInterface (const FFTJetInterface &)=delete
FFTJetInterface &	operator= (const FFTJetInterface &)=delete
	~FFTJetInterface () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
template<class Ptr >
void	checkConfig (const Ptr &ptr, const char *message)
void	discretizeEnergyFlow ()
	FFTJetInterface (const edm::ParameterSet &)
double	getEventScale () const
void	loadInputCollection (const edm::Event &)
bool	storeInSinglePrecision () const
const reco::Particle::Point &	vertexUsed () const

Protected Attributes
const AnomalousTower	anomalous
std::vector< unsigned >	candidateIndex
const bool	doPVCorrection
std::unique_ptr< fftjet::Grid2d< fftjetcms::Real > >	energyFlow
const std::vector< double >	etaDependentMagnutideFactors
std::vector< fftjetcms::VectorLike >	eventData
edm::Handle< reco::CandidateView >	inputCollection
const edm::InputTag	inputLabel
const JetType	jetType
const std::string	outputLabel
const edm::InputTag	srcPVs

Private Attributes
const double	completeEventScale
edm::EDGetTokenT< reco::CandidateView >	inputToken
const bool	insertCompleteEvent
edm::EDGetTokenT< reco::VertexCollection >	srcPVsToken
reco::Particle::Point	vertex_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 52 of file FFTJetInterface.h.

Constructor & Destructor Documentation

FFTJetInterface() [1/3]

fftjetcms::FFTJetInterface::FFTJetInterface ( )

delete

FFTJetInterface() [2/3]

fftjetcms::FFTJetInterface::FFTJetInterface ( const FFTJetInterface &  )

delete

~FFTJetInterface()

fftjetcms::FFTJetInterface::~FFTJetInterface ( )

inlineoverride

Definition at line 58 of file FFTJetInterface.h.

{}

FFTJetInterface() [3/3]

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

explicitprotected

Definition at line 17 of file FFTJetInterface.cc.

References fftjetcms::CALOJET, completeEventScale, doPVCorrection, Exception, inputLabel, inputToken, insertCompleteEvent, jetType, srcPVs, and srcPVsToken.

      : inputLabel(ps.getParameter<edm::InputTag>("src")),
        init_param(std::string, outputLabel),
        jetType(parseJetType(ps.getParameter<std::string>("jetType"))),
        init_param(bool, doPVCorrection),
        init_param(edm::InputTag, srcPVs),
        etaDependentMagnutideFactors(ps.getParameter<std::vector<double> >("etaDependentMagnutideFactors")),
        init_param(edm::ParameterSet, anomalous),
        init_param(bool, insertCompleteEvent),
        init_param(double, completeEventScale),
        vertex_(0.0, 0.0, 0.0) {
    if (insertCompleteEvent && completeEventScale <= 0.0)
      throw cms::Exception("FFTJetBadConfig") << "Bad scale for the complete event : must be positive" << std::endl;

    inputToken = consumes<reco::CandidateView>(inputLabel);

    if (doPVCorrection && jetType == CALOJET)
      srcPVsToken = consumes<reco::VertexCollection>(srcPVs);
  }

Member Function Documentation

checkConfig()

template<class Ptr >

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

inlineprotected

Definition at line 64 of file FFTJetInterface.h.

Referenced by FFTJetProducer::beginStream(), FFTJetEFlowSmoother::FFTJetEFlowSmoother(), FFTJetPatRecoProducer::FFTJetPatRecoProducer(), FFTJetPileupProcessor::FFTJetPileupProcessor(), and FFTJetProducer::FFTJetProducer().

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

discretizeEnergyFlow()

void FFTJetInterface::discretizeEnergyFlow ( )

protected

Definition at line 78 of file FFTJetInterface.cc.

References cms::cuda::assert(), energyFlow, etaDependentMagnutideFactors, eventData, g, mps_fire::i, and B2GTnPMonitor_cfi::item.

Referenced by FFTJetEFlowSmoother::produce(), FFTJetPileupProcessor::produce(), FFTJetPatRecoProducer::produce(), and FFTJetProducer::produce().

                                             {
    // It is a bug to call this function before defining the energy flow grid
    assert(energyFlow.get());

    fftjet::Grid2d<Real>& g(*energyFlow);
    g.reset();

    const unsigned nInputs = eventData.size();
    const VectorLike* inp = nInputs ? &eventData[0] : nullptr;
    for (unsigned i = 0; i < nInputs; ++i) {
      const VectorLike& item(inp[i]);
      g.fill(item.Eta(), item.Phi(), item.Et());
    }

    if (!etaDependentMagnutideFactors.empty()) {
      if (etaDependentMagnutideFactors.size() != g.nEta())
        throw cms::Exception("FFTJetBadConfig") << "ERROR in FFTJetInterface::discretizeEnergyFlow() :"
                                                   " number of elements in the \"etaDependentMagnutideFactors\""
                                                   " vector is inconsistent with the grid binning"
                                                << std::endl;
      g.scaleData(&etaDependentMagnutideFactors[0], etaDependentMagnutideFactors.size());
    }
  }

getEventScale()

double FFTJetInterface::getEventScale ( ) const

protected

Definition at line 37 of file FFTJetInterface.cc.

References completeEventScale, and insertCompleteEvent.

Referenced by FFTJetPatRecoProducer::buildDenseProduct(), FFTJetPatRecoProducer::buildSparseProduct(), FFTJetPatRecoProducer::FFTJetPatRecoProducer(), FFTJetProducer::loadSparseTreeData(), FFTJetEFlowSmoother::produce(), and FFTJetPatRecoProducer::produce().

{ return insertCompleteEvent ? completeEventScale : 0.0; }

loadInputCollection()

void FFTJetInterface::loadInputCollection ( const edm::Event &  iEvent )

protected

Definition at line 39 of file FFTJetInterface.cc.

References anomalous, cms::cuda::assert(), fftjetcms::CALOJET, candidateIndex, doPVCorrection, eventData, iEvent, inputCollection, inputToken, edm::isNotFinite(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, B2GTnPMonitor_cfi::item, jetType, heavyFlavorDQMFirstStep_cff::pvCollection, srcPVsToken, l1tHGCalTowerProducer_cfi::tower, and vertex_.

Referenced by FFTJetEFlowSmoother::produce(), FFTJetPileupProcessor::produce(), FFTJetPatRecoProducer::produce(), and FFTJetProducer::produce().

                                                                  {
    // Figure out if we are going to perform the vertex adjustment
    const bool adjustForVertex = doPVCorrection && jetType == CALOJET;

    // Figure out the vertex
    if (adjustForVertex) {
      edm::Handle<reco::VertexCollection> pvCollection;
      iEvent.getByToken(srcPVsToken, pvCollection);
      if (pvCollection->empty())
        vertex_ = reco::Particle::Point(0.0, 0.0, 0.0);
      else
        vertex_ = pvCollection->begin()->position();
    }

    // Get the input collection
    iEvent.getByToken(inputToken, inputCollection);

    // Create the set of 4-vectors needed by the algorithm
    eventData.clear();
    candidateIndex.clear();
    unsigned index = 0;
    const reco::CandidateView::const_iterator end = inputCollection->end();
    for (reco::CandidateView::const_iterator it = inputCollection->begin(); it != end; ++it, ++index) {
      const reco::Candidate& item(*it);
      if (anomalous(item))
        continue;
      if (edm::isNotFinite(item.pt()))
        continue;

      if (adjustForVertex) {
        const CaloTower& tower(dynamic_cast<const CaloTower&>(item));
        eventData.push_back(VectorLike(tower.p4(vertex_)));
      } else
        eventData.push_back(item.p4());
      candidateIndex.push_back(index);
    }
    assert(eventData.size() == candidateIndex.size());
  }

operator=()

FFTJetInterface& fftjetcms::FFTJetInterface::operator= ( const FFTJetInterface &  )

delete

storeInSinglePrecision()

bool FFTJetInterface::storeInSinglePrecision ( ) const

protected

Definition at line 15 of file FFTJetInterface.cc.

Referenced by FFTJetPatRecoProducer::buildDenseProduct(), FFTJetPatRecoProducer::buildSparseProduct(), FFTJetPatRecoProducer::FFTJetPatRecoProducer(), FFTJetProducer::FFTJetProducer(), FFTJetPatRecoProducer::produce(), and FFTJetProducer::produce().

{ return true; }

vertexUsed()

const reco::Particle::Point& fftjetcms::FFTJetInterface::vertexUsed ( ) const

inlineprotected

Definition at line 74 of file FFTJetInterface.h.

References vertex_.

Referenced by FFTJetProducer::writeJets().

{ return vertex_; }

Member Data Documentation

anomalous

const AnomalousTower fftjetcms::FFTJetInterface::anomalous

protected

Definition at line 95 of file FFTJetInterface.h.

Referenced by loadInputCollection().

candidateIndex

std::vector<unsigned> fftjetcms::FFTJetInterface::candidateIndex

protected

Definition at line 101 of file FFTJetInterface.h.

Referenced by FFTJetProducer::determineGriddedConstituents(), FFTJetProducer::determineVectorConstituents(), and loadInputCollection().

completeEventScale

const double fftjetcms::FFTJetInterface::completeEventScale

private

Definition at line 111 of file FFTJetInterface.h.

Referenced by FFTJetInterface(), and getEventScale().

doPVCorrection

const bool fftjetcms::FFTJetInterface::doPVCorrection

protected

Definition at line 86 of file FFTJetInterface.h.

Referenced by FFTJetInterface(), and loadInputCollection().

energyFlow

std::unique_ptr<fftjet::Grid2d<fftjetcms::Real> > fftjetcms::FFTJetInterface::energyFlow

protected

Definition at line 104 of file FFTJetInterface.h.

Referenced by FFTJetProducer::beginStream(), FFTJetProducer::buildGridAlg(), FFTJetEFlowSmoother::buildKernelConvolver(), FFTJetPileupProcessor::buildKernelConvolver(), FFTJetPatRecoProducer::buildKernelConvolver(), FFTJetPatRecoProducer::buildPeakFinder(), FFTJetProducer::determineGriddedConstituents(), discretizeEnergyFlow(), FFTJetEFlowSmoother::FFTJetEFlowSmoother(), FFTJetPatRecoProducer::FFTJetPatRecoProducer(), FFTJetPileupProcessor::FFTJetPileupProcessor(), FFTJetProducer::iterateJetReconstruction(), FFTJetPileupProcessor::loadFlatteningFactors(), FFTJetPileupProcessor::mixExtraGrid(), FFTJetEFlowSmoother::produce(), FFTJetPileupProcessor::produce(), FFTJetPatRecoProducer::produce(), FFTJetProducer::produce(), and FFTJetProducer::writeJets().

etaDependentMagnutideFactors

const std::vector<double> fftjetcms::FFTJetInterface::etaDependentMagnutideFactors

protected

Definition at line 92 of file FFTJetInterface.h.

Referenced by discretizeEnergyFlow().

eventData

std::vector<fftjetcms::VectorLike> fftjetcms::FFTJetInterface::eventData

protected

Definition at line 98 of file FFTJetInterface.h.

Referenced by FFTJetProducer::determineGriddedConstituents(), FFTJetProducer::determineVectorConstituents(), discretizeEnergyFlow(), FFTJetProducer::iterateJetReconstruction(), loadInputCollection(), and FFTJetProducer::produce().

inputCollection

edm::Handle<reco::CandidateView> fftjetcms::FFTJetInterface::inputCollection

protected

Definition at line 107 of file FFTJetInterface.h.

Referenced by FFTJetProducer::determineGriddedConstituents(), FFTJetProducer::determineVectorConstituents(), and loadInputCollection().

inputLabel

const edm::InputTag fftjetcms::FFTJetInterface::inputLabel

protected

Definition at line 77 of file FFTJetInterface.h.

Referenced by FFTJetInterface().

inputToken

edm::EDGetTokenT<reco::CandidateView> fftjetcms::FFTJetInterface::inputToken

private

Definition at line 114 of file FFTJetInterface.h.

Referenced by FFTJetInterface(), and loadInputCollection().

insertCompleteEvent

const bool fftjetcms::FFTJetInterface::insertCompleteEvent

private

Definition at line 110 of file FFTJetInterface.h.

Referenced by FFTJetInterface(), and getEventScale().

jetType

const JetType fftjetcms::FFTJetInterface::jetType

protected

Definition at line 83 of file FFTJetInterface.h.

Referenced by FFTJetInterface(), FFTJetProducer::FFTJetProducer(), and loadInputCollection().

outputLabel

const std::string fftjetcms::FFTJetInterface::outputLabel

protected

Definition at line 80 of file FFTJetInterface.h.

Referenced by FFTJetPatRecoProducer::buildDenseProduct(), FFTJetPatRecoProducer::buildSparseProduct(), FFTJetEFlowSmoother::FFTJetEFlowSmoother(), FFTJetPatRecoProducer::FFTJetPatRecoProducer(), FFTJetPileupProcessor::FFTJetPileupProcessor(), FFTJetProducer::FFTJetProducer(), FFTJetEFlowSmoother::produce(), FFTJetPileupProcessor::produce(), FFTJetPatRecoProducer::produce(), FFTJetProducer::saveResults(), and FFTJetProducer::writeJets().

srcPVs

const edm::InputTag fftjetcms::FFTJetInterface::srcPVs

protected

Definition at line 89 of file FFTJetInterface.h.

Referenced by FFTJetInterface().

srcPVsToken

edm::EDGetTokenT<reco::VertexCollection> fftjetcms::FFTJetInterface::srcPVsToken

private

Definition at line 115 of file FFTJetInterface.h.

Referenced by FFTJetInterface(), and loadInputCollection().

vertex_

reco::Particle::Point fftjetcms::FFTJetInterface::vertex_

private

Definition at line 112 of file FFTJetInterface.h.

Referenced by loadInputCollection(), and vertexUsed().