ThinInDetClustersAlg Class Reference

#include <ThinInDetClustersAlg.h>

Inheritance diagram for ThinInDetClustersAlg:

Collaboration diagram for ThinInDetClustersAlg:

Public Member Functions
	ThinInDetClustersAlg (const std::string &name, ISvcLocator *pSvcLocator)
	Constructor with parameters:
virtual	~ThinInDetClustersAlg ()=default
	Destructor:
virtual StatusCode	initialize () override
	Athena algorithm's initalize hook.
virtual StatusCode	execute (const EventContext &ctx) const override
	Athena algorithm's execute hook.
virtual StatusCode	finalize () override
	Athena algorithm's finalize hook.
StatusCode	initializeParser (const ExpressionParsing::SelectionArg< 1 > &selection_string)
StatusCode	finalizeParser ()
virtual bool	updateDataNeeds (const std::vector< const DataObjID * > &input_data_in, const std::vector< const DataObjID * > &output_data_in, std::vector< Gaudi::DataHandle * > &new_input_handles, std::vector< Gaudi::DataHandle * > &new_output_handles) override
	Inform a dynamic data consumer about all consumed or produced data.
virtual bool	renounceInput (const DataObjID &output_data_id) 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
StatusCode	_initializeParser (const ExpressionParsing::SelectionArg< 1 > &selection_string, ExpressionParsing::AppendLoaderFunc pre_xaod_loader=ExpressionParsing::NoLoaderFunc, ExpressionParsing::AppendLoaderFunc post_xaod_loader=ExpressionParsing::NoLoaderFunc)
StatusCode	_initializeParser (const ExpressionParsing::SelectionArg< NUM_PARSER > &selection_string, ExpressionParsing::AppendLoaderFunc pre_xaod_loader=ExpressionParsing::NoLoaderFunc, ExpressionParsing::AppendLoaderFunc post_xaod_loader=ExpressionParsing::NoLoaderFunc)
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.

Protected Attributes
ExpressionParsing::SGxAODProxyLoader *	m_xAODProxyLoader
	The expression parser.
std::unique_ptr< ExpressionParsing::IProxyLoader >	m_proxyLoaders
std::conditional< 1==1, std::unique_ptr< ExpressionParsing::ExpressionParser >, std::array< std::unique_ptr< ExpressionParsing::ExpressionParser >, 1 > >::type	m_parser
std::vector< std::string >	m_renounce

Private Types
enum	MeasurementType {   TrkState_unidentified = 0 , TrkState_Pixel = 1 , TrkState_SCT = 2 , TrkState_TRT = 3 ,   TrkState_Pseudo = 8 , TrkState_Vertex = 9 , TrkState_SpacePoint = 11 , TrkState_NumberOfMeasurementTypes =16 }
	Track state types Subset of enum MeasurementType from Athena: Tracking/TrkEvent/TrkEventPrimitives/TrkEventPrimitives/TrackStateDefs.h. More...
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	selectTrackHits (const xAOD::TrackParticleContainer &inputTrackParticles, const std::vector< bool > &inputMask, MeasurementType detTypeToSelect, std::vector< bool > &outputStatesMask, std::vector< bool > &outputMeasurementsMask) const
	Select TrackStateValidation and TrackMeasurementValidation objects that are used in the (thinned) track container.
StatusCode	filterTrackHits (MeasurementType detTypeToSelect, const xAOD::TrackParticleContainer &inputTrackParticles, const std::vector< bool > &inputMask, const SG::ThinningHandleKey< xAOD::TrackStateValidationContainer > &statesKey, const SG::ThinningHandleKey< xAOD::TrackMeasurementValidationContainer > &measurementsKey, std::atomic< unsigned int > &ntot_states, std::atomic< unsigned int > &npass_states, std::atomic< unsigned int > &npass_measurements, const EventContext &ctx) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
std::atomic< unsigned int >	m_ntot {}
	Counters and keys for xAOD::TrackParticle container.
std::atomic< unsigned int >	m_npass {}
BooleanProperty	m_thinPixelHitsOnTrack { this, "ThinPixelHitsOnTrack", false, ""}
	Thinning logic.
BooleanProperty	m_thinSCTHitsOnTrack { this, "ThinSCTHitsOnTrack", false, ""}
BooleanProperty	m_thinTRTHitsOnTrack { this, "ThinTRTHitsOnTrack", false, ""}
StringProperty	m_streamName { this, "StreamName", "", "Name of the stream being thinned" }
	Stream for object thinning selections.
