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 35 of file DTSegmentCand.h.

Member Typedef Documentation

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

Definition at line 38 of file DTSegmentCand.h.

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

Definition at line 40 of file DTSegmentCand.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 23 of file DTSegmentCand.cc.

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

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

Definition at line 28 of file DTSegmentCand.cc.

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

DTSegmentCand::~DTSegmentCand ( )

virtual

Destructor.

Definition at line 39 of file DTSegmentCand.cc.

                              {
}

Member Function Documentation

void DTSegmentCand::add ( AssPoint  newHit )

virtual

add hits to the hit list.

Definition at line 61 of file DTSegmentCand.cc.

References theHits.

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

                                       {
  theHits.insert(newHit);
}

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

virtual

Definition at line 65 of file DTSegmentCand.cc.

References theHits.

Referenced by counter.Counter::register().

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

virtual double DTSegmentCand::chi2 ( void  ) const

inlinevirtual

the chi2 (NOT chi2/NDOF) of the fit

Definition at line 64 of file DTSegmentCand.h.

References theChi2.

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

{ return theChi2; }

virtual double DTSegmentCand::chi2ndof ( ) const

inlinevirtual

the chi2/NDOF of the fit

Definition at line 67 of file DTSegmentCand.h.

References nHits(), and theChi2.

Referenced by DTMeantimerPatternReco::checkDoubleCandidates().

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

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 92 of file DTSegmentCand.cc.

References query::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;
}

virtual AlgebraicSymMatrix DTSegmentCand::covMatrix ( ) const

inlinevirtual

the covariance matrix

Definition at line 88 of file DTSegmentCand.h.

References theCovMatrix.

{return theCovMatrix; }

virtual LocalVector DTSegmentCand::direction ( ) const

inlinevirtual

Definition at line 85 of file DTSegmentCand.h.

References theDirection.

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

{ return theDirection;}

bool DTSegmentCand::good ( ) const

virtual

Reimplemented in DTSegmentExtendedCand.

Definition at line 112 of file DTSegmentCand.cc.

References chi2(), chi2max, edm::false, 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;
}

virtual const AssPointCont& DTSegmentCand::hits ( ) const

inlinevirtual

the used hits

Definition at line 125 of file DTSegmentCand.h.

References theHits.

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

{ return theHits;}

bool DTSegmentCand::hitsShareLayer ( ) const

virtual

Definition at line 123 of file DTSegmentCand.cc.

References i, j, and theHits.

Referenced by good().

{
  int layerN[20];
  int i=0;
  
  // we don't expect so many 1D hits, if such a segment arrives just drop it
  if (theHits.size()>20) return false;

  for(DTSegmentCand::AssPointCont::iterator assHit=theHits.begin();
      assHit!=theHits.end(); ++assHit) {
    layerN[i]=(*assHit).first->id().layerId().layer()+10*(*assHit).first->id().superlayerId().superlayer();
    i++;
  }

  for(int i=0;i<(int)theHits.size();i++){
    for(int j=0;j<i;j++){
      if(layerN[i] == layerN[j]) return true;
    }
  }

  return false;
}

virtual unsigned int DTSegmentCand::NDOF ( ) const

inlinevirtual

Definition at line 90 of file DTSegmentCand.h.

References nHits().

Referenced by good().

{ return nHits()-2; }

virtual unsigned int DTSegmentCand::nHits ( void  ) const

inlinevirtual

Reimplemented in DTSegmentExtendedCand.

Definition at line 61 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(); }

int DTSegmentCand::nLayers ( ) const

virtual

number of different layers with hits

Definition at line 146 of file DTSegmentCand.cc.

                                 {
  // TODO
  return 0;
}

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

virtual

number of shared hit pair with other segment candidate

Definition at line 74 of file DTSegmentCand.cc.

References hits(), query::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;
}

DTSegmentCand::operator DTChamberRecSegment2D *

(

)

const

convert this DTSegmentCand into a DTChamberRecSegment2D

Definition at line 181 of file DTSegmentCand.cc.

