InDet::TRT_TrajectoryElement_xk Class Reference

#include <TRT_TrajectoryElement_xk.h>

Collaboration diagram for InDet::TRT_TrajectoryElement_xk:

Public Member Functions
	TRT_TrajectoryElement_xk ()
	TRT_TrajectoryElement_xk (const TRT_TrajectoryElement_xk &)
	~TRT_TrajectoryElement_xk ()
TRT_TrajectoryElement_xk &	operator= (const TRT_TrajectoryElement_xk &)
const bool &	isCluster () const
const bool &	isBarrel () const
const int &	status () const
const int &	nlinks () const
const int &	bestlink () const
const double &	radius () const
const double &	z () const
const double &	radiusMin () const
const double &	radiusMax () const
const double &	dpositive () const
const double &	dnegative () const
const InDetDD::TRT_BaseElement *	detElement () const
const TRT_ExtensionDriftCircleLink_xk &	link (int i) const
void	set (const TRT_ID *, const Trk::IPatternParametersPropagator *, const Trk::IPatternParametersUpdator *, const Trk::IRIO_OnTrackCreator *, const Trk::IRIO_OnTrackCreator *, double)
void	set (const Trk::MagneticFieldProperties &, const AtlasFieldCacheCondObj *)
bool	initiateForPrecisionSeed (bool, const InDetDD::TRT_BaseElement *, InDet::TRT_DriftCircleCollection::const_iterator &, InDet::TRT_DriftCircleCollection::const_iterator &, std::pair< Amg::Vector3D, double > &, const double *, double)
bool	initiateForTRTSeed (bool, const InDetDD::TRT_BaseElement *, InDet::TRT_DriftCircleCollection::const_iterator &, InDet::TRT_DriftCircleCollection::const_iterator &, std::pair< Amg::Vector3D, double > &, const double *, double)
void	initiateLinksForPrecisionSeed (InDet::TRT_DriftCircleCollection::const_iterator &, InDet::TRT_DriftCircleCollection::const_iterator &, std::pair< Amg::Vector3D, double > &, const double *, double)
void	initiateLinksForTRTSeed (InDet::TRT_DriftCircleCollection::const_iterator &, InDet::TRT_DriftCircleCollection::const_iterator &, std::pair< Amg::Vector3D, double > &, const double *, double)
bool	boundaryTest (double, std::pair< Amg::Vector3D, double > &)
bool	buildForPrecisionSeed (double, double, bool &, bool &)
bool	buildForTRTSeed (double, double, bool &, bool &)
double	findCloseLink (double, double)
void	radiusCorrection (const double &)
const Trk::RIO_OnTrack *	rioOnTrack ()
std::unique_ptr< Trk::RIO_OnTrack >	rioOnTrackSimple () const
bool	trajectoryGlobalPosition (Amg::Vector3D &, double &)
bool	trackParametersEstimation (TRT_TrajectoryElement_xk *, TRT_TrajectoryElement_xk *, Trk::PatternTrackParameters &, double)
bool	trackParametersEstimation (TRT_TrajectoryElement_xk *, Trk::PatternTrackParameters &, double)
void	polarAngleEstimation (TRT_TrajectoryElement_xk *, Amg::Vector3D &, Amg::Vector3D &, double, double, double *) const
bool	addCluster (Trk::PatternTrackParameters &, Trk::PatternTrackParameters &, double &)
bool	propagate (Trk::PatternTrackParameters &, Trk::PatternTrackParameters &)

Protected Attributes
bool	m_barrel
bool	m_isCluster
int	m_status
int	m_bestlink
int	m_nlinks
double	m_z
double	m_zMin
double	m_zMax
double	m_radius
double	m_radiusMin
double	m_radiusMax
double	m_scale_error
double	m_dpositive
double	m_dnegative
TRT_ExtensionDriftCircleLink_xk	m_link [24]
const InDetDD::TRT_BaseElement *	m_detelement
const TRT_ID *	m_trtid
const Trk::IPatternParametersPropagator *	m_proptool
const Trk::IPatternParametersUpdator *	m_updatortool
const Trk::IRIO_OnTrackCreator *	m_riomakerD
const Trk::IRIO_OnTrackCreator *	m_riomakerN
Trk::MagneticFieldProperties	m_fieldprop
MagField::AtlasFieldCache	m_fieldCache

Detailed Description

Definition at line 35 of file TRT_TrajectoryElement_xk.h.

Constructor & Destructor Documentation

TRT_TrajectoryElement_xk() [1/2]

InDet::TRT_TrajectoryElement_xk::TRT_TrajectoryElement_xk ( )

inline

Definition at line 170 of file TRT_TrajectoryElement_xk.h.

    {
      m_isCluster   = false;
      m_status      =-1    ;
      m_riomakerD   = 0    ;
      m_riomakerN   = 0    ;
      m_proptool    = 0    ;
      m_updatortool = 0    ;
      m_trtid       = 0    ;
      m_scale_error = 2.   ;
      m_barrel      = false;
      m_bestlink    = 0    ;
      m_nlinks      = 0    ;
      m_z           = 0.   ;
      m_zMin        = 0.   ;
      m_zMax        = 0.   ;
      m_radius      = 0.   ;
      m_radiusMin   = 0.   ;
      m_radiusMax   = 0.   ;
      m_dpositive   = 0.   ;
      m_dnegative   = 0.   ;
      m_detelement  = 0    ;
    }

TRT_TrajectoryElement_xk() [2/2]

InDet::TRT_TrajectoryElement_xk::TRT_TrajectoryElement_xk ( const TRT_TrajectoryElement_xk & E )

inline

Definition at line 194 of file TRT_TrajectoryElement_xk.h.

    {
      (*this) = E;
    }

~TRT_TrajectoryElement_xk()

InDet::TRT_TrajectoryElement_xk::~TRT_TrajectoryElement_xk ( )

inline

Definition at line 228 of file TRT_TrajectoryElement_xk.h.

{}

Member Function Documentation

addCluster()

bool InDet::TRT_TrajectoryElement_xk::addCluster	(	Trk::PatternTrackParameters &	Ta,
		Trk::PatternTrackParameters &	Tb,
		double &	Xi2 )

Definition at line 614 of file TRT_TrajectoryElement_xk.cxx.

