|
ATLAS Offline Software
|
#include <AnalysisMuonThinningAlg.h>
|
| AnalysisMuonThinningAlg (const std::string &n, ISvcLocator *p) |
|
| ~AnalysisMuonThinningAlg ()=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 > &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 |
|
|
ToolHandle< CP::IMuonSelectionTool > | m_muonSelTool {this, "SelectionTool", "", "Configured instance of the MuonSelectionTool"} |
|
Gaudi::Property< int > | m_quality |
|
SG::ReadDecorHandleKeyArray< xAOD::MuonContainer > | m_MuonPassKeys |
|
SG::ReadDecorHandleKeyArray< xAOD::TrackParticleContainer > | m_TrkPassKeys |
|
Gaudi::Property< std::string > | m_streamName {this, "StreamName", "", "Name of the stream being thinned"} |
| Thinning of unneeded muon tracks. More...
|
|
SG::ThinningHandleKey< xAOD::MuonContainer > | m_muonKey {this, "MuonThinning", "Muons", "Name of the muon container behind"} |
|
SG::ThinningHandleKey< xAOD::TrackParticleContainer > | m_IdTrkKey |
|
SG::ThinningHandleKey< xAOD::TrackParticleContainer > | m_FwdIdTrkKey |
|
SG::ThinningHandleKey< xAOD::TrackParticleContainer > | m_MSTrkKey |
|
SG::ThinningHandleKey< xAOD::TrackParticleContainer > | m_METrkKey |
|
SG::ThinningHandleKey< xAOD::TrackParticleContainer > | m_MSOETrkKey |
|
SG::ThinningHandleKey< xAOD::TrackParticleContainer > | m_CmbTrkKey |
|
SG::ThinningHandleKey< xAOD::MuonSegmentContainer > | m_SegmentKey |
|
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 16 of file AnalysisMuonThinningAlg.h.
◆ StoreGateSvc_t
◆ AnalysisMuonThinningAlg()
DerivationFramework::AnalysisMuonThinningAlg::AnalysisMuonThinningAlg |
( |
const std::string & |
n, |
|
|
ISvcLocator * |
p |
|
) |
| |
◆ ~AnalysisMuonThinningAlg()
DerivationFramework::AnalysisMuonThinningAlg::~AnalysisMuonThinningAlg |
( |
| ) |
|
|
default |
◆ cardinality()
unsigned int AthReentrantAlgorithm::cardinality |
( |
| ) |
const |
|
overridevirtualinherited |
Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
Override this to return 0 for reentrant algorithms.
Override this to return 0 for reentrant algorithms.
Definition at line 55 of file AthReentrantAlgorithm.cxx.
◆ 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()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode DerivationFramework::AnalysisMuonThinningAlg::execute |
( |
const EventContext & |
ctx | ) |
const |
|
overridevirtual |
The Muon container thinning
Initialize the decorators whitelisting the muon to be written to the output
Cache the thinning decisions for each associated track particle
Lambda function to fill the thinning helpers above
Do not forget to apply the thinning to the SegmentContainers as well
Loop over the muon container
The muon is rejected by the selection tool & not marked as pass by the augmentation tools upstream
Ensure that all associated track particles are written to the output file
Proceed with the Segment containers
The next step is to find all track particles that are accepted by the upstream augmentation tools as well
The vector has been freshly created -> the container is unknown to the current tool
Put the tracks on the white list
Propagate the decisions
Definition at line 61 of file AnalysisMuonThinningAlg.cxx.
66 return StatusCode::FAILURE;
71 std::vector<MuonPassDecor> mu_passFlags;
75 KeepMap thin_decisions;
76 std::vector<std::unique_ptr<TrkThinning>> TrkContainers;
79 auto make_handle = [&,
this](
const TrkThinKey&
key) {
82 return StatusCode::SUCCESS;
85 TrkContainers.emplace_back(std::make_unique<TrkThinning>(
key, ctx));
87 std::unique_ptr<TrkThinning>& last = TrkContainers.back();
88 if (!last->isValid()) {
90 return StatusCode::FAILURE;
93 thin_decisions[cont].
resize(cont->
size(),
false);
94 return StatusCode::SUCCESS;
105 std::unique_ptr<SegmentThinning> segmentContainer =
107 std::vector<bool> keep_seg(segmentContainer ? segmentContainer->cptr()->size() : 0,
false);
114 <<
" quality: " <<
muon->quality() <<
", author: " <<
muon->author());
117 std::find_if(mu_passFlags.begin(), mu_passFlags.end(),
118 [
this,
muon](
const MuonPassDecor& decor) ->bool {
119 ATH_MSG_VERBOSE(
"Check decorator "<<decor.decorKey()<<
" "
120 <<SG::AuxTypeRegistry::instance().getName(decor.auxid())
121 <<
" "<< decor(*muon));
123 }) == mu_passFlags.end())
125 keep_muo[
muon->index()] =
true;
130 {xAOD::Muon::InnerDetectorTrackParticle, xAOD::Muon::MuonSpectrometerTrackParticle, xAOD::Muon::CombinedTrackParticle,
131 xAOD::Muon::ExtrapolatedMuonSpectrometerTrackParticle, xAOD::Muon::MSOnlyExtrapolatedMuonSpectrometerTrackParticle}) {
137 KeepMap::iterator itr = std::find_if(thin_decisions.begin(), thin_decisions.end(),
138 [
track](
const KeepPair& pair) { return pair.first == track->container(); });
139 if (itr == thin_decisions.end()) {
141 <<
" a valid associated container");
144 itr->second[
track->index()] =
true;
147 if (segmentContainer) {
148 for (
size_t s = 0;
s <
muon->nMuonSegments(); ++
s) {
150 if (seg) keep_seg[seg->
index()] =
true;
157 TrackPassDecor handle{
key, ctx};
160 std::vector<bool>& trk_thin = thin_decisions[trks];
163 if (!trks->
empty() && trk_thin.empty()) {
164 ATH_MSG_WARNING(
"The container " <<
key.fullKey() <<
" is not part of the current thinning scheme. Please check");
169 if (handle(*trk)) trk_thin[trk->index()] =
true;
175 for (std::unique_ptr<TrkThinning>& thin : TrkContainers) {
176 std::vector<bool>& dec = thin_decisions[thin->cptr()];
179 if (segmentContainer) { segmentContainer->keep(keep_seg); }
181 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()
const DataObjIDColl & AthReentrantAlgorithm::extraOutputDeps |
( |
| ) |
const |
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 79 of file AthReentrantAlgorithm.cxx.
86 return Algorithm::extraOutputDeps();
◆ filterPassed()
virtual bool AthReentrantAlgorithm::filterPassed |
( |
const EventContext & |
ctx | ) |
const |
|
inlinevirtualinherited |
◆ initialize()
StatusCode DerivationFramework::AnalysisMuonThinningAlg::initialize |
( |
| ) |
|
|
overridevirtual |
Definition at line 25 of file AnalysisMuonThinningAlg.cxx.
29 return StatusCode::FAILURE;
44 std::stringstream sstr{};
46 sstr<<
" *** "<<
key.fullKey()<<std::endl;
48 ATH_MSG_INFO(
"Accept muons with upstream flag "<<std::endl<<sstr.str());
51 std::stringstream sstr{};
53 sstr<<
" *** "<<
key.fullKey()<<std::endl;
55 ATH_MSG_INFO(
"Accept tracks with upstream flag "<<std::endl<<sstr.str());
58 return StatusCode::SUCCESS;
◆ 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()
bool AthReentrantAlgorithm::isClonable |
( |
| ) |
const |
|
overridevirtualinherited |
Specify if the algorithm is clonable.
Reentrant algorithms are clonable.
Reimplemented in Simulation::BeamEffectsAlg, InDet::SiTrackerSpacePointFinder, InDet::SCT_Clusterization, InDet::SiSPSeededTrackFinder, SCTRawDataProvider, InDet::GNNSeedingTrackMaker, SCT_PrepDataToxAOD, RoIBResultToxAOD, SCT_CablingCondAlgFromCoraCool, SCT_ReadCalibDataTestAlg, SCT_CablingCondAlgFromText, InDet::SiSPGNNTrackMaker, SCT_ReadCalibChipDataTestAlg, SCT_TestCablingAlg, SCT_ConfigurationConditionsTestAlg, ITkPixelCablingAlg, ITkStripCablingAlg, SCTEventFlagWriter, SCT_ConditionsSummaryTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_TdaqEnabledTestAlg, SCT_SiliconConditionsTestAlg, SCTSiLorentzAngleTestAlg, SCT_ByteStreamErrorsTestAlg, SCT_ConditionsParameterTestAlg, SCT_FlaggedConditionTestAlg, SCT_StripVetoTestAlg, SCT_RawDataToxAOD, and SCTSiPropertiesTestAlg.
Definition at line 44 of file AthReentrantAlgorithm.cxx.
◆ 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.
◆ renounce()
◆ renounceArray()
◆ setFilterPassed()
virtual void AthReentrantAlgorithm::setFilterPassed |
( |
bool |
state, |
|
|
const EventContext & |
ctx |
|
) |
| const |
|
inlinevirtualinherited |
◆ sysExecute()
StatusCode AthReentrantAlgorithm::sysExecute |
( |
const EventContext & |
ctx | ) |
|
|
overridevirtualinherited |
Execute an algorithm.
We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.
Definition at line 67 of file AthReentrantAlgorithm.cxx.
69 return Gaudi::Algorithm::sysExecute (ctx);
◆ sysInitialize()
StatusCode AthReentrantAlgorithm::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 InputMakerBase, and HypoBase.
Definition at line 96 of file AthReentrantAlgorithm.cxx.
107 if ( cs.retrieve().isFailure() ) {
109 return StatusCode::SUCCESS;
111 if (cs->regHandle(
this,*
h).isFailure()) {
112 sc = StatusCode::FAILURE;
113 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.
◆ updateVHKA()
◆ m_CmbTrkKey
Initial value:{this, "CmbTrkThinning", "CombinedMuonTrackParticles",
"Key to get rid of the unneeded MSOE tracks"}
Definition at line 52 of file AnalysisMuonThinningAlg.h.
◆ m_detStore
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthReentrantAlgorithm::m_extendedExtraObjects |
|
privateinherited |
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Empty if no symlinks were found.
Definition at line 153 of file AthReentrantAlgorithm.h.
◆ m_FwdIdTrkKey
Initial value:{
this, "IdTrkFwdThinning", "InDetForwardTrackParticles",
"Thin the unneeded forward ID tracks (Needed for Sillicon associated Forward)"}
Definition at line 43 of file AnalysisMuonThinningAlg.h.
◆ m_IdTrkKey
Initial value:{this, "IdTrkThinnig", "InDetTrackParticles",
"Thin the unneeded ID track particles associated with muons"}
Definition at line 41 of file AnalysisMuonThinningAlg.h.
◆ m_METrkKey
Initial value:{this, "METrkThinning", "ExtrapolatedMuonTrackParticles",
"Key to get rid of the unneeded ME tracks"}
Definition at line 48 of file AnalysisMuonThinningAlg.h.
◆ m_MSOETrkKey
Initial value:{this, "MSOETrkThinning", "MSOnlyExtrapolatedMuonTrackParticles",
"Key to get rid of the unneeded MSOE tracks"}
Definition at line 50 of file AnalysisMuonThinningAlg.h.
◆ m_MSTrkKey
Initial value:{this, "MSTrkThinning", "MuonSpectrometerTrackParticles",
"Key to get rid of the unneeded MS tracks"}
Definition at line 46 of file AnalysisMuonThinningAlg.h.
◆ m_muonKey
◆ m_MuonPassKeys
Initial value:{
this, "MuonPassFlags", {}, "Decorators to safe analysis muons needed for MCP studies "}
Definition at line 31 of file AnalysisMuonThinningAlg.h.
◆ m_muonSelTool
ToolHandle<CP::IMuonSelectionTool> DerivationFramework::AnalysisMuonThinningAlg::m_muonSelTool {this, "SelectionTool", "", "Configured instance of the MuonSelectionTool"} |
|
private |
◆ m_quality
Gaudi::Property<int> DerivationFramework::AnalysisMuonThinningAlg::m_quality |
|
private |
Initial value:
"Minimum working point that the muon has to satisfy in order to be survive"}
Definition at line 28 of file AnalysisMuonThinningAlg.h.
◆ m_SegmentKey
Initial value:{this, "SegmentThinning", "MuonSegments",
"Key to get rid of the unneeded segments"}
Definition at line 54 of file AnalysisMuonThinningAlg.h.
◆ m_streamName
Gaudi::Property<std::string> DerivationFramework::AnalysisMuonThinningAlg::m_streamName {this, "StreamName", "", "Name of the stream being thinned"} |
|
private |
◆ m_TrkPassKeys
Initial value:{
this, "TrkPassFlags", {}, "Decorator to safe tracks needed for MCP studies"}
Definition at line 34 of file AnalysisMuonThinningAlg.h.
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
JetConstituentVector::iterator iterator
SG::ThinningHandleKey< xAOD::TrackParticleContainer > m_MSOETrkKey
Gaudi::Property< std::string > m_streamName
Thinning of unneeded muon tracks.
SG::ThinningHandleKey< xAOD::MuonContainer > m_muonKey
SG::ThinningHandleKey< xAOD::TrackParticleContainer > m_CmbTrkKey
MuonContainer_v1 MuonContainer
Definition of the current "Muon container version".
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
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)
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
Class describing a MuonSegment.
SG::ThinningHandleKey< xAOD::TrackParticleContainer > m_IdTrkKey
SG::ReadDecorHandleKeyArray< xAOD::TrackParticleContainer > m_TrkPassKeys
Handle for requesting thinning for a data object.
virtual void setOwner(IDataHandleHolder *o)=0
SG::ThinningHandleKey< xAOD::MuonSegmentContainer > m_SegmentKey
SG::ReadDecorHandleKeyArray< xAOD::MuonContainer > m_MuonPassKeys
AthReentrantAlgorithm()
Default constructor:
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
ToolHandle< CP::IMuonSelectionTool > m_muonSelTool
::StatusCode StatusCode
StatusCode definition for legacy code.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
size_t index() const
Return the index of this element within its container.
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
void resize(size_type sz)
Resizes the collection to the specified number of elements.
SG::ThinningHandleKey< xAOD::TrackParticleContainer > m_METrkKey
virtual StatusCode sysInitialize() override
Override sysInitialize.
#define ATH_MSG_WARNING(x)
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
SG::ThinningHandleKey< xAOD::TrackParticleContainer > m_MSTrkKey
Class describing a TrackParticle.
Property holding a SG store/key/clid/attr name from which a ReadDecorHandle is made.
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>
SG::ThinningHandleKey< xAOD::TrackParticleContainer > m_FwdIdTrkKey
bool empty() const noexcept
Returns true if the collection is empty.