ActsTrk::TrackToTrackParticleCnvAlg Class Reference

#include <TrackToTrackParticleCnvAlg.h>

Inheritance diagram for ActsTrk::TrackToTrackParticleCnvAlg:

Public Member Functions
	TrackToTrackParticleCnvAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) const override
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
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
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
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
enum class	expressionStrategy { DontRecalculate , BeamLine , Vertex }
using	Stepper = Acts::EigenStepper<>
using	Navigator = Acts::Navigator
using	Propagator = Acts::Propagator<Stepper, Navigator>
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Acts::BoundTrackParameters	parametersAtPerigee (const EventContext &ctx, const typename ActsTrk::TrackContainer::ConstTrackProxy &track, const Acts::PerigeeSurface &perigee_surface) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Static Private Member Functions
static std::shared_ptr< Acts::PerigeeSurface >	makePerigeeSurface (const InDet::BeamSpotData *beamspotptr)
static std::shared_ptr< Acts::PerigeeSurface >	makePerigeeSurface (const xAOD::Vertex &)
static xAOD::ParticleHypothesis	convertParticleHypothesis (Acts::PdgParticle abs_pdg_id)
static void	initParticleHypothesisMap ()

Private Attributes
ToolHandle< ActsTrk::IExtrapolationTool >	m_extrapolationTool {this, "ExtrapolationTool", ""}
PublicToolHandle< ActsTrk::ITrackingGeometryTool >	m_trackingGeometryTool {this, "TrackingGeometryTool", "ActsTrackingGeometryTool"}
SG::ReadHandleKeyArray< ActsTrk::TrackContainer >	m_tracksContainerKey {this, "ACTSTracksLocation", {},"Track collection (ActsTrk variant)"}
SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot or empty." }
SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCacheCondObjInputKey {this,"AtlasFieldCacheCondObj","fieldCondObj", "Name of the Magnetic Field conditions object key" }
SG::ReadCondHandleKeyArray< InDetDD::SiDetectorElementCollection >	m_siDetEleCollKey {this, "SiDetectorElementCollections", {}, "Pixel and strip element collections to get geometry information about measurements."}
Gaudi::Property< std::vector< unsigned int > >	m_siDetEleCollToMeasurementType {this, "SiDetEleCollToMeasurementType",{}, "One value per si detector collection: Pixel = 1, Strip = 2"}
SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexHandle {this, "VertexContainerKey", "", "Name of the Primary Vertex Container"}
SG::WriteHandleKey< xAOD::TrackParticleContainer >	m_trackParticlesOutKey {this, "TrackParticlesOutKey","", "Name of the produced track particle collection" }
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer >	m_decorator_actsTracks {this, "ActsTrackLink", "actsTrack"}
Gaudi::Property< double >	m_paramExtrapolationParLimit {this, "ExtrapolationPathLimit",std::numeric_limits<double>::max(), "PathLimit for extrapolating track parameters." }
Gaudi::Property< bool >	m_firstAndLastParamOnly {this, "FirstAndLastParameterOnly",true, "Only convert the first and the last parameter." }
Gaudi::Property< bool >	m_computeExpectedLayerPattern {this, "ComputeExpectedLayerPattern",true, "Compute the expected layer pattern. CPU expensive" }
Gaudi::Property< bool >	m_expectIfPixelContributes {this, "expectIfPixelContribution",true, "Only expect pixel hits if there are pixel hits on track." }
Gaudi::Property< double >	m_pixelExpectLayerPathLimitInMM
Gaudi::Property< std::string >	m_perigeeExpression {this, "PerigeeExpression", "DontRecalculate"}
expressionStrategy	m_expression_strategy {expressionStrategy::BeamLine}
std::unique_ptr< Propagator >	m_propagator
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

Static Private Attributes
static std::vector< std::pair< Acts::PdgParticle, xAOD::ParticleHypothesis > > s_actsHypothesisToxAOD	ATLAS_THREAD_SAFE

Detailed Description

Definition at line 50 of file TrackToTrackParticleCnvAlg.h.

Member Typedef Documentation

Navigator

using ActsTrk::TrackToTrackParticleCnvAlg::Navigator = Acts::Navigator

private

Definition at line 63 of file TrackToTrackParticleCnvAlg.h.

Propagator

using ActsTrk::TrackToTrackParticleCnvAlg::Propagator = Acts::Propagator<Stepper, Navigator>

private

Definition at line 64 of file TrackToTrackParticleCnvAlg.h.

Stepper

using ActsTrk::TrackToTrackParticleCnvAlg::Stepper = Acts::EigenStepper<>

private

Definition at line 62 of file TrackToTrackParticleCnvAlg.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

expressionStrategy

enum class ActsTrk::TrackToTrackParticleCnvAlg::expressionStrategy

strongprivate

Enumerator
DontRecalculate
BeamLine
Vertex

Definition at line 113 of file TrackToTrackParticleCnvAlg.h.

{DontRecalculate, BeamLine, Vertex};

Constructor & Destructor Documentation

TrackToTrackParticleCnvAlg()

ActsTrk::TrackToTrackParticleCnvAlg::TrackToTrackParticleCnvAlg	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 124 of file TrackToTrackParticleCnvAlg.cxx.

      : AthReentrantAlgorithm(name, pSvcLocator)
  {
  }

Member Function Documentation

cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

convertParticleHypothesis()

xAOD::ParticleHypothesis ActsTrk::TrackToTrackParticleCnvAlg::convertParticleHypothesis ( Acts::PdgParticle abs_pdg_id )

inlinestaticprivate

