#include <MatchDCWithLine.h>

Collaboration diagram for TrkDriftCircleMath::MatchDCWithLine:

Public Types
enum	MatchStrategy { Road, Pull }

Public Member Functions
	MatchDCWithLine ()=default

	MatchDCWithLine (const Line &l, double deltaCut, MatchStrategy strategy, double tubeRadius)

void	set (const Line &l, double deltaCut, MatchStrategy strategy, double tubeRadius)

void	setTubeRadius (double radius)

const DCOnTrackVec &	match (const DCVec &dcs)

const DCOnTrackVec &	match (const DCOnTrackVec &dcs, const HitSelection *selection=nullptr, bool recoverMdtHits=true, bool usePreciseErrors=false)

unsigned int	deltas () const

unsigned int	hitsOutOfTime () const

unsigned int	hitsOnTrack () const

unsigned int	closeHits () const

unsigned int	showerHits () const

unsigned int	passedTubes () const

unsigned int	hitsMl1 () const

unsigned int	hitsMl2 () const

unsigned int	matchDifference () const
	returns the number of DCOnTrack that have a different status after the match, returns 0 if used with DCVec More...

Private Member Functions
void	reset ()

void	matchDC (DCOnTrack &dc, double &res, double &dline, bool forceOnTrack=false, bool forceOffTrack=false, bool usePreciseErrors=false)

bool	select (double residual, double error) const

Private Attributes
ResidualWithLine	m_resWithLine {}

double	m_tubeRadius {0.}

double	m_closeByCut {0.}

double	m_showerCut {0.}

double	m_deltaCut {0.}

MatchStrategy	m_strategy {}

unsigned int	m_deltas {0}

unsigned int	m_outOfTimes {0}

unsigned int	m_onTracks {0}

unsigned int	m_passedTubes {0}

unsigned int	m_ml1 {0}

unsigned int	m_ml2 {0}

unsigned int	m_closeHits {0}

unsigned int	m_showerHits {0}

unsigned int	m_matchDifference {0}

DCOnTrackVec	m_dcOnTrackVec

Detailed Description

Definition at line 16 of file MatchDCWithLine.h.

Member Enumeration Documentation

◆ MatchStrategy

enum TrkDriftCircleMath::MatchDCWithLine::MatchStrategy

Enumerator
Road
Pull

Definition at line 18 of file MatchDCWithLine.h.

18 { Road, Pull };

Constructor & Destructor Documentation

◆ MatchDCWithLine() [1/2]

TrkDriftCircleMath::MatchDCWithLine::MatchDCWithLine ( )

default

◆ MatchDCWithLine() [2/2]

TrkDriftCircleMath::MatchDCWithLine::MatchDCWithLine	(	const Line &	l,
		double	deltaCut,
		MatchStrategy	strategy,
		double	tubeRadius
	)

inline

Definition at line 22 of file MatchDCWithLine.h.

                                                                                                    :
             m_resWithLine{l}, m_deltaCut{deltaCut}, m_strategy{strategy} {
             m_dcOnTrackVec.reserve(50);
             setTubeRadius(tubeRadius);
         }

Member Function Documentation

◆ closeHits()

unsigned int TrkDriftCircleMath::MatchDCWithLine::closeHits ( ) const

inline

Definition at line 48 of file MatchDCWithLine.h.

48 { return m_closeHits; }

◆ deltas()

unsigned int TrkDriftCircleMath::MatchDCWithLine::deltas ( ) const

inline

Definition at line 45 of file MatchDCWithLine.h.

45 { return m_deltas; }

◆ hitsMl1()

unsigned int TrkDriftCircleMath::MatchDCWithLine::hitsMl1 ( ) const

inline

Definition at line 51 of file MatchDCWithLine.h.

51 { return m_ml1; }

◆ hitsMl2()

unsigned int TrkDriftCircleMath::MatchDCWithLine::hitsMl2 ( ) const

inline

Definition at line 52 of file MatchDCWithLine.h.

52 { return m_ml2; }

◆ hitsOnTrack()

unsigned int TrkDriftCircleMath::MatchDCWithLine::hitsOnTrack ( ) const

inline

Definition at line 47 of file MatchDCWithLine.h.

47 { return m_onTracks; }

◆ hitsOutOfTime()

unsigned int TrkDriftCircleMath::MatchDCWithLine::hitsOutOfTime ( ) const

inline

Definition at line 46 of file MatchDCWithLine.h.

46 { return m_outOfTimes; }

◆ match() [1/2]

const DCOnTrackVec & TrkDriftCircleMath::MatchDCWithLine::match	(	const DCOnTrackVec &	dcs,
		const HitSelection *	selection = `nullptr`,
		bool	recoverMdtHits = `true`,
		bool	usePreciseErrors = `false`
	)

