#include <DetIdAssociator.h>

Inheritance diagram for DetIdAssociator:

Classes
struct	MapRange

Public Types
typedef std::vector < GlobalPoint > ::const_iterator	const_iterator

enum	PropagationTarget { Barrel, ForwardEndcap, BackwardEndcap }

Public Member Functions
virtual void	buildMap ()
	make the look-up map More...

	DetIdAssociator (const int nPhi, const int nEta, const double etaBinSize)

double	etaBinSize () const
	look-up map bin size in eta dimension More...

virtual std::vector< DetId >	getCrossedDetIds (const std::set< DetId > &, const std::vector< GlobalPoint > &trajectory) const

virtual std::vector< DetId >	getCrossedDetIds (const std::set< DetId > &, const std::vector< SteppingHelixStateInfo > &trajectory, const double toleranceInSigmas=-1) const

virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &, const int iN=0) const

virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &direction, const unsigned int iNEtaPlus, const unsigned int iNEtaMinus, const unsigned int iNPhiPlus, const unsigned int iNPhiMinus) const

virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &direction, const MapRange &mapRange) const

virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &point, const double d=0) const

virtual std::set< DetId >	getDetIdsCloseToAPoint (const GlobalPoint &point, const double dThetaPlus, const double dThetaMinus, const double dPhiPlus, const double dPhiMinus) const

virtual std::set< DetId >	getDetIdsInACone (const std::set< DetId > &, const std::vector< GlobalPoint > &trajectory, const double dR) const

virtual const GeomDet *	getGeomDet (const DetId &) const =0

virtual int	iEta (const GlobalPoint &) const
	look-up map eta index More...

virtual int	iPhi (const GlobalPoint &) const
	look-up map phi index More...

virtual const char *	name () const =0

int	nEtaBins () const
	number of bins of the look-up map in eta dimension More...

int	nPhiBins () const
	number of bins of the look-up map in phi dimension More...

virtual bool	selectAllInACone (const double dR) const
	helper to see if getDetIdsInACone is useful More...

const FiducialVolume &	volume () const
	get active detector volume More...

virtual	~DetIdAssociator ()

Protected Member Functions
virtual void	check_setup () const

virtual bool	crossedElement (const GlobalPoint &, const GlobalPoint &, const DetId &, const double toleranceInSigmas=-1, const SteppingHelixStateInfo *=nullptr) const

virtual void	dumpMapContent (int, int) const

virtual void	dumpMapContent (int, int, int, int) const

void	fillSet (std::set< DetId > &set, unsigned int iEta, unsigned int iPhi) const

virtual std::pair < const_iterator, const_iterator >	getDetIdPoints (const DetId &, std::vector< GlobalPoint > &) const =0

virtual const unsigned int	getNumberOfSubdetectors () const

virtual GlobalPoint	getPosition (const DetId &) const =0

virtual void	getValidDetIds (unsigned int subDetectorIndex, std::vector< DetId > &) const =0

unsigned int	index (unsigned int iEta, unsigned int iPhi) const

virtual bool	insideElement (const GlobalPoint &, const DetId &) const =0

virtual bool	nearElement (const GlobalPoint &point, const DetId &id, const double distance) const

Protected Attributes
std::vector< DetId >	container_

const double	etaBinSize_

std::vector< std::pair < unsigned int, unsigned int > >	lookupMap_

double	maxEta_

double	minTheta_

const int	nEta_

const int	nPhi_

bool	theMapIsValid_

FiducialVolume	volume_

Detailed Description

\

Description: Abstract base class for 3D point -> std::set<DetId>

Implementation: A look up map of active detector elements in eta-phi space is built to speed up access to the detector element geometry as well as associated hits. The map is uniformly binned in eta and phi dimensions, which can be viewed as a histogram with every bin containing DetId of elements crossed in a given eta-phi window. It is very likely that a single DetId can be found in a few bins if it's geometrical size is bigger than eta-phi bin size.

The map is implemented as a double array. The first one has fixed size and points to the range of array elements in the second one.

Definition at line 45 of file DetIdAssociator.h.

Member Typedef Documentation

typedef std::vector<GlobalPoint>::const_iterator DetIdAssociator::const_iterator

Definition at line 54 of file DetIdAssociator.h.

Member Enumeration Documentation

enum DetIdAssociator::PropagationTarget

Enumerator
Barrel
ForwardEndcap
BackwardEndcap

Definition at line 47 of file DetIdAssociator.h.

47 { Barrel, ForwardEndcap, BackwardEndcap };

DetIdAssociator::Barrel

Definition: DetIdAssociator.h:47

DetIdAssociator::BackwardEndcap

