#include <Chi2ChargeMeasurementEstimator.h>

Inheritance diagram for Chi2ChargeMeasurementEstimator:

Public Member Functions
	Chi2ChargeMeasurementEstimator (double maxChi2, double nSigma, bool cutOnPixelCharge, bool cutOnStripCharge, double minGoodPixelCharge, double minGoodStripCharge, float pTChargeCutThreshold=100000.)

virtual Chi2ChargeMeasurementEstimator *	clone () const

virtual std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &, const TrackingRecHit &) const

Public Member Functions inherited from Chi2MeasurementEstimator
	Chi2MeasurementEstimator (double maxChi2, double nSigma=3.)

Public Member Functions inherited from Chi2MeasurementEstimatorBase
	Chi2MeasurementEstimatorBase (double maxChi2, double nSigma=3.)

double	chiSquaredCut () const

virtual bool	estimate (const TrajectoryStateOnSurface &ts, const Plane &plane) const

virtual Local2DVector	maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const Plane &plane) const

double	nSigmaCut () const

Public Member Functions inherited from MeasurementEstimator
virtual	~MeasurementEstimator ()

Private Member Functions
bool	checkCharge (const TrackingRecHit &aRecHit, int subdet, float chargeCut) const

bool	checkClusterCharge (const OmniClusterRef::ClusterStripRef cluster, float chargeCut) const

double	minGoodCharge (int subdet) const

float	sensorThickness (const DetId &detid) const

bool	thickSensors (const SiStripDetId &detid) const

Private Attributes
bool	cutOnPixelCharge_

bool	cutOnStripCharge_

double	minGoodPixelCharge_

double	minGoodStripCharge_

float	pTChargeCutThreshold_

Additional Inherited Members
Public Types inherited from MeasurementEstimator
typedef std::pair< bool, double >	HitReturnType

typedef Vector2DBase< float, LocalTag >	Local2DVector

typedef bool	SurfaceReturnType

Protected Member Functions inherited from Chi2MeasurementEstimatorBase
std::pair< bool, double >	returnIt (double est) const

Detailed Description

A Chi2 Measurement Estimator, checking also the charge of the cluster. Computhes the Chi^2 of a TrajectoryState with a RecHit or a Plane. The TrajectoryState must have errors. If the cluster passes the chi2 cut, the charge is checked. This cut can be bypassed for high-pt cut. Works for any RecHit dimension.

Author: todorov, cerati, speer

Definition at line 20 of file Chi2ChargeMeasurementEstimator.h.

Constructor & Destructor Documentation

Chi2ChargeMeasurementEstimator::Chi2ChargeMeasurementEstimator	(	double	maxChi2,
		double	nSigma,
		bool	cutOnPixelCharge,
		bool	cutOnStripCharge,
		double	minGoodPixelCharge,
		double	minGoodStripCharge,
		float	pTChargeCutThreshold = `100000.`
	)

inlineexplicit

Construct with cuts on chi2 and nSigma. The cut on Chi2 is used to define the acceptance of RecHits. The errors of the trajectory state are multiplied by nSigma to define acceptance of Plane and maximalLocalDisplacement.

Definition at line 28 of file Chi2ChargeMeasurementEstimator.h.

References pTChargeCutThreshold_.

Referenced by clone().

                                           : 
     Chi2MeasurementEstimator( maxChi2, nSigma), cutOnPixelCharge_(cutOnPixelCharge),
     cutOnStripCharge_(cutOnStripCharge), minGoodPixelCharge_(minGoodPixelCharge),
     minGoodStripCharge_(minGoodStripCharge) {
       if (pTChargeCutThreshold>=0.) pTChargeCutThreshold_=pTChargeCutThreshold;
       else pTChargeCutThreshold_=100000;
     }

Member Function Documentation

bool Chi2ChargeMeasurementEstimator::checkCharge	(	const TrackingRecHit &	aRecHit,
		int	subdet,
		float	chargeCut
	)		const

private

Definition at line 18 of file Chi2ChargeMeasurementEstimator.cc.

References checkClusterCharge(), OmniClusterRef::cluster_strip(), edm::hlt::Exception, BaseTrackerRecHit::firstClusterRef(), TrackingRecHit::getRTTI(), TrackingRecHit::hit(), likely, SiStripMatchedRecHit2D::monoClusterRef(), and SiStripMatchedRecHit2D::stereoClusterRef().

Referenced by estimate().

 {
   auto const & hit = aRecHit.hit();   // this works for both TRH and TTRH! 
 
   if (aRecHit.getRTTI() == 4) {
     const SiStripMatchedRecHit2D *matchHit = static_cast<const SiStripMatchedRecHit2D *>(hit);
     if likely(matchHit!=0) return (checkClusterCharge(matchHit->monoClusterRef().cluster_strip(), chargeCut)
       &&checkClusterCharge(matchHit->stereoClusterRef().cluster_strip(), chargeCut));
     else{
       const std::type_info &type = typeid(aRecHit);
       throw cms::Exception("Unknown RecHit Type") << "checkCharge/SiStripMatchedRecHit2D: Wrong recHit type: " << type.name();
     }
   } else {
     auto const & thit = static_cast<const BaseTrackerRecHit *>(hit);
     auto const & clus = thit->firstClusterRef();
     if (!clus.isPixel()) return checkClusterCharge(clus.cluster_strip(), chargeCut);
     else return (clus.cluster_pixel()->charge()>chargeCut);
   }
 }