Definition at line 45 of file MatchDCWithLine.cxx.

                                                                 {
         // clear cached information
         reset();
  
         DCOnTrackCit it = dcs.begin();
         DCOnTrackCit it_end = dcs.end();
         unsigned int index = 0;
         for (; it != it_end; ++it, ++index) {
             double dline = m_resWithLine.residual(it->position());
             double rhit = it->r() < 0. ? -it->r() : it->r();
             if (dline < 0.) dline *= -1.;
  
             double res = dline - rhit;
  
             DCOnTrack dcOnTrack(*it, res);
             dcOnTrack.errorTrack(it->errorTrack());
  
             // check whether any DCOnTracks change state from or to OnTrack with respect to the old selection
             if ((dcOnTrack.state() == DCOnTrack::OnTrack || it->state() == DCOnTrack::OnTrack) && dcOnTrack.state() != it->state())
                 ++m_matchDifference;
  
             if (selection) {
                 if ((*selection)[index] == 0) {
                     matchDC(dcOnTrack, res, dline, true, false, usePrecise);
                 } else if (recoverMdtHits) {
                     matchDC(dcOnTrack, res, dline, false, false, usePrecise);
                 } else {
                     matchDC(dcOnTrack, res, dline, false, true, usePrecise);
                 }
             } else {
                 matchDC(dcOnTrack, res, dline, false, false, usePrecise);
             }
  
             if (dcOnTrack.state() != DCOnTrack::OutsideRoad) m_dcOnTrackVec.push_back(dcOnTrack);
         }
         return m_dcOnTrackVec;
     }

◆ match() [2/2]

const DCOnTrackVec & TrkDriftCircleMath::MatchDCWithLine::match ( const DCVec & dcs )

add check to remove hits in the same tube check whether

this is not the first tube
the tubes have the same id

Definition at line 9 of file MatchDCWithLine.cxx.

                                                                {
         // clear cached information
         reset();
         DCCit it = dcs.begin();
         DCCit it_end = dcs.end();
         for (; it != it_end; ++it) {
             double dline = m_resWithLine.residual(it->position());
             double rhit = it->r() < 0. ? -it->r() : it->r();
             if (dline < 0.) dline *= -1.;
  
             double res = dline - rhit;
  
             DCOnTrack dcOnTrack(*it, res);
  
             matchDC(dcOnTrack, res, dline);
  
             if (dcOnTrack.state() != DCOnTrack::OutsideRoad) {
                 if (!m_dcOnTrackVec.empty() && m_dcOnTrackVec.back().id() == it->id()) {
                     // take the hit with the smallest residual
                     if (fabs(m_dcOnTrackVec.back().residual()) > fabs(res)) {
                         // exchange the two hits
                         m_dcOnTrackVec.pop_back();
                         m_dcOnTrackVec.push_back(dcOnTrack);
                     }
                 } else {
                     m_dcOnTrackVec.push_back(dcOnTrack);
                 }
             }
         }
         return m_dcOnTrackVec;
     }

◆ matchDC()

void TrkDriftCircleMath::MatchDCWithLine::matchDC	(	DCOnTrack &	dc,
		double &	res,
		double &	dline,
		bool	forceOnTrack = `false`,
		bool	forceOffTrack = `false`,
		bool	usePreciseErrors = `false`
	)

private

Definition at line 84 of file MatchDCWithLine.cxx.

                                                    {
         // check whether hit within range
         double error = usePrecise ? dcOnTrack.drPrecise() : dcOnTrack.dr();
         if (forceOnTrack || (!forceOffTrack && (dcOnTrack.driftState() == DriftCircle::InTime && select(res, error)))) {
             ++m_onTracks;
             ++m_passedTubes;
             if (dcOnTrack.id().ml() == 0)
                 ++m_ml1;
             else
                 ++m_ml2;
             dcOnTrack.state(DCOnTrack::OnTrack);
         } else {
             if (dline < m_tubeRadius) {
                 if (res < 0. || dcOnTrack.driftState() == DriftCircle::LateHit) {
                     dcOnTrack.state(DCOnTrack::OutOfTime);
                     ++m_outOfTimes;
                 } else {
                     dcOnTrack.state(DCOnTrack::Delta);
                     ++m_deltas;
                 }
                 ++m_passedTubes;
             } else {
                 if (dline < m_closeByCut) {
                     dcOnTrack.state(DCOnTrack::CloseDC);
                     ++m_closeHits;
                 } else {
                     dcOnTrack.state(DCOnTrack::OutsideRoad);
                 }
                 if (dline < m_showerCut) { ++m_showerHits; }
             }
         }
     }

◆ matchDifference()

unsigned int TrkDriftCircleMath::MatchDCWithLine::matchDifference ( ) const

inline

returns the number of DCOnTrack that have a different status after the match, returns 0 if used with DCVec

Definition at line 55 of file MatchDCWithLine.h.

55 { return m_matchDifference; }

◆ passedTubes()

unsigned int TrkDriftCircleMath::MatchDCWithLine::passedTubes ( ) const

inline

Definition at line 50 of file MatchDCWithLine.h.

50 { return m_passedTubes; }

◆ reset()

void TrkDriftCircleMath::MatchDCWithLine::reset ( )

private

