#include <PFHGCalRecHitCreator.h>

Inheritance diagram for PFHGCalRecHitCreator< DET, Layer, subdet >:

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

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

Public Member Functions inherited from PFRecHitCreatorBase
virtual void	init (const edm::EventSetup &es)

	PFRecHitCreatorBase ()

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

virtual	~PFRecHitCreatorBase ()=default

Protected Attributes
std::string	geometryInstance_

edm::EDGetTokenT< HGCRecHitCollection >	recHitToken_

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

Private Attributes
std::vector< const CaloCellGeometryHGCALAdapter * >	caloCells_

hgcal::RecHitTools	recHitTools_

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

Detailed Description

template<typename DET, PFLayer::Layer Layer, unsigned subdet>
class PFHGCalRecHitCreator< DET, Layer, subdet >

Definition at line 27 of file PFHGCalRecHitCreator.h.

Constructor & Destructor Documentation

template<typename DET , PFLayer::Layer Layer, unsigned subdet>

PFHGCalRecHitCreator< DET, Layer, subdet >::PFHGCalRecHitCreator	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	iC
	)

inline

Definition at line 30 of file PFHGCalRecHitCreator.h.

References edm::ConsumesCollector::consumes(), PFHGCalRecHitCreator< DET, Layer, subdet >::geometryInstance_, edm::ParameterSet::getParameter(), PFHGCalRecHitCreator< DET, Layer, subdet >::recHitToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                :
     PFRecHitCreatorBase(iConfig,iC)
     {
       recHitToken_ = iC.consumes<HGCRecHitCollection>(iConfig.getParameter<edm::InputTag>("src"));
       geometryInstance_ = iConfig.getParameter<std::string>("geometryInstance");
     }

Member Function Documentation

template<typename DET , PFLayer::Layer Layer, unsigned subdet>

void PFHGCalRecHitCreator< DET, Layer, subdet >::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 37 of file PFHGCalRecHitCreator.h.

References EnergyCorrector::c, PFHGCalRecHitCreator< DET, Layer, subdet >::caloCells_, Exception, relativeConstraints::geom, edm::EventSetup::get(), edm::Event::getByToken(), hgcal::RecHitTools::getEventSetup(), ecaldqm::getGeometry(), hgcal::RecHitTools::getPosition(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, mps_fire::i, keep, LogDebug, edm::ESHandle< T >::product(), PFRecHitCreatorBase::qualityTests_, TrackInfoProducer_cfi::rechits, PFHGCalRecHitCreator< DET, Layer, subdet >::recHitToken_, PFHGCalRecHitCreator< DET, Layer, subdet >::recHitTools_, edm::SortedCollection< T, SORT >::size(), and ntuplemaker::time.

                                                                                                                                                                               {
 
       // Setup RecHitTools to properly compute the position of the HGCAL Cells vie their DetIds
       recHitTools_.getEventSetup(iSetup);
 
       for (unsigned int i=0;i<qualityTests_.size();++i) {
     qualityTests_.at(i)->beginEvent(iEvent,iSetup);
       }
 
       edm::Handle<HGCRecHitCollection> recHitHandle;
       iEvent.getByToken(recHitToken_,recHitHandle);
       const HGCRecHitCollection& rechits = *recHitHandle;
 
       edm::ESHandle<CaloGeometry> geoHandle;
       iSetup.get<CaloGeometryRecord>().get(geoHandle);
       const CaloGeometry* geom = geoHandle.product();
 
       // Get rid of the content of the previous event
       for (auto c : caloCells_)
         delete c;
 
       caloCells_.clear();
       caloCells_.reserve(rechits.size());
 
       unsigned skipped_rechits = 0;
       for (const auto & hgrh : rechits) {
         const DET detid(hgrh.detid());
     
     if( subdet != detid.subdetId() ) {
       throw cms::Exception("IncorrectHGCSubdetector")
         << "subdet expected: " << subdet 
         << " subdet gotten: " << detid.subdetId() << std::endl;
     }
     
     double energy = hgrh.energy();
     double time = hgrh.time();  
     
     const CaloCellGeometry *thisCell;
     if( detid.det() == DetId::Hcal ) {
       thisCell = geom->getSubdetectorGeometry(detid.det(),detid.subdetId())->getGeometry(detid);
     } else {
       const auto* hg = static_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(detid.det(),detid.subdetId()));
       caloCells_.push_back(new CaloCellGeometryHGCALAdapter(static_cast<const FlatTrd*>(hg->getGeometry(detid)),
                            recHitTools_.getPosition(detid)));
       thisCell = caloCells_.back();
     }
 
     // find rechit geometry
     if(!thisCell) {
       LogDebug("PFHGCalRecHitCreator")
         <<"warning detid "<<detid.rawId()
         <<" not found in geometry"<<std::endl;
       ++skipped_rechits;
       continue;
     }
   
 
     reco::PFRecHit rh(thisCell, detid.rawId(),Layer,
                energy); 
     //  rh.setOriginalRecHit(edm::Ref<HGCRecHitCollection>(recHitHandle,i));
 
     
     bool rcleaned = false;
     bool keep=true;
 
     //Apply Q tests
     for (unsigned int i=0;i<qualityTests_.size();++i) {
       if (!qualityTests_.at(i)->test(rh,hgrh,rcleaned)) {
         keep = false;       
       }
     }
       
     if(keep) {
       rh.setTime(time);
       out->push_back(rh);
     }
     else if (rcleaned) 
       cleaned->push_back(rh);
       }
       edm::LogInfo("HGCalRecHitCreator") 
     <<  "Skipped " << skipped_rechits 
     << " out of " << rechits.size() << " rechits!" << std::endl;
       edm::LogInfo("HGCalRecHitCreator")
     << "Created " << out->size() << " PFRecHits!" << std::endl;
     }