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 &iC)
Public Member Functions inherited from PFRecHitCreatorBase
virtual void	init (const edm::EventSetup &es)
	PFRecHitCreatorBase ()
	PFRecHitCreatorBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
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
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

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

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

inline

Definition at line 28 of file PFHGCalRecHitCreator.h.

References edm::ConsumesCollector::consumes(), PFHGCalRecHitCreator< DET, Layer, det, subdet >::geometryInstance_, edm::ParameterSet::getParameter(), PFHGCalRecHitCreator< DET, Layer, det, subdet >::recHitToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

      : PFRecHitCreatorBase(iConfig, iC) {
    recHitToken_ = iC.consumes<HGCRecHitCollection>(iConfig.getParameter<edm::InputTag>("src"));
    geometryInstance_ = iConfig.getParameter<std::string>("geometryInstance");
  }

Member Function Documentation

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, edm::EventSetup::get(), edm::Event::getByToken(), hgcal::RecHitTools::getEventSetup(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), mps_fire::i, keep, LogDebug, or, edm::ESHandle< T >::product(), PFRecHitCreatorBase::qualityTests_, TrackInfoProducer_cfi::rechits, PFHGCalRecHitCreator< DET, Layer, det, subdet >::recHitToken_, PFHGCalRecHitCreator< DET, Layer, det, subdet >::recHitTools_, and ntuplemaker::time.

                                                           {
    // Setup RecHitTools to properly compute the position of the HGCAL Cells vie their DetIds
    recHitTools_.getEventSetup(iSetup);

    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;

    edm::ESHandle<CaloGeometry> geoHandle;
    iSetup.get<CaloGeometryRecord>().get(geoHandle);
    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

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

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

protected

Definition at line 100 of file PFHGCalRecHitCreator.h.

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

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

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

protected

Definition at line 99 of file PFHGCalRecHitCreator.h.

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

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

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

private

Definition at line 103 of file PFHGCalRecHitCreator.h.

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