{
 
  if(m_status <=0) return false;
 
  int n;
  if(m_status==2) {
 
    const TRT_DriftCircle* DS = m_link[m_bestlink].cluster();
    double dr = DS->localPosition().x(); if(m_link[m_bestlink].impact() < 0.) dr=-dr;
 
    Trk::DefinedParameter radius(dr,Trk::locX);
    Trk::LocalParameters  lp(radius);
 
    Trk::PatternTrackParameters To = Ta;
    bool Q = m_updatortool->addToState(Ta,lp,DS->localCovariance(),Tb,Xi2,n);
    if(Q && Xi2 < 15.) return Q;
    m_status = 1; Ta = To;
  }
 
  Trk::DefinedParameter  radius(0.,Trk::locX);
  Trk::LocalParameters   lp(radius);
 
  AmgSymMatrix(1) cov; cov<<1.33333;
  return m_updatortool->addToState(Ta,lp,cov,Tb,Xi2,n);
}

bestlink()

const int & InDet::TRT_TrajectoryElement_xk::bestlink ( ) const

inline

Definition at line 52 of file TRT_TrajectoryElement_xk.h.

{return m_bestlink  ;}

boundaryTest()

bool InDet::TRT_TrajectoryElement_xk::boundaryTest	(	double	dw,
		std::pair< Amg::Vector3D, double > &	gp )

Definition at line 90 of file TRT_TrajectoryElement_xk.cxx.

{
  if (!m_detelement){
    return false;
  }
 
  double x = gp.first.x();
  double y = gp.first.y();
  double z = gp.first.z();
 
  const Amg::Vector3D& C = m_detelement->center();
 
  // Test track position
  //
  const Trk::SurfaceBounds& surfBounds = m_detelement->bounds();
  if (surfBounds.type() == Trk::SurfaceBounds::Rectangle) {
    const Trk::RectangleBounds* rb =
      static_cast<const Trk::RectangleBounds*>(&surfBounds);
    // Barrel
    //
    m_barrel = true;
    m_z = z;
    m_zMin = C.z() - rb->halflengthY();
    m_zMax = C.z() + rb->halflengthY();
    double d = std::abs(z - C.z());
    if (d > rb->halflengthY() + dw) {
      return false;
    }
  } else if (surfBounds.type() == Trk::SurfaceBounds::Disc) {
    const Trk::DiscBounds* db =
      static_cast<const Trk::DiscBounds*>(&surfBounds);
    // Endcap
    //
    m_barrel = false;
    m_radius = std::sqrt(x * x + y * y);
    m_z = z;
    m_radiusMin = db->rMin();
    m_radiusMax = db->rMax();
 
    double d = m_radius - m_radiusMin;
    if (d < -dw) {
      return false;
    }
    d = m_radiusMax - m_radius;
    if (d < -dw) {
      return false;
    }
  }
  return true;
}

buildForPrecisionSeed()

bool InDet::TRT_TrajectoryElement_xk::buildForPrecisionSeed	(	double	a,
		double	b,
		bool &	useDriftTime,
		bool &	hole )

Definition at line 455 of file TRT_TrajectoryElement_xk.cxx.

{
 
  m_bestlink = -1;
  m_status   = -1;
  hole = false; if(!m_nlinks) return false;
 
  for(int l=0; l!=m_nlinks; ++l) {
 
    double v  = (a*m_link[l].way()+b)*m_link[l].way();
    double d  = m_link[l].distance()-v;
    double ad = std::abs(d);
 
    if(ad > 2.05) continue;
 
    if(!m_link[l].cluster()) {
 
      if(ad < 1.85) {
    m_status = 0; hole=true; m_bestlink=l; m_link[l].newImpactParameter(d);
      }
      return false;
    }
 
    m_status  = 1; m_bestlink=l; m_link[l].newImpactParameter(d);
    if(!useDriftTime) return true;
 
    double r = m_link[l].cluster()->localPosition()[Trk::driftRadius];
    double e = m_scale_error*std::sqrt(m_link[l].cluster()->localCovariance()(0,0));
 
    if(r > 2.05) r = 2.05;
    double r2 = r+e; if(r2 > 2.05) r2 = 2.05;
    double r1 = r-e; if(r1 >= r2) r1 = r2-2.*e;
 
    if(ad < r1 || ad > r2) useDriftTime = false; else m_status = 2;
 
    return true;
  }
  return false;
}

buildForTRTSeed()

bool InDet::TRT_TrajectoryElement_xk::buildForTRTSeed	(	double	a,
		double	b,
		bool &	useDriftTime,
		bool &	hole )

Definition at line 501 of file TRT_TrajectoryElement_xk.cxx.

{
 
  m_status   = -1;
  hole = false; if( m_bestlink < 0) return false;
  int l = m_bestlink;
 
  double v  = (a*m_link[l].way()+b)*m_link[l].way();
  double d  = m_link[l].distance()-v;
  double ad = std::abs(d);
 
  if(ad > 2.05) return false;
 
  if(!m_link[l].cluster()) {
 
    if(ad < 1.85) {
      m_status = 0; hole=true; m_bestlink=l; m_link[l].newImpactParameter(d);
    }
    return false;
  }
  m_status  = 1; m_link[l].newImpactParameter(d);
  if(!useDriftTime) return true;
 
  double r = m_link[l].cluster()->localPosition()[Trk::driftRadius];
  double e = m_scale_error*std::sqrt(m_link[l].cluster()->localCovariance()(0,0));
 
  if(r > 2.05) r = 2.05;
  double r2 = r+e; if(r2 > 2.05) r2 = 2.05;
  double r1 = r-e; if(r1 >= r2) r1 = r2-2.*e;
  if(ad < r1 || ad > r2) useDriftTime = false; else m_status = 2;
  return true;
}

detElement()

const InDetDD::TRT_BaseElement * InDet::TRT_TrajectoryElement_xk::detElement ( ) const

inline

Definition at line 59 of file TRT_TrajectoryElement_xk.h.

{return m_detelement;}

dnegative()

const double & InDet::TRT_TrajectoryElement_xk::dnegative ( ) const

inline

Definition at line 58 of file TRT_TrajectoryElement_xk.h.

{return m_dnegative ;}

dpositive()

const double & InDet::TRT_TrajectoryElement_xk::dpositive ( ) const

