Loading [MathJax]/extensions/tex2jax.js
 |
ATLAS Offline Software
|
Thin taus below a minimum pt threshold.
More...
#include <TauThinningAlg.h>
|
virtual StatusCode | initialize () override |
| Gaudi initialize method. More...
|
|
virtual StatusCode | execute (const EventContext &ctx) const override |
| Execute the algorithm. 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 |
|
|
Gaudi::Property< bool > | m_doVertexCorrection { this, "VertexCorrection", true, "Tau vertex correction" } |
|
StringProperty | m_streamName { this, "StreamName", "StreamAOD", "Name of the stream being thinned" } |
|
SG::ThinningHandleKey< xAOD::TauJetContainer > | m_taus { this, "Taus", "TauJets", "Tau container to thin" } |
|
SG::ThinningHandleKey< xAOD::TauTrackContainer > | m_tauTracks { this, "TauTracks", "TauTracks", "Tau track container to thin" } |
|
SG::ThinningHandleKey< xAOD::PFOContainer > | m_neutralPFOs { this, "TauNeutralPFOs", "TauNeutralParticleFlowObjects", "Tau neutral PFO container to thin" } |
|
SG::ThinningHandleKey< xAOD::CaloClusterContainer > | m_pi0clusters { this, "TauPi0Clusters", "TauPi0Clusters", "Tau pi0 cluster container to thin" } |
|
SG::ThinningHandleKey< CaloClusterCellLinkContainer > | m_pi0CellLinks { this, "TauPi0CellLinks", "TauPi0Clusters_links", "Tau pi0 cluster cell link container to thin" } |
|
SG::ThinningHandleKey< xAOD::ParticleContainer > | m_finalPi0s { this, "TauFinalPi0s", "TauFinalPi0s", "Tau final pi0 container to thin" } |
|
SG::ThinningHandleKey< xAOD::PFOContainer > | m_shotPFOs { this, "TauShotPFOs", "TauShotParticleFlowObjects", "Tau shot PFO container to thin" } |
|
SG::ThinningHandleKey< xAOD::CaloClusterContainer > | m_shotclusters { this, "TauShotClusters", "TauShotClusters", "Tau shot cluster container to thin" } |
|
SG::ThinningHandleKey< CaloClusterCellLinkContainer > | m_shotCellLinks { this, "TauShotCellLinks", "TauShotClusters_links", "Tau shot cluster cell link container to thin" } |
|
SG::ThinningHandleKey< xAOD::PFOContainer > | m_hadronicPFOs { this, "TauHadronicPFOs", "TauHadronicParticleFlowObjects", "Tau hadronic PFO container to thin" } |
|
SG::ThinningHandleKey< xAOD::VertexContainer > | m_secondaryVertices { this, "TauSecondaryVertices", "TauSecondaryVertices", "Tau secondary vertex container to thin" } |
|
SG::ThinningHandleKey< CaloCellContainer > | m_cells { this, "Cells", "AllCalo", "Cell container to thin" } |
|
SG::ThinningHandleKey< CaloClusterCellLinkContainer > | m_tauCellLinks { this, "TauCellLinks", "CaloCalTopoClusters_links", "Tau cluster cell link container to thin" } |
|
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 |
|
Thin taus below a minimum pt threshold.
Thinning algorithm that removes entries from all tau-related containers for tau candidates below a minimum pt threshold
Definition at line 25 of file TauThinningAlg.h.
◆ StoreGateSvc_t
◆ 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.
◆ 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 TauThinningAlg::execute |
( |
const EventContext & |
ctx | ) |
const |
|
overridevirtual |
Execute the algorithm.
- Parameters
-
ctx | Current event context. |
Definition at line 36 of file TauThinningAlg.cxx.
45 neutralPFOs.thinAll();
48 pi0clusters.thinAll();
51 pi0CellLinks.thinAll();
57 hadronicPFOs.thinAll();
60 secondaryVertices.thinAll();
66 tauCellLinks.thinAll();
71 std::optional<SG::ThinningHandle<xAOD::ParticleContainer> > finalPi0sOpt;
74 finalPi0sOpt->thinAll();
77 std::optional<SG::ThinningHandle<xAOD::CaloClusterContainer> > shotclustersOpt;
80 shotclustersOpt->thinAll();
83 std::optional<SG::ThinningHandle<CaloClusterCellLinkContainer> > shotCellLinksOpt;
86 shotCellLinksOpt->thinAll();
93 if (!acc_passThinning(*tau))
continue;
96 taus.keep(tau->index());
100 tauTracks.keep(
track->index());
109 TLorentzVector clusterP4 =
cluster->p4();
114 clusterP4 = vertexedCluster.p4();
117 if (clusterP4.DeltaR(tauAxis) > 0.2)
continue;
127 if(cellLinks_it != tauCellLinks->end()) {
128 size_t link_index =
std::distance(tauCellLinks->begin(), cellLinks_it);
129 tauCellLinks.keep(link_index);
140 if (
it.index() >=
cells->size()) {
149 for(
size_t i=0;
i<tau->nNeutralPFOs();
i++) {
151 neutralPFOs.keep(tau->neutralPFO(
i)->index());
155 pi0clusters.keep(
cluster->index());
160 if(cellLinks_it != pi0CellLinks->end()) {
161 size_t link_index =
std::distance(pi0CellLinks->begin(), cellLinks_it);
162 pi0CellLinks.keep(link_index);
172 if (
it.index() >=
cells->size()) {
182 for(
size_t i=0;
i<tau->nPi0s();
i++) {
183 finalPi0sOpt->keep(tau->pi0(
i)->index());
188 for(
size_t i=0;
i<tau->nShotPFOs();
i++) {
190 shotPFOs.keep(tau->shotPFO(
i)->index());
194 if (!cluster)
continue;
195 if (shotclustersOpt) {
196 shotclustersOpt->keep(
cluster->index());
201 if (shotCellLinksOpt) {
203 if(cellLinks_it != (*shotCellLinksOpt)->end()) {
204 size_t link_index =
std::distance((*shotCellLinksOpt)->begin(), cellLinks_it);
205 shotCellLinksOpt->keep(link_index);
216 if (
it.index() >=
cells->size()) {
225 for(
size_t i=0;
i<tau->nHadronicPFOs();
i++) {
226 hadronicPFOs.keep(tau->hadronicPFO(
i)->index());
230 if (tau->secondaryVertex() !=
nullptr) {
231 secondaryVertices.keep(tau->secondaryVertex()->index());
235 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.
◆ filterPassed()
◆ initialize()
StatusCode TauThinningAlg::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()
◆ 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()
◆ 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.
◆ updateVHKA()
◆ m_cells
◆ m_detStore
◆ m_doVertexCorrection
Gaudi::Property<bool> TauThinningAlg::m_doVertexCorrection { this, "VertexCorrection", true, "Tau vertex correction" } |
|
private |
◆ 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_finalPi0s
◆ m_hadronicPFOs
◆ m_neutralPFOs
◆ m_pi0CellLinks
◆ m_pi0clusters
◆ m_secondaryVertices
◆ m_shotCellLinks
◆ m_shotclusters
◆ m_shotPFOs
◆ m_streamName
StringProperty TauThinningAlg::m_streamName { this, "StreamName", "StreamAOD", "Name of the stream being thinned" } |
|
private |
◆ m_tauCellLinks
◆ m_taus
◆ m_tauTracks
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
Const iterator class for DataVector/DataList.
constexpr ParticleHypothesis particle[PARTICLEHYPOTHESES]
the array of masses
std::string find(const std::string &s)
return a remapped string
SG::ThinningHandleKey< xAOD::ParticleContainer > m_finalPi0s
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
SG::ThinningHandleKey< xAOD::PFOContainer > m_neutralPFOs
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)
const_iterator begin() const
const begin method
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
SG::ThinningHandleKey< xAOD::PFOContainer > m_shotPFOs
const Amg::Vector3D & position() const
Returns the 3-pos.
const std::string & key() const
Return the StoreGate ID for the referenced object.
Class providing the definition of the 4-vector interface.
Handle for requesting thinning for a data object.
StatusCode initialize(const std::string &stream, const std::string &qualifier, bool used=true)
Should be called during the initialize phase.
virtual void setOwner(IDataHandleHolder *o)=0
StringProperty m_streamName
Bookkeeping of cells that make up a cluster Simplified replacement for CaloCellLink,...
SG::ThinningHandleKey< CaloCellContainer > m_cells
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
virtual const DataObjIDColl & extraOutputDeps() const override
Return the list of extra output dependencies.
Description of a calorimeter cluster.
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
Gaudi::Property< bool > m_doVertexCorrection
::StatusCode StatusCode
StatusCode definition for legacy code.
SG::ThinningHandleKey< CaloClusterCellLinkContainer > m_tauCellLinks
Class describing a tau jet.
virtual StatusCode sysExecute(const EventContext &ctx) override
Execute an algorithm.
SG::ThinningHandleKey< xAOD::CaloClusterContainer > m_shotclusters
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
SG::ThinningHandleKey< CaloClusterCellLinkContainer > m_pi0CellLinks
SG::ThinningHandleKey< CaloClusterCellLinkContainer > m_shotCellLinks
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
SG::ThinningHandleKey< xAOD::TauTrackContainer > m_tauTracks
const_iterator end() const
const end method
Class describing a Vertex.
#define ATH_MSG_WARNING(x)
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
virtual StatusCode sysInitialize() override
Override sysInitialize.
Evaluate cluster kinematics with a different vertex / signal state.
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
SG::ThinningHandleKey< xAOD::CaloClusterContainer > m_pi0clusters
const_iterator to loop over cells belonging to a cluster
float distance(const Amg::Vector3D &p1, const Amg::Vector3D &p2)
calculates the distance between two point in 3D space
SG::ThinningHandleKey< xAOD::VertexContainer > m_secondaryVertices
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
SG::ThinningHandleKey< xAOD::PFOContainer > m_hadronicPFOs
SG::ThinningHandleKey< xAOD::TauJetContainer > m_taus