ActsTrk::TracccTrackConverterAlg Class Reference

Algorithm converting traccc tracks ( tracks + track states sdevice buffer) to ACTS tracks (host container). More...

#include <TracccTrackConverterAlg.h>

Inheritance diagram for ActsTrk::TracccTrackConverterAlg:

Public Member Functions
virtual StatusCode	initialize () override
	Function initializing the algorithm.
virtual StatusCode	execute (const EventContext &ctx) const override
	Function executing the algorithm.
virtual StatusCode	finalize () override
	Function finalizing the algorthm.
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
	Get filter decision:
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
	Set filter decision:
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	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
const xAOD::UncalibratedMeasurement &	makeSourceLink (const unsigned int &meas_index, std::span< const unsigned int > pixelMap, std::span< const unsigned int > stripMap, const xAOD::PixelClusterContainer &pixel_clusters, const xAOD::StripClusterContainer &strip_clusters) const
const Acts::Surface &	findActsSurface (const Acts::GeometryIdentifier &actsID) const
std::optional< Acts::BoundTrackParameters >	convertGlobalToActsParameters (traccc::bound_track_parameters< traccc::default_algebra > const &trkParams) const
template<typename state_t>
std::optional< Acts::BoundTrackParameters >	convertSmoothedToActsParameters (traccc::edm::track_state< state_t > const &state) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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
The input device resident cluster, measurement and cell collection names
{@
SG::ReadHandleKey< xAOD::PixelClusterContainer >	m_inputPixelClustersKey
SG::ReadHandleKey< xAOD::StripClusterContainer >	m_inputStripClustersKey
SG::ReadHandleKey< std::vector< unsigned int > >	m_inputMeasToPixelSPKey
SG::ReadHandleKey< std::vector< unsigned int > >	m_inputMeasToStripClKey
SG::ReadHandleKey< traccc_track_container::buffer >	m_inputTracksKey
The output host resident track container name
{@
SG::WriteHandleKey< ActsTrk::TrackContainer >	m_outputTracksKey
The host memory resource tool to use for memory allocations
ToolHandle< AthDevice::IMemoryResourceTool >	m_hostMR
The copy tool used for copying data from device
ToolHandle< AthDevice::ICopyTool >	m_copy
std::shared_ptr< const Acts::TrackingGeometry >	m_trackingGeometry
ServiceHandle< ActsTrk::ITrackingGeometrySvc >	m_trackingGeometrySvc {this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"}
std::map< Acts::GeometryIdentifier, const Acts::Surface * >	m_actsSurfaceMap
ActsTrk::MutableTrackContainerHandlesHelper	m_tracksBackendHandlesHelper
The detector description service providing the Athena<->Detray ID map
Gaudi::Property< std::string >	m_geoIdMappingObjectName {this, "GeoIdMapping", "", "ID mapping between the three detector description realms."}
const ActsTrk::GeometryIdMapping *	m_idMapping {nullptr}
Gaudi::Property< std::string >	m_hostDetectorObjectName
const traccc::host_detector *	m_hostDetector = nullptr
std::atomic< int >	m_nTracksIn = 0
	The object counters for debug prints in finalize method {.
std::atomic< int >	m_nTracksOut = 0
std::atomic< int >	m_nExcludedFitOutcome = 0
std::atomic< int >	m_nExcludedNoState = 0
std::atomic< int >	m_nExcludedBadNdf = 0
std::atomic< int >	m_nExcludedWeirdState = 0
std::atomic< int >	m_nExcludedWeirdGlobal = 0

Detailed Description

Algorithm converting traccc tracks ( tracks + track states sdevice buffer) to ACTS tracks (host container).

This algorithm retrieves the input device resident traccc track and track states collection from the event store, copies the data to host buffer and converts the traccc tracks to ActsTrk::Track. It works under the assumption that the other relevant traccc objects have already been converted to xAOD and are available on host in the event store and that the indices between the HOST and DEVICE resident containers have been recorded in a index map.

Author

Neža Ribarič neza..nosp@m.riba.nosp@m.ric@c.nosp@m.ern..nosp@m.ch

Definition at line 45 of file TracccTrackConverterAlg.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Function Documentation

convertGlobalToActsParameters()

std::optional< Acts::BoundTrackParameters > ActsTrk::TracccTrackConverterAlg::convertGlobalToActsParameters ( traccc::bound_track_parameters< traccc::default_algebra > const & trkParams ) const

private

Definition at line 126 of file TracccTrackConverterAlg.cxx.

{
    if (trkParams.bound_local()[0] == 0 && trkParams.bound_local()[1] == 0 &&
        trkParams.phi() == 0 && trkParams.theta() == 0 &&
        trkParams.qop() == 0 && trkParams.time() == 0) {
        // traccc reports degenerate states like this; treat as "no
        // parameters" rather than as a real (but zeroed-out) track.
        return std::nullopt;
    }
 
    const auto& detrayDetector = m_hostDetector->as<traccc::itk_detector>();
    const detray::tracking_surface detray_surface{detrayDetector, trkParams.surface_link()};
    const auto geo_id = detray_surface.source();
    const Acts::GeometryIdentifier acts_geom_id{geo_id};
 
    const Acts::Surface& surface = findActsSurface(acts_geom_id);
 
    Acts::BoundVector params;
    params << trkParams.bound_local()[0], trkParams.bound_local()[1],
        trkParams.phi(), trkParams.theta(), trkParams.qop(), trkParams.time();
 
    // ---- sanity-check units ----
    const double qop = trkParams.qop();
    const double p   = (qop != 0) ? 1.0 / std::abs(qop) : 0.0;
    const double pT  = p * std::sin(trkParams.theta());
    ATH_MSG_DEBUG("Global params: d0=" << trkParams.bound_local()[0]
                  << " z0=" << trkParams.bound_local()[1]
                  << " phi=" << trkParams.phi()
                  << " theta=" << trkParams.theta()
                  << " qop=" << qop
                  << " -> |p|=" << p << " pT=" << pT
                  << " charge=" << (qop > 0 ? +1 : -1)
                  << " time=" << trkParams.time());    
 
    // Traccc covariance has no time uncertainty; Acts::BoundMatrix does
    // (it's eBoundSize x eBoundSize == 6x6). For now leave the time row/column at the
    // Identity() values set above and only fill the spatial+qop values.
    static constexpr unsigned int kAtlasCovSize = 5; // d0, z0, phi, theta, qop (no time)
 
    Acts::BoundMatrix cov = Acts::BoundMatrix::Identity();
    const auto& atlasCov = trkParams.covariance();
    for (unsigned i = 0; i < kAtlasCovSize; ++i) {
        for (unsigned j = 0; j < kAtlasCovSize; ++j) {
            cov(i, j) = atlasCov[i][j];
        }
    }
 
    return Acts::BoundTrackParameters(surface.getSharedPtr(), params, cov,
                                      Acts::ParticleHypothesis::pion());
}

convertSmoothedToActsParameters()

template<typename state_t>

std::optional< Acts::BoundTrackParameters > ActsTrk::TracccTrackConverterAlg::convertSmoothedToActsParameters ( traccc::edm::track_state< state_t > const & state ) const

private

Definition at line 181 of file TracccTrackConverterAlg.cxx.

{
    const auto& atlasParam = state.smoothed_params();
 
    if (atlasParam.bound_local()[0] == 0 && atlasParam.bound_local()[1] == 0 &&
        atlasParam.phi() == 0 && atlasParam.theta() == 0 &&
        atlasParam.qop() == 0 && atlasParam.time() == 0) {
        // traccc reports degenerate states like this; treat as "no
        // parameters" rather than as a real (but zeroed-out) state.
        return std::nullopt;
    }
 
    const auto& detrayDetector = m_hostDetector->as<traccc::itk_detector>();
    const detray::tracking_surface detray_surface{detrayDetector, atlasParam.surface_link()};
    const auto geo_id = detray_surface.source();
    const Acts::GeometryIdentifier acts_geom_id{geo_id};
 
    const Acts::Surface& surface = findActsSurface(acts_geom_id);
 
    // d0/z0 come from the measurement's local position; the remaining bound
    // parameters come from the smoothed track state.
    Acts::BoundVector params;
    params << atlasParam.bound_local()[0], atlasParam.bound_local()[1],
        atlasParam.phi(), atlasParam.theta(), atlasParam.qop(),
        atlasParam.time();
 
    // ---- per-state momentum check for sanity ----
    const double qop = atlasParam.qop();
    const double p   = (qop != 0) ? 1.0 / std::abs(qop) : 0.0;
    const double pT  = p * std::sin(atlasParam.theta());
    ATH_MSG_DEBUG("Smoothed state: d0=" << atlasParam.bound_local()[0]
                  << " z0=" << atlasParam.bound_local()[1]
                  << " phi=" << atlasParam.phi()
                  << " theta=" << atlasParam.theta()
                  << " qop=" << qop
                  << " -> |p|=" << p << " pT=" << pT);    
 
    // Traccc covariance has no time uncertainty; Acts::BoundMatrix does
    // (it's eBoundSize x eBoundSize == 6x6). For now leave the time row/column at the
    // Identity() values set above and only fill the spatial+qop values.
    static constexpr unsigned int kAtlasCovSize = 5; // d0, z0, phi, theta, qop (no time)
 
    Acts::BoundMatrix cov = Acts::BoundMatrix::Identity();
    const auto& atlasCov = atlasParam.covariance();
    for (unsigned i = 0; i < kAtlasCovSize; ++i) {
        for (unsigned j = 0; j < kAtlasCovSize; ++j) {
            cov(i, j) = atlasCov[i][j];
        }
    }
 
    return Acts::BoundTrackParameters(surface.getSharedPtr(), params, cov,
                                      Acts::ParticleHypothesis::pion());
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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; }

execute()

StatusCode ActsTrk::TracccTrackConverterAlg::execute ( const EventContext & ctx ) const

overridevirtual

Function executing the algorithm.

Definition at line 236 of file TracccTrackConverterAlg.cxx.

{
    // ---- Retrieve HOST RESIDENT clusters (always needed) ----
    // These are made during the TracccMeasurementConverterAlg 
    auto pixel_clusters = SG::makeHandle(m_inputPixelClustersKey, ctx);
    ATH_CHECK(pixel_clusters.isValid());
 
    auto strip_clusters = SG::makeHandle(m_inputStripClustersKey, ctx);
    ATH_CHECK(strip_clusters.isValid());
 
    // ---- Retrieve mapping from traccc measurement index to pixel cluster/spacepoint index (always needed) ----
    // Because the SPs were created from xAOD pixel clusters on the host, the same traccc measurement index 
    // points to the xAOD pixel cluster and xAOD spacepoint in their respective containers
    auto pixelMeasMap = SG::makeHandle(m_inputMeasToPixelSPKey, ctx);
    ATH_CHECK(pixelMeasMap.isValid());
 
    auto stripMeasMap = SG::makeHandle(m_inputMeasToStripClKey, ctx);
    ATH_CHECK(stripMeasMap.isValid());
 
    if(stripMeasMap->size() == 0 || pixelMeasMap->size() == 0){
        ATH_MSG_FATAL("Maps relating traccc measurements to xAOD cluster containers are empty!!");
        return StatusCode::FAILURE;
    }
 
    // ---- Retrieve DEVICE resident traccc tracks and track states ----
    auto tracks = SG::makeHandle(m_inputTracksKey, ctx);
    ATH_CHECK(tracks.isValid());
 
    auto copy = m_copy->copy(ctx);
 
    traccc_track_container::buffer traccc_tracks_buffer;
    traccc_tracks_buffer.tracks =
        copy->to(tracks->tracks, m_hostMR->mr(),
            nullptr, vecmem::copy::type::device_to_host);
    traccc_tracks_buffer.states =
        copy->to(tracks->states, m_hostMR->mr(),
            nullptr, vecmem::copy::type::device_to_host);
        
    traccc_track_container::const_device traccc_tracks(
        traccc_tracks_buffer);  
 
    ATH_MSG_DEBUG("Read " << traccc_tracks.tracks.size() << " tracks from device.");
    m_nTracksIn += traccc_tracks.tracks.size();
 
    // -- Write HOST resident ACTS track container ----
    Acts::VectorTrackContainer trackBackend;
    Acts::VectorMultiTrajectory trackStateBackend;
    ActsTrk::MutableTrackContainer trackContainer(std::move(trackBackend),
                                                  std::move(trackStateBackend));
    
    // Bookkeeping for debug summary
    unsigned nExcludedFitOutcome = 0;
    unsigned nExcludedNoState = 0;
    unsigned nExcludedBadNdf = 0;
    unsigned nExcludedWeirdState = 0;
    unsigned nExcludedWeirdGlobal = 0;
 
    for (std::size_t i = 0; i < traccc_tracks.tracks.size(); ++i) {
 
        const traccc::edm::track track = traccc_tracks.tracks.at(i);
 
        const auto fitOutcome = track.fit_outcome();
        if (fitOutcome == traccc::track_fit_outcome::FAILURE_NON_POSITIVE_NDF ||
            fitOutcome == traccc::track_fit_outcome::FAILURE_NOT_ALL_SMOOTHED ||
            fitOutcome == traccc::track_fit_outcome::UNKNOWN) {
            ATH_MSG_DEBUG("Skipping track " << i << ": fit outcome "
                                            << static_cast<int>(fitOutcome));
            ++nExcludedFitOutcome;
            continue;
        }
 
        if (track.constituent_links().empty()) {
            ++nExcludedNoState;
            continue;
        }
 
        if (track.ndf() >
                static_cast<float>(std::numeric_limits<unsigned int>::max()) ||
            track.ndf() <
                static_cast<float>(std::numeric_limits<unsigned int>::min())) {
            ++nExcludedBadNdf;
            continue;
        }
 
        TrackValidity track_validity = kValid;
 
        auto actsTrack = trackContainer.makeTrack();
        actsTrack.chi2() = track.chi2();
        actsTrack.nDoF() = track.ndf(); 
 
        Acts::TrackStatePropMask const state_mask =
            Acts::TrackStatePropMask::Smoothed;
 
        bool firstState = true;
 
        for (const auto [linkType, stateIdx] : track.constituent_links()) {
            
            assert(linkType == traccc::edm::track_constituent_link::track_state);
 
            const auto& state = traccc_tracks.states.at(stateIdx);
            const auto& meas_index = state.measurement_index();
 
            auto trackState = actsTrack.appendTrackState(state_mask);
            trackState.typeFlags().setIsMeasurement();
 
            const std::optional<Acts::BoundTrackParameters> smoothed =
                convertSmoothedToActsParameters(state);
 
            if (!smoothed) {
                ATH_MSG_DEBUG("Track " << i << ": degenerate smoothed state, "
                                             "dropping track");
                track_validity = TrackValidity::kInvalidState;
                break;
            }
 
            const xAOD::UncalibratedMeasurement* sourceLink =
                &makeSourceLink(meas_index, *pixelMeasMap, *stripMeasMap, *pixel_clusters, *strip_clusters);
            
            trackState.setUncalibratedSourceLink(ActsTrk::detail::xAODUncalibMeasCalibrator::pack(sourceLink));
 
            // traccc does not yet do backpropagation, so there is no true
            // reference surface (perigee) for the track. We use the surface
            // of the first measurement instead and rely on back-propagation
            // during the later Acts -> xAOD conversion step to find the real
            // perigee.
            if (firstState) {
                const std::optional<Acts::BoundTrackParameters> global =
                    convertGlobalToActsParameters(track.params());
                if (!global) {
                    ATH_MSG_DEBUG("Track "
                                  << i
                                  << ": degenerate global parameters, "
                                     "dropping track");
                    track_validity = TrackValidity::kInvalidGlobalParams;
                    break;
                }
                actsTrack.parameters() = global->parameters();
                actsTrack.covariance() = *global->covariance();
                actsTrack.setReferenceSurface(
                    global->referenceSurface().getSharedPtr());
                firstState = false;
            }
 
            try {
                trackState.setReferenceSurface(
                    smoothed->referenceSurface().getSharedPtr());
                trackState.smoothed() = smoothed->parameters();
                trackState.smoothedCovariance() = *smoothed->covariance();
            } catch (const std::exception& e) {
                ATH_MSG_ERROR("Track " << i << ": failed to set track state ("
                                       << e.what() << ")");
            }
 
        }
 
        if (track_validity == TrackValidity::kInvalidState) {
            ++nExcludedWeirdState;
            trackContainer.removeTrack(actsTrack.index());
        } else if (track_validity == TrackValidity::kInvalidGlobalParams) {
            ++nExcludedWeirdGlobal;
            trackContainer.removeTrack(actsTrack.index());
        } else {
            // sanity checks
            const auto& pars = actsTrack.parameters();
            const double qop = pars[Acts::eBoundQOverP];
            const double p   = (qop != 0) ? 1.0 / std::abs(qop) : 0.0;
            const double pT  = p * std::sin(pars[Acts::eBoundTheta]);
            ATH_MSG_DEBUG("Track " << i << " ACCEPTED: nHits="
                          << track.constituent_links().size()
                          << " chi2=" << track.chi2()
                          << " ndf=" << track.ndf()
                          << " |p|=" << p << " GeV pT=" << pT << " GeV");
        }
 
    }
 
    ATH_MSG_DEBUG("Converted " << trackContainer.size() << " tracks to " << m_outputTracksKey.key() );
    ATH_MSG_DEBUG("excluded: "
                            << nExcludedFitOutcome << " (fit outcome), "
                            << nExcludedNoState << " (no state), "
                            << nExcludedBadNdf << " (bad ndf), "
                            << nExcludedWeirdState << " (weird state), "
                            << nExcludedWeirdGlobal << " (weird global params)");
    
    m_nExcludedFitOutcome += nExcludedFitOutcome;
    m_nExcludedNoState += nExcludedNoState;
    m_nExcludedBadNdf += nExcludedBadNdf;
    m_nExcludedWeirdState += nExcludedWeirdState;
    m_nExcludedWeirdGlobal += nExcludedWeirdGlobal;
    
    m_nTracksOut += trackContainer.size();
    
    Acts::ConstVectorTrackContainer constTrackBackend(
        std::move(trackContainer.container()));
    Acts::ConstVectorMultiTrajectory constTrackStateBackend(
        std::move(trackContainer.trackStateContainer()));
    auto constTrackContainer = std::make_unique<ActsTrk::TrackContainer>(
        std::move(constTrackBackend), std::move(constTrackStateBackend));
 
    SG::WriteHandle<ActsTrk::TrackContainer> outputHandle(
        m_outputTracksKey, ctx);
    ATH_CHECK(outputHandle.record(std::move(constTrackContainer)));                        
 
    return StatusCode::SUCCESS;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

extraOutputDeps()

const DataObjIDColl & AthCommonAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 89 of file AthCommonAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

filterPassed()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

finalize()

StatusCode ActsTrk::TracccTrackConverterAlg::finalize ( )

overridevirtual

Function finalizing the algorthm.

Definition at line 442 of file TracccTrackConverterAlg.cxx.

{
    ATH_MSG_DEBUG("Finalizing.");
 
    ATH_MSG_DEBUG("Received " << m_nTracksIn << " tracks, wrote " << m_nTracksOut);
    ATH_MSG_DEBUG("Excluded: "
                            << m_nExcludedFitOutcome << " (fit outcome), "
                            << m_nExcludedNoState << " (no state), "
                            << m_nExcludedBadNdf << " (bad ndf), "
                            << m_nExcludedWeirdState << " (weird state), "
                            << m_nExcludedWeirdGlobal << " (weird global params)");
 
    ATH_MSG_DEBUG("Successfully finalized");
    return StatusCode::SUCCESS;
}

findActsSurface()

const Acts::Surface & ActsTrk::TracccTrackConverterAlg::findActsSurface ( const Acts::GeometryIdentifier & actsID ) const

private

Definition at line 113 of file TracccTrackConverterAlg.cxx.

{
    const auto surfaceIt = m_actsSurfaceMap.find(actsID);
    if (surfaceIt == m_actsSurfaceMap.end()) {
        ATH_MSG_ERROR("No Acts surface known for Acts geometry id "
                      << actsID);
        throw std::domain_error("No Acts surface for this geometry id");
    }
    return *surfaceIt->second;
}

initialize()

StatusCode ActsTrk::TracccTrackConverterAlg::initialize ( )

overridevirtual

Function initializing the algorithm.

Definition at line 19 of file TracccTrackConverterAlg.cxx.

{
    ATH_MSG_DEBUG("Initializing.");
 
    ATH_CHECK(m_hostMR.retrieve());
    ATH_CHECK(m_copy.retrieve());
 
    ATH_CHECK(m_inputPixelClustersKey.initialize());
    ATH_CHECK(m_inputStripClustersKey.initialize());
    ATH_CHECK(m_inputMeasToPixelSPKey.initialize());
    ATH_CHECK(m_inputMeasToStripClKey.initialize());
 
    ATH_CHECK(m_inputTracksKey.initialize());
    ATH_CHECK(m_outputTracksKey.initialize());
 
    ATH_CHECK(detStore()->retrieve(m_idMapping, m_geoIdMappingObjectName.value()));
    ATH_CHECK(detStore()->retrieve(m_hostDetector, m_hostDetectorObjectName.value()));
 
    // Build the ACTS-surface <-> ACTS-id lookup map once,
    // up front, rather than re-deriving them per event.
    // Possibly this will be built during ACTS->Detray conversion
    // Possibly this is an issue when the tracking geometry changes between events
    if (!m_trackingGeometrySvc.empty()) {
        ATH_CHECK(m_trackingGeometrySvc.retrieve());
        m_trackingGeometry = m_trackingGeometrySvc->trackingGeometry();
        m_trackingGeometry->visitSurfaces([&](const Acts::Surface* surface) {
            if (!surface) return;
            // m_allsurfaces.push_back(surface->getSharedPtr());
            const auto *actsElement = getActsDetectorElement(surface);
            if (!actsElement){
                ATH_MSG_DEBUG("Could not find matching Acts detector element for surface with geometryId " << surface->geometryId());
                return;
            } 
            if(dynamic_cast<const InDetDD::HGTD_DetectorElement*>(actsElement->upstreamDetectorElement()) != nullptr) { 
                // skip HGTD surfaces for now
                // currently no time info is being used in traccc
                ATH_MSG_VERBOSE("Found HGTD surface with geometryId " << surface->geometryId());
                return;
            }
            const auto *geoElement = actsElement->upstreamDetectorElement();
            const auto *detElem = dynamic_cast<const InDetDD::SiDetectorElement*>(geoElement);
            if (!geoElement || !detElem) {
                ATH_MSG_DEBUG("Could not find matching Athena silicon detector element for surface with geometryId " << surface->geometryId());
                return;
            }
    
            const Acts::GeometryIdentifier actsID = surface->geometryId();
            const auto [it, inserted] = m_actsSurfaceMap.insert({actsID, surface});
            if (!inserted) {
                ATH_MSG_WARNING("ACTS id " << actsID
                                        << " maps to two Acts surfaces: "
                                        << it->second->geometryId() << " and "
                                        << surface->geometryId());
            }
    
        });
 
    }else {
        ATH_MSG_FATAL("TrackingGeometrySvc is not configured, cannot translate Traccc tracks to ACTS tracks");
        return StatusCode::FAILURE;
    }    
 
    ATH_CHECK(m_tracksBackendHandlesHelper.initialize(
        ActsTrk::prefixFromTrackContainerName(
            m_outputTracksKey.key())));
 
    ATH_MSG_DEBUG("Successfully initialized");
    return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

isClonable()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

                                           {
    return true;
  }

makeSourceLink()

const xAOD::UncalibratedMeasurement & ActsTrk::TracccTrackConverterAlg::makeSourceLink ( const unsigned int & meas_index, std::span< const unsigned int > pixelMap, std::span< const unsigned int > stripMap, const xAOD::PixelClusterContainer & pixel_clusters, const xAOD::StripClusterContainer & strip_clusters ) const

private

Definition at line 90 of file TracccTrackConverterAlg.cxx.

{
    const unsigned int pixel_index = (pixelMap)[meas_index];
    if(pixel_index != std::numeric_limits<unsigned int>::max()){
    
        return *pixel_clusters.at(pixel_index);
    }
 
    const unsigned int strip_index = (stripMap)[meas_index];
    if(strip_index != std::numeric_limits<unsigned int>::max()){
        return *strip_clusters.at(strip_index);
    }
 
    ATH_MSG_FATAL("Measurement index could not be found!!!");
    throw std::domain_error("No measurement index match found in xAOD map");
    
}

msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

virtual void AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, HypoBase, InputMakerBase, and PyAthena::Alg.

Definition at line 60 of file AthCommonAlgorithm.cxx.

                                                      {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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_actsSurfaceMap

std::map<Acts::GeometryIdentifier, const Acts::Surface*> ActsTrk::TracccTrackConverterAlg::m_actsSurfaceMap

private

Definition at line 110 of file TracccTrackConverterAlg.h.

m_copy

ToolHandle<AthDevice::ICopyTool> ActsTrk::TracccTrackConverterAlg::m_copy

private

Initial value:

{
      this, "CopyProviderTool", "", "Vecmem copy provider tool"}

Definition at line 105 of file TracccTrackConverterAlg.h.

                                       {
      this, "CopyProviderTool", "", "Vecmem copy provider tool"};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

m_geoIdMappingObjectName

Gaudi::Property<std::string> ActsTrk::TracccTrackConverterAlg::m_geoIdMappingObjectName {this, "GeoIdMapping", "", "ID mapping between the three detector description realms."}

private

Definition at line 115 of file TracccTrackConverterAlg.h.

{this, "GeoIdMapping", "", "ID mapping between the three detector description realms."};

m_hostDetector

const traccc::host_detector* ActsTrk::TracccTrackConverterAlg::m_hostDetector = nullptr

private

Definition at line 120 of file TracccTrackConverterAlg.h.

m_hostDetectorObjectName

Gaudi::Property<std::string> ActsTrk::TracccTrackConverterAlg::m_hostDetectorObjectName

private

Initial value:

{
        this, "HostDetectorName", "",
        "Detray host detector object"}

Definition at line 117 of file TracccTrackConverterAlg.h.

                                                   {
        this, "HostDetectorName", "",
        "Detray host detector object"};

m_hostMR

ToolHandle<AthDevice::IMemoryResourceTool> ActsTrk::TracccTrackConverterAlg::m_hostMR

private

Initial value:

{
    this, "HostMR", "", "Host memory resource tool"}

Definition at line 102 of file TracccTrackConverterAlg.h.

                                                   {
    this, "HostMR", "", "Host memory resource tool"};

m_idMapping

const ActsTrk::GeometryIdMapping* ActsTrk::TracccTrackConverterAlg::m_idMapping {nullptr}

private

Definition at line 116 of file TracccTrackConverterAlg.h.

{nullptr};

m_inputMeasToPixelSPKey

SG::ReadHandleKey<std::vector<unsigned int> > ActsTrk::TracccTrackConverterAlg::m_inputMeasToPixelSPKey

private

Initial value:

{
      this, "InputMeasToPixelSP", "",
      "Input mapping from traccc measurement index to pixel spacepoint index"}

Definition at line 83 of file TracccTrackConverterAlg.h.

                                                                  {
      this, "InputMeasToPixelSP", "",
      "Input mapping from traccc measurement index to pixel spacepoint index"};    

m_inputMeasToStripClKey

SG::ReadHandleKey<std::vector<unsigned int> > ActsTrk::TracccTrackConverterAlg::m_inputMeasToStripClKey

private

Initial value:

{
      this, "InputMeasToStripCl", "",
      "Input mapping from traccc measurement index to strip cluster index"}

Definition at line 86 of file TracccTrackConverterAlg.h.

                                                                  {
      this, "InputMeasToStripCl", "",
      "Input mapping from traccc measurement index to strip cluster index"};            

m_inputPixelClustersKey

SG::ReadHandleKey<xAOD::PixelClusterContainer> ActsTrk::TracccTrackConverterAlg::m_inputPixelClustersKey

private

Initial value:

{
      this, "InputPixelClusters", "",
      "Input xAOD pixel clusters converted from Traccc"}

Definition at line 77 of file TracccTrackConverterAlg.h.

                                                                    {
      this, "InputPixelClusters", "",
      "Input xAOD pixel clusters converted from Traccc"};

m_inputStripClustersKey

SG::ReadHandleKey<xAOD::StripClusterContainer> ActsTrk::TracccTrackConverterAlg::m_inputStripClustersKey

private

Initial value:

{
      this, "InputStripClusters", "",
      "Input xAOD strip clusters converted from Traccc"}

Definition at line 80 of file TracccTrackConverterAlg.h.

                                                                    {
      this, "InputStripClusters", "",
      "Input xAOD strip clusters converted from Traccc"};

m_inputTracksKey

SG::ReadHandleKey<traccc_track_container::buffer> ActsTrk::TracccTrackConverterAlg::m_inputTracksKey

private

Initial value:

{
      this, "InputTracks", "",
      "Input traccc track collection buffer"}

Definition at line 89 of file TracccTrackConverterAlg.h.

                                                                {
      this, "InputTracks", "",
      "Input traccc track collection buffer"}; 

m_nExcludedBadNdf

std::atomic<int> ActsTrk::TracccTrackConverterAlg::m_nExcludedBadNdf = 0

mutableprivate

Definition at line 127 of file TracccTrackConverterAlg.h.

m_nExcludedFitOutcome

std::atomic<int> ActsTrk::TracccTrackConverterAlg::m_nExcludedFitOutcome = 0

mutableprivate

Definition at line 125 of file TracccTrackConverterAlg.h.

m_nExcludedNoState

std::atomic<int> ActsTrk::TracccTrackConverterAlg::m_nExcludedNoState = 0

mutableprivate

Definition at line 126 of file TracccTrackConverterAlg.h.

m_nExcludedWeirdGlobal

std::atomic<int> ActsTrk::TracccTrackConverterAlg::m_nExcludedWeirdGlobal = 0

mutableprivate

Definition at line 129 of file TracccTrackConverterAlg.h.

m_nExcludedWeirdState

std::atomic<int> ActsTrk::TracccTrackConverterAlg::m_nExcludedWeirdState = 0

mutableprivate

Definition at line 128 of file TracccTrackConverterAlg.h.

m_nTracksIn

std::atomic<int> ActsTrk::TracccTrackConverterAlg::m_nTracksIn = 0

mutableprivate

The object counters for debug prints in finalize method {.

Definition at line 123 of file TracccTrackConverterAlg.h.

m_nTracksOut

std::atomic<int> ActsTrk::TracccTrackConverterAlg::m_nTracksOut = 0

mutableprivate

Definition at line 124 of file TracccTrackConverterAlg.h.

m_outputTracksKey

SG::WriteHandleKey<ActsTrk::TrackContainer> ActsTrk::TracccTrackConverterAlg::m_outputTracksKey

private

Initial value:

{
      this, "OutputTracks", "",
      "Output ACTS track container"}

Definition at line 96 of file TracccTrackConverterAlg.h.

                                                           {
      this, "OutputTracks", "",
      "Output ACTS track container"};

m_trackingGeometry

std::shared_ptr<const Acts::TrackingGeometry> ActsTrk::TracccTrackConverterAlg::m_trackingGeometry

private

Definition at line 108 of file TracccTrackConverterAlg.h.

m_trackingGeometrySvc

ServiceHandle<ActsTrk::ITrackingGeometrySvc> ActsTrk::TracccTrackConverterAlg::m_trackingGeometrySvc {this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"}

private

Definition at line 109 of file TracccTrackConverterAlg.h.

{this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"};    

m_tracksBackendHandlesHelper

ActsTrk::MutableTrackContainerHandlesHelper ActsTrk::TracccTrackConverterAlg::m_tracksBackendHandlesHelper

private

Initial value:

{
        this}

Definition at line 112 of file TracccTrackConverterAlg.h.

                                                                        {
        this};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• TracccTrackConverterAlg.h
• TracccTrackConverterAlg.cxx