#include <RecHitsSortedInPhi.h>

Classes
struct	HitLessPhi

class	HitWithPhi

Public Types
using	DoubleRange = std::array< int, 4 >

typedef BaseTrackerRecHit const *	Hit

typedef std::vector< HitWithPhi >::const_iterator	HitIter

typedef std::pair< HitIter, HitIter >	Range

Public Member Functions
Range	all () const

DoubleRange	doubleRange (float phiMin, float phiMax) const

bool	empty () const

GlobalPoint	gp (int i) const

float	gv (int i) const

std::vector< Hit >	hits () const

std::vector< Hit >	hits (float phiMin, float phiMax) const

void	hits (float phiMin, float phiMax, std::vector< Hit > &result) const

float	phi (int i) const

	RecHitsSortedInPhi (const std::vector< Hit > &hits, GlobalPoint const &origin, DetLayer const *il)

float	rv (int i) const

std::size_t	size () const

Range	unsafeRange (float phiMin, float phiMax) const

Static Public Member Functions
static void	copyResult (const Range &range, std::vector< Hit > &result)

Public Attributes
std::vector< float >	drphi

std::vector< float >	du

std::vector< float >	dv

bool	isBarrel

DetLayer const *	layer

std::vector< float >	lphi

std::vector< HitWithPhi >	theHits

GlobalPoint	theOrigin

std::vector< float >	u

std::vector< float >	v

std::vector< float >	x

std::vector< float >	y

std::vector< float >	z

Detailed Description

A RecHit container sorted in phi. Provides fast access for hits in a given phi window using binary search.

Definition at line 17 of file RecHitsSortedInPhi.h.

Member Typedef Documentation

◆ DoubleRange

using RecHitsSortedInPhi::DoubleRange = std::array<int, 4>

Definition at line 41 of file RecHitsSortedInPhi.h.

◆ Hit

typedef BaseTrackerRecHit const* RecHitsSortedInPhi::Hit

Definition at line 19 of file RecHitsSortedInPhi.h.

◆ HitIter

typedef std::vector<HitWithPhi>::const_iterator RecHitsSortedInPhi::HitIter

Definition at line 38 of file RecHitsSortedInPhi.h.

◆ Range

typedef std::pair<HitIter, HitIter> RecHitsSortedInPhi::Range

Definition at line 39 of file RecHitsSortedInPhi.h.

Constructor & Destructor Documentation

◆ RecHitsSortedInPhi()

RecHitsSortedInPhi::RecHitsSortedInPhi	(	const std::vector< Hit > &	hits,
		GlobalPoint const &	origin,
		DetLayer const *	il
	)

Definition at line 7 of file RecHitsSortedInPhi.cc.

     : layer(il),
       isBarrel(il->isBarrel()),
       x(hits.size()),
       y(hits.size()),
       z(hits.size()),
       drphi(hits.size()),
       u(hits.size()),
       v(hits.size()),
       du(hits.size()),
       dv(hits.size()),
       lphi(hits.size()) {
   // standard region have origin as 0,0,z (not true!!!!0
   // cosmic region never used here
   // assert(origin.x()==0 && origin.y()==0);
  
   theHits.reserve(hits.size());
   for (auto const& hp : hits)
     theHits.emplace_back(hp);
  
   std::sort(theHits.begin(), theHits.end(), HitLessPhi());
  
   for (unsigned int i = 0; i != theHits.size(); ++i) {
     auto const& h = *theHits[i].hit();
     auto const& gs = static_cast<BaseTrackerRecHit const&>(h).globalState();
     auto loc = gs.position - origin.basicVector();
     float lr = loc.perp();
     // float lr = gs.position.perp();
     float lz = gs.position.z();
     float dr = gs.errorR;
     float dz = gs.errorZ;
     // r[i] = gs.position.perp();
     // phi[i] = gs.position.barePhi();
     x[i] = gs.position.x();
     y[i] = gs.position.y();
     z[i] = lz;
     drphi[i] = gs.errorRPhi;
     u[i] = isBarrel ? lr : lz;
     v[i] = isBarrel ? lz : lr;
     du[i] = isBarrel ? dr : dz;
     dv[i] = isBarrel ? dz : dr;
     lphi[i] = loc.barePhi();
   }
 }

