#include <CSCSegment.h>

Inheritance diagram for CSCSegment:

Public Member Functions
double	chi2 () const override
	Chi2 of the segment fit. More...

CSCSegment *	clone () const override

CSCDetId	cscDetId () const

	CSCSegment ()
	Default constructor. More...

	CSCSegment (const std::vector< const CSCRecHit2D * > &proto_segment, LocalPoint origin, LocalVector direction, const AlgebraicSymMatrix &errors, double chi2)
	Constructor. More...

int	degreesOfFreedom () const override
	Degrees of freedom of the segment fit. More...

int	dimension () const override
	Dimension (in parameter space) More...

const std::vector< CSCSegment > &	duplicateSegments () const

bool	isME11a_duplicate () const

LocalVector	localDirection () const override
	Local direction. More...

LocalError	localDirectionError () const override
	Error on the local direction. More...

LocalPoint	localPosition () const override

LocalError	localPositionError () const override

int	nRecHits () const

AlgebraicVector	parameters () const override
	Parameters of the segment, for the track fit in the order (dx/dz, dy/dz, x, y ) More...

AlgebraicSymMatrix	parametersError () const override
	Covariance matrix of parameters() More...

void	print () const

AlgebraicMatrix	projectionMatrix () const override
	The projection matrix relates the trajectory state parameters to the segment parameters(). More...

std::vector< const TrackingRecHit * >	recHits () const override
	Access to component RecHits (if any) More...

std::vector< TrackingRecHit * >	recHits () override
	Non-const access to component RecHits (if any) More...

void	setDuplicateSegments (std::vector< CSCSegment * > &duplicates)

bool	sharesRecHits (const CSCSegment &anotherSegment, CSCRecHit2D::SharedInputType sharesInput) const

bool	sharesRecHits (const CSCSegment &anotherSegment) const

const std::vector< CSCRecHit2D > &	specificRecHits () const

bool	testSharesAllInSpecificRecHits (const std::vector< CSCRecHit2D > &specificRecHits_1, const std::vector< CSCRecHit2D > &specificRecHits_2, CSCRecHit2D::SharedInputType) const

float	time () const

	~CSCSegment () override
	Destructor. More...

Public Member Functions inherited from RecSegment
	RecSegment (DetId id)

	RecSegment (TrackingRecHit::id_type id=0)

	~RecSegment () override
	Destructor. More...

Public Member Functions inherited from TrackingRecHit
virtual bool	canImproveWithTrack () const

RecHitPointer	cloneForFit (const GeomDet &idet) const

virtual TrackingRecHit *	cloneHit () const

virtual RecHitPointer	cloneSH () const

const GeomDet *	det () const

virtual const GeomDetUnit *	detUnit () const

virtual float	errorGlobalR () const

virtual float	errorGlobalRPhi () const

virtual float	errorGlobalZ () const

DetId	geographicalId () const

virtual void	getKfComponents (KfComponentsHolder &holder) const

unsigned int	getRTTI () const

Type	getType () const

virtual GlobalPoint	globalPosition () const

virtual GlobalError	globalPositionError () const

virtual bool	hasPositionAndError () const
	to be redefined by daughter class More...

virtual TrackingRecHit const *	hit () const

bool	isValid () const

id_type	rawId () const

virtual void	recHitsV (std::vector< const TrackingRecHit * > &) const

virtual void	recHitsV (std::vector< TrackingRecHit * > &)

virtual void	setDet (const GeomDet &idet)

virtual bool	sharesInput (const TrackingRecHit *other, SharedInputType what) const

virtual const Surface *	surface () const

	TrackingRecHit (DetId id, Type type=valid)

	TrackingRecHit (id_type id=0, Type type=valid)

	TrackingRecHit (DetId id, unsigned int rt, Type type=valid)

	TrackingRecHit (const GeomDet &idet, Type type=valid)

	TrackingRecHit (const GeomDet &idet, unsigned int rt, Type type=valid)

	TrackingRecHit (const GeomDet &idet, TrackingRecHit const &rh)

virtual ConstRecHitContainer	transientHits () const

Type	type () const

virtual float	weight () const

virtual	~TrackingRecHit ()

Private Attributes
bool	aME11a_duplicate

double	theChi2

AlgebraicSymMatrix	theCovMatrix

std::vector< CSCRecHit2D >	theCSCRecHits

std::vector< CSCSegment >	theDuplicateSegments

