#include <CSCSegment.h>

Inheritance diagram for CSCSegment:

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

CSCSegment *	clone () const

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...

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

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

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

bool	isME11a_duplicate () const

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

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

LocalPoint	localPosition () const

LocalError	localPositionError () const

int	nRecHits () const

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

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

void	print () const

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

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

virtual std::vector < TrackingRecHit * >	recHits ()
	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

virtual	~CSCSegment ()
	Destructor. More...

Public Member Functions inherited from RecSegment
	RecSegment (DetId id)

	RecSegment (TrackingRecHit::id_type id=0)

virtual	~RecSegment ()
	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 }

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 26 of file CSCSegment.h.

Referenced by clone().

26 : 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 18 of file CSCSegment.cc.

References 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 ( )

virtual

Destructor.

Definition at line 28 of file CSCSegment.cc.

28 {}

Member Function Documentation

double CSCSegment::chi2 ( ) const

inlinevirtual

Chi2 of the segment fit.

Implements RecSegment.

Definition at line 57 of file CSCSegment.h.

References theChi2.

Referenced by CSCSegAlgoST::buildSegments(), CSCSegAlgoST::dumpSegment(), operator<<(), CSCSegAlgoST::prune_bad_hits(), CSCSegAlgoHitPruning::pruneBadHits(), and MuonTrackCSCChamberResidual::setSegmentResidual().

57 { return theChi2; };

CSCSegment::theChi2

double theChi2

Definition: CSCSegment.h:97

CSCSegment* CSCSegment::clone ( void ) const

inlinevirtual

Implements TrackingRecHit.

Definition at line 36 of file CSCSegment.h.

References CSCSegment().

36 { return new CSCSegment(*this); }

CSCSegment::CSCSegment

CSCSegment()

Default constructor.

Definition: CSCSegment.h:26

CSCDetId CSCSegment::cscDetId ( ) const

inline

Definition at line 69 of file CSCSegment.h.

References TrackingRecHit::geographicalId().

Referenced by FWCSCSegmentProxyBuilder::build(), DynamicTruncation::chooseLayers(), DynamicTruncation::fillDYTInfos(), and DynamicTruncation::preliminaryFit().

69 { return geographicalId(); }

TrackingRecHit::geographicalId

DetId geographicalId() const

Definition: TrackingRecHit.h:119

virtual int CSCSegment::degreesOfFreedom ( ) const

inlinevirtual

Degrees of freedom of the segment fit.

Implements RecSegment.

Definition at line 61 of file CSCSegment.h.

References nRecHits().

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

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

CSCSegment::nRecHits

int nRecHits() const

Definition: CSCSegment.h:67

virtual int CSCSegment::dimension ( ) const

inlinevirtual

Dimension (in parameter space)

Implements RecSegment.

Definition at line 59 of file CSCSegment.h.

59 { return 4; }

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

inline

Definition at line 75 of file CSCSegment.h.

References theDuplicateSegments.

Referenced by operator<<().

75 { return theDuplicateSegments; }

CSCSegment::theDuplicateSegments

std::vector< CSCSegment > theDuplicateSegments

Definition: CSCSegment.h:99

bool CSCSegment::isME11a_duplicate ( ) const

inline

Definition at line 73 of file CSCSegment.h.

References theDuplicateSegments.

73 { return (theDuplicateSegments.size() > 0 ? true : false); }

CSCSegment::theDuplicateSegments

std::vector< CSCSegment > theDuplicateSegments

Definition: CSCSegment.h:99

LocalVector CSCSegment::localDirection ( ) const

inlinevirtual

Local direction.

Implements RecSegment.

Definition at line 41 of file CSCSegment.h.

References theLocalDirection.

Referenced by FWCSCSegmentProxyBuilder::build(), CSCOfflineMonitor::doEfficiencies(), CSCValidation::doEfficiencies(), CSCSegAlgoST::dumpSegment(), operator<<(), and CSCSegmentValidation::plotResolution().

41 { return theLocalDirection; }

CSCSegment::theLocalDirection

LocalVector theLocalDirection

Definition: CSCSegment.h:95

LocalError CSCSegment::localDirectionError ( ) const

virtual

Error on the local direction.

Implements RecSegment.

Definition at line 51 of file CSCSegment.cc.

References theCovMatrix.

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

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

LocalPoint CSCSegment::localPosition ( ) const

inlinevirtual

Implements TrackingRecHit.

Definition at line 38 of file CSCSegment.h.

References theOrigin.

Referenced by FWCSCSegmentProxyBuilder::build(), CSCOfflineMonitor::doEfficiencies(), CSCValidation::doEfficiencies(), CSCSegAlgoST::dumpSegment(), operator<<(), and CSCSegmentValidation::plotResolution().

38 { return theOrigin; }

CSCSegment::theOrigin

LocalPoint theOrigin

Definition: CSCSegment.h:94

LocalError CSCSegment::localPositionError ( ) const

virtual

Implements TrackingRecHit.

Definition at line 47 of file CSCSegment.cc.

References theCovMatrix.

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

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

int CSCSegment::nRecHits ( ) const

inline

Definition at line 67 of file CSCSegment.h.

References theCSCRecHits.

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

67 { return theCSCRecHits.size(); }

CSCSegment::theCSCRecHits

std::vector< CSCRecHit2D > theCSCRecHits

Definition: CSCSegment.h:93

AlgebraicVector CSCSegment::parameters ( void ) const

virtual

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 query::result, theLocalDirection, theOrigin, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties(), and 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

inlinevirtual

Covariance matrix of parameters()

Implements TrackingRecHit.

Definition at line 48 of file CSCSegment.h.

References theCovMatrix.

Referenced by CSCSegAlgoST::dumpSegment(), and StateSegmentMatcher::StateSegmentMatcher().

48 { return theCovMatrix; }

CSCSegment::theCovMatrix

AlgebraicSymMatrix theCovMatrix

Definition: CSCSegment.h:96

void CSCSegment::print ( void ) const

Definition at line 186 of file CSCSegment.cc.

References gather_cfg::cout.

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

AlgebraicMatrix CSCSegment::projectionMatrix ( ) const

virtual

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

Implements TrackingRecHit.

Definition at line 82 of file CSCSegment.cc.

References theProjectionMatrix.

                                                    {
   return theProjectionMatrix;
 }

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

virtual

Access to component RecHits (if any)

Implements TrackingRecHit.

Definition at line 30 of file CSCSegment.cc.

References theCSCRecHits.

Referenced by MuonAssociatorByHits::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 ( )

virtual

Non-const access to component RecHits (if any)

Implements TrackingRecHit.

Definition at line 38 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 86 of file CSCSegment.cc.

References 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 133 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 138 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 65 of file CSCSegment.h.

References theCSCRecHits.

Referenced by CSCSegAlgoST::dumpSegment(), operator<<(), sharesRecHits(), and DynamicTruncation::updateWithCSChits().

65 { 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 95 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 149 of file CSCSegment.cc.

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

                              {
   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;
 }