inline

Definition at line 57 of file TRT_TrajectoryElement_xk.h.

{return m_dpositive ;}

findCloseLink()

double InDet::TRT_TrajectoryElement_xk::findCloseLink	(	double	a,
		double	b )

Definition at line 539 of file TRT_TrajectoryElement_xk.cxx.

{
  m_bestlink =    -1;
  double dm = 10000.;
 
  for(int l=0; l < m_nlinks; ++l) {
 
    double v  = (a*m_link[l].way()+b)*m_link[l].way();
    double d  = std::abs(m_link[l].distance()-v);
 
    if(d < dm) {dm=d;  m_bestlink=l;}
  }
  return dm;
}

initiateForPrecisionSeed()

bool InDet::TRT_TrajectoryElement_xk::initiateForPrecisionSeed	(	bool	st,
		const InDetDD::TRT_BaseElement *	de,
		InDet::TRT_DriftCircleCollection::const_iterator &	sb,
		InDet::TRT_DriftCircleCollection::const_iterator &	se,
		std::pair< Amg::Vector3D, double > &	gp,
		const double *	dir,
		double	width2 )

Definition at line 41 of file TRT_TrajectoryElement_xk.cxx.

{
  m_isCluster  = st;
  m_detelement = de;
 
  // Test boundary active region of detector elements
  //
  if(!boundaryTest(3.,gp)) return false;
 
  // Trajectory element links production
  //
  initiateLinksForPrecisionSeed(sb,se,gp,dir,width2);
  return true;
}

initiateForTRTSeed()

bool InDet::TRT_TrajectoryElement_xk::initiateForTRTSeed	(	bool	st,
		const InDetDD::TRT_BaseElement *	de,
		InDet::TRT_DriftCircleCollection::const_iterator &	sb,
		InDet::TRT_DriftCircleCollection::const_iterator &	se,
		std::pair< Amg::Vector3D, double > &	gp,
		const double *	dir,
		double	width2 )

Definition at line 65 of file TRT_TrajectoryElement_xk.cxx.

{
 
  m_isCluster  = st  ;
  m_detelement = de  ;
 
  // Test boundary active region of detector elements
  //
  bool boundary = boundaryTest(7.,gp);
 
  // Trajectory element links production
  //
  initiateLinksForTRTSeed(sb,se,gp,dir,width2);
  return boundary;
}

initiateLinksForPrecisionSeed()

void InDet::TRT_TrajectoryElement_xk::initiateLinksForPrecisionSeed	(	InDet::TRT_DriftCircleCollection::const_iterator &	sb,
		InDet::TRT_DriftCircleCollection::const_iterator &	se,
		std::pair< Amg::Vector3D, double > &	gp,
		const double *	dir,
		double	width2 )

Definition at line 147 of file TRT_TrajectoryElement_xk.cxx.

{
  m_status = -1;
  m_bestlink = -1;
  m_nlinks = 0;
  double x = gp.first.x();
  double y = gp.first.y();
  double z = gp.first.z();
 
  const Amg::Transform3D& T = m_detelement->surface().transform();
 
  double step = std::abs(dir[0] * T(0, 2) + dir[1] * T(1, 2) + dir[2] * T(2, 2));
  step > .05 ? step = 1. / step : step = 20.;
 
  int Nstraws = m_detelement->nStraws();
 
  if (m_detelement->type() == InDetDD::TRT_BaseElement::BARREL) {
 
    for (int ns = 0; ns != Nstraws; ++ns) {
 
      const Amg::Transform3D& t = m_detelement->strawTransform(ns);
 
      double xs = t(0, 3);
      double ys = t(1, 3);
      double dx = x - xs;
      double dy = y - ys;
      double d = dx * dx + dy * dy;
      if (d > width2) {
        continue;
      }
 
      double Az[3] = { t(0, 2), t(1, 2), t(2, 2) };
      double D = dir[0] * Az[0] + dir[1] * Az[1] + dir[2] * Az[2];
      double A = (1. - D) * (1. + D);
      double dz = z - t(2, 3);
      double S = (dx * (D * Az[0] - dir[0]) + dy * (D * Az[1] - dir[1]) +
                  dz * (D * Az[2] - dir[2])) /
                 A;
      dx += (dir[0] * S);
      dy += (dir[1] * S);
      dz += (dir[2] * S);
      double Bx = Az[1] * dir[2] - Az[2] * dir[1];
      double By = Az[2] * dir[0] - Az[0] * dir[2];
      double Bz = Az[0] * dir[1] - Az[1] * dir[0];
      double im =
        (dx * Bx + dy * By + dz * Bz) / std::sqrt(Bx * Bx + By * By + Bz * Bz);
      double zl = dx * Az[0] + dy * Az[1] + dz * Az[2];
      S += gp.second;
      d = std::abs(im);
      if (y * xs - x * ys > 0.) {
        d = -d;
      }
      m_link[m_nlinks].set(ns, d, im, zl, S);
      if (++m_nlinks == 24) {
        break;
      }
    }
  } else {
 
    for (int ns = 0; ns != Nstraws; ++ns) {
 
      const Amg::Transform3D& t = m_detelement->strawTransform(ns);
 
      double xs = t(0, 3);
      double ys = t(1, 3);
      double as = m_radius / std::sqrt(xs * xs + ys * ys);
      double dx = x - xs * as;
      double dy = y - ys * as;
      double d = dx * dx + dy * dy;
      if (d > width2) {
        continue;
      }
 
      double zs = t(2, 3);
      double Az[3] = { t(0, 2), t(1, 2), t(2, 2) };
      double D = dir[0] * Az[0] + dir[1] * Az[1] + dir[2] * Az[2];
      double A = (1. - D) * (1. + D);
      dx = x - xs;
      dy = y - ys;
      double dz = z - zs;
      double S = (dx * (D * Az[0] - dir[0]) + dy * (D * Az[1] - dir[1]) +
                  dz * (D * Az[2] - dir[2])) /
                 A;
      dx += (dir[0] * S);
      dy += (dir[1] * S);
      dz += (dir[2] * S);
      double Bx = Az[1] * dir[2] - Az[2] * dir[1];
      double By = Az[2] * dir[0] - Az[0] * dir[2];
      double Bz = Az[0] * dir[1] - Az[1] * dir[0];
      double im =
        (dx * Bx + dy * By + dz * Bz) / std::sqrt(Bx * Bx + By * By + Bz * Bz);
      double zl = dx * Az[0] + dy * Az[1] + dz * Az[2];
      S += gp.second;
      d = std::abs(im);
      if (y * xs - x * ys > 0.) {
        d = -d;
      }
      m_link[m_nlinks].set(ns, d, im, zl, S);
      if (++m_nlinks == 24) {
        break;
      }
    }
  }
  m_dpositive = 1000.;
  m_dnegative = -1000.;
  if (m_isCluster && m_nlinks) {
 
    bool nl = false;
    for (; sb != se; ++sb) {
 
      int ns = m_trtid->straw((*sb)->identify());
 
      for (int l = 0; l != m_nlinks; ++l) {
        if (ns != m_link[l].number()) {
          continue;
        }
        nl = true;
        m_link[l].set((*sb));
 
        double d = m_link[l].distance();
        if (d >= 0.) {
          if (d < m_dpositive)
            m_dpositive = d;
        } else {
          if (d > m_dnegative) {
            m_dnegative = d;
          }
        }
        break;
      }
    }
    if (!nl) {
      m_isCluster = false;
    }
  }
}

