PFHGCalRecHitCreator< DET, Layer, det, subdet > Class Template Reference

#include <PFHGCalRecHitCreator.h>

Inheritance diagram for PFHGCalRecHitCreator< DET, Layer, det, subdet >:

Public Member Functions
void	importRecHits (std::unique_ptr< reco::PFRecHitCollection > &out, std::unique_ptr< reco::PFRecHitCollection > &cleaned, const edm::Event &iEvent, const edm::EventSetup &iSetup) override
	PFHGCalRecHitCreator (const edm::ParameterSet &iConfig, edm::ConsumesCollector &cc)
Public Member Functions inherited from PFRecHitCreatorBase
virtual void	init (const edm::EventSetup &es)
	PFRecHitCreatorBase ()
	PFRecHitCreatorBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &cc)
virtual	~PFRecHitCreatorBase ()=default

Protected Attributes
std::string	geometryInstance_
edm::EDGetTokenT< HGCRecHitCollection >	recHitToken_
Protected Attributes inherited from PFRecHitCreatorBase
std::vector< std::unique_ptr< PFRecHitQTestBase > >	qualityTests_

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomToken_
hgcal::RecHitTools	recHitTools_

Additional Inherited Members
Protected Member Functions inherited from PFRecHitCreatorBase
void	beginEvent (const edm::Event &event, const edm::EventSetup &setup)

Detailed Description

template<typename DET, PFLayer::Layer Layer, DetId::Detector det, unsigned subdet>
class PFHGCalRecHitCreator< DET, Layer, det, subdet >

Definition at line 26 of file PFHGCalRecHitCreator.h.

Constructor & Destructor Documentation

PFHGCalRecHitCreator()

template<typename DET , PFLayer::Layer Layer, DetId::Detector det, unsigned subdet>

PFHGCalRecHitCreator< DET, Layer, det, subdet >::PFHGCalRecHitCreator ( const edm::ParameterSet &  iConfig, edm::ConsumesCollector &  cc  )

inline

Definition at line 28 of file PFHGCalRecHitCreator.h.

      : PFRecHitCreatorBase(iConfig, cc),
        recHitToken_(cc.consumes<HGCRecHitCollection>(iConfig.getParameter<edm::InputTag>("src"))),
        geometryInstance_(iConfig.getParameter<std::string>("geometryInstance")),
        geomToken_(cc.esConsumes()) {}

Member Function Documentation

importRecHits()

template<typename DET , PFLayer::Layer Layer, DetId::Detector det, unsigned subdet>

void PFHGCalRecHitCreator< DET, Layer, det, subdet >::importRecHits ( std::unique_ptr< reco::PFRecHitCollection > &  out, std::unique_ptr< reco::PFRecHitCollection > &  cleaned, const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

inlineoverridevirtual

Implements PFRecHitCreatorBase.

Definition at line 34 of file PFHGCalRecHitCreator.h.

References HCALHighEnergyHPDFilter_cfi::energy, Exception, relativeConstraints::geom, PFHGCalRecHitCreator< DET, Layer, det, subdet >::geomToken_, edm::EventSetup::getHandle(), mps_fire::i, iEvent, runTheMatrix::keep, LogDebug, or, MillePedeFileConverter_cfg::out, edm::ESHandle< T >::product(), PFRecHitCreatorBase::qualityTests_, HI_PhotonSkim_cff::rechits, PFHGCalRecHitCreator< DET, Layer, det, subdet >::recHitToken_, PFHGCalRecHitCreator< DET, Layer, det, subdet >::recHitTools_, hgcal::RecHitTools::setGeometry(), and protons_cff::time.

                                                           {
    // Setup RecHitTools to properly compute the position of the HGCAL Cells vie their DetIds
    edm::ESHandle<CaloGeometry> geoHandle = iSetup.getHandle(geomToken_);
    recHitTools_.setGeometry(*geoHandle);

    for (unsigned int i = 0; i < qualityTests_.size(); ++i) {
      qualityTests_.at(i)->beginEvent(iEvent, iSetup);
    }

    edm::Handle<HGCRecHitCollection> recHitHandle;
    iEvent.getByToken(recHitToken_, recHitHandle);
    const HGCRecHitCollection& rechits = *recHitHandle;

    const CaloGeometry* geom = geoHandle.product();

    unsigned skipped_rechits = 0;
    for (const auto& hgrh : rechits) {
      const DET detid(hgrh.detid());

      if (det != detid.det() or (subdet != 0 and subdet != detid.subdetId())) {
        throw cms::Exception("IncorrectHGCSubdetector")
            << "det expected: " << det << " det gotten: " << detid.det() << " ; "
            << "subdet expected: " << subdet << " subdet gotten: " << detid.subdetId() << std::endl;
      }

      double energy = hgrh.energy();
      double time = hgrh.time();

      auto thisCell = geom->getSubdetectorGeometry(det, subdet)->getGeometry(detid);

      // find rechit geometry
      if (!thisCell) {
        LogDebug("PFHGCalRecHitCreator") << "warning detid " << detid.rawId() << " not found in geometry" << std::endl;
        ++skipped_rechits;
        continue;
      }

      reco::PFRecHit rh(thisCell, detid.rawId(), Layer, energy);

      bool rcleaned = false;
      bool keep = true;

      //Apply Q tests
      for (unsigned int i = 0; i < qualityTests_.size(); ++i) {
        if (!qualityTests_.at(i)->test(rh, hgrh, rcleaned)) {
          keep = false;
        }
      }

      if (keep) {
        rh.setTime(time);
        out->push_back(rh);
      } else if (rcleaned)
        cleaned->push_back(rh);
    }
    edm::LogInfo("HGCalRecHitCreator") << "Skipped " << skipped_rechits << " out of " << rechits.size() << " rechits!"
                                       << std::endl;
    edm::LogInfo("HGCalRecHitCreator") << "Created " << out->size() << " PFRecHits!" << std::endl;
  }

Member Data Documentation

geometryInstance_

template<typename DET , PFLayer::Layer Layer, DetId::Detector det, unsigned subdet>

std::string PFHGCalRecHitCreator< DET, Layer, det, subdet >::geometryInstance_

protected

Definition at line 99 of file PFHGCalRecHitCreator.h.

geomToken_

template<typename DET , PFLayer::Layer Layer, DetId::Detector det, unsigned subdet>

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> PFHGCalRecHitCreator< DET, Layer, det, subdet >::geomToken_

private

Definition at line 103 of file PFHGCalRecHitCreator.h.

Referenced by PFHGCalRecHitCreator< DET, Layer, det, subdet >::importRecHits().

recHitToken_

template<typename DET , PFLayer::Layer Layer, DetId::Detector det, unsigned subdet>

edm::EDGetTokenT<HGCRecHitCollection> PFHGCalRecHitCreator< DET, Layer, det, subdet >::recHitToken_

protected

Definition at line 98 of file PFHGCalRecHitCreator.h.

Referenced by PFHGCalRecHitCreator< DET, Layer, det, subdet >::importRecHits().

recHitTools_

template<typename DET , PFLayer::Layer Layer, DetId::Detector det, unsigned subdet>

hgcal::RecHitTools PFHGCalRecHitCreator< DET, Layer, det, subdet >::recHitTools_

private

Definition at line 102 of file PFHGCalRecHitCreator.h.

Referenced by PFHGCalRecHitCreator< DET, Layer, det, subdet >::importRecHits().