LocalVector	theLocalDirection

LocalPoint	theOrigin

Additional Inherited Members
Public Types inherited from TrackingRecHit
typedef std::vector < ConstRecHitPointer >	ConstRecHitContainer

using	ConstRecHitPointer = std::shared_ptr< TrackingRecHit const >

typedef unsigned int	id_type

typedef std::vector < ConstRecHitPointer >	RecHitContainer

using	RecHitPointer = std::shared_ptr< TrackingRecHit const >

enum	SharedInputType { all, some }
	definition of equality via shared input More...

enum	Type { valid = 0, missing = 1, inactive = 2, bad = 3, missing_inner = 4, missing_outer = 5, inactive_inner = 6, inactive_outer = 7 }

Static Public Attributes inherited from TrackingRecHit
static const int	rttiShift = 24

static const int	typeMask = 0xf

Protected Member Functions inherited from TrackingRecHit
void	setId (id_type iid)

void	setRTTI (unsigned int rt)

void	setType (Type ttype)

Detailed Description

Describes a reconstructed track segment in the 6 layers of a CSC chamber. This is 4-dimensional since it has an origin (x,y) and a direction (x,y) in the local coordinate system of the chamber.

Author: Matteo Sani; Rick Wilkinson; Tim Cox

Definition at line 21 of file CSCSegment.h.

Constructor & Destructor Documentation

CSCSegment::CSCSegment ( )

inline

Default constructor.

Definition at line 24 of file CSCSegment.h.

Referenced by clone().

24 : theChi2(0.), aME11a_duplicate(false) {}

CSCSegment::aME11a_duplicate

bool aME11a_duplicate

Definition: CSCSegment.h:98

CSCSegment::theChi2

double theChi2

Definition: CSCSegment.h:97

CSCSegment::CSCSegment	(	const std::vector< const CSCRecHit2D * > &	proto_segment,
		LocalPoint	origin,
		LocalVector	direction,
		const AlgebraicSymMatrix &	errors,
		double	chi2
	)

Constructor.

Definition at line 15 of file CSCSegment.cc.

References mps_fire::i, and theCSCRecHits.

     : RecSegment(buildDetId(proto_segment.front()->cscDetId())),
       theOrigin(origin),
       theLocalDirection(direction),
       theCovMatrix(errors),
       theChi2(chi2),
       aME11a_duplicate(false) {
   for (unsigned int i = 0; i < proto_segment.size(); ++i)
     theCSCRecHits.push_back(*proto_segment[i]);
 }

CSCSegment::~CSCSegment ( )

override

Destructor.

Definition at line 30 of file CSCSegment.cc.

30 {}

Member Function Documentation

double CSCSegment::chi2 ( ) const

inlineoverridevirtual

Chi2 of the segment fit.

Implements RecSegment.

Definition at line 58 of file CSCSegment.h.

References theChi2.

Referenced by CSCSegAlgoST::buildSegments(), GEMCSCSegAlgoRR::buildSegments(), CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), operator<<(), and MuonTrackCSCChamberResidual::setSegmentResidual().

58 { return theChi2; };

CSCSegment::theChi2

double theChi2

Definition: CSCSegment.h:97

CSCSegment* CSCSegment::clone ( void ) const

inlineoverridevirtual

Implements TrackingRecHit.

Definition at line 37 of file CSCSegment.h.

References CSCSegment().

37 { return new CSCSegment(*this); }

CSCSegment::CSCSegment

CSCSegment()

Default constructor.

Definition: CSCSegment.h:24

CSCDetId CSCSegment::cscDetId ( ) const

inline

Definition at line 70 of file CSCSegment.h.

References TrackingRecHit::geographicalId().

Referenced by FWCSCSegmentProxyBuilder::build(), GEMCSCSegAlgoRR::buildSegments(), DynamicTruncation::chooseLayers(), DynamicTruncation::fillDYTInfos(), GEMCSCSegment::GEMCSCSegment(), CSCRecHitMatcher::globalPoint(), GEMEffByGEMCSCSegmentSource::isME11SegmentMatched(), and DynamicTruncation::preliminaryFit().

70 { return geographicalId(); }

TrackingRecHit::geographicalId

DetId geographicalId() const

Definition: TrackingRecHit.h:120

int CSCSegment::degreesOfFreedom ( ) const

inlineoverridevirtual

Degrees of freedom of the segment fit.

Implements RecSegment.