initiateLinksForTRTSeed()

void InDet::TRT_TrajectoryElement_xk::initiateLinksForTRTSeed	(	InDet::TRT_DriftCircleCollection::const_iterator &	sb,
		InDet::TRT_DriftCircleCollection::const_iterator &	se,
		std::pair< Amg::Vector3D, double > &	gp,
		const double *	dir,
		double	width2 )

Definition at line 293 of file TRT_TrajectoryElement_xk.cxx.

{
  m_status = -1;
  m_bestlink = -1;
  m_nlinks = 0;
  double x = gp.first.x();
  double y = gp.first.y();
  double z = gp.first.z();
 
  const Amg::Transform3D& T = m_detelement->surface().transform();
 
  double step =
    std::abs(dir[0] * T(0, 2) + dir[1] * T(1, 2) + dir[2] * T(2, 2));
  step > .05 ? step = 1. / step : step = 20.;
 
  int Nstraws = m_detelement->nStraws();
  if (m_detelement->type() == InDetDD::TRT_BaseElement::BARREL) {
 
    for (int ns = 0; ns != Nstraws; ++ns) {
 
      const Amg::Transform3D& t = m_detelement->strawTransform(ns);
 
      double xs = t(0, 3);
      double ys = t(1, 3);
      double dx = x - xs;
      double dy = y - ys;
      double d = dx * dx + dy * dy;
      if (d > width2){
        continue;
      }
 
      double Az[3] = { t(0, 2), t(1, 2), t(2, 2) };
      double D = dir[0] * Az[0] + dir[1] * Az[1] + dir[2] * Az[2];
      double A = (1. - D) * (1. + D);
      double dz = z - t(2, 3);
      double S = (dx * (D * Az[0] - dir[0]) + dy * (D * Az[1] - dir[1]) +
                  dz * (D * Az[2] - dir[2])) /
                 A;
      dx += (dir[0] * S);
      dy += (dir[1] * S);
      dz += (dir[2] * S);
      double Bx = Az[1] * dir[2] - Az[2] * dir[1];
      double By = Az[2] * dir[0] - Az[0] * dir[2];
      double Bz = Az[0] * dir[1] - Az[1] * dir[0];
      double im =
        (dx * Bx + dy * By + dz * Bz) / std::sqrt(Bx * Bx + By * By + Bz * Bz);
      double zl = dx * Az[0] + dy * Az[1] + dz * Az[2];
      S += gp.second;
      d = std::abs(im);
      if (y * xs - x * ys > 0.){
        d = -d;
      }
      m_link[m_nlinks].set(ns, d, im, zl, S);
      if (++m_nlinks == 24){
        break;
      }
    }
  } else {
 
    double ri = 1. / m_radius;
    double ax = x * ri;
    double ay = y * ri;
 
    for (int ns = 0; ns != Nstraws; ++ns) {
 
      const Amg::Transform3D& t = m_detelement->strawTransform(ns);
 
      double xs = t(0, 3);
      double ys = t(1, 3);
      double as = m_radius / std::sqrt(xs * xs + ys * ys);
      double dx = x - xs * as;
      double dy = y - ys * as;
      double d = dx * dx + dy * dy;
      if (d > width2){
        continue;
      }
 
      double zs = t(2, 3);
      double Az[3] = { t(0, 2), t(1, 2), t(2, 2) };
      double D = dir[0] * Az[0] + dir[1] * Az[1] + dir[2] * Az[2];
      double A = 1. / ((1. - D) * (1. + D));
      dx = x - xs;
      dy = y - ys;
      double dz = z - zs;
      double Dx = (D * Az[0] - dir[0]);
      double Dy = (D * Az[1] - dir[1]);
      double S = (dx * Dx + dy * Dy + dz * (D * Az[2] - dir[2])) * A;
      dx += (dir[0] * S);
      dy += (dir[1] * S);
      dz += (dir[2] * S);
      double Bx = Az[1] * dir[2] - Az[2] * dir[1];
      double By = Az[2] * dir[0] - Az[0] * dir[2];
      double Bz = Az[0] * dir[1] - Az[1] * dir[0];
      double B = 1. / std::sqrt(Bx * Bx + By * By + Bz * Bz);
      double im = (dx * Bx + dy * By + dz * Bz) * B;
      double zl = dx * Az[0] + dy * Az[1] + dz * Az[2];
 
      // d(im)/dr calculation for endcap
      //
      double dS = (ax * Dx + ay * Dy) * A;
      double sx = ax + dir[0] * dS;
      double sy = ay + dir[1] * dS;
      double sd = (sx * Bx + sy * By) * B;
      double sz = sx * Az[0] + sy * Az[1];
      S += gp.second;
      if (y * xs - x * ys > 0.) {
        d = -std::abs(im);
        sd = -std::abs(im + sd) - d;
      } else {
        d = std::abs(im);
        sd = std::abs(im + sd) - d;
      }
 
      m_link[m_nlinks].set(ns, d, im, zl, S, sd, sz);
      if (++m_nlinks == 24){
        break;
      }
    }
  }
 
  m_dpositive = 1000.;
  m_dnegative = -1000.;
  if (m_isCluster && m_nlinks) {
 
    bool nl = false;
    for (; sb != se; ++sb) {
 
      int ns = m_trtid->straw((*sb)->identify());
 
      for (int l = 0; l != m_nlinks; ++l) {
        if (ns != m_link[l].number())
          continue;
        nl = true;
        m_link[l].set((*sb));
 
        double d = m_link[l].distance();
        if (d >= 0.) {
          if (d < m_dpositive){
            m_dpositive = d;
          }
        } else {
          if (d > m_dnegative){
            m_dnegative = d;
          }
        }
        break;
      }
    }
    if (!nl){
      m_isCluster = false;
    }
  }
}

