#include <PFRecHitCaloNavigatorWithTime.h>

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

Public Member Functions
void	associateNeighbours (reco::PFRecHit &hit, std::auto_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::auto_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::auto_ptr< reco::PFRecHitCollection > &hits, edm::RefProd< reco::PFRecHitCollection > &refProd, short eta, short phi, short depth)

Protected Attributes
std::unique_ptr < ECALRecHitResolutionProvider >	_timeResolutionCalc

double	constantTerm2_

double	noiseLevel2_

double	noiseTerm2_

double	sigmaCut2_

std::unique_ptr< const T >	topology_

Detailed Description

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

Definition at line 26 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 28 of file PFRecHitCaloNavigatorWithTime.h.

                                                                 {
     noiseLevel2_ = pow(iConfig.getParameter<double>("noiseLevel"), 2);
     noiseTerm2_ = pow(iConfig.getParameter<double>("noiseTerm"), 2);
     constantTerm2_ = pow(iConfig.getParameter<double>("constantTerm"), 2);
     sigmaCut2_ = pow(iConfig.getParameter<double>("sigmaCut"), 2);
 
     _timeResolutionCalc.reset(NULL);
     if( iConfig.exists("timeResolutionCalc") ) {
       const edm::ParameterSet& timeResConf = 
         iConfig.getParameterSet("timeResolutionCalc");
         _timeResolutionCalc.reset(new ECALRecHitResolutionProvider(timeResConf));
     }
   }

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

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

inlinevirtual

Definition at line 42 of file PFRecHitCaloNavigatorWithTime.h.

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

PFRecHitCaloNavigatorWithTime::topology_

std::unique_ptr< const T > topology_

Definition: PFRecHitCaloNavigatorWithTime.h:149

Member Function Documentation

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

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

inlineprotected

Definition at line 121 of file PFRecHitCaloNavigatorWithTime.h.

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

                                                                                                                                                                          {
     double sigma2=10000.0;
     
     const reco::PFRecHit temp(id,PFLayer::NONE,0.0,math::XYZPoint(0,0,0),math::XYZVector(0,0,0),std::vector<math::XYZPoint>());
     auto found_hit = std::lower_bound(hits->begin(),hits->end(),
                       temp,
                       [](const reco::PFRecHit& a, 
                      const reco::PFRecHit& b){
                     return a.detId() < b.detId();
                       });
     if( found_hit != hits->end() && found_hit->detId() == id.rawId() ) {
       if (_timeResolutionCalc) {
         sigma2 = _timeResolutionCalc->timeResolution2(hit.energy()) + _timeResolutionCalc->timeResolution2(found_hit->energy());
       }
       else {
         const double hitEnergy = hit.energy();
         const double hitEnergy2 = hitEnergy*hitEnergy;
         const double fhEnergy = found_hit->energy();
         const double fhEnergy2 = fhEnergy*fhEnergy;
         sigma2 = noiseTerm2_*noiseLevel2_*(hitEnergy2+fhEnergy2)/(hitEnergy2*fhEnergy2) + 2*constantTerm2_;
       }
       const double deltaTime = hit.time()-found_hit->time();
       if(deltaTime*deltaTime/sigma2<sigmaCut2_) {
         hit.addNeighbour(eta,phi,0,reco::PFRecHitRef(refProd,std::distance(hits->begin(),found_hit)));
       }
     }
   }

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

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

inlinevirtual

Implements PFRecHitNavigatorBase.

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