ConversionBarrelEstimator Class Reference

#include <ConversionBarrelEstimator.h>

Inheritance diagram for ConversionBarrelEstimator:

Public Member Functions
ConversionBarrelEstimator *	clone () const override
	ConversionBarrelEstimator ()
	ConversionBarrelEstimator (float phiRangeMin, float phiRangeMax, float zRangeMin, float zRangeMax, double nSigma=3.)
bool	estimate (const TrajectoryStateOnSurface &ts, const Plane &plane) const override
std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &ts, const TrackingRecHit &hit) const override
Local2DVector	maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const Plane &plane) const override
double	nSigmaCut () const
Public Member Functions inherited from MeasurementEstimator
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
double	theNSigma
float	thePhiRangeMax
float	thePhiRangeMin
float	theZRangeMax
float	theZRangeMin

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

Detailed Description

Defines the search area in the barrel

Author

Nancy Marinelli, U. of Notre Dame, US

Definition at line 18 of file ConversionBarrelEstimator.h.

Constructor & Destructor Documentation

ConversionBarrelEstimator() [1/2]

ConversionBarrelEstimator::ConversionBarrelEstimator ( )

inline

Definition at line 20 of file ConversionBarrelEstimator.h.

{};

Referenced by clone().

ConversionBarrelEstimator() [2/2]

ConversionBarrelEstimator::ConversionBarrelEstimator ( float  phiRangeMin, float  phiRangeMax, float  zRangeMin, float  zRangeMax, double  nSigma = 3.  )

inline

Definition at line 21 of file ConversionBarrelEstimator.h.

      : thePhiRangeMin(phiRangeMin),
        thePhiRangeMax(phiRangeMax),
        theZRangeMin(zRangeMin),
        theZRangeMax(zRangeMax),
        theNSigma(nSigma) {
    //    std::cout << " ConversionBarrelEstimator CTOR " << std::endl;
  }

Member Function Documentation

clone()

ConversionBarrelEstimator* ConversionBarrelEstimator::clone ( void  ) const

inlineoverridevirtual

Implements MeasurementEstimator.

Definition at line 33 of file ConversionBarrelEstimator.h.

{ return new ConversionBarrelEstimator(*this); }

References ConversionBarrelEstimator().

estimate() [1/2]

bool ConversionBarrelEstimator::estimate ( const TrajectoryStateOnSurface &  ts, const Plane &  plane  ) const

overridevirtual

Returns true if the TrajectoryStateOnSurface is compatible with the Plane, false otherwise. The TrajectoryStateOnSurface must be on the plane.

Implements MeasurementEstimator.

estimate() [2/2]

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

                                                                                             {
  std::pair<bool, double> result;
 
  //std::cout << "  ConversionBarrelEstimator::estimate( const TrajectoryStateOnSurface& ts, const TransientTrackingRecHit& hit) " << std::endl;
 
  float tsPhi = ts.globalParameters().position().phi();
  GlobalPoint gp = hit.globalPosition();
  float rhPhi = gp.phi();
 
  // allow a z fudge of 2 sigma
  float dz = 2. * sqrt(hit.localPositionError().yy());
  float zDiff = ts.globalParameters().position().z() - gp.z();
  float phiDiff = tsPhi - rhPhi;
  if (phiDiff > pi)
    phiDiff -= twopi;
  if (phiDiff < -pi)
    phiDiff += twopi;
 
  // add the errors on the window and the point in quadrature
  float zrange = sqrt(theZRangeMax * theZRangeMax + dz * dz);
 
  /*
  std::cout << "  BarrelEstimator ts local error " <<ts.localError().positionError()  << " hit local error " << hit.localPositionError() << std::endl; 
  std::cout << "  BarrelEstimator:  RecHit at " << gp << " phi " << rhPhi << " eta " << gp.eta() <<  std::endl;
  std::cout << "  BarrelEstimator:  ts at " << ts.globalParameters().position() << " phi " <<ts.globalParameters().position().phi() << " eta " << ts.globalParameters().position().eta()<<  std::endl;
  std::cout << "                    zrange = +/-" << zrange << ", zDiff = " << zDiff << std::endl;
  std::cout << "                    thePhiRangeMin = " << thePhiRangeMin << ", thePhiRangeMax = " << thePhiRangeMax << ", phiDiff = " << phiDiff << std::endl;
  */
 
  if (phiDiff < thePhiRangeMax && phiDiff > thePhiRangeMin && zDiff < zrange && zDiff > -zrange) {
    //    std::cout << "      estimator returns 1 with phiDiff " << thePhiRangeMin << " < " << phiDiff << " < "
    //    << thePhiRangeMax << " and zDiff " << zDiff << " < " << zrange << std::endl;
    // std::cout << " YES " << phiDiff << " " << zDiff << std::endl;
    // std::cout << "                  => RECHIT ACCEPTED " << std::endl;
 
    result.first = true;
    result.second = phiDiff;
  } else {
    //     std::cout << "      estimator returns NOT ACCEPTED  with phiDiff " << thePhiRangeMin << " < " << phiDiff << " < "
    //<< thePhiRangeMax << " and zDiff " << zDiff << " < " << theZRangeMax+dz << std::endl;
 
    result.first = false;
    result.second = 0;
  }
 
  return result;
}

