PFRecHitCaloNavigatorWithTime< D, T, ownsTopo > Class Template Reference

#include <PFRecHitCaloNavigatorWithTime.h>

Inheritance diagram for PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >:

Public Member Functions
void	associateNeighbours (reco::PFRecHit &hit, std::unique_ptr< reco::PFRecHitCollection > &hits, edm::RefProd< reco::PFRecHitCollection > &refProd)
	PFRecHitCaloNavigatorWithTime (const edm::ParameterSet &iConfig)
virtual	~PFRecHitCaloNavigatorWithTime ()
Public Member Functions inherited from PFRecHitNavigatorBase
virtual void	beginEvent (const edm::EventSetup &)=0
	PFRecHitNavigatorBase ()
	PFRecHitNavigatorBase (const edm::ParameterSet &iConfig)
virtual	~PFRecHitNavigatorBase ()

Protected Member Functions
void	associateNeighbour (const DetId &id, reco::PFRecHit &hit, std::unique_ptr< reco::PFRecHitCollection > &hits, edm::RefProd< reco::PFRecHitCollection > &refProd, short eta, short phi)
Protected Member Functions inherited from PFRecHitNavigatorBase
void	associateNeighbour (const DetId &id, reco::PFRecHit &hit, std::unique_ptr< reco::PFRecHitCollection > &hits, edm::RefProd< reco::PFRecHitCollection > &refProd, short eta, short phi, short depth)

Protected Attributes
std::unique_ptr< CaloRecHitResolutionProvider >	_timeResolutionCalc
double	sigmaCut2_
std::unique_ptr< const T >	topology_

Additional Inherited Members
Public Types inherited from PFRecHitNavigatorBase
typedef std::unordered_map< unsigned, unsigned >	DetIdToHitIdx

Detailed Description

template<typename D, typename T, bool ownsTopo = true>
class PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >

Definition at line 27 of file PFRecHitCaloNavigatorWithTime.h.

Constructor & Destructor Documentation

template<typename D, typename T, bool ownsTopo = true>

PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >::PFRecHitCaloNavigatorWithTime ( const edm::ParameterSet &  iConfig )

inline

Definition at line 29 of file PFRecHitCaloNavigatorWithTime.h.

                                                                {
    sigmaCut2_ = pow(iConfig.getParameter<double>("sigmaCut"), 2);
    const edm::ParameterSet& timeResConf = iConfig.getParameterSet("timeResolutionCalc");
    _timeResolutionCalc.reset(new CaloRecHitResolutionProvider(timeResConf));
  }

template<typename D, typename T, bool ownsTopo = true>

virtual PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >::~PFRecHitCaloNavigatorWithTime ( )

inlinevirtual

Definition at line 36 of file PFRecHitCaloNavigatorWithTime.h.

{ if(!ownsTopo) { topology_.release(); } }

Member Function Documentation

template<typename D, typename T, bool ownsTopo = true>

void PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >::associateNeighbour ( const DetId &  id, reco::PFRecHit &  hit, std::unique_ptr< reco::PFRecHitCollection > &  hits, edm::RefProd< reco::PFRecHitCollection > &  refProd, short  eta, short  phi  )

inlineprotected

Definition at line 119 of file PFRecHitCaloNavigatorWithTime.h.

Referenced by PFRecHitCaloNavigatorWithTime< EEDetId, EcalEndcapTopology >::associateNeighbours().

                                                                                                                                                                           {
    double sigma2=10000.0;
    

    auto found_hit = std::lower_bound(hits->begin(),hits->end(),
                      id,
                      [](const reco::PFRecHit& a, 
                      DetId b){
                    if (DetId(a.detId()).det() == DetId::Hcal || b.det() == DetId::Hcal) return (HcalDetId)(a.detId()) < (HcalDetId)(b);

                    else return a.detId() < b.rawId();
                      });


    if (found_hit != hits->end() && ((found_hit->detId() == id.rawId()) ||
                     ((id.det()==DetId::Hcal) &&
                      ((HcalDetId)(found_hit->detId()) == (HcalDetId)(id))))) {
      sigma2 = _timeResolutionCalc->timeResolution2(hit.energy()) + _timeResolutionCalc->timeResolution2(found_hit->energy());
      const double deltaTime = hit.time()-found_hit->time();
      if(deltaTime*deltaTime<sigmaCut2_*sigma2) {
    hit.addNeighbour(eta,phi,0,found_hit-hits->begin());
      }
    }

  }

