BarrelMeasurementEstimator Class Reference

#include <BarrelMeasurementEstimator.h>

Inheritance diagram for BarrelMeasurementEstimator:

Public Member Functions
	BarrelMeasurementEstimator ()
	BarrelMeasurementEstimator (float phiMin, float phiMax, float zMin, float zMax)
virtual BarrelMeasurementEstimator *	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
MeasurementEstimator::Local2DVector	maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const
void	setPhiRange (float dummyphiMin, float dummyphiMax)
void	setZRange (float zmin, float zmax)
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	theZMax
float	theZMin

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

Detailed Description

Definition at line 31 of file BarrelMeasurementEstimator.h.

Constructor & Destructor Documentation

BarrelMeasurementEstimator::BarrelMeasurementEstimator ( )

inline

Definition at line 35 of file BarrelMeasurementEstimator.h.

Referenced by clone().

     {}

BarrelMeasurementEstimator::BarrelMeasurementEstimator ( float  phiMin, float  phiMax, float  zMin, float  zMax  )

inline

Definition at line 37 of file BarrelMeasurementEstimator.h.

     : thePhiMin(phiMin), thePhiMax(phiMax), theZMin(zMin), theZMax(zMax)
     {}

Member Function Documentation

virtual BarrelMeasurementEstimator* BarrelMeasurementEstimator::clone ( void  ) const

inlinevirtual

Implements MeasurementEstimator.

Definition at line 52 of file BarrelMeasurementEstimator.h.

References BarrelMeasurementEstimator().

     { return new BarrelMeasurementEstimator(*this) ; }

std::pair< bool, double > BarrelMeasurementEstimator::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 27 of file BarrelMeasurementEstimator.cc.

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

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

std::pair< bool, double > BarrelMeasurementEstimator::estimate ( const TrajectoryStateOnSurface &  ts, const GlobalPoint &  gp  ) const

virtual

Definition at line 35 of file BarrelMeasurementEstimator.cc.

References funct::abs(), TrajectoryStateOnSurface::globalParameters(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), pi, GlobalTrajectoryParameters::position(), thePhiMin, theZMax, theZMin, and PV3DBase< T, PVType, FrameType >::z().

                                                    {
  
  float myZ = gp.z();
  float myR = gp.perp();
  float zDiff = ts.globalParameters().position().z() - myZ;

  float myZmax =  theZMax;
  float myZmin =  theZMin;

  if(std::abs(myZ)<30. && myR>8.) {
    myZmax = 0.09;
    myZmin = -0.09;
  }

  if( zDiff >= myZmax || zDiff <= myZmin ) return std::pair<bool,double>(false,0.);

  float tsPhi = ts.globalParameters().position().phi();  

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

  if ( phiDiff < thePhiMax && phiDiff > thePhiMin ) {

    return std::pair<bool,double>(true,1.);
     } else {

    return std::pair<bool,double>(false,0.);
    }
}

std::pair< bool, double > BarrelMeasurementEstimator::estimate ( const GlobalPoint &  vprim, const TrajectoryStateOnSurface &  ts, const GlobalPoint &  gp  ) const

virtual

Definition at line 71 of file BarrelMeasurementEstimator.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::barePhi(), f, TrajectoryStateOnSurface::globalParameters(), 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 myR = gp.perp();
  float myZ = gp.z();
  float zDiff = myZ -ts.z() ;  
  float myZmax =  theZMax;
  float myZmin =  theZMin;
  if( (std::abs(myZ)<30.f) & (myR>8.f) )
    {
      myZmax = 0.09f;
      myZmin = -0.09f;
    }
  

  if( zDiff >= myZmax || zDiff <= myZmin ) return std::pair<bool,double>(false,0.);

  float rhPhi = gp.barePhi() ;
  float tsPhi = ts.barePhi();  
  float phiDiff = normalized_phi(rhPhi-tsPhi) ;

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

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

virtual

Definition at line 103 of file BarrelMeasurementEstimator.cc.

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

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


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

  Range trajZRange(trajPos.z() - std::abs(theZMin), trajPos.z() + std::abs(theZMax));
  Range trajPhiRange(trajPos.phi() - std::abs(thePhiMin), trajPos.phi() + std::abs(thePhiMax));

  if(rangesIntersect(trajZRange, plane.zSpan()) &&
     rangesIntersect(trajPhiRange, plane.phiSpan(), PhiLess()))
   {
    return true;
   }
  else
   {
    //     cout <<cout<<" barrel boundpl est returns false!!"<<endl;
    //     cout<<"BarrelMeasurementEstimator(estimate) :thePhiMin,thePhiMax, theZMin,theZMax "<<thePhiMin<<" "<<thePhiMax<<" "<< theZMin<<" "<<theZMax<<endl;
    //     cout<<" trajZRange "<<trajZRange.first<<" "<<trajZRange.second<<endl;
    //     cout<<" trajPhiRange "<<trajPhiRange.first<<" "<<trajPhiRange.second<<endl;
    //     cout<<" detZRange "<<detRange.zRange().first<<" "<<detRange.zRange().second<<endl;
    //     cout<<" detPhiRange "<<detRange.phiRange().first<<" "<<detRange.phiRange().second<<endl;
    //     cout<<" intersect z: "<<rangesIntersect(trajZRange, detRange.zRange())<<endl;
    //     cout<<" intersect phi: "<<rangesIntersect(trajPhiRange, detRange.phiRange(), PhiLess())<<endl;
    return false ;
   }

 }

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

Definition at line 136 of file BarrelMeasurementEstimator.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(99999,99999) ;
 }

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

inline

Definition at line 41 of file BarrelMeasurementEstimator.h.

References thePhiMax, and thePhiMin.

     { thePhiMin = dummyphiMin ; thePhiMax = dummyphiMax ; }

void BarrelMeasurementEstimator::setZRange ( float  zmin, float  zmax  )

inline

Definition at line 43 of file BarrelMeasurementEstimator.h.

References theZMax, theZMin, SiStripMonitorClusterAlca_cfi::zmax, and SiStripMonitorClusterAlca_cfi::zmin.

     { theZMin=zmin ; theZMax=zmax ; }

Member Data Documentation

float BarrelMeasurementEstimator::thePhiMax

private

Definition at line 61 of file BarrelMeasurementEstimator.h.

Referenced by setPhiRange().

float BarrelMeasurementEstimator::thePhiMin

private

Definition at line 60 of file BarrelMeasurementEstimator.h.

Referenced by estimate(), and setPhiRange().

float BarrelMeasurementEstimator::theZMax

private

Definition at line 63 of file BarrelMeasurementEstimator.h.

Referenced by estimate(), and setZRange().

float BarrelMeasurementEstimator::theZMin

private

Definition at line 62 of file BarrelMeasurementEstimator.h.

Referenced by estimate(), and setZRange().