#include <SeedAnalysisAlg.h>

Inheritance diagram for ActsTrk::SeedAnalysisAlg:

Collaboration diagram for ActsTrk::SeedAnalysisAlg:

Public Types
enum class	Environment_t { user = 0 , online , tier0 , tier0Raw , tier0ESD , AOD , altprod }
	Specifies the processing environment. More...
enum class	DataType_t { userDefined = 0 , monteCarlo , collisions , cosmics , heavyIonCollisions }
	Specifies what type of input data is being monitored. More...

Public Member Functions
	SeedAnalysisAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~SeedAnalysisAlg () override=default
virtual StatusCode	initialize () override
	initialize
virtual StatusCode	fillHistograms (const EventContext &ctx) const override
	adds event to the monitoring histograms
virtual StatusCode	execute (const EventContext &ctx) const override
	Applies filters and trigger requirements.
void	fill (const ToolHandle< GenericMonitoringTool > &groupHandle, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
	Fills a vector of variables to a group by reference.
void	fill (const ToolHandle< GenericMonitoringTool > &groupHandle, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
	Fills a vector of variables to a group by reference.
template<typename... T>
void	fill (const ToolHandle< GenericMonitoringTool > &groupHandle, T &&... variables) const
	Fills a variadic list of variables to a group by reference.
void	fill (const std::string &groupName, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
	Fills a vector of variables to a group by name.
void	fill (const std::string &groupName, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
	Fills a vector of variables to a group by name.
template<typename... T>
void	fill (const std::string &groupName, T &&... variables) const
	Fills a variadic list of variables to a group by name.
Environment_t	environment () const
	Accessor functions for the environment.
Environment_t	envStringToEnum (const std::string &str) const
	Convert the environment string from the python configuration to an enum object.
DataType_t	dataType () const
	Accessor functions for the data type.
DataType_t	dataTypeStringToEnum (const std::string &str) const
	Convert the data type string from the python configuration to an enum object.
const ToolHandle< GenericMonitoringTool > &	getGroup (const std::string &name) const
	Get a specific monitoring tool from the tool handle array.
const ToolHandle< Trig::TrigDecisionTool > &	getTrigDecisionTool () const
	Get the trigger decision tool member.
bool	trigChainsArePassed (const std::vector< std::string > &vTrigNames) const
	Check whether triggers are passed.
SG::ReadHandle< xAOD::EventInfo >	GetEventInfo (const EventContext &) const
	Return a ReadHandle for an EventInfo object (get run/event numbers, etc.).
virtual float	lbAverageInteractionsPerCrossing (const EventContext &ctx) const
	Calculate the average mu, i.e.
virtual float	lbInteractionsPerCrossing (const EventContext &ctx) const
	Calculate instantaneous number of interactions, i.e.
virtual float	lbAverageLuminosity (const EventContext &ctx) const
	Calculate average luminosity (in ub-1 s-1 => 10^30 cm-2 s-1).
virtual float	lbLuminosityPerBCID (const EventContext &ctx) const
	Calculate the instantaneous luminosity per bunch crossing.
virtual double	lbDuration (const EventContext &ctx) const
	Calculate the duration of the luminosity block (in seconds).
virtual float	lbAverageLivefraction (const EventContext &ctx) const
	Calculate the average luminosity livefraction.
virtual float	livefractionPerBCID (const EventContext &ctx) const
	Calculate the live fraction per bunch crossing ID.
virtual double	lbLumiWeight (const EventContext &ctx) const
	Calculate the average integrated luminosity multiplied by the live fraction.
virtual StatusCode	parseList (const std::string &line, std::vector< std::string > &result) const
	Parse a string into a vector.
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
	Get filter decision:
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
	Set filter decision:
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

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

Protected Attributes
ToolHandleArray< GenericMonitoringTool >	m_tools {this,"GMTools",{}}
	Array of Generic Monitoring Tools.
PublicToolHandle< Trig::TrigDecisionTool >	m_trigDecTool
	Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterTool >	m_DQFilterTools {this,"FilterTools",{}}
	Array of Data Quality filter tools.
SG::ReadCondHandleKey< LuminosityCondData >	m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
SG::ReadCondHandleKey< LBDurationCondData >	m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
SG::ReadCondHandleKey< TrigLiveFractionCondData >	m_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}
AthMonitorAlgorithm::Environment_t	m_environment
	Instance of the Environment_t enum.
AthMonitorAlgorithm::DataType_t	m_dataType
	Instance of the DataType_t enum.
Gaudi::Property< std::string >	m_environmentStr {this,"Environment","user"}
	Environment string pulled from the job option and converted to enum.
Gaudi::Property< std::string >	m_dataTypeStr {this,"DataType","userDefined"}
	DataType string pulled from the job option and converted to enum.
Gaudi::Property< std::string >	m_triggerChainString {this,"TriggerChain",""}
	Trigger chain string pulled from the job option and parsed into a vector.
std::vector< std::string >	m_vTrigChainNames
	Vector of trigger chain names parsed from trigger chain string.
Gaudi::Property< std::string >	m_fileKey {this,"FileKey",""}
	Internal Athena name for file.
Gaudi::Property< bool >	m_useLumi {this,"EnableLumi",false}
	Allows use of various luminosity functions.
Gaudi::Property< float >	m_defaultLBDuration {this,"DefaultLBDuration",60.}
	Default duration of one lumi block.
Gaudi::Property< int >	m_detailLevel {this,"DetailLevel",0}
	Sets the level of detail used in the monitoring.
SG::ReadHandleKey< xAOD::EventInfo >	m_EventInfoKey {this,"EventInfoKey","EventInfo"}
	Key for retrieving EventInfo from StoreGate.

Private Types
typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > >	MonVarVec_t
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	fillTruthHistograms (const EventContext &ctx, const ActsTrk::SeedContainer &seed_container, std::vector< int > &, std::vector< double > &) const
const Identifier	identify (const xAOD::PixelCluster &) const
const Identifier	identify (const xAOD::StripCluster &) const
void	matchParticleToSeedClusters (const PRD_MultiTruthCollection *prdTruth, const Identifier &id, std::map< int, int > &countMap) const
std::pair< int, double >	findSeedMajorityTruthParticle (const std::map< int, int > &countMap) const
std::array< float, 7 >	estimateParameters (const ActsTrk::Seed &seed, float pTPerHelixRadius) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ToolHandle< ActsTrk::ITrackParamsEstimationTool >	m_paramEstimationTool {this, "TrackParamsEstimationTool", "", "Track Param Estimation from Seeds"}
PublicToolHandle< IGeometryRealmConvTool >	m_geometryConvTool {this, "GeometryRealmConvTool", ""}
SG::ReadHandleKey< ActsTrk::SeedContainer >	m_inputSeedColletionKey {this, "InputSeedCollection", "", ""}
ServiceHandle< ActsTrk::ITrackingGeometrySvc >	m_trackingGeometrySvc {this, "TrackingGeometrySvc", "ActsTrackingGeometrySvc"}
SG::ReadHandleKey< PRD_MultiTruthCollection >	m_prdTruth {this, "ITkClustersTruth", "", ""}
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollection >	m_detEleCollKey {this, "DetectorElements", "", "Key of input SiDetectorElementCollection"}
SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"}
SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCondObjInputKey
Gaudi::Property< std::string >	m_monGroupName {this, "MonGroupName", "ActsSeedAnalysisAlg"}
Gaudi::Property< bool >	m_useTopSp {this, "useTopSp", false, "Use top SP. By default, use bottom SP."}
Gaudi::Property< bool >	m_usePixel {this, "UsePixel", true, ""}
ContextUtility	m_ctxProvider {this}
	Auxiliary class to access the calibration context.
std::string	m_name
std::unordered_map< std::string, size_t >	m_toolLookupMap
const ToolHandle< GenericMonitoringTool >	m_dummy
Gaudi::Property< bool >	m_enforceExpressTriggers
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 36 of file SeedAnalysisAlg.h.

Member Typedef Documentation

◆ MonVarVec_t

typedef std::vector<std::reference_wrapper<Monitored::IMonitoredVariable> > AthMonitorAlgorithm::MonVarVec_t

privateinherited

Definition at line 370 of file AthMonitorAlgorithm.h.

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ DataType_t

enum class AthMonitorAlgorithm::DataType_t

stronginherited

Specifies what type of input data is being monitored.

An enumeration of the different types of data the monitoring application may be running over. This can be used to select which histograms to produce, e.g. to prevent the production of colliding-beam histograms when running on cosmic-ray data. Strings of the same names may be given as jobOptions.

Enumerator
userDefined
monteCarlo
collisions
cosmics
heavyIonCollisions

Definition at line 194 of file AthMonitorAlgorithm.h.

                          {
        userDefined = 0, 
        monteCarlo, 
        collisions, 
        cosmics, 
        heavyIonCollisions, 
    };

◆ Environment_t

enum class AthMonitorAlgorithm::Environment_t

stronginherited

Specifies the processing environment.

The running environment may be used to select which histograms are produced, and where they are located in the output. For example, the output paths of the histograms are different for the "user", "online" and the various offline flags. Strings of the same names may be given as jobOptions.

Enumerator
user
online
tier0
tier0Raw
tier0ESD
AOD
altprod

Definition at line 175 of file AthMonitorAlgorithm.h.

                             {
        user = 0, 
        online, 
        tier0, 
        tier0Raw, 
        tier0ESD, 
        AOD, 
        altprod, 
    };

Constructor & Destructor Documentation

◆ SeedAnalysisAlg()

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

Definition at line 20 of file SeedAnalysisAlg.cxx.

21 : AthMonitorAlgorithm(name, pSvcLocator)

22 {}

AthMonitorAlgorithm::AthMonitorAlgorithm

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

Constructor.

Definition AthMonitorAlgorithm.cxx:8

◆ ~SeedAnalysisAlg()

virtual ActsTrk::SeedAnalysisAlg::~SeedAnalysisAlg ( )

overridevirtualdefault

Member Function Documentation

◆ dataType()

DataType_t AthMonitorAlgorithm::dataType ( ) const

inlineinherited

Accessor functions for the data type.

Returns: the current value of the class's DataType_t instance.

Definition at line 224 of file AthMonitorAlgorithm.h.

224{ return m_dataType; }

AthMonitorAlgorithm::m_dataType

AthMonitorAlgorithm::DataType_t m_dataType

Instance of the DataType_t enum.

Definition AthMonitorAlgorithm.h:356

◆ dataTypeStringToEnum()

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::dataTypeStringToEnum ( const std::string & str ) const

inherited

Convert the data type string from the python configuration to an enum object.

Returns: a value in the DataType_t enumeration which matches the input string.

Definition at line 144 of file AthMonitorAlgorithm.cxx.

                                                                                                    {
    // convert the string to all lowercase
    std::string lowerCaseStr = str;
    std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
 
    // check if it matches one of the enum choices
    if( lowerCaseStr == "userdefined" ) {
        return DataType_t::userDefined;
    } else if( lowerCaseStr == "montecarlo" ) {
        return DataType_t::monteCarlo;
    } else if( lowerCaseStr == "collisions" ) {
        return DataType_t::collisions;
    } else if( lowerCaseStr == "cosmics" ) {
        return DataType_t::cosmics;
    } else if( lowerCaseStr == "heavyioncollisions" ) {
        return DataType_t::heavyIonCollisions;
    } else { // otherwise, warn the user and return "userDefined"
        ATH_MSG_WARNING("AthMonitorAlgorithm::dataTypeStringToEnum(): Unknown data type "
            <<str<<", returning userDefined.");
        return DataType_t::userDefined;
    }
}

◆ 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.

95{ return m_detStore; }

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

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ environment()

Environment_t AthMonitorAlgorithm::environment ( ) const

inlineinherited

Accessor functions for the environment.

Returns: the current value of the class's Environment_t instance.

Definition at line 208 of file AthMonitorAlgorithm.h.

208{ return m_environment; }

AthMonitorAlgorithm::m_environment

AthMonitorAlgorithm::Environment_t m_environment

Instance of the Environment_t enum.

Definition AthMonitorAlgorithm.h:355

◆ envStringToEnum()

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::envStringToEnum ( const std::string & str ) const

inherited

Convert the environment string from the python configuration to an enum object.

Returns: a value in the Environment_t enumeration which matches the input string.

Definition at line 116 of file AthMonitorAlgorithm.cxx.

                                                                                                  {
    // convert the string to all lowercase
    std::string lowerCaseStr = str;
    std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
 
    // check if it matches one of the enum choices
    if( lowerCaseStr == "user" ) {
        return Environment_t::user;
    } else if( lowerCaseStr == "online" ) {
        return Environment_t::online;
    } else if( lowerCaseStr == "tier0" ) {
        return Environment_t::tier0;
    } else if( lowerCaseStr == "tier0raw" ) {
        return Environment_t::tier0Raw;
    } else if( lowerCaseStr == "tier0esd" ) {
        return Environment_t::tier0ESD;
    } else if( lowerCaseStr == "aod" ) {
        return Environment_t::AOD;
    } else if( lowerCaseStr == "altprod" ) {
        return Environment_t::altprod;
    } else { // otherwise, warn the user and return "user"
        ATH_MSG_WARNING("AthMonitorAlgorithm::envStringToEnum(): Unknown environment "
            <<str<<", returning user.");
        return Environment_t::user;
    }
}

◆ estimateParameters()

std::array< float, 7 > ActsTrk::SeedAnalysisAlg::estimateParameters	(	const ActsTrk::Seed &	seed,
		float	pTPerHelixRadius ) const

private

Definition at line 379 of file SeedAnalysisAlg.cxx.

  {
    auto extractCoordinates = 
      [] (const xAOD::SpacePoint* sp) -> std::array<float,4>
      {
    std::array<float, 4> coordinates {static_cast<float>(sp->x()), 
                      static_cast<float>(sp->y()), 
                      static_cast<float>(sp->z()), 
                      static_cast<float>(std::sqrt(sp->x()*sp->x() + sp->y()*sp->y()))};
    return coordinates;
      };
    
    auto extractQuantities = 
      [] (const std::array<float, 4>& sp,
      const std::array<float, 4>& spM,
      bool isBottom) -> std::array<float, 5>
      {
    auto& [xM, yM, zM, rM] = spM;
    auto& [xO, yO, zO, rO] = sp;
 
    float cosPhiM = xM / rM;
    float sinPhiM = yM / rM;
    float deltaX = xO - xM;
    float deltaY = yO - yM;
    float deltaZ = zO - zM;
    float x = deltaX * cosPhiM + deltaY * sinPhiM;
    float y = deltaY * cosPhiM - deltaX * sinPhiM;
    float iDeltaR2 = 1.f / (deltaX * deltaX + deltaY * deltaY);
    float iDeltaR = std::sqrt(iDeltaR2);
    int bottomFactor = int(not isBottom) - int(isBottom);
    float cot_theta = deltaZ * iDeltaR * bottomFactor;
 
    // cotTheta, Zo, iDeltaR, U, V
    std::array<float, 5> params =
    { 
      cot_theta, 
      zM - rM * cot_theta,
      iDeltaR,
      x * iDeltaR2,
        y * iDeltaR2
    };
 
    return params;
      };
 
    const auto& sps = seed.sp();
    const auto* bottom = sps[0];
    const auto* medium = sps[1];
    const auto* top = sps[2];
 
    auto coo_b = extractCoordinates(bottom);
    auto coo_m = extractCoordinates(medium);
    auto coo_t = extractCoordinates(top);
 
    // Compute the variables we need
    auto [cotThetaB, Zob, iDeltaRB, Ub, Vb] = extractQuantities(coo_b, coo_m, true);
    auto [cotThetaT, Zot, iDeltaRT, Ut, Vt] = extractQuantities(coo_t, coo_m, false);
 
    float squarediDeltaR2B = iDeltaRB*iDeltaRB;
    float squarediDeltaR2T = iDeltaRB*iDeltaRT;
    float squarediDeltaR = std::min(squarediDeltaR2B, squarediDeltaR2T);
 
    auto& [xB, yB, zB, rB] = coo_b;
    auto& [xM, yM, zM, rM] = coo_m;
    auto& [xT, yT, zT, rT] = coo_t;
 
    float ax = xM / rM;
    float ay = yM/ rM;
 
    float dxb = xM - xB;
    float dyb = yM - yB;
    float dzb = zM - zB;
    float xb = dxb * ax + dyb *ay;
    float yb = dyb * ax - dxb * ay; 
    float dxyb = xb * xb + yb * yb;                                                                                                                                                                                                                                            
    float dxt = xT - xM;
    float dyt = yT - yM;
    float dzt = zT - zM;
    float xt = dxt * ax + dyt *ay;
    float yt = dyt * ax - dxt * ay;
    float dxyt = xt * xt + yt * yt;
 
    float tzb = dzb * std::sqrt( 1.f/dxyb );
    float tzt = dzt * std::sqrt( 1.f/dxyt );
 
    float sTzb2 = std::sqrt(1.f + tzb*tzb);
 
    float dU = Ut - Ub;
    if (dU == 0.) {
      return {-1, -1, -1, -1, -1, -1, -1};
    }
 
    float A = (Vt - Vb) / dU;
    float S2 = 1.f + A * A;
    float B = Vb - A * Ub;
    float B2 = B * B;
    if (B2 == 0) B2 = 1e-8;
 
    // dzdr
    float dzdr_b = (zM - zB) / (rM - rB);
    float dzdr_t = (zT - zM) / (rT - rM);    
 
    // eta
    float cotThetaAvg2 = cotThetaB * cotThetaT;
    if (cotThetaAvg2 <= 0) {
      return {-1, -1, -1, -1, -1, -1, -1};
    }
    float theta = std::atan(1.f / std::sqrt(cotThetaAvg2));
    float eta = -std::log(std::tan(0.5f * theta));
 
    // pt
    float pt = pTPerHelixRadius * std::sqrt(S2 / B2) / 2.f;
 
    // d0
    float d0 = std::abs((A - B * rM) * rM);
 
    // curvature 
    // not used in current version of the code. We may want to use it later
    //    float curvature = B / std::sqrt(S2);
    float penalty = std::abs((tzb - tzt) / (squarediDeltaR * sTzb2));
 
    return {pt, theta, eta, d0, dzdr_b, dzdr_t, penalty};
  }

◆ 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.

85{ return m_evtStore; }

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

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ execute()

StatusCode AthMonitorAlgorithm::execute ( const EventContext & ctx ) const

overridevirtualinherited

Applies filters and trigger requirements.

Then, calls fillHistograms().

Parameters

ctx	event context for reentrant Athena call

Returns: StatusCode

Definition at line 77 of file AthMonitorAlgorithm.cxx.

                                                                       {
 
    // Checks that all of the  DQ filters are passed. If any one of the filters
    // fails, return SUCCESS code and do not fill the histograms with the event.
    for ( const auto& filterItr : m_DQFilterTools ) {
        if (!filterItr->accept(ctx)) {
            ATH_MSG_DEBUG("Event rejected due to filter tool.");
            return StatusCode::SUCCESS;
        }
    }
 
    // Trigger: If there is a decision tool and the chains fail, skip the event.
    if ( !m_trigDecTool.empty() && !trigChainsArePassed(m_vTrigChainNames) ) {
        ATH_MSG_DEBUG("Event rejected due to trigger filter.");
        return StatusCode::SUCCESS;
    }
 
    ATH_MSG_DEBUG("Event accepted!");
    return fillHistograms(ctx);
}

◆ extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

◆ extraOutputDeps()

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

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 89 of file AthCommonAlgorithm.cxx.

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

◆ fillHistograms()

StatusCode ActsTrk::SeedAnalysisAlg::fillHistograms ( const EventContext & ctx ) const

overridevirtual

adds event to the monitoring histograms

User will overwrite this function. Histogram booking is no longer done in C++. This function is called in execute once the filters are all passed.

Parameters

ctx	forwarded from execute

Returns: StatusCode

Implements AthMonitorAlgorithm.

Definition at line 48 of file SeedAnalysisAlg.cxx.

                                                                          {
    ATH_MSG_DEBUG( "Filling Histograms for " << name() << " ... " );
 
    // CONDS
    // Read the Beam Spot information
    SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey, ctx };
    ATH_CHECK( beamSpotHandle.isValid() );
    auto beamSpotData = beamSpotHandle.cptr();
    
    // Read the b-field information
    SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle { m_fieldCondObjInputKey, ctx };
    ATH_CHECK( readHandle.isValid() );
    const AtlasFieldCacheCondObj* fieldCondObj{ *readHandle };
    
    // Get the magnetic field
    // Using ACTS classes in order to be sure we are consistent
    Acts::MagneticFieldContext magFieldContext(fieldCondObj);
    // Beam Spot Position
    Acts::Vector2 beamPos( beamSpotData->beamPos()[ Amg::x ] * Acts::UnitConstants::mm,
               beamSpotData->beamPos()[ Amg::y ] * Acts::UnitConstants::mm );
    
    // Magnetic Field
    ATLASMagneticFieldWrapper magneticField;
    Acts::MagneticFieldProvider::Cache magFieldCache = magneticField.makeCache( magFieldContext );
    Acts::Vector3 bField = *magneticField.getField( Acts::Vector3(beamPos.x(), beamPos.y(), 0),
                            magFieldCache );
 
 
    // SEEDS
    ATH_MSG_DEBUG( "Reading input collection with key " << m_inputSeedColletionKey.key() );
    SG::ReadHandle< ActsTrk::SeedContainer > handle = SG::makeHandle( m_inputSeedColletionKey, ctx );
    ATH_CHECK( handle.isValid() );
    const ActsTrk::SeedContainer* seed_collection = handle.get();
    ATH_MSG_DEBUG( "Retrieved " << seed_collection->size() << " input elements from key " << m_inputSeedColletionKey.key() );
    
    auto monitor_nseed = Monitored::Scalar<int>("Nseed", seed_collection->size());
    fill(m_monGroupName.value(), monitor_nseed);
 
    // bottom  
    auto monitor_x1 = 
      Monitored::Collection("x1", *seed_collection,
                [] (const auto& seed) -> double
                { return seed.sp()[0]->x(); });
    auto monitor_y1 = 
      Monitored::Collection("y1", *seed_collection,
                [] (const auto& seed) -> double
                { return seed.sp()[0]->y(); });
    auto monitor_z1 = 
      Monitored::Collection("z1", *seed_collection,
                [] (const auto& seed) -> double
                { return seed.sp()[0]->z(); });
    auto monitor_r1 =
      Monitored::Collection("r1", *seed_collection,
                [] (const auto& seed) -> double
                {
                  const auto* sp = seed.sp()[0];
                  return std::sqrt(sp->x()*sp->x() + sp->y()*sp->y());
                });
 
    // middle
    auto monitor_x2 =
      Monitored::Collection("x2", *seed_collection,
                [] (const auto& seed) -> double
                { return seed.sp()[1]->x(); });
    auto monitor_y2 =
      Monitored::Collection("y2", *seed_collection,
                            [] (const auto& seed) -> double
                            { return seed.sp()[1]->y(); });
    auto monitor_z2 =
      Monitored::Collection("z2", *seed_collection,
                            [] (const auto& seed) -> double
                            { return seed.sp()[1]->z(); });
    auto monitor_r2 =
      Monitored::Collection("r2", *seed_collection,
                            [] (const auto& seed) -> double
                            {
                              const auto* sp = seed.sp()[1];
                              return std::sqrt(sp->x()*sp->x() + sp->y()*sp->y());
                            });
    
    // top
    auto monitor_x3 =
      Monitored::Collection("x3", *seed_collection,
                [] (const auto& seed) -> double
                { return seed.sp()[2]->x(); });
    auto monitor_y3 =
      Monitored::Collection("y3", *seed_collection,
                            [] (const auto& seed) -> double
                            { return seed.sp()[2]->y(); });
    auto monitor_z3 =
      Monitored::Collection("z3", *seed_collection,
                            [] (const auto& seed) -> double
                            { return seed.sp()[2]->z(); });
    auto monitor_r3 =
      Monitored::Collection("r3", *seed_collection,
                            [] (const auto& seed) -> double
                            {
                              const auto* sp = seed.sp()[2];
                              return std::sqrt(sp->x()*sp->x() + sp->y()*sp->y());
                            });
 
    std::vector< std::array<float, 7> > parametersCollection;
    parametersCollection.reserve(seed_collection->size());
 
    for (auto seed : *seed_collection) {
      parametersCollection.push_back( estimateParameters(seed, 300. * bField[2] / 1000.) );
    }
    
    auto monitor_param_pt = Monitored::Collection("pt", parametersCollection,
                          [] (const auto& params) -> float
                          { return params[0]; });
    auto monitor_param_theta = Monitored::Collection("theta", parametersCollection,
                             [] (const auto& params) -> float
                             { return params[1]; });
    auto monitor_param_eta = Monitored::Collection("eta", parametersCollection,
                           [] (const auto& params) -> float
                           { return params[2]; });
    auto monitor_param_d0 = Monitored::Collection("d0", parametersCollection,
                          [] (const auto& params) -> float
                          { return params[3]; });
 
    auto monitor_param_dzdr_b = Monitored::Collection("dzdr_b", parametersCollection,
                              [] (const auto& params) -> float
                              { return params[4]; });
    auto monitor_param_dzdr_t = Monitored::Collection("dzdr_t", parametersCollection,
                              [] (const auto& params) -> float
                              { return params[5]; });
 
    
    auto monitor_param_penalty = Monitored::Collection("penalty", parametersCollection,
                               [] (const auto& params) -> float
                               { return params[6]; });
      
 
    SG::ReadHandle<xAOD::EventInfo> eventInfo = SG::makeHandle(m_EventInfoKey, ctx);
    ATH_CHECK(eventInfo.isValid());
 
    auto monitor_event_number = Monitored::Scalar<long>("event_number", static_cast<long>(eventInfo->eventNumber()));
    auto monitor_actual_mu = Monitored::Scalar<float>("actual_mu", eventInfo->actualInteractionsPerCrossing());
 
    std::vector<int> vec_truthBarcode;
    std::vector<double> vec_truthProb;
    if (not m_prdTruth.empty())
      ATH_CHECK( fillTruthHistograms(ctx, *seed_collection, vec_truthBarcode, vec_truthProb) );
    auto monitor_truth_barcode = Monitored::Collection("truth_barcode", vec_truthBarcode);
    auto monitor_truth_prob = Monitored::Collection("truth_prob", vec_truthProb);
 
    fill(m_monGroupName.value(),
     monitor_x1, monitor_y1, monitor_z1, monitor_r1,
     monitor_x2, monitor_y2, monitor_z2, monitor_r2,
     monitor_x3, monitor_y3, monitor_z3, monitor_r3,
     monitor_param_pt, monitor_param_theta, monitor_param_eta, monitor_param_d0,
     monitor_param_dzdr_b, monitor_param_dzdr_t,
     monitor_param_penalty,
     monitor_event_number, monitor_actual_mu,
     monitor_truth_barcode, monitor_truth_prob);
        
    return StatusCode::SUCCESS;
  }

◆ fillTruthHistograms()

StatusCode ActsTrk::SeedAnalysisAlg::fillTruthHistograms	(	const EventContext &	ctx,
		const ActsTrk::SeedContainer &	seed_container,
		std::vector< int > &	truthBarCodeVec,
		std::vector< double > &	truthProbVec ) const

private

Definition at line 208 of file SeedAnalysisAlg.cxx.

  {
    ATH_MSG_DEBUG( "Filling Truth Histograms for " << name() << " ... " );
 
    SG::ReadHandle< PRD_MultiTruthCollection > prdTruthHandle = SG::makeHandle( m_prdTruth, ctx );
    ATH_CHECK(prdTruthHandle.isValid());
    const PRD_MultiTruthCollection* prdTruth = prdTruthHandle.get();
 
    SG::ReadCondHandle< InDetDD::SiDetectorElementCollection > detEleHandle( m_detEleCollKey, ctx );
    ATH_CHECK( detEleHandle.isValid() );
    const InDetDD::SiDetectorElementCollection& detElements = *detEleHandle.retrieve();
 
    // Read the b-field information
    SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle { m_fieldCondObjInputKey, ctx };
    ATH_CHECK( readHandle.isValid() );
    
    const AtlasFieldCacheCondObj* fieldCondObj{ *readHandle };
    ATH_CHECK( fieldCondObj != nullptr );
 
    Acts::MagneticFieldContext magFieldContext(fieldCondObj);
    const auto & geo_context = m_trackingGeometrySvc->getNominalContext();
    const Acts::CalibrationContext calContext = m_ctxProvider.getCalibrationContext(ctx);
 
    // utilities
    // Used for param estimation
    auto retrieveSurfaceFunction = 
      [this, &detElements] (const ActsTrk::Seed& seed, bool useTopSp) -> const Acts::Surface& 
      { 
          const xAOD::SpacePoint* sp = useTopSp ? seed.sp().back() : seed.sp().front();
          const InDetDD::SiDetectorElement* element = detElements.getDetectorElement(
                useTopSp ? sp->elementIdList().back() : sp->elementIdList().front());
          const Trk::Surface& atlas_surface = element->surface();
          return *m_geometryConvTool->convertSurfaceToActs(atlas_surface);
      };
 
 
    // computation
    std::vector<bool> vec_pass;
    vec_pass.reserve(seed_container.size());
 
    std::vector<double> estimated_pt;
    std::vector<double> estimated_eta;
    estimated_pt.reserve(seed_container.size());
    estimated_eta.reserve(seed_container.size());
 
    for (auto seed : seed_container) {
      const auto& [optTrackParams, estimationStatus] =
        m_paramEstimationTool->estimateTrackParameters(
                               seed,
                               m_useTopSp,
                               geo_context.context(),
                               magFieldContext,
                               calContext,
                               retrieveSurfaceFunction);
 
      if ( not optTrackParams.has_value() ) continue;
 
      const auto & param = optTrackParams.value();
      estimated_pt.push_back( param.transverseMomentum() );
      estimated_eta.push_back( -std::log( std::tan(0.5 * param.parameters()[Acts::eBoundTheta]) ) );
 
      std::map<int, int> truthHits;
      
      const auto& sps = seed.sp();
      for (const auto* sp : sps) {
    int number_of_clusters = m_usePixel ? 1 : 2;
    for (int cluster_number(0); cluster_number < number_of_clusters; cluster_number++) {
      const auto& els = sp->measurements();
      const auto* cluster = els[cluster_number];
      const xAOD::UncalibMeasType cluster_type = cluster->type();
      const Identifier id = cluster_type == xAOD::UncalibMeasType::PixelClusterType
        ? identify(*reinterpret_cast<const xAOD::PixelCluster*>(cluster))
        : identify(*reinterpret_cast<const xAOD::StripCluster*>(cluster));
      matchParticleToSeedClusters(prdTruth, id, truthHits);
    }
      }
 
      auto [barcode, prob] = findSeedMajorityTruthParticle(truthHits);
      truthBarCodeVec.push_back(barcode);
      truthProbVec.push_back(prob);
      vec_pass.push_back( barcode != 0 and prob > 0.5 );
    }
 
    auto monitor_estimated_pt = Monitored::Collection("estimated_pt", estimated_pt);
    auto monitor_estimated_eta = Monitored::Collection("estimated_eta", estimated_eta);
    auto monitor_pass = Monitored::Collection("passed", vec_pass);
 
    fill(m_monGroupName.value(),
     monitor_pass,
     monitor_estimated_pt, monitor_estimated_eta);
 
    return StatusCode::SUCCESS;
  }

◆ filterPassed()

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

inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

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

◆ findSeedMajorityTruthParticle()

std::pair< int, double > ActsTrk::SeedAnalysisAlg::findSeedMajorityTruthParticle ( const std::map< int, int > & countMap ) const

private

Definition at line 361 of file SeedAnalysisAlg.cxx.

                                                                                                            {
    int bestCount = 0;
    int bestBarcode = std::numeric_limits<int>::min();
  
    for (auto const& [barcode, count] : countMap) {
      if (count > bestCount) {
    bestCount = count;
    bestBarcode = barcode;
      }
    }
 
    // 3 spacepoints per seed, 1 (2) clusters per spacepoint for pixel (strip)
    double nTotal = m_usePixel ? 3. : 6.;
    double prob = bestCount / nTotal;
 
    return std::make_pair(bestBarcode, prob);
  }

◆ GetEventInfo()

SG::ReadHandle< xAOD::EventInfo > AthMonitorAlgorithm::GetEventInfo ( const EventContext & ctx ) const

inherited

Return a ReadHandle for an EventInfo object (get run/event numbers, etc.).

Parameters

ctx	EventContext for the event

Returns: a SG::ReadHandle<xAOD::EventInfo>

Definition at line 111 of file AthMonitorAlgorithm.cxx.

                                                                                             {
    return SG::ReadHandle<xAOD::EventInfo>(m_EventInfoKey, ctx);
}

◆ getGroup()

const ToolHandle< GenericMonitoringTool > & AthMonitorAlgorithm::getGroup ( const std::string & name ) const

inherited

Get a specific monitoring tool from the tool handle array.

Finds a specific GenericMonitoringTool instance from the list of monitoring tools (a ToolHandleArray). Throws a FATAL warning if the object found is empty.

Parameters

name	string name of the desired tool

Returns: reference to the desired monitoring tool

Definition at line 168 of file AthMonitorAlgorithm.cxx.

                                                                                                    {
    // get the pointer to the tool, and check that it exists
    auto idx = m_toolLookupMap.find(name);
    //Check if the tool exists in the map
    if (ATH_LIKELY(idx != m_toolLookupMap.end())) {
        return m_tools[idx->second];
    }
    else {
      //Check if the map is empty
     if (m_toolLookupMap.empty()) {
            ATH_MSG_FATAL("The m_toolLookupMap is empty. The tool " << name << " cannot be found in an empty map.");
            return m_dummy;
      } 
      //If the map is not empty and the tool was not found, print a fatal error
      if (!m_toolLookupMap.empty()) {
            std::string available = std::accumulate(
            m_toolLookupMap.begin(), m_toolLookupMap.end(), std::string(""),
            [](const std::string& s, auto h) { return s + "," + h.first; });
            ATH_MSG_FATAL("The tool " << name << " could not be found in the tool array of the "
                      << "monitoring algorithm " << m_name << ". This probably reflects a discrepancy between "
                      << "your python configuration and c++ filling code. Note: your available groups are {"
                      << available << "}.");
 
            return m_dummy;
      } 
      if (!isInitialized()) {
            ATH_MSG_FATAL(
            "It seems that the AthMonitorAlgorithm::initialize was not called "
            "in derived class initialize method, group name: " << name);
      } 
    }
    return m_dummy;
}

◆ getTrigDecisionTool()

const ToolHandle< Trig::TrigDecisionTool > & AthMonitorAlgorithm::getTrigDecisionTool ( ) const

inherited

Get the trigger decision tool member.

The trigger decision tool is used to check whether a specific trigger is passed by an event.

Returns: m_trigDecTool

Definition at line 203 of file AthMonitorAlgorithm.cxx.

                                                                                       {
    return m_trigDecTool;
}

◆ identify() [1/2]

const Identifier ActsTrk::SeedAnalysisAlg::identify ( const xAOD::PixelCluster & cluster ) const

private

Definition at line 305 of file SeedAnalysisAlg.cxx.

  {
    static const SG::ConstAccessor< ElementLink< InDet::PixelClusterCollection > > pixelLinkAcc("pixelClusterLink");
 
    // TO-DO -- AODs will not have this decoration, we'll need to provide a function for recomputing 
    // the identifier from local position
    if (not pixelLinkAcc.isAvailable (cluster))
      return Identifier();
 
    ElementLink<InDet::PixelClusterCollection> pixelLink = pixelLinkAcc(cluster);
    return (*pixelLink)->identify();
  }

◆ identify() [2/2]

const Identifier ActsTrk::SeedAnalysisAlg::identify ( const xAOD::StripCluster & cluster ) const

private

Definition at line 318 of file SeedAnalysisAlg.cxx.

  {
    static const SG::ConstAccessor< ElementLink< InDet::SCT_ClusterCollection > > stripLinkAcc("sctClusterLink");
 
    // TO-DO -- AODs will not have this decoration, we'll need to provide a function for recomputing 
    // the identifier from local position
    if (not stripLinkAcc.isAvailable (cluster))
      return Identifier();
 
    ElementLink<InDet::SCT_ClusterCollection> stripLink = stripLinkAcc(cluster);
    return (*stripLink)->identify();
  }

◆ initialize()

StatusCode ActsTrk::SeedAnalysisAlg::initialize ( )

overridevirtual

initialize

Returns: StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 24 of file SeedAnalysisAlg.cxx.

                                         {
    ATH_MSG_INFO("Initializing " << name() << " ...");
 
    ATH_CHECK( m_beamSpotKey.initialize() );
    ATH_CHECK( m_fieldCondObjInputKey.initialize() );
 
    ATH_CHECK( m_inputSeedColletionKey.initialize() );
 
    ATH_CHECK( m_prdTruth.initialize(not m_prdTruth.empty()));
    ATH_CHECK( m_detEleCollKey.initialize(not m_detEleCollKey.empty()) );
 
    ATH_CHECK( m_EventInfoKey.initialize() );
 
    ATH_CHECK( m_paramEstimationTool.retrieve(EnableTool{!m_prdTruth.empty()}) );
    ATH_CHECK( m_trackingGeometrySvc.retrieve() );
    ATH_CHECK( m_geometryConvTool.retrieve(EnableTool{!m_prdTruth.empty()}) );
    ATH_CHECK( m_ctxProvider.initialize() );
 
    ATH_MSG_DEBUG("Monitoring settings ...");
    ATH_MSG_DEBUG(m_monGroupName);
 
    return AthMonitorAlgorithm::initialize();
  }

◆ inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

◆ isClonable()

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

inlineoverridevirtualinherited

Specify if the algorithm is clonable.

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

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

Definition at line 68 of file AthCommonAlgorithm.h.

                                           {
    return true;
  }

◆ matchParticleToSeedClusters()

void ActsTrk::SeedAnalysisAlg::matchParticleToSeedClusters	(	const PRD_MultiTruthCollection *	prdTruth,
		const Identifier &	id,
		std::map< int, int > &	countMap ) const

private

Definition at line 332 of file SeedAnalysisAlg.cxx.

                                                              {
    auto n1 = prdTruth->count(id);
    if (n1 == 0) {
      int bc = 0;
      auto nBC = countMap.count(bc);
      if (nBC == 0) {
    countMap[bc] = 1;
      } else {
    countMap[bc] += 1;
      }
    } else {
      using iprdt = PRD_MultiTruthCollection::const_iterator;
      std::pair<iprdt, iprdt> range = prdTruth->equal_range(id);
      for (iprdt itr = range.first; itr != range.second; ++itr) {
    auto bc = itr->second.barcode();
    auto nBC = countMap.count(bc);
    if (nBC == 0) {
      countMap[bc] = 1;
    } else {
      countMap[bc] += 1;
    }
      }
    }
  }

◆ 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.

◆ parseList()

StatusCode AthMonitorAlgorithm::parseList	(	const std::string &	line,
		std::vector< std::string > &	result ) const

virtualinherited

Parse a string into a vector.

The input string is a single long string of all of the trigger names. It parses this string and turns it into a vector, where each element is one trigger or trigger category.

Parameters

line	The input string.
result	The parsed output vector of strings.

Returns: StatusCode

Definition at line 350 of file AthMonitorAlgorithm.cxx.

                                                                                                   {
    std::string item;
    std::stringstream ss(line);
 
    ATH_MSG_DEBUG( "AthMonitorAlgorithm::parseList()" );
 
    while ( std::getline(ss, item, ',') ) {
        std::stringstream iss(item); // remove whitespace
        iss >> item;
        result.push_back(item);
    }
 
    return StatusCode::SUCCESS;
}

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ setFilterPassed()

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

inlinevirtualinherited

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

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

◆ sysExecute()

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

overridevirtualinherited

Execute an algorithm.

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

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

◆ sysInitialize()

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

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

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

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

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

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

Definition at line 60 of file AthCommonAlgorithm.cxx.

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

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

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

◆ trigChainsArePassed()

bool AthMonitorAlgorithm::trigChainsArePassed ( const std::vector< std::string > & vTrigNames ) const

inherited

Check whether triggers are passed.

For the event, use the trigger decision tool to check that at least one of the triggers listed in the supplied vector is passed.

Parameters

vTrigNames List of trigger names.

Returns: If empty input, default to true. If at least one trigger is specified, returns whether at least one trigger was passed.

Definition at line 208 of file AthMonitorAlgorithm.cxx.

                                                                                            {
 
  
  // If no triggers were given, return true.
  if (vTrigNames.empty()) return true;
  
  
  // Trigger: Check if this Algorithm is being run as an Express Stream job.
  // Events are entering the express stream are chosen randomly, and by chain,
  // Hence an additional check should be aplied to see if the chain(s)
  // monitored here are responsible for the event being selected for
  // the express stream.
 
  const auto group =  m_trigDecTool->getChainGroup(vTrigNames);
  if (m_enforceExpressTriggers){  
    const auto passedBits = m_trigDecTool->isPassedBits(group);
    bool expressPass = passedBits & TrigDefs::Express_passed; //bitwise AND
    if(!expressPass) {
      return false;
    }
  }
  
  // monitor the event if any of the chains in the chain group passes the event.
  return group->isPassed();
  
}

◆ 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_beamSpotKey

SG::ReadCondHandleKey< InDet::BeamSpotData > ActsTrk::SeedAnalysisAlg::m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"}

private

Definition at line 74 of file SeedAnalysisAlg.h.

74{this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"};

◆ m_ctxProvider

ContextUtility ActsTrk::SeedAnalysisAlg::m_ctxProvider {this}

private

Auxiliary class to access the calibration context.

TODO: Could also be used for magnetic field and geometry?

Definition at line 86 of file SeedAnalysisAlg.h.

86{this};

◆ m_dataType

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::m_dataType

protectedinherited

Instance of the DataType_t enum.

Definition at line 356 of file AthMonitorAlgorithm.h.

◆ m_dataTypeStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_dataTypeStr {this,"DataType","userDefined"}

protectedinherited

DataType string pulled from the job option and converted to enum.

Definition at line 358 of file AthMonitorAlgorithm.h.

358{this,"DataType","userDefined"};

◆ m_defaultLBDuration

Gaudi::Property<float> AthMonitorAlgorithm::m_defaultLBDuration {this,"DefaultLBDuration",60.}

protectedinherited

Default duration of one lumi block.

Definition at line 365 of file AthMonitorAlgorithm.h.

365{this,"DefaultLBDuration",60.};

◆ m_detailLevel

Gaudi::Property<int> AthMonitorAlgorithm::m_detailLevel {this,"DetailLevel",0}

protectedinherited

Sets the level of detail used in the monitoring.

Definition at line 366 of file AthMonitorAlgorithm.h.

366{this,"DetailLevel",0};

◆ m_detEleCollKey

SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollection > ActsTrk::SeedAnalysisAlg::m_detEleCollKey {this, "DetectorElements", "", "Key of input SiDetectorElementCollection"}

private

Definition at line 72 of file SeedAnalysisAlg.h.

72{this, "DetectorElements", "", "Key of input SiDetectorElementCollection"};

◆ 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_DQFilterTools

ToolHandleArray<IDQFilterTool> AthMonitorAlgorithm::m_DQFilterTools {this,"FilterTools",{}}

protectedinherited

Array of Data Quality filter tools.

Definition at line 346 of file AthMonitorAlgorithm.h.

346{this,"FilterTools",{}};

◆ m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy

privateinherited

Definition at line 374 of file AthMonitorAlgorithm.h.

◆ m_enforceExpressTriggers

Gaudi::Property<bool> AthMonitorAlgorithm::m_enforceExpressTriggers

privateinherited

Initial value:

{this,
                          "EnforceExpressTriggers", false,
                          "Requires that matched triggers made the event enter the express stream"}

Definition at line 377 of file AthMonitorAlgorithm.h.

                                               {this,
                          "EnforceExpressTriggers", false,
                          "Requires that matched triggers made the event enter the express stream"};

◆ m_environment

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::m_environment

protectedinherited

Instance of the Environment_t enum.

Definition at line 355 of file AthMonitorAlgorithm.h.

◆ m_environmentStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_environmentStr {this,"Environment","user"}

protectedinherited

Environment string pulled from the job option and converted to enum.

Definition at line 357 of file AthMonitorAlgorithm.h.

357{this,"Environment","user"};

◆ m_EventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> AthMonitorAlgorithm::m_EventInfoKey {this,"EventInfoKey","EventInfo"}

protectedinherited

Key for retrieving EventInfo from StoreGate.

Definition at line 367 of file AthMonitorAlgorithm.h.

367{this,"EventInfoKey","EventInfo"};

◆ m_evtStore

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

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

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

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_fieldCondObjInputKey

SG::ReadCondHandleKey< AtlasFieldCacheCondObj > ActsTrk::SeedAnalysisAlg::m_fieldCondObjInputKey

private

Initial value:

{this, "AtlasFieldCacheCondObj", "fieldCondObj",

"Name of the Magnetic Field conditions object key"}

Definition at line 75 of file SeedAnalysisAlg.h.

75 {this, "AtlasFieldCacheCondObj", "fieldCondObj",

76 "Name of the Magnetic Field conditions object key"};

◆ m_fileKey

Gaudi::Property<std::string> AthMonitorAlgorithm::m_fileKey {this,"FileKey",""}

protectedinherited

Internal Athena name for file.

Definition at line 363 of file AthMonitorAlgorithm.h.

363{this,"FileKey",""};

◆ m_geometryConvTool

PublicToolHandle<IGeometryRealmConvTool> ActsTrk::SeedAnalysisAlg::m_geometryConvTool {this, "GeometryRealmConvTool", ""}

private

Definition at line 66 of file SeedAnalysisAlg.h.

66{this, "GeometryRealmConvTool", ""};

◆ m_inputSeedColletionKey

SG::ReadHandleKey< ActsTrk::SeedContainer > ActsTrk::SeedAnalysisAlg::m_inputSeedColletionKey {this, "InputSeedCollection", "", ""}

private

Definition at line 68 of file SeedAnalysisAlg.h.

68{this, "InputSeedCollection", "", ""};

◆ m_lbDurationDataKey

SG::ReadCondHandleKey<LBDurationCondData> AthMonitorAlgorithm::m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}

protectedinherited

Definition at line 350 of file AthMonitorAlgorithm.h.

351{this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"};

◆ m_lumiDataKey

SG::ReadCondHandleKey<LuminosityCondData> AthMonitorAlgorithm::m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}

protectedinherited

Definition at line 348 of file AthMonitorAlgorithm.h.

349{this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"};

◆ m_monGroupName

Gaudi::Property< std::string > ActsTrk::SeedAnalysisAlg::m_monGroupName {this, "MonGroupName", "ActsSeedAnalysisAlg"}

private

Definition at line 78 of file SeedAnalysisAlg.h.

79{this, "MonGroupName", "ActsSeedAnalysisAlg"};

◆ m_name

std::string AthMonitorAlgorithm::m_name

privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

◆ m_paramEstimationTool

ToolHandle< ActsTrk::ITrackParamsEstimationTool > ActsTrk::SeedAnalysisAlg::m_paramEstimationTool {this, "TrackParamsEstimationTool", "", "Track Param Estimation from Seeds"}

private

Definition at line 64 of file SeedAnalysisAlg.h.

64{this, "TrackParamsEstimationTool", "", "Track Param Estimation from Seeds"};

◆ m_prdTruth

SG::ReadHandleKey< PRD_MultiTruthCollection > ActsTrk::SeedAnalysisAlg::m_prdTruth {this, "ITkClustersTruth", "", ""}

private

Definition at line 71 of file SeedAnalysisAlg.h.

71{this, "ITkClustersTruth", "", ""};

◆ m_toolLookupMap

std::unordered_map<std::string, size_t> AthMonitorAlgorithm::m_toolLookupMap

privateinherited

Definition at line 372 of file AthMonitorAlgorithm.h.

◆ m_tools

ToolHandleArray<GenericMonitoringTool> AthMonitorAlgorithm::m_tools {this,"GMTools",{}}

protectedinherited

Array of Generic Monitoring Tools.

Definition at line 341 of file AthMonitorAlgorithm.h.

341{this,"GMTools",{}};

◆ m_trackingGeometrySvc

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

private

Definition at line 69 of file SeedAnalysisAlg.h.

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

◆ m_trigDecTool

PublicToolHandle<Trig::TrigDecisionTool> AthMonitorAlgorithm::m_trigDecTool

protectedinherited

Tool to tell whether a specific trigger is passed.

Definition at line 345 of file AthMonitorAlgorithm.h.

◆ m_triggerChainString

Gaudi::Property<std::string> AthMonitorAlgorithm::m_triggerChainString {this,"TriggerChain",""}

protectedinherited

Trigger chain string pulled from the job option and parsed into a vector.

Definition at line 360 of file AthMonitorAlgorithm.h.

360{this,"TriggerChain",""};

◆ m_trigLiveFractionDataKey

SG::ReadCondHandleKey<TrigLiveFractionCondData> AthMonitorAlgorithm::m_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}

protectedinherited

Definition at line 352 of file AthMonitorAlgorithm.h.

353{this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"};

◆ m_useLumi

Gaudi::Property<bool> AthMonitorAlgorithm::m_useLumi {this,"EnableLumi",false}

protectedinherited

Allows use of various luminosity functions.

Definition at line 364 of file AthMonitorAlgorithm.h.

364{this,"EnableLumi",false};

◆ m_usePixel

Gaudi::Property< bool > ActsTrk::SeedAnalysisAlg::m_usePixel {this, "UsePixel", true, ""}

private

Definition at line 82 of file SeedAnalysisAlg.h.

82{this, "UsePixel", true, ""};

◆ m_useTopSp

Gaudi::Property< bool > ActsTrk::SeedAnalysisAlg::m_useTopSp {this, "useTopSp", false, "Use top SP. By default, use bottom SP."}

private

Definition at line 81 of file SeedAnalysisAlg.h.

81{this, "useTopSp", false, "Use top SP. By default, use bottom SP."};

◆ 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_vTrigChainNames

std::vector<std::string> AthMonitorAlgorithm::m_vTrigChainNames

protectedinherited

Vector of trigger chain names parsed from trigger chain string.

Definition at line 361 of file AthMonitorAlgorithm.h.

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

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ MonVarVec_t

◆ StoreGateSvc_t

Member Enumeration Documentation

◆ DataType_t

◆ Environment_t

Constructor & Destructor Documentation

◆ SeedAnalysisAlg()

◆ ~SeedAnalysisAlg()

Member Function Documentation

◆ dataType()

◆ dataTypeStringToEnum()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ environment()

◆ envStringToEnum()

◆ estimateParameters()

◆ evtStore()

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ fillHistograms()

◆ fillTruthHistograms()

◆ filterPassed()

◆ findSeedMajorityTruthParticle()

◆ GetEventInfo()

◆ getGroup()

◆ getTrigDecisionTool()

◆ identify() [1/2]

◆ identify() [2/2]

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ matchParticleToSeedClusters()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ parseList()

◆ renounce()

◆ renounceArray()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ trigChainsArePassed()

◆ updateVHKA()

Member Data Documentation

◆ m_beamSpotKey

◆ m_ctxProvider

◆ m_dataType

◆ m_dataTypeStr

◆ m_defaultLBDuration

◆ m_detailLevel

◆ m_detEleCollKey

◆ m_detStore

◆ m_DQFilterTools

◆ m_dummy

◆ m_enforceExpressTriggers

◆ m_environment

◆ m_environmentStr

◆ m_EventInfoKey

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_fieldCondObjInputKey

◆ m_fileKey

◆ m_geometryConvTool

◆ m_inputSeedColletionKey

◆ m_lbDurationDataKey

◆ m_lumiDataKey

◆ m_monGroupName

◆ m_name