InDetProjHelper Class Reference

#include <InDetProjHelper.h>

Inheritance diagram for InDetProjHelper:

Collaboration diagram for InDetProjHelper:

Classes
class	Imp

Public Types
enum	PartsFlag {   NoParts = 0x00 , BarrelA = 0x01 , BarrelC = 0x02 , EndCapA = 0x04 ,   EndCapC = 0x08 }

Public Member Functions
virtual	~InDetProjHelper ()
InDetProjFlags::InDetProjPartsFlags	setParts (InDetProjFlags::InDetProjPartsFlags)
InDetProjFlags::InDetProjPartsFlags	parts () const
void	clipPath (const std::vector< Amg::Vector3D > &path, Amg::SetVectorVector3D &resulting_subpaths) const
void	clipPath (const std::vector< Amg::Vector3D > &path, Amg::SetVectorVector3D &resulting_subpaths_barrelA, Amg::SetVectorVector3D &resulting_subpaths_barrelC, Amg::SetVectorVector3D &resulting_subpaths_endcapA, Amg::SetVectorVector3D &resulting_subpaths_endcapC) const
void	projectPath (const std::vector< Amg::Vector3D > &path, Amg::SetVectorVector3D &resulting_projections) const
void	projectPath (const std::vector< Amg::Vector3D > &path, Amg::SetVectorVector3D &resulting_projections_barrelA, Amg::SetVectorVector3D &resulting_projections_barrelC, Amg::SetVectorVector3D &resulting_projections_endcapA, Amg::SetVectorVector3D &resulting_projections_endcapC) const
PartsFlags	touchedParts (const std::vector< Amg::Vector3D > &path) 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 InDetProjHelper *	createPixelHelper (IVP1System *sys=0)
static InDetProjHelper *	createSCTHelper (IVP1System *sys=0)
static InDetProjHelper *	createTRTHelper (IVP1System *sys=0)
static void	transformECPointToZPlane_specialZtoR (Amg::Vector3D &p, const double &planeZ, const double &planeRBegin, const double &endcapZBegin, const double &squeezeFactor)
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 Member Functions
	InDetProjHelper (double surfacethickness, double data_disttosurface_epsilon, double barrel_inner_radius, double barrel_outer_radius, double barrel_posneg_z, double endcap_surface_z, double endcap_surface_length, double endcap_inner_radius, double endcap_outer_radius, double endcap_zasr_innerradius, double endcap_zasr_endcapz_begin, double endcap_zasr_squeezefact, IVP1System *sys)

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 34 of file InDetProjHelper.h.

Member Enumeration Documentation

PartsFlag

enum InDetProjHelper::PartsFlag

Enumerator
NoParts
BarrelA
BarrelC
EndCapA
EndCapC

Definition at line 66 of file InDetProjHelper.h.

{ NoParts = 0x00, BarrelA = 0x01, BarrelC = 0x02, EndCapA = 0x04, EndCapC = 0x08 };

Constructor & Destructor Documentation

~InDetProjHelper()

InDetProjHelper::~InDetProjHelper ( )

virtual

Definition at line 190 of file InDetProjHelper.cxx.

{
  delete m_d;
}

InDetProjHelper()

InDetProjHelper::InDetProjHelper ( double surfacethickness, double data_disttosurface_epsilon, double barrel_inner_radius, double barrel_outer_radius, double barrel_posneg_z, double endcap_surface_z, double endcap_surface_length, double endcap_inner_radius, double endcap_outer_radius, double endcap_zasr_innerradius, double endcap_zasr_endcapz_begin, double endcap_zasr_squeezefact, IVP1System * sys )

private

Definition at line 151 of file InDetProjHelper.cxx.

  : VP1HelperClassBase(sys,"InDetProjHelper"), m_d(new Imp)
{
  m_d->theclass = this;
 
  m_d->surfacethickness = surfacethickness;
  m_d->data_disttosurface_epsilon = data_disttosurface_epsilon;
  m_d->barrel_inner_radius = barrel_inner_radius;
  m_d->barrel_outer_radius = barrel_outer_radius;
  m_d->barrel_posneg_z = barrel_posneg_z;
  m_d->endcap_surface_z = endcap_surface_z;
  m_d->endcap_surface_length = endcap_surface_length;
  m_d->endcap_inner_radius = endcap_inner_radius;
  m_d->endcap_outer_radius = endcap_outer_radius;
  m_d->endcap_zasr_innerradius = endcap_zasr_innerradius;
  m_d->endcap_zasr_endcapz_begin = endcap_zasr_endcapz_begin;
  m_d->endcap_zasr_squeezefact = endcap_zasr_squeezefact;
 
  m_d->parts = InDetProjFlags::NoProjections;
  m_d->covercyl_zmin = 0.0;
  m_d->covercyl_zmax = 0.0;
  m_d->covercyl_rmin = 0.0;
  m_d->covercyl_rmax = 0.0;
 
}