bool Chi2ChargeMeasurementEstimator::checkClusterCharge	(	const OmniClusterRef::ClusterStripRef	cluster,
		float	chargeCut
	)		const

private

Definition at line 12 of file Chi2ChargeMeasurementEstimator.cc.

Referenced by checkCharge().

 {
   int clusCharge=accumulate( cluster->amplitudes().begin(), cluster->amplitudes().end(), uint16_t(0));
   return (clusCharge>chargeCut);
 }

virtual Chi2ChargeMeasurementEstimator* Chi2ChargeMeasurementEstimator::clone ( void ) const

inlinevirtual

Reimplemented from Chi2MeasurementEstimator.

Definition at line 43 of file Chi2ChargeMeasurementEstimator.h.

References Chi2ChargeMeasurementEstimator().

                                                         {
     return new Chi2ChargeMeasurementEstimator(*this);
   }

std::pair< bool, double > Chi2ChargeMeasurementEstimator::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.

Reimplemented from Chi2MeasurementEstimator.

Definition at line 56 of file Chi2ChargeMeasurementEstimator.cc.

References checkCharge(), cutOnPixelCharge_, cutOnStripCharge_, DetId::det(), cond::rpcobgas::detid, LocalTrajectoryParameters::dxdz(), LocalTrajectoryParameters::dydz(), Chi2MeasurementEstimator::estimate(), TrackingRecHit::geographicalId(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::localParameters(), minGoodCharge(), PV3DBase< T, PVType, FrameType >::perp(), pTChargeCutThreshold_, sensorThickness(), mathSSE::sqrt(), and DetId::subdetId().

                                                         {
 
 
 
   std::pair<bool,double> estimateResult = Chi2MeasurementEstimator::estimate(tsos, aRecHit);
   if ( !estimateResult.first || (!(cutOnStripCharge_||cutOnPixelCharge_)) || 
     (tsos.globalMomentum ().perp()>pTChargeCutThreshold_)) return estimateResult;
 
   SiStripDetId detid = aRecHit.geographicalId(); 
   uint32_t subdet = detid.subdetId();
   if((detid.det()==1) && (((subdet>2)&&cutOnStripCharge_) || ((subdet<3)&&cutOnPixelCharge_)))   {
 
     float theDxdz = tsos.localParameters().dxdz();
     float theDydz = tsos.localParameters().dydz();
     float chargeCut = (float) minGoodCharge(subdet) * sensorThickness(detid)*(sqrt(1. + theDxdz*theDxdz + theDydz*theDydz));
 
     if (!checkCharge(aRecHit, subdet, chargeCut)) return HitReturnType(false,0.);
     else return estimateResult;
   }
 
   return estimateResult;
 }

double Chi2ChargeMeasurementEstimator::minGoodCharge ( int subdet ) const

inlineprivate

Definition at line 52 of file Chi2ChargeMeasurementEstimator.h.

References minGoodPixelCharge_, and minGoodStripCharge_.

Referenced by estimate().

52 {return (subdet>2?minGoodStripCharge_:minGoodPixelCharge_);}

Chi2ChargeMeasurementEstimator::minGoodPixelCharge_

double minGoodPixelCharge_

Definition: Chi2ChargeMeasurementEstimator.h:49

Chi2ChargeMeasurementEstimator::minGoodStripCharge_

double minGoodStripCharge_

Definition: Chi2ChargeMeasurementEstimator.h:50

float Chi2ChargeMeasurementEstimator::sensorThickness ( const DetId & detid ) const

private

Definition at line 41 of file Chi2ChargeMeasurementEstimator.cc.

References cond::rpcobgas::detid, SiStripDetId::moduleGeometry(), DetId::subdetId(), SiStripDetId::TIB, SiStripDetId::TOB, SiStripDetId::W5, SiStripDetId::W6, and SiStripDetId::W7.

Referenced by estimate().

 {
   if (detid.subdetId()>=SiStripDetId::TIB) {
     SiStripDetId siStripDetId = detid(); 
     if (siStripDetId.subdetId()==SiStripDetId::TOB) return 0.047;
     if (siStripDetId.moduleGeometry()==SiStripDetId::W5 || siStripDetId.moduleGeometry()==SiStripDetId::W6 ||
         siStripDetId.moduleGeometry()==SiStripDetId::W7)
     return 0.047;
     return 0.029; // so it is TEC ring 1-4 or TIB or TOB;
   } else if (detid.subdetId()==1) return 0.0285;
   else return 0.027;
 }