Definition: DetIdAssociator.h:47

DetIdAssociator::ForwardEndcap

Definition: DetIdAssociator.h:47

Constructor & Destructor Documentation

DetIdAssociator::DetIdAssociator	(	const int	nPhi,
		const int	nEta,
		const double	etaBinSize
	)

Definition at line 24 of file DetIdAssociator.cc.

References etaBinSize_, Exception, funct::exp(), maxEta_, minTheta_, nEta_, and nPhi_.

     : nPhi_(nPhi),
       nEta_(nEta),
       lookupMap_(nPhi_ * nEta_, std::pair<unsigned int, unsigned int>(0, 0)),
       theMapIsValid_(false),
       etaBinSize_(etaBinSize) {
   if (nEta_ <= 0 || nPhi_ <= 0)
     throw cms::Exception("FatalError") << "incorrect look-up map size. Cannot initialize such a map.";
   maxEta_ = etaBinSize_ * nEta_ / 2;
   minTheta_ = 2 * atan(exp(-maxEta_));
 }

virtual DetIdAssociator::~DetIdAssociator ( )

inlinevirtual

Definition at line 57 of file DetIdAssociator.h.

57 {}

Member Function Documentation

void DetIdAssociator::buildMap ( )

virtual

make the look-up map

Definition at line 138 of file DetIdAssociator.cc.