StringProperty	m_selectionString { this, "SelectionString", "InDetTrackParticles.pt>10*GeV", "Selection string for each TrackParticle container" }
	Expressions for object thinning selections Default InDetTrackParticles (also accomplished through flag.Tracking.thinPixelClustersSelectionStrings defaults).
SG::ThinningHandleKey< xAOD::TrackParticleContainer >	m_inDetSGKey { this, "InDetTrackParticlesKey", "InDetTrackParticles", "" }
	SGKey for TrackParticleContainer.
std::atomic< unsigned int >	m_ntot_pix_states {}
	Counters and keys for xAOD::TrackStateValidation and xAOD::TrackMeasurementValidation containers.
std::atomic< unsigned int >	m_npass_pix_states {}
SG::ThinningHandleKey< xAOD::TrackStateValidationContainer >	m_statesPixSGKey { this, "InDetTrackStatesPixKey", "PixelMSOSs", "" }
std::atomic< unsigned int >	m_ntot_pix_measurements {}
std::atomic< unsigned int >	m_npass_pix_measurements {}
SG::ThinningHandleKey< xAOD::TrackMeasurementValidationContainer >	m_measurementsPixSGKey { this, "InDetTrackMeasurementsPixKey", "PixelClusters", "" }
std::atomic< unsigned int >	m_ntot_sct_states {}
std::atomic< unsigned int >	m_npass_sct_states {}
SG::ThinningHandleKey< xAOD::TrackStateValidationContainer >	m_statesSctSGKey { this, "InDetTrackStatesSctKey", "SCT_MSOSs", "" }
std::atomic< unsigned int >	m_ntot_sct_measurements {}
std::atomic< unsigned int >	m_npass_sct_measurements {}
SG::ThinningHandleKey< xAOD::TrackMeasurementValidationContainer >	m_measurementsSctSGKey { this, "InDetTrackMeasurementsSctKey", "SCT_Clusters", "" }
std::atomic< unsigned int >	m_ntot_trt_states {}
std::atomic< unsigned int >	m_npass_trt_states {}
SG::ThinningHandleKey< xAOD::TrackStateValidationContainer >	m_statesTrtSGKey { this, "InDetTrackStatesTrtKey", "TRT_MSOSs", "" }
std::atomic< unsigned int >	m_ntot_trt_measurements {}
std::atomic< unsigned int >	m_npass_trt_measurements {}
SG::ThinningHandleKey< xAOD::TrackMeasurementValidationContainer >	m_measurementsTrtSGKey { this, "InDetTrackMeasurementsTrtKey", "TRT_DriftCircles", "" }
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollection >	m_SCTDetEleCollKey {this, "SCTDetEleCollKey", "SCT_DetectorElementCollection", "Key of SiDetectorElementCollection for SCT"}
	For P->T converter of SCT_Clusters.
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

Detailed Description

Definition at line 32 of file ThinInDetClustersAlg.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

MeasurementType

enum ThinInDetClustersAlg::MeasurementType

private

Track state types Subset of enum MeasurementType from Athena: Tracking/TrkEvent/TrkEventPrimitives/TrkEventPrimitives/TrackStateDefs.h.

Enumerator
TrkState_unidentified
TrkState_Pixel
TrkState_SCT
TrkState_TRT
TrkState_Pseudo
TrkState_Vertex
TrkState_SpacePoint
TrkState_NumberOfMeasurementTypes

Definition at line 112 of file ThinInDetClustersAlg.h.

                       {
    TrkState_unidentified = 0,
    TrkState_Pixel      = 1,
    TrkState_SCT        = 2,
    TrkState_TRT        = 3,
    TrkState_Pseudo     = 8,
    TrkState_Vertex     = 9,
    TrkState_SpacePoint = 11,
    TrkState_NumberOfMeasurementTypes=16
  };

Constructor & Destructor Documentation

ThinInDetClustersAlg()

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

Constructor with parameters:

Definition at line 32 of file ThinInDetClustersAlg.cxx.

                           : 
 
  ExpressionParserUser< ::AthReentrantAlgorithm>( name, pSvcLocator )
{
}

~ThinInDetClustersAlg()

virtual ThinInDetClustersAlg::~ThinInDetClustersAlg ( )

virtualdefault

Destructor:

Member Function Documentation

