BarrelMeasurementEstimator Class Reference

#include <BarrelMeasurementEstimator.h>

Inheritance diagram for BarrelMeasurementEstimator:

Public Member Functions
	BarrelMeasurementEstimator ()
	BarrelMeasurementEstimator (float phiMin, float phiMax, float zMin, float zMax)
BarrelMeasurementEstimator *	clone () 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
MeasurementEstimator::Local2DVector	maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const BoundPlane &plane) const override
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

BarrelMeasurementEstimator* BarrelMeasurementEstimator::clone ( void  ) const

inlineoverridevirtual

Implements MeasurementEstimator.

Definition at line 52 of file BarrelMeasurementEstimator.h.

References BarrelMeasurementEstimator(), and maximalLocalDisplacement().

     { return new BarrelMeasurementEstimator(*this) ; }

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

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

Referenced by estimate(), and setZRange().

                                                                  {
  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 36 of file BarrelMeasurementEstimator.cc.

References funct::abs(), estimate(), TrajectoryStateOnSurface::globalParameters(), normalizedPhi(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), 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 = normalizedPhi(tsPhi - rhPhi);

  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 70 of file BarrelMeasurementEstimator.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::barePhi(), estimate(), f, TrajectoryStateOnSurface::globalParameters(), runTauDisplay::gp, normalizedPhi(), PV3DBase< T, PVType, FrameType >::perp(), GlobalTrajectoryParameters::position(), thePhiMax, thePhiMin, theZMax, theZMin, 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 = normalizedPhi(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

override

Definition at line 102 of file BarrelMeasurementEstimator.cc.

References funct::abs(), funct::false, TrajectoryStateOnSurface::globalParameters(), maximalLocalDisplacement(), Geom::phiLess(), GlobalTrajectoryParameters::position(), rangesIntersect(), thePhiMax, thePhiMin, theZMax, theZMin, x, and y.

 {
  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(), [](auto x,auto y){ return Geom::phiLess(x, y);}))
   {
    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

override

Definition at line 135 of file BarrelMeasurementEstimator.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().

 {
  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 estimate(), runTauDisplay::gp, theZMax, and theZMin.

     { theZMin=zmin ; theZMax=zmax ; }

Member Data Documentation

float BarrelMeasurementEstimator::thePhiMax

private

Definition at line 61 of file BarrelMeasurementEstimator.h.

Referenced by estimate(), and 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().