#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 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.

48 { Barrel, ForwardEndcap, BackwardEndcap };

DetIdAssociator::Barrel

Definition: DetIdAssociator.h:48

DetIdAssociator::BackwardEndcap

Definition: DetIdAssociator.h:48

DetIdAssociator::ForwardEndcap

Definition: DetIdAssociator.h:48

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.

58 {}

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_25ns14e33_v1_cff::etaMax, HLT_25ns14e33_v1_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 134 of file DetIdAssociator.h.

Referenced by getCrossedDetIds().

138 { 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 109 of file DetIdAssociator.h.

References etaBinSize_.

109 { return etaBinSize_;};

DetIdAssociator::etaBinSize_

const double etaBinSize_

Definition: DetIdAssociator.h:157

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(), cmsBatch::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(), M_PI, maxEta_, and nearElement().

 {
    if ( dR > 2*M_PI && dR > maxEta_ ) 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 129 of file DetIdAssociator.h.

Referenced by buildMap().

129 { 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 143 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 107 of file DetIdAssociator.h.

References nEta_.

107 { return nEta_;}

DetIdAssociator::nEta_

const int nEta_

Definition: DetIdAssociator.h:150

int DetIdAssociator::nPhiBins ( ) const

inline

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

Definition at line 105 of file DetIdAssociator.h.

References nPhi_.

105 { return nPhi_;}

DetIdAssociator::nPhi_

const int nPhi_

Definition: DetIdAssociator.h:149

virtual void DetIdAssociator::setConditions ( const DetIdAssociatorRecord & )

inlinevirtual

Reimplemented in MuonDetIdAssociator.

Definition at line 118 of file DetIdAssociator.h.

118 {};

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 155 of file DetIdAssociator.h.

Referenced by buildMap(), and fillSet().

const double DetIdAssociator::etaBinSize_

protected

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

Referenced by buildMap(), and fillSet().

double DetIdAssociator::maxEta_

protected

Definition at line 158 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and getDetIdsInACone().

double DetIdAssociator::minTheta_

protected

Definition at line 159 of file DetIdAssociator.h.

Referenced by DetIdAssociator(), and getDetIdsCloseToAPoint().

const int DetIdAssociator::nEta_

protected

Definition at line 150 of file DetIdAssociator.h.

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

const int DetIdAssociator::nPhi_

protected

Definition at line 149 of file DetIdAssociator.h.

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

bool DetIdAssociator::theMapIsValid_

protected

Definition at line 156 of file DetIdAssociator.h.

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

FiducialVolume DetIdAssociator::volume_

protected

Definition at line 164 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