_initializeParser() [1/2]

StatusCode ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::_initializeParser ( const ExpressionParsing::SelectionArg< 1 > & selection_string, ExpressionParsing::AppendLoaderFunc pre_xaod_loader = ExpressionParsing::NoLoaderFunc, ExpressionParsing::AppendLoaderFunc post_xaod_loader = ExpressionParsing::NoLoaderFunc )

protectedinherited

_initializeParser() [2/2]

StatusCode ExpressionParserUserBase< ::AthReentrantAlgorithm, NUM_PARSER >::_initializeParser ( const ExpressionParsing::SelectionArg< NUM_PARSER > & selection_string, ExpressionParsing::AppendLoaderFunc pre_xaod_loader = ExpressionParsing::NoLoaderFunc, ExpressionParsing::AppendLoaderFunc post_xaod_loader = ExpressionParsing::NoLoaderFunc )

protectedinherited

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

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

overridevirtual

Athena algorithm's execute hook.

Definition at line 116 of file ThinInDetClustersAlg.cxx.

{
 
  // Retrieve main TrackParticle collection
  SG::ThinningHandle<xAOD::TrackParticleContainer> importedTrackParticles(m_inDetSGKey, ctx);
  ATH_CHECK( importedTrackParticles.isValid() );

  if (m_thinPixelHitsOnTrack && !m_measurementsPixSGKey.empty()) {
    SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> importedMeasurements(m_measurementsPixSGKey, ctx);
    ATH_CHECK( importedMeasurements.isValid() );
    unsigned int size_measurements = importedMeasurements->size();
    m_ntot_pix_measurements += size_measurements;
  }
  if (m_thinSCTHitsOnTrack && !m_measurementsSctSGKey.empty()) {
    SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> importedMeasurements(m_measurementsSctSGKey, ctx);
    ATH_CHECK( importedMeasurements.isValid() );
    unsigned int size_measurements = importedMeasurements->size();
    m_ntot_sct_measurements += size_measurements;
  }
  if (m_thinTRTHitsOnTrack && !m_measurementsTrtSGKey.empty()) {
    SG::ReadHandle<xAOD::TrackMeasurementValidationContainer> importedMeasurements(m_measurementsTrtSGKey, ctx);
    ATH_CHECK( importedMeasurements.isValid() );
    unsigned int size_measurements = importedMeasurements->size();
    m_ntot_trt_measurements += size_measurements;
  }
 
  // Check the event contains tracks
  unsigned int nTracks = importedTrackParticles->size();
  if (nTracks==0) return StatusCode::SUCCESS;
 
  // Set up a mask with the same entries as the full TrackParticle collection
  std::vector<bool> mask;
  mask.assign(nTracks,false); // default: don't keep any tracks
  m_ntot += nTracks;
 
  // Execute the text parser and update the mask
  if (m_parser) {
    std::vector<int> entries =  m_parser->evaluateAsVector();
    unsigned int nEntries = entries.size();
    // check the sizes are compatible
    if (nTracks != nEntries ) {
      ATH_MSG_ERROR("Sizes incompatible! Are you sure your selection string used ID TrackParticles?");
      return StatusCode::FAILURE;
    }
    // set mask
    for (unsigned int i=0; i<nTracks; ++i) {
      if (entries[i]==1) {
        mask[i]=true; // Set track mask
        m_npass += 1; // Count passing tracks.
      }
    }
  }
 
  if (m_thinPixelHitsOnTrack) {
    ATH_CHECK( filterTrackHits (
                                TrkState_Pixel,
                                *importedTrackParticles,
                                mask,
                                m_statesPixSGKey,
                                m_measurementsPixSGKey,
                                m_ntot_pix_states,
                                m_npass_pix_states,
                                m_npass_pix_measurements,
                                ctx) );
  }
  if (m_thinSCTHitsOnTrack) {
    ATH_CHECK( filterTrackHits (
                                TrkState_SCT,
                                *importedTrackParticles,
                                mask,
                                m_statesSctSGKey,
                                m_measurementsSctSGKey,
                                m_ntot_sct_states,
                                m_npass_sct_states,
                                m_npass_sct_measurements,
                                ctx) );
  }
  if (m_thinTRTHitsOnTrack) {
    ATH_CHECK( filterTrackHits (
                                TrkState_TRT,
                                *importedTrackParticles,
                                mask,
                                m_statesTrtSGKey,
                                m_measurementsTrtSGKey,
                                m_ntot_trt_states,
                                m_npass_trt_states,
                                m_npass_trt_measurements,
                                ctx) );
  }
  
  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();
  }