References funct::abs(), FiducialVolume::addActivePoint(), cms::cuda::assert(), newFWLiteAna::bin, check_setup(), container_, FiducialVolume::determinInnerDimensions(), etaBinSize_, Puppi_cff::etaMax, Puppi_cff::etaMin, Exception, getDetIdPoints(), getNumberOfSubdetectors(), getPosition(), getValidDetIds(), iEta(), index(), DetIdInfo::info(), iPhi(), edm::isNotFinite(), LogTrace, lookupMap_, FiducialVolume::maxR(), FiducialVolume::maxZ(), FiducialVolume::minR(), FiducialVolume::minZ(), eostools::move(), name(), nEta_, nPhi_, L1TMuonDQMOffline_cfi::phiMax, L1TMuonDQMOffline_cfi::phiMin, theMapIsValid_, trackerHitRTTI::vector, and volume_.

                                {
   // the map is built in two steps to create a clean memory footprint
   // 0) determine how many elements each bin has
   // 1) fill the container map
   check_setup();
   LogTrace("TrackAssociator") << "building map for " << name() << "\n";
   // clear the map
   if (nEta_ <= 0 || nPhi_ <= 0)
     throw cms::Exception("FatalError") << "incorrect look-up map size. Cannot build such a map.";
   std::vector<GlobalPoint> pointBuffer;
   std::vector<DetId> detIdBuffer;
   unsigned int numberOfSubDetectors = getNumberOfSubdetectors();
   unsigned totalNumberOfElementsInTheContainer(0);
   for (unsigned int step = 0; step < 2; ++step) {
     unsigned int numberOfDetIdsOutsideEtaRange = 0;
     unsigned int numberOfDetIdsActive = 0;
     LogTrace("TrackAssociator") << "Step: " << step;
     for (unsigned int subDetectorIndex = 0; subDetectorIndex < numberOfSubDetectors; ++subDetectorIndex) {
       getValidDetIds(subDetectorIndex, detIdBuffer);
       // This std::move prevents modification
       std::vector<DetId> const& validIds = std::move(detIdBuffer);
       LogTrace("TrackAssociator") << "Number of valid DetIds for subdetector: " << subDetectorIndex << " is "
                                   << validIds.size();
       for (std::vector<DetId>::const_iterator id_itr = validIds.begin(); id_itr != validIds.end(); id_itr++) {
         std::pair<const_iterator, const_iterator> points = getDetIdPoints(*id_itr, pointBuffer);
         LogTrace("TrackAssociatorVerbose") << "Found " << points.second - points.first
                                            << " global points to describe geometry of DetId: " << id_itr->rawId();
         int etaMax(-1);
         int etaMin(-1);
         int phiMax(-1);
         int phiMin(-1);
         // this is a bit overkill, but it should be 100% proof (when debugged :)
         for (std::vector<GlobalPoint>::const_iterator iter = points.first; iter != points.second; iter++) {
           LogTrace("TrackAssociatorVerbose")
               << "\tpoint (rho,phi,z): " << iter->perp() << ", " << iter->phi() << ", " << iter->z();
           // FIX ME: this should be a fatal error
           if (edm::isNotFinite(iter->mag()) || iter->mag() > 1e5) {  //Detector parts cannot be 1 km away or be NaN
             edm::LogWarning("TrackAssociator")
                 << "Critical error! Bad detector unit geometry:\n\tDetId:" << id_itr->rawId()
                 << "\t mag(): " << iter->mag() << "\n"
                 << DetIdInfo::info(*id_itr, nullptr) << "\nSkipped the element";
             continue;
           }
           volume_.addActivePoint(*iter);
           int ieta = iEta(*iter);
           int iphi = iPhi(*iter);
           if (ieta < 0 || ieta >= nEta_) {
             LogTrace("TrackAssociator") << "Out of range: DetId:" << id_itr->rawId() << "\t (ieta,iphi): " << ieta
                                         << "," << iphi << "\n"
                                         << "Point: " << *iter << "\t(eta,phi): " << (*iter).eta() << ","
                                         << (*iter).phi() << "\n center: " << getPosition(*id_itr);
             continue;
           }
           if (phiMin < 0) {
             // first element
             etaMin = ieta;
             etaMax = ieta;
             phiMin = iphi;
             phiMax = iphi;
           } else {
             // check for discontinuity in phi
             int deltaMin = abs(phiMin - iphi);
             int deltaMax = abs(phiMax - iphi);
             // assume that no single detector element has more than 3.1416 coverage in phi
             if (deltaMin > nPhi_ / 2 && phiMin < nPhi_ / 2)
               phiMin += nPhi_;
             if (deltaMax > nPhi_ / 2) {
               if (phiMax < nPhi_ / 2)
                 phiMax += nPhi_;
               else
                 iphi += nPhi_;
             }
             assert(iphi >= 0);
             if (etaMin > ieta)
               etaMin = ieta;
             if (etaMax < ieta)
               etaMax = ieta;
             if (phiMin > iphi)
               phiMin = iphi;
             if (phiMax < iphi)
               phiMax = iphi;
           }
         }
         if (etaMax < 0 || phiMax < 0 || etaMin >= nEta_ || phiMin >= nPhi_) {
           LogTrace("TrackAssociatorVerbose")
               << "Out of range or no geometry: DetId:" << id_itr->rawId() << "\n\teta (min,max): " << etaMin << ","
               << etaMax << "\n\tphi (min,max): " << phiMin << "," << phiMax << "\nTower id: " << id_itr->rawId()
               << "\n";
           numberOfDetIdsOutsideEtaRange++;
           continue;
         }
         numberOfDetIdsActive++;
 
         LogTrace("TrackAssociatorVerbose") << "DetId (ieta_min,ieta_max,iphi_min,iphi_max): " << id_itr->rawId() << ", "
                                            << etaMin << ", " << etaMax << ", " << phiMin << ", " << phiMax;
         for (int ieta = etaMin; ieta <= etaMax; ieta++)
           for (int iphi = phiMin; iphi <= phiMax; iphi++)
             if (step == 0) {
               lookupMap_.at(index(ieta, iphi % nPhi_)).second++;
               totalNumberOfElementsInTheContainer++;
             } else {
               container_.at(lookupMap_.at(index(ieta, iphi % nPhi_)).first) = *id_itr;
               lookupMap_.at(index(ieta, iphi % nPhi_)).first--;
               totalNumberOfElementsInTheContainer--;
             }
       }
     }
     LogTrace("TrackAssociator") << "Number of elements outside the allowed range ( |eta|>" << nEta_ / 2 * etaBinSize_
                                 << "): " << numberOfDetIdsOutsideEtaRange << "\n";
     LogTrace("TrackAssociator") << "Number of active DetId's mapped: " << numberOfDetIdsActive << "\n";
     if (step == 0) {
       // allocate
       container_.resize(totalNumberOfElementsInTheContainer);
       // fill the range index in the lookup map to point to the last element in the range
       unsigned int index(0);
       for (std::vector<std::pair<unsigned int, unsigned int> >::iterator bin = lookupMap_.begin();
            bin != lookupMap_.end();
            ++bin) {
         if (bin->second == 0)
           continue;
         index += bin->second;
         bin->first = index - 1;
       }
     }
   }
   if (totalNumberOfElementsInTheContainer != 0)
     throw cms::Exception("FatalError")
         << "Look-up map filled incorrectly. Structural problem. Get in touch with the developer.";
   volume_.determinInnerDimensions();
   edm::LogVerbatim("TrackAssociator") << "Fiducial volume for " << name()
                                       << " (minR, maxR, minZ, maxZ): " << volume_.minR() << ", " << volume_.maxR()
                                       << ", " << volume_.minZ() << ", " << volume_.maxZ();
   theMapIsValid_ = true;
 }

void DetIdAssociator::check_setup ( ) const

protectedvirtual

