File indexing completed on 2025-11-04 10:09:36
0001 
0002 
0003 #ifndef EDM4HEP_MutableReconstructedParticle_H
0004 #define EDM4HEP_MutableReconstructedParticle_H
0005 
0006 #include "edm4hep/ReconstructedParticleObj.h"
0007 
0008 #include "edm4hep/ReconstructedParticle.h"
0009 
0010 #include "edm4hep/Cluster.h"
0011 #include "edm4hep/CovMatrix4f.h"
0012 #include "edm4hep/Track.h"
0013 #include "edm4hep/Vector3f.h"
0014 #include <edm4hep/Constants.h>
0015 #include "podio/RelationRange.h"
0016 #include <cstdint>
0017 #include <vector>
0018 
0019 #include "podio/utilities/MaybeSharedPtr.h"
0020 
0021 #include <cstdint>
0022 
0023 #if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
0024 #include "nlohmann/json_fwd.hpp"
0025 #endif
0026 
0027 
0028 namespace edm4hep {
0029 class ReconstructedParticleCollection;
0030 class Vertex;
0031 class MutableVertex;
0032 }
0033 
0034 
0035 namespace edm4hep {
0036 
0037 
0038 
0039 
0040 
0041 
0042 class MutableReconstructedParticle {
0043 
0044   friend class ReconstructedParticleCollection;
0045   friend class ReconstructedParticleMutableCollectionIterator;
0046   friend class ReconstructedParticle;
0047 
0048 public:
0049   using object_type = ReconstructedParticle;
0050   using collection_type = ReconstructedParticleCollection;
0051 
0052   
0053   MutableReconstructedParticle();
0054 
0055   
0056   MutableReconstructedParticle(const std::int32_t PDG, const float energy, const edm4hep::Vector3f& momentum, const edm4hep::Vector3f& referencePoint, const float charge, const float mass, const float goodnessOfPID, const edm4hep::CovMatrix4f& covMatrix);
0057 
0058   
0059   MutableReconstructedParticle(const MutableReconstructedParticle& other) = default;
0060 
0061   
0062   MutableReconstructedParticle& operator=(MutableReconstructedParticle other) &; 
0063   MutableReconstructedParticle& operator=(MutableReconstructedParticle other) && = delete; 
0064 
0065   
0066   
0067   MutableReconstructedParticle clone(bool cloneRelations=true) const;
0068 
0069   
0070   ~MutableReconstructedParticle() = default;
0071 
0072 
0073 public:
0074 
0075   
0076   std::int32_t getPDG() const;
0077 
0078   
0079   float getEnergy() const;
0080 
0081   
0082   const edm4hep::Vector3f& getMomentum() const;
0083 
0084   
0085   const edm4hep::Vector3f& getReferencePoint() const;
0086 
0087   
0088   float getCharge() const;
0089 
0090   
0091   float getMass() const;
0092 
0093   
0094   float getGoodnessOfPID() const;
0095 
0096   
0097   const edm4hep::CovMatrix4f& getCovMatrix() const;
0098 
0099 
0100   
0101   const edm4hep::Vertex getDecayVertex() const;
0102 
0103   
0104   void setPDG(const std::int32_t PDG);
0105   
0106   std::int32_t& getPDG();
0107   
0108   [[deprecated("use getPDG instead")]]
0109   std::int32_t& PDG();
0110 
0111   
0112   void setEnergy(const float energy);
0113   
0114   float& getEnergy();
0115   
0116   [[deprecated("use getEnergy instead")]]
0117   float& energy();
0118 
0119   
0120   void setMomentum(const edm4hep::Vector3f& momentum);
0121   
0122   edm4hep::Vector3f& getMomentum();
0123   
0124   [[deprecated("use getMomentum instead")]]
0125   edm4hep::Vector3f& momentum();
0126 
0127   
0128   void setReferencePoint(const edm4hep::Vector3f& referencePoint);
0129   
0130   edm4hep::Vector3f& getReferencePoint();
0131   
0132   [[deprecated("use getReferencePoint instead")]]
0133   edm4hep::Vector3f& referencePoint();
0134 
0135   
0136   void setCharge(const float charge);
0137   
0138   float& getCharge();
0139   
0140   [[deprecated("use getCharge instead")]]
0141   float& charge();
0142 
0143   
0144   void setMass(const float mass);
0145   
0146   float& getMass();
0147   
0148   [[deprecated("use getMass instead")]]
0149   float& mass();
0150 
0151   
0152   void setGoodnessOfPID(const float goodnessOfPID);
0153   
0154   float& getGoodnessOfPID();
0155   
0156   [[deprecated("use getGoodnessOfPID instead")]]
0157   float& goodnessOfPID();
0158 
0159   
0160   void setCovMatrix(const edm4hep::CovMatrix4f& covMatrix);
0161   
0162   edm4hep::CovMatrix4f& getCovMatrix();
0163   
0164   [[deprecated("use getCovMatrix instead")]]
0165   edm4hep::CovMatrix4f& covMatrix();
0166 
0167 
0168   
0169   void setDecayVertex(const edm4hep::Vertex& value);
0170 
0171   void addToClusters(const edm4hep::Cluster&);
0172   std::size_t clusters_size() const;
0173   edm4hep::Cluster getClusters(std::size_t) const;
0174   std::vector<edm4hep::Cluster>::const_iterator clusters_begin() const;
0175   std::vector<edm4hep::Cluster>::const_iterator clusters_end() const;
0176   podio::RelationRange<edm4hep::Cluster> getClusters() const;
0177   void addToTracks(const edm4hep::Track&);
0178   std::size_t tracks_size() const;
0179   edm4hep::Track getTracks(std::size_t) const;
0180   std::vector<edm4hep::Track>::const_iterator tracks_begin() const;
0181   std::vector<edm4hep::Track>::const_iterator tracks_end() const;
0182   podio::RelationRange<edm4hep::Track> getTracks() const;
0183   void addToParticles(const edm4hep::ReconstructedParticle&);
0184   std::size_t particles_size() const;
0185   edm4hep::ReconstructedParticle getParticles(std::size_t) const;
0186   std::vector<edm4hep::ReconstructedParticle>::const_iterator particles_begin() const;
0187   std::vector<edm4hep::ReconstructedParticle>::const_iterator particles_end() const;
0188   podio::RelationRange<edm4hep::ReconstructedParticle> getParticles() const;
0189 
0190  bool isCompound() const { return particles_size() > 0 ;}
0191  [[deprecated("use setPDG instead")]]
0192  int32_t getType() const { return getPDG(); }
0193  
0194  float getCovMatrix(edm4hep::FourMomCoords dimI, edm4hep::FourMomCoords dimJ) const { return getCovMatrix().getValue(dimI, dimJ); }
0195  
0196 
0197  
0198  
0199  [[deprecated("use setPDG instead")]]
0200  void setType(int32_t pdg) { setPDG(pdg); }
0201  
0202  void setCovMatrix(float value, edm4hep::FourMomCoords dimI, edm4hep::FourMomCoords dimJ) { getCovMatrix().setValue(value, dimI, dimJ); }
0203  
0204 
0205   
0206   bool isAvailable() const;
0207   
0208   void unlink() { m_obj = podio::utils::MaybeSharedPtr<ReconstructedParticleObj>{nullptr}; }
0209 
0210   bool operator==(const MutableReconstructedParticle& other) const { return m_obj == other.m_obj; }
0211   bool operator==(const ReconstructedParticle& other) const;
0212 
0213   bool operator!=(const MutableReconstructedParticle& other) const { return !(*this == other); }
0214   bool operator!=(const ReconstructedParticle& other) const { return !(*this == other); }
0215 
0216   
0217   bool operator<(const MutableReconstructedParticle& other) const { return podio::detail::getOrderKey(*this) < podio::detail::getOrderKey(other); }
0218 
0219   podio::ObjectID id() const { return getObjectID(); }
0220 
0221   const podio::ObjectID getObjectID() const;
0222 
0223   friend std::hash<MutableReconstructedParticle>;
0224 
0225   friend void swap(MutableReconstructedParticle& a, MutableReconstructedParticle& b) {
0226     using std::swap;
0227     swap(a.m_obj, b.m_obj); 
0228   }
0229 
0230 private:
0231   
0232   explicit MutableReconstructedParticle(podio::utils::MaybeSharedPtr<ReconstructedParticleObj> obj);
0233 
0234   podio::utils::MaybeSharedPtr<ReconstructedParticleObj> m_obj{nullptr};
0235 };
0236 
0237 #if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
0238 void to_json(nlohmann::json& j, const MutableReconstructedParticle& value);
0239 #endif
0240 
0241 
0242 } 
0243 
0244 
0245 
0246 template<>
0247 struct std::hash<edm4hep::MutableReconstructedParticle> {
0248   std::size_t operator()(const edm4hep::MutableReconstructedParticle& obj) const {
0249     return std::hash<edm4hep::ReconstructedParticleObj*>{}(obj.m_obj.get());
0250   }
0251 };
0252 
0253 
0254 #endif