PFHcalRecHitCreator< Digi, Geometry, Layer, Detector > Class Template Reference

#include <PFHcalRecHitCreator.h>

Inheritance diagram for PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >:

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

Protected Attributes
int	hoDepth_
edm::EDGetTokenT < edm::SortedCollection< Digi > >	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 Digi, typename Geometry, PFLayer::Layer Layer, int Detector>
class PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >

Definition at line 21 of file PFHcalRecHitCreator.h.

Constructor & Destructor Documentation

template<typename Digi , typename Geometry , PFLayer::Layer Layer, int Detector>

PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::PFHcalRecHitCreator ( const edm::ParameterSet &  iConfig, edm::ConsumesCollector &  iC  )

inline

Definition at line 24 of file PFHcalRecHitCreator.h.

References edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), and PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::recHitToken_.

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

Member Function Documentation

template<typename Digi , typename Geometry , PFLayer::Layer Layer, int Detector>

void PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::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 30 of file PFHcalRecHitCreator.h.

References assert(), PFRecHitCreatorBase::beginEvent(), HcalDetId::depth(), HLT_FULL_cff::depth, align::Detector, relval_parameters_module::energy, edm::EventSetup::get(), edm::Event::getByToken(), CaloCellGeometry::getCorners(), CaloCellGeometry::getPosition(), DetId::Hcal, HcalOther, keep, point, position, PFRecHitCreatorBase::qualityTests_, DetId::rawId(), PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::recHitToken_, reco::PFRecHit::setTime(), EZArrayFL< T >::size(), DetId::subdetId(), x, y, and z.

                                                                                                                                                                 {


      beginEvent(iEvent,iSetup);

      edm::Handle<edm::SortedCollection<Digi> > recHitHandle;

      edm::ESHandle<CaloGeometry> geoHandle;
      iSetup.get<CaloGeometryRecord>().get(geoHandle);
  
      // get the ecal geometry
      const CaloSubdetectorGeometry *gTmp = 
    geoHandle->getSubdetectorGeometry(DetId::Hcal, Detector);

      const Geometry *hcalGeo =dynamic_cast< const Geometry* > (gTmp);

      iEvent.getByToken(recHitToken_,recHitHandle);
      for( const auto& erh : *recHitHandle ) {      
    const HcalDetId& detid = (HcalDetId)erh.detid();
    HcalSubdetector esd=(HcalSubdetector)detid.subdetId();
    
    //since hbhe are together kill other detector
    if (esd !=Detector && Detector != HcalOther  ) 
      continue;


    double energy = erh.energy();
    double time = erh.time();
    int depth =detid.depth();
      
    math::XYZVector position;
    math::XYZVector axis;
    
    const CaloCellGeometry *thisCell;
    thisCell= hcalGeo->getGeometry(detid);
  
    // find rechit geometry
    if(!thisCell) {
      edm::LogError("PFHcalRecHitCreator")
        <<"warning detid "<<detid.rawId()
        <<" not found in geometry"<<std::endl;
      continue;
    }

    auto const point  = thisCell->getPosition();
    position.SetCoordinates ( point.x(),
                  point.y(),
                  point.z() );
  



    reco::PFRecHit rh( detid.rawId(),Layer,
               energy, 
               position.x(), position.y(), position.z(), 
               0,0,0);
    rh.setTime(time); //Mike: This we will use later
    rh.setDepth(depth);

    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

template<typename Digi , typename Geometry , PFLayer::Layer Layer, int Detector>

int PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::hoDepth_

protected

Definition at line 121 of file PFHcalRecHitCreator.h.

template<typename Digi , typename Geometry , PFLayer::Layer Layer, int Detector>

edm::EDGetTokenT<edm::SortedCollection<Digi> > PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::recHitToken_

protected

Definition at line 120 of file PFHcalRecHitCreator.h.

Referenced by PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::importRecHits(), and PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::PFHcalRecHitCreator().