|
| 1 | +#ifndef DataFormats_L1Scouting_OrbitFlatTable_h |
| 2 | +#define DataFormats_L1Scouting_OrbitFlatTable_h |
| 3 | + |
| 4 | +/** |
| 5 | + * A cross-breed of a FlatTable and an OrbitCollection |
| 6 | + */ |
| 7 | + |
| 8 | +#include "DataFormats/NanoAOD/interface/FlatTable.h" |
| 9 | + |
| 10 | +#include <cstdint> |
| 11 | +#include <vector> |
| 12 | +#include <string> |
| 13 | +#include <type_traits> |
| 14 | + |
| 15 | +namespace l1ScoutingRun3 { |
| 16 | + |
| 17 | + class OrbitFlatTable : public nanoaod::FlatTable { |
| 18 | + public: |
| 19 | + static constexpr unsigned int NBX = 3564; |
| 20 | + |
| 21 | + OrbitFlatTable() : nanoaod::FlatTable(), bxOffsets_(orbitBufferSize_ + 1, 0) {} |
| 22 | + |
| 23 | + OrbitFlatTable(std::vector<unsigned> bxOffsets, |
| 24 | + const std::string &name, |
| 25 | + bool singleton = false, |
| 26 | + bool extension = false) |
| 27 | + : nanoaod::FlatTable(bxOffsets.back(), name, singleton, extension), bxOffsets_(bxOffsets) { |
| 28 | + if (bxOffsets.size() != orbitBufferSize_ + 1) { |
| 29 | + throw cms::Exception("LogicError") << "Mismatch between bxOffsets.size() " << bxOffsets.size() |
| 30 | + << " and orbitBufferSize_ + 1" << (orbitBufferSize_ + 1); |
| 31 | + } |
| 32 | + } |
| 33 | + |
| 34 | + ~OrbitFlatTable() {} |
| 35 | + |
| 36 | + using FlatTable::nRows; |
| 37 | + using FlatTable::size; |
| 38 | + |
| 39 | + /// number of rows for single BX |
| 40 | + unsigned int nRows(unsigned bx) const { |
| 41 | + if (bx >= orbitBufferSize_) |
| 42 | + throwBadBx(bx); |
| 43 | + return bxOffsets_[bx + 1] - bxOffsets_[bx]; |
| 44 | + }; |
| 45 | + unsigned int size(unsigned bx) const { return nRows(bx); } |
| 46 | + |
| 47 | + /// get a column by index (const) |
| 48 | + template <typename T> |
| 49 | + auto columnData(unsigned int column) const { |
| 50 | + return nanoaod::FlatTable::columnData<T>(column); |
| 51 | + } |
| 52 | + |
| 53 | + /// get a column by index and bx (const) |
| 54 | + template <typename T> |
| 55 | + auto columnData(unsigned int column, unsigned bx) const { |
| 56 | + if (bx >= orbitBufferSize_) |
| 57 | + throwBadBx(bx); |
| 58 | + auto begin = beginData<T>(column); |
| 59 | + return std::span(begin + bxOffsets_[bx], begin + bxOffsets_[bx + 1]); |
| 60 | + } |
| 61 | + |
| 62 | + /// get a column by index (non-const) |
| 63 | + template <typename T> |
| 64 | + auto columnData(unsigned int column) { |
| 65 | + return nanoaod::FlatTable::columnData<T>(column); |
| 66 | + } |
| 67 | + |
| 68 | + /// get a column by index and bx (non-const) |
| 69 | + template <typename T> |
| 70 | + auto columnData(unsigned int column, unsigned bx) { |
| 71 | + if (bx >= orbitBufferSize_) |
| 72 | + throwBadBx(bx); |
| 73 | + auto begin = beginData<T>(column); |
| 74 | + return std::span(begin + bxOffsets_[bx], begin + bxOffsets_[bx + 1]); |
| 75 | + } |
| 76 | + |
| 77 | + /// get a column value for singleton (const) |
| 78 | + template <typename T> |
| 79 | + const auto &columValue(unsigned int column, unsigned bx) const { |
| 80 | + if (!singleton()) |
| 81 | + throw cms::Exception("LogicError", "columnValue works only for singleton tables"); |
| 82 | + if (bx >= orbitBufferSize_ || bxOffsets_[bx + 1] == bxOffsets_[bx]) |
| 83 | + throwBadBx(bx); |
| 84 | + auto begin = beginData<T>(column); |
| 85 | + return *(begin + bxOffsets_[bx]); |
| 86 | + } |
| 87 | + |
| 88 | + private: |
| 89 | + std::vector<unsigned> bxOffsets_; |
| 90 | + |
| 91 | + // there are 3564 BX in one orbtit [1,3564], one extra |
| 92 | + // count added to keep first entry of the vector |
| 93 | + static constexpr int orbitBufferSize_ = NBX + 1; |
| 94 | + |
| 95 | + [[noreturn]] void throwBadBx(unsigned bx) const { |
| 96 | + throw cms::Exception("OrbitFlatTable") << "Trying to access bad bx " << bx; |
| 97 | + } |
| 98 | + }; |
| 99 | + |
| 100 | +} // namespace l1ScoutingRun3 |
| 101 | + |
| 102 | +#endif |
0 commit comments