DetIdAssociator Class Reference

#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...
virtual void	setConditions (const DetIdAssociatorRecord &)
virtual void	setGeometry (const DetIdAssociatorRecord &)=0
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 *=0) 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 46 of file DetIdAssociator.h.

Member Typedef Documentation

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

Definition at line 55 of file DetIdAssociator.h.

Member Enumeration Documentation

enum DetIdAssociator::PropagationTarget

Enumerator
Barrel
ForwardEndcap
BackwardEndcap

Definition at line 48 of file DetIdAssociator.h.

{ Barrel, ForwardEndcap, BackwardEndcap };

Constructor & Destructor Documentation

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

Definition at line 25 of file DetIdAssociator.cc.

References etaBinSize_, Exception, create_public_lumi_plots::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 58 of file DetIdAssociator.h.

{}

Member Function Documentation

void DetIdAssociator::buildMap ( )

virtual

make the look-up map

Definition at line 142 of file DetIdAssociator.cc.

References funct::abs(), FiducialVolume::addActivePoint(), assert(), newFWLiteAna::bin, check_setup(), container_, FiducialVolume::determinInnerDimensions(), etaBinSize_, HLT_25ns10e33_v2_cff::etaMax, HLT_25ns10e33_v2_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_, relval_parameters_module::step, theMapIsValid_, 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,0 )
                          << "\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 377 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 *  = 0  ) const

inlineprotectedvirtual

Reimplemented in CaloDetIdAssociator.

Definition at line 138 of file DetIdAssociator.h.

Referenced by getCrossedDetIds().

                                                             { return false; }

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

protectedvirtual

Definition at line 324 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(), i, and 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 113 of file DetIdAssociator.h.

References etaBinSize_.

{ return etaBinSize_;};

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

protected

Definition at line 385 of file DetIdAssociator.cc.

References container_, 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 284 of file DetIdAssociator.cc.

References check_setup(), crossedElement(), 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 303 of file DetIdAssociator.cc.

References check_setup(), crossedElement(), 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 35 of file DetIdAssociator.cc.

References gen::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(), i, iEta(), iPhi(), 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 359 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 94 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(), dqm-mbProfile::log, LogTrace, M_PI, maxEta, benchmark_cfg::minEta, minTheta_, gen::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 268 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 EcalDetIdAssociator, and HcalDetIdAssociator.

Definition at line 133 of file DetIdAssociator.h.

Referenced by buildMap().

{ 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, EcalDetIdAssociator, HcalDetIdAssociator, HODetIdAssociator, and PreshowerDetIdAssociator.

Referenced by buildMap().

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

virtual

look-up map eta index

Definition at line 131 of file DetIdAssociator.cc.

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

Referenced by buildMap(), and getDetIdsCloseToAPoint().

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

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

inlineprotected

Definition at line 147 of file DetIdAssociator.h.

References iPhi(), and nPhi_.

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

                                                                  {
     return iEta*nPhi_+iPhi;
   }

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 136 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, EcalDetIdAssociator, HcalDetIdAssociator, HODetIdAssociator, and PreshowerDetIdAssociator.

Referenced by cuy.divideElement::__init__(), cuy.plotElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), cuy.graphElement::__init__(), config.CFG::__str__(), buildMap(), and validation.Sample::digest().

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

protectedvirtual

Definition at line 367 of file DetIdAssociator.cc.

References SiPixelRawToDigiRegional_cfi::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 111 of file DetIdAssociator.h.

References nEta_.

{ return nEta_;}

int DetIdAssociator::nPhiBins ( ) const

inline

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

Definition at line 109 of file DetIdAssociator.h.

References nPhi_.

{ return nPhi_;}

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

inlinevirtual

helper to see if getDetIdsInACone is useful

Definition at line 86 of file DetIdAssociator.h.

References M_PI, and maxEta_.

Referenced by getDetIdsInACone().

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

virtual void DetIdAssociator::setConditions ( const DetIdAssociatorRecord &  )

inlinevirtual

Reimplemented in MuonDetIdAssociator.

Definition at line 122 of file DetIdAssociator.h.

{};

virtual void DetIdAssociator::setGeometry ( const DetIdAssociatorRecord &  )

pure virtual

Implemented in CaloDetIdAssociator, and MuonDetIdAssociator.

const FiducialVolume & DetIdAssociator::volume

(

void

)

const

get active detector volume

Definition at line 353 of file DetIdAssociator.cc.

References Exception, theMapIsValid_, and volume_.

Referenced by BetaCalculatorECAL::calcEcalDeposit().

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

Member Data Documentation

std::vector<DetId> DetIdAssociator::container_

protected

Definition at line 159 of file DetIdAssociator.h.

Referenced by buildMap(), and fillSet().

const double DetIdAssociator::etaBinSize_

protected

Definition at line 161 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 158 of file DetIdAssociator.h.

Referenced by buildMap(), and fillSet().

double DetIdAssociator::maxEta_

protected

Definition at line 162 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and selectAllInACone().

double DetIdAssociator::minTheta_

protected

Definition at line 163 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and getDetIdsCloseToAPoint().

const int DetIdAssociator::nEta_

protected

Definition at line 154 of file DetIdAssociator.h.

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

const int DetIdAssociator::nPhi_

protected

Definition at line 153 of file DetIdAssociator.h.

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

bool DetIdAssociator::theMapIsValid_

protected

Definition at line 160 of file DetIdAssociator.h.

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

FiducialVolume DetIdAssociator::volume_

protected

Definition at line 168 of file DetIdAssociator.h.

Referenced by buildMap(), and volume().