TrackPropagationHelper Class Reference

#include <TrackPropagationHelper.h>

Inheritance diagram for TrackPropagationHelper:

Collaboration diagram for TrackPropagationHelper:

Classes
class	Imp

Public Member Functions
	TrackPropagationHelper (IVP1System *sys)
	~TrackPropagationHelper ()
bool	makePointsNeutral (std::vector< Amg::Vector3D > &points, const Trk::Track *)
bool	makePointsCharged (std::vector< Amg::Vector3D > &points, const Trk::Track *, Trk::IExtrapolator *extrapolator, Trk::ParticleHypothesis hypo=Trk::nonInteracting, bool useMEOT=false, const Trk::Volume *volume=0)
std::vector< Trk::PlaneSurface > &	getExtrapolationSurfaces () const
bool	showExtrapolationSurfaces () const
void	setSystemBasePointer (IVP1System *sys)
void	message (const QString &) const
void	message (const QStringList &, const QString &addtoend="") const
void	message (const QString &addtostart, const QStringList &, const QString &addtoend="") const
void	messageDebug (const QString &) const
void	messageDebug (const QStringList &, const QString &addtoend="") const
void	messageDebug (const QString &addtostart, const QStringList &, const QString &addtoend="") const
void	messageVerbose (const QString &) const
void	messageVerbose (const QStringList &, const QString &addtoend="") const
void	messageVerbose (const QString &addtostart, const QStringList &, const QString &addtoend="") const
IVP1System *	systemBase () const
const QString &	helperClassName () const

Static Public Member Functions
static bool	verbose ()
static void	warnUndeletedInstances ()
static QString	str (const QString &s)
static QString	str (const char *c)
static QString	str (const bool b)
static QString	str (const QColor &)
static QString	str (const SbColor &)
static QString	str (const VP1Interval &)
static QString	str (const SbVec2d &)
static QString	str (const SbVec2f &)
static QString	str (const SbVec2s &)
static QString	str (const SbVec3d &)
static QString	str (const SbVec3f &)
static QString	str (const SbVec3s &)
static QString	str (const SbVec4d &)
static QString	str (const SbVec4f &)
static QString	str (short int n)
static QString	str (unsigned short int n)
static QString	str (long n)
static QString	str (ulong n)
static QString	str (int n)
static QString	str (uint n)
static QString	str (qlonglong n)
static QString	str (qulonglong n)
static QString	str (const double &d)
static QString	str (const float &f)
static QString	str (const void *)
template<class T>
static QString	str (const T *t)
template<class T>
static QString	str (const QFlags< T > &f)
template<class T>
static QString	str (const HepGeom::BasicVector3D< T > &t)
static QString	str (const Amg::Vector3D &t)
template<class T>
static QString	str (const QList< T > &t)

Protected Member Functions
void	setHelperClassName (const QString &n)

Private Attributes
Imp *	m_d
QString	m_helpername
IVP1System *	m_system

Static Private Attributes
static const bool	s_vp1verbose = VP1QtUtils::environmentVariableIsOn("VP1_VERBOSE_OUTPUT")

Detailed Description

Definition at line 40 of file TrackPropagationHelper.h.

Constructor & Destructor Documentation

TrackPropagationHelper()

TrackPropagationHelper::TrackPropagationHelper

(

IVP1System *

sys

)

Definition at line 140 of file TrackPropagationHelper.cxx.

  : VP1HelperClassBase(sys,"TrackPropagationHelper"), m_d(new Imp(this))
{
}

~TrackPropagationHelper()

TrackPropagationHelper::~TrackPropagationHelper

(

)

Definition at line 146 of file TrackPropagationHelper.cxx.

{
  delete m_d;
}

Member Function Documentation

getExtrapolationSurfaces()

std::vector< Trk::PlaneSurface > & TrackPropagationHelper::getExtrapolationSurfaces

(

)

const

Definition at line 97 of file TrackPropagationHelper.cxx.

                                                                                   {
  return m_d->surfaces;
}

helperClassName()

const QString & VP1HelperClassBase::helperClassName ( ) const

inlineinherited

Definition at line 51 of file VP1HelperClassBase.h.