References funct::cos(), position, 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!
}

DTSegmentCand::operator DTSLRecSegment2D *

(

)

const

convert this DTSegmentCand into a DTRecSegment2D

Definition at line 151 of file DTSegmentCand.cc.

References position, 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);
}

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

virtual

less operator based on nHits and chi2

Definition at line 56 of file DTSegmentCand.cc.

References chi2(), and nHits().

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

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

virtual

equality operator based on position, direction, chi2 and nHits

Definition at line 43 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;
}

virtual LocalPoint DTSegmentCand::position ( ) const

inlinevirtual

Definition at line 82 of file DTSegmentCand.h.

References thePosition.

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

{ return thePosition; }

void DTSegmentCand::removeHit ( AssPoint  hit )

virtual

remove hit from the candidate

Definition at line 70 of file DTSegmentCand.cc.

References theHits.

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

                                             {
  theHits.erase(badHit);
}

virtual void DTSegmentCand::setChi2 ( double &  chi2 )

inlinevirtual

set chi2

Definition at line 106 of file DTSegmentCand.h.

References chi2(), and theChi2.

Referenced by DTSegmentUpdator::fit().

{ theChi2 = chi2; }

virtual void DTSegmentCand::setCovMatrix ( AlgebraicSymMatrix &  cov )

inlinevirtual

set the cov matrix

Definition at line 119 of file DTSegmentCand.h.

References theCovMatrix.

Referenced by DTSegmentUpdator::fit().

{ theCovMatrix = cov; }

virtual void DTSegmentCand::setDirection ( LocalVector &  dir )

inlinevirtual

set direction

Definition at line 96 of file DTSegmentCand.h.

References dir, and theDirection.

Referenced by DTSegmentUpdator::fit().

{ theDirection = dir; }

virtual void DTSegmentCand::setPosition ( LocalPoint &  pos )

inlinevirtual

set position

Definition at line 93 of file DTSegmentCand.h.

References thePosition.

Referenced by DTSegmentUpdator::fit().

{ thePosition=pos; }

virtual void DTSegmentCand::sett0 ( double &  t0 )

inlinevirtual

set t0

Definition at line 109 of file DTSegmentCand.h.

References t0(), and thet0.

Referenced by DTSegmentUpdator::fit().

{ thet0 = t0; }

const DTSuperLayer* DTSegmentCand::superLayer ( ) const

inline

the super layer on which relies

Definition at line 79 of file DTSegmentCand.h.

References theSL.

Referenced by DTCombinatorialExtendedPatternReco::extendCandidates().

{return theSL;}

virtual double DTSegmentCand::t0 ( ) const

inlinevirtual

the t0 of the segment

Definition at line 70 of file DTSegmentCand.h.

References thet0.

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

{ return thet0; }

Member Data Documentation

const double DTSegmentCand::chi2max =20.

staticprotected

Definition at line 157 of file DTSegmentCand.h.

Referenced by good().

const unsigned int DTSegmentCand::nHitsMin =3

staticprotected

Definition at line 158 of file DTSegmentCand.h.

Referenced by good().

double DTSegmentCand::theChi2

private

Definition at line 146 of file DTSegmentCand.h.

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

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 152 of file DTSegmentCand.h.

Referenced by covMatrix(), and setCovMatrix().

LocalVector DTSegmentCand::theDirection

private

Definition at line 145 of file DTSegmentCand.h.

Referenced by direction(), and setDirection().

AssPointCont DTSegmentCand::theHits

private

Definition at line 154 of file DTSegmentCand.h.

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

LocalPoint DTSegmentCand::thePosition

private

Definition at line 144 of file DTSegmentCand.h.

Referenced by position(), and setPosition().

const DTSuperLayer* DTSegmentCand::theSL

private

Definition at line 143 of file DTSegmentCand.h.

Referenced by superLayer().

double DTSegmentCand::thet0

private

Definition at line 147 of file DTSegmentCand.h.

Referenced by sett0(), and t0().