Definition at line 127 of file MatchDCWithLine.cxx.

                                 {
         // clear cache
         m_dcOnTrackVec.clear();
  
         m_outOfTimes = 0;
         m_deltas = 0;
         m_onTracks = 0;
         m_passedTubes = 0;
         m_ml1 = 0;
         m_ml2 = 0;
         m_closeHits = 0;
         m_showerHits = 0;
         m_matchDifference = 0;
     }

◆ select()

bool TrkDriftCircleMath::MatchDCWithLine::select	(	double	residual,
		double	error
	)		const

private

Definition at line 118 of file MatchDCWithLine.cxx.

                                                                     {
         if (residual < 0.) residual *= -1.;
         switch (m_strategy) {
             case Road: return residual < m_deltaCut;
             case Pull: return residual < m_deltaCut * error;
             default: return false;
         }
     }

◆ set()

void TrkDriftCircleMath::MatchDCWithLine::set	(	const Line &	l,
		double	deltaCut,
		MatchStrategy	strategy,
		double	tubeRadius
	)

inline

Definition at line 28 of file MatchDCWithLine.h.

                                                                                             {
             m_resWithLine.set(l);
             m_deltaCut = deltaCut;
             m_strategy = strategy;
             setTubeRadius(tubeRadius);
         }

◆ setTubeRadius()

void TrkDriftCircleMath::MatchDCWithLine::setTubeRadius ( double radius )

inline

Definition at line 35 of file MatchDCWithLine.h.

                                           {
             m_tubeRadius = radius;
             m_closeByCut = 2 * radius;
             m_showerCut = 4 * radius;
         }

◆ showerHits()

unsigned int TrkDriftCircleMath::MatchDCWithLine::showerHits ( ) const

inline

Definition at line 49 of file MatchDCWithLine.h.

49 { return m_showerHits; }

Member Data Documentation

◆ m_closeByCut

double TrkDriftCircleMath::MatchDCWithLine::m_closeByCut {0.}

private

Definition at line 65 of file MatchDCWithLine.h.

◆ m_closeHits

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_closeHits {0}

private

Definition at line 75 of file MatchDCWithLine.h.

◆ m_dcOnTrackVec

DCOnTrackVec TrkDriftCircleMath::MatchDCWithLine::m_dcOnTrackVec

private

Definition at line 80 of file MatchDCWithLine.h.

◆ m_deltaCut

double TrkDriftCircleMath::MatchDCWithLine::m_deltaCut {0.}

private

Definition at line 67 of file MatchDCWithLine.h.

◆ m_deltas

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_deltas {0}

private

Definition at line 69 of file MatchDCWithLine.h.

◆ m_matchDifference

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_matchDifference {0}

private

Definition at line 77 of file MatchDCWithLine.h.

◆ m_ml1

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_ml1 {0}

private

Definition at line 73 of file MatchDCWithLine.h.

◆ m_ml2

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_ml2 {0}

private

Definition at line 74 of file MatchDCWithLine.h.

◆ m_onTracks

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_onTracks {0}

private

Definition at line 71 of file MatchDCWithLine.h.

◆ m_outOfTimes

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_outOfTimes {0}

private

Definition at line 70 of file MatchDCWithLine.h.

◆ m_passedTubes

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_passedTubes {0}

private

Definition at line 72 of file MatchDCWithLine.h.

◆ m_resWithLine

ResidualWithLine TrkDriftCircleMath::MatchDCWithLine::m_resWithLine {}

private

Definition at line 63 of file MatchDCWithLine.h.

◆ m_showerCut

double TrkDriftCircleMath::MatchDCWithLine::m_showerCut {0.}

private

Definition at line 66 of file MatchDCWithLine.h.

◆ m_showerHits

unsigned int TrkDriftCircleMath::MatchDCWithLine::m_showerHits {0}

private

Definition at line 76 of file MatchDCWithLine.h.

◆ m_strategy

MatchStrategy TrkDriftCircleMath::MatchDCWithLine::m_strategy {}

private

Definition at line 68 of file MatchDCWithLine.h.

◆ m_tubeRadius

double TrkDriftCircleMath::MatchDCWithLine::m_tubeRadius {0.}

private

Definition at line 64 of file MatchDCWithLine.h.

The documentation for this class was generated from the following files:

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ MatchStrategy

Constructor & Destructor Documentation

◆ MatchDCWithLine() [1/2]

◆ MatchDCWithLine() [2/2]

Member Function Documentation

◆ closeHits()

◆ deltas()

◆ hitsMl1()

◆ hitsMl2()

◆ hitsOnTrack()

◆ hitsOutOfTime()

◆ match() [1/2]

◆ match() [2/2]

◆ matchDC()

◆ matchDifference()

◆ passedTubes()

◆ reset()

◆ select()

◆ set()

◆ setTubeRadius()

◆ showerHits()

Member Data Documentation

◆ m_closeByCut

◆ m_closeHits

◆ m_dcOnTrackVec

◆ m_deltaCut

◆ m_deltas

◆ m_matchDifference

◆ m_ml1

◆ m_ml2

◆ m_onTracks

◆ m_outOfTimes

◆ m_passedTubes

◆ m_resWithLine

◆ m_showerCut

◆ m_showerHits

◆ m_strategy

◆ m_tubeRadius