ActsTrk::ActsInspectTruthContentAlg Class Reference

#include <ActsInspectTruthContentAlg.h>

Inheritance diagram for ActsTrk::ActsInspectTruthContentAlg:

Public Types
enum class	EStatClusters : std::size_t {   kNTotal , kNClustersFromPrimaries , kNClustersWithNoBarcode , kNClustersWith1Contribution ,   kNClustersWith1ValidContribution , kNClustersWith2Contribution , kNClustersWith2ValidContribution , kNClustersWith3Contribution ,   kNClustersWith3ValidContribution , kNClustersWith200kBarcode , kNStat }
enum class	EStatSeeds : std::size_t {   kNTotal , nKSeedsWith0Matches , nKSeedsWith1Matches , nKSeedsWith2Matches ,   nKSeedsSame2Matches , nKSeedsWith3Matches , nKSeedsSame3Matches , nKSeedsWith4Matches ,   nKSeedsSame4Matches , nKSeedsWith5Matches , nKSeedsSame5Matches , nKSeedsWith6Matches ,   nKSeedsSame6Matches , kNStat }
enum class	EStatTracks : std::size_t {   kNTotal , kNFullMatch , kNPerfectMatch , kNTracks0Holes ,   kNTracks1Holes , kNTracks2Holes , kNTracks3Holes , kNTracks0Outliers ,   kNTracks1Outliers , kNTracks2Outliers , kNTracks3Outliers , kNStat }
enum class	SeedType : std::size_t {   PPP , SSS , PPS , PSS ,   Others , nTypes }
enum class	TrackType : std::size_t { Main , nTypes }

Public Member Functions
	ActsInspectTruthContentAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~ActsInspectTruthContentAlg () override=default
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) const override
virtual StatusCode	finalize () 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
using	cluster_stat_t
using	seed_stat_t
using	track_stat_t
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	fillStatClusters (const xAOD::UncalibratedMeasurementContainer &container, const ActsTrk::MeasurementToTruthParticleAssociation &truth, cluster_stat_t &stat) const
StatusCode	fillStatSeeds (const ActsTrk::SeedContainer &seeds, std::array< const ActsTrk::MeasurementToTruthParticleAssociation *, s_nClusterTypes > &truths, seed_stat_t &stat) const
StatusCode	fillStatTracks (const ActsTrk::TrackContainer &tracks, std::array< const ActsTrk::MeasurementToTruthParticleAssociation *, s_nClusterTypes > &truths, track_stat_t &trackStat, cluster_stat_t &onTrackClusterStat) const
template<typename row_t, typename coll_t, typename stat_t>
StatusCode	printStatTables (const std::string &objectCollectionName, const stat_t &stat) const
template<typename stat_t>
StatusCode	copyStatTable (const stat_t &contextual, stat_t &global) const
SeedType	deduceSeedType (const ActsTrk::Seed &) const
std::string	to_string (xAOD::UncalibMeasType type) const
std::string	to_string (ActsInspectTruthContentAlg::SeedType type) const
std::string	to_string (ActsInspectTruthContentAlg::TrackType type) const
std::string	to_label (ActsInspectTruthContentAlg::EStatClusters type) const
std::string	to_label (ActsInspectTruthContentAlg::EStatSeeds type) const
std::string	to_label (ActsInspectTruthContentAlg::EStatTracks type) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
SG::ReadHandleKeyArray< xAOD::UncalibratedMeasurementContainer >	m_clusters {this, "Clusters", {}}
SG::ReadHandleKeyArray< ActsTrk::SeedContainer >	m_seeds {this, "Seeds", {}}
SG::ReadHandleKeyArray< ActsTrk::MeasurementToTruthParticleAssociation >	m_associationMap_key {this,"TruthAssociationMaps", {}}
SG::ReadHandleKeyArray< ActsTrk::TrackContainer >	m_tracks {this, "Tracks", {}}
std::mutex m_mutex	ATLAS_THREAD_SAFE {}
cluster_stat_t m_clusterStat	ATLAS_THREAD_SAFE {}
seed_stat_t m_seedStat	ATLAS_THREAD_SAFE {}
std::vector< std::pair< std::string, track_stat_t > > m_trackStat	ATLAS_THREAD_SAFE {}
std::vector< std::pair< std::string, cluster_stat_t > > m_onTrack_clusterStat	ATLAS_THREAD_SAFE {}
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 constexpr std::size_t	s_nClusterTypes = static_cast<std::size_t>(xAOD::UncalibMeasType::nTypes)
static constexpr std::size_t	s_nSeedTypes = static_cast<std::size_t>(SeedType::nTypes)

