#include <DTSegmentCand.h>

Inheritance diagram for DTSegmentCand:

Classes
struct	AssPointLessZ

Public Types
typedef std::pair< std::shared_ptr< DTHitPairForFit >, DTEnums::DTCellSide >	AssPoint

typedef std::set< AssPoint, AssPointLessZ >	AssPointCont

Public Member Functions
virtual void	add (AssPoint newHit)
	add hits to the hit list. More...

virtual void	add (std::shared_ptr< DTHitPairForFit > hit, DTEnums::DTCellSide code)

virtual double	chi2 () const
	the chi2 (NOT chi2/NDOF) of the fit More...

virtual double	chi2ndof () const
	the chi2/NDOF of the fit More...

virtual AssPointCont	conflictingHitPairs (const DTSegmentCand &seg) const

virtual AlgebraicSymMatrix	covMatrix () const
	the covariance matrix More...

virtual LocalVector	direction () const

	DTSegmentCand (AssPointCont &hits, const DTSuperLayer *sl)
	Constructor. More...

	DTSegmentCand (const AssPointCont &hits, LocalPoint &position, LocalVector &direction, double chi2, const AlgebraicSymMatrix &covMat, const DTSuperLayer *sl)

virtual bool	good () const

virtual const AssPointCont &	hits () const
	the used hits More...

virtual bool	hitsShareLayer () const

virtual unsigned int	NDOF () const

virtual unsigned int	nHits () const

virtual int	nLayers () const
	number of different layers with hits More...

virtual int	nSharedHitPairs (const DTSegmentCand &seg) const
	number of shared hit pair with other segment candidate More...

	operator DTChamberRecSegment2D * () const
	convert this DTSegmentCand into a DTChamberRecSegment2D More...

	operator DTSLRecSegment2D * () const
	convert this DTSegmentCand into a DTRecSegment2D More...

virtual bool	operator< (const DTSegmentCand &seg)
	less operator based on nHits and chi2 More...

virtual bool	operator== (const DTSegmentCand &seg)
	equality operator based on position, direction, chi2 and nHits More...

virtual LocalPoint	position () const

virtual void	removeHit (AssPoint hit)
	remove hit from the candidate More...

virtual void	setChi2 (double &chi2)
	set chi2 More...

virtual void	setCovMatrix (AlgebraicSymMatrix &cov)
	set the cov matrix More...

virtual void	setDirection (LocalVector &dir)
	set direction More...

virtual void	setPosition (LocalPoint &pos)
	set position More...

virtual void	sett0 (double &t0)
	set t0 More...

const DTSuperLayer *	superLayer () const
	the super layer on which relies More...

virtual double	t0 () const
	the t0 of the segment More...

virtual	~DTSegmentCand ()
	Destructor. More...

Static Protected Attributes
static const double	chi2max = 20.

static const unsigned int	nHitsMin = 3

Private Attributes
double	theChi2

AlgebraicSymMatrix	theCovMatrix

LocalVector	theDirection

AssPointCont	theHits

LocalPoint	thePosition

const DTSuperLayer *	theSL

double	thet0

Detailed Description

A Candidate for a DT segment. It's used by the algorithm to build segments and store relative information. It must be transformed into a DTSegment for further use.

Author: Stefano Lacaprara - INFN Legnaro stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t; Riccardo Bellan - INFN TO ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 34 of file DTSegmentCand.h.

Member Typedef Documentation

◆ AssPoint

typedef std::pair<std::shared_ptr<DTHitPairForFit>, DTEnums::DTCellSide> DTSegmentCand::AssPoint

Definition at line 36 of file DTSegmentCand.h.

◆ AssPointCont

typedef std::set<AssPoint, AssPointLessZ> DTSegmentCand::AssPointCont

Definition at line 38 of file DTSegmentCand.h.

Constructor & Destructor Documentation

◆ DTSegmentCand() [1/2]

DTSegmentCand::DTSegmentCand	(	AssPointCont &	hits,
		const DTSuperLayer *	sl
	)

Constructor.

Definition at line 23 of file DTSegmentCand.cc.

23 : theSL(sl), theChi2(-1.), theHits(hits) {}

◆ DTSegmentCand() [2/2]

DTSegmentCand::DTSegmentCand	(	const AssPointCont &	hits,
		LocalPoint &	position,
		LocalVector &	direction,
		double	chi2,
		const AlgebraicSymMatrix &	covMat,
		const DTSuperLayer *	sl
	)

Definition at line 25 of file DTSegmentCand.cc.