{ return m_helpername; }

makePointsCharged()

bool TrackPropagationHelper::makePointsCharged	(	std::vector< Amg::Vector3D > &	points,
		const Trk::Track *	track,
		Trk::IExtrapolator *	extrapolator,
		Trk::ParticleHypothesis	hypo = Trk::nonInteracting,
		bool	useMEOT = false,
		const Trk::Volume *	volume = 0 )

Definition at line 240 of file TrackPropagationHelper.cxx.

{
  if (VP1Msg::verbose())
    messageVerbose("makePointsCharged start with hypo="+str(hypo)+", useMEOT="+str(useMEOT)+", volume=" +str(volume));
  // if (volume) std::cout<<volume->volumeBounds()<<std::endl;
    //  Clear input points, check/sanitise input and initialise  //
 
  points.clear();
  if (!extrapolator) {
    message("makePointsCharged ERROR: Null extrapolator tool provided!");
    return false;
  }
  if (!track) {
    message("makePointsCharged ERROR: Called with null track pointer");
    return false;
  }
  const unsigned npars = track->trackParameters()->size();
  if (npars==0) {
    message("makePointsCharged ERROR: No TrackParameters on track!");
    return false;
  }
  if (npars==1)
    return m_d->makePointsCharged_SinglePar(points,track,extrapolator,hypo);
 
  points.reserve(npars);//At least we need this.
  const EventContext& ctx = Gaudi::Hive::currentContext();
    //Add a point for each parameter, and add extra points between them where appropriate.
  Trk::TrackStates::const_iterator tsos_iter = track->trackStateOnSurfaces()->begin();
  Trk::TrackStates::const_iterator tsos_end = track->trackStateOnSurfaces()->end();
  const Trk::TrackParameters* prevpar(nullptr);
  const Trk::TrackParameters* trackParam(nullptr);
  bool problems(false);
  for (; tsos_iter != tsos_end; ++tsos_iter) {
    if (!m_d->tracksanity->isSafe(*tsos_iter))
      continue;
    if ((*tsos_iter)->measurementOnTrack()==nullptr && ( (*tsos_iter)->materialEffectsOnTrack()&&!useMEOT ) )
      continue;
    trackParam = (*tsos_iter)->trackParameters();
    if (!m_d->tracksanity->isSafe(trackParam))
      continue;
    if (!prevpar) {
        //first time.
      prevpar = trackParam;
      points.push_back(prevpar->position());
      continue;
    }
    if (!m_d->addPointsBetweenParameters_Charged(points,track,prevpar,trackParam,extrapolator,hypo))
      problems = true;
    points.push_back(trackParam->position());
    prevpar = trackParam;
  }
 
  if (problems)
    messageDebug("WARNING: Problems encountered adding point(s) between track parameters");
 
// restrict to ID tracks for now.
  if (volume && trackParam && !m_d->outsideIDVolume(trackParam->position())) {
    messageVerbose("Extending to Volume");
    //get individual surfaces
 
    // TODO - optimise this!
    const std::vector<std::unique_ptr<Trk::Surface>> bsurfs =
      const_cast<Trk::VolumeBounds&>(volume->volumeBounds())
        .decomposeToSurfaces(volume->transform());
 
    if (!bsurfs.empty()){
      messageVerbose("Has this many surfaces:"+str(bsurfs.size()));
 
      auto bSurfsIt = bsurfs.begin();
      for (;bSurfsIt!= bsurfs.end(); ++bSurfsIt){
 
        messageVerbose("Extrap value:"+str((extrapolator)));
        messageVerbose("trackParam:"+str((trackParam)));
        const Trk::TrackParameters*  trackPar = extrapolator->extrapolate(ctx,
                                                                          *trackParam,
                                                                          **bSurfsIt,
                                                                          Trk::alongMomentum,true,hypo).release(); // change this to extrapolate current param to surface.
 
        if (trackPar){
          messageVerbose("Extrapolation succeeded");
 
          if (!m_d->addPointsBetweenParameters_Charged(points,track,trackParam,trackPar,extrapolator,hypo))
            messageDebug("WARNING: Problems encountered adding point(s) between track parameters in extending to Volume");
        }
      }
    } else {
      messageDebug("WARNING: Problems encountered getting boundary surfaces from Volume");
    }
  }
  if (VP1Msg::verbose())
    messageVerbose("makePointsCharged end with "+str(points.size())+"points");
  return true;
}