Detailed Description

Definition at line 18 of file ActsInspectTruthContentAlg.h.

Member Typedef Documentation

cluster_stat_t

using ActsTrk::ActsInspectTruthContentAlg::cluster_stat_t

private

Initial value:

 std::array<
      std::array<std::size_t, s_nClusterTypes>,
      static_cast<std::size_t>(EStatClusters::kNStat)>

Definition at line 92 of file ActsInspectTruthContentAlg.h.

seed_stat_t

using ActsTrk::ActsInspectTruthContentAlg::seed_stat_t

private

Initial value:

 std::array<
      std::array<std::size_t, s_nSeedTypes>,
      static_cast<std::size_t>(EStatSeeds::kNStat)>

Definition at line 95 of file ActsInspectTruthContentAlg.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

track_stat_t

using ActsTrk::ActsInspectTruthContentAlg::track_stat_t

private

Initial value:

 std::array<
      std::array<std::size_t, 1>,
      static_cast<std::size_t>(EStatTracks::kNStat)>

Definition at line 98 of file ActsInspectTruthContentAlg.h.

Member Enumeration Documentation

EStatClusters

enum class ActsTrk::ActsInspectTruthContentAlg::EStatClusters : std::size_t

strong

Enumerator
kNTotal
kNClustersFromPrimaries
kNClustersWithNoBarcode
kNClustersWith1Contribution
kNClustersWith1ValidContribution
kNClustersWith2Contribution
kNClustersWith2ValidContribution
kNClustersWith3Contribution
kNClustersWith3ValidContribution
kNClustersWith200kBarcode
kNStat

Definition at line 29 of file ActsInspectTruthContentAlg.h.

                             : std::size_t {
      kNTotal,
      kNClustersFromPrimaries,
      kNClustersWithNoBarcode,
      kNClustersWith1Contribution,
      kNClustersWith1ValidContribution,
      kNClustersWith2Contribution,
      kNClustersWith2ValidContribution,
      kNClustersWith3Contribution,
      kNClustersWith3ValidContribution,
      kNClustersWith200kBarcode,
      kNStat
    };

EStatSeeds

enum class ActsTrk::ActsInspectTruthContentAlg::EStatSeeds : std::size_t

strong

Enumerator
kNTotal
nKSeedsWith0Matches
nKSeedsWith1Matches
nKSeedsWith2Matches
nKSeedsSame2Matches
nKSeedsWith3Matches
nKSeedsSame3Matches
nKSeedsWith4Matches
nKSeedsSame4Matches
nKSeedsWith5Matches
nKSeedsSame5Matches
nKSeedsWith6Matches
nKSeedsSame6Matches
kNStat

Definition at line 43 of file ActsInspectTruthContentAlg.h.

                          : std::size_t {
      kNTotal,
      nKSeedsWith0Matches,
      nKSeedsWith1Matches,
      nKSeedsWith2Matches,
      nKSeedsSame2Matches,
      nKSeedsWith3Matches,
      nKSeedsSame3Matches,
      nKSeedsWith4Matches,
      nKSeedsSame4Matches,
      nKSeedsWith5Matches,
      nKSeedsSame5Matches,
      nKSeedsWith6Matches,
      nKSeedsSame6Matches,
      kNStat
    };

EStatTracks

enum class ActsTrk::ActsInspectTruthContentAlg::EStatTracks : std::size_t

strong

Enumerator
kNTotal
kNFullMatch
kNPerfectMatch
kNTracks0Holes
kNTracks1Holes
kNTracks2Holes
kNTracks3Holes
kNTracks0Outliers
kNTracks1Outliers
kNTracks2Outliers
kNTracks3Outliers
kNStat

Definition at line 60 of file ActsInspectTruthContentAlg.h.

                           : std::size_t {
      kNTotal,
      kNFullMatch,
      kNPerfectMatch,
      kNTracks0Holes,
      kNTracks1Holes,
      kNTracks2Holes,
      kNTracks3Holes,
      kNTracks0Outliers,
      kNTracks1Outliers,
      kNTracks2Outliers,
      kNTracks3Outliers,
      kNStat
    };

SeedType

enum class ActsTrk::ActsInspectTruthContentAlg::SeedType : std::size_t

