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