#include <PFEcalBarrelRecHitCreator.h>

Inheritance diagram for PFEcalBarrelRecHitCreator:

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

void	init (const edm::EventSetup &es) override

	PFEcalBarrelRecHitCreator (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)

Public Member Functions inherited from PFRecHitCreatorBase
	PFRecHitCreatorBase ()

	PFRecHitCreatorBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)

virtual	~PFRecHitCreatorBase ()=default

Protected Member Functions
bool	isHighInterest (const EBDetId &detid)

EcalTrigTowerDetId	readOutUnitOf (const EBDetId &detid) const

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

Protected Attributes
edm::EDGetTokenT< EcalRecHitCollection >	recHitToken_

edm::Handle< EBSrFlagCollection >	srFlagHandle_

edm::EDGetTokenT< EBSrFlagCollection >	srFlagToken_

const EcalTrigTowerConstituentsMap *	triggerTowerMap_

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

Detailed Description

Definition at line 28 of file PFEcalBarrelRecHitCreator.h.

Constructor & Destructor Documentation

◆ PFEcalBarrelRecHitCreator()

PFEcalBarrelRecHitCreator::PFEcalBarrelRecHitCreator	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	iC
	)

inline

Definition at line 30 of file PFEcalBarrelRecHitCreator.h.

       : PFRecHitCreatorBase(iConfig, iC) {
     recHitToken_ = iC.consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("src"));
     auto srF = iConfig.getParameter<edm::InputTag>("srFlags");
     if (not srF.label().empty())
       srFlagToken_ = iC.consumes<EBSrFlagCollection>(srF);
     triggerTowerMap_ = nullptr;
   }

References edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), recHitToken_, srFlagToken_, and triggerTowerMap_.

Member Function Documentation

◆ importRecHits()

void PFEcalBarrelRecHitCreator::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 39 of file PFEcalBarrelRecHitCreator.h.

                                                            {
     beginEvent(iEvent, iSetup);
  
     edm::Handle<EcalRecHitCollection> recHitHandle;
  
     edm::ESHandle<CaloGeometry> geoHandle;
     iSetup.get<CaloGeometryRecord>().get(geoHandle);
  
     bool useSrF = false;
     if (not srFlagToken_.isUninitialized()) {
       iEvent.getByToken(srFlagToken_, srFlagHandle_);
       useSrF = true;
     }
  
     // get the ecal geometry
     const CaloSubdetectorGeometry* gTmp = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
  
     const EcalBarrelGeometry* ecalGeo = dynamic_cast<const EcalBarrelGeometry*>(gTmp);
  
     iEvent.getByToken(recHitToken_, recHitHandle);
     for (const auto& erh : *recHitHandle) {
       const DetId& detid = erh.detid();
       auto energy = erh.energy();
       auto time = erh.time();
       auto flags = erh.flagsBits();
       bool hi = (useSrF ? isHighInterest(detid) : true);
  
       const auto thisCell = ecalGeo->getGeometry(detid);
  
       // find rechit geometry
       if (!thisCell) {
         throw cms::Exception("PFEcalBarrelRecHitCreator") << "detid " << detid.rawId() << "not found in geometry";
       }
  
       out->emplace_back(thisCell, detid.rawId(), PFLayer::ECAL_BARREL, energy, flags);
  
       auto& rh = out->back();
  
       bool rcleaned = false;
       bool keep = true;
  
       //Apply Q tests
       for (const auto& qtest : qualityTests_) {
         if (!qtest->test(rh, erh, rcleaned, hi)) {
           keep = false;
         }
       }
  
       if (keep) {
         rh.setTime(time);
         rh.setDepth(1);
       } else {
         if (rcleaned)
           cleaned->push_back(std::move(out->back()));
         out->pop_back();
       }
     }
   }

References PFRecHitCreatorBase::beginEvent(), DetId::Ecal, PFLayer::ECAL_BARREL, EcalBarrel, HCALHighEnergyHPDFilter_cfi::energy, Exception, HLT_2018_cff::flags, edm::EventSetup::get(), get, CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), iEvent, isHighInterest(), keep, eostools::move(), MillePedeFileConverter_cfg::out, PFRecHitCreatorBase::qualityTests_, DetId::rawId(), recHitToken_, srFlagHandle_, srFlagToken_, and ntuplemaker::time.

◆ init()

void PFEcalBarrelRecHitCreator::init ( const edm::EventSetup & es )

inlineoverridevirtual

Reimplemented from PFRecHitCreatorBase.

Definition at line 101 of file PFEcalBarrelRecHitCreator.h.

                                               {
     edm::ESHandle<EcalTrigTowerConstituentsMap> hTriggerTowerMap;
     es.get<IdealGeometryRecord>().get(hTriggerTowerMap);
     triggerTowerMap_ = hTriggerTowerMap.product();
   }

References edm::EventSetup::get(), get, edm::ESHandle< T >::product(), and triggerTowerMap_.

◆ isHighInterest()

bool PFEcalBarrelRecHitCreator::isHighInterest ( const EBDetId & detid )

inlineprotected

Definition at line 108 of file PFEcalBarrelRecHitCreator.h.

                                             {
     bool result = false;
     auto srf = srFlagHandle_->find(readOutUnitOf(detid));
     if (srf == srFlagHandle_->end())
       return false;
     else
       result = ((srf->value() & ~EcalSrFlag::SRF_FORCED_MASK) == EcalSrFlag::SRF_FULL);
     return result;
   }