makePointsNeutral()

bool TrackPropagationHelper::makePointsNeutral	(	std::vector< Amg::Vector3D > &	points,
		const Trk::Track *	track )

Definition at line 152 of file TrackPropagationHelper.cxx.

{
  if (VP1Msg::verbose())
    messageVerbose("makePointsNeutral start");
 
  points.clear();
  if (!track) {
    message("makePointsNeutral called with null track pointer");
    return false;
  }
  const unsigned npars = track->trackParameters()->size();
  if (npars==0) {
    message("makePointsNeutral Error: No TrackParameters on track!");
    return false;
  }
  if (npars==1)
    return m_d->makePointsNeutral_SinglePar(points,track);
 
  points.reserve(track->trackParameters()->size());
 
 
  Trk::TrackStates::const_iterator tsos_iter = track->trackStateOnSurfaces()->begin();
  Trk::TrackStates::const_iterator tsos_end = track->trackStateOnSurfaces()->end();
  bool problems(false);
  for (; tsos_iter != tsos_end; ++tsos_iter) {
    if (!m_d->tracksanity->isSafe(*tsos_iter)) {
      problems = true;
      continue;
    }
    //if ((*tsos_iter)->materialEffectsOnTrack())
    //  continue;
    const Trk::TrackParameters* trackParam = (*tsos_iter)->trackParameters();
    if (!m_d->tracksanity->isSafe(trackParam)) {
      problems = true;
      continue;
    }
    points.push_back(trackParam->position());
  }
  if (problems) {
    messageDebug("makePointsNeutral WARNING: Track had one or more track parameters which is unsafe to use in job.");
    if (points.size()<2) {
      messageDebug("makePointsNeutral ERROR: Track did not have at least two safe parameters.");
      points.clear();
      return false;
    }
  }
 
  if (VP1Msg::verbose())
    messageVerbose("makePointsNeutral_SinglePar end");
  return true;
}

message() [1/3]

void VP1HelperClassBase::message ( const QString & str ) const

inherited

Definition at line 49 of file VP1HelperClassBase.cxx.

{
  if (m_helpername.isEmpty()) {
    if (m_system)
      m_system->message(str);
    else
      std::cout<<VP1Msg::prefix_msg()<<" [nameless helper class]: "<<str.toStdString()<<std::endl;
  } else {
    if (m_system)
      m_system->message("["+m_helpername+"] " + str);
    else
      std::cout<<VP1Msg::prefix_msg()<<" ["<<m_helpername.toStdString()<<"]: "<<str.toStdString()<<std::endl;
  }
}

message() [2/3]

void VP1HelperClassBase::message ( const QString & addtostart, const QStringList & l, const QString & addtoend = "" ) const

inherited

Definition at line 131 of file VP1HelperClassBase.cxx.

{
  if (addtostart.isEmpty()) {
    message(l,addtoend);
    return;
  }
  if (addtoend.isEmpty()) {
    for (const QString& s :  l)
      message(addtostart+s);
  } else {
    for (const QString& s :  l)
      message(addtostart+s+addtoend);
  }
}

message() [3/3]

void VP1HelperClassBase::message ( const QStringList & l, const QString & addtoend = "" ) const

inherited

Definition at line 91 of file VP1HelperClassBase.cxx.

{
  if (addtoend.isEmpty()) {
    for (const QString& s :  l)
      message(s);
  } else {
    for (const QString& s :  l)
      message(s+addtoend);
  }
}

messageDebug() [1/3]

void VP1HelperClassBase::messageDebug ( const QString & str ) const

inherited

Definition at line 65 of file VP1HelperClassBase.cxx.

