FixedGridRhoProducerFastjetFromRecHit Class Reference

Inheritance diagram for FixedGridRhoProducerFastjetFromRecHit:

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

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

Private Member Functions
std::array< double, 4 >	getHitP4 (const DetId &detId, const double hitE, const CaloGeometry &caloGeometry) const
bool	passedEcalNoiseCut (const EcalRecHit &hit, const EcalPFRecHitThresholds &thresholds) const
bool	passedHcalNoiseCut (const HBHERecHit &hit) const
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
fastjet::GridMedianBackgroundEstimator	bge_
const edm::ESGetToken < CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
const edm::EDGetTokenT < EcalRecHitCollection >	ebRecHitsTag_
const edm::ESGetToken < EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd >	ecalPFRecHitThresholdsToken_
const edm::EDGetTokenT < EcalRecHitCollection >	eeRecHitsTag_
const EgammaHcalIsolation::arrayHB	eThresHB_
const EgammaHcalIsolation::arrayHE	eThresHE_
const edm::EDGetTokenT < HBHERecHitCollection >	hbheRecHitsTag_
const bool	skipECAL_
const bool	skipHCAL_

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 31 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Constructor & Destructor Documentation

FixedGridRhoProducerFastjetFromRecHit::FixedGridRhoProducerFastjetFromRecHit ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 60 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References DeDxTools::esConsumes().

    : bge_(iConfig.getParameter<double>("maxRapidity"), iConfig.getParameter<double>("gridSpacing")),
      hbheRecHitsTag_(consumes(iConfig.getParameter<edm::InputTag>("hbheRecHitsTag"))),
      ebRecHitsTag_(consumes(iConfig.getParameter<edm::InputTag>("ebRecHitsTag"))),
      eeRecHitsTag_(consumes(iConfig.getParameter<edm::InputTag>("eeRecHitsTag"))),
      eThresHB_(iConfig.getParameter<EgammaHcalIsolation::arrayHB>("eThresHB")),
      eThresHE_(iConfig.getParameter<EgammaHcalIsolation::arrayHE>("eThresHE")),
      skipHCAL_(iConfig.getParameter<bool>("skipHCAL")),
      skipECAL_(iConfig.getParameter<bool>("skipECAL")),
      ecalPFRecHitThresholdsToken_{esConsumes()},
      caloGeometryToken_{esConsumes()} {
  if (skipHCAL_ && skipECAL_) {
    throw cms::Exception("FixedGridRhoProducerFastjetFromRecHit")
        << "skipHCAL and skipECAL both can't be True. Please make at least one of them False.";
  }
  produces<double>();
}

FixedGridRhoProducerFastjetFromRecHit::~FixedGridRhoProducerFastjetFromRecHit ( )

overridedefault

Member Function Documentation

void FixedGridRhoProducerFastjetFromRecHit::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 78 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

                                                                                                       {
  edm::ParameterSetDescription desc;
  //We use raw recHits and not the PF recHits, because in Phase1 muon and egamma HLT,
  //PF recHit collection is regional, while the full raw recHit collection is available.
  desc.add<edm::InputTag>("hbheRecHitsTag", edm::InputTag("hltHbhereco"));
  desc.add<edm::InputTag>("ebRecHitsTag", edm::InputTag("hltEcalRecHit", "EcalRecHitsEB"));
  desc.add<edm::InputTag>("eeRecHitsTag", edm::InputTag("hltEcalRecHit", "EcalRecHitsEE"));
  desc.add<bool>("skipHCAL", false);
  desc.add<bool>("skipECAL", false);
  //eThresHB/HE are from RecoParticleFlow/PFClusterProducer/python/particleFlowRecHitHBHE_cfi.py
  desc.add<std::vector<double> >("eThresHB", {0.1, 0.2, 0.3, 0.3});
  desc.add<std::vector<double> >("eThresHE", {0.1, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2});
  desc.add<double>("maxRapidity", 2.5);
  desc.add<double>("gridSpacing", 0.55);
  descriptions.addWithDefaultLabel(desc);
}

std::array< double, 4 > FixedGridRhoProducerFastjetFromRecHit::getHitP4 ( const DetId &  detId, const double  hitE, const CaloGeometry &  caloGeometry  ) const

private

Definition at line 137 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References funct::cos(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), funct::sin(), and GeomDetEnumerators::subDetGeom.

Referenced by produce().

                                                                                                              {
  const CaloSubdetectorGeometry *subDetGeom = caloGeometry.getSubdetectorGeometry(detId);
  const auto &gpPos = subDetGeom->getGeometry(detId)->repPos();
  const double thispt = hitE / cosh(gpPos.eta());
  const double thispx = thispt * cos(gpPos.phi());
  const double thispy = thispt * sin(gpPos.phi());
  const double thispz = thispt * sinh(gpPos.eta());
  return std::array<double, 4>{{thispx, thispy, thispz, hitE}};
}

bool FixedGridRhoProducerFastjetFromRecHit::passedEcalNoiseCut ( const EcalRecHit &  hit, const EcalPFRecHitThresholds &  thresholds  ) const

private