Reimplemented in MuonDetIdAssociator, and CaloDetIdAssociator.

Definition at line 371 of file DetIdAssociator.cc.

References etaBinSize_, Exception, nEta_, and nPhi_.

Referenced by buildMap(), CaloDetIdAssociator::check_setup(), MuonDetIdAssociator::check_setup(), getCrossedDetIds(), getDetIdsCloseToAPoint(), and getDetIdsInACone().

                                         {
   if (nEta_ == 0)
     throw cms::Exception("FatalError") << "Number of eta bins is not set.\n";
   if (nPhi_ == 0)
     throw cms::Exception("FatalError") << "Number of phi bins is not set.\n";
   // if (ivProp_==0) throw cms::Exception("FatalError") << "Track propagator is not defined\n";
   if (etaBinSize_ == 0)
     throw cms::Exception("FatalError") << "Eta bin size is not set.\n";
 }

virtual bool DetIdAssociator::crossedElement	(	const GlobalPoint &	,
		const GlobalPoint &	,
		const DetId &	,
		const double	toleranceInSigmas = `-1`,
		const SteppingHelixStateInfo *	= `nullptr`
	)		const

inlineprotectedvirtual

Reimplemented in CaloDetIdAssociator.

Definition at line 130 of file DetIdAssociator.h.

Referenced by getCrossedDetIds().

                                                                              {
     return false;
   }

void DetIdAssociator::dumpMapContent	(	int	ieta,
		int	iphi
	)		const

protectedvirtual

Definition at line 327 of file DetIdAssociator.cc.

References fillSet(), getDetIdPoints(), LogTrace, nPhi_, and point.

Referenced by dumpMapContent().

                                                              {
   if (!(ieta >= 0 && ieta < nEta_ && iphi >= 0)) {
     edm::LogWarning("TrackAssociator") << "ieta or iphi is out of range. Skipped.";
     return;
   }
 
   std::vector<GlobalPoint> pointBuffer;
   std::set<DetId> set;
   fillSet(set, ieta, iphi % nPhi_);
   LogTrace("TrackAssociator") << "Map content for cell (ieta,iphi): " << ieta << ", " << iphi % nPhi_;
   for (std::set<DetId>::const_iterator itr = set.begin(); itr != set.end(); itr++) {
     LogTrace("TrackAssociator") << "\tDetId " << itr->rawId() << ", geometry (x,y,z,rho,eta,phi):";
     std::pair<const_iterator, const_iterator> points = getDetIdPoints(*itr, pointBuffer);
     for (std::vector<GlobalPoint>::const_iterator point = points.first; point != points.second; point++)
       LogTrace("TrackAssociator") << "\t\t" << point->x() << ", " << point->y() << ", " << point->z() << ", "
                                   << point->perp() << ", " << point->eta() << ", " << point->phi();
   }
 }

void DetIdAssociator::dumpMapContent	(	int	ieta_min,
		int	ieta_max,
		int	iphi_min,
		int	iphi_max
	)		const

protectedvirtual

Definition at line 346 of file DetIdAssociator.cc.

References dumpMapContent(), mps_fire::i, and dqmiolumiharvest::j.

                                                                                                  {
   for (int i = ieta_min; i <= ieta_max; i++)
     for (int j = iphi_min; j <= iphi_max; j++)
       dumpMapContent(i, j);
 }

double DetIdAssociator::etaBinSize ( ) const

inline

look-up map bin size in eta dimension

Definition at line 108 of file DetIdAssociator.h.

References etaBinSize_.

108 { return etaBinSize_; };

DetIdAssociator::etaBinSize_

const double etaBinSize_

Definition: DetIdAssociator.h:151

void DetIdAssociator::fillSet	(	std::set< DetId > &	set,
		unsigned int	iEta,
		unsigned int	iPhi
	)		const

protected

Definition at line 381 of file DetIdAssociator.cc.

References container_, mps_fire::i, index(), lookupMap_, and findQualityFiles::size.

Referenced by dumpMapContent(), and getDetIdsCloseToAPoint().

                                                                                             {
   unsigned int i = index(iEta, iPhi);
   unsigned int i0 = lookupMap_.at(i).first;
   unsigned int size = lookupMap_.at(i).second;
   for (i = i0; i < i0 + size; ++i)
     set.insert(container_.at(i));
 }

std::vector< DetId > DetIdAssociator::getCrossedDetIds	(	const std::set< DetId > &	inset,
		const std::vector< GlobalPoint > &	trajectory
	)		const

virtual