Member Function Documentation

clipPath() [1/2]

void InDetProjHelper::clipPath	(	const std::vector< Amg::Vector3D > &	path,
		Amg::SetVectorVector3D &	resulting_subpaths ) const

Definition at line 274 of file InDetProjHelper.cxx.

{
  clipPath(path,resulting_subpaths,resulting_subpaths,resulting_subpaths,resulting_subpaths);
}

clipPath() [2/2]

void InDetProjHelper::clipPath	(	const std::vector< Amg::Vector3D > &	path,
		Amg::SetVectorVector3D &	resulting_subpaths_barrelA,
		Amg::SetVectorVector3D &	resulting_subpaths_barrelC,
		Amg::SetVectorVector3D &	resulting_subpaths_endcapA,
		Amg::SetVectorVector3D &	resulting_subpaths_endcapC ) const

Definition at line 281 of file InDetProjHelper.cxx.

{
  if (VP1Msg::verbose())
    messageVerbose("clipPath(..) called. Input path has "+QString::number(path.size())+" points.");
 
  resulting_subpaths_barrelA.clear();
  resulting_subpaths_barrelC.clear();
  resulting_subpaths_endcapA.clear();
  resulting_subpaths_endcapC.clear();
 
  //Fixme: If verbose - perform sanity check of input data (check for NAN's).
  if (m_d->parts == InDetProjFlags::NoProjections ) {
    if (VP1Msg::verbose())
      messageVerbose("All projections currently off.");
    return;
  }
  if ( path.size()<2 ) {
    if (VP1Msg::verbose())
      messageVerbose("Input path too short.");
    return;
  }
 
  // Find the clipped path's in all of the enabled detector parts.
 
  //For efficiency, we first clip the path to the smallest
  //axis-aligned cylinder containing all of the projective volumes
  Amg::SetVectorVector3D paths_clipped;
  m_d->clipPathToHollowCylinder( path, paths_clipped,
                   m_d->covercyl_rmin, m_d->covercyl_rmax,
                   m_d->covercyl_zmin, m_d->covercyl_zmax );
 
  if (paths_clipped.empty()) {
    if (VP1Msg::verbose())
      messageVerbose("Path entirely outside clip volumes.");
    return;
  }
 
  const bool enabled_brlA = m_d->parts & InDetProjFlags::Barrel_AllPos;
  const bool enabled_brlC = m_d->parts & InDetProjFlags::Barrel_AllNeg;
  const bool enabled_ecA = m_d->parts & InDetProjFlags::EndCap_AllPos;
  const bool enabled_ecC = m_d->parts & InDetProjFlags::EndCap_AllNeg;
 
  //Special case: If exactly one of the four parts is enabled, we already have our result:
  if ( ( (enabled_brlA?1:0) + (enabled_brlC?1:0) + (enabled_ecA?1:0) + (enabled_ecC?1:0) ) == 1 ) {
    if (enabled_brlA) {
      resulting_subpaths_barrelA = paths_clipped;
      if (VP1Msg::verbose())
    messageVerbose("clipPath(..) only brlA enabled. Returning.");
      return;
    }
    if (enabled_brlC) {
      resulting_subpaths_barrelC = paths_clipped;
      if (VP1Msg::verbose())
    messageVerbose("clipPath(..) only brlC enabled. Returning.");
      return;
    }
    if (enabled_ecA) {
      resulting_subpaths_endcapA = paths_clipped;
      if (VP1Msg::verbose())
    messageVerbose("clipPath(..) only ecA enabled. Returning.");
      return;
    }
    if (enabled_ecC) {
      resulting_subpaths_endcapC = paths_clipped;
      if (VP1Msg::verbose())
    messageVerbose("clipPath(..) only ecC enabled. Returning.");
      return;
    }
  }
 
 
  //For each of the segments, we then find its clipped parts inside
  //the four detector volumes: BarrelA, BarrelC, EndCapA, EndCapC.
  //   Amg::SetVectorVector3D paths_brlA, paths_brlC, paths_ecA,paths_ecC;
  Amg::SetVectorVector3D::const_iterator it, itE(paths_clipped.end());
  for (it = paths_clipped.begin();it!=itE;++it) {
    if ( enabled_brlA )
      m_d->clipPathToHollowCylinder( *it, resulting_subpaths_barrelA, m_d->barrel_inner_radius, m_d->barrel_outer_radius, 0, m_d->barrel_posneg_z );
    if ( enabled_brlC )
      m_d->clipPathToHollowCylinder( *it, resulting_subpaths_barrelC, m_d->barrel_inner_radius, m_d->barrel_outer_radius, - m_d->barrel_posneg_z, 0 );
    if ( enabled_ecA )
      m_d->clipPathToHollowCylinder( *it, resulting_subpaths_endcapA, m_d->endcap_inner_radius, m_d->endcap_outer_radius,
                 m_d->endcap_surface_z - m_d->endcap_surface_length * 0.5, m_d->endcap_surface_z + m_d->endcap_surface_length * 0.5 );
    if ( enabled_ecC )
      m_d->clipPathToHollowCylinder( *it, resulting_subpaths_endcapC, m_d->endcap_inner_radius, m_d->endcap_outer_radius,
                 - m_d->endcap_surface_z - m_d->endcap_surface_length * 0.5, - m_d->endcap_surface_z + m_d->endcap_surface_length * 0.5 );
  }
 
  messageVerbose("clipPath(..) end.");
  //Fixme: If verbose: sanity check on output!
}

