HLTJetTimingProducer< T > Class Template Reference

This produces timing and associated ecal cell information for calo jets. More...

#include <HLTJetTimingProducer.h>

Inheritance diagram for HLTJetTimingProducer< T >:

Public Member Functions
	HLTJetTimingProducer (const edm::ParameterSet &)
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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	jetTimeFromEcalCells (const T &, const edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit >> &, const CaloGeometry &, float &, float &, uint &)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
const bool	barrelJets_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
const double	ecalCellEnergyThresh_
const double	ecalCellTimeErrorThresh_
const double	ecalCellTimeThresh_
const edm::EDGetTokenT< edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit > > >	ecalRecHitsEBToken_
const edm::EDGetTokenT< edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit > > >	ecalRecHitsEEToken_
const bool	endcapJets_
const edm::EDGetTokenT< std::vector< T > >	jetInputToken_
const double	matchingRadius2_

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

template<typename T>
class HLTJetTimingProducer< T >

This produces timing and associated ecal cell information for calo jets.

Author

Matthew Citron

Definition at line 28 of file HLTJetTimingProducer.h.

Constructor & Destructor Documentation

HLTJetTimingProducer()

template<typename T >

HLTJetTimingProducer< T >::HLTJetTimingProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 60 of file HLTJetTimingProducer.h.

References edm::ParameterSet::getParameter().

    : caloGeometryToken_(esConsumes()),
      jetInputToken_{consumes<std::vector<T>>(iConfig.getParameter<edm::InputTag>("jets"))},
      ecalRecHitsEBToken_{consumes<edm::SortedCollection<EcalRecHit, edm::StrictWeakOrdering<EcalRecHit>>>(
          iConfig.getParameter<edm::InputTag>("ebRecHitsColl"))},
      ecalRecHitsEEToken_{consumes<edm::SortedCollection<EcalRecHit, edm::StrictWeakOrdering<EcalRecHit>>>(
          iConfig.getParameter<edm::InputTag>("eeRecHitsColl"))},
      barrelJets_{iConfig.getParameter<bool>("barrelJets")},
      endcapJets_{iConfig.getParameter<bool>("endcapJets")},
      ecalCellEnergyThresh_{iConfig.getParameter<double>("ecalCellEnergyThresh")},
      ecalCellTimeThresh_{iConfig.getParameter<double>("ecalCellTimeThresh")},
      ecalCellTimeErrorThresh_{iConfig.getParameter<double>("ecalCellTimeErrorThresh")},
      matchingRadius2_{std::pow(iConfig.getParameter<double>("matchingRadius"), 2)} {
  produces<edm::ValueMap<float>>("");
  produces<edm::ValueMap<unsigned int>>("jetCellsForTiming");
  produces<edm::ValueMap<float>>("jetEcalEtForTiming");
}

Member Function Documentation

fillDescriptions()

template<typename T >

void HLTJetTimingProducer< T >::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 161 of file HLTJetTimingProducer.h.