template<typename D, typename T, bool ownsTopo = true>

void PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >::associateNeighbours ( reco::PFRecHit &  hit, std::unique_ptr< reco::PFRecHitCollection > &  hits, edm::RefProd< reco::PFRecHitCollection > &  refProd  )

inlinevirtual

Implements PFRecHitNavigatorBase.

Definition at line 39 of file PFRecHitCaloNavigatorWithTime.h.

                                                                                                                                      {
      DetId detid( hit.detId() );
      
      CaloNavigator<D> navigator(detid, topology_.get());
      
      DetId N(0);
      DetId E(0);
      DetId S(0);
      DetId W(0);
      DetId NW(0);
      DetId NE(0);
      DetId SW(0);
      DetId SE(0);


      N=navigator.north();  
      associateNeighbour(N,hit,hits,refProd,0,1);


      if (N !=DetId(0)) {
    NE=navigator.east();
      }
      else 
    {
      navigator.home();
      E=navigator.east();
      NE=navigator.north();
    }
      associateNeighbour(NE,hit,hits,refProd,1,1);
      navigator.home();

      S = navigator.south();
      associateNeighbour(S,hit,hits,refProd,0,-1);
      
      if (S !=DetId(0)) {
    SW = navigator.west();
      } else {
    navigator.home();
    W=navigator.west();
    SW=navigator.south();
      }
      associateNeighbour(SW,hit,hits,refProd,-1,-1);
      navigator.home();

      E = navigator.east();
      associateNeighbour(E,hit,hits,refProd,1,0);
      
      if (E !=DetId(0)) {
    SE = navigator.south();
      } else {
    navigator.home();
    S=navigator.south();
    SE=navigator.east();
      }
      associateNeighbour(SE,hit,hits,refProd,1,-1);
      navigator.home();


      W = navigator.west();
      associateNeighbour(W,hit,hits,refProd,-1,0);

      if (W !=DetId(0)) {
    NW = navigator.north();
      } else {
    navigator.home();
    N=navigator.north();
    NW=navigator.west();
      }
      associateNeighbour(NW,hit,hits,refProd,-1,1);
  }

Member Data Documentation

template<typename D, typename T, bool ownsTopo = true>

std::unique_ptr<CaloRecHitResolutionProvider> PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >::_timeResolutionCalc

protected

Definition at line 116 of file PFRecHitCaloNavigatorWithTime.h.

Referenced by PFRecHitCaloNavigatorWithTime< EEDetId, EcalEndcapTopology >::PFRecHitCaloNavigatorWithTime().

template<typename D, typename T, bool ownsTopo = true>

double PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >::sigmaCut2_

protected

Definition at line 114 of file PFRecHitCaloNavigatorWithTime.h.

Referenced by PFRecHitCaloNavigatorWithTime< EEDetId, EcalEndcapTopology >::PFRecHitCaloNavigatorWithTime().

template<typename D, typename T, bool ownsTopo = true>

std::unique_ptr<const T> PFRecHitCaloNavigatorWithTime< D, T, ownsTopo >::topology_

protected

Definition at line 115 of file PFRecHitCaloNavigatorWithTime.h.

Referenced by PFRecHitCaloNavigatorWithTime< EEDetId, EcalEndcapTopology >::associateNeighbours(), and PFRecHitCaloNavigatorWithTime< EEDetId, EcalEndcapTopology >::~PFRecHitCaloNavigatorWithTime().