File indexing completed on 2026-05-04 08:36:57
0001
0002
0003 #ifndef EDM4EIC_TrackCollection_H
0004 #define EDM4EIC_TrackCollection_H
0005
0006
0007 #include "edm4eic/Track.h"
0008 #include "edm4eic/MutableTrack.h"
0009 #include "edm4eic/TrackObj.h"
0010 #include "edm4eic/TrackCollectionData.h"
0011
0012
0013 #include "podio/ICollectionProvider.h"
0014 #include "podio/CollectionBase.h"
0015 #include "podio/detail/Pythonizations.h"
0016
0017 #if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
0018 #include "nlohmann/json_fwd.hpp"
0019 #endif
0020
0021 #include <string_view>
0022 #include <vector>
0023 #include <algorithm>
0024 #include <ostream>
0025 #include <mutex>
0026 #include <memory>
0027 #include <cstddef>
0028
0029 namespace podio {
0030 struct RelationNames;
0031 }
0032
0033 namespace edm4eic {
0034
0035
0036
0037 class TrackCollectionIterator {
0038 public:
0039 using value_type = Track;
0040 using difference_type = ptrdiff_t;
0041 using reference = Track;
0042 using pointer = Track*;
0043 using iterator_category = std::input_iterator_tag;
0044
0045
0046 using iterator_concept = std::random_access_iterator_tag;
0047
0048 TrackCollectionIterator(size_t index, const TrackObjPointerContainer* collection) : m_index(index), m_object(podio::utils::MaybeSharedPtr<TrackObj>{nullptr}), m_collection(collection) {}
0049 TrackCollectionIterator() = default;
0050
0051 TrackCollectionIterator(const TrackCollectionIterator&) = default;
0052 TrackCollectionIterator(TrackCollectionIterator&&) = default;
0053 TrackCollectionIterator& operator=(const TrackCollectionIterator&) = default;
0054 TrackCollectionIterator& operator=(TrackCollectionIterator&&) = default;
0055 ~TrackCollectionIterator() = default;
0056
0057 auto operator<=>(const TrackCollectionIterator& other) const {
0058 return m_index <=> other.m_index;
0059 }
0060
0061 bool operator==(const TrackCollectionIterator& x) const {
0062 return m_index == x.m_index;
0063 }
0064
0065 reference operator*() const;
0066 pointer operator->();
0067 TrackCollectionIterator& operator++();
0068 TrackCollectionIterator operator++(int);
0069 TrackCollectionIterator& operator--();
0070 TrackCollectionIterator operator--(int);
0071 TrackCollectionIterator& operator+=(difference_type n);
0072 TrackCollectionIterator operator+(difference_type n) const;
0073 friend TrackCollectionIterator operator+(difference_type n, const TrackCollectionIterator& it);
0074 TrackCollectionIterator& operator-=(difference_type n);
0075 TrackCollectionIterator operator-(difference_type n) const;
0076 reference operator[](difference_type n) const;
0077 difference_type operator-(const TrackCollectionIterator& other) const;
0078
0079 private:
0080 size_t m_index{0};
0081 Track m_object { podio::utils::MaybeSharedPtr<TrackObj>{nullptr} };
0082 const TrackObjPointerContainer* m_collection{nullptr};
0083 };
0084
0085
0086 class TrackMutableCollectionIterator {
0087 public:
0088 using value_type = Track;
0089 using difference_type = ptrdiff_t;
0090 using reference = MutableTrack;
0091 using pointer = MutableTrack*;
0092 using iterator_category = std::input_iterator_tag;
0093
0094
0095 using iterator_concept = std::random_access_iterator_tag;
0096
0097 TrackMutableCollectionIterator(size_t index, const TrackObjPointerContainer* collection) : m_index(index), m_object(podio::utils::MaybeSharedPtr<TrackObj>{nullptr}), m_collection(collection) {}
0098 TrackMutableCollectionIterator() = default;
0099
0100 TrackMutableCollectionIterator(const TrackMutableCollectionIterator&) = default;
0101 TrackMutableCollectionIterator(TrackMutableCollectionIterator&&) = default;
0102 TrackMutableCollectionIterator& operator=(const TrackMutableCollectionIterator&) = default;
0103 TrackMutableCollectionIterator& operator=(TrackMutableCollectionIterator&&) = default;
0104 ~TrackMutableCollectionIterator() = default;
0105
0106 auto operator<=>(const TrackMutableCollectionIterator& other) const {
0107 return m_index <=> other.m_index;
0108 }
0109
0110 bool operator==(const TrackMutableCollectionIterator& x) const {
0111 return m_index == x.m_index;
0112 }
0113
0114 reference operator*() const;
0115 pointer operator->();
0116 TrackMutableCollectionIterator& operator++();
0117 TrackMutableCollectionIterator operator++(int);
0118 TrackMutableCollectionIterator& operator--();
0119 TrackMutableCollectionIterator operator--(int);
0120 TrackMutableCollectionIterator& operator+=(difference_type n);
0121 TrackMutableCollectionIterator operator+(difference_type n) const;
0122 friend TrackMutableCollectionIterator operator+(difference_type n, const TrackMutableCollectionIterator& it);
0123 TrackMutableCollectionIterator& operator-=(difference_type n);
0124 TrackMutableCollectionIterator operator-(difference_type n) const;
0125 reference operator[](difference_type n) const;
0126 difference_type operator-(const TrackMutableCollectionIterator& other) const;
0127
0128 private:
0129 size_t m_index{0};
0130 MutableTrack m_object { podio::utils::MaybeSharedPtr<TrackObj>{nullptr} };
0131 const TrackObjPointerContainer* m_collection{nullptr};
0132 };
0133
0134
0135
0136
0137
0138 class TrackCollection : public podio::CollectionBase {
0139 public:
0140 using value_type = Track;
0141 using mutable_type = MutableTrack;
0142 using const_iterator = TrackCollectionIterator;
0143 using iterator = TrackMutableCollectionIterator;
0144 using difference_type = ptrdiff_t;
0145 using size_type = size_t;
0146 using const_reverse_iterator = std::reverse_iterator<const_iterator>;
0147 using reverse_iterator = std::reverse_iterator<iterator>;
0148
0149 TrackCollection() = default;
0150 TrackCollection(TrackCollectionData&& data, bool isSubsetColl);
0151
0152 TrackCollection(const TrackCollection& ) = delete;
0153 TrackCollection& operator=(const TrackCollection& ) = delete;
0154 TrackCollection(TrackCollection&&) = default;
0155 TrackCollection& operator=(TrackCollection&&) = default;
0156
0157
0158 ~TrackCollection() override;
0159
0160 constexpr static std::string_view typeName = "edm4eic::TrackCollection";
0161 constexpr static std::string_view valueTypeName = "edm4eic::Track";
0162 constexpr static std::string_view dataTypeName = "edm4eic::TrackData";
0163
0164 void clear() final;
0165
0166
0167 static void __cppyy_pythonize__(PyObject* klass, const std::string& name){
0168 podio::detail::pythonizations::pythonize_subscript(klass, name);
0169 }
0170
0171
0172 void print(std::ostream& os=std::cout, bool flush=true) const final;
0173
0174
0175 MutableTrack create();
0176
0177
0178
0179 template<typename... Args>
0180 MutableTrack create(Args&&... args);
0181
0182
0183 std::size_t size() const final;
0184
0185
0186 std::size_t max_size() const final;
0187
0188
0189 bool empty() const final;
0190
0191
0192 const std::string_view getTypeName() const final { return typeName; }
0193
0194 const std::string_view getValueTypeName() const final { return valueTypeName; }
0195
0196 const std::string_view getDataTypeName() const final { return dataTypeName; }
0197
0198 podio::SchemaVersionT getSchemaVersion() const final;
0199
0200 bool isSubsetCollection() const final {
0201 return m_isSubsetColl;
0202 }
0203
0204 void setSubsetCollection(bool setSubset=true) final;
0205
0206
0207 Track operator[](std::size_t index) const;
0208
0209 MutableTrack operator[](std::size_t index);
0210
0211 Track at(std::size_t index) const;
0212
0213 MutableTrack at(std::size_t index);
0214
0215
0216
0217 void push_back(const MutableTrack& object);
0218
0219 void push_back(const Track& object);
0220
0221 void prepareForWrite() const final;
0222 void prepareAfterRead() final;
0223 bool setReferences(const podio::ICollectionProvider* collectionProvider) final;
0224
0225
0226 podio::CollectionWriteBuffers getBuffers() final;
0227
0228 void setID(uint32_t ID) final {
0229 m_collectionID = ID;
0230 if (!m_isSubsetColl) {
0231 std::for_each(m_storage.entries.begin(), m_storage.entries.end(),
0232 [ID] (TrackObj* obj) { obj->id = {obj->id.index, static_cast<uint32_t>(ID)}; }
0233 );
0234 }
0235 }
0236
0237 uint32_t getID() const final {
0238 return m_collectionID;
0239 }
0240
0241
0242 bool hasID() const final {
0243 return getID() != static_cast<uint32_t>(podio::ObjectID::untracked) &&
0244 getID() != static_cast<uint32_t>(podio::ObjectID::invalid);
0245 }
0246
0247 [[deprecated("isValid will be removed, use hasID() if you want to check if it has an ID, otherwise assume the "
0248 "collection is valid")]]
0249 bool isValid() const final {
0250 return hasID();
0251 }
0252
0253 size_t getDatamodelRegistryIndex() const final;
0254
0255
0256 iterator begin() {
0257 return iterator(0, &m_storage.entries);
0258 }
0259 const_iterator begin() const {
0260 return const_iterator(0, &m_storage.entries);
0261 }
0262 const_iterator cbegin() const {
0263 return begin();
0264 }
0265 iterator end() {
0266 return iterator(m_storage.entries.size(), &m_storage.entries);
0267 }
0268 const_iterator end() const {
0269 return const_iterator(m_storage.entries.size(), &m_storage.entries);
0270 }
0271 const_iterator cend() const {
0272 return end();
0273 }
0274
0275 reverse_iterator rbegin() {
0276 return reverse_iterator(end());
0277 }
0278 const_reverse_iterator rbegin() const {
0279 return const_reverse_iterator(end());
0280 }
0281 const_reverse_iterator crbegin() const {
0282 return rbegin();
0283 }
0284 reverse_iterator rend() {
0285 return reverse_iterator(begin());
0286 }
0287 const_reverse_iterator rend() const {
0288 return const_reverse_iterator(begin());
0289 }
0290 const_reverse_iterator crend() const {
0291 return rend();
0292 }
0293
0294
0295 std::vector<std::int32_t> type(const size_t nElem = 0) const;
0296 std::vector<edm4hep::Vector3f> position(const size_t nElem = 0) const;
0297 std::vector<edm4hep::Vector3f> momentum(const size_t nElem = 0) const;
0298 std::vector<edm4eic::Cov6f> positionMomentumCovariance(const size_t nElem = 0) const;
0299 std::vector<float> time(const size_t nElem = 0) const;
0300 std::vector<float> timeError(const size_t nElem = 0) const;
0301 std::vector<float> charge(const size_t nElem = 0) const;
0302 std::vector<float> chi2(const size_t nElem = 0) const;
0303 std::vector<std::uint32_t> ndf(const size_t nElem = 0) const;
0304 std::vector<std::int32_t> pdg(const size_t nElem = 0) const;
0305
0306 private:
0307
0308
0309
0310 friend class TrackCollectionData;
0311
0312 mutable bool m_isPrepared{false};
0313 bool m_isSubsetColl{false};
0314 uint32_t m_collectionID{static_cast<uint32_t>(podio::ObjectID::untracked)};
0315 mutable std::unique_ptr<std::mutex> m_storageMtx{std::make_unique<std::mutex>()};
0316 mutable TrackCollectionData m_storage{};
0317 };
0318
0319 std::ostream& operator<<(std::ostream& o, const TrackCollection& v);
0320
0321 template<typename... Args>
0322 MutableTrack TrackCollection::create(Args&&... args) {
0323 if (m_isSubsetColl) {
0324 throw std::logic_error("Cannot create new elements on a subset collection");
0325 }
0326 auto obj = new TrackObj({static_cast<int>(m_storage.entries.size()), m_collectionID}, {std::forward<Args>(args)...});
0327 m_storage.entries.push_back(obj);
0328
0329
0330 obj->m_measurements = new std::vector<edm4eic::Measurement2D>();
0331 obj->m_tracks = new std::vector<edm4eic::Track>();
0332 m_storage.createRelations(obj);
0333 return MutableTrack(podio::utils::MaybeSharedPtr(obj));
0334 }
0335
0336 #if defined(PODIO_JSON_OUTPUT) && !defined(__CLING__)
0337 void to_json(nlohmann::json& j, const TrackCollection& collection);
0338 #endif
0339
0340 }
0341
0342
0343
0344
0345
0346 #if defined(__clang__)
0347 #pragma clang diagnostic push
0348 #pragma clang diagnostic ignored "-Wunknown-warning-option"
0349 #pragma clang diagnostic ignored "-Wdeprecated-redundant-constexpr-static-def"
0350 #pragma clang diagnostic ignored "-Wdeprecated"
0351 constexpr std::string_view edm4eic::TrackCollection::typeName;
0352 constexpr std::string_view edm4eic::TrackCollection::valueTypeName;
0353 constexpr std::string_view edm4eic::TrackCollection::dataTypeName;
0354 #pragma clang diagnostic pop
0355 #elif defined(__GNUC__)
0356 #pragma GCC diagnostic push
0357 #pragma GCC diagnostic ignored "-Wdeprecated"
0358 constexpr std::string_view edm4eic::TrackCollection::typeName;
0359 constexpr std::string_view edm4eic::TrackCollection::valueTypeName;
0360 constexpr std::string_view edm4eic::TrackCollection::dataTypeName;
0361 #pragma GCC diagnostic pop
0362 #endif
0363
0364
0365 #endif