References PVValHelper::dz, TrajectoryStateOnSurface::globalParameters(), runTauDisplay::gp, PV3DBase< T, PVType, FrameType >::phi(), pi, GlobalTrajectoryParameters::position(), mps_fire::result, mathSSE::sqrt(), thePhiRangeMin, theZRangeMax, PV3DBase< T, PVType, FrameType >::z(), and TkClusParameters_cff::zrange.

maximalLocalDisplacement()

MeasurementEstimator::Local2DVector ConversionBarrelEstimator::maximalLocalDisplacement ( const TrajectoryStateOnSurface &  ts, const Plane &  plane  ) const

overridevirtual

Returns the size of the compatibility region around the local position of the TrajectoryStateOnSurface along the directions of local x and y axis. The TrajectoryStateOnSurface must be on the plane. This method allows to limit the search for compatible detectors or RecHits. The MeasurementEstimator should not return "true" for any RecHit or Plane which is entirely outside of the compatibility region defined by maximalLocalDisplacement().

Implements MeasurementEstimator.

Definition at line 78 of file ConversionBarrelEstimator.cc.

                                                                       {
  /* 
  if ( ts.hasError() ) {
    LocalError le = ts.localError().positionError();
    std::cout << "  ConversionBarrelEstimator::maximalLocalDisplacent local error " << sqrt(le.xx()) << " " << sqrt(le.yy()) << " nSigma " << nSigmaCut() << " sqrt(le.xx())*nSigmaCut() " << sqrt(le.xx())*nSigmaCut()  << "  sqrt(le.yy())*nSigmaCut() " <<  sqrt(le.yy())*nSigmaCut() << std::endl;
    return Local2DVector( sqrt(le.xx())*nSigmaCut(), sqrt(le.yy())*nSigmaCut());
  }
 
  else return Local2DVector(9999,9999);
  */
  return Local2DVector(9999, 9999);
}

nSigmaCut()

double ConversionBarrelEstimator::nSigmaCut ( ) const

inline

Definition at line 37 of file ConversionBarrelEstimator.h.

{ return theNSigma; }

References theNSigma.

Member Data Documentation

theNSigma

double ConversionBarrelEstimator::theNSigma

private

Definition at line 44 of file ConversionBarrelEstimator.h.

Referenced by nSigmaCut().

thePhiRangeMax

float ConversionBarrelEstimator::thePhiRangeMax

private

Definition at line 41 of file ConversionBarrelEstimator.h.

thePhiRangeMin

float ConversionBarrelEstimator::thePhiRangeMin

private

Definition at line 40 of file ConversionBarrelEstimator.h.

Referenced by estimate().

theZRangeMax

float ConversionBarrelEstimator::theZRangeMax

private

Definition at line 43 of file ConversionBarrelEstimator.h.

Referenced by estimate().

theZRangeMin

float ConversionBarrelEstimator::theZRangeMin

private

Definition at line 42 of file ConversionBarrelEstimator.h.