isBarrel()

const bool & InDet::TRT_TrajectoryElement_xk::isBarrel ( ) const

inline

Definition at line 49 of file TRT_TrajectoryElement_xk.h.

{return m_barrel    ;}

isCluster()

const bool & InDet::TRT_TrajectoryElement_xk::isCluster ( ) const

inline

Definition at line 48 of file TRT_TrajectoryElement_xk.h.

{return m_isCluster ;}

link()

const TRT_ExtensionDriftCircleLink_xk & InDet::TRT_TrajectoryElement_xk::link ( int i ) const

inline

Definition at line 60 of file TRT_TrajectoryElement_xk.h.

{return m_link[i]    ;}

nlinks()

const int & InDet::TRT_TrajectoryElement_xk::nlinks ( ) const

inline

Definition at line 51 of file TRT_TrajectoryElement_xk.h.

{return m_nlinks    ;}

operator=()

TRT_TrajectoryElement_xk & InDet::TRT_TrajectoryElement_xk::operator= ( const TRT_TrajectoryElement_xk & E )

inline

Definition at line 200 of file TRT_TrajectoryElement_xk.h.

    {
      m_status      = E.m_status     ;
      m_isCluster   = E.m_isCluster  ;
      m_barrel      = E.m_barrel     ;
      m_bestlink    = E.m_bestlink   ;
      m_nlinks      = E.m_nlinks     ;
      m_detelement  = E.m_detelement ;
      m_trtid       = E.m_trtid      ;
      m_proptool    = E.m_proptool   ;
      m_updatortool = E.m_updatortool;
      m_fieldprop   = E.m_fieldprop  ;
      m_riomakerD   = E.m_riomakerD  ;
      m_riomakerN   = E.m_riomakerN  ;
      m_z           = E.m_z          ;
      m_zMin        = E.m_zMin       ;
      m_zMax        = E.m_zMax       ;
      m_radius      = E.m_radius     ;
      m_radiusMin   = E.m_radiusMin  ;
      m_radiusMax   = E.m_radiusMax  ;
      m_scale_error = E.m_scale_error;
      m_dpositive   = E.m_dpositive  ;
      m_dnegative   = E.m_dnegative  ;
      for(int i=0; i!=m_nlinks; ++i) {m_link[i]=E.m_link[i];}
      return(*this);
    }

polarAngleEstimation()

void InDet::TRT_TrajectoryElement_xk::polarAngleEstimation	(	TRT_TrajectoryElement_xk *	E,
		Amg::Vector3D &	G0,
		Amg::Vector3D &	G1,
		double	C,
		double	VZ,
		double *	Tp ) const

Definition at line 881 of file TRT_TrajectoryElement_xk.cxx.

{
  double dx = G1[0]-G0[0]           ;
  double dy = G1[1]-G0[1]           ;
  double dz = G1[2]-G0[2]           ;
  double dr = std::sqrt(dx*dx+dy*dy)     ;
  double rc = dr*C                  ;
  Tp[0]     = dz/(dr*(1.+.04*rc*rc));
 
  double r0 = G0[0]*G0[0]+G0[1]*G0[1];
  double r1 = G1[0]*G1[0]+G1[1]*G1[1];
 
  if( r0 > r1) {
 
    Tp[0]  =-Tp[0];
    if     (  m_barrel )   { //---------------------------->BB
 
      double dZ = m_zMax-m_zMin              ;
      double dT = (dZ+2.*VZ)                 ;
      Tp[1]     = dZ*dZ*.1                   ;
      Tp[2]     = dT*dT*.1/r0                ;
    }
    else if(E->m_barrel)   { //---------------------------->EB
 
      double T2 = Tp[0]*Tp[0]                ;
      Tp[1]     = 400.                       ;
      Tp[2]     = (VZ*VZ)/((1.+T2)*r0)       ;
    }
    else                   { //---------------------------->EE
 
      double T2 = Tp[0]*Tp[0]                ;
      Tp[1]     = 400.                       ;
      Tp[2]     = (VZ*VZ)/((1.+T2)*r0)       ;
    }
  }
  else {
    if     (E->m_barrel)   { //---------------------------->BB
 
      double dZ = m_zMax-m_zMin              ;
      double dT = (E->m_zMax-E->m_zMin+2.*VZ);
      Tp[1]     = dZ*dZ*.1                   ;
      Tp[2]     = dT*dT*.1/r1                ;
    }
    else if(  m_barrel )   { //---------------------------->BE
 
      double T2 = Tp[0]*Tp[0]                ;
      double dr = std::sqrt(r1)-std::sqrt(r0)          ;
      Tp[1]     = dr*dr*VZ*VZ/r1             ;
      Tp[2]     = (VZ*VZ)/((1.+T2)*r1)       ;
 
    }
    else                   { //---------------------------->EE
 
      double T2 = Tp[0]*Tp[0]                ;
      Tp[1]     = 400.                       ;
      Tp[2]     = (VZ*VZ)/((1.+T2)*r1)       ;
    }
  }
}

propagate()

bool InDet::TRT_TrajectoryElement_xk::propagate	(	Trk::PatternTrackParameters &	Ta,
		Trk::PatternTrackParameters &	Tb )

Definition at line 600 of file TRT_TrajectoryElement_xk.cxx.

{
  if(m_bestlink<0) return false;
 
  return m_proptool->propagate
    (Gaudi::Hive::currentContext(),
     Ta,m_link[m_bestlink].surface(),Tb,Trk::anyDirection,m_fieldprop,Trk::pion);
}

radius()

