#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 MapRange &mapRange) 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 &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

◆ const_iterator

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

Definition at line 54 of file DetIdAssociator.h.

Member Enumeration Documentation

◆ PropagationTarget

enum DetIdAssociator::PropagationTarget

Enumerator
Barrel
ForwardEndcap
BackwardEndcap

Definition at line 47 of file DetIdAssociator.h.

47 { Barrel, ForwardEndcap, BackwardEndcap };

Constructor & Destructor Documentation

◆ DetIdAssociator()

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

Definition at line 24 of file DetIdAssociator.cc.

     : 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_));
 }

References etaBinSize_, Exception, JetChargeProducer_cfi::exp, maxEta_, minTheta_, nEta_, and nPhi_.

◆ ~DetIdAssociator()

virtual DetIdAssociator::~DetIdAssociator ( )

inlinevirtual

Definition at line 57 of file DetIdAssociator.h.

57 {}

Member Function Documentation

◆ buildMap()

void DetIdAssociator::buildMap ( )

virtual

make the look-up map

Definition at line 138 of file DetIdAssociator.cc.

                                {
   // 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;
 }

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

◆ check_setup()

void DetIdAssociator::check_setup ( ) const

protectedvirtual

Reimplemented in MuonDetIdAssociator, and CaloDetIdAssociator.

Definition at line 371 of file DetIdAssociator.cc.

                                         {
   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";
 }

References etaBinSize_, Exception, nEta_, and nPhi_.

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

◆ crossedElement()

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.

                                                                              {
     return false;
   }

Referenced by getCrossedDetIds().

◆ dumpMapContent() [1/2]

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

protectedvirtual

Definition at line 327 of file DetIdAssociator.cc.

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