DetIds crossed by the track tolerance is the radius of the trajectory used for matching -1 is default and represent the case with no uncertainty on the trajectory direction. It's the fastest option
DetIds crossed by the track, ordered according to the order that they were crossed by the track flying outside the detector

Definition at line 290 of file DetIdAssociator.cc.

References check_setup(), crossedElement(), mps_fire::i, and convertSQLitetoXML_cfg::output.

                                                                                                      {
   check_setup();
   std::vector<DetId> output;
   std::set<DetId> ids(inset);
   for (unsigned int i = 0; i + 1 < trajectory.size(); ++i) {
     std::set<DetId>::const_iterator id_iter = ids.begin();
     while (id_iter != ids.end()) {
       if (crossedElement(trajectory[i], trajectory[i + 1], *id_iter)) {
         output.push_back(*id_iter);
         ids.erase(id_iter++);
       } else
         id_iter++;
     }
   }
   return output;
 }

std::vector< DetId > DetIdAssociator::getCrossedDetIds	(	const std::set< DetId > &	inset,
		const std::vector< SteppingHelixStateInfo > &	trajectory,
		const double	toleranceInSigmas = `-1`
	)		const

virtual

Definition at line 308 of file DetIdAssociator.cc.

References check_setup(), crossedElement(), mps_fire::i, convertSQLitetoXML_cfg::output, and position.

                                                                                    {
   check_setup();
   std::vector<DetId> output;
   std::set<DetId> ids(inset);
   for (unsigned int i = 0; i + 1 < trajectory.size(); ++i) {
     std::set<DetId>::const_iterator id_iter = ids.begin();
     while (id_iter != ids.end()) {
       if (crossedElement(trajectory[i].position(), trajectory[i + 1].position(), *id_iter, tolerance, &trajectory[i])) {
         output.push_back(*id_iter);
         ids.erase(id_iter++);
       } else
         id_iter++;
     }
   }
   return output;
 }

virtual std::pair<const_iterator, const_iterator> DetIdAssociator::getDetIdPoints	(	const DetId &	,
		std::vector< GlobalPoint > &
	)		const

protectedpure virtual

Implemented in MuonDetIdAssociator, and CaloDetIdAssociator.

Referenced by buildMap(), and dumpMapContent().

std::set< DetId > DetIdAssociator::getDetIdsCloseToAPoint	(	const GlobalPoint &	direction,
		const int	iN = `0`
	)		const

virtual

Preselect DetIds close to a point on the inner surface of the detector. "iN" is a number of the adjacent bins of the map to retrieve

Definition at line 36 of file DetIdAssociator.cc.

References dqmiodumpmetadata::n.

Referenced by getDetIdsCloseToAPoint().

                                                                                                       {
   unsigned int n = 0;
   if (iN > 0)
     n = iN;
   return getDetIdsCloseToAPoint(direction, n, n, n, n);
 }

std::set< DetId > DetIdAssociator::getDetIdsCloseToAPoint	(	const GlobalPoint &	direction,
		const unsigned int	iNEtaPlus,
		const unsigned int	iNEtaMinus,
		const unsigned int	iNPhiPlus,
		const unsigned int	iNPhiMinus
	)		const

virtual

Definition at line 43 of file DetIdAssociator.cc.

References check_setup(), PV3DBase< T, PVType, FrameType >::eta(), Exception, fillSet(), mps_fire::i, iEta(), iPhi(), dqmiolumiharvest::j, LogTrace, nEta_, nPhi_, PV3DBase< T, PVType, FrameType >::phi(), and theMapIsValid_.

                                                                                              {
   std::set<DetId> set;
   check_setup();
   if (!theMapIsValid_)
     throw cms::Exception("FatalError") << "map is not valid.";
   LogTrace("TrackAssociator") << "(iNEtaPlus, iNEtaMinus, iNPhiPlus, iNPhiMinus): " << iNEtaPlus << ", " << iNEtaMinus
                               << ", " << iNPhiPlus << ", " << iNPhiMinus;
   LogTrace("TrackAssociator") << "point (eta,phi): " << direction.eta() << "," << direction.phi();
   int ieta = iEta(direction);
   int iphi = iPhi(direction);
   LogTrace("TrackAssociator") << "(ieta,iphi): " << ieta << "," << iphi << "\n";
   if (ieta >= 0 && ieta < nEta_ && iphi >= 0 && iphi < nPhi_) {
     fillSet(set, ieta, iphi);
     // dumpMapContent(ieta,iphi);
     // check if any neighbor bin is requested
     if (iNEtaPlus + iNEtaMinus + iNPhiPlus + iNPhiMinus > 0) {
       LogTrace("TrackAssociator") << "Add neighbors (ieta,iphi): " << ieta << "," << iphi;
       // eta
       int maxIEta = ieta + iNEtaPlus;
       int minIEta = ieta - iNEtaMinus;
       if (maxIEta >= nEta_)
         maxIEta = nEta_ - 1;
       if (minIEta < 0)
         minIEta = 0;
       // phi
       int maxIPhi = iphi + iNPhiPlus;
       int minIPhi = iphi - iNPhiMinus;
       if (maxIPhi - minIPhi >= nPhi_) {  // all elements in phi
         minIPhi = 0;
         maxIPhi = nPhi_ - 1;
       }
       if (minIPhi < 0) {
         minIPhi += nPhi_;
         maxIPhi += nPhi_;
       }
       LogTrace("TrackAssociator") << "\tieta (min,max): " << minIEta << "," << maxIEta;
       LogTrace("TrackAssociator") << "\tiphi (min,max): " << minIPhi << "," << maxIPhi << "\n";
       // dumpMapContent(minIEta,maxIEta,minIPhi,maxIPhi);
       for (int i = minIEta; i <= maxIEta; i++)
         for (int j = minIPhi; j <= maxIPhi; j++) {
           if (i == ieta && j == iphi)
             continue;  // already in the set
           fillSet(set, i, j % nPhi_);
         }
     }
   }
   return set;
 }