Definition at line 62 of file CSCSegment.h.

References nRecHits().

Referenced by CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), operator<<(), and MuonTrackCSCChamberResidual::setSegmentResidual().

62 { return 2 * nRecHits() - 4; }

CSCSegment::nRecHits

int nRecHits() const

Definition: CSCSegment.h:68

int CSCSegment::dimension ( ) const

inlineoverridevirtual

Dimension (in parameter space)

Implements RecSegment.

Definition at line 60 of file CSCSegment.h.

60 { return 4; }

const std::vector<CSCSegment>& CSCSegment::duplicateSegments ( ) const

inline

Definition at line 76 of file CSCSegment.h.

References theDuplicateSegments.

Referenced by operator<<().

76 { return theDuplicateSegments; }

CSCSegment::theDuplicateSegments

std::vector< CSCSegment > theDuplicateSegments

Definition: CSCSegment.h:99

bool CSCSegment::isME11a_duplicate ( ) const

inline

Definition at line 74 of file CSCSegment.h.

References theDuplicateSegments.

74 { return (!theDuplicateSegments.empty() ? true : false); }

CSCSegment::theDuplicateSegments

std::vector< CSCSegment > theDuplicateSegments

Definition: CSCSegment.h:99

LocalVector CSCSegment::localDirection ( ) const

inlineoverridevirtual

Local direction.

Implements RecSegment.

Definition at line 42 of file CSCSegment.h.

References theLocalDirection.

Referenced by CSCRecHitMatcher::areCSCSegmentsSame(), FWCSCSegmentProxyBuilder::build(), GEMCSCSegAlgoRR::chainHitsToSegm(), CSCOfflineMonitor::doEfficiencies(), CSCValidation::doEfficiencies(), CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), operator<<(), and CSCSegmentValidation::plotResolution().

42 { return theLocalDirection; }

CSCSegment::theLocalDirection

LocalVector theLocalDirection

Definition: CSCSegment.h:95

LocalError CSCSegment::localDirectionError ( ) const

overridevirtual

Error on the local direction.

Implements RecSegment.

Definition at line 52 of file CSCSegment.cc.

References theCovMatrix.

Referenced by CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), operator<<(), and CSCSegmentValidation::plotResolution().

                                                  {
   return LocalError(theCovMatrix[0][0], theCovMatrix[0][1], theCovMatrix[1][1]);
 }

LocalPoint CSCSegment::localPosition ( ) const

inlineoverridevirtual

Implements TrackingRecHit.

Definition at line 39 of file CSCSegment.h.

References theOrigin.

Referenced by L1MuonRecoTreeProducer::analyze(), CSCRecHitMatcher::areCSCSegmentsSame(), FWCSCSegmentProxyBuilder::build(), GEMCSCSegAlgoRR::chainHitsToSegm(), CSCOfflineMonitor::doEfficiencies(), CSCValidation::doEfficiencies(), CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), CSCRecHitMatcher::globalPoint(), GEMEffByGEMCSCSegmentSource::isME11SegmentMatched(), operator<<(), and CSCSegmentValidation::plotResolution().

39 { return theOrigin; }

CSCSegment::theOrigin

LocalPoint theOrigin

Definition: CSCSegment.h:94

LocalError CSCSegment::localPositionError ( ) const

overridevirtual

Implements TrackingRecHit.

Definition at line 48 of file CSCSegment.cc.

References theCovMatrix.

Referenced by CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), operator<<(), and CSCSegmentValidation::plotResolution().

                                                 {
   return LocalError(theCovMatrix[2][2], theCovMatrix[2][3], theCovMatrix[3][3]);
 }

int CSCSegment::nRecHits ( ) const

inline

Definition at line 68 of file CSCSegment.h.

References theCSCRecHits.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoDF::buildSegments(), degreesOfFreedom(), and MuonTrackCSCChamberResidual::setSegmentResidual().

68 { return theCSCRecHits.size(); }

CSCSegment::theCSCRecHits

std::vector< CSCRecHit2D > theCSCRecHits

Definition: CSCSegment.h:93

AlgebraicVector CSCSegment::parameters ( void ) const

overridevirtual

Parameters of the segment, for the track fit in the order (dx/dz, dy/dz, x, y )

Implements TrackingRecHit.

Definition at line 56 of file CSCSegment.cc.