References fillSet(), getDetIdPoints(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, LogTrace, nPhi_, point, and HLT_2018_cff::points.

Referenced by dumpMapContent().

◆ dumpMapContent() [2/2]

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.

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

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

◆ etaBinSize()

double DetIdAssociator::etaBinSize ( ) const

inline

look-up map bin size in eta dimension

Definition at line 108 of file DetIdAssociator.h.

108 { return etaBinSize_; };

References etaBinSize_.

◆ fillSet()

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

protected

Definition at line 381 of file DetIdAssociator.cc.

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

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

Referenced by dumpMapContent(), and getDetIdsCloseToAPoint().

◆ getCrossedDetIds() [1/2]

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.

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

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

◆ getCrossedDetIds() [2/2]

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.

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

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

◆ getDetIdPoints()

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().

◆ getDetIdsCloseToAPoint() [1/5]

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.

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

References dqmiodumpmetadata::n.

Referenced by getDetIdsCloseToAPoint().

◆ getDetIdsCloseToAPoint() [2/5]

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

virtual

Definition at line 358 of file DetIdAssociator.cc.

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

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

◆ getDetIdsCloseToAPoint() [3/5]

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.

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

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

◆ getDetIdsCloseToAPoint() [4/5]

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.

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

References ztail::d, getDetIdsCloseToAPoint(), and point.

◆ getDetIdsCloseToAPoint() [5/5]

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.

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

References funct::abs(), etaBinSize_, getDetIdsCloseToAPoint(), dqm-mbProfile::log, LogTrace, M_PI, maxEta, CombinatorialSeedGeneratorForBeamHalo_cfi::maxTheta, EgHLTOffEleSelection_cfi::minEta, minTheta_, dqmiodumpmetadata::n, nPhi_, point, and funct::tan().

◆ getDetIdsInACone()

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.

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

References check_setup(), HGC3DClusterGenMatchSelector_cfi::dR, nearElement(), and selectAllInACone().

◆ getGeomDet()

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

pure virtual

Implemented in MuonDetIdAssociator, and CaloDetIdAssociator.

Referenced by MuonResidualsFromTrack::MuonResidualsFromTrack().

◆ getNumberOfSubdetectors()

virtual const unsigned int DetIdAssociator::getNumberOfSubdetectors ( ) const

inlineprotectedvirtual

Reimplemented in HcalDetIdAssociator, and EcalDetIdAssociator.

Definition at line 125 of file DetIdAssociator.h.

125 { return 1; }

Referenced by buildMap().

◆ getPosition()

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

protectedpure virtual

Implemented in MuonDetIdAssociator, and CaloDetIdAssociator.

Referenced by buildMap(), and nearElement().

◆ getValidDetIds()

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

protectedpure virtual

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

Referenced by buildMap().

◆ iEta()

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

virtual

look-up map eta index

Definition at line 132 of file DetIdAssociator.cc.

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

References etaBinSize_, createfilelist::int, nEta_, and point.

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

◆ index()

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

inlineprotected

Definition at line 139 of file DetIdAssociator.h.

139 { return iEta * nPhi_ + iPhi; }

References iEta(), iPhi(), and nPhi_.

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

◆ insideElement()

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

protectedpure virtual

Implemented in MuonDetIdAssociator, and CaloDetIdAssociator.

◆ iPhi()

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

virtual

look-up map phi index

Definition at line 134 of file DetIdAssociator.cc.

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

References createfilelist::int, M_PI, nPhi_, and point.

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

◆ name()

virtual const char* DetIdAssociator::name ( ) const

pure virtual

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

Referenced by config.CFG::__str__(), buildMap(), and validation.Sample::digest().

◆ nearElement()

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

protectedvirtual

Definition at line 363 of file DetIdAssociator.cc.

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

References SiPixelRawToDigiRegional_cfi::deltaPhi, HLT_2018_cff::distance, PV3DBase< T, PVType, FrameType >::eta(), getPosition(), M_PI, PV3DBase< T, PVType, FrameType >::phi(), and point.

Referenced by getDetIdsInACone().

◆ nEtaBins()

int DetIdAssociator::nEtaBins ( ) const

inline

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

Definition at line 106 of file DetIdAssociator.h.

106 { return nEta_; }

References nEta_.

◆ nPhiBins()

int DetIdAssociator::nPhiBins ( ) const

inline

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

Definition at line 104 of file DetIdAssociator.h.

104 { return nPhi_; }

References nPhi_.

◆ selectAllInACone()

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

inlinevirtual

helper to see if getDetIdsInACone is useful

Definition at line 82 of file DetIdAssociator.h.

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

References HGC3DClusterGenMatchSelector_cfi::dR, M_PI, and maxEta_.

Referenced by getDetIdsInACone().

◆ volume()

const FiducialVolume & DetIdAssociator::volume ( void ) const

get active detector volume

Definition at line 352 of file DetIdAssociator.cc.

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

References Exception, theMapIsValid_, and volume_.

Referenced by TrackExtrapolator::produce().

Member Data Documentation

◆ container_

std::vector<DetId> DetIdAssociator::container_

protected

Definition at line 149 of file DetIdAssociator.h.

Referenced by buildMap(), and fillSet().

◆ etaBinSize_

const double DetIdAssociator::etaBinSize_

protected

Definition at line 151 of file DetIdAssociator.h.

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

◆ lookupMap_

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

protected

Definition at line 148 of file DetIdAssociator.h.

Referenced by buildMap(), and fillSet().

◆ maxEta_

double DetIdAssociator::maxEta_

protected

Definition at line 152 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and selectAllInACone().

◆ minTheta_

double DetIdAssociator::minTheta_

protected

Definition at line 153 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and getDetIdsCloseToAPoint().

◆ nEta_

const int DetIdAssociator::nEta_

protected

Definition at line 144 of file DetIdAssociator.h.

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

◆ nPhi_

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().

◆ theMapIsValid_

bool DetIdAssociator::theMapIsValid_

protected

Definition at line 150 of file DetIdAssociator.h.

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

◆ 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

◆ const_iterator

Member Enumeration Documentation

◆ PropagationTarget

Constructor & Destructor Documentation

◆ DetIdAssociator()

◆ ~DetIdAssociator()

Member Function Documentation

◆ buildMap()

◆ check_setup()

◆ crossedElement()

◆ dumpMapContent() [1/2]

◆ dumpMapContent() [2/2]

◆ etaBinSize()

◆ fillSet()

◆ getCrossedDetIds() [1/2]

◆ getCrossedDetIds() [2/2]

◆ getDetIdPoints()

◆ getDetIdsCloseToAPoint() [1/5]

◆ getDetIdsCloseToAPoint() [2/5]

◆ getDetIdsCloseToAPoint() [3/5]

◆ getDetIdsCloseToAPoint() [4/5]

◆ getDetIdsCloseToAPoint() [5/5]

◆ getDetIdsInACone()

◆ getGeomDet()

◆ getNumberOfSubdetectors()

◆ getPosition()

◆ getValidDetIds()

◆ iEta()

◆ index()

◆ insideElement()

◆ iPhi()

◆ name()

◆ nearElement()

◆ nEtaBins()

◆ nPhiBins()

◆ selectAllInACone()

◆ volume()

Member Data Documentation

◆ container_

◆ etaBinSize_

◆ lookupMap_

◆ maxEta_

◆ minTheta_

◆ nEta_

◆ nPhi_

◆ theMapIsValid_

◆ volume_