![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
#include <SecVertexMergingTool.h>
|
virtual StatusCode | initialize () override |
|
virtual StatusCode | finalize () override |
| EndOfInitialize. More...
|
|
| SecVertexMergingTool (const std::string &t, const std::string &n, const IInterface *p) |
| constructor More...
|
|
virtual | ~SecVertexMergingTool () |
| destructor More...
|
|
virtual std::pair< xAOD::VertexContainer *, xAOD::VertexAuxContainer * > | mergeVertexContainer (const xAOD::VertexContainer &MyVxCont) const override |
| Merging
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 |
|
- Author
- Maria Ilaria Besana, June 2014
Changes:
David Shope david.nosp@m..ric.nosp@m.hard..nosp@m.shop.nosp@m.e@cer.nosp@m.n.ch (2016-04-26)
EDM Migration to xAOD - remove method using VxCandidate
Definition at line 41 of file SecVertexMergingTool.h.
◆ StoreGateSvc_t
◆ SecVertexMergingTool()
Trk::SecVertexMergingTool::SecVertexMergingTool |
( |
const std::string & |
t, |
|
|
const std::string & |
n, |
|
|
const IInterface * |
p |
|
) |
| |
◆ ~SecVertexMergingTool()
Trk::SecVertexMergingTool::~SecVertexMergingTool |
( |
| ) |
|
|
virtualdefault |
◆ checkCompatibility()
Definition at line 242 of file SecVertexMergingTool.cxx.
249 AmgSymMatrix(3) vErrs = v1->covariancePosition() +
v2->covariancePosition() ;
250 vErrs = vErrs.inverse().eval();
254 sigma = sqrt( vdif.dot( vErrs * vdif ) ) ;
257 sigma = vdif(0)*vdif(0)*vErrs(0,0) + vdif(1)*vdif(1)*vErrs(1,1) + 2*vdif(0)*vdif(1)*vErrs(0,1) ;
262 sigma = vdif(2)*sqrt( vErrs(2,2) );
◆ 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]
◆ 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
◆ finalize()
StatusCode Trk::SecVertexMergingTool::finalize |
( |
| ) |
|
|
overridevirtual |
◆ initialize()
StatusCode Trk::SecVertexMergingTool::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.
◆ interfaceID()
static const InterfaceID& Trk::IVertexMergingTool::interfaceID |
( |
| ) |
|
|
inlinestaticinherited |
◆ mergeVertexContainer()
Merging
SG::Decorator<std::vector<float> > mDecor_trkWght(
"trkWeight" ) ; SG::Decorator<float> mDecor_trkDOE( "trkDistOverError" ) ; SG::Decorator<float> mDecor_direction( "MomentaDirection" ); SG::Decorator< float > mDecor_HitsFilter( "radiiPattern" );
Implements Trk::IVertexMergingTool.
Definition at line 50 of file SecVertexMergingTool.cxx.
60 NewContainer->setStore(auxNewContainer);
78 bool moreDeco = mAcc_momdir.isAvailable(*MyVxCont.
front());
92 std::vector<bool> remerged(MyVxCont.
size(),
false);
104 std::vector<const xAOD::TrackParticle*> combinedTracks;
105 std::vector<ElementLink<xAOD::TrackParticleContainer>> tpLinks1 =
107 if (!tpLinks1.empty()) {
108 for (
const auto& tp_EL : tpLinks1) {
110 combinedTracks.push_back(trk);
113 unsigned int Nj = Ni + 1;
114 bool newmerge =
false;
122 combinedTracks.clear();
124 if (tpLinks1.empty())
126 for (
const auto& tp_EL : tpLinks1) {
128 combinedTracks.push_back(trk);
141 const std::vector<ElementLink<xAOD::TrackParticleContainer>>
143 if (tpLinks2.empty())
146 for (
const auto& tp_EL : tpLinks2) {
148 combinedTracks.push_back(trk);
169 float pt1 = sqrt(mAcc_sumPt2(*vx));
170 float pt2 = sqrt(mAcc_sumPt2(*mergeCand));
174 0.6 * pt1 / (pt1 + pt2) + 0.4 * ntrk1 / (ntrk1 + ntrk2);
176 0.6 * pt2 / (pt1 + pt2) + 0.4 * ntrk2 / (ntrk1 + ntrk2);
180 float mas1 = mAcc_mass(*vx);
181 float mas2 = mAcc_mass(*mergeCand);
182 float e1 = mAcc_energy(*vx);
183 float e2 = mAcc_energy(*mergeCand);
184 int inNtrk1 = mAcc_intrk(*vx);
185 int inNtrk2 = mAcc_intrk(*mergeCand);
188 float md1 = 0., md2 = 0., hf1 = 0., hf2 = 0.;
189 std::vector<float> trkW1, trkW2, doe1, doe2;
191 doe1 = mAcc_trkdoe(*vx);
192 doe2 = mAcc_trkdoe(*mergeCand);
193 doe2.insert(doe2.end(), doe1.begin(), doe1.end());
194 md1 = mAcc_momdir(*vx);
195 md2 = mAcc_momdir(*mergeCand);
196 hf1 = mAcc_radpat(*vx);
197 hf2 = mAcc_radpat(*mergeCand);
198 trkW1 = mAcc_trkwt(*vx);
199 trkW2 = mAcc_trkwt(*mergeCand);
200 trkW2.insert(trkW2.end(), trkW1.begin(), trkW1.end());
201 ntrks = mAcc_numtav(*vx) + mAcc_numtav(*mergeCand);
214 mAcc_trkdoe(*vx) = doe2;
215 mAcc_momdir(*vx) = wght1 * md1 + wght2 * md2;
216 mAcc_radpat(*vx) = wght1 * hf1 + wght2 * hf2;
217 mAcc_trkwt(*vx) = trkW2;
218 mAcc_numtav(*vx) = ntrks;
221 mDecor_sumPt2(*vx) = pt1 * pt1 + pt2 * pt2;
222 mDecor_mass(*vx) = wght1 * mas1 + wght2 * mas2;
223 mDecor_energy(*vx) = wght1 *
e1 + wght2 *
e2;
224 mDecor_nrobbed(*vx) = 0;
225 mDecor_intrk(*vx) = (
int)(wght1 * inNtrk1 + wght1 * inNtrk2);
237 return std::make_pair(NewContainer, auxNewContainer);
◆ 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()
◆ 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_Compatidime
int Trk::SecVertexMergingTool::m_Compatidime |
|
private |
◆ m_detStore
◆ m_evtStore
◆ m_iVertexFitter
◆ m_minDist
float Trk::SecVertexMergingTool::m_minDist |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
Const iterator class for DataVector/DataList.
size_t nTrackParticles() const
Get the number of tracks associated with this vertex.
Temporary container used until we have I/O for AuxStoreInternal.
Vertex_v1 Vertex
Define the latest version of the vertex class.
double e1(const xAOD::CaloCluster &cluster)
return the uncorrected cluster energy in 1st sampling
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
const TrackParticleLinks_t & trackParticleLinks() const
Get all the particles associated with the vertex.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
const Amg::Vector3D & position() const
Returns the 3-pos.
VertexContainer_v1 VertexContainer
Definition of the current "Vertex container version".
VxType::VertexType vertexType() const
The type of the vertex.
virtual void setOwner(IDataHandleHolder *o)=0
AmgSymMatrix(5) &GXFTrackState
void setVertexType(VxType::VertexType vType)
Set the type of the vertex.
VertexAuxContainer_v1 VertexAuxContainer
Definition of the current jet auxiliary container.
const T * front() const
Access the first element in the collection as an rvalue.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
value_type push_back(value_type pElem)
Add an element to the end of the collection.
Eigen::Matrix< double, 3, 1 > Vector3D
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
Class describing a Vertex.
double e2(const xAOD::CaloCluster &cluster)
return the uncorrected cluster energy in 2nd sampling
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Class describing a TrackParticle.
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_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.