![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
Implementation for the photon vertex selection tool.
More...
#include <PhotonVertexSelectionTool.h>
|
enum | FailType {
Unkown = -99,
NoFail = 0,
NoVxCont = 1,
NoEventInfo = 2,
FailPointing = 3,
FailEgamVect = 4,
NoGdCandidate = 5,
MatchedTrack = 6
} |
| Declare the interface that the class provides. More...
|
|
enum | yyVtxType { Unknown = -1,
ConvTrack = 0,
SiConvTrack = 1,
NoSiTracks = 2
} |
|
|
| PhotonVertexSelectionTool (const std::string &name) |
|
virtual | ~PhotonVertexSelectionTool () |
|
virtual void | print () const =0 |
| Print the state of the tool. More...
|
|
virtual void | print () const |
| Print the state of the tool. More...
|
|
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 | sysInitialize () override |
| Perform system initialization for an algorithm. 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 |
|
|
virtual StatusCode | initialize () |
| Function initialising the tool. More...
|
|
|
StatusCode | decorateInputs (const xAOD::EgammaContainer &egammas, FailType *failType=nullptr) const |
| Given a list of photons, decorate vertex container with MVA variables. More...
|
|
StatusCode | getVertex (const xAOD::EgammaContainer &egammas, const xAOD::Vertex *&vertex, bool ignoreConv=false) const |
| Given a list of photons, return the most likely vertex based on MVA likelihood. More...
|
|
std::vector< std::pair< const xAOD::Vertex *, float > > | getVertex (const xAOD::EgammaContainer &egammas, bool ignoreConv=false, bool noDecorate=false, yyVtxType *vtxCase=nullptr, FailType *failType=nullptr) const |
| Given a list of photons, return the MLPs of all vertices in the event. More...
|
|
int | getCase () const |
| Return the last case treated: More...
|
|
const xAOD::Vertex * | getPrimaryVertexFromConv (const xAOD::PhotonContainer *photons) const |
| Get possible vertex directly associated with photon conversions. More...
|
|
|
TLorentzVector | getEgammaVector (const xAOD::EgammaContainer *egammas, FailType &failType) const |
| Get combined 4-vector of photon container. More...
|
|
StatusCode | getVertexImp (const xAOD::EgammaContainer &egammas, const xAOD::Vertex *&vertex, bool ignoreConv, bool noDecorate, std::vector< std::pair< const xAOD::Vertex *, float > > &, yyVtxType &, FailType &) const |
| Given a list of photons, return the MLPs of all vertices in the event. More...
|
|
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...
|
|
Implementation for the photon vertex selection tool.
Takes a list of photons (for example, to two leading photons) and the most likely primary vertex, based on an MVA.
- Author
- Christopher Meyer chris.nosp@m..mey.nosp@m.er@ce.nosp@m.rn.c.nosp@m.h
-
Bruno Lenzi bruno.nosp@m..len.nosp@m.zi@ce.nosp@m.rn.c.nosp@m.h
Definition at line 34 of file PhotonVertexSelectionTool.h.
◆ StoreGateSvc_t
◆ FailType
Declare the interface that the class provides.
enum indicating where the tool has failed
Enumerator |
---|
Unkown | |
NoFail | |
NoVxCont | |
NoEventInfo | |
FailPointing | |
FailEgamVect | |
NoGdCandidate | |
MatchedTrack | |
Definition at line 33 of file IPhotonVertexSelectionTool.h.
◆ yyVtxType
◆ PhotonVertexSelectionTool()
CP::PhotonVertexSelectionTool::PhotonVertexSelectionTool |
( |
const std::string & |
name | ) |
|
◆ ~PhotonVertexSelectionTool()
CP::PhotonVertexSelectionTool::~PhotonVertexSelectionTool |
( |
| ) |
|
|
virtualdefault |
◆ 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]
◆ decorateInputs()
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ 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
◆ getCase()
int CP::PhotonVertexSelectionTool::getCase |
( |
| ) |
const |
|
inlinevirtual |
◆ getEgammaVector()
Get combined 4-vector of photon container.
Definition at line 360 of file PhotonVertexSelectionTool.cxx.
362 TLorentzVector
v, v1;
367 failType = FailType::FailEgamVect;
370 cluster =
egamma->caloCluster();
371 if (cluster ==
nullptr) {
372 ATH_MSG_WARNING(
"No cluster associated to egamma, not adding to 4-vector.");
◆ getKey()
Get the (hashed) key of an object that is in the event store.
This is a bit of a special one. StoreGateSvc
and xAOD::TEvent
both provide ways for getting the SG::sgkey_t
key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.
In order to allow tools to efficiently perform this operation, they can use this helper function.
- See also
- asg::AsgTool::getName
- Parameters
-
ptr | The bare pointer to the object that the event store should know about |
- Returns
- The hashed key of the object in the store. If not found, an invalid (zero) key.
Definition at line 119 of file AsgTool.cxx.
121 #ifdef XAOD_STANDALONE
124 return evtStore()->event()->getKey( ptr );
127 return (
proxy ==
nullptr ? 0 :
proxy->sgkey() );
128 #endif // XAOD_STANDALONE
◆ getName()
const std::string & asg::AsgTool::getName |
( |
const void * |
ptr | ) |
const |
|
inherited |
Get the name of an object that is / should be in the event store.
This is a bit of a special one. StoreGateSvc
and xAOD::TEvent
both provide ways for getting the std::string
name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.
In order to allow tools to efficiently perform this operation, they can use this helper function.
- See also
- asg::AsgTool::getKey
- Parameters
-
ptr | The bare pointer to the object that the event store should know about |
- Returns
- The string name of the object in the store. If not found, an empty string.
Definition at line 106 of file AsgTool.cxx.
108 #ifdef XAOD_STANDALONE
111 return evtStore()->event()->getName( ptr );
114 static const std::string
dummy =
"";
116 #endif // XAOD_STANDALONE
◆ getPrimaryVertexFromConv()
Get possible vertex directly associated with photon conversions.
Implements CP::IPhotonVertexSelectionTool.
Definition at line 306 of file PhotonVertexSelectionTool.cxx.
308 if (photons ==
nullptr) {
309 ATH_MSG_WARNING(
"Passed nullptr photon container, returning nullptr vertex from getPrimaryVertexFromConv");
313 std::vector<const xAOD::Vertex*> vertices;
316 size_t NumberOfTracks = 0;
322 for (
const auto *
photon: *photons) {
323 conversionVertex =
photon->vertex();
324 if (conversionVertex ==
nullptr)
continue;
327 for (
size_t i = 0;
i < NumberOfTracks; ++
i) {
330 if (gsfTp ==
nullptr)
continue;
335 if (
tp ==
nullptr)
continue;
338 if (
primary ==
nullptr)
continue;
342 if (
std::find(vertices.begin(), vertices.end(),
primary) == vertices.end()) {
350 if (!vertices.empty()) {
351 if (vertices.size() > 1)
352 ATH_MSG_WARNING(
"Photons associated to different vertices! Returning lead photon association.");
◆ getProperty()
template<class T >
const T* asg::AsgTool::getProperty |
( |
const std::string & |
name | ) |
const |
|
inherited |
Get one of the tool's properties.
◆ getVertex() [1/2]
Given a list of photons, return the MLPs of all vertices in the event.
Implements CP::IPhotonVertexSelectionTool.
Definition at line 150 of file PhotonVertexSelectionTool.cxx.
153 std::vector<std::pair<const xAOD::Vertex*, float> > vertexMLP;
155 FailType failType = FailType::NoFail;
156 if (
getVertexImp( egammas,
vertex, ignoreConv, noDecorate, vertexMLP, vtxCase, failType ).isSuccess()) {
159 if(vtxCasePtr!=
nullptr)
160 *vtxCasePtr = vtxCase;
161 if(failTypePtr!=
nullptr)
162 *failTypePtr = failType;
◆ getVertex() [2/2]
◆ getVertexImp()
Given a list of photons, return the MLPs of all vertices in the event.
Definition at line 178 of file PhotonVertexSelectionTool.cxx.
192 return StatusCode::FAILURE;
197 if (!ignoreConv && photons) {
199 if (prime_vertex !=
nullptr) {
200 vtxCase = yyVtxType::ConvTrack;
201 fail = FailType::MatchedTrack;
202 vertexMLP.emplace_back(prime_vertex, 0.);
203 return StatusCode::SUCCESS;
207 if (
fail != FailType::NoFail){
211 return StatusCode::SUCCESS;
220 vtxCase = yyVtxType::NoSiTracks;
221 if (!ignoreConv && photons) {
222 for (
const auto *
photon: *photons) {
226 return StatusCode::FAILURE;
232 vtxCase = yyVtxType::SiConvTrack;
247 float mlp = 0.0, mlp_max = -99999.0;
257 sumPt = log10((sumPtA)(*
vertex));
258 sumPt2 = log10((sumPt2A)(*
vertex));
262 " log(sumPt2): " << sumPt2 <<
273 TString mvaMethod(
"MLP method");
274 mlp =
reader->EvaluateMVA(mvaInput,mvaMethod);
278 vertexMLP.emplace_back(
vertex, mlp);
288 if (mlp_max <= -9999.0) {
289 ATH_MSG_DEBUG(
"No good vertex candidates from pointing, returning hardest vertex.");
291 fail = FailType::NoGdCandidate;
297 return StatusCode::SUCCESS;
◆ initialize()
StatusCode CP::PhotonVertexSelectionTool::initialize |
( |
| ) |
|
|
virtual |
Function initialising the tool.
Reimplemented from asg::AsgTool.
Definition at line 70 of file PhotonVertexSelectionTool.cxx.
72 ATH_MSG_INFO(
"Initializing PhotonVertexSelectionTool...");
79 std::vector<std::string> var_names = {
"deltaZ := TMath::Min(abs(PrimaryVerticesAuxDyn.z-zCommon)/zCommonError,20)",
80 "deltaPhi := abs(deltaPhi(PrimaryVerticesAuxDyn.phi,egamma_phi))" ,
81 "logSumpt := log10(PrimaryVerticesAuxDyn.sumPt)" ,
82 "logSumpt2 := log10(PrimaryVerticesAuxDyn.sumPt2)" };
84 auto *mva1 =
new TMVA::Reader(var_names,
"!Silent:Color");
86 m_mva1 = std::unique_ptr<TMVA::Reader>( mva1 );
88 auto mva2 = std::make_unique<TMVA::Reader>(var_names,
"!Silent:Color");
90 m_mva2 = std::unique_ptr<TMVA::Reader>( std::move(mva2) );
95 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.
◆ msg() [1/2]
◆ msg() [2/2]
◆ msg_level_name()
const std::string & asg::AsgTool::msg_level_name |
( |
| ) |
const |
|
inherited |
A deprecated function for getting the message level's name.
Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:
MSG::name( msg().level() )
This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.
- Returns
- The string name of the current minimum message level that's printed
Definition at line 101 of file AsgTool.cxx.
◆ 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.
◆ print() [1/2]
void asg::AsgTool::print |
( |
| ) |
const |
|
virtualinherited |
Print the state of the tool.
Implements asg::IAsgTool.
Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ print() [2/2]
virtual void asg::IAsgTool::print |
( |
| ) |
const |
|
pure virtualinherited |
Print the state of the tool.
Implemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, asg::AsgTool, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
◆ renounce()
◆ renounceArray()
◆ sortMLP()
◆ sysInitialize()
◆ 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_configFileCase1
std::string CP::PhotonVertexSelectionTool::m_configFileCase1 |
|
private |
◆ m_configFileCase2
std::string CP::PhotonVertexSelectionTool::m_configFileCase2 |
|
private |
◆ m_convPtCut
float CP::PhotonVertexSelectionTool::m_convPtCut |
|
private |
◆ m_derivationPrefix
std::string CP::PhotonVertexSelectionTool::m_derivationPrefix |
|
private |
◆ m_detStore
◆ m_eventInfo
◆ m_evtStore
◆ m_mva1
std::unique_ptr<TMVA::Reader> CP::PhotonVertexSelectionTool::m_mva1 |
|
private |
◆ m_mva2
std::unique_ptr<TMVA::Reader> CP::PhotonVertexSelectionTool::m_mva2 |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vertexContainer
◆ m_vertexContainerName
std::string CP::PhotonVertexSelectionTool::m_vertexContainerName |
|
private |
◆ m_vhka
The documentation for this class was generated from the following files:
float getVertexSumPt(const xAOD::Vertex *vertex, int power=1, bool useAux=true)
Loop over track particles associated with vertex and return scalar sum of pT^power in GeV (from auxda...
size_t nTrackParticles() const
Get the number of tracks associated with this vertex.
std::string find(const std::string &s)
return a remapped string
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
setSAddress setEtaMS setDirPhiMS setDirZMS setBarrelRadius setEndcapAlpha setEndcapRadius setInterceptInner setEtaMap setEtaBin setIsTgcFailure setDeltaPt deltaPhi
std::pair< float, float > getZCommonAndError(const xAOD::EventInfo *eventInfo, const xAOD::EgammaContainer *egammas, float convPtCut=2e3)
Return zCommon and zCommonError.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
float phiBE(const unsigned layer) const
Get the phi in one layer of the EM Calo.
Helper class to provide constant type-safe access to aux data.
#define ATH_MSG_VERBOSE(x)
virtual void setOwner(IDataHandleHolder *o)=0
virtual double e() const
energy
float etaBE(const unsigned layer) const
Get the eta in one layer of the EM Calo.
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Description of a calorimeter cluster.
const xAOD::Vertex * getHardestVertex(const xAOD::VertexContainer *vertices)
Return vertex with highest sum pT^2.
Helper class to provide type-safe access to aux data.
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
const TrackParticle * trackParticle(size_t i) const
Get the pointer to a given track that was used in vertex reco.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
TLorentzVector getVertexMomentum(const xAOD::Vertex *vertex, bool useAux=true, const std::string &derivationPrefix="")
Return vector sum of tracks associated with vertex (from auxdata if available and useAux = true)
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
const xAOD::TrackParticle * getOriginalTrackParticleFromGSF(const xAOD::TrackParticle *trkPar)
Helper function for getting the "Original" Track Particle (i.e before GSF) via the GSF Track Particle...
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Class describing a Vertex.
void sort(typename std::reverse_iterator< DataModel_detail::iterator< DVL > > beg, typename std::reverse_iterator< DataModel_detail::iterator< DVL > > end, const Compare &comp)
Specialization of sort for DataVector/List.
#define ATH_MSG_WARNING(x)
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
bool passConvSelection(const xAOD::Photon *photon, float convPtCut=2e3)
Check if photon is converted, and tracks have Si hits and pass selection.
Class describing a TrackParticle.
reader
read the goodrunslist xml file(s)
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
const xAOD::Vertex * getVertexFromTrack(const xAOD::TrackParticle *track, const xAOD::VertexContainer *vertices)