31 : theSL(sl), thePosition(position), theDirection(direction), theChi2(chi2), theCovMatrix(covMat), theHits(hits) {}

◆ ~DTSegmentCand()

DTSegmentCand::~DTSegmentCand ( )

virtual

Destructor.

Definition at line 34 of file DTSegmentCand.cc.

34 {}

Member Function Documentation

◆ add() [1/2]

void DTSegmentCand::add ( AssPoint newHit )

virtual

add hits to the hit list.

Definition at line 58 of file DTSegmentCand.cc.

58 { theHits.insert(newHit); }

References theHits.

Referenced by DTMeantimerPatternReco::addHits(), and counter.Counter::register().

◆ add() [2/2]

void DTSegmentCand::add	(	std::shared_ptr< DTHitPairForFit >	hit,
		DTEnums::DTCellSide	code
	)

virtual

Definition at line 60 of file DTSegmentCand.cc.

                                                                                   {
   AssPoint newHit(hit, code);
   theHits.insert(newHit);
 }

References theHits.

Referenced by counter.Counter::register().

◆ chi2()

virtual double DTSegmentCand::chi2 ( void ) const

inlinevirtual

the chi2 (NOT chi2/NDOF) of the fit

Definition at line 61 of file DTSegmentCand.h.

61 { return theChi2; }

References theChi2.

Referenced by DTMeantimerPatternReco::addHits(), DTMeantimerPatternReco::fitWithT0(), good(), operator<(), operator<<(), operator==(), and setChi2().

◆ chi2ndof()

virtual double DTSegmentCand::chi2ndof ( ) const

inlinevirtual

the chi2/NDOF of the fit

Definition at line 64 of file DTSegmentCand.h.

64 { return theChi2 / (nHits() - 2.); }

References nHits(), and theChi2.

Referenced by DTMeantimerPatternReco::checkDoubleCandidates().

◆ conflictingHitPairs()

DTSegmentCand::AssPointCont DTSegmentCand::conflictingHitPairs ( const DTSegmentCand & seg ) const

virtual

return the hits shared with other segment and with confliction L/R assignment

Definition at line 82 of file DTSegmentCand.cc.

                                                                                            {
   AssPointCont result;
   const AssPointCont& hits2 = seg.theHits;
  
   //  if (nSharedHitPairs(seg)==0) return result;
  
   AssPointCont::const_iterator hitBegin2 = hits2.begin(), hitEnd2 = hits2.end();
   for (AssPointCont::const_iterator hit = theHits.begin(), hitEnd = theHits.end(); hit != hitEnd; ++hit) {
     for (AssPointCont::const_iterator hit2 = hitBegin2; hit2 != hitEnd2; ++hit2) {
       if ((*(*hit).first) == (*(*hit2).first) && (*hit).second != (*hit2).second) {
         result.insert(*hit);
         continue;
       }
     }
   }
   return result;
 }

References mps_fire::result, and theHits.

◆ covMatrix()

virtual AlgebraicSymMatrix DTSegmentCand::covMatrix ( ) const

inlinevirtual

the covariance matrix

Definition at line 85 of file DTSegmentCand.h.

85 { return theCovMatrix; }

References theCovMatrix.

◆ direction()

virtual LocalVector DTSegmentCand::direction ( ) const

inlinevirtual

Definition at line 82 of file DTSegmentCand.h.

82 { return theDirection; }

References theDirection.

Referenced by operator<<(), and operator==().

◆ good()

bool DTSegmentCand::good ( ) const

virtual

Reimplemented in DTSegmentExtendedCand.

Definition at line 100 of file DTSegmentCand.cc.

                                {
   // std::cout << NDOF() << "  " << chi2()/NDOF() << "   " << nHits() << std::endl;
   if (NDOF() == 0)
     return false;
   if (chi2() / NDOF() > chi2max || nHits() < nHitsMin)
     return false;
  
   if (nHits() == nHitsMin && hitsShareLayer())
     return false;
  
   return true;
 }

References chi2(), chi2max, hitsShareLayer(), NDOF(), nHits(), and nHitsMin.

Referenced by DTMeantimerPatternReco::addHits(), and DTCombinatorialPatternReco::buildSegments().

◆ hits()

virtual const AssPointCont& DTSegmentCand::hits ( void ) const

inlinevirtual

the used hits

Definition at line 122 of file DTSegmentCand.h.

122 { return theHits; }

References theHits.

Referenced by DTSegmentUpdator::fit(), and nSharedHitPairs().

