DTSegmentCand Class Reference

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

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

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

{}

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.

References theHits.

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

{ theHits.insert(newHit); }

add() [2/2]

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

virtual

Definition at line 60 of file DTSegmentCand.cc.

References theHits.

Referenced by counter.Counter::register().

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

chi2()

virtual double DTSegmentCand::chi2 ( void  ) const

inlinevirtual

the chi2 (NOT chi2/NDOF) of the fit

Definition at line 61 of file DTSegmentCand.h.

References theChi2.

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

{ return theChi2; }

chi2ndof()

virtual double DTSegmentCand::chi2ndof ( ) const

inlinevirtual

the chi2/NDOF of the fit

Definition at line 64 of file DTSegmentCand.h.

References nHits(), and theChi2.

Referenced by DTMeantimerPatternReco::checkDoubleCandidates().

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

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.

References mps_fire::result, and theHits.

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

covMatrix()

virtual AlgebraicSymMatrix DTSegmentCand::covMatrix ( ) const

inlinevirtual

the covariance matrix

Definition at line 85 of file DTSegmentCand.h.

References theCovMatrix.

{ return theCovMatrix; }

direction()

virtual LocalVector DTSegmentCand::direction ( ) const

inlinevirtual

Definition at line 82 of file DTSegmentCand.h.

References theDirection.

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

{ return theDirection; }

good()

bool DTSegmentCand::good ( ) const

virtual

Reimplemented in DTSegmentExtendedCand.

Definition at line 100 of file DTSegmentCand.cc.

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

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

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

hits()

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

inlinevirtual

the used hits

Definition at line 122 of file DTSegmentCand.h.

References theHits.

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

{ return theHits; }

hitsShareLayer()

bool DTSegmentCand::hitsShareLayer ( ) const

virtual

Definition at line 113 of file DTSegmentCand.cc.

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

Referenced by good().

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

NDOF()

virtual unsigned int DTSegmentCand::NDOF ( ) const

inlinevirtual

Definition at line 87 of file DTSegmentCand.h.

References nHits().

Referenced by good().

{ return nHits() - 2; }

nHits()

virtual unsigned int DTSegmentCand::nHits ( void  ) const

inlinevirtual

Reimplemented in DTSegmentExtendedCand.

Definition at line 58 of file DTSegmentCand.h.

References theHits.

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

{ return theHits.size(); }

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.

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

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

operator DTChamberRecSegment2D *()

DTSegmentCand::operator DTChamberRecSegment2D *

(

)

const

convert this DTSegmentCand into a DTChamberRecSegment2D

Definition at line 161 of file DTSegmentCand.cc.

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

                                                     {
  // 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!
}

operator DTSLRecSegment2D *()

DTSegmentCand::operator DTSLRecSegment2D *

(

)

const

convert this DTSegmentCand into a DTRecSegment2D

Definition at line 138 of file DTSegmentCand.cc.

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

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

operator<()

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

virtual

less operator based on nHits and chi2

Definition at line 52 of file DTSegmentCand.cc.

References chi2(), and nHits().

                                                      {
  if (nHits() == seg.nHits())
    return (chi2() > seg.chi2());
  return (nHits() < seg.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.

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

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

position()

virtual LocalPoint DTSegmentCand::position ( ) const

inlinevirtual

Definition at line 79 of file DTSegmentCand.h.

References thePosition.

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

{ return thePosition; }

removeHit()

void DTSegmentCand::removeHit ( AssPoint  hit )

virtual

remove hit from the candidate

Definition at line 65 of file DTSegmentCand.cc.

References theHits.

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

{ theHits.erase(badHit); }

setChi2()

virtual void DTSegmentCand::setChi2 ( double &  chi2 )

inlinevirtual

set chi2

Definition at line 103 of file DTSegmentCand.h.

References chi2(), and theChi2.

Referenced by DTSegmentUpdator::fit().

{ theChi2 = chi2; }

setCovMatrix()

virtual void DTSegmentCand::setCovMatrix ( AlgebraicSymMatrix &  cov )

inlinevirtual

set the cov matrix

Definition at line 116 of file DTSegmentCand.h.

References theCovMatrix.

Referenced by DTSegmentUpdator::fit().

{ theCovMatrix = cov; }

setDirection()

virtual void DTSegmentCand::setDirection ( LocalVector &  dir )

inlinevirtual

set direction

Definition at line 93 of file DTSegmentCand.h.

References DeadROC_duringRun::dir, and theDirection.

Referenced by DTSegmentUpdator::fit().

{ theDirection = dir; }

setPosition()

virtual void DTSegmentCand::setPosition ( LocalPoint &  pos )

inlinevirtual

set position

Definition at line 90 of file DTSegmentCand.h.

References thePosition.

Referenced by DTSegmentUpdator::fit().

{ thePosition = pos; }

sett0()

virtual void DTSegmentCand::sett0 ( double &  t0 )

inlinevirtual

set t0

Definition at line 106 of file DTSegmentCand.h.

References t0(), and thet0.

Referenced by DTSegmentUpdator::fit().

{ thet0 = t0; }

superLayer()

const DTSuperLayer* DTSegmentCand::superLayer ( ) const

inline

the super layer on which relies

Definition at line 76 of file DTSegmentCand.h.

References theSL.

Referenced by DTCombinatorialExtendedPatternReco::extendCandidates().

{ return theSL; }

t0()

virtual double DTSegmentCand::t0 ( ) const

inlinevirtual

the t0 of the segment

Definition at line 67 of file DTSegmentCand.h.

References thet0.

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

{ return thet0; }

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