strong

Enumerator
PPP
SSS
PPS
PSS
Others
nTypes

Definition at line 75 of file ActsInspectTruthContentAlg.h.

                        : std::size_t {
      PPP,
      SSS,
      PPS,
      PSS,
      Others,
      nTypes
    };

TrackType

enum class ActsTrk::ActsInspectTruthContentAlg::TrackType : std::size_t

strong

Enumerator
Main
nTypes

Definition at line 84 of file ActsInspectTruthContentAlg.h.

                         : std::size_t {
      Main,
      nTypes
    };

Constructor & Destructor Documentation

ActsInspectTruthContentAlg()

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

Definition at line 18 of file ActsInspectTruthContentAlg.cxx.

    : AthReentrantAlgorithm(name, pSvcLocator)
  {}

~ActsInspectTruthContentAlg()

virtual ActsTrk::ActsInspectTruthContentAlg::~ActsInspectTruthContentAlg ( )

overridevirtualdefault

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

copyStatTable()

template<typename stat_t>

StatusCode ActsTrk::ActsInspectTruthContentAlg::copyStatTable ( const stat_t & contextual, stat_t & global ) const

private

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

deduceSeedType()

ActsInspectTruthContentAlg::SeedType ActsTrk::ActsInspectTruthContentAlg::deduceSeedType ( const ActsTrk::Seed & seed ) const

private

Definition at line 278 of file ActsInspectTruthContentAlg.cxx.

                                                                                                             {
    assert(seed.sp().size() == 3ul);
    const auto& bottom = seed.sp().at(0);
    const auto& middle = seed.sp().at(1);
    const auto& top = seed.sp().at(2);
    xAOD::UncalibMeasType bottom_type = bottom->measurements().front()->type();
    xAOD::UncalibMeasType middle_type = middle->measurements().front()->type();
    xAOD::UncalibMeasType top_type = top->measurements().front()->type();
 
    if (bottom_type == xAOD::UncalibMeasType::PixelClusterType and
    middle_type == xAOD::UncalibMeasType::PixelClusterType and
    top_type == xAOD::UncalibMeasType::PixelClusterType) {
      return SeedType::PPP;
    } else if (bottom_type == xAOD::UncalibMeasType::PixelClusterType and
           middle_type == xAOD::UncalibMeasType::PixelClusterType and
           top_type == xAOD::UncalibMeasType::StripClusterType) {
      return SeedType::PPS;
    } else if (bottom_type == xAOD::UncalibMeasType::PixelClusterType and
           middle_type == xAOD::UncalibMeasType::StripClusterType and
           top_type == xAOD::UncalibMeasType::StripClusterType) {
      return SeedType::PSS;
    } else if (bottom_type == xAOD::UncalibMeasType::StripClusterType and
           middle_type == xAOD::UncalibMeasType::StripClusterType and
           top_type == xAOD::UncalibMeasType::StripClusterType) {
      return SeedType::SSS;
    } else {
      return SeedType::Others;
    }
  }

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

overridevirtual