References edm::ConfigurationDescriptions::add(), defaultModuleLabel(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                         {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("jets", edm::InputTag(""));
  desc.add<bool>("barrelJets", false);
  desc.add<bool>("endcapJets", false);
  desc.add<double>("ecalCellEnergyThresh", 0.5);
  desc.add<double>("ecalCellTimeThresh", 12.5);
  desc.add<double>("ecalCellTimeErrorThresh", 100.);
  desc.add<double>("matchingRadius", 0.4);
  desc.add<edm::InputTag>("ebRecHitsColl", edm::InputTag("hltEcalRecHit", "EcalRecHitsEB"));
  desc.add<edm::InputTag>("eeRecHitsColl", edm::InputTag("hltEcalRecHit", "EcalRecHitsEE"));
  descriptions.add(defaultModuleLabel<HLTJetTimingProducer<T>>(), desc);
}

jetTimeFromEcalCells()

template<typename T >

void HLTJetTimingProducer< T >::jetTimeFromEcalCells ( const T &  jet, const edm::SortedCollection< EcalRecHit, edm::StrictWeakOrdering< EcalRecHit >> &  ecalRecHits, const CaloGeometry &  caloGeometry, float &  weightedTimeCell, float &  totalEmEnergyCell, uint &  nCells  )

private

Definition at line 79 of file HLTJetTimingProducer.h.

References funct::abs(), reco::deltaR2(), CaloGeometry::getPosition(), metsig::jet, EcalRecHit::kDiWeird, EcalRecHit::kLeadingEdgeRecovered, EcalRecHit::kPoorReco, EcalRecHit::kSaturated, EcalRecHit::kWeird, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::nCells, and funct::sin().

                  {
  for (auto const& ecalRH : ecalRecHits) {
    if (ecalRH.checkFlag(EcalRecHit::kSaturated) || ecalRH.checkFlag(EcalRecHit::kLeadingEdgeRecovered) ||
        ecalRH.checkFlag(EcalRecHit::kPoorReco) || ecalRH.checkFlag(EcalRecHit::kWeird) ||
        ecalRH.checkFlag(EcalRecHit::kDiWeird))
      continue;
    if (ecalRH.energy() < ecalCellEnergyThresh_)
      continue;
    if (ecalRH.timeError() <= 0. || ecalRH.timeError() > ecalCellTimeErrorThresh_)
      continue;
    if (std::abs(ecalRH.time()) > ecalCellTimeThresh_)
      continue;
    auto const pos = caloGeometry.getPosition(ecalRH.detid());
    if (reco::deltaR2(jet, pos) > matchingRadius2_)
      continue;
    weightedTimeCell += ecalRH.time() * ecalRH.energy() * sin(pos.theta());
    totalEmEnergyCell += ecalRH.energy() * sin(pos.theta());
    nCells++;
  }
  if (totalEmEnergyCell > 0) {
    weightedTimeCell /= totalEmEnergyCell;
  }
}

produce()

template<typename T >

void HLTJetTimingProducer< T >::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 110 of file HLTJetTimingProducer.h.

References HLT_2023v12_cff::ecalRecHitsEB, HLT_2023v12_cff::ecalRecHitsEE, edm::helper::Filler< Map >::fill(), edm::EventSetup::getData(), iEvent, edm::helper::Filler< Map >::insert(), metsig::jet, PDWG_EXODelayedJetMET_cff::jetCellsForTiming, PDWG_EXODelayedJetMET_cff::jetEcalEtForTiming, PDWG_EXODelayedJetMET_cff::jets, eostools::move(), and ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::nCells.

                                                                                   {
  auto const& caloGeometry = iSetup.getData(caloGeometryToken_);
  auto const jets = iEvent.getHandle(jetInputToken_);
  auto const& ecalRecHitsEB = iEvent.get(ecalRecHitsEBToken_);
  auto const& ecalRecHitsEE = iEvent.get(ecalRecHitsEEToken_);

  std::vector<float> jetTimings;
  std::vector<unsigned int> jetCellsForTiming;
  std::vector<float> jetEcalEtForTiming;

  jetTimings.reserve(jets->size());
  jetEcalEtForTiming.reserve(jets->size());
  jetCellsForTiming.reserve(jets->size());

  for (auto const& jet : *jets) {
    float weightedTimeCell = 0;
    float totalEmEnergyCell = 0;
    unsigned int nCells = 0;
    if (barrelJets_)
      jetTimeFromEcalCells(jet, ecalRecHitsEB, caloGeometry, weightedTimeCell, totalEmEnergyCell, nCells);
    if (endcapJets_) {
      weightedTimeCell *= totalEmEnergyCell;
      jetTimeFromEcalCells(jet, ecalRecHitsEE, caloGeometry, weightedTimeCell, totalEmEnergyCell, nCells);
    }

    // If there is at least one ecal cell passing selection, calculate timing
    jetTimings.emplace_back(totalEmEnergyCell > 0 ? weightedTimeCell : -50);
    jetEcalEtForTiming.emplace_back(totalEmEnergyCell);
    jetCellsForTiming.emplace_back(nCells);
  }

  std::unique_ptr<edm::ValueMap<float>> jetTimings_out(new edm::ValueMap<float>());
  edm::ValueMap<float>::Filler jetTimings_filler(*jetTimings_out);
  jetTimings_filler.insert(jets, jetTimings.begin(), jetTimings.end());
  jetTimings_filler.fill();
  iEvent.put(std::move(jetTimings_out), "");

  std::unique_ptr<edm::ValueMap<float>> jetEcalEtForTiming_out(new edm::ValueMap<float>());
  edm::ValueMap<float>::Filler jetEcalEtForTiming_filler(*jetEcalEtForTiming_out);
  jetEcalEtForTiming_filler.insert(jets, jetEcalEtForTiming.begin(), jetEcalEtForTiming.end());
  jetEcalEtForTiming_filler.fill();
  iEvent.put(std::move(jetEcalEtForTiming_out), "jetEcalEtForTiming");

  std::unique_ptr<edm::ValueMap<unsigned int>> jetCellsForTiming_out(new edm::ValueMap<unsigned int>());
  edm::ValueMap<unsigned int>::Filler jetCellsForTiming_filler(*jetCellsForTiming_out);
  jetCellsForTiming_filler.insert(jets, jetCellsForTiming.begin(), jetCellsForTiming.end());
  jetCellsForTiming_filler.fill();
  iEvent.put(std::move(jetCellsForTiming_out), "jetCellsForTiming");
}

Member Data Documentation

barrelJets_

template<typename T >

const bool HLTJetTimingProducer< T >::barrelJets_

private

Definition at line 49 of file HLTJetTimingProducer.h.

caloGeometryToken_

template<typename T >

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HLTJetTimingProducer< T >::caloGeometryToken_

private

Definition at line 42 of file HLTJetTimingProducer.h.

ecalCellEnergyThresh_

template<typename T >

const double HLTJetTimingProducer< T >::ecalCellEnergyThresh_

private

Definition at line 53 of file HLTJetTimingProducer.h.

ecalCellTimeErrorThresh_

template<typename T >

const double HLTJetTimingProducer< T >::ecalCellTimeErrorThresh_

private

Definition at line 55 of file HLTJetTimingProducer.h.

ecalCellTimeThresh_

template<typename T >

const double HLTJetTimingProducer< T >::ecalCellTimeThresh_

private

Definition at line 54 of file HLTJetTimingProducer.h.

ecalRecHitsEBToken_

template<typename T >

const edm::EDGetTokenT<edm::SortedCollection<EcalRecHit, edm::StrictWeakOrdering<EcalRecHit> > > HLTJetTimingProducer< T >::ecalRecHitsEBToken_

private

Definition at line 45 of file HLTJetTimingProducer.h.

ecalRecHitsEEToken_

template<typename T >

const edm::EDGetTokenT<edm::SortedCollection<EcalRecHit, edm::StrictWeakOrdering<EcalRecHit> > > HLTJetTimingProducer< T >::ecalRecHitsEEToken_

private

Definition at line 46 of file HLTJetTimingProducer.h.

endcapJets_

template<typename T >

const bool HLTJetTimingProducer< T >::endcapJets_

private

Definition at line 50 of file HLTJetTimingProducer.h.

jetInputToken_

template<typename T >

const edm::EDGetTokenT<std::vector<T> > HLTJetTimingProducer< T >::jetInputToken_

private

Definition at line 44 of file HLTJetTimingProducer.h.

matchingRadius2_

template<typename T >

const double HLTJetTimingProducer< T >::matchingRadius2_

private

Definition at line 56 of file HLTJetTimingProducer.h.