CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
ForwardMeasurementEstimator Class Reference

#include <ForwardMeasurementEstimator.h>

Inheritance diagram for ForwardMeasurementEstimator:
MeasurementEstimator

Public Member Functions

ForwardMeasurementEstimatorclone () const override
 
std::pair< bool, double > estimate (const TrajectoryStateOnSurface &ts, const TrackingRecHit &hit) const override
 
virtual std::pair< bool, double > estimate (const TrajectoryStateOnSurface &ts, const GlobalPoint &gp) const
 
virtual std::pair< bool, double > estimate (const GlobalPoint &vprim, const TrajectoryStateOnSurface &ts, const GlobalPoint &gp) const
 
bool estimate (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const override
 
 ForwardMeasurementEstimator ()
 
 ForwardMeasurementEstimator (float phiMin, float phiMax, float rMin, float rMax)
 
MeasurementEstimator::Local2DVector maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const override
 
void setPhiRange (float dummyphiMin, float dummyphiMax)
 
void setRRange (float rmin, float rmax)
 
void setRRangeI (float rmin, float rmax)
 
- Public Member Functions inherited from MeasurementEstimator
virtual SurfaceReturnType estimate (const TrajectoryStateOnSurface &ts, const Plane &plane) const =0
 
virtual Local2DVector maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const Plane &plane) const =0
 
float maxSagitta () const
 
 MeasurementEstimator ()
 
 MeasurementEstimator (float maxSag, float minToll, float mpt)
 
float minPt2ForHitRecoveryInGluedDet () const
 
float minTolerance2 () const
 
virtual bool preFilter (const TrajectoryStateOnSurface &, OpaquePayload const &) const
 
virtual ~MeasurementEstimator ()
 

Private Attributes

float thePhiMax
 
float thePhiMin
 
float theRMax
 
float theRMaxI
 
float theRMin
 
float theRMinI
 

Additional Inherited Members

- Public Types inherited from MeasurementEstimator
using HitReturnType = std::pair< bool, double >
 
using Local2DVector = Vector2DBase< float, LocalTag >
 
using SurfaceReturnType = bool
 

Detailed Description

Description: Class defining the search area in the forward disks in the pixel match, ported from ORCA

Implementation: <Notes on="" implementation>="">

Description: MeasurementEstimator for Pixel Endcap, ported from ORCA

Implementation: <Notes on="" implementation>="">

Definition at line 26 of file ForwardMeasurementEstimator.h.

Constructor & Destructor Documentation

ForwardMeasurementEstimator::ForwardMeasurementEstimator ( )
inline

Definition at line 31 of file ForwardMeasurementEstimator.h.

Referenced by clone().

32  {}
ForwardMeasurementEstimator::ForwardMeasurementEstimator ( float  phiMin,
float  phiMax,
float  rMin,
float  rMax 
)
inline

Member Function Documentation

ForwardMeasurementEstimator* ForwardMeasurementEstimator::clone ( void  ) const
inlineoverridevirtual
std::pair< bool, double > ForwardMeasurementEstimator::estimate ( const TrajectoryStateOnSurface ts,
const TrackingRecHit hit 
) const
overridevirtual

Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible. The TrajectoryStateOnSurface must be on the same Surface as the RecHit. For an estimator where there is no value computed, e.g. fixed window estimator, only the first(bool) part is of interest.

Implements MeasurementEstimator.

Definition at line 26 of file ForwardMeasurementEstimator.cc.

References TrackingRecHit::det(), runTauDisplay::gp, TrackingRecHit::localPosition(), GeomDet::surface(), and Surface::toGlobal().

Referenced by estimate(), and setRRangeI().

27  {
28  LocalPoint lp = hit.localPosition();
29  GlobalPoint gp = hit.det()->surface().toGlobal( lp);
30  return estimate(ts,gp);
31 }
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:106
std::pair< bool, double > estimate(const TrajectoryStateOnSurface &ts, const TrackingRecHit &hit) const override
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:42
const GeomDet * det() const
virtual LocalPoint localPosition() const =0
std::pair< bool, double > ForwardMeasurementEstimator::estimate ( const TrajectoryStateOnSurface ts,
const GlobalPoint gp 
) const
virtual

