#include <PFEcalEndcapRecHitCreator.h>

Inheritance diagram for PFEcalEndcapRecHitCreator:

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

	PFEcalEndcapRecHitCreator (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 EEDetId &detid)

EcalScDetId	readOutUnitOf (const EEDetId &detid) const

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

Protected Attributes
const EcalElectronicsMapping *	elecMap_

const EcalTrigTowerConstituentsMap *	eTTmap_

edm::EDGetTokenT< EcalRecHitCollection >	recHitToken_

edm::Handle< EESrFlagCollection >	srFlagHandle_

edm::EDGetTokenT< EESrFlagCollection >	srFlagToken_

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

Detailed Description

Definition at line 30 of file PFEcalEndcapRecHitCreator.h.

Constructor & Destructor Documentation

◆ PFEcalEndcapRecHitCreator()

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

inline

Definition at line 32 of file PFEcalEndcapRecHitCreator.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<EESrFlagCollection>(srF);
     elecMap_ = nullptr;
   }

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

Member Function Documentation

◆ importRecHits()

void PFEcalEndcapRecHitCreator::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 PFEcalEndcapRecHitCreator.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, EcalEndcap);
  
     const EcalEndcapGeometry* ecalGeo = dynamic_cast<const EcalEndcapGeometry*>(gTmp);
  
     iEvent.getByToken(recHitToken_, recHitHandle);
     for (const auto& erh : *recHitHandle) {
       const DetId& detid = erh.detid();
       auto energy = erh.energy();
       auto time = erh.time();
  
       bool hi = (useSrF ? isHighInterest(detid) : true);
  
       std::shared_ptr<const CaloCellGeometry> thisCell = ecalGeo->getGeometry(detid);
  
       // find rechit geometry
       if (!thisCell) {
         throw cms::Exception("PFEcalEndcapRecHitCreator") << "detid " << detid.rawId() << "not found in geometry";
       }
  
       out->emplace_back(thisCell, detid.rawId(), PFLayer::ECAL_ENDCAP, energy);
  
       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_ENDCAP, EcalEndcap, HCALHighEnergyHPDFilter_cfi::energy, Exception, 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 PFEcalEndcapRecHitCreator::init ( const edm::EventSetup & es )

inlineoverridevirtual

Reimplemented from PFRecHitCreatorBase.

Definition at line 103 of file PFEcalEndcapRecHitCreator.h.

                                               {
     edm::ESHandle<EcalElectronicsMapping> ecalmapping;
     es.get<EcalMappingRcd>().get(ecalmapping);
     elecMap_ = ecalmapping.product();
   }

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

◆ isHighInterest()

bool PFEcalEndcapRecHitCreator::isHighInterest ( const EEDetId & detid )

inlineprotected

Definition at line 110 of file PFEcalEndcapRecHitCreator.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;
   }

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

Referenced by importRecHits().

◆ readOutUnitOf()

EcalScDetId PFEcalEndcapRecHitCreator::readOutUnitOf ( const EEDetId & detid ) const

inlineprotected

Definition at line 120 of file PFEcalEndcapRecHitCreator.h.

                                                         {
     const EcalElectronicsId& EcalElecId = elecMap_->getElectronicsId(detid);
     int iDCC = EcalElecId.dccId();
     int iDccChan = EcalElecId.towerId();
     const bool ignoreSingle = true;
     const std::vector<EcalScDetId> id = elecMap_->getEcalScDetId(iDCC, iDccChan, ignoreSingle);
     return !id.empty() ? id[0] : EcalScDetId();
   }