Definition at line 80 of file ActsInspectTruthContentAlg.cxx.

                                                                              {
    ATH_MSG_DEBUG( "Executing " << name() << " ..." );
    
    // contextual stat collectors
    cluster_stat_t clusterStat {};
    seed_stat_t seedStat {};
    
    std::array<const ActsTrk::MeasurementToTruthParticleAssociation*, s_nClusterTypes> truths {};
    // Fill Stat info for Cluster collection(s)
    for (std::size_t i(0); i<m_clusters.size(); ++i) {
      ATH_MSG_DEBUG( "Retrieving cluster collection with key " << m_clusters.at(i).key() );
      SG::ReadHandle<xAOD::UncalibratedMeasurementContainer> clusterHandle = SG::makeHandle( m_clusters.at(i), ctx );
      ATH_CHECK( clusterHandle.isValid() );
      const xAOD::UncalibratedMeasurementContainer *clusters = clusterHandle.cptr();
 
      if (clusters->empty()) continue;
      xAOD::UncalibMeasType elementType = clusters->front()->type();
    
      ATH_MSG_DEBUG( "Retrieving Measurement to Truth Particle Association map with key: " << m_associationMap_key.at(i).key() );
      SG::ReadHandle<ActsTrk::MeasurementToTruthParticleAssociation> truthHandle = SG::makeHandle( m_associationMap_key.at(i), ctx );
      ATH_CHECK( truthHandle.isValid() );
      truths[static_cast<std::size_t>(elementType)] = truthHandle.cptr();
 
      // Check truth and clusters are compatible
      ATH_CHECK( truths[static_cast<std::size_t>(elementType)]->isCompatibleWith(clusters) );
 
      ATH_CHECK( fillStatClusters(*clusters,
                  *truths[static_cast<std::size_t>(elementType)],
                  clusterStat) );
    }
 
    ATH_CHECK( copyStatTable(clusterStat, m_clusterStat) );
    
    for (std::size_t i(0); i<m_seeds.size(); ++i) {
      ATH_MSG_DEBUG( "Retrieving seed collection with key: " << m_seeds.at(i).key() );
      SG::ReadHandle<ActsTrk::SeedContainer> seedHandle = SG::makeHandle( m_seeds.at(i), ctx );
      ATH_CHECK( seedHandle.isValid() );
      const ActsTrk::SeedContainer* seeds = seedHandle.cptr();
 
      ATH_CHECK( fillStatSeeds(*seeds,
                   truths,
                   seedStat) );
    }
 
    ATH_CHECK( copyStatTable(seedStat, m_seedStat) );
 
    for (std::size_t i(0); i<m_tracks.size(); ++i) {
      ATH_MSG_DEBUG( "Retrieving tracks collection with key " << m_tracks.at(i).key() );
      SG::ReadHandle<ActsTrk::TrackContainer> trackHandle = SG::makeHandle( m_tracks.at(i), ctx );
      ATH_CHECK( trackHandle.isValid() );
      const ActsTrk::TrackContainer* tracks = trackHandle.cptr();
 
      cluster_stat_t onTrack_clusterStat {};
      track_stat_t trackStat {};
      ATH_CHECK( fillStatTracks(*tracks,             
                truths,
                trackStat,
                onTrack_clusterStat) );
      
      ATH_CHECK( copyStatTable(trackStat, m_trackStat.at(i).second) );
      ATH_CHECK( copyStatTable(onTrack_clusterStat, m_onTrack_clusterStat.at(i).second) );
    }
 
    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();
}

fillStatClusters()

StatusCode ActsTrk::ActsInspectTruthContentAlg::fillStatClusters ( const xAOD::UncalibratedMeasurementContainer & container, const ActsTrk::MeasurementToTruthParticleAssociation & truth, typename ActsInspectTruthContentAlg::cluster_stat_t & stat ) const

private

Definition at line 147 of file ActsInspectTruthContentAlg.cxx.

  {
    ATH_MSG_DEBUG( "Checking truth for clusters ..." );
    for (const xAOD::UncalibratedMeasurement* meas : container) {
      std::size_t clusterTypeIndex = static_cast<std::size_t>(meas->type());
      ++stat[to_underlying(EStatClusters::kNTotal)][clusterTypeIndex];
      const auto& tps = truth.at(meas->index());
 
      // get number of (valid) contributions
      if (tps.empty()) {
        ++stat[to_underlying(EStatClusters::kNClustersWithNoBarcode)][clusterTypeIndex];
        continue;
      }
      
      // Check all barcodes are from primary particles
      bool allValidParticles = true;
      for (const auto* tp : tps) {
    if ( not HepMC::is_simulation_particle(*tp) ) continue;
    allValidParticles = false;
    break;
      }
      
      if (tps.size() == 1) {
    ++stat[to_underlying(EStatClusters::kNClustersWith1Contribution)][clusterTypeIndex];
    if (allValidParticles) ++stat[to_underlying(EStatClusters::kNClustersWith1ValidContribution)][clusterTypeIndex];
      }
      else if (tps.size() == 2) {
    ++stat[to_underlying(EStatClusters::kNClustersWith2Contribution)][clusterTypeIndex];
    if (allValidParticles) ++stat[to_underlying(EStatClusters::kNClustersWith2ValidContribution)][clusterTypeIndex];
      }
      else {
    ++stat[to_underlying(EStatClusters::kNClustersWith3Contribution)][clusterTypeIndex];
    if (allValidParticles) ++stat[to_underlying(EStatClusters::kNClustersWith3ValidContribution)][clusterTypeIndex];
      }
      
      // get main contribution
      bool hasContributionFromPrimaryParticle = false;
      for (const auto* tp : tps) {
    if ( HepMC::is_simulation_particle(*tp) ) continue;
    hasContributionFromPrimaryParticle = true;
    break;
      }
 
      if (not hasContributionFromPrimaryParticle) ++stat[to_underlying(EStatClusters::kNClustersWith200kBarcode)][clusterTypeIndex];
      else ++stat[to_underlying(EStatClusters::kNClustersFromPrimaries)][clusterTypeIndex];
    }
    
    return StatusCode::SUCCESS;
  }