filterTrackHits()

StatusCode ThinInDetClustersAlg::filterTrackHits ( MeasurementType detTypeToSelect, const xAOD::TrackParticleContainer & inputTrackParticles, const std::vector< bool > & inputMask, const SG::ThinningHandleKey< xAOD::TrackStateValidationContainer > & statesKey, const SG::ThinningHandleKey< xAOD::TrackMeasurementValidationContainer > & measurementsKey, std::atomic< unsigned int > & ntot_states, std::atomic< unsigned int > & npass_states, std::atomic< unsigned int > & npass_measurements, const EventContext & ctx ) const

private

Definition at line 211 of file ThinInDetClustersAlg.cxx.

{
  std::vector<bool> maskStates;
  std::vector<bool> maskMeasurements;
 
  selectTrackHits (inputTrackParticles, inputMask, detTypeToSelect,
                   maskStates, maskMeasurements);
 
  auto count = [] (const std::vector<bool>& m)
    { return std::count (m.begin(), m.end(), true); };
  npass_states += count (maskStates);
  npass_measurements += count (maskMeasurements);
 
  if (!statesKey.empty()) {
    SG::ThinningHandle<xAOD::TrackStateValidationContainer> importedStates(statesKey, ctx);
    ATH_CHECK( importedStates.isValid() );
    unsigned int size_states = importedStates->size();
    if (size_states == 0) {
      ATH_MSG_WARNING("States container is empty: " << statesKey.key());
    }
    else {
      ntot_states += size_states;
      if (maskStates.size() > size_states) {
        ATH_MSG_ERROR("States mask size mismatch: mask size (" << maskStates.size() << ") > number of states (" << size_states << ").");
        return StatusCode::FAILURE;
      }
      maskStates.resize (size_states);
      importedStates.keep (maskStates);
    }
  }
 
  if (!measurementsKey.empty()) {
    SG::ThinningHandle<xAOD::TrackMeasurementValidationContainer> importedMeasurements(measurementsKey, ctx);
    ATH_CHECK( importedMeasurements.isValid() );
    unsigned int size_measurements = importedMeasurements->size();
    if (size_measurements == 0) {
      ATH_MSG_WARNING("Measurements container is empty: " << measurementsKey.key());
    }
    else {
      if (maskMeasurements.size() > size_measurements) {
        ATH_MSG_ERROR("Measurements mask size mismatch: mask size (" << maskMeasurements.size() << ") > number of measurements (" << size_measurements << ").");
        return StatusCode::FAILURE;
      }
      maskMeasurements.resize (size_measurements);
      importedMeasurements.keep (maskMeasurements);
    }
  }
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode ThinInDetClustersAlg::finalize ( )

overridevirtual

Athena algorithm's finalize hook.

Definition at line 84 of file ThinInDetClustersAlg.cxx.

{
  ATH_MSG_VERBOSE("finalize() ...");
  ATH_MSG_INFO("Processed "<< m_ntot <<" tracks, "<< m_npass<< " had their hits retained ");
  if (m_thinPixelHitsOnTrack) {
    ATH_MSG_INFO("Pixel state objects thinning, Total / Passed (Efficiency): "
                 << m_ntot_pix_states << " / " << m_npass_pix_states
                 << " (" << (m_ntot_pix_states == 0 ? 0 : static_cast<float>(m_npass_pix_states) / m_ntot_pix_states) << ")");
    ATH_MSG_INFO("Total pixel measurements objects from *all* track collections:: " << m_ntot_pix_measurements);
    ATH_MSG_INFO("Pixel measurements objects from this track collection passing selection: " << m_npass_pix_measurements);
  }
  if (m_thinSCTHitsOnTrack) {
    ATH_MSG_INFO("SCT state objects thinning, Total / Passed (Efficiency): "
                 << m_ntot_sct_states << " / " << m_npass_sct_states
                 << " (" << (m_ntot_sct_states == 0 ? 0 : static_cast<float>(m_npass_sct_states) / m_ntot_sct_states) << ")");
    ATH_MSG_INFO("Total SCT measurements objects from *all* track collections:: " << m_ntot_sct_measurements);
    ATH_MSG_INFO("SCT measurements objects from this track collection passing selection: " << m_npass_sct_measurements);
  }
  if (m_thinTRTHitsOnTrack) {
    ATH_MSG_INFO("TRT state objects thinning, Total / Passed (Efficiency): "
                 << m_ntot_trt_states << " / " << m_npass_trt_states
                 << " (" << (m_ntot_trt_states == 0 ? 0 : static_cast<float>(m_npass_trt_states) / m_ntot_trt_states) << ")");
    ATH_MSG_INFO("Total TRT measurements objects from *all* track collections:: " << m_ntot_trt_measurements);
    ATH_MSG_INFO("TRT measurements objects from this track collection passing selection: " << m_npass_trt_measurements);
  }
 
  ATH_CHECK( finalizeParser() );
 
  return StatusCode::SUCCESS;
}

finalizeParser()

StatusCode ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::finalizeParser ( )

inherited

initialize()

StatusCode ThinInDetClustersAlg::initialize ( )

overridevirtual

Athena algorithm's initalize hook.

Definition at line 42 of file ThinInDetClustersAlg.cxx.

{
 
  // Set up the text-parsing machinery for thinning the tracks directly according to user cuts
  if (m_selectionString.empty()) {
    ATH_MSG_FATAL("No inner detector track selection strings provided!");
    return StatusCode::FAILURE;
  }
  ATH_MSG_INFO("Track thinning selection string: " << m_selectionString);
  ATH_CHECK(initializeParser(m_selectionString));
 
  ATH_CHECK( m_inDetSGKey.initialize (m_streamName) );
  ATH_MSG_INFO("Using " << m_inDetSGKey << "as the source collection for inner detector track particles");
 
  if (m_thinPixelHitsOnTrack) {
    ATH_MSG_INFO("Pixel states collection as source for thinning: " << m_statesPixSGKey.key());
    ATH_MSG_INFO("Pixel measurements collection as source for thinning: " << m_measurementsPixSGKey.key());
  }
  if (m_thinSCTHitsOnTrack) {
    ATH_MSG_INFO("SCT states collection as source for thinning: " << m_statesSctSGKey.key());
    ATH_MSG_INFO("SCT measurements collection as source for thinning: " << m_measurementsSctSGKey.key());
  }
  if (m_thinTRTHitsOnTrack) {
    ATH_MSG_INFO("TRT states collection as source for thinning: " << m_statesTrtSGKey.key());
    ATH_MSG_INFO("TRT measurements collection as source for thinning: " << m_measurementsTrtSGKey.key());
  }
 
  ATH_MSG_DEBUG("Initializing measurement ThinningHandleKeys.");
  ATH_CHECK( m_measurementsPixSGKey.initialize (m_streamName, m_thinPixelHitsOnTrack && !m_measurementsPixSGKey.empty()) );
  ATH_CHECK( m_measurementsSctSGKey.initialize (m_streamName, m_thinSCTHitsOnTrack && !m_measurementsSctSGKey.empty()) );
  ATH_CHECK( m_measurementsTrtSGKey.initialize (m_streamName, m_thinTRTHitsOnTrack && !m_measurementsTrtSGKey.empty()) );
 
  ATH_MSG_DEBUG("Initializing MSOS ThinningHandleKeyArrays.");
  ATH_CHECK( m_statesPixSGKey.initialize (m_streamName, m_thinPixelHitsOnTrack && !m_statesPixSGKey.empty()) );
  ATH_CHECK( m_statesSctSGKey.initialize (m_streamName, m_thinSCTHitsOnTrack && !m_statesSctSGKey.empty()) );
  ATH_CHECK( m_statesTrtSGKey.initialize (m_streamName, m_thinTRTHitsOnTrack && !m_statesTrtSGKey.empty()) );
 
  ATH_CHECK(m_SCTDetEleCollKey.initialize( !m_SCTDetEleCollKey.key().empty() ));
 
  return StatusCode::SUCCESS;
}

initializeParser()

StatusCode ExpressionParserUser< ::AthReentrantAlgorithm, 1 >::initializeParser ( const ExpressionParsing::SelectionArg< 1 > & selection_string )

inherited

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 HiveAlgBase, 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.

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

renounceInput()

virtual bool ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::renounceInput ( const DataObjID & output_data_id )

overridevirtualinherited

selectTrackHits()

void ThinInDetClustersAlg::selectTrackHits ( const xAOD::TrackParticleContainer & inputTrackParticles, const std::vector< bool > & inputMask, MeasurementType detTypeToSelect, std::vector< bool > & outputStatesMask, std::vector< bool > & outputMeasurementsMask ) const

private

Select TrackStateValidation and TrackMeasurementValidation objects that are used in the (thinned) track container.

Definition at line 271 of file ThinInDetClustersAlg.cxx.

{
  using StatesOnTrack = std::vector<ElementLink<xAOD::TrackStateValidationContainer>>;
 
  // loop over track particles, consider only the ones pre-selected by the mask
  int trkIndex=-1;
  for (const xAOD::TrackParticle* trkIt : inputTrackParticles) {
    trkIndex++;
 
    if (not inputMask[trkIndex]) continue;
 
    // loop over the TrackStateValidation objects, and add them to the outputStatesMask
    static const SG::AuxElement::ConstAccessor< StatesOnTrack > trackStateAcc("Reco_msosLink");
    if( ! trackStateAcc.isAvailable( *trkIt ) ) {
      ATH_MSG_INFO("Cannot find TrackState link from xAOD::TrackParticle. Skipping track.");
      continue;
    }
    const StatesOnTrack& measurementsOnTrack = trackStateAcc(*trkIt);
 
    for( const ElementLink<xAOD::TrackStateValidationContainer>& trkState_el : measurementsOnTrack) {
      if (not trkState_el.isValid()) {
        ATH_MSG_INFO("Cannot find a valid link to TrackStateValidation object for track index: " << trkIndex);
        continue; //not a valid link
      }
      if ((*trkState_el)->detType() != detTypeToSelect) {
        ATH_MSG_VERBOSE("Discarding TrackState as not of correct type " << detTypeToSelect);
        continue;
      }
      if (trkState_el.index() >= outputStatesMask.size()) {
        outputStatesMask.resize (trkState_el.index()+1);
      }
      outputStatesMask[trkState_el.index()] = true;
 
      // get the corresponding TrackMeasurementValidation object, if any, and add it to the outputMeasurementsMask
      const ElementLink<xAOD::TrackMeasurementValidationContainer> trkMeasurement_el = (*trkState_el)->trackMeasurementValidationLink();
      if (not trkMeasurement_el.isValid()) {
        ATH_MSG_VERBOSE("Cannot find a valid link to TrackMeasurementValidation object from track state for track index: " << trkIndex
                        << ", trackState index: " << trkState_el.index());
        continue; //not a valid link
      }
      if (*trkMeasurement_el == nullptr) {
        ATH_MSG_VERBOSE("Invalid pointer to TrackMeasurementValidation object from track state for track index: " << trkIndex
                        << ", trackState index: " << trkState_el.index());
        continue; //not linking to a valid object -- is it necessary?
      }
      if (trkMeasurement_el.index() >= outputMeasurementsMask.size()) {
        outputMeasurementsMask.resize (trkMeasurement_el.index()+1);
      }
      outputMeasurementsMask[trkMeasurement_el.index()] = true;
    }
  } // end loop over xAOD::TrackParticle container
}

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.

updateDataNeeds()

virtual bool ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::updateDataNeeds ( const std::vector< const DataObjID * > & input_data_in, const std::vector< const DataObjID * > & output_data_in, std::vector< Gaudi::DataHandle * > & new_input_handles, std::vector< Gaudi::DataHandle * > & new_output_handles )

overridevirtualinherited

Inform a dynamic data consumer about all consumed or produced data.

Parameters

input_data_in	the input data object IDs of all algorithms and tools
output_data_in	the output data object IDs of all algorithms and tools
new_input_handles	the new input data handles added by this call
new_output_handles	the new output data handles added by this call

Returns

true if the consumer did update its handles, false otherwise

updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

m_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_inDetSGKey

SG::ThinningHandleKey<xAOD::TrackParticleContainer> ThinInDetClustersAlg::m_inDetSGKey { this, "InDetTrackParticlesKey", "InDetTrackParticles", "" }

private

SGKey for TrackParticleContainer.

Definition at line 73 of file ThinInDetClustersAlg.h.

{ this, "InDetTrackParticlesKey", "InDetTrackParticles", "" };

m_measurementsPixSGKey

SG::ThinningHandleKey<xAOD::TrackMeasurementValidationContainer> ThinInDetClustersAlg::m_measurementsPixSGKey { this, "InDetTrackMeasurementsPixKey", "PixelClusters", "" }

private

Definition at line 84 of file ThinInDetClustersAlg.h.

{ this, "InDetTrackMeasurementsPixKey", "PixelClusters", "" };

m_measurementsSctSGKey

SG::ThinningHandleKey<xAOD::TrackMeasurementValidationContainer> ThinInDetClustersAlg::m_measurementsSctSGKey { this, "InDetTrackMeasurementsSctKey", "SCT_Clusters", "" }

private

Definition at line 94 of file ThinInDetClustersAlg.h.

{ this, "InDetTrackMeasurementsSctKey", "SCT_Clusters", "" };

m_measurementsTrtSGKey

SG::ThinningHandleKey<xAOD::TrackMeasurementValidationContainer> ThinInDetClustersAlg::m_measurementsTrtSGKey { this, "InDetTrackMeasurementsTrtKey", "TRT_DriftCircles", "" }

private

Definition at line 104 of file ThinInDetClustersAlg.h.

{ this, "InDetTrackMeasurementsTrtKey", "TRT_DriftCircles", "" };

m_npass

std::atomic<unsigned int> ThinInDetClustersAlg::m_npass {}

mutableprivate

Definition at line 53 of file ThinInDetClustersAlg.h.

{};

m_npass_pix_measurements

std::atomic<unsigned int> ThinInDetClustersAlg::m_npass_pix_measurements {}

mutableprivate

Definition at line 83 of file ThinInDetClustersAlg.h.

{};

m_npass_pix_states

std::atomic<unsigned int> ThinInDetClustersAlg::m_npass_pix_states {}

mutableprivate

Definition at line 78 of file ThinInDetClustersAlg.h.

{};

m_npass_sct_measurements

std::atomic<unsigned int> ThinInDetClustersAlg::m_npass_sct_measurements {}

mutableprivate

Definition at line 93 of file ThinInDetClustersAlg.h.

{};

m_npass_sct_states

std::atomic<unsigned int> ThinInDetClustersAlg::m_npass_sct_states {}

mutableprivate

Definition at line 88 of file ThinInDetClustersAlg.h.

{};

m_npass_trt_measurements

std::atomic<unsigned int> ThinInDetClustersAlg::m_npass_trt_measurements {}

mutableprivate

Definition at line 103 of file ThinInDetClustersAlg.h.

{};

m_npass_trt_states

std::atomic<unsigned int> ThinInDetClustersAlg::m_npass_trt_states {}

mutableprivate

Definition at line 98 of file ThinInDetClustersAlg.h.

{};

m_ntot

std::atomic<unsigned int> ThinInDetClustersAlg::m_ntot {}

mutableprivate

Counters and keys for xAOD::TrackParticle container.

Definition at line 52 of file ThinInDetClustersAlg.h.

{};

m_ntot_pix_measurements

std::atomic<unsigned int> ThinInDetClustersAlg::m_ntot_pix_measurements {}

mutableprivate

Definition at line 82 of file ThinInDetClustersAlg.h.

{};

m_ntot_pix_states

std::atomic<unsigned int> ThinInDetClustersAlg::m_ntot_pix_states {}

mutableprivate

Counters and keys for xAOD::TrackStateValidation and xAOD::TrackMeasurementValidation containers.

Definition at line 77 of file ThinInDetClustersAlg.h.

{};

m_ntot_sct_measurements

std::atomic<unsigned int> ThinInDetClustersAlg::m_ntot_sct_measurements {}

mutableprivate

Definition at line 92 of file ThinInDetClustersAlg.h.

{};

m_ntot_sct_states

std::atomic<unsigned int> ThinInDetClustersAlg::m_ntot_sct_states {}

mutableprivate

Definition at line 87 of file ThinInDetClustersAlg.h.

{};

m_ntot_trt_measurements

std::atomic<unsigned int> ThinInDetClustersAlg::m_ntot_trt_measurements {}

mutableprivate

Definition at line 102 of file ThinInDetClustersAlg.h.

{};

m_ntot_trt_states

std::atomic<unsigned int> ThinInDetClustersAlg::m_ntot_trt_states {}

mutableprivate

Definition at line 97 of file ThinInDetClustersAlg.h.

{};

m_parser

std::conditional<1==1,std::unique_ptr<ExpressionParsing::ExpressionParser>,std::array<std::unique_ptr<ExpressionParsing::ExpressionParser>,1>>::type ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::m_parser

protectedinherited

Definition at line 100 of file ExpressionParserUser.h.

m_proxyLoaders

std::unique_ptr<ExpressionParsing::IProxyLoader> ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::m_proxyLoaders

protectedinherited

Definition at line 97 of file ExpressionParserUser.h.

m_renounce

std::vector<std::string> ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::m_renounce

protectedinherited

Definition at line 102 of file ExpressionParserUser.h.

m_SCTDetEleCollKey

SG::ReadCondHandleKey<InDetDD::SiDetectorElementCollection> ThinInDetClustersAlg::m_SCTDetEleCollKey {this, "SCTDetEleCollKey", "SCT_DetectorElementCollection", "Key of SiDetectorElementCollection for SCT"}

private

For P->T converter of SCT_Clusters.

Definition at line 108 of file ThinInDetClustersAlg.h.

{this, "SCTDetEleCollKey", "SCT_DetectorElementCollection", "Key of SiDetectorElementCollection for SCT"};

m_selectionString

StringProperty ThinInDetClustersAlg::m_selectionString { this, "SelectionString", "InDetTrackParticles.pt>10*GeV", "Selection string for each TrackParticle container" }

private

Expressions for object thinning selections Default InDetTrackParticles (also accomplished through flag.Tracking.thinPixelClustersSelectionStrings defaults).

Definition at line 69 of file ThinInDetClustersAlg.h.

{ this, "SelectionString", "InDetTrackParticles.pt>10*GeV", "Selection string for each TrackParticle container" };

m_statesPixSGKey

SG::ThinningHandleKey<xAOD::TrackStateValidationContainer> ThinInDetClustersAlg::m_statesPixSGKey { this, "InDetTrackStatesPixKey", "PixelMSOSs", "" }

private

Definition at line 79 of file ThinInDetClustersAlg.h.

{ this, "InDetTrackStatesPixKey", "PixelMSOSs", "" }; // original

m_statesSctSGKey

SG::ThinningHandleKey<xAOD::TrackStateValidationContainer> ThinInDetClustersAlg::m_statesSctSGKey { this, "InDetTrackStatesSctKey", "SCT_MSOSs", "" }

private

Definition at line 89 of file ThinInDetClustersAlg.h.

{ this, "InDetTrackStatesSctKey", "SCT_MSOSs", "" };

m_statesTrtSGKey

SG::ThinningHandleKey<xAOD::TrackStateValidationContainer> ThinInDetClustersAlg::m_statesTrtSGKey { this, "InDetTrackStatesTrtKey", "TRT_MSOSs", "" }

private

Definition at line 99 of file ThinInDetClustersAlg.h.

{ this, "InDetTrackStatesTrtKey", "TRT_MSOSs", "" };

m_streamName

StringProperty ThinInDetClustersAlg::m_streamName { this, "StreamName", "", "Name of the stream being thinned" }

private

Stream for object thinning selections.

Definition at line 63 of file ThinInDetClustersAlg.h.

{ this, "StreamName", "", "Name of the stream being thinned" };

m_thinPixelHitsOnTrack

BooleanProperty ThinInDetClustersAlg::m_thinPixelHitsOnTrack { this, "ThinPixelHitsOnTrack", false, ""}

private

Thinning logic.

Definition at line 55 of file ThinInDetClustersAlg.h.

{ this, "ThinPixelHitsOnTrack", false, ""};

m_thinSCTHitsOnTrack

BooleanProperty ThinInDetClustersAlg::m_thinSCTHitsOnTrack { this, "ThinSCTHitsOnTrack", false, ""}

private

Definition at line 57 of file ThinInDetClustersAlg.h.

{ this, "ThinSCTHitsOnTrack", false, ""};

m_thinTRTHitsOnTrack

BooleanProperty ThinInDetClustersAlg::m_thinTRTHitsOnTrack { this, "ThinTRTHitsOnTrack", false, ""}

private

Definition at line 59 of file ThinInDetClustersAlg.h.

{ this, "ThinTRTHitsOnTrack", false, ""};

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.

m_xAODProxyLoader

ExpressionParsing::SGxAODProxyLoader* ExpressionParserUserBase< ::AthReentrantAlgorithm, 1 >::m_xAODProxyLoader

protectedinherited

The expression parser.

Definition at line 96 of file ExpressionParserUser.h.

---

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

• ThinInDetClustersAlg.h
• ThinInDetClustersAlg.cxx