std::set< DetId > DetIdAssociator::getDetIdsCloseToAPoint	(	const GlobalPoint &	direction,
		const MapRange &	mapRange
	)		const

virtual

Definition at line 358 of file DetIdAssociator.cc.

References DetIdAssociator::MapRange::dPhiMinus, DetIdAssociator::MapRange::dPhiPlus, DetIdAssociator::MapRange::dThetaMinus, DetIdAssociator::MapRange::dThetaPlus, and getDetIdsCloseToAPoint().

                                                                                                                   {
   return getDetIdsCloseToAPoint(
       direction, mapRange.dThetaPlus, mapRange.dThetaMinus, mapRange.dPhiPlus, mapRange.dPhiMinus);
 }

std::set< DetId > DetIdAssociator::getDetIdsCloseToAPoint	(	const GlobalPoint &	point,
		const double	d = `0`
	)		const

virtual

Preselect DetIds close to a point on the inner surface of the detector. "d" defines the allowed range in theta-phi space:

theta is in [point.theta()-d, point.theta()+d]
phi is in [point.phi()-d, point.phi()+d]

Definition at line 96 of file DetIdAssociator.cc.

References getDetIdsCloseToAPoint().

                                                                                                     {
   return getDetIdsCloseToAPoint(point, d, d, d, d);
 }

std::set< DetId > DetIdAssociator::getDetIdsCloseToAPoint	(	const GlobalPoint &	point,
		const double	dThetaPlus,
		const double	dThetaMinus,
		const double	dPhiPlus,
		const double	dPhiMinus
	)		const

virtual

theta is in [point.theta()-dThetaMinus, point.theta()+dThetaPlus]
phi is in [point.phi()-dPhiMinus, point.phi()+dPhiPlus]

Definition at line 100 of file DetIdAssociator.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::eta(), etaBinSize_, getDetIdsCloseToAPoint(), log, LogTrace, M_PI, maxEta, minTheta_, dqmiodumpmetadata::n, nPhi_, funct::tan(), and PV3DBase< T, PVType, FrameType >::theta().

                                                                                       {
   LogTrace("TrackAssociator") << "(dThetaPlus,dThetaMinus,dPhiPlus,dPhiMinus): " << dThetaPlus << ", " << dThetaMinus
                               << ", " << dPhiPlus << ", " << dPhiMinus;
   unsigned int n = 0;
   if (dThetaPlus < 0 || dThetaMinus < 0 || dPhiPlus < 0 || dPhiMinus < 0)
     return getDetIdsCloseToAPoint(point, n, n, n, n);
   // check that region of interest overlaps with the look-up map
   double maxTheta = point.theta() + dThetaPlus;
   if (maxTheta > M_PI - minTheta_)
     maxTheta = M_PI - minTheta_;
   double minTheta = point.theta() - dThetaMinus;
   if (minTheta < minTheta_)
     minTheta = minTheta_;
   if (maxTheta < minTheta_ || minTheta > M_PI - minTheta_)
     return std::set<DetId>();
 
   // take into account non-linear dependence of eta from
   // theta in regions with large |eta|
   double minEta = -log(tan(maxTheta / 2));
   double maxEta = -log(tan(minTheta / 2));
   unsigned int iNEtaPlus = abs(int((maxEta - point.eta()) / etaBinSize_));
   unsigned int iNEtaMinus = abs(int((point.eta() - minEta) / etaBinSize_));
   unsigned int iNPhiPlus = abs(int(dPhiPlus / (2 * M_PI) * nPhi_));
   unsigned int iNPhiMinus = abs(int(dPhiMinus / (2 * M_PI) * nPhi_));
   // add one more bin in each direction to guaranty that we don't miss anything
   return getDetIdsCloseToAPoint(point, iNEtaPlus + 1, iNEtaMinus + 1, iNPhiPlus + 1, iNPhiMinus + 1);
 }