fillStatSeeds()

StatusCode ActsTrk::ActsInspectTruthContentAlg::fillStatSeeds ( const ActsTrk::SeedContainer & seeds, std::array< const ActsTrk::MeasurementToTruthParticleAssociation *, s_nClusterTypes > & truths, seed_stat_t & stat ) const

private

Definition at line 200 of file ActsInspectTruthContentAlg.cxx.

  {
    ATH_MSG_DEBUG( "Checking truth for seeds ..." );
    for (std::size_t i(0); i<seeds.size(); ++i) {
      ActsTrk::Seed seed = seeds.at(i);
 
      int nMatches = 0;
      int nMeasurements = 0;
      std::unordered_map<std::size_t, int> particleIds {};
 
      std::size_t seedType = to_underlying(deduceSeedType(seed));
      ++stat[to_underlying(EStatSeeds::kNTotal)][seedType];
 
      const auto& sps = seed.sp();
      for ( const xAOD::SpacePoint* sp : sps ) {
    const auto& measurements = sp->measurements();
    for (const xAOD::UncalibratedMeasurement* meas : measurements ) {
      ++nMeasurements;
      
      std::size_t clusterTypeIndex = to_underlying(meas->type());
      const ActsTrk::MeasurementToTruthParticleAssociation* truth = truths.at(clusterTypeIndex);
      auto tps = truth->at(meas->index());
 
      if (tps.empty()) continue;
      
      bool contributionOnlyFromSimulationParticles = true;
      for (const auto* tp : tps) {
        if ( HepMC::is_simulation_particle(*tp) ) continue;
        contributionOnlyFromSimulationParticles = false;
 
        std::size_t pid = HepMC::uniqueID(tp);
        particleIds.try_emplace( pid, 0 );
        ++particleIds[pid];
      }
      if (contributionOnlyFromSimulationParticles) continue;
 
      ++nMatches;
    } // loop on measurements
      } // loop on space points
 
      // check if we have measurements associated to the same particle id
      // if the number of entries for the same id (the key) is the same as the
      // number of measurements in the seed, then we have a match
      bool isFromSameParticle = false;
      for (const auto [pid, nEntries] : particleIds) {
    if (nEntries != nMeasurements) continue;
    isFromSameParticle = true;
    break;
      }
 
      
      if (nMatches == 0) {
    ++stat[to_underlying(EStatSeeds::nKSeedsWith0Matches)][seedType];
      } else if (nMatches == 1) {
    ++stat[to_underlying(EStatSeeds::nKSeedsWith1Matches)][seedType];
      } else if (nMatches == 2) {
    ++stat[to_underlying(EStatSeeds::nKSeedsWith2Matches)][seedType];
    if (isFromSameParticle) ++stat[to_underlying(EStatSeeds::nKSeedsSame2Matches)][seedType];
      } else if (nMatches == 3) {
        ++stat[to_underlying(EStatSeeds::nKSeedsWith3Matches)][seedType];
    if (isFromSameParticle) ++stat[to_underlying(EStatSeeds::nKSeedsSame3Matches)][seedType];
      } else if (nMatches == 4) {
        ++stat[to_underlying(EStatSeeds::nKSeedsWith4Matches)][seedType];
    if (isFromSameParticle) ++stat[to_underlying(EStatSeeds::nKSeedsSame4Matches)][seedType];
      } else if (nMatches == 5) {
        ++stat[to_underlying(EStatSeeds::nKSeedsWith5Matches)][seedType];
    if (isFromSameParticle) ++stat[to_underlying(EStatSeeds::nKSeedsSame5Matches)][seedType];
      } else if (nMatches == 6) {
        ++stat[to_underlying(EStatSeeds::nKSeedsWith6Matches)][seedType];
    if (isFromSameParticle) ++stat[to_underlying(EStatSeeds::nKSeedsSame6Matches)][seedType];
      }                  
    } // loop on seed
    
    return StatusCode::SUCCESS;
  }

fillStatTracks()

StatusCode ActsTrk::ActsInspectTruthContentAlg::fillStatTracks ( const ActsTrk::TrackContainer & tracks, std::array< const ActsTrk::MeasurementToTruthParticleAssociation *, s_nClusterTypes > & truths, track_stat_t & trackStat, cluster_stat_t & onTrackClusterStat ) const

private

Definition at line 309 of file ActsInspectTruthContentAlg.cxx.

                                                               {
    ATH_MSG_DEBUG( "Checking truth for tracks ..." );
    for (const auto track : tracks) {
      ++trackStat[to_underlying(EStatTracks::kNTotal)][0];
      std::size_t nHoles = track.nHoles();
      if (nHoles == 0) ++trackStat[to_underlying(EStatTracks::kNTracks0Holes)][0];
      else if (nHoles == 1) ++trackStat[to_underlying(EStatTracks::kNTracks1Holes)][0];
      else if (nHoles == 2) ++trackStat[to_underlying(EStatTracks::kNTracks2Holes)][0];
      else ++trackStat[to_underlying(EStatTracks::kNTracks3Holes)][0];
 
      std::size_t nOutliers = track.nOutliers();
      if (nOutliers == 0) ++trackStat[to_underlying(EStatTracks::kNTracks0Outliers)][0];
      else if (nOutliers == 1) ++trackStat[to_underlying(EStatTracks::kNTracks1Outliers)][0];
      else if (nOutliers == 2) ++trackStat[to_underlying(EStatTracks::kNTracks2Outliers)][0];
      else ++trackStat[to_underlying(EStatTracks::kNTracks3Outliers)][0];
 
 
      bool AllValids = true;
      int nConsideredMeasurements = 0;
      std::unordered_map<std::size_t, int> particleIds {};
      
      // on track clusters
      track.container()
    .trackStateContainer().visitBackwards(track.tipIndex(),
                          [&truths, &onTrackStat,
                           &AllValids,
                           &nConsideredMeasurements, &particleIds]
                          (const auto& state) {
                        auto flags = state.typeFlags();
                        if (not flags.test(Acts::TrackStateFlag::MeasurementFlag) and
                            not flags.test(Acts::TrackStateFlag::OutlierFlag)) return;
                        ++nConsideredMeasurements;
                        
                        // get cluster
                        auto sl = state.getUncalibratedSourceLink().template get<ATLASUncalibSourceLink>();
                        assert( sl != nullptr);
                        const xAOD::UncalibratedMeasurement &meas = getUncalibratedMeasurement(sl);
                        
                        std::size_t clusterTypeIndex = to_underlying(meas.type());
                        ++onTrackStat[to_underlying(EStatClusters::kNTotal)][clusterTypeIndex];
                        
                        const ActsTrk::MeasurementToTruthParticleAssociation* truth = truths[clusterTypeIndex];
                        const auto& tps = truth->at(meas.index());
 
                        if (tps.empty()) {
                          ++onTrackStat[to_underlying(EStatClusters::kNClustersWithNoBarcode)][clusterTypeIndex];
                          AllValids = false;
                          return;
                        }
                        
                        bool allBarcodesValid = true;
                        bool hasContributionFromPrimaryParticle = false;
                        for (const auto* tp : tps) {
                          if ( not HepMC::is_simulation_particle(*tp) ) {
                            hasContributionFromPrimaryParticle = true;
                            continue;
                          }
                          allBarcodesValid = false;
                          break;
                        }
                        if (not hasContributionFromPrimaryParticle) AllValids = false;
                        
                        if (tps.size() == 1) {
                          ++onTrackStat[to_underlying(EStatClusters::kNClustersWith1Contribution)][clusterTypeIndex];
                          if (allBarcodesValid) ++onTrackStat[to_underlying(EStatClusters::kNClustersWith1ValidContribution)][clusterTypeIndex];
                        }
                        else if (tps.size() == 2) {
                          ++onTrackStat[to_underlying(EStatClusters::kNClustersWith2Contribution)][clusterTypeIndex];
                          if (allBarcodesValid) ++onTrackStat[to_underlying(EStatClusters::kNClustersWith2ValidContribution)][clusterTypeIndex];
                        }
                        else {
                          ++onTrackStat[to_underlying(EStatClusters::kNClustersWith3Contribution)][clusterTypeIndex];
                          if (allBarcodesValid) ++onTrackStat[to_underlying(EStatClusters::kNClustersWith3ValidContribution)][clusterTypeIndex];
                        }
 
                        bool contributionOnlyFromSimulationParticles = true;
                        for (const auto* tp : tps) {
                          if ( HepMC::is_simulation_particle(*tp) ) continue;
                          contributionOnlyFromSimulationParticles = false;
                          
                          std::size_t pid = HepMC::uniqueID(tp);
                          particleIds.try_emplace( pid, 0 );
                          ++particleIds[pid];
                        }
                        
                        // get main contribution
                        if (contributionOnlyFromSimulationParticles) ++onTrackStat[to_underlying(EStatClusters::kNClustersWith200kBarcode)][clusterTypeIndex];
                        else ++onTrackStat[to_underlying(EStatClusters::kNClustersFromPrimaries)][clusterTypeIndex];
                          }); // llop on states
 
      bool AllSameBarcode = false;
      for ( const auto [pid, nEntries] : particleIds ) {
    if (nEntries != nConsideredMeasurements) continue;
    AllSameBarcode = true;
    break;
      }
      
      if (AllValids) ++trackStat[to_underlying(EStatTracks::kNFullMatch)][0];
      if (AllSameBarcode) ++trackStat[to_underlying(EStatTracks::kNPerfectMatch)][0];
    } // loop on tracks    
 
    return StatusCode::SUCCESS;
  }

filterPassed()

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

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

finalize()

StatusCode ActsTrk::ActsInspectTruthContentAlg::finalize ( )

overridevirtual

Definition at line 46 of file ActsInspectTruthContentAlg.cxx.

                                                  {
    ATH_MSG_INFO( "Finalizing " << name() << " ..." );
    ATH_MSG_INFO( "Statistics from Seed check with truth info:" );
 
    if ( printStatTables<ActsInspectTruthContentAlg::EStatClusters, xAOD::UncalibMeasType>("Clusters", m_clusterStat).isFailure() ) {
      ATH_MSG_FATAL("Problem dumping Cluster truth information");
      return StatusCode::FAILURE;
    }
 
    if ( printStatTables<ActsInspectTruthContentAlg::EStatSeeds, ActsInspectTruthContentAlg::SeedType>("Seeds", m_seedStat).isFailure() ) {
      ATH_MSG_FATAL("Problem dumping Seed truth information");
      return StatusCode::FAILURE;
    }
 
    for (const auto& [trackCollectionName, onTrackStat] : m_onTrack_clusterStat) {
      std::string reportName = "On Track Clusters (" + trackCollectionName + ")";
      if ( printStatTables<ActsInspectTruthContentAlg::EStatClusters, xAOD::UncalibMeasType>(reportName, onTrackStat).isFailure() ) {
    ATH_MSG_FATAL("Problem dumping On Track Cluster truth info (" << trackCollectionName << ")");
    return StatusCode::FAILURE;
      }
    }
    
    for (const auto& [trackCollectionName, trackStat] : m_trackStat) {
      std::string reportName = "Track (" + trackCollectionName + ")";
      if ( printStatTables<ActsInspectTruthContentAlg::EStatTracks, ActsInspectTruthContentAlg::TrackType>(reportName, trackStat).isFailure() ) {
        ATH_MSG_FATAL("Problem dumping Track truth info (" << trackCollectionName << ")");
        return StatusCode::FAILURE;
      }
    }
    
    return StatusCode::SUCCESS;
  }

initialize()

StatusCode ActsTrk::ActsInspectTruthContentAlg::initialize ( )

overridevirtual

Definition at line 23 of file ActsInspectTruthContentAlg.cxx.

                                                    {
    ATH_MSG_INFO( "Initializing " << name() << " ..." );
 
    ATH_CHECK( m_clusters.initialize() );
    ATH_CHECK( m_seeds.initialize() );
    ATH_CHECK( m_associationMap_key.initialize() );
    ATH_CHECK( m_tracks.initialize() );
    
    if (m_clusters.size() != m_associationMap_key.size()) {
      ATH_MSG_ERROR("Inconsistent sizes of Clusters and TruthAssociationMaps");
      return StatusCode::FAILURE;
    }
 
    for (const auto& trackKey : m_tracks) {
      std::string containerName = trackKey.key();
      m_onTrack_clusterStat.push_back( std::make_pair(containerName, cluster_stat_t()) );
      m_trackStat.push_back( std::make_pair(containerName, track_stat_t()) );
    }
 
    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()

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

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.

printStatTables()

template<typename row_t, typename coll_t, typename stat_t>

StatusCode ActsTrk::ActsInspectTruthContentAlg::printStatTables ( const std::string & objectCollectionName, const stat_t & stat ) const

private

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.

to_label() [1/3]

std::string ActsTrk::ActsInspectTruthContentAlg::to_label ( ActsInspectTruthContentAlg::EStatClusters type ) const

inlineprivate

to_label() [2/3]

std::string ActsTrk::ActsInspectTruthContentAlg::to_label ( ActsInspectTruthContentAlg::EStatSeeds type ) const

inlineprivate

to_label() [3/3]

std::string ActsTrk::ActsInspectTruthContentAlg::to_label ( ActsInspectTruthContentAlg::EStatTracks type ) const

inlineprivate

to_string() [1/3]

std::string ActsTrk::ActsInspectTruthContentAlg::to_string ( ActsInspectTruthContentAlg::SeedType type ) const

inlineprivate

to_string() [2/3]

std::string ActsTrk::ActsInspectTruthContentAlg::to_string ( ActsInspectTruthContentAlg::TrackType type ) const

inlineprivate

to_string() [3/3]

std::string ActsTrk::ActsInspectTruthContentAlg::to_string ( xAOD::UncalibMeasType type ) const

inlineprivate

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 [1/5]

std::mutex m_mutex ActsTrk::ActsInspectTruthContentAlg::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 130 of file ActsInspectTruthContentAlg.h.

{};

ATLAS_THREAD_SAFE [2/5]

cluster_stat_t m_clusterStat ActsTrk::ActsInspectTruthContentAlg::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 131 of file ActsInspectTruthContentAlg.h.

{};

ATLAS_THREAD_SAFE [3/5]

seed_stat_t m_seedStat ActsTrk::ActsInspectTruthContentAlg::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 132 of file ActsInspectTruthContentAlg.h.

{};

ATLAS_THREAD_SAFE [4/5]

std::vector<std::pair<std::string, track_stat_t> > m_trackStat ActsTrk::ActsInspectTruthContentAlg::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 133 of file ActsInspectTruthContentAlg.h.

{};

ATLAS_THREAD_SAFE [5/5]

std::vector<std::pair<std::string, cluster_stat_t> > m_onTrack_clusterStat ActsTrk::ActsInspectTruthContentAlg::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 134 of file ActsInspectTruthContentAlg.h.

{};

m_associationMap_key

SG::ReadHandleKeyArray<ActsTrk::MeasurementToTruthParticleAssociation> ActsTrk::ActsInspectTruthContentAlg::m_associationMap_key {this,"TruthAssociationMaps", {}}

private

Definition at line 126 of file ActsInspectTruthContentAlg.h.

{this,"TruthAssociationMaps", {}};

m_clusters

SG::ReadHandleKeyArray<xAOD::UncalibratedMeasurementContainer> ActsTrk::ActsInspectTruthContentAlg::m_clusters {this, "Clusters", {}}

private

Definition at line 124 of file ActsInspectTruthContentAlg.h.

{this, "Clusters", {}};

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 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_seeds

SG::ReadHandleKeyArray<ActsTrk::SeedContainer> ActsTrk::ActsInspectTruthContentAlg::m_seeds {this, "Seeds", {}}

private

Definition at line 125 of file ActsInspectTruthContentAlg.h.

{this, "Seeds", {}};

m_tracks

SG::ReadHandleKeyArray<ActsTrk::TrackContainer> ActsTrk::ActsInspectTruthContentAlg::m_tracks {this, "Tracks", {}}

private

Definition at line 127 of file ActsInspectTruthContentAlg.h.

{this, "Tracks", {}};

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.

s_nClusterTypes

std::size_t ActsTrk::ActsInspectTruthContentAlg::s_nClusterTypes = static_cast<std::size_t>(xAOD::UncalibMeasType::nTypes)

staticconstexprprivate

Definition at line 90 of file ActsInspectTruthContentAlg.h.

s_nSeedTypes

std::size_t ActsTrk::ActsInspectTruthContentAlg::s_nSeedTypes = static_cast<std::size_t>(SeedType::nTypes)

staticconstexprprivate

Definition at line 91 of file ActsInspectTruthContentAlg.h.

---

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

• ActsInspectTruthContentAlg.h
• ActsInspectTruthContentAlg.cxx