Definition at line 125 of file TrackToTrackParticleCnvAlg.h.

                                                                                                                     {
   std::vector<std::pair<Acts::PdgParticle, xAOD::ParticleHypothesis> >::const_iterator
      iter = std::find_if(s_actsHypothesisToxAOD.begin(),
                          s_actsHypothesisToxAOD.end(),
                          [abs_pdg_id](const std::pair<Acts::PdgParticle, xAOD::ParticleHypothesis> &elm) {
                             return abs_pdg_id == elm.first;
                          });
   return  (iter != s_actsHypothesisToxAOD.end() ? iter->second : xAOD::noHypothesis);
}

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::TrackToTrackParticleCnvAlg::execute ( const EventContext & ctx ) const

overridevirtual

Definition at line 201 of file TrackToTrackParticleCnvAlg.cxx.

  {
    SG::WriteHandle<xAOD::TrackParticleContainer> wh_track_particles( m_trackParticlesOutKey, ctx);
    if (wh_track_particles.record(std::make_unique<xAOD::TrackParticleContainer>(),
                                  std::make_unique<xAOD::TrackParticleAuxContainer>()).isFailure()) {
       ATH_MSG_ERROR("Failed to record track particle container with key " << m_trackParticlesOutKey.key() );
       return StatusCode::FAILURE;
    }
 
    SG::WriteDecorHandle<xAOD::TrackParticleContainer, ElementLink<ActsTrk::TrackContainer>> trackLink(m_decorator_actsTracks, ctx);
 
    xAOD::TrackParticleContainer *track_particles = wh_track_particles.ptr();
 
    const InDet::BeamSpotData *beamspot_data {nullptr};
    const xAOD::VertexContainer *vertexContainer {nullptr};
    const xAOD::Vertex* primaryVertex {nullptr};
 
    if (m_expression_strategy == expressionStrategy::BeamLine) {
      SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle = SG::makeHandle( m_beamSpotKey, ctx );
      ATH_CHECK(beamSpotHandle.isValid());
      beamspot_data = beamSpotHandle.cptr();
    } 
 
    if (m_expression_strategy == expressionStrategy::Vertex) {
      SG::ReadHandle<xAOD::VertexContainer> vertexHandle = SG::makeHandle( m_vertexHandle, ctx );
      ATH_CHECK( vertexHandle.isValid() );
      vertexContainer = vertexHandle.cptr();
      if (vertexContainer->size() == 0) {
    ATH_MSG_ERROR("Retrieved an empty vertex container. This is totally wrong!");
    return StatusCode::FAILURE;
      }
      
      for(const xAOD::Vertex* vtx : *vertexContainer) {
         if(vtx->vertexType() == xAOD::VxType::PriVtx) {
       primaryVertex = vtx;
           break;
         }
      }
 
      if (not primaryVertex) {
    ATH_MSG_WARNING("Requested to compute track particles wrt primary vertex, but no primary vertex is found. Using dummy vertex");
    primaryVertex = vertexContainer->front();
      }
    }
 
    std::size_t nTracks = 0ul;
    std::vector<const ActsTrk::TrackContainer *> trackContainers;
    for (const SG::ReadHandleKey<ActsTrk::TrackContainer>& handleKey : m_tracksContainerKey) {
      SG::ReadHandle<ActsTrk::TrackContainer> handle = SG::makeHandle( handleKey, ctx );
      ATH_CHECK(handle.isValid());
      trackContainers.push_back( handle.cptr() );
      nTracks += trackContainers.back()->size();
    }
 
    // Fast Insertion Trick
    std::vector<xAOD::TrackParticle*> toAddParticles;
    toAddParticles.reserve(nTracks);
    for (std::size_t i(0); i<nTracks; ++i) {
      toAddParticles.push_back( new xAOD::TrackParticle() );
    }
    track_particles->insert(track_particles->end(),
                toAddParticles.begin(),
                toAddParticles.end());
    
    SG::ReadCondHandle<AtlasFieldCacheCondObj> fieldHandle = SG::makeHandle( m_fieldCacheCondObjInputKey, ctx );
    ATH_CHECK(fieldHandle.isValid());
    const AtlasFieldCacheCondObj *field_cond_data = fieldHandle.cptr();
    MagField::AtlasFieldCache fieldCache;
    field_cond_data->getInitializedCache(fieldCache);
 
    const GeometryContext &gctx = m_trackingGeometryTool->getNominalGeometryContext();
    std::shared_ptr<Acts::PerigeeSurface> perigee_surface {nullptr};
    if (m_expression_strategy == expressionStrategy::BeamLine) {
      perigee_surface = makePerigeeSurface(beamspot_data);
    } else if (m_expression_strategy == expressionStrategy::Vertex) {
      perigee_surface = makePerigeeSurface(*primaryVertex);
    }
 
    std::array<const InDetDD::SiDetectorElementCollection *,ActsTrk::detail::to_underlying(xAOD::UncalibMeasType::nTypes)> siDetEleColl {};
    for (unsigned int idx=0; idx <m_siDetEleCollToMeasurementType.size(); ++idx ) {
      SG::ReadCondHandle<InDetDD::SiDetectorElementCollection> detHandle = SG::makeHandle( m_siDetEleCollKey[idx], ctx );
      ATH_CHECK(detHandle.isValid());
      siDetEleColl[m_siDetEleCollToMeasurementType[idx] ] = detHandle.cptr();
    }
 
    static const std::array<unsigned short, ActsTrk::detail::to_underlying(xAOD::UncalibMeasType::nTypes)>
      measurementToSummaryType ATLAS_THREAD_SAFE (makeMeasurementToSummaryTypeMap());
 
 
    
    // re-used temporaries
    std::vector<float> tmp_cov_vector;
    std::vector<ActsTrk::TrackStateBackend::ConstTrackStateProxy::IndexType > tmp_param_state_idx;
    tmp_param_state_idx.reserve(30);
    Amg::Vector3D magnFieldVect;
    std::vector<std::vector<float>> parametersVec;
    ActsTrk::detail::HitSummaryData hitInfo;
 
    unsigned int converted_track_states=0;
 
    std::size_t particleCounter = 0ul;
    using namespace Acts::UnitLiterals;
    for (const ActsTrk::TrackContainer *tracksContainer : trackContainers) {
 
      bool precalculatedLayerPattern = detail::ExpectedLayerPatternHelper::exists(*tracksContainer);
 
      for (const typename ActsTrk::TrackContainer::ConstTrackProxy track : *tracksContainer) {
    xAOD::TrackParticle *track_particle = track_particles->at(particleCounter++);
    
    // convert defining parameters
    // @TODO add support for other modes available in the legacy converter : wrt a vertex, origin, beamspot ?
  Acts::BoundTrackParameters perigeeParam = [&] {
        if (m_expression_strategy == expressionStrategy::DontRecalculate) {
          // If the strategy is "DontRecalculate", we will take the reference surface as is 
          // from the track finding without modification. Consult track finding configuration to
          // find out what that is.
          return track.createParametersAtReference();
        }
        else {
          return parametersAtPerigee(ctx, track, *perigee_surface);
        }
  }();
 
  Acts::BoundVector boundParams = perigeeParam.parameters();
  track_particle->setDefiningParameters(boundParams[Acts::eBoundLoc0],
                                        boundParams[Acts::eBoundLoc1],
                                        boundParams[Acts::eBoundPhi],
                                        boundParams[Acts::eBoundTheta],
                                        boundParams[Acts::eBoundQOverP] * 1_MeV);
 
    if (perigeeParam.covariance().has_value()) {
          // only use the 5x5 sub-matrix of the full covariance matrix
          lowerTriangleToVectorScaleLastRow(perigeeParam.covariance().value(),tmp_cov_vector,5, 1_MeV);
          track_particle->setDefiningParametersCovMatrixVec(tmp_cov_vector);
    }
    // optional beam tilt
    if (beamspot_data) {
      track_particle->setBeamlineTiltX(beamspot_data->beamTilt(0));
      track_particle->setBeamlineTiltY(beamspot_data->beamTilt(1));
    }
    
    // fit info, quality
    track_particle->setFitQuality(track.chi2(), track.nDoF());
    track_particle->setPatternRecognitionInfo( (1ul << xAOD::SiSPSeededFinder) );
    track_particle->setTrackFitter(xAOD::KalmanFitter);
    
    const Acts::ParticleHypothesis &hypothesis = track.particleHypothesis();
    track_particle->setParticleHypothesis(convertParticleHypothesis( hypothesis.absolutePdg() ));
    constexpr float inv_1_MeV = 1/1_MeV;
    // gather track state indices for parameter conversion
    // @TODO add support for muons
    
    // xAOD::UncalibMeasType::underlying_type is expected to be the number of UncalibMeasTypes
    std::array<std::array<uint8_t, ActsTrk::detail::to_underlying(ActsTrk::detail::HitCategory::N)>,
           ActsTrk::detail::to_underlying(xAOD::UncalibMeasType::nTypes)> specialHitCounts{};
    
    ActsTrk::detail::SumOfValues chi2_stat;
    gatherTrackSummaryData(*tracksContainer,
                   track,
                   siDetEleColl,
                   measurementToSummaryType,
                   chi2_stat,
                   hitInfo,
                   tmp_param_state_idx,
                   specialHitCounts);
    
    // Muon
    //    MdtDriftCircleType = 3
    //    RpcStripType = 4,
    //    TgcStripType = 5,
    //    MMClusterType = 6,
    //    sTgcStripType = 7,
    
    // pixel summaries
    std::array< std::tuple< uint8_t, uint8_t, uint8_t, bool >, 4> copy_summary {
          std::make_tuple(static_cast<uint8_t>(ActsTrk::detail::HitSummaryData::pixelTotal),
                          static_cast<uint8_t>(xAOD::numberOfContribPixelLayers),
                          static_cast<uint8_t>(xAOD::numberOfPixelHits),
                          false),
      
          std::make_tuple(static_cast<uint8_t>(ActsTrk::detail::HitSummaryData::pixelBarrelFlat),
                          static_cast<uint8_t>(xAOD::numberOfContribPixelBarrelFlatLayers),
                          static_cast<uint8_t>(xAOD::numberOfPixelBarrelFlatHits),
                          true),
 
          std::make_tuple(static_cast<uint8_t>(ActsTrk::detail::HitSummaryData::pixelBarrelInclined),
                          static_cast<uint8_t>(xAOD::numberOfContribPixelBarrelInclinedLayers),
                          static_cast<uint8_t>(xAOD::numberOfPixelBarrelInclinedHits),
                          true),
 
          std::make_tuple(static_cast<uint8_t>(ActsTrk::detail::HitSummaryData::pixelEndcap),
                          static_cast<uint8_t>(xAOD::numberOfContribPixelEndcap),
                          static_cast<uint8_t>(xAOD::numberOfPixelEndcapHits),\
                          true) };
    
    for (auto [src_region, dest_xaod_summary_layer, dest_xaod_summary_hits, add_outlier] : copy_summary ) {
          setSummaryValue(*track_particle,
                          hitInfo.contributingLayers( static_cast<ActsTrk::detail::HitSummaryData::DetectorRegion>(src_region)),
                          static_cast<xAOD::SummaryType>(dest_xaod_summary_layer));
          setSummaryValue(*track_particle,
              hitInfo.contributingHits(static_cast<ActsTrk::detail::HitSummaryData::DetectorRegion>(src_region))
                          + ( add_outlier
                              ? hitInfo.contributingOutlierHits(static_cast<ActsTrk::detail::HitSummaryData::DetectorRegion>(src_region))
                              : 0),
                          static_cast<xAOD::SummaryType>(dest_xaod_summary_hits));
    }
    setSummaryValue(*track_particle,
                        hitInfo.sum<ActsTrk::detail::HitSummaryData::Hit>(ActsTrk::detail::HitSummaryData::pixelEndcap,0)
            +hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelEndcap,0),
            xAOD::numberOfInnermostPixelLayerEndcapHits);
    setSummaryValue(*track_particle,
            hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelEndcap,0),
            xAOD::numberOfInnermostPixelLayerEndcapOutliers);
    setSummaryValue(*track_particle,
                        hitInfo.sum<ActsTrk::detail::HitSummaryData::Hit>(ActsTrk::detail::HitSummaryData::pixelEndcap,1)
            +hitInfo.sum<ActsTrk::detail::HitSummaryData::Hit>(ActsTrk::detail::HitSummaryData::pixelEndcap,2)
            +hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelEndcap,1)
            +hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelEndcap,2),
            xAOD::numberOfNextToInnermostPixelLayerEndcapHits);
    setSummaryValue(*track_particle,
                        hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelEndcap,1)
            +hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelEndcap,2),
            xAOD::numberOfNextToInnermostPixelLayerEndcapOutliers);
    setSummaryValue(*track_particle,
            hitInfo.contributingOutlierHits(ActsTrk::detail::HitSummaryData::pixelTotal),
            xAOD::numberOfPixelOutliers);
    setSummaryValue(*track_particle,
            specialHitCounts[ActsTrk::detail::to_underlying(xAOD::UncalibMeasType::PixelClusterType)][ActsTrk::detail::HitCategory::Hole],
            xAOD::numberOfPixelHoles);
    setSummaryValue(*track_particle,
            hitInfo.sum<ActsTrk::detail::HitSummaryData::SharedHit>(ActsTrk::detail::HitSummaryData::pixelEndcap, 0),
            xAOD::numberOfInnermostPixelLayerSharedEndcapHits);
    setSummaryValue(*track_particle,
            hitInfo.sum<ActsTrk::detail::HitSummaryData::SharedHit>(ActsTrk::detail::HitSummaryData::pixelEndcap, 1)
            +hitInfo.sum<ActsTrk::detail::HitSummaryData::SharedHit>(ActsTrk::detail::HitSummaryData::pixelEndcap, 2),
            xAOD::numberOfNextToInnermostPixelLayerSharedEndcapHits);
    setSummaryValue(*track_particle,
            hitInfo.contributingSharedHits(ActsTrk::detail::HitSummaryData::pixelTotal),
            xAOD::numberOfPixelSharedHits);
 
 
  std::array<unsigned int,4> expect_layer_pattern{};
  if (precalculatedLayerPattern) {
    // We have a pre-calculated layer pattern from track finding, use as is
    expect_layer_pattern = detail::ExpectedLayerPatternHelper::get(track);
  }
  else {
    // Only check if computeExpectedLayerPattern is true. TODO:: move this computation to the track finding to avoid calling propagator steps here.
    // Do not expect pixel hits if there are not contributing pixel hits in the flat barrel and expectIfPixelContributes is true
    expect_layer_pattern = (m_computeExpectedLayerPattern.value()
                  && (!m_expectIfPixelContributes.value()
                      || hitInfo.contributingLayers(ActsTrk::detail::HitSummaryData::pixelTotal))
                  ? detail::expectedLayerPattern(ctx,
                          *m_extrapolationTool,
                          perigeeParam,
                          m_pixelExpectLayerPathLimitInMM.value() * Acts::UnitConstants::mm)
                  : std::array<unsigned int,4> {0u,0u, 0u,0u} );
  }
    
    // @TODO consider end-caps  for inner most pixel hits ?
    setSummaryValue(*track_particle,
            static_cast<uint8_t>((expect_layer_pattern[0] & (1<<0)) != 0 ),
            xAOD::expectInnermostPixelLayerHit);
    setSummaryValue(*track_particle,
            static_cast<uint8_t>((expect_layer_pattern[0] & (1<<1)) != 0 ),
            xAOD::expectNextToInnermostPixelLayerHit);
    setSummaryValue(*track_particle,
            static_cast<unsigned int >(hitInfo.sum<ActsTrk::detail::HitSummaryData::Hit>(ActsTrk::detail::HitSummaryData::pixelBarrelFlat,0)),
            xAOD::numberOfInnermostPixelLayerHits);
    setSummaryValue(*track_particle,
            static_cast<unsigned int >(hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelBarrelFlat,0)),
            xAOD::numberOfInnermostPixelLayerOutliers);
    setSummaryValue(*track_particle,
            static_cast<unsigned int >(hitInfo.sum<ActsTrk::detail::HitSummaryData::Hit>(ActsTrk::detail::HitSummaryData::pixelBarrelFlat,1)),
            xAOD::numberOfNextToInnermostPixelLayerHits);
    setSummaryValue(*track_particle,
            static_cast<unsigned int >(hitInfo.sum<ActsTrk::detail::HitSummaryData::Outlier>(ActsTrk::detail::HitSummaryData::pixelBarrelFlat,1)),
            xAOD::numberOfNextToInnermostPixelLayerOutliers);
    setSummaryValue(*track_particle,
            static_cast<unsigned int >(hitInfo.sum<ActsTrk::detail::HitSummaryData::SharedHit>(ActsTrk::detail::HitSummaryData::pixelBarrelFlat,0)),
            xAOD::numberOfInnermostPixelLayerSharedHits);
    setSummaryValue(*track_particle,
            static_cast<unsigned int >(hitInfo.sum<ActsTrk::detail::HitSummaryData::SharedHit>(ActsTrk::detail::HitSummaryData::pixelBarrelFlat,1)),
            xAOD::numberOfNextToInnermostPixelLayerSharedHits);
    
    // Strip summaries
    setSummaryValue(*track_particle,
            hitInfo.contributingHits( ActsTrk::detail::HitSummaryData::stripTotal ),
            xAOD::numberOfSCTHits);
    setSummaryValue(*track_particle,
            hitInfo.contributingOutlierHits( ActsTrk::detail::HitSummaryData::stripTotal ),
            xAOD::numberOfSCTOutliers);
    setSummaryValue(*track_particle,
            hitInfo.contributingSharedHits( ActsTrk::detail::HitSummaryData::stripTotal ),
            xAOD::numberOfSCTSharedHits);
    setSummaryValue(*track_particle,
            specialHitCounts[ActsTrk::detail::to_underlying(xAOD::UncalibMeasType::StripClusterType)][ActsTrk::detail::HitCategory::Hole],
            xAOD::numberOfSCTHoles);
    
    double biased_chi2_variance = chi2_stat.biasedVariance();
    setSummaryValue(*track_particle,
            static_cast<uint8_t> (biased_chi2_variance>0.
                          ? std::min(static_cast<unsigned int>(std::sqrt(biased_chi2_variance) * 100),255u)
                          : 0u),
            xAOD::standardDeviationOfChi2OS);
    
    setSummaryValue(*track_particle,
                        hitInfo.contributingOutlierHits( ActsTrk::detail::HitSummaryData::pixelTotal )
            +hitInfo.contributingOutlierHits( ActsTrk::detail::HitSummaryData::stripTotal ),
            xAOD::numberOfOutliersOnTrack);
    
    
    // @TODO slect states for which parameters are stored
    if (m_firstAndLastParamOnly && tmp_param_state_idx.size()>2) {
          tmp_param_state_idx[1]=tmp_param_state_idx.back();
          tmp_param_state_idx.erase(tmp_param_state_idx.begin()+2,tmp_param_state_idx.end());
    }
 
    // store track parameters and covariances for slected states
    parametersVec.clear();
    parametersVec.reserve(tmp_param_state_idx.size());
    
    for(std::vector<ActsTrk::TrackStateBackend::ConstTrackStateProxy::IndexType>::const_reverse_iterator
              idx_iter = tmp_param_state_idx.rbegin();
        idx_iter != tmp_param_state_idx.rend();
        ++idx_iter) {
          //       for(ActsTrk::TrackStateBackend::ConstTrackStateProxy::IndexType idx : tmp_param_state_idx) {
          ActsTrk::TrackStateBackend::ConstTrackStateProxy
        state = tracksContainer->trackStateContainer().getTrackState(*idx_iter);
          const Acts::BoundTrackParameters actsParam = track.createParametersFromState(state);
      
          Acts::Vector3 position = actsParam.position(gctx.context());
          Acts::Vector3 momentum = actsParam.momentum();
      
          // scaling from Acts momentume units (GeV) to Athena Units (MeV)
          for (unsigned int i=0; i<momentum.rows(); ++i) {
             momentum(i) *= inv_1_MeV;
          }
 
 
          if (actsParam.covariance()) {
             Acts::MagneticFieldContext mfContext = m_extrapolationTool->getMagneticFieldContext(ctx);
             Acts::GeometryContext tgContext = gctx.context();
 
             magnFieldVect.setZero();
             fieldCache.getField(position.data(), magnFieldVect.data());
             // scaling from Athena magnetic field units kT to Acts units T
             {
                using namespace Acts::UnitLiterals;
                magnFieldVect *= 1000_T;
             }
 
             auto curvilinear_cov_result = ActsTrk::detail::convertActsBoundCovToCurvilinearParam(tgContext, actsParam, magnFieldVect, hypothesis);
             if (curvilinear_cov_result.has_value()) {
                Acts::BoundSquareMatrix &curvilinear_cov = curvilinear_cov_result.value();
 
                // convert q/p components from GeV (Acts) to MeV (Athena)
                for (unsigned int col_i=0; col_i<4; ++col_i) {
                   curvilinear_cov(col_i,4) *= 1_MeV;
                   curvilinear_cov(4,col_i) *= 1_MeV;
                }
                curvilinear_cov(4,4) *= (1_MeV * 1_MeV);
 
                std::size_t param_idx = parametersVec.size();
                // only use the 5x5 sub-matrix of the full covariance matrix
                lowerTriangleToVector(curvilinear_cov,tmp_cov_vector,5);
                if (tmp_cov_vector.size() != 15) {
                   ATH_MSG_ERROR("Invalid size of lower triangle cov " << tmp_cov_vector.size() <<  " != 15" 
                                 << " input matrix : " << curvilinear_cov.rows() << " x " << curvilinear_cov.cols() );
                }
                track_particle->setTrackParameterCovarianceMatrix(param_idx, tmp_cov_vector);
             }
          }
          parametersVec.emplace_back(std::vector<float>{
          static_cast<float>(position[0]),static_cast<float>(position[1]),static_cast<float>(position[2]),
          static_cast<float>(momentum[0]),static_cast<float>(momentum[1]),static_cast<float>(momentum[2]) });
          ++converted_track_states;
      
      
    }
    for (const std::vector<float> &param : parametersVec) {
          if (param.size() != 6) {
        ATH_MSG_ERROR("Invalid size of param element " << param.size() <<  " != 6" );
          }
    }
    
    track_particle->setTrackParameters(parametersVec);
    
    // add element to link to the correspond track
    trackLink(*track_particle)
          = ElementLink<ActsTrk::TrackContainer>( tracksContainer,
                                                  track.index() );
    ATH_CHECK( (trackLink(*track_particle)).isValid() );
      }
    }
    ATH_MSG_DEBUG( "Converted " << nTracks << " acts tracks into " << track_particles->size()
           << " track particles with parameters for " << converted_track_states << " track states.");
    
    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 & AthCommonReentrantAlgorithm< 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 94 of file AthCommonReentrantAlgorithm.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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

initialize()

StatusCode ActsTrk::TrackToTrackParticleCnvAlg::initialize ( )

overridevirtual

Definition at line 130 of file TrackToTrackParticleCnvAlg.cxx.

  {
     std::vector<std::string> supportedStrategies {"BeamLine", "Vertex"};
     bool isAllowedStrategy = false;
     for (const std::string& strategy : supportedStrategies) {
       if (m_perigeeExpression != strategy) continue;
       isAllowedStrategy = true;
       break;
     }
     ATH_MSG_DEBUG("- perigeeExpression: " << m_perigeeExpression.value());
     if (not isAllowedStrategy) {
       ATH_MSG_ERROR("Wrong configuration of the Track to Track Particle Cnv algorithm: perigeeExpression is not supported");
       return StatusCode::FAILURE;
     }
 
     if (m_perigeeExpression == "BeamLine") m_expression_strategy = expressionStrategy::BeamLine;
     else if (m_perigeeExpression == "Vertex") m_expression_strategy = expressionStrategy::Vertex;
     else if (m_perigeeExpression == "DontRecalculate") m_expression_strategy = expressionStrategy::DontRecalculate;
     else return StatusCode::FAILURE;
    
     ATH_CHECK(m_trackingGeometryTool.retrieve());
     ATH_CHECK( m_tracksContainerKey.initialize() );
     ATH_CHECK( m_trackParticlesOutKey.initialize() );
     ATH_CHECK( m_beamSpotKey.initialize(m_expression_strategy == expressionStrategy::BeamLine) );
     ATH_CHECK( m_vertexHandle.initialize(m_expression_strategy == expressionStrategy::Vertex) );
     ATH_CHECK( m_fieldCacheCondObjInputKey.initialize() );
 
     m_decorator_actsTracks = m_trackParticlesOutKey.key() + "." + m_decorator_actsTracks.key();
     ATH_CHECK(m_decorator_actsTracks.initialize());
     
     ATH_CHECK( m_extrapolationTool.retrieve() ); // for extrapolation to beamline
 
     // propagator for conversion to curvilnear parameters
     {
        auto logger = makeActsAthenaLogger(this, "Prop");
 
        Navigator::Config cfg{m_trackingGeometryTool->trackingGeometry()};
        cfg.resolvePassive = false;
        cfg.resolveMaterial = true;
        cfg.resolveSensitive = true;
        auto navigtor_logger = logger->cloneWithSuffix("Navigator");
        m_propagator = std::make_unique<Propagator>(Stepper(std::make_shared<ATLASMagneticFieldWrapper>()),
                                                    Navigator(cfg,std::move(navigtor_logger)),
                                                    std::move(logger));
     }
 
     // for layer/module-type information per hit
     ATH_CHECK( m_siDetEleCollKey.initialize() );
     if (m_siDetEleCollToMeasurementType.size() == m_siDetEleCollKey.size()) {
        unsigned int collection_idx=0;
        for (int type : m_siDetEleCollToMeasurementType ) {
           if (type <1 || type >2) {
              ATH_MSG_ERROR("Invalid measurement type (" << type << ") given for collection "  << collection_idx << " : "
                            << m_siDetEleCollKey[collection_idx].key()
                            << ". Expected 1 for pixel, 2 for strips.");
              return StatusCode::FAILURE;
           }
           ++collection_idx;
        }
     }
     else {
        ATH_MSG_ERROR("Expected exactly one value in SiDetEleCollToMeasurementType per SiDetectorElementCollection. But got "
                      << m_siDetEleCollToMeasurementType.size() << " instead of " << m_siDetEleCollKey.size() << ".");
        return StatusCode::FAILURE;
     }
 
     initParticleHypothesisMap();
     
     return StatusCode::SUCCESS;
  }

initParticleHypothesisMap()

void ActsTrk::TrackToTrackParticleCnvAlg::initParticleHypothesisMap ( )

staticprivate

Definition at line 110 of file TrackToTrackParticleCnvAlg.cxx.

                                                             {
      if (s_actsHypothesisToxAOD.empty()) {
         s_actsHypothesisToxAOD.reserve(7);
         s_actsHypothesisToxAOD.push_back( std::make_pair( Acts::eElectron , xAOD::electron) );
         s_actsHypothesisToxAOD.push_back( std::make_pair( Acts::eMuon , xAOD::muon) );
         s_actsHypothesisToxAOD.push_back( std::make_pair( Acts::ePionPlus , xAOD::pion) );
         s_actsHypothesisToxAOD.push_back( std::make_pair( Acts::eProton , xAOD::proton) );
         s_actsHypothesisToxAOD.push_back( std::make_pair( Acts::ePionZero , xAOD::pi0) );
         s_actsHypothesisToxAOD.push_back( std::make_pair( Acts::eNeutron , xAOD::neutron) );
         s_actsHypothesisToxAOD.push_back( std::make_pair( Acts::eGamma , xAOD::photon) );
     }
   }

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()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, 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, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

makePerigeeSurface() [1/2]

std::shared_ptr< Acts::PerigeeSurface > ActsTrk::TrackToTrackParticleCnvAlg::makePerigeeSurface ( const InDet::BeamSpotData * beamspotptr )

staticprivate

Definition at line 602 of file TrackToTrackParticleCnvAlg.cxx.

                                                                                                                         {
     // @from TrackToVertex::trackAtBeamline
     Acts::Vector3 beamspot(0., 0., 0.);
     float tiltx = 0.0;
     float tilty = 0.0;
     if (beamspot_data) {
        beamspot = Acts::Vector3(beamspot_data->beamVtx().position());
        tiltx =  beamspot_data->beamTilt(0);
        tilty =  beamspot_data->beamTilt(1);
     }
     Acts::Translation3 translation(beamspot);
     Acts::Transform3 transform( translation * Acts::RotationMatrix3::Identity() );
     transform *= Acts::AngleAxis3(tilty, Acts::Vector3(0.,1.,0.));
     transform *= Acts::AngleAxis3(tiltx, Acts::Vector3(1.,0.,0.));
     return Acts::Surface::makeShared<Acts::PerigeeSurface>(transform);
  }

makePerigeeSurface() [2/2]

std::shared_ptr< Acts::PerigeeSurface > ActsTrk::TrackToTrackParticleCnvAlg::makePerigeeSurface ( const xAOD::Vertex & vertex )

staticprivate

Definition at line 619 of file TrackToTrackParticleCnvAlg.cxx.

                                                                                                           {
    Acts::Translation3 translation(Acts::Vector3(vertex.position()));
    Acts::Transform3 transform( translation * Acts::RotationMatrix3::Identity() );
    return Acts::Surface::makeShared<Acts::PerigeeSurface>(transform);    
  }

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.

parametersAtPerigee()

Acts::BoundTrackParameters ActsTrk::TrackToTrackParticleCnvAlg::parametersAtPerigee ( const EventContext & ctx, const typename ActsTrk::TrackContainer::ConstTrackProxy & track, const Acts::PerigeeSurface & perigee_surface ) const

private

Definition at line 625 of file TrackToTrackParticleCnvAlg.cxx.

                                                                                                         {
     const Acts::BoundTrackParameters trackParam = track.createParametersAtReference();
 
     std::optional<const Acts::BoundTrackParameters>
        perigeeParam = m_extrapolationTool->propagate(ctx,
                                                      trackParam,
                                                      perigee_surface,
                                                      Acts::Direction::Backward(), // @TODO try forward if backward fails ?
                                                      m_paramExtrapolationParLimit.value());
     if (!perigeeParam.has_value()) {
        ATH_MSG_WARNING( "Failed to extrapolate to perigee, started from \n" << trackParam << " " << trackParam.referenceSurface().name() );
 
        return trackParam;
     }
     else {
        ATH_MSG_DEBUG( "Succeeded to extrapolate to perigee ");
        return perigeeParam.value();
     }
  }

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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

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

sysExecute()

StatusCode AthCommonReentrantAlgorithm< 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.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonReentrantAlgorithm< 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 HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.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