std::set< DetId > DetIdAssociator::getDetIdsInACone	(	const std::set< DetId > &	inset,
		const std::vector< GlobalPoint > &	trajectory,
		const double	dR
	)		const

virtual

Find DetIds that satisfy given requirements

inside eta-phi cone of radius dR

Definition at line 273 of file DetIdAssociator.cc.

References check_setup(), nearElement(), and selectAllInACone().

                                                                          {
   if (selectAllInACone(dR))
     return inset;
   check_setup();
   std::set<DetId> outset;
   for (std::set<DetId>::const_iterator id_iter = inset.begin(); id_iter != inset.end(); id_iter++)
     for (std::vector<GlobalPoint>::const_iterator point_iter = trajectory.begin(); point_iter != trajectory.end();
          point_iter++)
       if (nearElement(*point_iter, *id_iter, dR)) {
         outset.insert(*id_iter);
         break;
       }
   return outset;
 }

virtual const GeomDet* DetIdAssociator::getGeomDet ( const DetId & ) const

pure virtual

Implemented in MuonDetIdAssociator, and CaloDetIdAssociator.

virtual const unsigned int DetIdAssociator::getNumberOfSubdetectors ( ) const

inlineprotectedvirtual

Reimplemented in HcalDetIdAssociator, and EcalDetIdAssociator.

Definition at line 125 of file DetIdAssociator.h.

Referenced by buildMap().

125 { return 1; }

virtual GlobalPoint DetIdAssociator::getPosition ( const DetId & ) const

protectedpure virtual

Implemented in MuonDetIdAssociator, and CaloDetIdAssociator.

Referenced by buildMap(), and nearElement().

virtual void DetIdAssociator::getValidDetIds	(	unsigned int	subDetectorIndex,
		std::vector< DetId > &
	)		const

protectedpure virtual

Implemented in MuonDetIdAssociator, CaloDetIdAssociator, HcalDetIdAssociator, EcalDetIdAssociator, HODetIdAssociator, and PreshowerDetIdAssociator.

Referenced by buildMap().

int DetIdAssociator::iEta ( const GlobalPoint & point ) const

virtual

look-up map eta index

Definition at line 132 of file DetIdAssociator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), etaBinSize_, and nEta_.

Referenced by buildMap(), and getDetIdsCloseToAPoint().

132 { return int(point.eta() / etaBinSize_ + nEta_ / 2); }

DetIdAssociator::etaBinSize_

const double etaBinSize_

Definition: DetIdAssociator.h:151

PV3DBase::eta

T eta() const

Definition: PV3DBase.h:73

DetIdAssociator::nEta_

const int nEta_

Definition: DetIdAssociator.h:144

unsigned int DetIdAssociator::index	(	unsigned int	iEta,
		unsigned int	iPhi
	)		const

inlineprotected

Definition at line 139 of file DetIdAssociator.h.

References iPhi(), and nPhi_.

Referenced by BeautifulSoup.PageElement::_invert(), buildMap(), and fillSet().

139 { return iEta * nPhi_ + iPhi; }

DetIdAssociator::nPhi_

const int nPhi_

Definition: DetIdAssociator.h:143

DetIdAssociator::iEta

virtual int iEta(const GlobalPoint &) const

look-up map eta index

Definition: DetIdAssociator.cc:132

DetIdAssociator::iPhi

virtual int iPhi(const GlobalPoint &) const

look-up map phi index

Definition: DetIdAssociator.cc:134

virtual bool DetIdAssociator::insideElement	(	const GlobalPoint &	,
		const DetId &
	)		const

protectedpure virtual

Implemented in MuonDetIdAssociator, and CaloDetIdAssociator.

int DetIdAssociator::iPhi ( const GlobalPoint & point ) const

virtual

look-up map phi index

Definition at line 134 of file DetIdAssociator.cc.

References M_PI, nPhi_, and PV3DBase< T, PVType, FrameType >::phi().

Referenced by buildMap(), getDetIdsCloseToAPoint(), and index().

                                                         {
   return int((double(point.phi()) + M_PI) / (2 * M_PI) * nPhi_);
 }