const double & InDet::TRT_TrajectoryElement_xk::radius ( ) const

inline

Definition at line 53 of file TRT_TrajectoryElement_xk.h.

{return m_radius    ;}

radiusCorrection()

void InDet::TRT_TrajectoryElement_xk::radiusCorrection ( const double & dr )

inline

Definition at line 230 of file TRT_TrajectoryElement_xk.h.

    {
      if(m_bestlink >= 0) m_link[m_bestlink].radiusCorrection(dr);
    }

radiusMax()

const double & InDet::TRT_TrajectoryElement_xk::radiusMax ( ) const

inline

Definition at line 56 of file TRT_TrajectoryElement_xk.h.

{return m_radiusMax ;}

radiusMin()

const double & InDet::TRT_TrajectoryElement_xk::radiusMin ( ) const

inline

Definition at line 55 of file TRT_TrajectoryElement_xk.h.

{return m_radiusMin ;}

rioOnTrack()

const Trk::RIO_OnTrack * InDet::TRT_TrajectoryElement_xk::rioOnTrack

(

)

Definition at line 559 of file TRT_TrajectoryElement_xk.cxx.

{
  if(m_bestlink < 0 || m_status<=0) return nullptr;
 
  int l = m_bestlink;
 
  const Trk::StraightLineSurface* line = static_cast<const Trk::StraightLineSurface*>
    (&(m_link[l].cluster()->detectorElement())->surface(m_link[l].cluster()->identify()));
  Trk::AtaStraightLine Tp(m_link[l].impact(),m_link[l].zlocal(),1.,1.,1.,*line);
  if(m_status==2)
    return m_riomakerD->correct(*m_link[l].cluster(),Tp, Gaudi::Hive::currentContext());
  else
    return m_riomakerN->correct(*m_link[l].cluster(),Tp, Gaudi::Hive::currentContext());
}

rioOnTrackSimple()

std::unique_ptr< Trk::RIO_OnTrack > InDet::TRT_TrajectoryElement_xk::rioOnTrackSimple

(

)

const

Definition at line 579 of file TRT_TrajectoryElement_xk.cxx.

{
  if(m_bestlink < 0 || m_status<=0) return nullptr;
 
  int l = m_bestlink;
  Amg::Vector3D dir(1.,0.,0.);
  Trk::DefinedParameter radius(0.,Trk::locX);
  Amg::MatrixX cov(1,1);
  cov(0,0) = 1.;
  if(m_status==2) {
    return std::make_unique<InDet::TRT_DriftCircleOnTrack>(m_link[l].cluster(),Trk::LocalParameters(radius),
                                                           std::move(cov),0,0.,dir,Trk::DECIDED);
  }
  return  std::make_unique<InDet::TRT_DriftCircleOnTrack>(m_link[l].cluster(),Trk::LocalParameters(radius),
                                                          std::move(cov),0,0.,dir,Trk::NODRIFTTIME);
}

set() [1/2]

void InDet::TRT_TrajectoryElement_xk::set	(	const Trk::MagneticFieldProperties &	mf,
		const AtlasFieldCacheCondObj *	fieldCondObj )

Definition at line 30 of file TRT_TrajectoryElement_xk.cxx.

{
  m_fieldprop    = mf;
  fieldCondObj->getInitializedCache (m_fieldCache);
}

set() [2/2]

void InDet::TRT_TrajectoryElement_xk::set	(	const TRT_ID *	m,
		const Trk::IPatternParametersPropagator *	pr,
		const Trk::IPatternParametersUpdator *	up,
		const Trk::IRIO_OnTrackCreator *	riod,
		const Trk::IRIO_OnTrackCreator *	rion,
		double	scale )

Definition at line 14 of file TRT_TrajectoryElement_xk.cxx.

{
  m_trtid       = m               ;
  m_riomakerD   = riod            ;
  m_riomakerN   = rion            ;
  m_proptool    = pr              ;
  m_updatortool = up              ;
  m_scale_error = scale           ;
}

status()

const int & InDet::TRT_TrajectoryElement_xk::status ( ) const

inline

Definition at line 50 of file TRT_TrajectoryElement_xk.h.

{return m_status    ;}

trackParametersEstimation() [1/2]

bool InDet::TRT_TrajectoryElement_xk::trackParametersEstimation	(	TRT_TrajectoryElement_xk *	E1,
		Trk::PatternTrackParameters &	Tp,
		double	ZvHW )

Definition at line 790 of file TRT_TrajectoryElement_xk.cxx.

{
 
  if(!E1) return false;
 
  Amg::Vector3D Gp[3];
  double     Wa[3];
  if   (!    trajectoryGlobalPosition(Gp[1],Wa[1])) return false;
  if   (!E1->trajectoryGlobalPosition(Gp[2],Wa[2])) return false;
 
  double dx = Gp[2][0]-Gp[1][0];
  double dy = Gp[2][1]-Gp[1][1];
  Gp[0][0]  = 0.;
  Gp[0][1]  = 0.;
  Gp[0][2]  = Gp[1][2]-(Gp[2][2]-Gp[1][2])*std::sqrt((Gp[1][0]*Gp[1][0]+Gp[1][1]*Gp[1][1])/(dx*dx+dy*dy));
  Wa[0]     = 0.;
 
  double x0 = Gp[0][0]               ;
  double y0 = Gp[0][1]               ;
  double x1 = Gp[1][0]-x0            ;
  double y1 = Gp[1][1]-y0            ;
  double x2 = Gp[2][0]-x0            ;
  double y2 = Gp[2][1]-y0            ;
  double r1 = std::sqrt(x1*x1+y1*y1)      ;
  double u1 = 1./r1                  ;
  double r2 = 1./(x2*x2+y2*y2)       ;
  double a  = x1*u1                  ;
  double b  = y1*u1                  ;
  double u2 = (a*x2+b*y2)*r2         ;
  double v2 = (a*y2-b*x2)*r2         ;
  double A  = v2/(u2-u1)             ;
  double B  = 2.*(v2-A*u2)           ;
  double C  = B/std::sqrt(1.+A*A)         ;
 
  double f,TP[3];
 
  f = atan2( b+a*A, a-b*A); polarAngleEstimation(E1,Gp[1],Gp[2],C,ZvHW,TP);
 
  // Track parameters estimation
  //
  double P[5] = {0.,Gp[0][2],f,atan2(1.,TP[0]),.0001};
 
  double pos0[3]; pos0[0]=Gp[0][0]; pos0[1]=Gp[0][1]; pos0[2]=Gp[0][2];
  double pos1[3]; pos1[0]=Gp[1][0]; pos1[1]=Gp[1][1]; pos1[2]=Gp[1][2];
  double pos2[3]; pos2[0]=Gp[2][0]; pos2[1]=Gp[2][1]; pos2[2]=Gp[2][2];
  double H0  [3];
  double H1  [3];
  double H2  [3];
  m_fieldCache.getFieldZR    (pos0,H0);
  m_fieldCache.getFieldZR    (pos1,H1);
  m_fieldCache.getFieldZR    (pos2,H2);
 
  double Hz = .333333*(H0[2]+H1[2]+H2[2]);
 
  // If magnetic field exist
  //
  double Cm =          1.;
  double T2 = TP[0]*TP[0];
  if(std::abs(Hz)>1.e-9) {Cm = 1./(300.*Hz*std::sqrt(1.+T2)); P[4] = -C*Cm;}
 
  // Covariance of track parameters estimation
  //
  double wa = 1./Wa[2];
 
  const double dS = 4.           ;
  double df       = dS*wa        ;
  double dp       = 8.*df*wa*Cm  ;
  double c0       = .1           ;
  double c1       = ZvHW*ZvHW*.1 ;
  double c2       = df*df*(1.+T2);
  double c3       = TP[2]        ;
  double c4       = dp*dp        ;
  double V[15] ={c0,
         0.,c1,
         0.,0.,c2,
         0.,0.,0.,c3,
         0.,0.,0.,0.,c4};
 
  Tp.setParametersWithCovariance(nullptr,P,V);
  return true;
}

