StateSegmentMatcher Class Reference

#include <StateSegmentMatcher.h>

Public Member Functions
	StateSegmentMatcher (TrajectoryStateOnSurface const &, CSCSegment const &, LocalError const &)
	StateSegmentMatcher (TrajectoryStateOnSurface const &, DTRecSegment4D const &, LocalError const &)
double	value ()

Private Member Functions
void	setAPE2d (LocalError const &apeLoc)
void	setAPE4d (LocalError const &apeLoc)

Private Attributes
AlgebraicSymMatrix44	ape
AlgebraicSymMatrix22	ape_2d
double	estValue
AlgebraicSymMatrix44	m1
AlgebraicSymMatrix22	m1_2d
AlgebraicSymMatrix44	m2
AlgebraicSymMatrix22	m2_2d
bool	match2D
AlgebraicVector4	v1
AlgebraicVector2	v1_2d
AlgebraicVector4	v2
AlgebraicVector2	v2_2d

Detailed Description

Definition at line 68 of file StateSegmentMatcher.h.

Constructor & Destructor Documentation

StateSegmentMatcher() [1/2]

StateSegmentMatcher::StateSegmentMatcher	(	TrajectoryStateOnSurface const &	tsos,
		CSCSegment const &	cscseg4d,
		LocalError const &	apeLoc
	)

Definition at line 76 of file StateSegmentMatcher.cc.

References Tsos4D::errorMatrix(), callgraph::m2, CSCSegment::parameters(), CSCSegment::parametersError(), and Tsos4D::paramVector().

                                                                   {
  setAPE4d(apeLoc);
  match2D = false;
  AlgebraicVector cscseg = cscseg4d.parameters();
  v1[0] = cscseg[0];
  v1[1] = cscseg[1];
  v1[2] = cscseg[2];
  v1[3] = cscseg[3];

  AlgebraicSymMatrix rhErr_vect = cscseg4d.parametersError();
  m1(0, 0) = rhErr_vect(1, 1);
  m1(0, 1) = rhErr_vect(1, 2);
  m1(0, 2) = rhErr_vect(1, 3);
  m1(0, 3) = rhErr_vect(1, 4);
  m1(1, 0) = rhErr_vect(2, 1);
  m1(1, 1) = rhErr_vect(2, 2);
  m1(1, 2) = rhErr_vect(2, 3);
  m1(1, 3) = rhErr_vect(2, 4);
  m1(2, 0) = rhErr_vect(3, 1);
  m1(2, 1) = rhErr_vect(3, 2);
  m1(2, 2) = rhErr_vect(3, 3);
  m1(2, 3) = rhErr_vect(3, 4);
  m1(3, 0) = rhErr_vect(4, 1);
  m1(3, 1) = rhErr_vect(4, 2);
  m1(3, 2) = rhErr_vect(4, 3);
  m1(3, 3) = rhErr_vect(4, 4);

  Tsos4D tsos4d = Tsos4D(tsos);
  v2 = tsos4d.paramVector();
  m2 = tsos4d.errorMatrix();
}

StateSegmentMatcher() [2/2]

StateSegmentMatcher::StateSegmentMatcher	(	TrajectoryStateOnSurface const &	tsos,
		DTRecSegment4D const &	dtseg4d,
		LocalError const &	apeLoc
	)

Definition at line 16 of file StateSegmentMatcher.cc.

References Tsos4D::errorMatrix(), Tsos2DPhi::errorMatrix(), Tsos2DZed::errorMatrix(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), callgraph::m2, DTRecSegment4D::parameters(), DTRecSegment4D::parametersError(), Tsos4D::paramVector(), Tsos2DPhi::paramVector(), and Tsos2DZed::paramVector().

                                                                   {
  if (dtseg4d.hasPhi() && dtseg4d.hasZed()) {
    setAPE4d(apeLoc);
    match2D = false;
    AlgebraicVector dtseg = dtseg4d.parameters();
    v1[0] = dtseg[0];
    v1[1] = dtseg[1];
    v1[2] = dtseg[2];
    v1[3] = dtseg[3];

    AlgebraicSymMatrix rhErr_vect = dtseg4d.parametersError();
    m1(0, 0) = rhErr_vect(1, 1);
    m1(0, 1) = rhErr_vect(1, 2);
    m1(0, 2) = rhErr_vect(1, 3);
    m1(0, 3) = rhErr_vect(1, 4);
    m1(1, 0) = rhErr_vect(2, 1);
    m1(1, 1) = rhErr_vect(2, 2);
    m1(1, 2) = rhErr_vect(2, 3);
    m1(1, 3) = rhErr_vect(2, 4);
    m1(2, 0) = rhErr_vect(3, 1);
    m1(2, 1) = rhErr_vect(3, 2);
    m1(2, 2) = rhErr_vect(3, 3);
    m1(2, 3) = rhErr_vect(3, 4);
    m1(3, 0) = rhErr_vect(4, 1);
    m1(3, 1) = rhErr_vect(4, 2);
    m1(3, 2) = rhErr_vect(4, 3);
    m1(3, 3) = rhErr_vect(4, 4);

    Tsos4D tsos4d = Tsos4D(tsos);
    v2 = tsos4d.paramVector();
    m2 = tsos4d.errorMatrix();
  } else {
    setAPE2d(apeLoc);
    match2D = true;
    AlgebraicVector dtseg = dtseg4d.parameters();
    v1_2d[0] = dtseg[0];
    v1_2d[1] = dtseg[1];

    AlgebraicSymMatrix rhErr_vect = dtseg4d.parametersError();
    m1_2d(0, 0) = rhErr_vect(1, 1);
    m1_2d(0, 1) = rhErr_vect(1, 2);
    m1_2d(1, 0) = rhErr_vect(2, 1);
    m1_2d(1, 1) = rhErr_vect(2, 2);

    if (dtseg4d.hasPhi()) {
      Tsos2DPhi tsos2d = Tsos2DPhi(tsos);
      v2_2d = tsos2d.paramVector();
      m2_2d = tsos2d.errorMatrix();
    }

    if (dtseg4d.hasZed()) {
      Tsos2DZed tsos2d = Tsos2DZed(tsos);
      v2_2d = tsos2d.paramVector();
      m2_2d = tsos2d.errorMatrix();
    }
  }
}