{  
  if (!VP1Msg::debug())
    return;
  std::string sysstring(m_system ? " in "+m_system->name().toStdString() : std::string(""));
  if (m_helpername.isEmpty()) {
    std::cout<<VP1Msg::prefix_debug()<<" [helper"<<sysstring<<"]: "<<str.toStdString()<<std::endl;
  } else {
    std::cout<<VP1Msg::prefix_debug()<<" ["<<m_helpername.toStdString()<<sysstring<<"]: "<<str.toStdString()<<std::endl;
  }
}

messageDebug() [2/3]

void VP1HelperClassBase::messageDebug ( const QString & addtostart, const QStringList & l, const QString & addtoend = "" ) const

inherited

Definition at line 147 of file VP1HelperClassBase.cxx.

{
  if (!VP1Msg::debug())
    return;
  if (addtostart.isEmpty()) {
    messageDebug(l,addtoend);
    return;
  }
  if (addtoend.isEmpty()) {
    for (const QString& s :  l)
      messageDebug(addtostart+s);
  } else {
    for (const QString& s :  l)
      messageDebug(addtostart+s+addtoend);
  }
}

messageDebug() [3/3]

void VP1HelperClassBase::messageDebug ( const QStringList & l, const QString & addtoend = "" ) const

inherited

Definition at line 103 of file VP1HelperClassBase.cxx.

{
  if (!VP1Msg::debug())
    return;
  if (addtoend.isEmpty()) {
    for (const QString& s :  l)
      messageDebug(s);
  } else {
    for (const QString& s :  l)
      messageDebug(s+addtoend);
  }
}

messageVerbose() [1/3]

void VP1HelperClassBase::messageVerbose ( const QString & str ) const

inherited

Definition at line 78 of file VP1HelperClassBase.cxx.

{
  if (!VP1Msg::verbose())
    return;
  std::string sysstring(m_system ? " in "+m_system->name().toStdString() : std::string(""));
  if (m_helpername.isEmpty()) {
    std::cout<<VP1Msg::prefix_verbose()<<" [helper"<<sysstring<<"]: "<<str.toStdString()<<std::endl;
  } else {
    std::cout<<VP1Msg::prefix_verbose()<<" ["<<m_helpername.toStdString()<<sysstring<<"]: "<<str.toStdString()<<std::endl;
  }
}

messageVerbose() [2/3]

void VP1HelperClassBase::messageVerbose ( const QString & addtostart, const QStringList & l, const QString & addtoend = "" ) const

inherited

Definition at line 165 of file VP1HelperClassBase.cxx.

{
  if (!VP1Msg::verbose())
    return;
  if (addtostart.isEmpty()) {
    messageVerbose(l,addtoend);
    return;
  }
  if (addtoend.isEmpty()) {
    for (const QString& s :  l)
      messageVerbose(addtostart+s);
  } else {
    for (const QString& s :  l)
      messageVerbose(addtostart+s+addtoend);
  }
}

messageVerbose() [3/3]

void VP1HelperClassBase::messageVerbose ( const QStringList & l, const QString & addtoend = "" ) const

inherited

Definition at line 117 of file VP1HelperClassBase.cxx.

{
  if (!VP1Msg::verbose())
    return;
  if (addtoend.isEmpty()) {
    for (const QString& s :  l)
      messageVerbose(s);
  } else {
    for (const QString& s :  l)
      messageVerbose(s+addtoend);
  }
}

setHelperClassName()

void VP1HelperClassBase::setHelperClassName ( const QString & n )

inlineprotectedinherited

Definition at line 59 of file VP1HelperClassBase.h.

{ m_helpername = n; }

setSystemBasePointer()

void VP1HelperClassBase::setSystemBasePointer ( IVP1System * sys )

inherited

Definition at line 197 of file VP1HelperClassBase.cxx.

{
  m_system = sys;
  if (VP1Msg::verbose()){
    vp1helperclassbase_instanceMap[this] = (m_system?m_system->name():QString(""));
  }
}

showExtrapolationSurfaces()

bool TrackPropagationHelper::showExtrapolationSurfaces

(

)

const

Definition at line 379 of file TrackPropagationHelper.cxx.

                                                             {
  return m_d->showExtrapSurfaces;
}

str() [1/30]

QString VP1String::str ( const Amg::Vector3D & t )