Definition at line 158 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References EcalRecHit::detid(), and EcalRecHit::energy().

Referenced by produce().

                                                                                                               {
  return (hit.energy() > thresholds[hit.detid()]);
}

bool FixedGridRhoProducerFastjetFromRecHit::passedHcalNoiseCut ( const HBHERecHit &  hit ) const

private

Definition at line 150 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References HcalDetId::depth(), CaloRecHit::energy(), eThresHB_, eThresHE_, HcalBarrel, HcalEndcap, and HBHERecHit::id().

Referenced by produce().

                                                                                          {
  const auto thisDetId = hit.id();
  const auto thisDepth = thisDetId.depth();
  return (thisDetId.subdet() == HcalBarrel && hit.energy() > eThresHB_[thisDepth - 1]) ||
         (thisDetId.subdet() == HcalEndcap && hit.energy() > eThresHE_[thisDepth - 1]);
}

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

overrideprivate

Definition at line 97 of file FixedGridRhoProducerFastjetFromRecHit.cc.

References bge_, caloGeometryToken_, ebRecHitsTag_, ecalPFRecHitThresholdsToken_, eeRecHitsTag_, edm::Event::get(), edm::EventSetup::getData(), getHitP4(), ootPhotons_cff::hbheRecHits, hbheRecHitsTag_, hit::id, PixelMapPlotter::inputs, passedEcalNoiseCut(), passedHcalNoiseCut(), edm::Event::put(), skipECAL_, skipHCAL_, and lowPtGsfElectronSeeds_cfi::thresholds().

                                                                                                 {
  std::vector<fastjet::PseudoJet> inputs;
  auto const &thresholds = iSetup.getData(ecalPFRecHitThresholdsToken_);
  auto const &caloGeometry = iSetup.getData(caloGeometryToken_);

  if (!skipHCAL_) {
    auto const &hbheRecHits = iEvent.get(hbheRecHitsTag_);
    inputs.reserve(inputs.size() + hbheRecHits.size());
    for (const auto &hit : hbheRecHits) {
      if (passedHcalNoiseCut(hit)) {
        const auto &hitp4 = getHitP4(hit.id(), hit.energy(), caloGeometry);
        inputs.emplace_back(fastjet::PseudoJet(hitp4[0], hitp4[1], hitp4[2], hitp4[3]));
      }
    }
  }

  if (!skipECAL_) {
    auto const &ebRecHits = iEvent.get(ebRecHitsTag_);
    inputs.reserve(inputs.size() + ebRecHits.size());
    for (const auto &hit : ebRecHits) {
      if (passedEcalNoiseCut(hit, thresholds)) {
        const auto &hitp4 = getHitP4(hit.id(), hit.energy(), caloGeometry);
        inputs.emplace_back(fastjet::PseudoJet(hitp4[0], hitp4[1], hitp4[2], hitp4[3]));
      }
    }

    auto const &eeRecHits = iEvent.get(eeRecHitsTag_);
    inputs.reserve(inputs.size() + eeRecHits.size());
    for (const auto &hit : eeRecHits) {
      if (passedEcalNoiseCut(hit, thresholds)) {
        const auto &hitp4 = getHitP4(hit.id(), hit.energy(), caloGeometry);
        inputs.emplace_back(fastjet::PseudoJet(hitp4[0], hitp4[1], hitp4[2], hitp4[3]));
      }
    }
  }

  bge_.set_particles(inputs);
  iEvent.put(std::make_unique<double>(bge_.rho()));
}

Member Data Documentation

fastjet::GridMedianBackgroundEstimator FixedGridRhoProducerFastjetFromRecHit::bge_

private

Definition at line 43 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> FixedGridRhoProducerFastjetFromRecHit::caloGeometryToken_

private

Definition at line 57 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().

const edm::EDGetTokenT<EcalRecHitCollection> FixedGridRhoProducerFastjetFromRecHit::ebRecHitsTag_

private

Definition at line 45 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().

const edm::ESGetToken<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd> FixedGridRhoProducerFastjetFromRecHit::ecalPFRecHitThresholdsToken_

private

Definition at line 56 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().

const edm::EDGetTokenT<EcalRecHitCollection> FixedGridRhoProducerFastjetFromRecHit::eeRecHitsTag_

private

Definition at line 46 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().

const EgammaHcalIsolation::arrayHB FixedGridRhoProducerFastjetFromRecHit::eThresHB_

private

Definition at line 48 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by passedHcalNoiseCut().

const EgammaHcalIsolation::arrayHE FixedGridRhoProducerFastjetFromRecHit::eThresHE_

private

Definition at line 49 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by passedHcalNoiseCut().

const edm::EDGetTokenT<HBHERecHitCollection> FixedGridRhoProducerFastjetFromRecHit::hbheRecHitsTag_

private

Definition at line 44 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().

const bool FixedGridRhoProducerFastjetFromRecHit::skipECAL_

private

Definition at line 54 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().

const bool FixedGridRhoProducerFastjetFromRecHit::skipHCAL_

private

Definition at line 53 of file FixedGridRhoProducerFastjetFromRecHit.cc.

Referenced by produce().