Member Function Documentation

setAPE2d()

void StateSegmentMatcher::setAPE2d ( LocalError const &  apeLoc )

inlineprivate

Definition at line 96 of file StateSegmentMatcher.h.

References ape_2d, and LocalError::xx().

                                          {
    ape_2d[0][0] = 0;            //sigma (dx/dz)
    ape_2d[1][1] = apeLoc.xx();  //sigma (x)
    ape_2d[0][1] = 0;            //cov(dx/dz,x)
  };

setAPE4d()

void StateSegmentMatcher::setAPE4d ( LocalError const &  apeLoc )

inlineprivate

Definition at line 87 of file StateSegmentMatcher.h.

References ape, LocalError::xx(), and LocalError::yy().

                                          {
    ape[0][0] = 0;            //sigma (dx/dz)
    ape[1][1] = 0;            //sigma (dy/dz)
    ape[2][2] = apeLoc.xx();  //sigma (x)
    ape[3][3] = apeLoc.yy();  //sigma (y)
    ape[0][2] = 0;            //cov(dx/dz,x)
    ape[1][3] = 0;            //cov(dy/dz,y)
  };

value()

double StateSegmentMatcher::value

(

void

)

Definition at line 110 of file StateSegmentMatcher.cc.

References callgraph::m2.

Referenced by DynamicTruncation::testCSCstation(), and DynamicTruncation::testDTstation().

                                  {
  if (match2D) {
    AlgebraicVector2 v3(v1_2d - v2_2d);
    AlgebraicSymMatrix22 m3(m1_2d + m2_2d + ape_2d);
    bool m3i = !m3.Invert();
    if (m3i) {
      return 1e7;
    } else {
      estValue = ROOT::Math::Similarity(v3, m3);
    }

  } else {
    AlgebraicVector4 v3(v1 - v2);
    AlgebraicSymMatrix44 m3(m1 + m2 + ape);
    bool m3i = !m3.Invert();
    if (m3i) {
      return 1e7;
    } else {
      estValue = ROOT::Math::Similarity(v3, m3);
    }
  }
  return estValue;
}

Member Data Documentation

ape

AlgebraicSymMatrix44 StateSegmentMatcher::ape

private

Definition at line 81 of file StateSegmentMatcher.h.

Referenced by setAPE4d().

ape_2d

AlgebraicSymMatrix22 StateSegmentMatcher::ape_2d

private

Definition at line 83 of file StateSegmentMatcher.h.

Referenced by setAPE2d().

estValue

double StateSegmentMatcher::estValue

private

Definition at line 85 of file StateSegmentMatcher.h.

m1

AlgebraicSymMatrix44 StateSegmentMatcher::m1

private

Definition at line 81 of file StateSegmentMatcher.h.

m1_2d

AlgebraicSymMatrix22 StateSegmentMatcher::m1_2d

private

Definition at line 83 of file StateSegmentMatcher.h.

m2

AlgebraicSymMatrix44 StateSegmentMatcher::m2

private

Definition at line 81 of file StateSegmentMatcher.h.

m2_2d

AlgebraicSymMatrix22 StateSegmentMatcher::m2_2d

private

Definition at line 83 of file StateSegmentMatcher.h.

match2D

bool StateSegmentMatcher::match2D

private

Definition at line 84 of file StateSegmentMatcher.h.

v1

AlgebraicVector4 StateSegmentMatcher::v1

private

Definition at line 80 of file StateSegmentMatcher.h.

v1_2d

AlgebraicVector2 StateSegmentMatcher::v1_2d

private

Definition at line 82 of file StateSegmentMatcher.h.

v2

AlgebraicVector4 StateSegmentMatcher::v2

private

Definition at line 80 of file StateSegmentMatcher.h.

v2_2d

AlgebraicVector2 StateSegmentMatcher::v2_2d

private

Definition at line 82 of file StateSegmentMatcher.h.