inlinestaticinherited

Definition at line 98 of file VP1String.h.

{ return "("+str(t.x())+", "+str(t.y())+", "+str(t.z())+")"; }

str() [2/30]

QString VP1String::str ( const bool b )

inlinestaticinherited

Definition at line 53 of file VP1String.h.

{ return b?"True":"False"; }

str() [3/30]

QString VP1String::str ( const char * c )

inlinestaticinherited

Definition at line 50 of file VP1String.h.

{ return c; }

str() [4/30]

QString VP1String::str ( const double & d )

inlinestaticinherited

Definition at line 81 of file VP1String.h.

{ return QString::number(d); }

str() [5/30]

QString VP1String::str ( const float & f )

inlinestaticinherited

Definition at line 82 of file VP1String.h.

{ return QString::number(f); }

str() [6/30]

template<class T>

QString VP1String::str ( const HepGeom::BasicVector3D< T > & t )

inlinestaticinherited

Definition at line 95 of file VP1String.h.

{ return "("+str(t.x())+", "+str(t.y())+", "+str(t.z())+")"; }

str() [7/30]

QString VP1String::str ( const QColor & c )

staticinherited

Definition at line 30 of file VP1String.cxx.

{
  return c.isValid() ? c.name() : "Invalid";
}

str() [8/30]

template<class T>

QString VP1String::str ( const QFlags< T > & f )

inlinestaticinherited

Definition at line 91 of file VP1String.h.

{ return "0x"+QString::number(f, 16).toUpper().rightJustified(8,'0'); }

str() [9/30]

template<class T>

QString VP1String::str ( const QList< T > & t )

inlinestaticinherited

Definition at line 102 of file VP1String.h.

{  return "QList of size"+QString::number(t.size()); }

str() [10/30]

QString VP1String::str ( const QString & s )

inlinestaticinherited

Definition at line 49 of file VP1String.h.

{ return s; }

str() [11/30]

QString VP1String::str ( const SbColor & c )

staticinherited

Definition at line 36 of file VP1String.cxx.

{
  return str(VP1QtInventorUtils::sbcol2qcol(c));
}

str() [12/30]

QString VP1String::str ( const SbVec2d & v )

staticinherited

Definition at line 61 of file VP1String.cxx.

{ double x,y; v.getValue(x,y); return "("+str(x)+", "+str(y)+")"; }

str() [13/30]

QString VP1String::str ( const SbVec2f & v )

staticinherited

Definition at line 62 of file VP1String.cxx.

{ float  x,y; v.getValue(x,y); return "("+str(x)+", "+str(y)+")"; }

str() [14/30]

QString VP1String::str ( const SbVec2s & v )

staticinherited

Definition at line 63 of file VP1String.cxx.

{ short  x,y; v.getValue(x,y); return "("+str(x)+", "+str(y)+")"; }

str() [15/30]

QString VP1String::str ( const SbVec3d & v )

staticinherited

Definition at line 64 of file VP1String.cxx.

{ double x,y,z; v.getValue(x,y,z); return "("+str(x)+", "+str(y)+", "+str(z)+")"; }

str() [16/30]

QString VP1String::str ( const SbVec3f & v )

staticinherited

Definition at line 65 of file VP1String.cxx.

{ float  x,y,z; v.getValue(x,y,z); return "("+str(x)+", "+str(y)+", "+str(z)+")"; }

str() [17/30]

QString VP1String::str ( const SbVec3s & v )

staticinherited

Definition at line 66 of file VP1String.cxx.

{ short  x,y,z; v.getValue(x,y,z); return "("+str(x)+", "+str(y)+", "+str(z)+")"; }

str() [18/30]

QString VP1String::str ( const SbVec4d & v )

staticinherited

Definition at line 67 of file VP1String.cxx.

{ double x,y,z,t; v.getValue(x,y,z,t); return "("+str(x)+", "+str(y)+", "+str(z)+", "+str(t)+")"; }

str() [19/30]

QString VP1String::str ( const SbVec4f & v )

staticinherited

Definition at line 68 of file VP1String.cxx.

