#include <CaloDualConeSelector.h>

Public Member Functions
	CaloDualConeSelector (double dRmin, double dRmax, const CaloGeometry *geom)

	CaloDualConeSelector (double dRmin, double dRmax, const CaloGeometry *geom, DetId::Detector detector, int subdet=0)

void	selectCallback (const GlobalPoint &p, const edm::SortedCollection< T > &inputCollection, std::function< void(const T &)> callback)

void	selectCallback (double eta, double phi, const edm::SortedCollection< T > &inputCollection, std::function< void(const T &)> callback)

Private Attributes
double	deltaRmax_

double	deltaRmin_

DetId::Detector	detector_

const CaloGeometry *	geom_

int	subdet_

Detailed Description

template<class T>
class CaloDualConeSelector< T >

Author: J. Mans - Minnesota

Definition at line 14 of file CaloDualConeSelector.h.

Constructor & Destructor Documentation

◆ CaloDualConeSelector() [1/2]

template<class T>

CaloDualConeSelector< T >::CaloDualConeSelector	(	double	dRmin,
		double	dRmax,
		const CaloGeometry *	geom
	)

inline

Definition at line 16 of file CaloDualConeSelector.h.

17 : geom_(geom), deltaRmin_(dRmin), deltaRmax_(dRmax), detector_(DetId::Detector(0)), subdet_(0) {}

◆ CaloDualConeSelector() [2/2]

template<class T>

CaloDualConeSelector< T >::CaloDualConeSelector	(	double	dRmin,
		double	dRmax,
		const CaloGeometry *	geom,
		DetId::Detector	detector,
		int	subdet = `0`
	)

inline

Definition at line 19 of file CaloDualConeSelector.h.

20 : geom_(geom), deltaRmin_(dRmin), deltaRmax_(dRmax), detector_(detector), subdet_(subdet) {}

Member Function Documentation

◆ selectCallback() [1/2]

template<class T>

void CaloDualConeSelector< T >::selectCallback	(	const GlobalPoint &	p,
		const edm::SortedCollection< T > &	inputCollection,
		std::function< void(const T &)>	callback
	)

inline

Definition at line 30 of file CaloDualConeSelector.h.

                                                                    {
     // TODO: handle default setting of detector_ (loops over subdet)
     // TODO: heuristics of when it is better to loop over inputCollection instead (small # hits)
     for (int subdet = subdet_; subdet <= 7 && (subdet_ == 0 || subdet_ == subdet); subdet++) {
       const CaloSubdetectorGeometry* sdg = geom_->getSubdetectorGeometry(detector_, subdet);
       if (sdg != nullptr) {
         // get the list of detids within range (from geometry)
         CaloSubdetectorGeometry::DetIdSet dis_excl = sdg->getCells(p, deltaRmin_);
         CaloSubdetectorGeometry::DetIdSet dis_all = sdg->getCells(p, deltaRmax_);
         // use set operations to determine detids in annulus
         CaloSubdetectorGeometry::DetIdSet dis;
         std::set_difference(
             dis_all.begin(), dis_all.end(), dis_excl.begin(), dis_excl.end(), std::inserter(dis, dis.begin()));
  
         // loop over detids...
         typename edm::SortedCollection<T>::const_iterator j, je = inputCollection.end();
  
         for (CaloSubdetectorGeometry::DetIdSet::iterator i = dis.begin(); i != dis.end(); i++) {
           if (i->subdetId() != subdet)
             continue;  // possible for HCAL where the same geometry object handles all the detectors
           j = inputCollection.find(*i);
           if (j != je)
             callback(*j);
         }
       }
     }
   }

◆ selectCallback() [2/2]

template<class T>

void CaloDualConeSelector< T >::selectCallback	(	double	eta,
		double	phi,
		const edm::SortedCollection< T > &	inputCollection,
		std::function< void(const T &)>	callback
	)

inline

Definition at line 22 of file CaloDualConeSelector.h.

                                                                    {
     GlobalPoint p(GlobalPoint::Cylindrical(1, phi, tanh(eta)));
     return selectCallback(p, inputCollection, callback);
   }