virtual const char* DetIdAssociator::name ( ) const

pure virtual

Implemented in MuonDetIdAssociator, CaloDetIdAssociator, HcalDetIdAssociator, EcalDetIdAssociator, HODetIdAssociator, and PreshowerDetIdAssociator.

Referenced by config.CFG::__str__(), buildMap(), geometryComparison.GeometryComparison::createScript(), validation.Sample::digest(), primaryVertexResolution.PrimaryVertexResolution::getRepMap(), primaryVertexValidation.PrimaryVertexValidation::getRepMap(), and zMuMuValidation.ZMuMuValidation::getRepMap().

bool DetIdAssociator::nearElement	(	const GlobalPoint &	point,
		const DetId &	id,
		const double	distance
	)		const

protectedvirtual

Definition at line 363 of file DetIdAssociator.cc.

References srCondWrite_cfg::deltaPhi, PV3DBase< T, PVType, FrameType >::eta(), getPosition(), M_PI, and PV3DBase< T, PVType, FrameType >::phi().

Referenced by getDetIdsInACone().

                                                                                                         {
   GlobalPoint center = getPosition(id);
   double deltaPhi(fabs(point.phi() - center.phi()));
   if (deltaPhi > M_PI)
     deltaPhi = fabs(deltaPhi - M_PI * 2.);
   return (point.eta() - center.eta()) * (point.eta() - center.eta()) + deltaPhi * deltaPhi < distance * distance;
 }

int DetIdAssociator::nEtaBins ( ) const

inline

number of bins of the look-up map in eta dimension

Definition at line 106 of file DetIdAssociator.h.

References nEta_.

106 { return nEta_; }

DetIdAssociator::nEta_

const int nEta_

Definition: DetIdAssociator.h:144

int DetIdAssociator::nPhiBins ( ) const

inline

number of bins of the look-up map in phi dimension

Definition at line 104 of file DetIdAssociator.h.

References nPhi_.

104 { return nPhi_; }

DetIdAssociator::nPhi_

const int nPhi_

Definition: DetIdAssociator.h:143

virtual bool DetIdAssociator::selectAllInACone ( const double dR ) const

inlinevirtual

helper to see if getDetIdsInACone is useful

Definition at line 82 of file DetIdAssociator.h.

References M_PI, and maxEta_.

Referenced by getDetIdsInACone().

82 { return dR > 2 * M_PI && dR > maxEta_; }

DetIdAssociator::maxEta_

double maxEta_

Definition: DetIdAssociator.h:152

M_PI

#define M_PI

Definition: BXVectorInputProducer.cc:50

const FiducialVolume & DetIdAssociator::volume ( void ) const

get active detector volume

Definition at line 352 of file DetIdAssociator.cc.

References Exception, theMapIsValid_, and volume_.

                                                     {
   if (!theMapIsValid_)
     throw cms::Exception("FatalError") << "map is not valid.";
   return volume_;
 }

Member Data Documentation

std::vector<DetId> DetIdAssociator::container_

protected

Definition at line 149 of file DetIdAssociator.h.

Referenced by buildMap(), and fillSet().

const double DetIdAssociator::etaBinSize_

protected

Definition at line 151 of file DetIdAssociator.h.

Referenced by buildMap(), check_setup(), DetIdAssociator(), etaBinSize(), getDetIdsCloseToAPoint(), and iEta().

std::vector<std::pair<unsigned int, unsigned int> > DetIdAssociator::lookupMap_

protected

Definition at line 148 of file DetIdAssociator.h.

Referenced by buildMap(), and fillSet().

double DetIdAssociator::maxEta_

protected

Definition at line 152 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and selectAllInACone().

double DetIdAssociator::minTheta_

protected

Definition at line 153 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and getDetIdsCloseToAPoint().

const int DetIdAssociator::nEta_

protected

Definition at line 144 of file DetIdAssociator.h.

Referenced by buildMap(), check_setup(), DetIdAssociator(), getDetIdsCloseToAPoint(), iEta(), and nEtaBins().

const int DetIdAssociator::nPhi_

protected

Definition at line 143 of file DetIdAssociator.h.

Referenced by buildMap(), check_setup(), DetIdAssociator(), dumpMapContent(), getDetIdsCloseToAPoint(), index(), iPhi(), and nPhiBins().

bool DetIdAssociator::theMapIsValid_

protected

Definition at line 150 of file DetIdAssociator.h.

Referenced by buildMap(), getDetIdsCloseToAPoint(), and volume().

FiducialVolume DetIdAssociator::volume_

protected

Definition at line 158 of file DetIdAssociator.h.

Referenced by buildMap(), and volume().

Classes

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation