![]() |
ATLAS Offline Software
|
#include <CombinatorialNSWSeedFinderAlg.h>
Public Member Functions | |
virtual | ~CombinatorialNSWSeedFinderAlg ()=default |
virtual StatusCode | initialize () override |
virtual StatusCode | execute (const EventContext &ctx) const override |
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, V, H > &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 |
Protected Member Functions | |
void | renounceArray (SG::VarHandleKeyArray &handlesArray) |
remove all handles from I/O resolution More... | |
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. More... | |
Private Types | |
enum | StripOrient { StripOrient::U, StripOrient::V, StripOrient::X, StripOrient::P, StripOrient::C, StripOrient::Unknown } |
Enumeration to classify the orientation of a NSW strip More... | |
enum | HitWindow { HitWindow::tooLow, HitWindow::inside, HitWindow::tooHigh } |
To fastly check whether a hit is roughly compatible with a muon trajectory a narrow corridor is opened from the estimated beamspot to the first tested hit in the seed finding. More... | |
using | HitVec = SpacePointPerLayerSplitter::HitVec |
using | HitLayVec = SpacePointPerLayerSplitter::HitLayVec |
using | HitLaySpan_t = std::vector< std::reference_wrapper< const HitVec > > |
Abbrivation of the space comprising multiple hit vectors without copy. More... | |
using | UsedHitMarker_t = std::vector< std::vector< char > > |
Abbrivation of the container book keeping whether a hit is used or not. More... | |
using | UsedHitSpan_t = std::vector< std::reference_wrapper< std::vector< char > >> |
Abbrivation of the container to pass a subset of markers wtihout copy. More... | |
using | InitialSeed_t = std::array< const SpacePoint *, 4 > |
Abbrivation of the More... | |
using | InitialSeedVec_t = std::vector< InitialSeed_t > |
Vector of initial seeds. More... | |
typedef ServiceHandle< StoreGateSvc > | StoreGateSvc_t |
Private Member Functions | |
StripOrient | classifyStrip (const SpacePoint &spacePoint) const |
Determines the orientation of the strip space point. More... | |
UsedHitMarker_t | emptyBookKeeper (const HitLayVec &sortedSp) const |
Constructs an empty HitMarker from the split space points. More... | |
HitWindow | hitFromIPCorridor (const SpacePoint &testHit, const Amg::Vector3D &beamSpotPos, const Amg::Vector3D &dirEstUp, const Amg::Vector3D &dirEstDn) const |
The hit is above the predefined corridor. More... | |
void | constructPrelimnarySeeds (const Amg::Vector3D &beamSpot, const HitLaySpan_t &combinatoricLayers, const UsedHitSpan_t &usedHits, InitialSeedVec_t &outVec) const |
Construct a set of prelimnary seeds from the selected combinatoric layers. More... | |
std::unique_ptr< SegmentSeed > | buildSegmentSeed (const InitialSeed_t &initialSeed, const AmgSymMatrix(2)&bMatrix, const HoughMaximum &max, const HitLaySpan_t &extensionLayers, const UsedHitSpan_t &usedHits) const |
Build the final seed from the initial seed hits and then attempt to append hits from the complementary layers onto the seed. More... | |
void | markHitsAsUsed (const SegmentSeed &seed, const HitLayVec &allSortHits, UsedHitMarker_t &usedHitMarker) const |
HitVec | extendHits (const Amg::Vector3D &startPos, const Amg::Vector3D &direction, const HitLaySpan_t &extensionLayers, const UsedHitSpan_t &usedHits) const |
std::vector< std::unique_ptr< SegmentSeed > > | findSeedsFromMaximum (const HoughMaximum &max, const ActsGeometryContext &gctx) const |
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &) |
specialization for handling Gaudi::Property<SG::VarHandleKey> More... | |
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &) |
specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More... | |
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &) |
specialization for handling Gaudi::Property<SG::VarHandleBase> More... | |
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &) |
specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More... | |
Private Attributes | |
SG::ReadHandleKey< EtaHoughMaxContainer > | m_etaKey {this, "CombinatorialReadKey", "MuonHoughNswMaxima"} |
SG::WriteHandleKey< SegmentSeedContainer > | m_writeKey {this, "CombinatorialPhiWriteKey", "MuonHoughNswSegmentSeeds"} |
SG::ReadHandleKey< ActsGeometryContext > | m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"} |
ServiceHandle< Muon::IMuonIdHelperSvc > | m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"} |
DoubleProperty | m_windowTheta {this, "thetaWindow", 0.5 * Gaudi::Units::deg} |
DoubleProperty | m_minPullThreshold {this, "maxPull", 5.} |
ToolHandle< MuonValR4::IPatternVisualizationTool > | m_visionTool {this, "VisualizationTool", ""} |
Pattern visualization tool. More... | |
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 28 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Abbrivation of the space comprising multiple hit vectors without copy.
Definition at line 54 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Definition at line 52 of file CombinatorialNSWSeedFinderAlg.h.
Definition at line 50 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Abbrivation of the
Definition at line 60 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Vector of initial seeds.
Definition at line 62 of file CombinatorialNSWSeedFinderAlg.h.
|
privateinherited |
Definition at line 388 of file AthCommonDataStore.h.
|
private |
Abbrivation of the container book keeping whether a hit is used or not.
Definition at line 56 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Abbrivation of the container to pass a subset of markers wtihout copy.
Definition at line 58 of file CombinatorialNSWSeedFinderAlg.h.
|
strongprivate |
To fastly check whether a hit is roughly compatible with a muon trajectory a narrow corridor is opened from the estimated beamspot to the first tested hit in the seed finding.
Hits in subsequent layers need to be within this corridor in order to be considered for seed construction. The HitWindow is the output classification of such a corridor test.
Enumerator | |
---|---|
tooLow | |
inside | The hit is below the predefined corridor. |
tooHigh | The hit is inside the defined window and hence an initial candidate. |
Definition at line 71 of file CombinatorialNSWSeedFinderAlg.h.
|
strongprivate |
Enumeration to classify the orientation of a NSW strip
Enumerator | |
---|---|
U | |
V | Stereo strips with positive angle. |
X | Stereo strips with negative angle. |
P | Ordinary eta strips. |
C | Single phi measurements. |
Unknown | Combined 2D space point (sTGC wire + strip / sTgc pad) |
Definition at line 39 of file CombinatorialNSWSeedFinderAlg.h.
|
virtualdefault |
|
private |
Build the final seed from the initial seed hits and then attempt to append hits from the complementary layers onto the seed.
initialSeed | Reference to the hit quadruplet that may form a seed |
bMatrix | Kernel matrix calculated from the layer configuration to construct the seed |
max | Refrence to the eta maximum from which the segment seed is constructed |
extensionLayers | Reference to the hits on the remaining layers of the detector |
usedExtensionHits | Refrence to the book keeper of which of the hits on the extension was already used |
Definition at line 284 of file CombinatorialNSWSeedFinderAlg.cxx.
|
overridevirtualinherited |
Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
Override this to return 0 for reentrant algorithms.
Override this to return 0 for reentrant algorithms.
Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.
|
private |
Determines the orientation of the strip space point.
Definition at line 85 of file CombinatorialNSWSeedFinderAlg.cxx.
|
private |
Construct a set of prelimnary seeds from the selected combinatoric layers.
Quadruplets of hits, one from each layer, are formed if they are all within the the corridor as described above
beamSpot | Position of the beam spot in the sector's frame |
combinatoricLayers | Quadruplet of four hit vectors from which the hits are retrieved |
usedHits | Mask marking hits that were already successfully added to a seed |
outVec | Reference to the output vector where the initial seeds are stored. The vector is cleared at the beginning and capacity is allocated accordingly |
Assign enough memory to the vector
The hit is alrady in a good seed. Don't consider again
Construct the beamspot to first hit connection to guestimate the angle
Apply cut window on theta of the seed.
Definition at line 182 of file CombinatorialNSWSeedFinderAlg.cxx.
|
inlineprivateinherited |
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
|
inlineprivateinherited |
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
|
inlineprivateinherited |
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
|
inlineprivateinherited |
specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>
Definition at line 199 of file AthCommonDataStore.h.
|
inlineinherited |
Declare a new Gaudi property.
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.
|
inlineinherited |
Declare a new Gaudi property.
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.
|
inlineinherited |
Definition at line 259 of file AthCommonDataStore.h.
|
inlineinherited |
Declare a new Gaudi property.
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.
|
inlineinherited |
Declare a new Gaudi property.
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.
|
inlineinherited |
Definition at line 145 of file AthCommonDataStore.h.
|
inlineinherited |
The standard StoreGateSvc/DetectorStore
Returns (kind of) a pointer to the StoreGateSvc
.
Definition at line 95 of file AthCommonDataStore.h.
|
private |
Constructs an empty HitMarker from the split space points.
sortedSp | List of space points sorted by layer |
Definition at line 76 of file CombinatorialNSWSeedFinderAlg.cxx.
|
inlineinherited |
The standard StoreGateSvc
(event store) Returns (kind of) a pointer to the StoreGateSvc
.
Definition at line 85 of file AthCommonDataStore.h.
|
inlineinherited |
The standard StoreGateSvc
(event store) Returns (kind of) a pointer to the StoreGateSvc
.
Definition at line 90 of file AthCommonDataStore.h.
|
overridevirtual |
Definition at line 432 of file CombinatorialNSWSeedFinderAlg.cxx.
|
private |
Definition at line 231 of file CombinatorialNSWSeedFinderAlg.cxx.
|
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
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.
|
inlinevirtualinherited |
Definition at line 96 of file AthCommonReentrantAlgorithm.h.
|
private |
Definition at line 322 of file CombinatorialNSWSeedFinderAlg.cxx.
|
inlineprivate |
The hit is above the predefined corridor.
Tests whether a hit is inside the corridor defined by line connecting the centre of the first candidate hit in the seed and the beam spot. The theta angle is varied by m_windowTheta to define the lower & upper direction etimate. The function tests whether the strip then crosses the corridor.
testHit | Reference to the hit to test |
beamSpotPos | Position of the beam spot serving as starting point |
dirEstUp | Direction vector defining the upper limit of the corridor |
dirEstDn | Direction vector defining the lower limit of the corridor |
Calculate the strip edges
Check whether the both edges are below the lower estimated muon arrival
Analogous check for the upper edge
No extrapolation needed
Analogous check for the upper edge
Definition at line 105 of file CombinatorialNSWSeedFinderAlg.cxx.
|
overridevirtual |
Definition at line 60 of file CombinatorialNSWSeedFinderAlg.cxx.
|
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.
|
overridevirtualinherited |
Specify if the algorithm is clonable.
Reentrant algorithms are clonable.
Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.
|
private |
That's ultra slow & should be revised
Definition at line 414 of file CombinatorialNSWSeedFinderAlg.cxx.
|
inlineinherited |
Definition at line 24 of file AthCommonMsg.h.
|
inlineinherited |
Definition at line 27 of file AthCommonMsg.h.
|
inlineinherited |
Definition at line 30 of file AthCommonMsg.h.
|
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.
|
inlineprotectedinherited |
Definition at line 380 of file AthCommonDataStore.h.
|
inlineprotectedinherited |
remove all handles from I/O resolution
Definition at line 364 of file AthCommonDataStore.h.
|
inlinevirtualinherited |
Definition at line 100 of file AthCommonReentrantAlgorithm.h.
|
overridevirtualinherited |
Execute an algorithm.
We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.
Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.
|
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 InputMakerBase, and HypoBase.
Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.
|
overridevirtualinherited |
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
|
inlineinherited |
Definition at line 308 of file AthCommonDataStore.h.
|
privateinherited |
Pointer to StoreGate (detector store by default)
Definition at line 393 of file AthCommonDataStore.h.
|
private |
Definition at line 120 of file CombinatorialNSWSeedFinderAlg.h.
|
privateinherited |
Pointer to StoreGate (event store by default)
Definition at line 390 of file AthCommonDataStore.h.
|
privateinherited |
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Empty if no symlinks were found.
Definition at line 114 of file AthCommonReentrantAlgorithm.h.
|
private |
Definition at line 126 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Definition at line 129 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Definition at line 140 of file CombinatorialNSWSeedFinderAlg.h.
|
privateinherited |
Definition at line 399 of file AthCommonDataStore.h.
|
privateinherited |
Definition at line 398 of file AthCommonDataStore.h.
|
private |
Pattern visualization tool.
Definition at line 143 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Definition at line 137 of file CombinatorialNSWSeedFinderAlg.h.
|
private |
Definition at line 123 of file CombinatorialNSWSeedFinderAlg.h.