References mps_fire::result, theLocalDirection, theOrigin, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by StateSegmentMatcher::StateSegmentMatcher().

                                              {
   // For consistency with DT and what we require for the TrackingRecHit interface,
   // the order of the parameters in the returned vector should be (dx/dz, dy/dz, x, z)
 
   AlgebraicVector result(4);
 
   result[0] = theLocalDirection.x() / theLocalDirection.z();
   result[1] = theLocalDirection.y() / theLocalDirection.z();
   result[2] = theOrigin.x();
   result[3] = theOrigin.y();
 
   return result;
 }

AlgebraicSymMatrix CSCSegment::parametersError ( ) const

inlineoverridevirtual

Covariance matrix of parameters()

Implements TrackingRecHit.

Definition at line 49 of file CSCSegment.h.

References theCovMatrix.

Referenced by CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), and StateSegmentMatcher::StateSegmentMatcher().

49 { return theCovMatrix; }

CSCSegment::theCovMatrix

AlgebraicSymMatrix theCovMatrix

Definition: CSCSegment.h:96

void CSCSegment::print ( void ) const

Definition at line 179 of file CSCSegment.cc.

References gather_cfg::cout.

179 { std::cout << *this << std::endl; }

gather_cfg.cout

tuple cout

Definition: gather_cfg.py:144

AlgebraicMatrix CSCSegment::projectionMatrix ( ) const

overridevirtual

The projection matrix relates the trajectory state parameters to the segment parameters().

Implements TrackingRecHit.

Definition at line 81 of file CSCSegment.cc.

References theProjectionMatrix.

81 { return theProjectionMatrix; }

theProjectionMatrix

static const AlgebraicMatrix theProjectionMatrix

Definition: CSCSegment.cc:79

std::vector< const TrackingRecHit * > CSCSegment::recHits ( ) const

overridevirtual

Access to component RecHits (if any)

Implements TrackingRecHit.

Definition at line 32 of file CSCSegment.cc.

References theCSCRecHits.

Referenced by MuonAssociatorByHitsHelper::getMatchedIds(), TrackerMuonHitExtractor::getMuonHits(), and MuonTrackProducer::produce().

                                                            {
   std::vector<const TrackingRecHit*> pointersOfRecHits;
   for (std::vector<CSCRecHit2D>::const_iterator irh = theCSCRecHits.begin(); irh != theCSCRecHits.end(); ++irh) {
     pointersOfRecHits.push_back(&(*irh));
   }
   return pointersOfRecHits;
 }

std::vector< TrackingRecHit * > CSCSegment::recHits ( )

overridevirtual

Non-const access to component RecHits (if any)

Implements TrackingRecHit.

Definition at line 40 of file CSCSegment.cc.

References theCSCRecHits.

                                                {
   std::vector<TrackingRecHit*> pointersOfRecHits;
   for (std::vector<CSCRecHit2D>::iterator irh = theCSCRecHits.begin(); irh != theCSCRecHits.end(); ++irh) {
     pointersOfRecHits.push_back(&(*irh));
   }
   return pointersOfRecHits;
 }

void CSCSegment::setDuplicateSegments ( std::vector< CSCSegment * > & duplicates )

Definition at line 83 of file CSCSegment.cc.

References mps_fire::i, and theDuplicateSegments.

                                                                         {
   theDuplicateSegments.clear();
   for (unsigned int i = 0; i < duplicates.size(); ++i) {
     theDuplicateSegments.push_back(*duplicates[i]);
     //avoid copying duplicates of duplicates of duplicates...
     theDuplicateSegments.back().theDuplicateSegments.resize(0);
   }
 }

bool CSCSegment::sharesRecHits	(	const CSCSegment &	anotherSegment,
		CSCRecHit2D::SharedInputType	sharesInput
	)		const

Definition at line 129 of file CSCSegment.cc.

References TrackingRecHit::sharesInput(), specificRecHits(), testSharesAllInSpecificRecHits(), and theCSCRecHits.

                                                                                                              {
   return testSharesAllInSpecificRecHits(theCSCRecHits, anotherSegment.specificRecHits(), sharesInput);
 }

bool CSCSegment::sharesRecHits ( const CSCSegment & anotherSegment ) const

Definition at line 134 of file CSCSegment.cc.

