#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	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	iC
	)

inline

Definition at line 33 of file PFEcalEndcapRecHitCreator.h.

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

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

Member Function Documentation

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 43 of file PFEcalEndcapRecHitCreator.h.

References PFRecHitCreatorBase::beginEvent(), DetId::Ecal, PFLayer::ECAL_ENDCAP, EcalEndcap, randomXiThetaGunProducer_cfi::energy, Exception, edm::EventSetup::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), isHighInterest(), keep, eostools::move(), PFRecHitCreatorBase::qualityTests_, DetId::rawId(), recHitToken_, srFlagHandle_, srFlagToken_, and ntuplemaker::time.

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

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

inlineoverridevirtual

Reimplemented from PFRecHitCreatorBase.

Definition at line 106 of file PFEcalEndcapRecHitCreator.h.

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

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

bool PFEcalEndcapRecHitCreator::isHighInterest ( const EEDetId & detid )

inlineprotected

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

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

inlineprotected

Definition at line 125 of file PFEcalEndcapRecHitCreator.h.

References EcalElectronicsId::dccId(), elecMap_, EcalElectronicsMapping::getEcalScDetId(), EcalElectronicsMapping::getElectronicsId(), and EcalElectronicsId::towerId().

Referenced by isHighInterest().

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