trackParametersEstimation() [2/2]

bool InDet::TRT_TrajectoryElement_xk::trackParametersEstimation	(	TRT_TrajectoryElement_xk *	E1,
		TRT_TrajectoryElement_xk *	E2,
		Trk::PatternTrackParameters &	Tp,
		double	ZvHW )

Definition at line 685 of file TRT_TrajectoryElement_xk.cxx.

{
 
  if(!E1) return false;
 
  Amg::Vector3D Gp[3];
  double     Wa[3];
 
  if   (!    trajectoryGlobalPosition(Gp[0],Wa[0])) return false;
  if   (!E1->trajectoryGlobalPosition(Gp[1],Wa[1])) return false;
 
  int mode; Wa[0] > Wa[1] ? mode=0 : mode=1;
 
  if(E2) {
    if(!E2->trajectoryGlobalPosition(Gp[2],Wa[2])) return false;
  }
  else    {
    double dx = Gp[1][0]-Gp[0][0];
    double dy = Gp[1][1]-Gp[0][1];
    Gp[2][0]  = 0.;
    Gp[2][1]  = 0.;
    Gp[2][2]  = Gp[0][2]-(Gp[1][2]-Gp[0][2])*std::sqrt((Gp[0][0]*Gp[0][0]+Gp[0][1]*Gp[0][1])/(dx*dx+dy*dy));
    Wa[2]     = 0.;
    mode      = 2 ;
  }
 
  double x0 = Gp[0][0]               ;
  double y0 = Gp[0][1]               ;
  double x1 = Gp[1][0]-x0            ;
  double y1 = Gp[1][1]-y0            ;
  double x2 = Gp[2][0]-x0            ;
  double y2 = Gp[2][1]-y0            ;
  double r1 = std::sqrt(x1*x1+y1*y1)      ;
  double u1 = 1./r1                  ;
  double r2 = 1./(x2*x2+y2*y2)       ;
  double a  = x1*u1                  ;
  double b  = y1*u1                  ;
  double u2 = (a*x2+b*y2)*r2         ;
  double v2 = (a*y2-b*x2)*r2         ;
  double A  = v2/(u2-u1)             ;
  double B  = 2.*(v2-A*u2)           ;
  double C  = B/std::sqrt(1.+A*A)         ;
 
  double f,TP[3];
 
  if(mode==0) {
    f = atan2(-b-a*A,-a+b*A); C=-C; polarAngleEstimation(E2,Gp[0],Gp[2],C,ZvHW,TP);
  }
  else if(mode==1) {
    f = atan2( b+a*A, a-b*A);       polarAngleEstimation(E2,Gp[0],Gp[2],C,ZvHW,TP);
  }
  else             {
    f = atan2( b+a*A, a-b*A);       polarAngleEstimation(E1,Gp[0],Gp[1],C,ZvHW,TP);
  }
 
  // Track parameters estimation
  //
  double P[5] = {m_link[m_bestlink].impact(),m_link[m_bestlink].zlocal(),f,atan2(1.,TP[0]),.0001};
 
  double pos0[3]; pos0[0]=Gp[0][0]; pos0[1]=Gp[0][1]; pos0[2]=Gp[0][2];
  double pos1[3]; pos1[0]=Gp[1][0]; pos1[1]=Gp[1][1]; pos1[2]=Gp[1][2];
  double pos2[3]; pos2[0]=Gp[2][0]; pos2[1]=Gp[2][1]; pos2[2]=Gp[2][2];
  double H0  [3];
  double H1  [3];
  double H2  [3];
  m_fieldCache.getFieldZR    (pos0,H0);
  m_fieldCache.getFieldZR    (pos1,H1);
  m_fieldCache.getFieldZR    (pos2,H2);
 
  double Hz = .333333*(H0[2]+H1[2]+H2[2]);
 
  // If magnetic field exist
  //
  double Cm = 1.         ;
  double T2 = TP[0]*TP[0];
  if(std::abs(Hz)>1.e-9) {Cm = 1./(300.*Hz*std::sqrt(1.+T2)); P[4] = -C*Cm;}
 
  // Covariance of track parameters estimation
  //
  double wa; mode != 2 ? wa = 1./(Wa[0]-Wa[2]) : wa = 1./Wa[1];
  const double dS = 4.         ;
 
  double df       = dS*wa        ;
  double dp       = 8.*df*wa*Cm  ;
  double c0       = 4.           ;
  double c1       = TP[1]        ;
  double c2       = df*df*(1.+T2);
  double c3       = TP[2]        ;
  double c4       = dp*dp        ;
  double V[15] ={c0,
         0.,c1,
         0.,0.,c2,
         0.,0.,0.,c3,
         0.,0.,0.,0.,c4};
 
  Tp.setParametersWithCovariance(&m_link[m_bestlink].surface(),P,V);
  return true;
}

