#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) override

	PFRecHitCaloNavigatorWithTime (const edm::ParameterSet &iConfig, edm::ConsumesCollector &cc)

	~PFRecHitCaloNavigatorWithTime () override

Public Member Functions inherited from PFRecHitNavigatorBase
virtual void	init (const edm::EventSetup &)=0

	PFRecHitNavigatorBase ()=default

	PFRecHitNavigatorBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &cc)

virtual	~PFRecHitNavigatorBase ()=default

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 29 of file PFRecHitCaloNavigatorWithTime.h.

Constructor & Destructor Documentation

◆ PFRecHitCaloNavigatorWithTime()

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

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

inline

Definition at line 31 of file PFRecHitCaloNavigatorWithTime.h.

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

◆ ~PFRecHitCaloNavigatorWithTime()

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

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

inlineoverride

Definition at line 37 of file PFRecHitCaloNavigatorWithTime.h.

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

Member Function Documentation

◆ associateNeighbour()

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 116 of file PFRecHitCaloNavigatorWithTime.h.

                                      {
     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());
       }
     }
   }

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

◆ associateNeighbours()

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
	)

inlineoverridevirtual

Implements PFRecHitNavigatorBase.

Definition at line 43 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);
   }