◆ hitsShareLayer()

bool DTSegmentCand::hitsShareLayer ( ) const

virtual

Definition at line 113 of file DTSegmentCand.cc.

                                          {
   const unsigned int hitsSize = theHits.size();
   // we don't expect so many 1D hits, if such a segment arrives just drop it
   if (hitsSize > 20)
     return false;
  
   int layerN[hitsSize];
   unsigned int i = 0;
   for (DTSegmentCand::AssPointCont::iterator assHit = theHits.begin(); assHit != theHits.end(); ++assHit) {
     layerN[i] = (*assHit).first->id().layerId().layer() + 10 * (*assHit).first->id().superlayerId().superlayer();
     for (unsigned int j = 0; j < i; j++) {
       if (layerN[i] == layerN[j])
         return true;
     }
     i++;
   }
  
   return false;
 }

References mps_fire::i, dqmiolumiharvest::j, and theHits.

Referenced by good().

◆ NDOF()

virtual unsigned int DTSegmentCand::NDOF ( ) const

inlinevirtual

Definition at line 87 of file DTSegmentCand.h.

87 { return nHits() - 2; }

References nHits().

Referenced by good().

◆ nHits()

virtual unsigned int DTSegmentCand::nHits ( void ) const

inlinevirtual

Reimplemented in DTSegmentExtendedCand.

Definition at line 58 of file DTSegmentCand.h.

58 { return theHits.size(); }

References theHits.

Referenced by DTMeantimerPatternReco::addHits(), DTMeantimerPatternReco::checkDoubleCandidates(), chi2ndof(), DTMeantimerPatternReco::fitWithT0(), DTSegmentExtendedCand::good(), good(), NDOF(), DTSegmentExtendedCand::nHits(), operator<(), and operator==().

◆ nLayers()

int DTSegmentCand::nLayers ( ) const

virtual

number of different layers with hits

Definition at line 133 of file DTSegmentCand.cc.

                                  {
   // TODO
   return 0;
 }

◆ nSharedHitPairs()

int DTSegmentCand::nSharedHitPairs ( const DTSegmentCand & seg ) const

virtual

number of shared hit pair with other segment candidate

Definition at line 67 of file DTSegmentCand.cc.

                                                                  {
   int result = 0;
  
   for (AssPointCont::const_iterator hit = theHits.begin(); hit != theHits.end(); ++hit) {
     for (AssPointCont::const_iterator hit2 = seg.hits().begin(); hit2 != seg.hits().end(); ++hit2) {
       //  if(result) return result ; // TODO, uncomm this line or move it in another func
       if ((*(*hit).first) == (*(*hit2).first)) {
         ++result;
         continue;
       }
     }
   }
   return result;
 }

References hits(), mps_fire::result, and theHits.

◆ operator DTChamberRecSegment2D *()

DTSegmentCand::operator DTChamberRecSegment2D * ( ) const

convert this DTSegmentCand into a DTChamberRecSegment2D

Definition at line 161 of file DTSegmentCand.cc.

                                                      {
   // input position and direction are in sl frame, while must be stored in
   // chamber one: so I have to extrapolate the position (along the direction) to
   // the chamber reference plane.
  
   LocalPoint posInCh = theSL->chamber()->toLocal(theSL->toGlobal(position()));
   LocalVector dirInCh = theSL->chamber()->toLocal(theSL->toGlobal(direction()));
  
   LocalPoint pos = posInCh + dirInCh * posInCh.z() / cos(dirInCh.theta());
  
   double seg2DChi2 = chi2();
   AlgebraicSymMatrix seg2DCovMatrix = covMatrix();
  
   std::vector<DTRecHit1D> hits1D;
   for (DTSegmentCand::AssPointCont::iterator assHit = theHits.begin(); assHit != theHits.end(); ++assHit) {
     GlobalPoint hitGlobalPos = theSL->toGlobal((*assHit).first->localPosition((*assHit).second));
  
     LocalPoint hitPosInLayer = theSL->chamber()
                                    ->superLayer((*assHit).first->id().superlayerId())
                                    ->layer((*assHit).first->id().layerId())
                                    ->toLocal(hitGlobalPos);
  
     DTRecHit1D hit(((*assHit).first)->id(),
                    (*assHit).second,
                    ((*assHit).first)->digiTime(),
                    hitPosInLayer,
                    ((*assHit).first)->localPositionError());
     hits1D.push_back(hit);
   }
  
   return new DTChamberRecSegment2D(theSL->chamber()->id(), pos, dirInCh, seg2DCovMatrix, seg2DChi2, hits1D);
  
   // chamber and Phi SLs' frame are oriented in the same way, only a transaltion,
   // so the covariance matrix is the same!
 }