trajectoryGlobalPosition()

bool InDet::TRT_TrajectoryElement_xk::trajectoryGlobalPosition	(	Amg::Vector3D &	G,
		double &	WAY )

Definition at line 646 of file TRT_TrajectoryElement_xk.cxx.

{
  int l = m_bestlink;
  if( l<0 ) return false;
 
  const Amg::Transform3D& t = m_detelement->strawTransform(m_link[l].number());
 
  double Az[3] = {t(0,2),t(1,2),t(2,2)};
  double Rc[3] = {t(0,3),t(1,3),t(2,3)};
 
  double Bx,By;
 
  if(std::abs(Az[2]) > .7) {              // Barrel
 
    double Ri = 1./std::sqrt(Rc[0]*Rc[0]+Rc[1]*Rc[1]);
    Bx =-Az[2]*Rc[1]*Ri; By = Az[2]*Rc[0]*Ri;
  }
  else if(Rc[2] > 0. ) {             // Positive endcap
 
    Bx = Az[1];          By =-Az[0];
  }
  else                 {             // Negative endcap
 
    Bx =-Az[1];          By = Az[0];
  }
  double zl = m_link[l].zlocal();
  double im = m_link[l].impact();
  G[0]      = zl*Az[0]+Bx*im+Rc[0];
  G[1]      = zl*Az[1]+By*im+Rc[1];
  G[2]      = zl*Az[2]      +Rc[2];
  WAY       = m_link[l].way()     ;
  return true;
}

z()

const double & InDet::TRT_TrajectoryElement_xk::z ( ) const

inline

Definition at line 54 of file TRT_TrajectoryElement_xk.h.

{return m_z         ;}

Member Data Documentation

m_barrel

bool InDet::TRT_TrajectoryElement_xk::m_barrel

protected

Definition at line 128 of file TRT_TrajectoryElement_xk.h.

m_bestlink

int InDet::TRT_TrajectoryElement_xk::m_bestlink

protected

Definition at line 131 of file TRT_TrajectoryElement_xk.h.

m_detelement

const InDetDD::TRT_BaseElement* InDet::TRT_TrajectoryElement_xk::m_detelement

protected

Definition at line 143 of file TRT_TrajectoryElement_xk.h.

m_dnegative

double InDet::TRT_TrajectoryElement_xk::m_dnegative

protected

Definition at line 141 of file TRT_TrajectoryElement_xk.h.

m_dpositive

double InDet::TRT_TrajectoryElement_xk::m_dpositive

protected

Definition at line 140 of file TRT_TrajectoryElement_xk.h.

m_fieldCache

MagField::AtlasFieldCache InDet::TRT_TrajectoryElement_xk::m_fieldCache

protected

Definition at line 150 of file TRT_TrajectoryElement_xk.h.

m_fieldprop

Trk::MagneticFieldProperties InDet::TRT_TrajectoryElement_xk::m_fieldprop

protected

Definition at line 149 of file TRT_TrajectoryElement_xk.h.

m_isCluster

bool InDet::TRT_TrajectoryElement_xk::m_isCluster

protected

Definition at line 129 of file TRT_TrajectoryElement_xk.h.

m_link

TRT_ExtensionDriftCircleLink_xk InDet::TRT_TrajectoryElement_xk::m_link[24]

protected

Definition at line 142 of file TRT_TrajectoryElement_xk.h.

m_nlinks

int InDet::TRT_TrajectoryElement_xk::m_nlinks

protected

Definition at line 132 of file TRT_TrajectoryElement_xk.h.

m_proptool

const Trk::IPatternParametersPropagator* InDet::TRT_TrajectoryElement_xk::m_proptool

protected

Definition at line 145 of file TRT_TrajectoryElement_xk.h.

m_radius

double InDet::TRT_TrajectoryElement_xk::m_radius

protected

Definition at line 136 of file TRT_TrajectoryElement_xk.h.

m_radiusMax

double InDet::TRT_TrajectoryElement_xk::m_radiusMax

protected

Definition at line 138 of file TRT_TrajectoryElement_xk.h.

m_radiusMin

double InDet::TRT_TrajectoryElement_xk::m_radiusMin

protected

Definition at line 137 of file TRT_TrajectoryElement_xk.h.

m_riomakerD

const Trk::IRIO_OnTrackCreator* InDet::TRT_TrajectoryElement_xk::m_riomakerD

protected

Definition at line 147 of file TRT_TrajectoryElement_xk.h.

m_riomakerN

const Trk::IRIO_OnTrackCreator* InDet::TRT_TrajectoryElement_xk::m_riomakerN

protected

Definition at line 148 of file TRT_TrajectoryElement_xk.h.

m_scale_error

double InDet::TRT_TrajectoryElement_xk::m_scale_error

protected

Definition at line 139 of file TRT_TrajectoryElement_xk.h.

m_status

int InDet::TRT_TrajectoryElement_xk::m_status

protected

Definition at line 130 of file TRT_TrajectoryElement_xk.h.

m_trtid

const TRT_ID* InDet::TRT_TrajectoryElement_xk::m_trtid

protected

Definition at line 144 of file TRT_TrajectoryElement_xk.h.

m_updatortool

const Trk::IPatternParametersUpdator* InDet::TRT_TrajectoryElement_xk::m_updatortool

protected

Definition at line 146 of file TRT_TrajectoryElement_xk.h.

m_z

double InDet::TRT_TrajectoryElement_xk::m_z

protected

Definition at line 133 of file TRT_TrajectoryElement_xk.h.

m_zMax

double InDet::TRT_TrajectoryElement_xk::m_zMax

protected

Definition at line 135 of file TRT_TrajectoryElement_xk.h.

m_zMin

double InDet::TRT_TrajectoryElement_xk::m_zMin

protected

Definition at line 134 of file TRT_TrajectoryElement_xk.h.

---

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

• TRT_TrajectoryElement_xk.h
• TRT_TrajectoryElement_xk.cxx