DTSegmentCand Class Reference

#include <DTSegmentCand.h>

Inheritance diagram for DTSegmentCand:

Classes
struct	AssPointLessZ

Public Types
typedef std::pair < DTHitPairForFit *, DTEnums::DTCellSide >	AssPoint
typedef std::set< AssPoint, AssPointLessZ >	AssPointCont

Public Member Functions
virtual void	add (DTHitPairForFit *hit, DTEnums::DTCellSide code)
	add hits to the hit list. More...
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 (AssPointCont hits, LocalPoint &position, LocalVector &direction, double chi2, AlgebraicSymMatrix covMat, const DTSuperLayer *sl)
virtual bool	good () const
virtual 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...
const DTSuperLayer *	superLayer () const
	the super layer on which relies More...
virtual	~DTSegmentCand ()
	Destructor. More...

Static Protected Attributes
static double	chi2max =20.
static unsigned int	nHitsMin =3

Private Attributes
double	theChi2
AlgebraicSymMatrix	theCovMatrix
LocalVector	theDirection
AssPointCont	theHits
LocalPoint	thePosition
const DTSuperLayer *	theSL

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.

Date:

2010/01/21 15:33:32

Revision:

1.9

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

Member Typedef Documentation

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

Definition at line 39 of file DTSegmentCand.h.

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

Definition at line 41 of file DTSegmentCand.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 25 of file DTSegmentCand.cc.

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

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

Definition at line 30 of file DTSegmentCand.cc.

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

DTSegmentCand::~DTSegmentCand ( )

virtual

Destructor.

Definition at line 41 of file DTSegmentCand.cc.

                              {
}

Member Function Documentation

void DTSegmentCand::add ( DTHitPairForFit *  hit, DTEnums::DTCellSide  code  )

virtual

add hits to the hit list.

Definition at line 63 of file DTSegmentCand.cc.

References theHits.

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

virtual double DTSegmentCand::chi2 ( void  ) const

inlinevirtual

the chi2 (NOT chi2/NDOF) of the fit

Definition at line 65 of file DTSegmentCand.h.

References theChi2.

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

{return theChi2; }

virtual double DTSegmentCand::chi2ndof ( ) const

inlinevirtual

the chi2/NDOF of the fit

Definition at line 68 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 90 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 86 of file DTSegmentCand.h.

References theCovMatrix.

{return theCovMatrix; }

virtual LocalVector DTSegmentCand::direction ( ) const

inlinevirtual

Definition at line 83 of file DTSegmentCand.h.

References theDirection.

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

{ return theDirection;}

bool DTSegmentCand::good ( ) const

virtual

Reimplemented in DTSegmentExtendedCand.

Definition at line 110 of file DTSegmentCand.cc.

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

Referenced by DTCombinatorialPatternReco::buildSegments(), and DTSegmentUpdator::fit().

{
  if(NDOF() == 0) return false;
  if(chi2()/NDOF() > chi2max || nHits() < nHitsMin) return false;

  if(nHits() == nHitsMin && hitsShareLayer()) return false ;

  return true;
}

virtual AssPointCont DTSegmentCand::hits ( ) const

inlinevirtual

the used hits

Definition at line 119 of file DTSegmentCand.h.

References theHits.

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

{ return theHits;}

bool DTSegmentCand::hitsShareLayer ( ) const

virtual

Definition at line 120 of file DTSegmentCand.cc.

References i, j, and theHits.

Referenced by good().

{
  std::vector<int> layerN;

  for(DTSegmentCand::AssPointCont::iterator assHit=theHits.begin();
      assHit!=theHits.end(); ++assHit) {
    layerN.push_back((*assHit).first->id().layerId().layer());

    //std::cout << (*assHit).first->id().layerId().layer() << std::endl;
  }

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

  return false;
}

virtual unsigned int DTSegmentCand::NDOF ( ) const

inlinevirtual

Definition at line 88 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 62 of file DTSegmentCand.h.

References theHits.

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

{ return theHits.size(); }

int DTSegmentCand::nLayers ( ) const

virtual

number of different layers with hits

Definition at line 140 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 71 of file DTSegmentCand.cc.

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

                                                                {
  int result=0;
  AssPointCont hitsCont = seg.hits();

  for (AssPointCont::const_iterator hit=theHits.begin(); 
       hit!=theHits.end() ; ++hit) {
    for (AssPointCont::const_iterator hit2=hitsCont.begin();
         hit2!=hitsCont.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 175 of file DTSegmentCand.cc.

References funct::cos(), pos, 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 145 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 58 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 45 of file DTSegmentCand.cc.

References chi2(), direction(), 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 80 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 67 of file DTSegmentCand.cc.

References theHits.

Referenced by DTSegmentCleaner::solveConflict().

                                             {
  theHits.erase(badHit);
}

virtual void DTSegmentCand::setChi2 ( double &  chi2 )

inlinevirtual

set chi2

Definition at line 103 of file DTSegmentCand.h.

References chi2(), and theChi2.

Referenced by DTMeantimerPatternReco::addHits(), and DTSegmentUpdator::fit().

{ theChi2 = chi2 ;}

virtual void DTSegmentCand::setCovMatrix ( AlgebraicSymMatrix &  cov )

inlinevirtual

set the cov matrix

Definition at line 113 of file DTSegmentCand.h.

References theCovMatrix.

Referenced by DTSegmentUpdator::fit().

{ theCovMatrix = cov; }

virtual void DTSegmentCand::setDirection ( LocalVector &  dir )

inlinevirtual

set direction

Definition at line 94 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 91 of file DTSegmentCand.h.

References pos, and thePosition.

Referenced by DTSegmentUpdator::fit().

{ thePosition=pos; }

const DTSuperLayer* DTSegmentCand::superLayer ( ) const

inline

the super layer on which relies

Definition at line 77 of file DTSegmentCand.h.

References theSL.

Referenced by DTCombinatorialExtendedPatternReco::extendCandidates().

{return theSL;}

Member Data Documentation

double DTSegmentCand::chi2max =20.

staticprotected

Definition at line 151 of file DTSegmentCand.h.

Referenced by good().

unsigned int DTSegmentCand::nHitsMin =3

staticprotected

Definition at line 152 of file DTSegmentCand.h.

Referenced by good().

double DTSegmentCand::theChi2

private

Definition at line 140 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 145 of file DTSegmentCand.h.

Referenced by covMatrix(), and setCovMatrix().

LocalVector DTSegmentCand::theDirection

private

Definition at line 139 of file DTSegmentCand.h.

Referenced by direction(), and setDirection().

AssPointCont DTSegmentCand::theHits

private

Definition at line 148 of file DTSegmentCand.h.

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

LocalPoint DTSegmentCand::thePosition

private

Definition at line 138 of file DTSegmentCand.h.

Referenced by position(), and setPosition().

const DTSuperLayer* DTSegmentCand::theSL

private

Definition at line 137 of file DTSegmentCand.h.

Referenced by superLayer().