RecHitsSortedInPhi Class Reference

#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 (float phiMin, float phiMax) const
void	hits (float phiMin, float phiMax, std::vector< Hit > &result) const
std::vector< Hit >	hits () 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

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

Definition at line 41 of file RecHitsSortedInPhi.h.

typedef BaseTrackerRecHit const* RecHitsSortedInPhi::Hit

Definition at line 19 of file RecHitsSortedInPhi.h.

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

Definition at line 38 of file RecHitsSortedInPhi.h.

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

Definition at line 39 of file RecHitsSortedInPhi.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file RecHitsSortedInPhi.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), runTheMatrix::const, runTauDisplay::dr, drphi, du, dv, PVValHelper::dz, h, mps_fire::i, isBarrel, lphi, Basic3DVector< T >::perp(), theHits, u, v, x, y, and z.

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

Member Function Documentation

Range RecHitsSortedInPhi::all ( ) const

inline

Definition at line 84 of file RecHitsSortedInPhi.h.

References theHits.

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

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

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

inlinestatic

Definition at line 112 of file RecHitsSortedInPhi.h.

References mps_fire::i.

Referenced by hits().

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

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

Definition at line 52 of file RecHitsSortedInPhi.cc.

References Geom::fpi(), Geom::ftwoPi(), Geom::pi(), diffTwoXMLs::r1, diffTwoXMLs::r2, theHits, and unsafeRange().

Referenced by HitPairGeneratorFromLayerPair::doublets().

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

bool RecHitsSortedInPhi::empty ( void  ) const

inline

Definition at line 45 of file RecHitsSortedInPhi.h.

References theHits.

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

{ return theHits.empty(); }

GlobalPoint RecHitsSortedInPhi::gp ( int  i ) const

inline

Definition at line 90 of file RecHitsSortedInPhi.h.

References x, y, and z.

Referenced by gv().

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

float RecHitsSortedInPhi::gv ( int  i ) const

inline

Definition at line 88 of file RecHitsSortedInPhi.h.

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

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

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

Definition at line 93 of file RecHitsSortedInPhi.cc.

References hits(), and mps_fire::result.

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

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

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

Definition at line 76 of file RecHitsSortedInPhi.cc.

References copyResult(), Geom::fpi(), Geom::ftwoPi(), Geom::pi(), and unsafeRange().

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

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

inline

Definition at line 76 of file RecHitsSortedInPhi.h.

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

Referenced by hits().

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

float RecHitsSortedInPhi::phi ( int  i ) const

inline

Definition at line 87 of file RecHitsSortedInPhi.h.

References mps_fire::i, and theHits.

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

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

float RecHitsSortedInPhi::rv ( int  i ) const

inline

Definition at line 89 of file RecHitsSortedInPhi.h.

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

Referenced by HitPairGeneratorFromLayerPair::doublets().

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

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

inline

Definition at line 46 of file RecHitsSortedInPhi.h.

References theHits.

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

{ return theHits.size(); }

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

Definition at line 99 of file RecHitsSortedInPhi.cc.

References cuda_std::lower_bound(), theHits, and cuda_std::upper_bound().

Referenced by doubleRange(), and hits().

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

Member Data Documentation

std::vector<float> RecHitsSortedInPhi::drphi

Definition at line 103 of file RecHitsSortedInPhi.h.

Referenced by HitPairGeneratorFromLayerPair::doublets(), and RecHitsSortedInPhi().

std::vector<float> RecHitsSortedInPhi::du

Definition at line 108 of file RecHitsSortedInPhi.h.

Referenced by HitPairGeneratorFromLayerPair::doublets(), and RecHitsSortedInPhi().

std::vector<float> RecHitsSortedInPhi::dv

Definition at line 109 of file RecHitsSortedInPhi.h.

Referenced by HitPairGeneratorFromLayerPair::doublets(), and RecHitsSortedInPhi().

bool RecHitsSortedInPhi::isBarrel

Definition at line 98 of file RecHitsSortedInPhi.h.

Referenced by HitPairGeneratorFromLayerPair::doublets(), gv(), MultiHitGeneratorFromChi2::hitSets(), RecHitsSortedInPhi(), and rv().

DetLayer const* RecHitsSortedInPhi::layer

Definition at line 97 of file RecHitsSortedInPhi.h.

Referenced by geometryXMLparser.DTAlignable::index(), and geometryXMLparser.CSCAlignable::index().

std::vector<float> RecHitsSortedInPhi::lphi

Definition at line 110 of file RecHitsSortedInPhi.h.

Referenced by RecHitsSortedInPhi().

std::vector<HitWithPhi> RecHitsSortedInPhi::theHits

Definition at line 95 of file RecHitsSortedInPhi.h.

Referenced by all(), doubleRange(), HitPairGeneratorFromLayerPair::doublets(), empty(), hits(), phi(), RecHitsSortedInPhi(), size(), and unsafeRange().

GlobalPoint RecHitsSortedInPhi::theOrigin

Definition at line 93 of file RecHitsSortedInPhi.h.

std::vector<float> RecHitsSortedInPhi::u

Definition at line 106 of file RecHitsSortedInPhi.h.

Referenced by RecHitsSortedInPhi(), and rv().

std::vector<float> RecHitsSortedInPhi::v

Definition at line 107 of file RecHitsSortedInPhi.h.

Referenced by RecHitsSortedInPhi(), and rv().

std::vector<float> RecHitsSortedInPhi::x

Definition at line 100 of file RecHitsSortedInPhi.h.

Referenced by svgfig.Curve.Sample::__repr__(), svgfig.Ellipse::__repr__(), HitPairGeneratorFromLayerPair::doublets(), gp(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r(), ntupleDataFormat._HitObject::r3D(), and RecHitsSortedInPhi().

std::vector<float> RecHitsSortedInPhi::y

Definition at line 101 of file RecHitsSortedInPhi.h.

Referenced by svgfig.Ellipse::__repr__(), HitPairGeneratorFromLayerPair::doublets(), gp(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r(), ntupleDataFormat._HitObject::r3D(), and RecHitsSortedInPhi().

std::vector<float> RecHitsSortedInPhi::z

Definition at line 102 of file RecHitsSortedInPhi.h.

Referenced by HitPairGeneratorFromLayerPair::doublets(), gp(), gv(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r3D(), and RecHitsSortedInPhi().