createPixelHelper()

InDetProjHelper * InDetProjHelper::createPixelHelper ( IVP1System * sys = 0 )

static

Definition at line 32 of file InDetProjHelper.cxx.

{
  return new InDetProjHelper( InDetProjParams::surfacethickness(),
                  InDetProjParams::pixel_data_disttosurface_epsilon(),
                  InDetProjParams::pixel_barrel_inner_radius(),
                  InDetProjParams::pixel_barrel_outer_radius(),
                  InDetProjParams::pixel_barrel_posneg_z(),
                  InDetProjParams::pixel_endcap_surface_z(),
                  InDetProjParams::pixel_endcap_surface_length(),
                  InDetProjParams::pixel_endcap_inner_radius(),
                  InDetProjParams::pixel_endcap_outer_radius(),
                  InDetProjParams::pixel_endcap_zasr_innerradius(),
                  InDetProjParams::pixel_endcap_zasr_endcapz_begin(),
                  InDetProjParams::pixel_endcap_zasr_squeezefact(),
                  system );
}

createSCTHelper()

InDetProjHelper * InDetProjHelper::createSCTHelper ( IVP1System * sys = 0 )

static

Definition at line 50 of file InDetProjHelper.cxx.

{
  return new InDetProjHelper( InDetProjParams::surfacethickness(),
                  InDetProjParams::sct_data_disttosurface_epsilon(),
                  InDetProjParams::sct_barrel_inner_radius(),
                  InDetProjParams::sct_barrel_outer_radius(),
                  InDetProjParams::sct_barrel_posneg_z(),
                  InDetProjParams::sct_endcap_surface_z(),
                  InDetProjParams::sct_endcap_surface_length(),
                  InDetProjParams::sct_endcap_inner_radius(),
                  InDetProjParams::sct_endcap_outer_radius(),
                  InDetProjParams::sct_endcap_zasr_innerradius(),
                  InDetProjParams::sct_endcap_zasr_endcapz_begin(),
                  InDetProjParams::sct_endcap_zasr_squeezefact(),
                  system );
}

createTRTHelper()

InDetProjHelper * InDetProjHelper::createTRTHelper ( IVP1System * sys = 0 )

static

Definition at line 68 of file InDetProjHelper.cxx.

{
  return new InDetProjHelper( InDetProjParams::surfacethickness(),
                  InDetProjParams::trt_data_disttosurface_epsilon(),
                  InDetProjParams::trt_barrel_inner_radius(),
                  InDetProjParams::trt_barrel_outer_radius(),
                  InDetProjParams::trt_barrel_posneg_z(),
                  InDetProjParams::trt_endcap_surface_z(),
                  InDetProjParams::trt_endcap_surface_length(),
                  InDetProjParams::trt_endcap_inner_radius(),
                  InDetProjParams::trt_endcap_outer_radius(),
                  InDetProjParams::trt_endcap_zasr_innerradius(),
                  InDetProjParams::trt_endcap_zasr_endcapz_begin(),
                  InDetProjParams::trt_endcap_zasr_squeezefact(),
                  system );
}