References CSCRecHit2D::someStrips, CSCRecHit2D::someWires, specificRecHits(), testSharesAllInSpecificRecHits(), and theCSCRecHits.

                                                                      {
   if (testSharesAllInSpecificRecHits(theCSCRecHits, anotherSegment.specificRecHits(), CSCRecHit2D::someWires) &&
       testSharesAllInSpecificRecHits(theCSCRecHits, anotherSegment.specificRecHits(), CSCRecHit2D::someStrips)) {
     return true;
   } else {
     return false;
   }
 }

const std::vector<CSCRecHit2D>& CSCSegment::specificRecHits ( ) const

inline

Definition at line 66 of file CSCSegment.h.

References theCSCRecHits.

Referenced by GEMCSCSegAlgoRR::buildSegments(), GEMCSCSegAlgoRR::chainHitsToSegm(), GEMCSCSegment::cscRecHits(), CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), GEMCSCSegment::GEMCSCSegment(), GEMCSCSegment::nRecHits(), operator<<(), GEMCSCSegment::recHits(), sharesRecHits(), and DynamicTruncation::updateWithCSChits().

66 { return theCSCRecHits; }

CSCSegment::theCSCRecHits

std::vector< CSCRecHit2D > theCSCRecHits

Definition: CSCSegment.h:93

bool CSCSegment::testSharesAllInSpecificRecHits	(	const std::vector< CSCRecHit2D > &	specificRecHits_1,
		const std::vector< CSCRecHit2D > &	specificRecHits_2,
		CSCRecHit2D::SharedInputType	sharesInput
	)		const

Definition at line 92 of file CSCSegment.cc.

Referenced by sharesRecHits().

                                                                                               {
   const std::vector<CSCRecHit2D>* rhContainer_1 = &specificRecHits_1;
   const std::vector<CSCRecHit2D>* rhContainer_2 = &specificRecHits_2;
   if (specificRecHits_1.size() > specificRecHits_2.size()) {
     rhContainer_2 = &specificRecHits_1;
     rhContainer_1 = &specificRecHits_2;
   }
   //
   bool shareConditionPassed = true;
   for (std::vector<CSCRecHit2D>::const_iterator itRH = rhContainer_1->begin(); itRH != rhContainer_1->end(); ++itRH) {
     const CSCRecHit2D* firstRecHit = &(*itRH);
     bool sharedHit = false;
     for (std::vector<CSCRecHit2D>::const_iterator itRH2 = rhContainer_2->begin(); itRH2 != rhContainer_2->end();
          ++itRH2) {
       if (itRH2->sharesInput(firstRecHit, sharesInput)) {
         sharedHit = true;
         break;
       }
     }
     if (!sharedHit) {
       shareConditionPassed = false;
       break;
     }
   }
   return shareConditionPassed;
 }

float CSCSegment::time ( ) const

Definition at line 144 of file CSCSegment.cc.

References change_name::diff, reco::castor::maxDiff(), N, theCSCRecHits, CSCRecHit2D::tpeak(), and CSCRecHit2D::wireTime().

Referenced by CSCSegAlgoST::dumpSegment(), CSCSegAlgoDF::dumpSegment(), CSCSegAlgoSK::dumpSegment(), CSCSegAlgoTC::dumpSegment(), and GEMEffByGEMCSCSegmentSource::isME11SegmentMatched().

                              {
   float averageTime = 0;
   std::vector<float> wireTimes;
   for (std::vector<CSCRecHit2D>::const_iterator itRH = theCSCRecHits.begin(); itRH != theCSCRecHits.end(); ++itRH) {
     const CSCRecHit2D* recHit = &(*itRH);
     averageTime += recHit->tpeak();
     averageTime += recHit->wireTime();
     wireTimes.push_back(recHit->wireTime());
   }
   averageTime = averageTime / (2 * theCSCRecHits.size());
 
   //The wire times have a long tail that has to be pruned.  The strip times (tpeak) are fine
   bool modified = true;
   while (modified) {
     modified = false;
     double maxDiff = -1;
     std::vector<float>::iterator maxHit;
     for (std::vector<float>::iterator itWT = wireTimes.begin(); itWT != wireTimes.end(); ++itWT) {
       float diff = fabs(*itWT - averageTime);
       if (diff > maxDiff) {
         maxDiff = diff;
         maxHit = itWT;
       }
     }
     if (maxDiff > 26) {
       int N = theCSCRecHits.size() + wireTimes.size();
       averageTime = (averageTime * N - (*maxHit)) / (N - 1);
       wireTimes.erase(maxHit);
       modified = true;
     }
   }
   return averageTime;
 }