ForwardMeasurementEstimator Class Reference

#include <ForwardMeasurementEstimator.h>

Inheritance diagram for ForwardMeasurementEstimator:

Public Member Functions
virtual ForwardMeasurementEstimator *	clone () const
virtual std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &ts, const TrackingRecHit &hit) const
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
virtual bool	estimate (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const
	ForwardMeasurementEstimator ()
	ForwardMeasurementEstimator (float phiMin, float phiMax, float rMin, float rMax)
MeasurementEstimator::Local2DVector	maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const
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().

     {}

ForwardMeasurementEstimator::ForwardMeasurementEstimator ( float  phiMin, float  phiMax, float  rMin, float  rMax  )

inline

Definition at line 33 of file ForwardMeasurementEstimator.h.

     : thePhiMin(phiMin), thePhiMax( phiMax), theRMin(rMin), theRMax(rMax)
     {}

Member Function Documentation

virtual ForwardMeasurementEstimator* ForwardMeasurementEstimator::clone ( void  ) const

inlinevirtual

Implements MeasurementEstimator.

Definition at line 50 of file ForwardMeasurementEstimator.h.

References ForwardMeasurementEstimator().

     { return new ForwardMeasurementEstimator(*this) ; }

std::pair< bool, double > ForwardMeasurementEstimator::estimate ( const TrajectoryStateOnSurface &  ts, const TrackingRecHit &  hit  ) const

virtual

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().

                                                                   {
  LocalPoint lp = hit.localPosition();
  GlobalPoint gp = hit.det()->surface().toGlobal( lp);
  return estimate(ts,gp);
}

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(), f, TrajectoryStateOnSurface::globalParameters(), PV3DBase< T, PVType, FrameType >::perp(), pi, GlobalTrajectoryParameters::position(), thePhiMax, thePhiMin, theRMax, theRMaxI, theRMin, theRMinI, and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                           {

  float tsR = ts.globalParameters().position().perp();
  float rhR = gp.perp();
  float rDiff = tsR - rhR;
  float rMin = theRMin;
  float rMax = theRMax;
  float myZ = gp.z();
  if( (std::abs(myZ)> 70.f)  &  (std::abs(myZ)<170.f)) {
    rMin = theRMinI;
    rMax = theRMaxI;
  }
  if( rDiff >= rMax || rDiff <= rMin ) return std::pair<bool,double>(false,0.);
  
  float tsPhi = ts.globalParameters().position().barePhi();
  float rhPhi = gp.barePhi();
  
  float myPhimin = thePhiMin;
  float myPhimax = thePhiMax;

  float phiDiff = tsPhi - rhPhi;
  if (phiDiff > pi) phiDiff -= twopi;
  if (phiDiff < -pi) phiDiff += twopi;

  if ( (phiDiff < myPhimax) & (phiDiff > myPhimin) ) {
    return std::pair<bool,double>(true,1.);
  } else {
    return std::pair<bool,double>(false,0.);
  }
}

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(), f, TrajectoryStateOnSurface::globalParameters(), runTauDisplay::gp, normalized_phi(), PV3DBase< T, PVType, FrameType >::perp(), GlobalTrajectoryParameters::position(), and PV3DBase< T, PVType, FrameType >::z().

 {
  GlobalVector ts = absolute_ts.globalParameters().position() - vprim ;
  GlobalVector gp = absolute_gp - vprim ;

  float rhR = gp.perp();
  float tsR = ts.perp();
  float rDiff = rhR - tsR;
  float rMin = theRMin;
  float rMax = theRMax;
  float myZ = gp.z();
  if( (std::abs(myZ)> 70.f) &  (std::abs(myZ)<170.f) ) {
    rMin = theRMinI;
    rMax = theRMaxI;
  }
  
  if( (rDiff >= rMax) | (rDiff <= rMin) ) return std::pair<bool,double>(false,0.);
  
  float tsPhi = ts.barePhi();
  float rhPhi = gp.barePhi();  
  
  float myPhimin = thePhiMin;
  float myPhimax = thePhiMax;  

  float phiDiff = normalized_phi(rhPhi - tsPhi) ;  

  if ( phiDiff < myPhimax && phiDiff > myPhimin )
   { return std::pair<bool,double>(true,1.) ; }
  else
   { return std::pair<bool,double>(false,0.) ; }
}

bool ForwardMeasurementEstimator::estimate ( const TrajectoryStateOnSurface &  ts, const BoundPlane &  plane  ) const

virtual

Definition at line 100 of file ForwardMeasurementEstimator.cc.

References funct::abs(), edm::false, TrajectoryStateOnSurface::globalParameters(), GlobalTrajectoryParameters::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, rangesIntersect(), and funct::true.

 {
  typedef std::pair<float,float> Range ;

  GlobalPoint trajPos(ts.globalParameters().position());

  float r1 = 0.;
  float r2 = 40.;

  Range trajRRange(trajPos.perp() - r1, trajPos.perp() + r2);
  Range trajPhiRange(trajPos.phi() - std::abs(thePhiMin), trajPos.phi() + std::abs(thePhiMax));

  if(rangesIntersect(trajRRange, plane.rSpan()) &&
     rangesIntersect(trajPhiRange, plane.phiSpan(), PhiLess()))
   { return true ; }
  else
   { return false ; }
}

MeasurementEstimator::Local2DVector ForwardMeasurementEstimator::maximalLocalDisplacement	(	const TrajectoryStateOnSurface &	ts,
		const BoundPlane &	plane
	)		const

Definition at line 122 of file ForwardMeasurementEstimator.cc.

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

 {
  float nSigmaCut = 3.;
  if ( ts.hasError())
   {
    LocalError le = ts.localError().positionError();
    return Local2DVector( std::sqrt(le.xx())*nSigmaCut, std::sqrt(le.yy())*nSigmaCut);
   }
  else
    return Local2DVector(999999,999999) ;
 }

void ForwardMeasurementEstimator::setPhiRange ( float  dummyphiMin, float  dummyphiMax  )

inline

Definition at line 37 of file ForwardMeasurementEstimator.h.

References thePhiMax, and thePhiMin.

     { thePhiMin = dummyphiMin ; thePhiMax = dummyphiMax ; }

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

inline

Definition at line 39 of file ForwardMeasurementEstimator.h.

References theRMax, and theRMin.

     { theRMin = rmin ; theRMax = rmax ; }

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

inline

Definition at line 41 of file ForwardMeasurementEstimator.h.

References theRMaxI, and theRMinI.

     { theRMinI = rmin ; theRMaxI = rmax ; }

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().