#include <PFHGCalRecHitCreator.h>

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

Public Member Functions
void	importRecHits (std::auto_ptr< reco::PFRecHitCollection > &out, std::auto_ptr< reco::PFRecHitCollection > &cleaned, const edm::Event &iEvent, const edm::EventSetup &iSetup)

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

Public Member Functions inherited from PFRecHitCreatorBase
	PFRecHitCreatorBase ()

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

Protected Attributes
std::string	geometryInstance_

edm::EDGetTokenT < HGCRecHitCollection >	recHitToken_

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

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, ForwardSubdetector subdet>
class PFHGCalRecHitCreator< DET, Layer, subdet >

Definition at line 23 of file PFHGCalRecHitCreator.h.

Constructor & Destructor Documentation

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

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

inline

Definition at line 26 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, ForwardSubdetector subdet>

void PFHGCalRecHitCreator< DET, Layer, subdet >::importRecHits	(	std::auto_ptr< reco::PFRecHitCollection > &	out,
		std::auto_ptr< reco::PFRecHitCollection > &	cleaned,
		const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

inlinevirtual

Implements PFRecHitCreatorBase.

Definition at line 33 of file PFHGCalRecHitCreator.h.

References assert(), CaloRecHit::detid(), CaloRecHit::energy(), relval_parameters_module::energy, Exception, PFHGCalRecHitCreator< DET, Layer, subdet >::geometryInstance_, edm::EventSetup::get(), edm::Event::getByToken(), HGCalGeometry::getCorners(), HGCalGeometry::getGeometry(), HGCalGeometry::getPosition(), i, keep, LogDebug, eostools::move(), position, PFRecHitCreatorBase::qualityTests_, HI_PhotonSkim_cff::rechits, PFHGCalRecHitCreator< DET, Layer, subdet >::recHitToken_, reco::PFRecHit::setOriginalRecHit(), edm::SortedCollection< T, SORT >::size(), CaloRecHit::time(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                                                                  {
 
       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<HGCalGeometry> geoHandle;
       iSetup.get<IdealGeometryRecord>().get(geometryInstance_,geoHandle);
       const HGCalGeometry& hgcGeo = *geoHandle;
       
       unsigned skipped_rechits = 0;
       for (unsigned int i=0;i<rechits.size();++i) {
     const HGCRecHit& hgrh = rechits[i];
     const DET detid(hgrh.detid());
     
     if( subdet != detid.subdet() ) {
       throw cms::Exception("IncorrectHGCSubdetector")
         << "subdet expected: " << subdet 
         << " subdet gotten: " << detid.subdet() << std::endl;
     }
     
     double energy = hgrh.energy();
     double time = hgrh.time();  
     
     const FlatTrd *thisCell = 
       static_cast<const FlatTrd*>(hgcGeo.getGeometry(detid));   
 
     // find rechit geometry
     if(!thisCell) {
       LogDebug("PFHGCalRecHitCreator")
         <<"warning detid "<<detid.rawId()
         <<" not found in geometry"<<std::endl;
       ++skipped_rechits;
       continue;
     }
   
     const GlobalPoint position( std::move( hgcGeo.getPosition( detid ) ) );
     //std::cout << "geometry cell position: " << position << std::endl;
 
     reco::PFRecHit rh( detid.rawId(),Layer,
                energy, 
                position.x(), position.y(), position.z(), 
                0, 0, 0 ); 
     
     rh.setOriginalRecHit(edm::Ref<HGCRecHitCollection>(recHitHandle,i));
 
     const HGCalGeometry::CornersVec corners( std::move( hgcGeo.getCorners( detid ) ) );
     assert( corners.size() == 8 );
 
     rh.setNECorner( corners[0].x(), corners[0].y(),  corners[0].z() );
     rh.setSECorner( corners[1].x(), corners[1].y(),  corners[1].z() );
     rh.setSWCorner( corners[2].x(), corners[2].y(),  corners[2].z() );
     rh.setNWCorner( corners[3].x(), corners[3].y(),  corners[3].z() );
     
     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;
     }