References PV3DBase< T, PVType, FrameType >::basicVector(), flavorHistoryFilter_cfi::dr, drphi, du, dv, PVValHelper::dz, hits(), trackingPlots::hp, mps_fire::i, isBarrel, lphi, Basic3DVector< T >::perp(), jetUpdater_cfi::sort, theHits, u, v, x, y, and z.

Member Function Documentation

◆ all()

Range RecHitsSortedInPhi::all ( ) const

inline

Definition at line 84 of file RecHitsSortedInPhi.h.

84 { return Range(theHits.begin(), theHits.end()); }

References theHits.

Referenced by HitQuadrupletGeneratorFromLayerPairForPhotonConversion::hitPairs(), MultiHitGeneratorFromChi2::hitSets(), and LayerHitMapCache::operator()().

◆ copyResult()

static void RecHitsSortedInPhi::copyResult	(	const Range &	range,
		std::vector< Hit > &	result
	)

inlinestatic

Definition at line 112 of file RecHitsSortedInPhi.h.

                                                                      {
     result.reserve(result.size() + (range.second - range.first));
     for (HitIter i = range.first; i != range.second; i++)
       result.push_back(i->hit());
   }

References mps_fire::i, FastTimerService_cff::range, and mps_fire::result.

Referenced by hits().

◆ doubleRange()

RecHitsSortedInPhi::DoubleRange RecHitsSortedInPhi::doubleRange	(	float	phiMin,
		float	phiMax
	)		const

Definition at line 52 of file RecHitsSortedInPhi.cc.

                                                                                               {
   Range r1, r2;
   if (phiMin < phiMax) {
     if (phiMin < -Geom::fpi()) {
       r1 = unsafeRange(phiMin + Geom::ftwoPi(), Geom::fpi());
       r2 = unsafeRange(-Geom::fpi(), phiMax);
     } else if (phiMax > Geom::pi()) {
       r1 = unsafeRange(phiMin, Geom::fpi());
       r2 = unsafeRange(-Geom::fpi(), phiMax - Geom::ftwoPi());
     } else {
       r1 = unsafeRange(phiMin, phiMax);
       r2 = Range(theHits.begin(), theHits.begin());
     }
   } else {
     r1 = unsafeRange(phiMin, Geom::fpi());
     r2 = unsafeRange(-Geom::fpi(), phiMax);
   }
  
   return (DoubleRange){{int(r1.first - theHits.begin()),
                         int(r1.second - theHits.begin()),
                         int(r2.first - theHits.begin()),
                         int(r2.second - theHits.begin())}};
 }

References Geom::fpi(), Geom::ftwoPi(), createfilelist::int, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, Geom::pi(), diffTwoXMLs::r1, diffTwoXMLs::r2, theHits, and unsafeRange().

Referenced by HitPairGeneratorFromLayerPair::doublets().

◆ empty()

bool RecHitsSortedInPhi::empty ( void ) const

inline

Definition at line 45 of file RecHitsSortedInPhi.h.

45 { return theHits.empty(); }

References theHits.

Referenced by HitPairGeneratorFromLayerPair::doublets(), HitQuadrupletGeneratorFromLayerPairForPhotonConversion::hitPairs(), and HitPairGeneratorFromLayerPairForPhotonConversion::hitPairs().

◆ gp()

GlobalPoint RecHitsSortedInPhi::gp ( int i ) const

inline

Definition at line 90 of file RecHitsSortedInPhi.h.

90 { return GlobalPoint(x[i], y[i], z[i]); }

References mps_fire::i, x, y, and z.

Referenced by gv().

◆ gv()

float RecHitsSortedInPhi::gv ( int i ) const

inline

Definition at line 88 of file RecHitsSortedInPhi.h.