References hltPixelTracks_cff::chi2, funct::cos(), position, bscTrigger_cfi::theHits, PV3DBase< T, PVType, FrameType >::theta(), toLocal(), and PV3DBase< T, PVType, FrameType >::z().

◆ operator DTSLRecSegment2D *()

DTSegmentCand::operator DTSLRecSegment2D * ( ) const

convert this DTSegmentCand into a DTRecSegment2D

Definition at line 138 of file DTSegmentCand.cc.

                                                 {
   LocalPoint seg2Dposition = position();
   LocalVector seg2DDirection = direction();
   double seg2DChi2 = chi2();
   AlgebraicSymMatrix seg2DCovMatrix = covMatrix();
  
   std::vector<DTRecHit1D> hits1D;
   for (DTSegmentCand::AssPointCont::iterator assHit = theHits.begin(); assHit != theHits.end(); ++assHit) {
     GlobalPoint hitGlobalPos = theSL->toGlobal((*assHit).first->localPosition((*assHit).second));
  
     LocalPoint hitPosInLayer = theSL->layer((*assHit).first->id().layerId())->toLocal(hitGlobalPos);
  
     DTRecHit1D hit(((*assHit).first)->id(),
                    (*assHit).second,
                    ((*assHit).first)->digiTime(),
                    hitPosInLayer,
                    ((*assHit).first)->localPositionError());
     hits1D.push_back(hit);
   }
  
   return new DTSLRecSegment2D(theSL->id(), seg2Dposition, seg2DDirection, seg2DCovMatrix, seg2DChi2, hits1D);
 }

References hltPixelTracks_cff::chi2, position, bscTrigger_cfi::theHits, and toLocal().

◆ operator<()

bool DTSegmentCand::operator< ( const DTSegmentCand & seg )

virtual

less operator based on nHits and chi2

Definition at line 52 of file DTSegmentCand.cc.

                                                       {
   if (nHits() == seg.nHits())
     return (chi2() > seg.chi2());
   return (nHits() < seg.nHits());
 }

References chi2(), and nHits().

◆ operator==()

bool DTSegmentCand::operator== ( const DTSegmentCand & seg )

virtual

equality operator based on position, direction, chi2 and nHits

Definition at line 37 of file DTSegmentCand.cc.

                                                        {
   static const double epsilon = 0.00001;
   if (nHits() != seg.nHits())
     return false;
   if (fabs(chi2() - seg.chi2()) > epsilon)
     return false;
   if (fabs(position().x() - seg.position().x()) > epsilon || fabs(position().y() - seg.position().y()) > epsilon ||
       fabs(position().z() - seg.position().z()) > epsilon)
     return false;
   if (fabs(direction().x() - seg.direction().x()) > epsilon || fabs(direction().y() - seg.direction().y()) > epsilon ||
       fabs(direction().z() - seg.direction().z()) > epsilon)
     return false;
   return true;
 }