Definition at line 33 of file ForwardMeasurementEstimator.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::barePhi(), estimate(), f, TrajectoryStateOnSurface::globalParameters(), PV3DBase< T, PVType, FrameType >::perp(), pi, GlobalTrajectoryParameters::position(), hfnoseDigiStudy_cfi::rMax, hfnoseDigiStudy_cfi::rMin, thePhiMax, thePhiMin, theRMax, theRMaxI, theRMin, theRMinI, and PV3DBase< T, PVType, FrameType >::z().

33  {
34 
35  float tsR = ts.globalParameters().position().perp();
36  float rhR = gp.perp();
37  float rDiff = tsR - rhR;
38  float rMin = theRMin;
39  float rMax = theRMax;
40  float myZ = gp.z();
41  if( (std::abs(myZ)> 70.f) & (std::abs(myZ)<170.f)) {
42  rMin = theRMinI;
43  rMax = theRMaxI;
44  }
45  if( rDiff >= rMax || rDiff <= rMin ) return std::pair<bool,double>(false,0.);
46 
47  float tsPhi = ts.globalParameters().position().barePhi();
48  float rhPhi = gp.barePhi();
49 
50  float myPhimin = thePhiMin;
51  float myPhimax = thePhiMax;
52 
53  float phiDiff = tsPhi - rhPhi;
54  if (phiDiff > pi) phiDiff -= twopi;
55  if (phiDiff < -pi) phiDiff += twopi;
56 
57  if ( (phiDiff < myPhimax) & (phiDiff > myPhimin) ) {
58  return std::pair<bool,double>(true,1.);
59  } else {
60  return std::pair<bool,double>(false,0.);
61  }
62 }
T perp() const
Definition: PV3DBase.h:72
const Double_t pi
T barePhi() const
Definition: PV3DBase.h:68
T z() const
Definition: PV3DBase.h:64
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
double f[11][100]
const GlobalTrajectoryParameters & globalParameters() const
std::pair< bool, double > ForwardMeasurementEstimator::estimate ( const GlobalPoint vprim,
const TrajectoryStateOnSurface ts,
const GlobalPoint gp 
) const
virtual

Definition at line 65 of file ForwardMeasurementEstimator.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::barePhi(), estimate(), f, TrajectoryStateOnSurface::globalParameters(), runTauDisplay::gp, normalized_phi(), PV3DBase< T, PVType, FrameType >::perp(), GlobalTrajectoryParameters::position(), hfnoseDigiStudy_cfi::rMax, hfnoseDigiStudy_cfi::rMin, thePhiMax, thePhiMin, theRMax, theRMaxI, theRMin, theRMinI, and PV3DBase< T, PVType, FrameType >::z().

68  {
69  GlobalVector ts = absolute_ts.globalParameters().position() - vprim ;
70  GlobalVector gp = absolute_gp - vprim ;
71 
72  float rhR = gp.perp();
73  float tsR = ts.perp();
74  float rDiff = rhR - tsR;
75  float rMin = theRMin;
76  float rMax = theRMax;
77  float myZ = gp.z();
78  if( (std::abs(myZ)> 70.f) & (std::abs(myZ)<170.f) ) {
79  rMin = theRMinI;
80  rMax = theRMaxI;
81  }
82 
83  if( (rDiff >= rMax) | (rDiff <= rMin) ) return std::pair<bool,double>(false,0.);
84 
85  float tsPhi = ts.barePhi();
86  float rhPhi = gp.barePhi();
87 
88  float myPhimin = thePhiMin;
89  float myPhimax = thePhiMax;
90 
91  float phiDiff = normalized_phi(rhPhi - tsPhi) ;
92 
93  if ( phiDiff < myPhimax && phiDiff > myPhimin )
94  { return std::pair<bool,double>(true,1.) ; }
95  else
96  { return std::pair<bool,double>(false,0.) ; }
97 }
T perp() const
Definition: PV3DBase.h:72
T barePhi() const
Definition: PV3DBase.h:68
RealType normalized_phi(RealType phi)
T z() const
Definition: PV3DBase.h:64
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
double f[11][100]
const GlobalTrajectoryParameters & globalParameters() const
bool ForwardMeasurementEstimator::estimate ( const TrajectoryStateOnSurface ts,
const BoundPlane plane 
) const
override

