#include <PFHBHERecHitCreator.h>

Inheritance diagram for PFHBHERecHitCreator:

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

	PFHBHERecHitCreator (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
edm::EDGetTokenT< edm::SortedCollection< HBHERecHit > >	recHitToken_

Protected Attributes inherited from PFRecHitCreatorBase
std::vector< std::unique_ptr< PFRecHitQTestBase > >	qualityTests_

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomToken_

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

Detailed Description

Definition at line 19 of file PFHBHERecHitCreator.h.

Constructor & Destructor Documentation

◆ PFHBHERecHitCreator()

PFHBHERecHitCreator::PFHBHERecHitCreator	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	cc
	)

inline

Definition at line 21 of file PFHBHERecHitCreator.h.

       : PFRecHitCreatorBase(iConfig, cc),
         recHitToken_(cc.consumes<edm::SortedCollection<HBHERecHit> >(iConfig.getParameter<edm::InputTag>("src"))),
         geomToken_(cc.esConsumes()) {}

Member Function Documentation

◆ importRecHits()

void PFHBHERecHitCreator::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 26 of file PFHBHERecHitCreator.h.

References PFRecHitCreatorBase::beginEvent(), TauDecayModes::dec, hcalRecHitTable_cff::depth, HcalDetId::depth(), hcalRecHitTable_cff::energy, geomToken_, CaloSubdetectorGeometry::getGeometry(), edm::EventSetup::getHandle(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, PFLayer::HCAL_BARREL1, PFLayer::HCAL_ENDCAP, HcalBarrel, HcalEndcap, iEvent, runTheMatrix::keep, nano_mu_digi_cff::layer, eostools::move(), MillePedeFileConverter_cfg::out, PFRecHitCreatorBase::qualityTests_, DetId::rawId(), recHitToken_, reco::PFRecHit::setTime(), DetId::subdetId(), and hcalRecHitTable_cff::time.

                                                            {
     beginEvent(iEvent, iSetup);
 
     edm::Handle<edm::SortedCollection<HBHERecHit> > recHitHandle;
 
     edm::ESHandle<CaloGeometry> geoHandle = iSetup.getHandle(geomToken_);
 
     // get the ecal geometry
     const CaloSubdetectorGeometry* hcalBarrelGeo = geoHandle->getSubdetectorGeometry(DetId::Hcal, HcalBarrel);
     const CaloSubdetectorGeometry* hcalEndcapGeo = geoHandle->getSubdetectorGeometry(DetId::Hcal, HcalEndcap);
 
     iEvent.getByToken(recHitToken_, recHitHandle);
     for (const auto& erh : *recHitHandle) {
       const HcalDetId detid = erh.idFront();
       HcalSubdetector esd = (HcalSubdetector)detid.subdetId();
 
       auto energy = erh.energy();
       auto time = erh.time();
       auto depth = detid.depth();
 
       std::shared_ptr<const CaloCellGeometry> thisCell = nullptr;
       PFLayer::Layer layer = PFLayer::HCAL_BARREL1;
       switch (esd) {
         case HcalBarrel:
           thisCell = hcalBarrelGeo->getGeometry(detid);
           layer = PFLayer::HCAL_BARREL1;
           break;
 
         case HcalEndcap:
           thisCell = hcalEndcapGeo->getGeometry(detid);
           layer = PFLayer::HCAL_ENDCAP;
           break;
         default:
           break;
       }
 
       // find rechit geometry
       if (!thisCell) {
         edm::LogError("PFHBHERecHitCreator") << "warning detid " << std::hex << detid.rawId() << std::dec << " "
                                              << detid << " not found in geometry" << std::endl;
         continue;
       }
 
       reco::PFRecHit rh(thisCell, detid.rawId(), layer, energy);
       rh.setTime(time);  //Mike: This we will use later
       rh.setDepth(depth);
 
       bool rcleaned = false;
       bool keep = true;
 
       //Apply Q tests
       for (const auto& qtest : qualityTests_) {
         if (!qtest->test(rh, erh, rcleaned)) {
           keep = false;
         }
       }
 
       if (keep) {
         out->push_back(std::move(rh));
       } else if (rcleaned)
         cleaned->push_back(std::move(rh));
     }
   }