ATLAS_THREAD_SAFE

std::vector<std::pair<Acts::PdgParticle, xAOD::ParticleHypothesis> > s_actsHypothesisToxAOD ActsTrk::TrackToTrackParticleCnvAlg::ATLAS_THREAD_SAFE

staticprivate

Definition at line 118 of file TrackToTrackParticleCnvAlg.h.

m_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> ActsTrk::TrackToTrackParticleCnvAlg::m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot or empty." }

private

Definition at line 80 of file TrackToTrackParticleCnvAlg.h.

{this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot or empty." };

m_computeExpectedLayerPattern

Gaudi::Property<bool> ActsTrk::TrackToTrackParticleCnvAlg::m_computeExpectedLayerPattern {this, "ComputeExpectedLayerPattern",true, "Compute the expected layer pattern. CPU expensive" }

private

Definition at line 102 of file TrackToTrackParticleCnvAlg.h.

{this, "ComputeExpectedLayerPattern",true, "Compute the expected layer pattern. CPU expensive" };

m_decorator_actsTracks

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> ActsTrk::TrackToTrackParticleCnvAlg::m_decorator_actsTracks {this, "ActsTrackLink", "actsTrack"}

private

Definition at line 95 of file TrackToTrackParticleCnvAlg.h.

{this, "ActsTrackLink", "actsTrack"};

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_expectIfPixelContributes

Gaudi::Property<bool> ActsTrk::TrackToTrackParticleCnvAlg::m_expectIfPixelContributes {this, "expectIfPixelContribution",true, "Only expect pixel hits if there are pixel hits on track." }

private

Definition at line 104 of file TrackToTrackParticleCnvAlg.h.

{this, "expectIfPixelContribution",true, "Only expect pixel hits if there are pixel hits on track." };

m_expression_strategy

expressionStrategy ActsTrk::TrackToTrackParticleCnvAlg::m_expression_strategy {expressionStrategy::BeamLine}

private

Definition at line 114 of file TrackToTrackParticleCnvAlg.h.

{expressionStrategy::BeamLine};

m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_extrapolationTool

ToolHandle<ActsTrk::IExtrapolationTool> ActsTrk::TrackToTrackParticleCnvAlg::m_extrapolationTool {this, "ExtrapolationTool", ""}

private

Definition at line 72 of file TrackToTrackParticleCnvAlg.h.

{this, "ExtrapolationTool", ""};

m_fieldCacheCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> ActsTrk::TrackToTrackParticleCnvAlg::m_fieldCacheCondObjInputKey {this,"AtlasFieldCacheCondObj","fieldCondObj", "Name of the Magnetic Field conditions object key" }

private

Definition at line 82 of file TrackToTrackParticleCnvAlg.h.

{this,"AtlasFieldCacheCondObj","fieldCondObj", "Name of the Magnetic Field conditions object key" };

m_firstAndLastParamOnly

Gaudi::Property<bool> ActsTrk::TrackToTrackParticleCnvAlg::m_firstAndLastParamOnly {this, "FirstAndLastParameterOnly",true, "Only convert the first and the last parameter." }

private

Definition at line 100 of file TrackToTrackParticleCnvAlg.h.

{this, "FirstAndLastParameterOnly",true, "Only convert the first and the last parameter." };

m_paramExtrapolationParLimit

Gaudi::Property<double> ActsTrk::TrackToTrackParticleCnvAlg::m_paramExtrapolationParLimit {this, "ExtrapolationPathLimit",std::numeric_limits<double>::max(), "PathLimit for extrapolating track parameters." }

private

Definition at line 98 of file TrackToTrackParticleCnvAlg.h.

{this, "ExtrapolationPathLimit",std::numeric_limits<double>::max(), "PathLimit for extrapolating track parameters." }; // @TODO (unit?mm?)

m_perigeeExpression

Gaudi::Property<std::string> ActsTrk::TrackToTrackParticleCnvAlg::m_perigeeExpression {this, "PerigeeExpression", "DontRecalculate"}

private

Definition at line 111 of file TrackToTrackParticleCnvAlg.h.

{this, "PerigeeExpression", "DontRecalculate"};

m_pixelExpectLayerPathLimitInMM

Gaudi::Property<double> ActsTrk::TrackToTrackParticleCnvAlg::m_pixelExpectLayerPathLimitInMM

private

Initial value:

{this, "PixelExpectLayerPathLimitInMM",1000,
        "PathLimit for extrapolating to get the expected pixel layer pattern in mm." }

Definition at line 107 of file TrackToTrackParticleCnvAlg.h.

       {this, "PixelExpectLayerPathLimitInMM",1000,
        "PathLimit for extrapolating to get the expected pixel layer pattern in mm." };

m_propagator

std::unique_ptr<Propagator> ActsTrk::TrackToTrackParticleCnvAlg::m_propagator

private

Definition at line 116 of file TrackToTrackParticleCnvAlg.h.

m_siDetEleCollKey

SG::ReadCondHandleKeyArray<InDetDD::SiDetectorElementCollection> ActsTrk::TrackToTrackParticleCnvAlg::m_siDetEleCollKey {this, "SiDetectorElementCollections", {}, "Pixel and strip element collections to get geometry information about measurements."}

private

Definition at line 85 of file TrackToTrackParticleCnvAlg.h.

{this, "SiDetectorElementCollections", {}, "Pixel and strip element collections to get geometry information about measurements."};

m_siDetEleCollToMeasurementType

Gaudi::Property<std::vector<unsigned int> > ActsTrk::TrackToTrackParticleCnvAlg::m_siDetEleCollToMeasurementType {this, "SiDetEleCollToMeasurementType",{}, "One value per si detector collection: Pixel = 1, Strip = 2"}

private

Definition at line 87 of file TrackToTrackParticleCnvAlg.h.

{this, "SiDetEleCollToMeasurementType",{}, "One value per si detector collection: Pixel = 1, Strip = 2"};

m_trackingGeometryTool

PublicToolHandle<ActsTrk::ITrackingGeometryTool> ActsTrk::TrackToTrackParticleCnvAlg::m_trackingGeometryTool {this, "TrackingGeometryTool", "ActsTrackingGeometryTool"}

private

Definition at line 75 of file TrackToTrackParticleCnvAlg.h.

{this, "TrackingGeometryTool", "ActsTrackingGeometryTool"};

m_trackParticlesOutKey

SG::WriteHandleKey<xAOD::TrackParticleContainer> ActsTrk::TrackToTrackParticleCnvAlg::m_trackParticlesOutKey {this, "TrackParticlesOutKey","", "Name of the produced track particle collection" }

private

Definition at line 92 of file TrackToTrackParticleCnvAlg.h.

{this, "TrackParticlesOutKey","", "Name of the produced track particle collection" };

m_tracksContainerKey

SG::ReadHandleKeyArray<ActsTrk::TrackContainer> ActsTrk::TrackToTrackParticleCnvAlg::m_tracksContainerKey {this, "ACTSTracksLocation", {},"Track collection (ActsTrk variant)"}

private

Definition at line 78 of file TrackToTrackParticleCnvAlg.h.

{this, "ACTSTracksLocation", {},"Track collection (ActsTrk variant)"};

m_varHandleArraysDeclared

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

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vertexHandle

SG::ReadHandleKey<xAOD::VertexContainer> ActsTrk::TrackToTrackParticleCnvAlg::m_vertexHandle {this, "VertexContainerKey", "", "Name of the Primary Vertex Container"}

private

Definition at line 90 of file TrackToTrackParticleCnvAlg.h.

{this, "VertexContainerKey", "", "Name of the Primary Vertex Container"};

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:

• TrackToTrackParticleCnvAlg.h
• TrackToTrackParticleCnvAlg.cxx