Definition at line 100 of file ForwardMeasurementEstimator.cc.

References funct::abs(), funct::false, TrajectoryStateOnSurface::globalParameters(), maximalLocalDisplacement(), Geom::phiLess(), GlobalTrajectoryParameters::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, rangesIntersect(), thePhiMax, thePhiMin, funct::true, x, and y.

102  {
103  typedef std::pair<float,float> Range ;
104 
105  GlobalPoint trajPos(ts.globalParameters().position());
106 
107  float r1 = 0.;
108  float r2 = 40.;
109 
110  Range trajRRange(trajPos.perp() - r1, trajPos.perp() + r2);
111  Range trajPhiRange(trajPos.phi() - std::abs(thePhiMin), trajPos.phi() + std::abs(thePhiMax));
112 
113  if(rangesIntersect(trajRRange, plane.rSpan()) &&
114  rangesIntersect(trajPhiRange, plane.phiSpan(), [](auto x,auto y){ return Geom::phiLess(x, y);}))
115  { return true ; }
116  else
117  { return false ; }
118 }
PixelRecoRange< float > Range
bool rangesIntersect(const Range &a, const Range &b)
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
bool phiLess(float phi1, float phi2)
Definition: VectorUtil.h:23
const GlobalTrajectoryParameters & globalParameters() const
MeasurementEstimator::Local2DVector ForwardMeasurementEstimator::maximalLocalDisplacement ( const TrajectoryStateOnSurface ts,
const BoundPlane plane 
) const
override

Definition at line 122 of file ForwardMeasurementEstimator.cc.

References TrajectoryStateOnSurface::hasError(), TrajectoryStateOnSurface::localError(), d0_phi_analyzer_cff::nSigmaCut, LocalTrajectoryError::positionError(), mathSSE::sqrt(), LocalError::xx(), and LocalError::yy().

Referenced by clone(), and estimate().

124  {
125  float nSigmaCut = 3.;
126  if ( ts.hasError())
127  {
128  LocalError le = ts.localError().positionError();
129  return Local2DVector( std::sqrt(le.xx())*nSigmaCut, std::sqrt(le.yy())*nSigmaCut);
130  }
131  else
132  return Local2DVector(999999,999999) ;
133  }
float xx() const
Definition: LocalError.h:24
LocalError positionError() const
float yy() const
Definition: LocalError.h:26
Vector2DBase< float, LocalTag > Local2DVector
T sqrt(T t)
Definition: SSEVec.h:18
const LocalTrajectoryError & localError() const
void ForwardMeasurementEstimator::setPhiRange ( float  dummyphiMin,
float  dummyphiMax 
)
inline

Definition at line 37 of file ForwardMeasurementEstimator.h.

References thePhiMax, and thePhiMin.

void ForwardMeasurementEstimator::setRRange ( float  rmin,
float  rmax 
)
inline

Definition at line 39 of file ForwardMeasurementEstimator.h.

References theRMax, and theRMin.

void ForwardMeasurementEstimator::setRRangeI ( float  rmin,
float  rmax 
)
inline

Member Data Documentation

float ForwardMeasurementEstimator::thePhiMax
private

Definition at line 59 of file ForwardMeasurementEstimator.h.

Referenced by estimate(), and setPhiRange().

float ForwardMeasurementEstimator::thePhiMin
private

Definition at line 58 of file ForwardMeasurementEstimator.h.

Referenced by estimate(), and setPhiRange().

float ForwardMeasurementEstimator::theRMax
private

Definition at line 61 of file ForwardMeasurementEstimator.h.

Referenced by estimate(), and setRRange().

float ForwardMeasurementEstimator::theRMaxI
private

Definition at line 63 of file ForwardMeasurementEstimator.h.

Referenced by estimate(), and setRRangeI().

float ForwardMeasurementEstimator::theRMin
private

Definition at line 60 of file ForwardMeasurementEstimator.h.

Referenced by estimate(), and setRRange().

float ForwardMeasurementEstimator::theRMinI
private

Definition at line 62 of file ForwardMeasurementEstimator.h.

Referenced by estimate(), and setRRangeI().