#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 &cc)

Public Member Functions inherited from PFRecHitCreatorBase
	PFRecHitCreatorBase ()

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

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_

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

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	towerToken_

Detailed Description

Definition at line 28 of file PFEcalBarrelRecHitCreator.h.

Constructor & Destructor Documentation

◆ PFEcalBarrelRecHitCreator()

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

inline

Definition at line 30 of file PFEcalBarrelRecHitCreator.h.

References gpuPixelDoublets::cc, edm::ParameterSet::getParameter(), and srFlagToken_.

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

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 41 of file PFEcalBarrelRecHitCreator.h.

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

                                                            {
     beginEvent(iEvent, iSetup);
 
     edm::Handle<EcalRecHitCollection> recHitHandle;
 
     edm::ESHandle<CaloGeometry> geoHandle = iSetup.getHandle(geomToken_);
 
     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();
       }
     }
   }

◆ init()

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

inlineoverridevirtual

Reimplemented from PFRecHitCreatorBase.

Definition at line 102 of file PFEcalBarrelRecHitCreator.h.

References edm::EventSetup::getData(), towerToken_, and triggerTowerMap_.

102 { triggerTowerMap_ = &es.getData(towerToken_); }

PFEcalBarrelRecHitCreator::towerToken_

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord > towerToken_

Definition: PFEcalBarrelRecHitCreator.h:127

edm::EventSetup::getData

T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)

Definition: EventSetup.h:119

PFEcalBarrelRecHitCreator::triggerTowerMap_

const EcalTrigTowerConstituentsMap * triggerTowerMap_

Definition: PFEcalBarrelRecHitCreator.h:121

◆ isHighInterest()

bool PFEcalBarrelRecHitCreator::isHighInterest ( const EBDetId & detid )

inlineprotected

Definition at line 105 of file PFEcalBarrelRecHitCreator.h.

References readOutUnitOf(), mps_fire::result, EcalSrFlag::SRF_FULL, and srFlagHandle_.

Referenced by importRecHits().

                                             {
     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;
   }