Loading [MathJax]/extensions/tex2jax.js
 |
ATLAS Offline Software
|
#include <SeedAnalysisAlg.h>
|
| enum | Environment_t {
Environment_t::user = 0,
Environment_t::online,
Environment_t::tier0,
Environment_t::tier0Raw,
Environment_t::tier0ESD,
Environment_t::AOD,
Environment_t::altprod
} |
| | Specifies the processing environment. More...
|
| |
| enum | DataType_t {
DataType_t::userDefined = 0,
DataType_t::monteCarlo,
DataType_t::collisions,
DataType_t::cosmics,
DataType_t::heavyIonCollisions
} |
| | Specifies what type of input data is being monitored. More...
|
| |
|
| | SeedAnalysisAlg (const std::string &name, ISvcLocator *pSvcLocator) |
| |
| virtual | ~SeedAnalysisAlg () override=default |
| |
| virtual StatusCode | initialize () override |
| | initialize More...
|
| |
| virtual StatusCode | fillHistograms (const EventContext &ctx) const override |
| | adds event to the monitoring histograms More...
|
| |
| virtual StatusCode | execute (const EventContext &ctx) const override |
| | Applies filters and trigger requirements. More...
|
| |
| 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. More...
|
| |
| 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. More...
|
| |
| template<typename... T> |
| void | fill (const ToolHandle< GenericMonitoringTool > &groupHandle, T &&... variables) const |
| | Fills a variadic list of variables to a group by reference. More...
|
| |
| 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. More...
|
| |
| 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. More...
|
| |
| template<typename... T> |
| void | fill (const std::string &groupName, T &&... variables) const |
| | Fills a variadic list of variables to a group by name. More...
|
| |
| Environment_t | environment () const |
| | Accessor functions for the environment. More...
|
| |
| Environment_t | envStringToEnum (const std::string &str) const |
| | Convert the environment string from the python configuration to an enum object. More...
|
| |
| DataType_t | dataType () const |
| | Accessor functions for the data type. More...
|
| |
| DataType_t | dataTypeStringToEnum (const std::string &str) const |
| | Convert the data type string from the python configuration to an enum object. More...
|
| |
| const ToolHandle< GenericMonitoringTool > & | getGroup (const std::string &name) const |
| | Get a specific monitoring tool from the tool handle array. More...
|
| |
| const ToolHandle< Trig::TrigDecisionTool > & | getTrigDecisionTool () const |
| | Get the trigger decision tool member. More...
|
| |
| bool | trigChainsArePassed (const std::vector< std::string > &vTrigNames) const |
| | Check whether triggers are passed. More...
|
| |
| SG::ReadHandle< xAOD::EventInfo > | GetEventInfo (const EventContext &) const |
| | Return a ReadHandle for an EventInfo object (get run/event numbers, etc.) More...
|
| |
| virtual float | lbAverageInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate the average mu, i.e. More...
|
| |
| virtual float | lbInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate instantaneous number of interactions, i.e. More...
|
| |
| virtual float | lbAverageLuminosity (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate average luminosity (in ub-1 s-1 => 10^30 cm-2 s-1). More...
|
| |
| virtual float | lbLuminosityPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate the instantaneous luminosity per bunch crossing. More...
|
| |
| virtual double | lbDuration (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate the duration of the luminosity block (in seconds) More...
|
| |
| virtual float | lbAverageLivefraction (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate the average luminosity livefraction. More...
|
| |
| virtual float | livefractionPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate the live fraction per bunch crossing ID. More...
|
| |
| virtual double | lbLumiWeight (const EventContext &ctx=Gaudi::Hive::currentContext()) const |
| | Calculate the average integrated luminosity multiplied by the live fraction. More...
|
| |
| virtual StatusCode | parseList (const std::string &line, std::vector< std::string > &result) const |
| | Parse a string into a vector. More...
|
| |
| virtual StatusCode | sysInitialize () override |
| | Override sysInitialize. More...
|
| |
| virtual bool | isClonable () const override |
| | Specify if the algorithm is clonable. More...
|
| |
| virtual unsigned int | cardinality () const override |
| | Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
|
| |
| virtual StatusCode | sysExecute (const EventContext &ctx) override |
| | Execute an algorithm. More...
|
| |
| virtual const DataObjIDColl & | extraOutputDeps () const override |
| | Return the list of extra output dependencies. More...
|
| |
| 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. More...
|
| |
| const ServiceHandle< StoreGateSvc > & | evtStore () const |
| | The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
|
| |
| const ServiceHandle< StoreGateSvc > & | detStore () const |
| | The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
|
| |
| virtual StatusCode | sysStart () override |
| | Handle START transition. More...
|
| |
| virtual std::vector< Gaudi::DataHandle * > | inputHandles () const override |
| | Return this algorithm's input handles. More...
|
| |
| virtual std::vector< Gaudi::DataHandle * > | outputHandles () const override |
| | Return this algorithm's output handles. More...
|
| |
| Gaudi::Details::PropertyBase & | declareProperty (Gaudi::Property< T > &t) |
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &) |
| | Declare a new Gaudi property. More...
|
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &) |
| | Declare a new Gaudi property. More...
|
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &) |
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &) |
| | Declare a new Gaudi property. More...
|
| |
| Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc="none") |
| | Declare a new Gaudi property. More...
|
| |
| void | updateVHKA (Gaudi::Details::PropertyBase &) |
| |
| MsgStream & | msg () const |
| |
| MsgStream & | msg (const MSG::Level lvl) const |
| |
| bool | msgLvl (const MSG::Level lvl) const |
| |
|
| ToolHandleArray< GenericMonitoringTool > | m_tools {this,"GMTools",{}} |
| | Array of Generic Monitoring Tools. More...
|
| |
| PublicToolHandle< Trig::TrigDecisionTool > | m_trigDecTool {this, "TrigDecisionTool",""} |
| | Tool to tell whether a specific trigger is passed. More...
|
| |
| ToolHandleArray< IDQFilterTool > | m_DQFilterTools {this,"FilterTools",{}} |
| | Array of Data Quality filter tools. More...
|
| |
| 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. More...
|
| |
| AthMonitorAlgorithm::DataType_t | m_dataType |
| | Instance of the DataType_t enum. More...
|
| |
| Gaudi::Property< std::string > | m_environmentStr {this,"Environment","user"} |
| | Environment string pulled from the job option and converted to enum. More...
|
| |
| Gaudi::Property< std::string > | m_dataTypeStr {this,"DataType","userDefined"} |
| | DataType string pulled from the job option and converted to enum. More...
|
| |
| Gaudi::Property< std::string > | m_triggerChainString {this,"TriggerChain",""} |
| | Trigger chain string pulled from the job option and parsed into a vector. More...
|
| |
| std::vector< std::string > | m_vTrigChainNames |
| | Vector of trigger chain names parsed from trigger chain string. More...
|
| |
| Gaudi::Property< std::string > | m_fileKey {this,"FileKey",""} |
| | Internal Athena name for file. More...
|
| |
| Gaudi::Property< bool > | m_useLumi {this,"EnableLumi",false} |
| | Allows use of various luminosity functions. More...
|
| |
| Gaudi::Property< float > | m_defaultLBDuration {this,"DefaultLBDuration",60.} |
| | Default duration of one lumi block. More...
|
| |
| Gaudi::Property< int > | m_detailLevel {this,"DetailLevel",0} |
| | Sets the level of detail used in the monitoring. More...
|
| |
| SG::ReadHandleKey< xAOD::EventInfo > | m_EventInfoKey {this,"EventInfoKey","EventInfo"} |
| | Key for retrieving EventInfo from StoreGate. More...
|
| |
|
| 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 > &hndl, const SG::VarHandleKeyType &) |
| | specialization for handling Gaudi::Property<SG::VarHandleKey> More...
|
| |
| Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &) |
| | specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
|
| |
| Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &) |
| | specialization for handling Gaudi::Property<SG::VarHandleBase> More...
|
| |
| Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &) |
| | specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
|
| |
|
| ToolHandle< ActsTrk::ITrackParamsEstimationTool > | m_paramEstimationTool {this, "TrackParamsEstimationTool", "", "Track Param Estimation from Seeds"} |
| |
| PublicToolHandle< IActsTrackingGeometryTool > | m_trackingGeometryTool {this, "TrackingGeometryTool", "ActsTrackingGeometryTool"} |
| |
| ToolHandle< ActsTrk::IActsToTrkConverterTool > | m_ATLASConverterTool {this, "ATLASConverterTool", "ActsToTrkConverterTool", "Acts to ATLAS converter tool"} |
| |
| SG::ReadHandleKey< ActsTrk::SeedContainer > | m_inputSeedColletionKey {this, "InputSeedCollection", "", ""} |
| |
| 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, ""} |
| |
| 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. More...
|
| |
| StoreGateSvc_t | m_evtStore |
| | Pointer to StoreGate (event store by default) More...
|
| |
| StoreGateSvc_t | m_detStore |
| | Pointer to StoreGate (detector store by default) More...
|
| |
| std::vector< SG::VarHandleKeyArray * > | m_vhka |
| |
| bool | m_varHandleArraysDeclared |
| |
Definition at line 37 of file SeedAnalysisAlg.h.
◆ MonVarVec_t
◆ StoreGateSvc_t
◆ DataType_t
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 191 of file AthMonitorAlgorithm.h.
◆ Environment_t
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 172 of file AthMonitorAlgorithm.h.
◆ SeedAnalysisAlg()
| ActsTrk::SeedAnalysisAlg::SeedAnalysisAlg |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
pSvcLocator |
|
) |
| |
◆ ~SeedAnalysisAlg()
| virtual ActsTrk::SeedAnalysisAlg::~SeedAnalysisAlg |
( |
| ) |
|
|
overridevirtualdefault |
◆ cardinality()
Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
Override this to return 0 for reentrant algorithms.
Override this to return 0 for reentrant algorithms.
Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.
◆ 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 221 of file AthMonitorAlgorithm.h.
◆ dataTypeStringToEnum()
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 140 of file AthMonitorAlgorithm.cxx.
142 std::string lowerCaseStr =
str;
146 if( lowerCaseStr ==
"userdefined" ) {
148 }
else if( lowerCaseStr ==
"montecarlo" ) {
150 }
else if( lowerCaseStr ==
"collisions" ) {
152 }
else if( lowerCaseStr ==
"cosmics" ) {
154 }
else if( lowerCaseStr ==
"heavyioncollisions" ) {
157 ATH_MSG_WARNING(
"AthMonitorAlgorithm::dataTypeStringToEnum(): Unknown data type "
158 <<
str<<
", returning userDefined.");
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| hndl | Object holding the property value. |
| doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| hndl | Object holding the property value. |
| doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 221 of file AthCommonDataStore.h.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| property | Object holding the property value. |
| doc | Documentation string for the property. |
This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.
Definition at line 333 of file AthCommonDataStore.h.
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
| name | Name of the property. |
| property | Object holding the property value. |
| doc | Documentation string for the property. |
This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.
Definition at line 352 of file AthCommonDataStore.h.
◆ declareProperty() [6/6]
◆ detStore()
◆ environment()
Accessor functions for the environment.
- Returns
- the current value of the class's Environment_t instance.
Definition at line 205 of file AthMonitorAlgorithm.h.
◆ envStringToEnum()
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 112 of file AthMonitorAlgorithm.cxx.
114 std::string lowerCaseStr =
str;
118 if( lowerCaseStr ==
"user" ) {
120 }
else if( lowerCaseStr ==
"online" ) {
122 }
else if( lowerCaseStr ==
"tier0" ) {
124 }
else if( lowerCaseStr ==
"tier0raw" ) {
126 }
else if( lowerCaseStr ==
"tier0esd" ) {
128 }
else if( lowerCaseStr ==
"aod" ) {
130 }
else if( lowerCaseStr ==
"altprod" ) {
133 ATH_MSG_WARNING(
"AthMonitorAlgorithm::envStringToEnum(): Unknown environment "
134 <<
str<<
", returning user.");
◆ estimateParameters()
| std::array< float, 7 > ActsTrk::SeedAnalysisAlg::estimateParameters |
( |
const ActsTrk::Seed & |
seed, |
|
|
float |
pTPerHelixRadius |
|
) |
| const |
|
private |
Definition at line 380 of file SeedAnalysisAlg.cxx.
383 auto extractCoordinates =
386 std::array<float, 4> coordinates {
static_cast<float>(sp->x()),
387 static_cast<float>(sp->y()),
388 static_cast<float>(sp->z()),
389 static_cast<float>(std::sqrt(sp->x()*sp->x() + sp->y()*sp->y()))};
393 auto extractQuantities =
394 [] (
const std::array<float, 4>& sp,
395 const std::array<float, 4>& spM,
396 bool isBottom) -> std::array<float, 5>
398 auto& [xM, yM, zM, rM] = spM;
399 auto& [xO, yO, zO, rO] = sp;
401 float cosPhiM = xM / rM;
402 float sinPhiM = yM / rM;
409 float iDeltaR = std::sqrt(iDeltaR2);
411 float cot_theta =
deltaZ * iDeltaR * bottomFactor;
414 std::array<float, 5>
params =
426 const auto& sps =
seed.sp();
427 const auto* bottom = sps[0];
428 const auto* medium = sps[1];
429 const auto*
top = sps[2];
431 auto coo_b = extractCoordinates(bottom);
432 auto coo_m = extractCoordinates(medium);
433 auto coo_t = extractCoordinates(
top);
436 auto [cotThetaB, Zob, iDeltaRB, Ub, Vb] = extractQuantities(coo_b, coo_m,
true);
437 auto [cotThetaT, Zot, iDeltaRT, Ut, Vt] = extractQuantities(coo_t, coo_m,
false);
439 float squarediDeltaR2B = iDeltaRB*iDeltaRB;
440 float squarediDeltaR2T = iDeltaRB*iDeltaRT;
441 float squarediDeltaR =
std::min(squarediDeltaR2B, squarediDeltaR2T);
443 auto& [xB, yB, zB, rB] = coo_b;
444 auto& [xM, yM, zM, rM] = coo_m;
445 auto& [xT, yT, zT, rT] = coo_t;
453 float xb = dxb *
ax + dyb *
ay;
454 float yb = dyb *
ax - dxb *
ay;
455 float dxyb =
xb *
xb + yb * yb;
459 float xt = dxt *
ax + dyt *
ay;
460 float yt = dyt *
ax - dxt *
ay;
463 float tzb = dzb * std::sqrt( 1./dxyb );
464 float tzt = dzt * std::sqrt( 1./dxyt );
466 float sTzb2 = std::sqrt(1 + tzb*tzb);
470 return {-1, -1, -1, -1, -1, -1, -1};
473 float A = (Vt - Vb) / dU;
474 float S2 = 1. +
A *
A;
475 float B = Vb -
A * Ub;
477 if (B2 == 0) B2 = 1
e-8;
480 float dzdr_b = (zM - zB) / (rM - rB);
481 float dzdr_t = (zT - zM) / (rT - rM);
484 float cotThetaAvg2 = cotThetaB * cotThetaT;
485 if (cotThetaAvg2 <= 0) {
486 return {-1, -1, -1, -1, -1, -1, -1};
492 float pt = pTPerHelixRadius * std::sqrt(
S2 / B2) / 2.;
495 float d0 = std::abs((
A - B * rM) * rM);
500 float penalty = std::abs((tzb - tzt) / (squarediDeltaR * sTzb2));
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ 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 73 of file AthMonitorAlgorithm.cxx.
78 if (!filterItr->accept()) {
80 return StatusCode::SUCCESS;
87 return StatusCode::SUCCESS;
◆ extraDeps_update_handler()
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()
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.
◆ 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 49 of file SeedAnalysisAlg.cxx.
56 auto beamSpotData = beamSpotHandle.cptr();
65 Acts::MagneticFieldContext magFieldContext(fieldCondObj);
72 Acts::MagneticFieldProvider::Cache magFieldCache = magneticField.
makeCache( magFieldContext );
73 Acts::Vector3 bField = *magneticField.
getField( Acts::Vector3(beamPos.x(), beamPos.y(), 0),
90 [] (
const auto* seed) ->
double
91 {
return seed->sp()[0]->
x(); });
94 [] (
const auto* seed) ->
double
95 {
return seed->sp()[0]->
y(); });
98 [] (
const auto* seed) ->
double
99 {
return seed->sp()[0]->
z(); });
102 [] (
const auto* seed) ->
double
104 const auto* sp =
seed->sp()[0];
105 return std::sqrt(sp->x()*sp->x() + sp->y()*sp->y());
111 [] (
const auto* seed) ->
double
112 {
return seed->sp()[1]->
x(); });
115 [] (
const auto* seed) ->
double
116 {
return seed->sp()[1]->
y(); });
119 [] (
const auto* seed) ->
double
120 {
return seed->sp()[1]->
z(); });
123 [] (
const auto* seed) ->
double
125 const auto* sp =
seed->sp()[1];
126 return std::sqrt(sp->x()*sp->x() + sp->y()*sp->y());
132 [] (
const auto* seed) ->
double
133 {
return seed->sp()[2]->
x(); });
136 [] (
const auto* seed) ->
double
137 {
return seed->sp()[2]->
y(); });
140 [] (
const auto* seed) ->
double
141 {
return seed->sp()[2]->
z(); });
144 [] (
const auto* seed) ->
double
146 const auto* sp =
seed->sp()[2];
147 return std::sqrt(sp->x()*sp->x() + sp->y()*sp->y());
150 std::vector< std::array<float, 7> > parametersCollection;
151 parametersCollection.reserve(seed_collection->
size());
153 for (
const auto* seed : *seed_collection) {
154 parametersCollection.push_back(
estimateParameters(*seed, 300. * bField[2] / 1000.) );
158 [] (
const auto&
params) ->
float
161 [] (
const auto&
params) ->
float
164 [] (
const auto&
params) ->
float
167 [] (
const auto&
params) ->
float
171 [] (
const auto&
params) ->
float
174 [] (
const auto&
params) ->
float
179 [] (
const auto&
params) ->
float
189 std::vector<int> vec_truthBarcode;
190 std::vector<double> vec_truthProb;
197 monitor_x1, monitor_y1, monitor_z1, monitor_r1,
198 monitor_x2, monitor_y2, monitor_z2, monitor_r2,
199 monitor_x3, monitor_y3, monitor_z3, monitor_r3,
200 monitor_param_pt, monitor_param_theta, monitor_param_eta, monitor_param_d0,
201 monitor_param_dzdr_b, monitor_param_dzdr_t,
202 monitor_param_penalty,
203 monitor_event_number, monitor_actual_mu,
204 monitor_truth_barcode, monitor_truth_prob);
206 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 209 of file SeedAnalysisAlg.cxx.
231 Acts::MagneticFieldContext magFieldContext(fieldCondObj);
237 auto retrieveSurfaceFunction =
238 [
this, &detElements] (
const ActsTrk::Seed&
seed,
bool useTopSp) ->
const Acts::Surface&
249 std::vector<bool> vec_pass;
250 vec_pass.reserve(seed_container.
size());
252 std::vector<double> estimated_pt;
253 std::vector<double> estimated_eta;
254 estimated_pt.reserve(seed_container.
size());
255 estimated_eta.reserve(seed_container.
size());
257 for (
const auto* seed : seed_container) {
258 std::optional<Acts::BoundTrackParameters> optTrackParams =
262 geo_context.context(),
264 retrieveSurfaceFunction);
266 if ( not optTrackParams.has_value() )
continue;
268 const auto param = optTrackParams.value();
269 estimated_pt.push_back( param.transverseMomentum() );
270 estimated_eta.push_back( -
std::log(
std::tan(0.5 * param.parameters()[Acts::eBoundTheta]) ) );
272 std::map<int, int> truthHits;
274 const auto& sps =
seed->sp();
275 for (
const auto* sp : sps) {
277 for (
int cluster_number(0); cluster_number < number_of_clusters; cluster_number++) {
279 const auto*
cluster = els[cluster_number];
289 truthBarCodeVec.push_back(
barcode);
290 truthProbVec.push_back(
prob);
300 monitor_estimated_pt, monitor_estimated_eta);
302 return StatusCode::SUCCESS;
◆ filterPassed()
◆ findSeedMajorityTruthParticle()
| std::pair< int, double > ActsTrk::SeedAnalysisAlg::findSeedMajorityTruthParticle |
( |
const std::map< int, int > & |
countMap | ) |
const |
|
private |
◆ GetEventInfo()
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 107 of file AthMonitorAlgorithm.cxx.
◆ getGroup()
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 164 of file AthMonitorAlgorithm.cxx.
176 if (!isInitialized()) {
178 "It seems that the AthMonitorAlgorithm::initialize was not called "
179 "in derived class initialize method");
183 [](
const std::string&
s,
auto h) { return s +
"," + h.first; });
184 ATH_MSG_FATAL(
"The tool " <<
name <<
" could not be found in the tool array of the "
185 <<
"monitoring algorithm " <<
m_name <<
". This probably reflects a discrepancy between "
186 <<
"your python configuration and c++ filling code. Note: your available groups are {"
187 << available <<
"}.");
◆ getTrigDecisionTool()
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 194 of file AthMonitorAlgorithm.cxx.
◆ identify() [1/2]
Definition at line 305 of file SeedAnalysisAlg.cxx.
311 if (not pixelLinkAcc.isAvailable (cluster))
315 return (*pixelLink)->identify();
◆ identify() [2/2]
Definition at line 318 of file SeedAnalysisAlg.cxx.
324 if (not stripLinkAcc.isAvailable (cluster))
328 return (*stripLink)->identify();
◆ initialize()
| StatusCode ActsTrk::SeedAnalysisAlg::initialize |
( |
| ) |
|
|
overridevirtual |
◆ inputHandles()
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()
◆ matchParticleToSeedClusters()
Definition at line 332 of file SeedAnalysisAlg.cxx.
335 auto n1 = prdTruth->count(
id);
338 auto nBC = countMap.count(bc);
345 using iprdt = PRD_MultiTruthCollection::const_iterator;
346 std::pair<iprdt, iprdt>
range = prdTruth->equal_range(
id);
347 for (iprdt itr =
range.first; itr !=
range.second; ++itr) {
348 auto bc = itr->second.barcode();
349 auto nBC = countMap.count(bc);
◆ msg() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ outputHandles()
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 341 of file AthMonitorAlgorithm.cxx.
343 std::stringstream
ss(
line);
347 while ( std::getline(
ss,
item,
',') ) {
348 std::stringstream iss(
item);
353 return StatusCode::SUCCESS;
◆ renounce()
◆ renounceArray()
◆ setFilterPassed()
◆ sysExecute()
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.
◆ sysInitialize()
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 InputMakerBase, and HypoBase.
Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.
110 if (
sc.isFailure()) {
118 if ( cs.retrieve().isFailure() ) {
120 return StatusCode::SUCCESS;
122 if (cs->regHandle(
this,*
h).isFailure()) {
123 sc = StatusCode::FAILURE;
124 ATH_MSG_ERROR(
"unable to register WriteCondHandle " <<
h->fullKey()
◆ sysStart()
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 199 of file AthMonitorAlgorithm.cxx.
203 if (vTrigNames.empty())
return true;
215 bool expressPass = passedBits & TrigDefs::Express_passed;
222 return group->isPassed();
◆ updateVHKA()
◆ m_ATLASConverterTool
◆ m_beamSpotKey
◆ m_dataType
◆ m_dataTypeStr
| Gaudi::Property<std::string> AthMonitorAlgorithm::m_dataTypeStr {this,"DataType","userDefined"} |
|
protectedinherited |
◆ m_defaultLBDuration
| Gaudi::Property<float> AthMonitorAlgorithm::m_defaultLBDuration {this,"DefaultLBDuration",60.} |
|
protectedinherited |
◆ m_detailLevel
| Gaudi::Property<int> AthMonitorAlgorithm::m_detailLevel {this,"DetailLevel",0} |
|
protectedinherited |
◆ m_detEleCollKey
◆ m_detStore
◆ m_DQFilterTools
| ToolHandleArray<IDQFilterTool> AthMonitorAlgorithm::m_DQFilterTools {this,"FilterTools",{}} |
|
protectedinherited |
◆ m_dummy
◆ 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 372 of file AthMonitorAlgorithm.h.
◆ m_environment
◆ 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 352 of file AthMonitorAlgorithm.h.
◆ m_EventInfoKey
◆ m_evtStore
◆ m_extendedExtraObjects
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_fieldCondObjInputKey
Initial value:{this, "AtlasFieldCacheCondObj", "fieldCondObj",
"Name of the Magnetic Field conditions object key"}
Definition at line 75 of file SeedAnalysisAlg.h.
◆ m_fileKey
| Gaudi::Property<std::string> AthMonitorAlgorithm::m_fileKey {this,"FileKey",""} |
|
protectedinherited |
◆ m_inputSeedColletionKey
◆ m_lbDurationDataKey
◆ m_lumiDataKey
◆ m_monGroupName
| Gaudi::Property< std::string > ActsTrk::SeedAnalysisAlg::m_monGroupName {this, "MonGroupName", "ActsSeedAnalysisAlg"} |
|
private |
◆ m_name
| std::string AthMonitorAlgorithm::m_name |
|
privateinherited |
◆ m_paramEstimationTool
◆ m_prdTruth
◆ m_toolLookupMap
| std::unordered_map<std::string, size_t> AthMonitorAlgorithm::m_toolLookupMap |
|
privateinherited |
◆ m_tools
◆ m_trackingGeometryTool
◆ m_trigDecTool
◆ 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 355 of file AthMonitorAlgorithm.h.
◆ m_trigLiveFractionDataKey
◆ m_useLumi
| Gaudi::Property<bool> AthMonitorAlgorithm::m_useLumi {this,"EnableLumi",false} |
|
protectedinherited |
◆ m_usePixel
| Gaudi::Property< bool > ActsTrk::SeedAnalysisAlg::m_usePixel {this, "UsePixel", true, ""} |
|
private |
◆ m_useTopSp
| Gaudi::Property< bool > ActsTrk::SeedAnalysisAlg::m_useTopSp {this, "useTopSp", false, "Use top SP. By default, use bottom SP."} |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vhka
◆ m_vTrigChainNames
| std::vector<std::string> AthMonitorAlgorithm::m_vTrigChainNames |
|
protectedinherited |
The documentation for this class was generated from the following files:
SG::ReadHandleKey< PRD_MultiTruthCollection > m_prdTruth
PublicToolHandle< Trig::TrigDecisionTool > m_trigDecTool
Tool to tell whether a specific trigger is passed.
SG::ReadCondHandleKey< InDetDD::SiDetectorElementCollection > m_detEleCollKey
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
A PRD is mapped onto all contributing particles.
bool accumulate(AccumulateMap &map, std::vector< module_t > const &modules, FPGATrackSimMatrixAccumulator const &acc)
Accumulates an accumulator (e.g.
std::vector< std::string > m_vTrigChainNames
Vector of trigger chain names parsed from trigger chain string.
StatusCode fillTruthHistograms(const EventContext &ctx, const ActsTrk::SeedContainer &seed_container, std::vector< int > &, std::vector< double > &) const
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)
SG::ReadHandleKey< xAOD::EventInfo > m_EventInfoKey
Key for retrieving EventInfo from StoreGate.
Trk::Surface & surface()
Element Surface.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
std::unordered_map< std::string, size_t > m_toolLookupMap
AthMonitorAlgorithm::Environment_t m_environment
Instance of the Environment_t enum.
ICaloAffectedTool is abstract interface for tools checking if 4 mom is in calo affected region.
Helper class to provide constant type-safe access to aux data.
bool const RAWDATA *ch2 const
bool empty() const
Test if the key is blank.
bool trigChainsArePassed(const std::vector< std::string > &vTrigNames) const
Check whether triggers are passed.
MagneticFieldProvider::Cache makeCache(const Acts::MagneticFieldContext &mctx) const override
ValuesCollection< T > Collection(std::string name, const T &collection)
Declare a monitored (double-convertible) collection.
virtual void setOwner(IDataHandleHolder *o)=0
const ToolHandle< GenericMonitoringTool > m_dummy
const std::vector< DetectorIDHashType > & elementIdList() const
Returns the IdentifierHash of the spacepoint (corresponds to the detector element IdentifierHash)
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
std::pair< int, double > findSeedMajorityTruthParticle(const std::map< int, int > &countMap) const
Gaudi::Property< bool > m_useTopSp
virtual const DataObjIDColl & extraOutputDeps() const override
Return the list of extra output dependencies.
StripCluster_v1 StripCluster
Define the version of the strip cluster class.
ToolHandleArray< IDQFilterTool > m_DQFilterTools
Array of Data Quality filter tools.
ToolHandle< ActsTrk::ITrackParamsEstimationTool > m_paramEstimationTool
bool isBottom(const T &p)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
::StatusCode StatusCode
StatusCode definition for legacy code.
const_pointer_type get() const
Dereference the pointer, but don't cache anything.
void tolower(std::string &s)
virtual StatusCode fillHistograms(const EventContext &ctx) const =0
adds event to the monitoring histograms
virtual StatusCode sysExecute(const EventContext &ctx) override
Execute an algorithm.
ToolHandle< ActsTrk::IActsToTrkConverterTool > m_ATLASConverterTool
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
SG::ReadCondHandleKey< AtlasFieldCacheCondObj > m_fieldCondObjInputKey
SG::ReadCondHandleKey< InDet::BeamSpotData > m_beamSpotKey
Acts::Seed< xAOD::SpacePoint, 3ul > Seed
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
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.
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
const std::vector< const xAOD::UncalibratedMeasurement * > & measurements() const
Returns the index of the measurements.
Gaudi::Property< std::string > m_monGroupName
void matchParticleToSeedClusters(const PRD_MultiTruthCollection *prdTruth, const Identifier &id, std::map< int, int > &countMap) const
virtual bool isValid() override final
Can the handle be successfully dereferenced?
AthMonitorAlgorithm::DataType_t m_dataType
Instance of the DataType_t enum.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
ElementLink implementation for ROOT usage.
ToolHandleArray< GenericMonitoringTool > m_tools
Array of Generic Monitoring Tools.
std::array< float, 7 > estimateParameters(const ActsTrk::Seed &seed, float pTPerHelixRadius) const
Gaudi::Property< bool > m_enforceExpressTriggers
StatusCode initialize(bool used=true)
virtual StatusCode initialize() override
initialize
Gaudi::Property< bool > m_usePixel
Acts::Result< Acts::Vector3 > getField(const Acts::Vector3 &position, Acts::MagneticFieldProvider::Cache &gcache) const override
const Identifier identify(const xAOD::PixelCluster &) const
#define ATH_MSG_WARNING(x)
UncalibMeasType
Define the type of the uncalibrated measurement.
struct TBPatternUnitContext S2
AthMonitorAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
virtual StatusCode sysInitialize() override
Override sysInitialize.
Declare a monitored scalar variable.
SG::ReadHandleKey< ActsTrk::SeedContainer > m_inputSeedColletionKey
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
size_type size() const noexcept
Returns the number of elements in the collection.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
const SiDetectorElement * getDetectorElement(const IdentifierHash &hash) const
PublicToolHandle< IActsTrackingGeometryTool > m_trackingGeometryTool
1.8.18