{ float  x,y,z,t; v.getValue(x,y,z,t); return "("+str(x)+", "+str(y)+", "+str(z)+", "+str(t)+")"; }

str() [20/30]

template<class T>

QString VP1String::str ( const T * t )

inlinestaticinherited

Definition at line 87 of file VP1String.h.

{ return str(static_cast<const void* >(t)); }

str() [21/30]

QString VP1String::str ( const void * p )

staticinherited

Definition at line 48 of file VP1String.cxx.

{
  if (p) {
    std::ostringstream s;
    s << p;
    // Explicitly naming QString here avoids a cppcheck warning.
    return QString (s.str().c_str());
  } else {
    return "NULL";
  }
}

str() [22/30]

QString VP1String::str ( const VP1Interval & i )

staticinherited

Definition at line 42 of file VP1String.cxx.

{
  return i.toString();
}

str() [23/30]

QString VP1String::str ( int n )

inlinestaticinherited

Definition at line 77 of file VP1String.h.

{ return QString::number(n); }

str() [24/30]

QString VP1String::str ( long n )

inlinestaticinherited

Definition at line 75 of file VP1String.h.

{ return QString::number(n); }

str() [25/30]

QString VP1String::str ( qlonglong n )

inlinestaticinherited

Definition at line 79 of file VP1String.h.

{ return QString::number(n); }

str() [26/30]

QString VP1String::str ( qulonglong n )

inlinestaticinherited

Definition at line 80 of file VP1String.h.

{ return QString::number(n); }

str() [27/30]

QString VP1String::str ( short int n )

inlinestaticinherited

Definition at line 73 of file VP1String.h.

{ return QString::number(n); }

str() [28/30]

QString VP1String::str ( uint n )

inlinestaticinherited

Definition at line 78 of file VP1String.h.

{ return QString::number(n); }

str() [29/30]

QString VP1String::str ( ulong n )

inlinestaticinherited

Definition at line 76 of file VP1String.h.

{ return QString::number(n); }

str() [30/30]

QString VP1String::str ( unsigned short int n )

inlinestaticinherited

Definition at line 74 of file VP1String.h.

{ return QString::number(n); }

systemBase()

IVP1System * VP1HelperClassBase::systemBase ( ) const

inlineinherited

Definition at line 50 of file VP1HelperClassBase.h.

{ return m_system; }

verbose()

bool VP1HelperClassBase::verbose ( )

inlinestaticinherited

Definition at line 32 of file VP1HelperClassBase.h.

{ return s_vp1verbose; } // Returns true if env var VP1_VERBOSE_OUTPUT=1

warnUndeletedInstances()

void VP1HelperClassBase::warnUndeletedInstances ( )

staticinherited

Definition at line 183 of file VP1HelperClassBase.cxx.

{
  if (vp1helperclassbase_instanceMap.empty())
    return;
  std::cout << "WARNING: Detected "<<vp1helperclassbase_instanceMap.size()<<" undeleted helper class instances:"<<std::endl;
 
  std::map<VP1HelperClassBase*,QString>::iterator it,itE(vp1helperclassbase_instanceMap.end());
  for (it = vp1helperclassbase_instanceMap.begin();it!=itE;++it) {
    std::cout << "  ==> "<<it->first<<": "<<it->first->m_helpername.toStdString()
          << (it->second.isEmpty()?QString(""):" (in system "+it->second+")").toStdString()<<std::endl;
  }
}

Member Data Documentation

m_d

Imp* TrackPropagationHelper::m_d

private

Definition at line 58 of file TrackPropagationHelper.h.

m_helpername

QString VP1HelperClassBase::m_helpername

privateinherited

Definition at line 67 of file VP1HelperClassBase.h.

m_system

IVP1System* VP1HelperClassBase::m_system

privateinherited

Definition at line 68 of file VP1HelperClassBase.h.

s_vp1verbose

const bool VP1HelperClassBase::s_vp1verbose = VP1QtUtils::environmentVariableIsOn("VP1_VERBOSE_OUTPUT")

staticprivateinherited

Definition at line 69 of file VP1HelperClassBase.h.

---

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

• TrackPropagationHelper.h
• TrackPropagationHelper.cxx