TRTTrackHoleSearchTool Class Reference

#include <TRTTrackHoleSearchTool.h>

Inheritance diagram for TRTTrackHoleSearchTool:

Collaboration diagram for TRTTrackHoleSearchTool:

Public Member Functions
	TRTTrackHoleSearchTool (const std::string &type, const std::string &name, const IInterface *parent)
StatusCode	initialize ()
StatusCode	finalize ()
void	countHoles (const Trk::Track &track, std::vector< int > &information, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, partHyp Output: Changes in information This method fills the fields relevant to the hole counts in the vector information.
const Trk::TrackStates *	getHolesOnTrack (const Trk::Track &track, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, parthyp Return: A DataVector containing pointers to TrackStateOnSurfaces which each represent an identified hole on the track.
const Trk::Track *	getTrackWithHoles (const Trk::Track &track, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes The parthyp argument is relevant for the extrapolation steps in the hole search.
const Trk::Track *	getTrackWithHolesAndOutliers (const Trk::Track &track, const Trk::ParticleHypothesis partHyp=Trk::pion) const
	Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes or outliers The parthyp argument is relevant for the extrapolation steps in the hole search.
	DeclareInterfaceID (ITrackHoleSearchTool, 1, 0)
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
int	extrapolateBetweenHits (const Trk::TrackParameters *start_parameters, const Trk::Surface &end_surf, Trk::TrackStates *holes, const Trk::ParticleHypothesis partHyp=Trk::pion) const
void	dump_bad_straw_log () const
Trk::TrackStates::const_iterator	find_first_trt_hit (const Trk::TrackStates &track_states) const
Trk::TrackStates::const_iterator	find_last_hit_before_trt (const Trk::TrackStates &track_states) const
const Trk::Track *	addHolesToTrack (const Trk::Track &track, const Trk::TrackStates *holes) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ToolHandle< Trk::IExtrapolator >	m_extrapolator {this, "extrapolator", "Trk::Extrapolator"}
ServiceHandle< IInDetConditionsSvc >	m_conditions_svc {this, "conditions_svc", "TRT_ConditionsSummarySvc"}
BooleanProperty	m_use_conditions_svc {this, "use_conditions_svc", true}
FloatProperty	m_outer_radius {this, "outer_radius", 1075.0*CLHEP::mm}
FloatProperty	m_max_z {this, "max_z", 2715.0*CLHEP::mm}
IntegerProperty	m_max_trailing_holes {this, "max_trailing_holes", 1}
BooleanProperty	m_begin_at_first_trt_hit {this, "begin_at_first_trt_hit", false}
BooleanProperty	m_end_at_last_trt_hit {this, "end_at_last_trt_hit", false}
BooleanProperty	m_bcheck {this, "bcheck", false}
BooleanProperty	m_do_dump_bad_straw_log {this, "do_dump_bad_straw_log", false}
FloatProperty	m_locR_cut {this, "locR_cut", -1.}
FloatProperty	m_locR_sigma_cut {this, "locR_sigma_cut", -1.}
FloatProperty	m_locZ_cut {this, "locZ_cut", 5.0*CLHEP::mm}
std::atomic_bool	m_has_been_called {false}
const TRT_ID *	m_TRT_ID {nullptr}
Trk::CylinderSurface *	m_trt_outer_surf {nullptr}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 37 of file TRTTrackHoleSearchTool.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

TRTTrackHoleSearchTool()

TRTTrackHoleSearchTool::TRTTrackHoleSearchTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Definition at line 39 of file TRTTrackHoleSearchTool.cxx.

    : AthAlgTool(type, name, parent)
{
    declareInterface<ITrackHoleSearchTool>(this);
}

Member Function Documentation

addHolesToTrack()

const Trk::Track * TRTTrackHoleSearchTool::addHolesToTrack ( const Trk::Track & track, const Trk::TrackStates * holes ) const

private

Definition at line 422 of file TRTTrackHoleSearchTool.cxx.

{
  ATH_MSG_DEBUG("TRTTrackHoleSearchTool::addHolesToTrack");
  /*
    This method was basically coppied from here:
    http://alxr.usatlas.bnl.gov/lxr-stb4/source/atlas/InnerDetector/InDetRecTools/InDetTrackHoleSearch/src/InDetTrackHoleSearchTool.cxx#931
  */
 
  // get states from track
  auto tsos = std::make_unique<Trk::TrackStates>();
  for (const auto *it : *track.trackStateOnSurfaces()) {
    // veto old holes
    if (!it->type(Trk::TrackStateOnSurface::Hole)) {
      tsos->push_back(new Trk::TrackStateOnSurface(*it));
    }
  }
 
  // if we have no holes on the old track and no holes found by search, then we
  // just copy the track
  if (track.trackStateOnSurfaces()->size() == tsos->size() && holes->empty()) {
    // create copy of track
    const Trk::Track* new_track = new Trk::Track(
      track.info(),
      std::move(tsos),
      track.fitQuality() ? track.fitQuality()->uniqueClone() : nullptr);
    return new_track;
  }
 
  // add new holes
  tsos->insert(tsos->end(), holes->begin(), holes->end());
 
  // sort
  const Trk::TrackParameters* perigee = track.perigeeParameters();
  if (!perigee)
    perigee = (*(track.trackStateOnSurfaces()->begin()))->trackParameters();
 
  if (perigee) {
    Trk::TrackStateOnSurfaceComparisonFunction CompFunc(perigee->momentum());
 
    if (fabs(perigee->parameters()[Trk::qOverP]) > 0.002) {
      /* invest n*(logN)**2 sorting time for lowPt, coping with a possibly
         not 100% transitive comparison functor. */
      if (msgLvl(MSG::DEBUG)) {
        msg() << "sorting vector with stable_sort" << endmsg;
      }
      std::stable_sort(tsos->begin(), tsos->end(), CompFunc);
    } else {
      tsos->sort(CompFunc); // respects DV object ownership
    }
  }
 
  // create copy of track
  const Trk::Track* new_track =
    new Trk::Track(track.info(),
                   std::move(tsos),
                   track.fitQuality() ? track.fitQuality()->uniqueClone() : nullptr);
 
  return new_track;
}

countHoles()

void TRTTrackHoleSearchTool::countHoles ( const Trk::Track & track, std::vector< int > & information, const Trk::ParticleHypothesis partHyp = Trk::pion ) const

virtual

Input : track, partHyp Output: Changes in information This method fills the fields relevant to the hole counts in the vector information.

These fields should be initialised to 0. The relevant indices are specified by the enumeration in Tracking/TrkEvent/TrkTrackSummary. If problems occur, the information counters are reset to -1 flagging them as not set. The parthyp argument is relevant for the extrapolation steps in the hole search.

Implements Trk::ITrackHoleSearchTool.

Definition at line 76 of file TRTTrackHoleSearchTool.cxx.

                                                                                   {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::countHoles" );
    const Trk::TrackStates* holes = getHolesOnTrack(track, partHyp);
    if (holes) {
        information[Trk::numberOfTRTHoles] = holes->size();
        delete holes;
    }   else {
        information[Trk::numberOfTRTHoles] = -1;
    }
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

DeclareInterfaceID()

Trk::ITrackHoleSearchTool::DeclareInterfaceID ( ITrackHoleSearchTool , 1 , 0  )

inherited

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

dump_bad_straw_log()

void TRTTrackHoleSearchTool::dump_bad_straw_log ( ) const

private

Definition at line 354 of file TRTTrackHoleSearchTool.cxx.

                                                      {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::dump_bad_straw_log" );
    std::ofstream out("TRT_ConditionsSummarySvc_bad_straws.log");
    out << "# id  barrel_ec  phi_module  layer_or_wheel  straw_layer  straw" << std::endl;
    for(std::vector<Identifier>::const_iterator it = m_TRT_ID->straw_layer_begin(); it != m_TRT_ID->straw_layer_end(); ++it) {
        for(int i=0; i<= m_TRT_ID->straw_max(*it); ++i) {
            Identifier id = m_TRT_ID->straw_id(*it, i);
            if(!m_conditions_svc->isGood(id)) {
                out << id.getString()
                    << std::setw(3) << m_TRT_ID->barrel_ec(id)
                    << std::setw(4) << m_TRT_ID->phi_module(id)
                    << std::setw(4) << m_TRT_ID->layer_or_wheel(id)
                    << std::setw(4) << m_TRT_ID->straw_layer(id)
                    << std::setw(4) << m_TRT_ID->straw(id)
                    << std::endl;
            }
        }
    }
    out.close();
}

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

extrapolateBetweenHits()

int TRTTrackHoleSearchTool::extrapolateBetweenHits ( const Trk::TrackParameters * start_parameters, const Trk::Surface & end_surf, Trk::TrackStates * holes, const Trk::ParticleHypothesis partHyp = Trk::pion ) const

private

Definition at line 226 of file TRTTrackHoleSearchTool.cxx.

                                                                                               {
    int hole_count = 0;
    // initialize previous id
    const Trk::Surface& start_surf = start_parameters->associatedSurface();
 
    Identifier previous_id = start_surf.associatedDetectorElementIdentifier();
 
    // get end id
    const Identifier end_id = end_surf.associatedDetectorElementIdentifier();
 
    // look for holes
    std::vector<std::unique_ptr<Trk::TrackParameters> > steps =
          m_extrapolator->extrapolateStepwise(Gaudi::Hive::currentContext(),
                                              *start_parameters, 
                                              end_surf, 
                                              Trk::alongMomentum, m_bcheck.value(),
                                              partHyp);
 
    if(steps.empty()) {
        ATH_MSG_DEBUG("extrapolateBetweenHits: extrapolateStepwise returned null");
    } else {
        // loop over parameters from extrapolation
        // note: the last element in the vector is always the track parameters at the destination surface
                for(std::vector<std::unique_ptr<Trk::TrackParameters> >::const_iterator step = steps.begin(); step != steps.end()-1; ++step) {
            // check for surface
            const Trk::Surface& surf = (*step)->associatedSurface();
            
            // get id
            Identifier id = surf.associatedDetectorElementIdentifier();
 
            // check that the surface is a TRT straw
            if(!m_TRT_ID->is_trt(id)) {
                ATH_MSG_DEBUG("extrapolateBetweenHits: surf is not a TRT straw.  Skipping.");
                continue;
            }
 
            // check that we are not still on the previous id
            if(id == previous_id) {
                ATH_MSG_DEBUG("extrapolateBetweenHits: id == previous_id");
                continue;
            }
 
            // check that we haven't hit the destination surface
            if(id == end_id) {
                ATH_MSG_DEBUG("extrapolateBetweenHits: id == end_id");
                continue;
            }
 
            // ignore id 0xffffffff from m_trt_outer_surf
            if(id == 0xffffffff) {
                ATH_MSG_DEBUG("extrapolateBetweenHits: id == 0xffffffff. Skipping.");
                continue;
            }
 
            // don't count bad straws as holes
            if( (m_use_conditions_svc)&&(!m_conditions_svc->isGood(id)) ) {
                ATH_MSG_DEBUG("extrapolateBetweenHits: ConditionsSvc says this straw is bad. Skipping.");
                continue;
            }
 
            // check fiducial |locR| by sigma
            const float locR = (*step)->parameters()[Trk::locR];
            if(m_locR_sigma_cut > 0.0) {
                //const Trk::MeasuredTrackParameters* meas = dynamic_cast< const Trk::MeasuredTrackParameters* >(*step);
                //if(meas)
                //{
 
                const AmgSymMatrix(5)* merr = (*step)->covariance();
 
                if(merr){
                    //const float locR_error = (*merr)(Trk::locR,Trk::locR);
                    
                    const float locR_error = Amg::error(*merr,Trk::locR);
 
                    // hole must be within the straw by some tolerance
                    if( (2.0*CLHEP::mm - fabs(locR))/locR_error < m_locR_sigma_cut ) {
                        continue;
                    }
                } else {
                    ATH_MSG_WARNING("extrapolateBetweenHits: Track parameters failed to dynamic_cast< const Trk::MeasuredTrackParameters* >.");
                    continue;
                }
            }
            // check fiducial |locR| by value
            if(m_locR_cut > 0.0) {
                if( fabs(locR) > m_locR_cut )
                    continue;
            }
            // check fiducial |locZ| by value
            if(m_locZ_cut > 0.0) {
                const float locZ = (*step)->parameters()[Trk::locZ];
                const Trk::CylinderBounds* cylb = dynamic_cast<const Trk::CylinderBounds *>(&surf.bounds());
                if(cylb) {
                    const float halfz=cylb->halflengthZ();
                    if( fabs(locZ) > halfz - m_locZ_cut ) continue;
                } else {
                    ATH_MSG_WARNING("extrapolateBetweenHits: Surface failed to dynamic_cast to Trk::CylinderBounds. Skipping.");
                    continue;
                }
            }
 
            // if we've gotten here, it is a legitimate hole
            ATH_MSG_DEBUG("      HOLE Found! Identifier = " << id.getString() << " indicates:\n"
                          << "            is_pixel =       " << m_TRT_ID->is_pixel(id) << "\n"
                          << "            is_sct =         " << m_TRT_ID->is_sct(id) << "\n"
                          << "            is_trt =         " << m_TRT_ID->is_trt(id) << "\n"
                          << "            barrel_ec =      " << m_TRT_ID->barrel_ec(id) << "\n"
                          << "            phi_module =     " << m_TRT_ID->phi_module(id) << "\n"
                          << "            layer_or_wheel = " << m_TRT_ID->layer_or_wheel(id) << "\n"
                          << "            straw_layer =    " << m_TRT_ID->straw_layer(id) << "\n"
                          << "            straw =          " << m_TRT_ID->straw(id) << "\n"
                          << "            rho = " << (*step)->position().perp() );
            
            std::bitset<Trk::TrackStateOnSurface::NumberOfTrackStateOnSurfaceTypes> typePattern;
            typePattern.set(Trk::TrackStateOnSurface::Hole);
            holes->push_back( new Trk::TrackStateOnSurface(nullptr, (*step)->uniqueClone(), nullptr, typePattern) );
            hole_count++;
            previous_id = id;
        } // end loop over parameters from extrapolation
    }
    return hole_count;
}

finalize()

StatusCode TRTTrackHoleSearchTool::finalize

(

)

Definition at line 68 of file TRTTrackHoleSearchTool.cxx.

                                            {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::finalize()" );
    delete m_trt_outer_surf;
    return StatusCode::SUCCESS;
}

find_first_trt_hit()

Trk::TrackStates::const_iterator TRTTrackHoleSearchTool::find_first_trt_hit ( const Trk::TrackStates & track_states ) const

private

Definition at line 378 of file TRTTrackHoleSearchTool.cxx.

                                                                                   {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::find_first_trt_hit" );
    // finds the first TRT hit or returns end
    Trk::TrackStates::const_iterator track_state = track_states.begin();
    for(; track_state != track_states.end(); ++track_state) {
        // skip track states that are not measurements
        if(!(*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
            ATH_MSG_VERBOSE( "    TrackStateOnSurface is not of type Trk::TrackStateOnSurface::Measurement." );
            continue;
        }
        if( dynamic_cast<const InDet::TRT_DriftCircleOnTrack*>( (*track_state)->measurementOnTrack() )  ) {
            ATH_MSG_DEBUG("    Found first TRT hit.");
            break;
        }
    }
    return track_state;
}

find_last_hit_before_trt()

Trk::TrackStates::const_iterator TRTTrackHoleSearchTool::find_last_hit_before_trt ( const Trk::TrackStates & track_states ) const

private

Definition at line 399 of file TRTTrackHoleSearchTool.cxx.

                                                                                         {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::find_last_hit_before_trt" );
    // start at first TRT hit
    Trk::TrackStates::const_iterator track_state = find_first_trt_hit(track_states);
    // if there is no room to step backwards in the hits, return end
    if(track_states.size() < 2 || track_state == track_states.begin()) {
        return track_states.end();
    }
    --track_state; // step back and look for last measurement before the TRT hit
    for(; track_state != track_states.begin(); --track_state) {
        if((*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
            break;
        }
    }
    // if the loop ended without finding a measurement, return end
    if(!(*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
        return track_states.end();
    }
    return track_state;
}

getHolesOnTrack()

const Trk::TrackStates * TRTTrackHoleSearchTool::getHolesOnTrack ( const Trk::Track & track, const Trk::ParticleHypothesis partHyp = Trk::pion ) const

virtual

Input : track, parthyp Return: A DataVector containing pointers to TrackStateOnSurfaces which each represent an identified hole on the track.

The parthyp argument is relevant for the extrapolation steps in the hole search. Attention: This is a factory, ownership of the return vector is passed to the calling method.

Implements Trk::ITrackHoleSearchTool.

Definition at line 91 of file TRTTrackHoleSearchTool.cxx.

                                                                   {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::getHolesOnTrack" );
 
    // write out list of bad straws for debugging purposes.
    // only dump log on first call of this function.
    if (m_do_dump_bad_straw_log && !m_has_been_called && m_use_conditions_svc) {
        dump_bad_straw_log();
    }
    m_has_been_called = true;
 
    if (track.perigeeParameters()) {
        ATH_MSG_VERBOSE( "  This track has perigee: \n"
                         << "    pT =  " << track.perigeeParameters()->pT()/CLHEP::GeV << " CLHEP::GeV\n"
                         << "    eta = " << track.perigeeParameters()->eta() << "\n"
                         << "    phi = " << track.perigeeParameters()->parameters()[Trk::phi0] << "\n"
                         << "    d0 =  " << track.perigeeParameters()->parameters()[Trk::d0] << "\n"
                         << "    z0 =  " << track.perigeeParameters()->parameters()[Trk::z0] );
    }
 
    // get TrackStateOnSurfaces
    const Trk::TrackStates* track_states = track.trackStateOnSurfaces();
    if (track_states) {
        ATH_MSG_DEBUG( "  This track has " << track_states->size() << " track states on surface." );
    } else {
        ATH_MSG_WARNING( "  This track has null track states on surface. Returning 0." );
        return nullptr;
    }
 
    if (track_states->size() < 2) {
        ATH_MSG_WARNING( "  Fewer than 2 TrackStatesOnSurface. Returning 0." );
        return nullptr;
    }
 
    // set beginning point of extrapolation
    Trk::TrackStates::const_iterator beginning_track_state;
    if (m_begin_at_first_trt_hit) {
        // begin at first TRT hit
        beginning_track_state = find_first_trt_hit(*track_states);
    } else {
        // begin at last Si hit
        beginning_track_state = find_last_hit_before_trt(*track_states);
        if (beginning_track_state == track_states->end()) {
            ATH_MSG_WARNING( "  beginning_track_state == track_states->end().  There must be no Si hits.\n"
                             << "  Will try to begin at the first TRT hit." );
            beginning_track_state = find_first_trt_hit(*track_states);
        }
    }
 
    if (beginning_track_state == track_states->end()) {
        ATH_MSG_WARNING( "  beginning_track_state == track_states->end(). No where to extrapolate to. Returning 0." );
        return nullptr;
    }
 
    // to be returned:
    Trk::TrackStates* holes = new Trk::TrackStates;
 
    Trk::TrackStates::const_iterator track_state = beginning_track_state;
    const Trk::TrackParameters* start_parameters = (*track_state)->trackParameters();
    ++track_state;
    // loop over TrackStateOnSurfaces (destination surfaces for extrapolation)
    for (; track_state != track_states->end(); ++track_state) {
        // skip track states that are not measurements
        if (!(*track_state)->type(Trk::TrackStateOnSurface::Measurement)) {
            ATH_MSG_VERBOSE( "    TrackStateOnSurface is not of type Trk::TrackStateOnSurface::Measurement." );
            continue;
        }
 
        const Trk::TrackParameters* end_parameters = (*track_state)->trackParameters();
        if(end_parameters) {
            ATH_MSG_VERBOSE( "    TrackStateOnSurface: (x, y, z) = (" 
                             << end_parameters->position().x() << ", "
                             << end_parameters->position().y() << ", "
                             << end_parameters->position().z() << "); (rho, eta, phi) = (" 
                             << end_parameters->position().perp() << ", "
                             << end_parameters->position().eta() << ", "
                             << end_parameters->position().phi() << ")");
        } else {
            ATH_MSG_WARNING( "    TrackStateOnSurface has no TrackParameters." );
            continue;
        }
 
        const Trk::Surface& end_surf = end_parameters->associatedSurface();
        if (!start_parameters) {
            start_parameters = end_parameters;
        } else {         
            extrapolateBetweenHits(start_parameters, end_surf, holes, partHyp);
            start_parameters = end_parameters;
        }
    } // end loop over TrackStateOnSurfaces
 
    if( !m_end_at_last_trt_hit ) {
        // final extrapolation to the edge of TRT to check for holes in outermost straws
        int trailing_hole_count = extrapolateBetweenHits(start_parameters, *m_trt_outer_surf, holes, partHyp);
 
        // remove trailing holes
        if(trailing_hole_count > m_max_trailing_holes) {
            ATH_MSG_DEBUG("There are " << trailing_hole_count << " trailing holes removed.");
            for(int i=0; i < trailing_hole_count; i++) {
                holes->pop_back();
                /*
                  Note: DataVector::pop_back() deletes pointers if it owns them.
                */
            }
        }
    }
 
    return holes;
}

getTrackWithHoles()

const Trk::Track * TRTTrackHoleSearchTool::getTrackWithHoles ( const Trk::Track & track, const Trk::ParticleHypothesis partHyp = Trk::pion ) const

virtual

Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes The parthyp argument is relevant for the extrapolation steps in the hole search.

Attention: This is a factory, ownership of the return track is passed to the calling method.

Implements Trk::ITrackHoleSearchTool.

Definition at line 204 of file TRTTrackHoleSearchTool.cxx.

                                                                                                       {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::getTrackWithHoles" );
    const Trk::TrackStates* holes = getHolesOnTrack(track, partHyp);
    const Trk::Track* new_track = addHolesToTrack(track, holes);
    delete holes;
    return new_track;
}

getTrackWithHolesAndOutliers()

const Trk::Track * TRTTrackHoleSearchTool::getTrackWithHolesAndOutliers ( const Trk::Track & track, const Trk::ParticleHypothesis partHyp = Trk::pion ) const

virtual

Input : track, parthyp Return: A pointer to a new Trk::Track which containes the information of the input track plus the tsos of the identified holes or outliers The parthyp argument is relevant for the extrapolation steps in the hole search.

Attention: This is a factory, ownership of the return track is passed to the calling method.

Implements Trk::ITrackHoleSearchTool.

Definition at line 215 of file TRTTrackHoleSearchTool.cxx.

                                                                                                                  {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::getTrackWithHolesAndOutliers" );
    return getTrackWithHoles(track, partHyp);
}

initialize()

StatusCode TRTTrackHoleSearchTool::initialize

(

)

Definition at line 47 of file TRTTrackHoleSearchTool.cxx.

                                              {
    ATH_MSG_DEBUG( "TRTTrackHoleSearchTool::initialize()" );
 
    // retrieve extrapolator
    ATH_CHECK(m_extrapolator.retrieve());
 
    // retrieve TRT_ID
    ATH_CHECK(detStore()->retrieve(m_TRT_ID, "TRT_ID"));
 
    // retrieve ConditionsSummarySvc
    ATH_CHECK(m_conditions_svc.retrieve());
 
    m_trt_outer_surf = new Trk::CylinderSurface(Amg::Transform3D(Amg::Transform3D::Identity()),
     m_outer_radius, m_max_z);
    // note: HepGeom::Translate3D is deleted by Trk::Surface destructor
 
    return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_bcheck

BooleanProperty TRTTrackHoleSearchTool::m_bcheck {this, "bcheck", false}

private

Definition at line 97 of file TRTTrackHoleSearchTool.h.

{this, "bcheck", false};

m_begin_at_first_trt_hit

BooleanProperty TRTTrackHoleSearchTool::m_begin_at_first_trt_hit {this, "begin_at_first_trt_hit", false}

private

Definition at line 95 of file TRTTrackHoleSearchTool.h.

{this, "begin_at_first_trt_hit", false}; // if not, extrapolate from last Si hit

m_conditions_svc

ServiceHandle<IInDetConditionsSvc> TRTTrackHoleSearchTool::m_conditions_svc {this, "conditions_svc", "TRT_ConditionsSummarySvc"}

private

Definition at line 87 of file TRTTrackHoleSearchTool.h.

{this, "conditions_svc", "TRT_ConditionsSummarySvc"};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_do_dump_bad_straw_log

BooleanProperty TRTTrackHoleSearchTool::m_do_dump_bad_straw_log {this, "do_dump_bad_straw_log", false}

private

Definition at line 98 of file TRTTrackHoleSearchTool.h.

{this, "do_dump_bad_straw_log", false};

m_end_at_last_trt_hit

BooleanProperty TRTTrackHoleSearchTool::m_end_at_last_trt_hit {this, "end_at_last_trt_hit", false}

private

Definition at line 96 of file TRTTrackHoleSearchTool.h.

{this, "end_at_last_trt_hit", false}; // if not, continue hole search to the edge of the TRT

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extrapolator

ToolHandle<Trk::IExtrapolator> TRTTrackHoleSearchTool::m_extrapolator {this, "extrapolator", "Trk::Extrapolator"}

private

Definition at line 85 of file TRTTrackHoleSearchTool.h.

{this, "extrapolator", "Trk::Extrapolator"};

m_has_been_called

std::atomic_bool TRTTrackHoleSearchTool::m_has_been_called {false}

mutableprivate

Definition at line 105 of file TRTTrackHoleSearchTool.h.

{false};

m_locR_cut

FloatProperty TRTTrackHoleSearchTool::m_locR_cut {this, "locR_cut", -1.}

private

Definition at line 99 of file TRTTrackHoleSearchTool.h.

{this, "locR_cut", -1.}; // 1.4*CLHEP::mm // negative means no cut

m_locR_sigma_cut

FloatProperty TRTTrackHoleSearchTool::m_locR_sigma_cut {this, "locR_sigma_cut", -1.}

private

Definition at line 100 of file TRTTrackHoleSearchTool.h.

{this, "locR_sigma_cut", -1.};

m_locZ_cut

FloatProperty TRTTrackHoleSearchTool::m_locZ_cut {this, "locZ_cut", 5.0*CLHEP::mm}

private

Definition at line 101 of file TRTTrackHoleSearchTool.h.

{this, "locZ_cut", 5.0*CLHEP::mm};

m_max_trailing_holes

IntegerProperty TRTTrackHoleSearchTool::m_max_trailing_holes {this, "max_trailing_holes", 1}

private

Definition at line 94 of file TRTTrackHoleSearchTool.h.

{this, "max_trailing_holes", 1}; // only used if not end_at_last_trt_hit

m_max_z

FloatProperty TRTTrackHoleSearchTool::m_max_z {this, "max_z", 2715.0*CLHEP::mm}

private

Definition at line 93 of file TRTTrackHoleSearchTool.h.

{this, "max_z", 2715.0*CLHEP::mm};

m_outer_radius

FloatProperty TRTTrackHoleSearchTool::m_outer_radius {this, "outer_radius", 1075.0*CLHEP::mm}

private

Definition at line 91 of file TRTTrackHoleSearchTool.h.

{this, "outer_radius", 1075.0*CLHEP::mm};

m_TRT_ID

const TRT_ID* TRTTrackHoleSearchTool::m_TRT_ID {nullptr}

private

Definition at line 106 of file TRTTrackHoleSearchTool.h.

{nullptr};

m_trt_outer_surf

Trk::CylinderSurface* TRTTrackHoleSearchTool::m_trt_outer_surf {nullptr}

private

Definition at line 107 of file TRTTrackHoleSearchTool.h.

{nullptr};

m_use_conditions_svc

BooleanProperty TRTTrackHoleSearchTool::m_use_conditions_svc {this, "use_conditions_svc", true}

private

Definition at line 89 of file TRTTrackHoleSearchTool.h.

{this, "use_conditions_svc", true};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• TRTTrackHoleSearchTool.h
• TRTTrackHoleSearchTool.cxx