References chi2(), direction(), geometryDiff::epsilon, nHits(), position(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ position()

virtual LocalPoint DTSegmentCand::position ( ) const

inlinevirtual

Definition at line 79 of file DTSegmentCand.h.

79 { return thePosition; }

References thePosition.

Referenced by operator<<(), and operator==().

◆ removeHit()

void DTSegmentCand::removeHit ( AssPoint hit )

virtual

remove hit from the candidate

Definition at line 65 of file DTSegmentCand.cc.

65 { theHits.erase(badHit); }

References theHits.

Referenced by DTMeantimerPatternReco::addHits(), and DTSegmentCleaner::solveConflict().

◆ setChi2()

virtual void DTSegmentCand::setChi2 ( double & chi2 )

inlinevirtual

set chi2

Definition at line 103 of file DTSegmentCand.h.

103 { theChi2 = chi2; }

References chi2(), and theChi2.

Referenced by DTSegmentUpdator::fit().

◆ setCovMatrix()

virtual void DTSegmentCand::setCovMatrix ( AlgebraicSymMatrix & cov )

inlinevirtual

set the cov matrix

Definition at line 116 of file DTSegmentCand.h.

116 { theCovMatrix = cov; }

References theCovMatrix.

Referenced by DTSegmentUpdator::fit().

◆ setDirection()

virtual void DTSegmentCand::setDirection ( LocalVector & dir )

inlinevirtual

set direction

Definition at line 93 of file DTSegmentCand.h.

93 { theDirection = dir; }

References DeadROC_duringRun::dir, and theDirection.

Referenced by DTSegmentUpdator::fit().

◆ setPosition()

virtual void DTSegmentCand::setPosition ( LocalPoint & pos )

inlinevirtual

set position

Definition at line 90 of file DTSegmentCand.h.

90 { thePosition = pos; }

References thePosition.

Referenced by DTSegmentUpdator::fit().

◆ sett0()

virtual void DTSegmentCand::sett0 ( double & t0 )

inlinevirtual

set t0

Definition at line 106 of file DTSegmentCand.h.

106 { thet0 = t0; }

References t0(), and thet0.

Referenced by DTSegmentUpdator::fit().

◆ superLayer()

const DTSuperLayer* DTSegmentCand::superLayer ( ) const

inline

the super layer on which relies

Definition at line 76 of file DTSegmentCand.h.

76 { return theSL; }

References theSL.

Referenced by DTCombinatorialExtendedPatternReco::extendCandidates().

◆ t0()

virtual double DTSegmentCand::t0 ( ) const

inlinevirtual

the t0 of the segment

Definition at line 67 of file DTSegmentCand.h.

67 { return thet0; }

References thet0.

Referenced by DTMeantimerPatternReco::addHits(), DTMeantimerPatternReco::fitWithT0(), operator<<(), and sett0().

Member Data Documentation

◆ chi2max

const double DTSegmentCand::chi2max = 20.

staticprotected

Definition at line 151 of file DTSegmentCand.h.

Referenced by good().

◆ nHitsMin

const unsigned int DTSegmentCand::nHitsMin = 3

staticprotected

Definition at line 152 of file DTSegmentCand.h.

Referenced by good().

◆ theChi2

double DTSegmentCand::theChi2

private

Definition at line 140 of file DTSegmentCand.h.

Referenced by chi2(), chi2ndof(), and setChi2().

◆ theCovMatrix

AlgebraicSymMatrix DTSegmentCand::theCovMatrix

private

mat[1][1]=sigma (dx/dz) mat[2][2]=sigma (x) mat[1][2]=cov(dx/dz,x)

Definition at line 146 of file DTSegmentCand.h.

Referenced by covMatrix(), and setCovMatrix().

◆ theDirection

LocalVector DTSegmentCand::theDirection

private

Definition at line 139 of file DTSegmentCand.h.

Referenced by direction(), and setDirection().

◆ theHits

AssPointCont DTSegmentCand::theHits

private

Definition at line 148 of file DTSegmentCand.h.

Referenced by add(), conflictingHitPairs(), hits(), hitsShareLayer(), nHits(), nSharedHitPairs(), and removeHit().

◆ thePosition

LocalPoint DTSegmentCand::thePosition

private

Definition at line 138 of file DTSegmentCand.h.

Referenced by position(), and setPosition().

◆ theSL

const DTSuperLayer* DTSegmentCand::theSL

private

Definition at line 137 of file DTSegmentCand.h.

Referenced by superLayer().

◆ thet0

double DTSegmentCand::thet0

private

Definition at line 141 of file DTSegmentCand.h.

Referenced by sett0(), and t0().

Classes

Public Types

Public Member Functions

Static Protected Attributes

Private Attributes

Detailed Description

Member Typedef Documentation

◆ AssPoint

◆ AssPointCont

Constructor & Destructor Documentation

◆ DTSegmentCand() [1/2]

◆ DTSegmentCand() [2/2]

◆ ~DTSegmentCand()

Member Function Documentation

◆ add() [1/2]

◆ add() [2/2]

◆ chi2()

◆ chi2ndof()

◆ conflictingHitPairs()

◆ covMatrix()

◆ direction()

◆ good()

◆ hits()

◆ hitsShareLayer()

◆ NDOF()

◆ nHits()

◆ nLayers()

◆ nSharedHitPairs()

◆ operator DTChamberRecSegment2D *()

◆ operator DTSLRecSegment2D *()

◆ operator<()

◆ operator==()

◆ position()

◆ removeHit()

◆ setChi2()

◆ setCovMatrix()

◆ setDirection()

◆ setPosition()

◆ sett0()

◆ superLayer()

◆ t0()

Member Data Documentation

◆ chi2max

◆ nHitsMin

◆ theChi2

◆ theCovMatrix

◆ theDirection

◆ theHits

◆ thePosition

◆ theSL

◆ thet0