File indexing completed on 2025-07-01 08:35:22
0001
0002
0003 #ifndef EDM4HEP_ClusterCollection_H
0004 #define EDM4HEP_ClusterCollection_H
0005
0006
0007 #include "edm4hep/Cluster.h"
0008 #include "edm4hep/ClusterCollectionData.h"
0009 #include "edm4hep/ClusterObj.h"
0010 #include "edm4hep/MutableCluster.h"
0011
0012
0013 #include "podio/CollectionBase.h"
0014 #include "podio/ICollectionProvider.h"
0015
0016 #if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
0017 #include "nlohmann/json_fwd.hpp"
0018 #endif
0019
0020 #include <algorithm>
0021 #include <cstddef>
0022 #include <memory>
0023 #include <mutex>
0024 #include <ostream>
0025 #include <string_view>
0026 #include <vector>
0027
0028 namespace podio {
0029 struct RelationNames;
0030 }
0031
0032 namespace edm4hep {
0033
0034 class ClusterCollectionIterator {
0035 public:
0036 ClusterCollectionIterator(size_t index, const ClusterObjPointerContainer* collection) :
0037 m_index(index), m_object(podio::utils::MaybeSharedPtr<ClusterObj>{nullptr}), m_collection(collection) {
0038 }
0039
0040 ClusterCollectionIterator(const ClusterCollectionIterator&) = delete;
0041 ClusterCollectionIterator& operator=(const ClusterCollectionIterator&) = delete;
0042
0043 bool operator!=(const ClusterCollectionIterator& x) const {
0044 return m_index != x.m_index;
0045 }
0046
0047 bool operator==(const ClusterCollectionIterator& x) const {
0048 return m_index == x.m_index;
0049 }
0050
0051 Cluster operator*();
0052 Cluster* operator->();
0053 ClusterCollectionIterator& operator++();
0054
0055 private:
0056 size_t m_index;
0057 Cluster m_object;
0058 const ClusterObjPointerContainer* m_collection;
0059 };
0060
0061 class ClusterMutableCollectionIterator {
0062 public:
0063 ClusterMutableCollectionIterator(size_t index, const ClusterObjPointerContainer* collection) :
0064 m_index(index), m_object(podio::utils::MaybeSharedPtr<ClusterObj>{nullptr}), m_collection(collection) {
0065 }
0066
0067 ClusterMutableCollectionIterator(const ClusterMutableCollectionIterator&) = delete;
0068 ClusterMutableCollectionIterator& operator=(const ClusterMutableCollectionIterator&) = delete;
0069
0070 bool operator!=(const ClusterMutableCollectionIterator& x) const {
0071 return m_index != x.m_index;
0072 }
0073
0074 bool operator==(const ClusterMutableCollectionIterator& x) const {
0075 return m_index == x.m_index;
0076 }
0077
0078 MutableCluster operator*();
0079 MutableCluster* operator->();
0080 ClusterMutableCollectionIterator& operator++();
0081
0082 private:
0083 size_t m_index;
0084 MutableCluster m_object;
0085 const ClusterObjPointerContainer* m_collection;
0086 };
0087
0088
0089
0090
0091 class ClusterCollection : public podio::CollectionBase {
0092 public:
0093 using value_type = Cluster;
0094 using const_iterator = ClusterCollectionIterator;
0095 using iterator = ClusterMutableCollectionIterator;
0096 using difference_type = ptrdiff_t;
0097 using size_type = size_t;
0098
0099 ClusterCollection();
0100 ClusterCollection(ClusterCollectionData&& data, bool isSubsetColl);
0101
0102 ClusterCollection(const ClusterCollection&) = delete;
0103 ClusterCollection& operator=(const ClusterCollection&) = delete;
0104 ClusterCollection(ClusterCollection&&) = default;
0105 ClusterCollection& operator=(ClusterCollection&&) = default;
0106
0107
0108 ~ClusterCollection() override;
0109
0110 constexpr static auto typeName = "edm4hep::ClusterCollection";
0111 constexpr static auto valueTypeName = "edm4hep::Cluster";
0112 constexpr static auto dataTypeName = "edm4hep::ClusterData";
0113
0114 void clear() final;
0115
0116
0117 void print(std::ostream& os = std::cout, bool flush = true) const final;
0118
0119
0120 ClusterCollection* operator->() {
0121 return static_cast<ClusterCollection*>(this);
0122 }
0123
0124
0125 MutableCluster create();
0126
0127
0128
0129 template <typename... Args>
0130 MutableCluster create(Args&&... args);
0131
0132
0133 std::size_t size() const final;
0134
0135
0136 std::size_t max_size() const final;
0137
0138
0139 bool empty() const final;
0140
0141
0142 const std::string_view getTypeName() const final {
0143 return typeName;
0144 }
0145
0146 const std::string_view getValueTypeName() const final {
0147 return valueTypeName;
0148 }
0149
0150 const std::string_view getDataTypeName() const final {
0151 return dataTypeName;
0152 }
0153
0154 podio::SchemaVersionT getSchemaVersion() const final;
0155
0156 bool isSubsetCollection() const final {
0157 return m_isSubsetColl;
0158 }
0159
0160 void setSubsetCollection(bool setSubset = true) final;
0161
0162
0163 Cluster operator[](std::size_t index) const;
0164
0165 MutableCluster operator[](std::size_t index);
0166
0167 Cluster at(std::size_t index) const;
0168
0169 MutableCluster at(std::size_t index);
0170
0171
0172 void push_back(const MutableCluster& object);
0173
0174 void push_back(const Cluster& object);
0175
0176 void prepareForWrite() const final;
0177 void prepareAfterRead() final;
0178 bool setReferences(const podio::ICollectionProvider* collectionProvider) final;
0179
0180
0181 podio::CollectionWriteBuffers getBuffers() final;
0182
0183 void setID(uint32_t ID) final {
0184 m_collectionID = ID;
0185 if (!m_isSubsetColl) {
0186 std::for_each(m_storage.entries.begin(), m_storage.entries.end(), [ID](ClusterObj* obj) {
0187 obj->id = {obj->id.index, static_cast<uint32_t>(ID)};
0188 });
0189 }
0190 m_isValid = true;
0191 }
0192
0193 uint32_t getID() const final {
0194 return m_collectionID;
0195 }
0196
0197 bool isValid() const final {
0198 return m_isValid;
0199 }
0200
0201 size_t getDatamodelRegistryIndex() const final;
0202
0203
0204 iterator begin() {
0205 return iterator(0, &m_storage.entries);
0206 }
0207 const_iterator begin() const {
0208 return const_iterator(0, &m_storage.entries);
0209 }
0210 const_iterator cbegin() const {
0211 return begin();
0212 }
0213 iterator end() {
0214 return iterator(m_storage.entries.size(), &m_storage.entries);
0215 }
0216 const_iterator end() const {
0217 return const_iterator(m_storage.entries.size(), &m_storage.entries);
0218 }
0219 const_iterator cend() const {
0220 return end();
0221 }
0222
0223 std::vector<std::int32_t> type(const size_t nElem = 0) const;
0224 std::vector<float> energy(const size_t nElem = 0) const;
0225 std::vector<float> energyError(const size_t nElem = 0) const;
0226 std::vector<edm4hep::Vector3f> position(const size_t nElem = 0) const;
0227 std::vector<edm4hep::CovMatrix3f> positionError(const size_t nElem = 0) const;
0228 std::vector<float> iTheta(const size_t nElem = 0) const;
0229 std::vector<float> phi(const size_t nElem = 0) const;
0230 std::vector<edm4hep::Vector3f> directionError(const size_t nElem = 0) const;
0231
0232 private:
0233
0234
0235
0236 friend class ClusterCollectionData;
0237
0238 bool m_isValid{false};
0239 mutable bool m_isPrepared{false};
0240 bool m_isSubsetColl{false};
0241 uint32_t m_collectionID{0};
0242 mutable std::unique_ptr<std::mutex> m_storageMtx{nullptr};
0243 mutable ClusterCollectionData m_storage{};
0244 };
0245
0246 std::ostream& operator<<(std::ostream& o, const ClusterCollection& v);
0247
0248 template <typename... Args>
0249 MutableCluster ClusterCollection::create(Args&&... args) {
0250 if (m_isSubsetColl) {
0251 throw std::logic_error("Cannot create new elements on a subset collection");
0252 }
0253 const int size = m_storage.entries.size();
0254 auto obj = new ClusterObj({size, m_collectionID}, {std::forward<Args>(args)...});
0255 m_storage.entries.push_back(obj);
0256
0257
0258 obj->m_clusters = new std::vector<edm4hep::Cluster>();
0259 obj->m_hits = new std::vector<edm4hep::CalorimeterHit>();
0260 obj->m_shapeParameters = new std::vector<float>();
0261 obj->m_subdetectorEnergies = new std::vector<float>();
0262 m_storage.createRelations(obj);
0263 return MutableCluster(podio::utils::MaybeSharedPtr(obj));
0264 }
0265
0266 #if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
0267 void to_json(nlohmann::json& j, const ClusterCollection& collection);
0268 #endif
0269
0270 }
0271
0272 #endif