PFPSRecHitCreator Class Reference

#include <PFPSRecHitCreator.h>

Inheritance diagram for PFPSRecHitCreator:

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)
	PFPSRecHitCreator (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
Public Member Functions inherited from PFRecHitCreatorBase
	PFRecHitCreatorBase ()
	PFRecHitCreatorBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)

Protected Attributes
edm::EDGetTokenT < EcalRecHitCollection >	recHitToken_
Protected Attributes inherited from PFRecHitCreatorBase
std::vector< std::unique_ptr < PFRecHitQTestBase > >	qualityTests_

Detailed Description

Definition at line 27 of file PFPSRecHitCreator.h.

Constructor & Destructor Documentation

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

inline

Definition at line 30 of file PFPSRecHitCreator.h.

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

                                                                            :
    PFRecHitCreatorBase(iConfig,iC)
    {
      recHitToken_ = iC.consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("src"));
    }

Member Function Documentation

void PFPSRecHitCreator::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 36 of file PFPSRecHitCreator.h.

References cond::rpcobgas::detid, DetId::Ecal, EcalPreshower, relval_parameters_module::energy, edm::hlt::Exception, edm::EventSetup::get(), edm::Event::getByToken(), CaloCellGeometry::getCorners(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), i, keep, PFLayer::NONE, position, PFLayer::PS1, PFLayer::PS2, PFRecHitCreatorBase::qualityTests_, recHitToken_, reco::PFRecHit::setNECorner(), reco::PFRecHit::setNWCorner(), reco::PFRecHit::setSECorner(), reco::PFRecHit::setSWCorner(), EZArrayFL< T >::size(), x, PV3DBase< T, PVType, FrameType >::x(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                                                                 {

      for (unsigned int i=0;i<qualityTests_.size();++i) {
    qualityTests_.at(i)->beginEvent(iEvent,iSetup);
      }


      edm::Handle<EcalRecHitCollection> recHitHandle;
      edm::ESHandle<CaloGeometry> geoHandle;
      iSetup.get<CaloGeometryRecord>().get(geoHandle);
  
      // get the ecal geometry
      const CaloSubdetectorGeometry *psGeometry = 
    geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);

      iEvent.getByToken(recHitToken_,recHitHandle);
      for( const auto& erh : *recHitHandle ) {      
    ESDetId detid(erh.detid());
    double energy = erh.energy();


    math::XYZVector position;

    PFLayer::Layer layer = PFLayer::NONE;
    
    switch( detid.plane() ) {
    case 1:
      layer = PFLayer::PS1;
      break;
    case 2:
      layer = PFLayer::PS2;
      break;
    default:
      throw cms::Exception("PFRecHitBadInput")
        <<"incorrect preshower plane !! plane number "
        <<detid.plane()<<std::endl;
    }
 

    
    const CaloCellGeometry *thisCell;
    thisCell= psGeometry->getGeometry(detid);
  
    // find rechit geometry
    if(!thisCell) {
      edm::LogError("PFPSRecHitCreator")
        <<"warning detid "<<detid.rawId()
        <<" not found in geometry"<<std::endl;
      continue;
    }
  
    position.SetCoordinates ( thisCell->getPosition().x(),
                  thisCell->getPosition().y(),
                  thisCell->getPosition().z() );
  
    reco::PFRecHit rh( detid.rawId(),layer,
               energy, 
               position.x(), position.y(), position.z(), 
               0.0,0.0,0.0);


    
    const CaloCellGeometry::CornersVec& corners = thisCell->getCorners();
    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( const auto& qtest : qualityTests_ ) {
      if (!qtest->test(rh,erh,rcleaned)) {
        keep = false;       
      }
    }
    
    if(keep) {
      out->push_back(rh);
    }
    else if (rcleaned) 
      cleaned->push_back(rh);
      }
    }

Member Data Documentation

edm::EDGetTokenT<EcalRecHitCollection> PFPSRecHitCreator::recHitToken_

protected

Definition at line 128 of file PFPSRecHitCreator.h.

Referenced by importRecHits(), and PFPSRecHitCreator().