helperClassName()

const QString & VP1HelperClassBase::helperClassName ( ) const

inlineinherited

Definition at line 51 of file VP1HelperClassBase.h.

{ return m_helpername; }

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);
  }
}

parts()

InDetProjFlags::InDetProjPartsFlags InDetProjHelper::parts

(

)

const

Definition at line 268 of file InDetProjHelper.cxx.

{
  return m_d->parts;
}

projectPath() [1/2]

void InDetProjHelper::projectPath	(	const std::vector< Amg::Vector3D > &	path,
		Amg::SetVectorVector3D &	resulting_projections ) const

Definition at line 819 of file InDetProjHelper.cxx.

{
  projectPath(path,resulting_projs,resulting_projs,resulting_projs,resulting_projs);
}

projectPath() [2/2]

void InDetProjHelper::projectPath	(	const std::vector< Amg::Vector3D > &	path,
		Amg::SetVectorVector3D &	resulting_projections_barrelA,
		Amg::SetVectorVector3D &	resulting_projections_barrelC,
		Amg::SetVectorVector3D &	resulting_projections_endcapA,
		Amg::SetVectorVector3D &	resulting_projections_endcapC ) const

Definition at line 826 of file InDetProjHelper.cxx.

{
  if (VP1Msg::verbose())
    messageVerbose("projectPath(..) called. Input path has "+QString::number(path.size())+" points.");
 
  resulting_projections_barrelA.clear();
  resulting_projections_barrelC.clear();
  resulting_projections_endcapA.clear();
  resulting_projections_endcapC.clear();
 
  //Fixme: If verbose - perform sanity check of input data (check for NAN's).
  if (m_d->parts == InDetProjFlags::NoProjections ) {
    if (VP1Msg::verbose())
      messageVerbose("All projections currently off.");
    return;
  }
  if ( path.size()<2 ) {
    if (VP1Msg::verbose())
      messageVerbose("Input path too short.");
    return;
  }
 
  // ===> First we must find the clipped path's in all of the enabled detector parts.
 
  Amg::SetVectorVector3D paths_brlA, paths_brlC, paths_ecA,paths_ecC;
  clipPath( path,paths_brlA, paths_brlC, paths_ecA,paths_ecC);
 
  // ===> Then we project those.
 
  //Fixme: The dependence on surface thickness and epsilon below is very preliminary.
 
  const double eps = m_d->data_disttosurface_epsilon;
  const double endcapeps(-5*SYSTEM_OF_UNITS::mm);//fixme hardcoding..
 
  Amg::SetVectorVector3D::const_iterator it,itE;
 
  if (m_d->parts & InDetProjFlags::Barrel_AllPos) {
    itE = paths_brlA.end();
    if ( m_d->parts & InDetProjFlags::BarrelCentral )
      for ( it = paths_brlA.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane( *it, resulting_projections_barrelA, 0.5*m_d->surfacethickness+eps );
    if ( m_d->parts & InDetProjFlags::BarrelPositive )
      for ( it = paths_brlA.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane( *it, resulting_projections_barrelA, m_d->barrel_posneg_z - eps );
  }
  if ( m_d->parts & InDetProjFlags::Barrel_AllNeg ) {
    itE = paths_brlC.end();
    if ( m_d->parts & InDetProjFlags::BarrelCentral )
      for ( it = paths_brlC.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane( *it, resulting_projections_barrelC, - 0.5*m_d->surfacethickness - eps);
    if ( m_d->parts & InDetProjFlags::BarrelNegative )
      for ( it = paths_brlC.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane( *it, resulting_projections_barrelC, - m_d->barrel_posneg_z );
  }
  if ( m_d->parts & InDetProjFlags::EndCap_AllPos ) {
    itE = paths_ecA.end();
    if ( m_d->parts & InDetProjFlags::EndCapInnerPositive )
      for ( it = paths_ecA.begin(); it!=itE; ++it )
    m_d->projectPathToInfiniteCylinder( *it, resulting_projections_endcapA, m_d->endcap_inner_radius + eps+endcapeps );
    if ( m_d->parts & InDetProjFlags::EndCapOuterPositive )
      for ( it = paths_ecA.begin(); it!=itE; ++it )
    m_d->projectPathToInfiniteCylinder( *it, resulting_projections_endcapA, m_d->endcap_outer_radius + eps+endcapeps );
    //Fixme: Make sure to use the same parameters here as in PRDHandle_TRT.cxx:
    if ( m_d->parts & InDetProjFlags::TRT_EndCapZToRCentral )
      for ( it = paths_ecA.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane_specialZtoR( *it, resulting_projections_endcapA,
                        0.5*m_d->surfacethickness + eps  );
    //Fixme: Make sure to use the same parameters here as in PRDHandle_TRT.cxx:
    if ( m_d->parts & InDetProjFlags::TRT_EndCapZToRPositive )
      for ( it = paths_ecA.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane_specialZtoR( *it, resulting_projections_endcapA,
                        m_d->barrel_posneg_z - 0.5*m_d->surfacethickness - eps /*fixme: +- epsilon??*/  );
  }
  if ( m_d->parts & InDetProjFlags::EndCap_AllNeg ) {
    itE = paths_ecC.end();
    if ( m_d->parts & InDetProjFlags::EndCapInnerNegative )
      for ( it = paths_ecC.begin(); it!=itE; ++it )
    m_d->projectPathToInfiniteCylinder( *it, resulting_projections_endcapC, m_d->endcap_inner_radius + eps+endcapeps );
    if ( m_d->parts & InDetProjFlags::EndCapOuterNegative )
      for ( it = paths_ecC.begin(); it!=itE; ++it )
    m_d->projectPathToInfiniteCylinder( *it, resulting_projections_endcapC, m_d->endcap_outer_radius + eps+endcapeps );
    //Fixme: Make sure to use the same parameters here as in PRDHandle_TRT.cxx:
    if ( m_d->parts & InDetProjFlags::TRT_EndCapZToRCentral )
      for ( it = paths_ecC.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane_specialZtoR( *it, resulting_projections_endcapC,
                        - 0.5*m_d->surfacethickness - eps  );
    //Fixme: Make sure to use the same parameters here as in PRDHandle_TRT.cxx:
    if ( m_d->parts & InDetProjFlags::TRT_EndCapZToRNegative )
      for ( it = paths_ecC.begin(); it!=itE; ++it )
    m_d->projectPathToZPlane_specialZtoR( *it, resulting_projections_endcapC,
                        - m_d->barrel_posneg_z + 0.5*m_d->surfacethickness + eps/*fixme: +- epsilon??*/  );
  }
 
}

setHelperClassName()

void VP1HelperClassBase::setHelperClassName ( const QString & n )

inlineprotectedinherited

Definition at line 59 of file VP1HelperClassBase.h.

{ m_helpername = n; }

setParts()

InDetProjFlags::InDetProjPartsFlags InDetProjHelper::setParts

(

InDetProjFlags::InDetProjPartsFlags

newparts

)

Definition at line 196 of file InDetProjHelper.cxx.

{
  if ( m_d->parts==newparts )
    return m_d->parts;
  InDetProjFlags::InDetProjPartsFlags oldparts = m_d->parts;
  m_d->parts = newparts;
 
  //Update parameters of smallest cylinder covering all enabled clip volumes.
  if (m_d->parts == InDetProjFlags::NoProjections) {
    m_d->covercyl_zmin = 0.0;
    m_d->covercyl_zmax = 0.0;
    m_d->covercyl_rmin = 0.0;
    m_d->covercyl_rmax = 0.0;
    return oldparts;
  }
 
  bool no_ec_neg = !( m_d->parts & InDetProjFlags::EndCap_AllNeg );
  bool no_ec_pos = !( m_d->parts & InDetProjFlags::EndCap_AllPos );
  bool no_brl_neg = !( m_d->parts & InDetProjFlags::Barrel_AllNeg );
  bool no_brl_pos = !( m_d->parts & InDetProjFlags::Barrel_AllPos );
  bool barrel = m_d->parts & InDetProjFlags::Barrel_All;
  bool endcap = m_d->parts & InDetProjFlags::EndCap_All;
 
  m_d->covercyl_zmin = - m_d->endcap_surface_z - 0.5*m_d->endcap_surface_length;
  if ( no_ec_neg ) {
    m_d->covercyl_zmin = - m_d->barrel_posneg_z;
    if ( no_brl_neg ) {
      m_d->covercyl_zmin = 0.0;
      if ( no_brl_pos ) {
    m_d->covercyl_zmin = m_d->barrel_posneg_z;
    if ( no_ec_pos )
      m_d->covercyl_zmin = m_d->endcap_surface_z + 0.5*m_d->endcap_surface_length + 1.0e99;
      }
    }
  }
  m_d->covercyl_zmax = m_d->endcap_surface_z + 0.5*m_d->endcap_surface_length;
  if ( no_ec_pos ) {
    m_d->covercyl_zmax = m_d->barrel_posneg_z;
    if ( no_brl_pos ) {
      m_d->covercyl_zmax = 0.0;
      if ( no_brl_neg ) {
    m_d->covercyl_zmax = - m_d->barrel_posneg_z;
    if ( no_ec_neg )
      m_d->covercyl_zmax = - m_d->endcap_surface_z - 0.5*m_d->endcap_surface_length - 1.0e99;
      }
    }
  }
  if ( m_d->covercyl_zmin >= m_d->covercyl_zmax )
    m_d->covercyl_zmin = m_d->covercyl_zmax = 0;
 
  if ( barrel && endcap ) {
    m_d->covercyl_rmin = std::min(m_d->barrel_inner_radius,m_d->endcap_inner_radius);
    m_d->covercyl_rmax = std::max(m_d->barrel_outer_radius,m_d->endcap_outer_radius);
  } else {
    if (barrel) {
      m_d->covercyl_rmin = m_d->barrel_inner_radius;
      m_d->covercyl_rmax = m_d->barrel_outer_radius;
    } else if (endcap) {
      m_d->covercyl_rmin = m_d->endcap_inner_radius;
      m_d->covercyl_rmax = m_d->endcap_outer_radius;
    } else {
      message("Unforeseen execution path encountered.");
      m_d->covercyl_rmin = 0;
      m_d->covercyl_rmax = 0;
    }
  }
  if ( m_d->covercyl_rmin >= m_d->covercyl_rmax )
    m_d->covercyl_rmin = m_d->covercyl_rmax = 0;
  return oldparts;
}

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(""));
  }
}

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; }

touchedParts()

InDetProjHelper::PartsFlags InDetProjHelper::touchedParts

(

const std::vector< Amg::Vector3D > &

path

)

const

Definition at line 926 of file InDetProjHelper.cxx.

{
  if (VP1Msg::verbose())
    messageVerbose("touchedParts(..) called. Input path has "+QString::number(path.size())+" points.");
  PartsFlags touchedparts = NoParts;
  if ( m_d->touchesHollowCylinder(path,m_d->barrel_inner_radius, m_d->barrel_outer_radius, 0, m_d->barrel_posneg_z) )
    touchedparts |= BarrelA;
  if ( m_d->touchesHollowCylinder(path,m_d->barrel_inner_radius, m_d->barrel_outer_radius, - m_d->barrel_posneg_z, 0) )
    touchedparts |= BarrelC;
  if ( m_d->touchesHollowCylinder(path,m_d->endcap_inner_radius, m_d->endcap_outer_radius,
                m_d->endcap_surface_z - m_d->endcap_surface_length * 0.5, m_d->endcap_surface_z + m_d->endcap_surface_length * 0.5 ) )
    touchedparts |= EndCapA;
  if ( m_d->touchesHollowCylinder(path, m_d->endcap_inner_radius, m_d->endcap_outer_radius,
                - m_d->endcap_surface_z - m_d->endcap_surface_length * 0.5, - m_d->endcap_surface_z + m_d->endcap_surface_length * 0.5) )
    touchedparts |= EndCapC;
  return touchedparts;
}

transformECPointToZPlane_specialZtoR()

void InDetProjHelper::transformECPointToZPlane_specialZtoR ( Amg::Vector3D & p, const double & planeZ, const double & planeRBegin, const double & endcapZBegin, const double & squeezeFactor )

static

Definition at line 780 of file InDetProjHelper.cxx.

{
    if ( p.x()==0.0 && p.y()==0.0 ) {
      VP1Msg::message("InDetProjHelper::transformECPointToZPlane_specialZtoR ERROR: "
              "Point has x==0 and y==0. Ambiguous projection of point.");
//      p.setX(1.0);
      p.x() = 1.0;
    }
    const double r = planeRBegin + (fabs(p.z())-endcapZBegin)/squeezeFactor;
    const double s = r / sqrt( p.x()*p.x()+p.y()*p.y() );
//    p.setX(p.x()*s);
//    p.setY(p.y()*s);
//    p.setZ(planeZ);
    p.x() = p.x()*s;
    p.y() = p.y()*s;
    p.z() = planeZ;
}

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* InDetProjHelper::m_d

private

Definition at line 97 of file InDetProjHelper.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:

• InDetProjHelper.h
• InDetProjHelper.cxx