88 { return isBarrel ? z[i] : gp(i).perp(); } // global v

References gp(), mps_fire::i, isBarrel, PV3DBase< T, PVType, FrameType >::perp(), and z.

◆ hits() [1/3]

std::vector<Hit> RecHitsSortedInPhi::hits ( void ) const

inline

Definition at line 76 of file RecHitsSortedInPhi.h.

                               {
     std::vector<Hit> result;
     result.reserve(theHits.size());
     for (HitIter i = theHits.begin(); i != theHits.end(); i++)
       result.push_back(i->hit());
     return result;
   }

References mps_fire::i, mps_fire::result, and theHits.

Referenced by hits(), and RecHitsSortedInPhi().

◆ hits() [2/3]

std::vector< RecHitsSortedInPhi::Hit > RecHitsSortedInPhi::hits	(	float	phiMin,
		float	phiMax
	)		const

Definition at line 93 of file RecHitsSortedInPhi.cc.

                                                                                           {
   std::vector<Hit> result;
   hits(phiMin, phiMax, result);
   return result;
 }

References hits(), AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, and mps_fire::result.

Referenced by HitQuadrupletGeneratorFromLayerPairForPhotonConversion::hitPairs(), HitPairGeneratorFromLayerPairForPhotonConversion::hitPairs(), PixelTripletNoTipGenerator::hitTriplets(), and PixelTripletLowPtGenerator::hitTriplets().

◆ hits() [3/3]

void RecHitsSortedInPhi::hits	(	float	phiMin,
		float	phiMax,
		std::vector< Hit > &	result
	)		const

Definition at line 76 of file RecHitsSortedInPhi.cc.

                                                                                       {
   if (phiMin < phiMax) {
     if (phiMin < -Geom::fpi()) {
       copyResult(unsafeRange(phiMin + Geom::ftwoPi(), Geom::fpi()), result);
       copyResult(unsafeRange(-Geom::fpi(), phiMax), result);
     } else if (phiMax > Geom::pi()) {
       copyResult(unsafeRange(phiMin, Geom::fpi()), result);
       copyResult(unsafeRange(-Geom::fpi(), phiMax - Geom::ftwoPi()), result);
     } else {
       copyResult(unsafeRange(phiMin, phiMax), result);
     }
   } else {
     copyResult(unsafeRange(phiMin, Geom::fpi()), result);
     copyResult(unsafeRange(-Geom::fpi(), phiMax), result);
   }
 }

References copyResult(), Geom::fpi(), Geom::ftwoPi(), AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, Geom::pi(), mps_fire::result, and unsafeRange().

◆ phi()

float RecHitsSortedInPhi::phi ( int i ) const

inline

Definition at line 87 of file RecHitsSortedInPhi.h.

87 { return theHits[i].phi(); }

References mps_fire::i, and theHits.

Referenced by Particle.Particle::__str__(), and ntupleDataFormat.Track::phiPull().

◆ rv()

float RecHitsSortedInPhi::rv ( int i ) const

inline

Definition at line 89 of file RecHitsSortedInPhi.h.

89 { return isBarrel ? u[i] : v[i]; } // dispaced r

References mps_fire::i, isBarrel, u, and v.

Referenced by HitPairGeneratorFromLayerPair::doublets().

◆ size()

std::size_t RecHitsSortedInPhi::size ( void ) const

inline

Definition at line 46 of file RecHitsSortedInPhi.h.

46 { return theHits.size(); }

References theHits.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), and HitPairGeneratorFromLayerPair::doublets().

◆ unsafeRange()

RecHitsSortedInPhi::Range RecHitsSortedInPhi::unsafeRange	(	float	phiMin,
		float	phiMax
	)		const

Definition at line 99 of file RecHitsSortedInPhi.cc.

                                                                                         {
   auto low = std::lower_bound(theHits.begin(), theHits.end(), HitWithPhi(phiMin), HitLessPhi());
   return Range(low, std::upper_bound(low, theHits.end(), HitWithPhi(phiMax), HitLessPhi()));
 }