Project
Loading...
Searching...
No Matches
ArrowTypes.h
Go to the documentation of this file.
1// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
2// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
3// All rights not expressly granted are reserved.
4//
5// This software is distributed under the terms of the GNU General Public
6// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
7//
8// In applying this license CERN does not waive the privileges and immunities
9// granted to it by virtue of its status as an Intergovernmental Organization
10// or submit itself to any jurisdiction.
11
12#ifndef O2_FRAMEWORK_ARROWTYPES_H
13#define O2_FRAMEWORK_ARROWTYPES_H
14#include <arrow/table.h>
15#include "Framework/Traits.h"
16#include "arrow/type_fwd.h"
17#include <span>
18
19namespace o2::soa
20{
21struct ArrowRange {
24
25 bool operator!=(ArrowRange const& other) const
26 {
27 return (offset != other.offset) || (size != other.size);
28 }
29};
30
32 std::shared_ptr<arrow::Table> tablePtr = nullptr;
34
35 ArrowTableRef() = default;
36 ArrowTableRef(std::shared_ptr<arrow::Table> table)
37 : tablePtr{table},
38 range{0, table->num_rows()}
39 {
40 }
41 ArrowTableRef(std::shared_ptr<arrow::Table> table, ArrowRange range_)
42 : tablePtr{table},
43 range{range_}
44 {
45 }
46
48 {
49 return {tablePtr, {0, 0}};
50 }
51
53 {
54 return {tablePtr, newRange};
55 }
56
57 std::shared_ptr<arrow::Table> const& operator->() const
58 {
59 return tablePtr;
60 }
61};
62
63template <typename T>
65};
66template <>
68 using type = arrow::BooleanArray;
69};
70template <>
72 using type = arrow::Int8Array;
73};
74template <>
76 using type = arrow::UInt8Array;
77};
78template <>
79struct arrow_array_for<int16_t> {
80 using type = arrow::Int16Array;
81};
82template <>
84 using type = arrow::UInt16Array;
85};
86template <>
88 using type = arrow::Int32Array;
89};
90template <>
92 using type = arrow::Int64Array;
93};
94template <>
96 using type = arrow::UInt32Array;
97};
98template <>
100 using type = arrow::UInt64Array;
101};
102template <>
104 using type = arrow::FloatArray;
105};
106template <>
107struct arrow_array_for<double> {
108 using type = arrow::DoubleArray;
109};
110template <>
111struct arrow_array_for<std::span<std::byte>> {
112 using type = arrow::BinaryViewArray;
113};
114template <int N>
116 using type = arrow::FixedSizeListArray;
118};
119template <int N>
120struct arrow_array_for<int[N]> {
121 using type = arrow::FixedSizeListArray;
122 using value_type = int;
123};
124template <int N>
125struct arrow_array_for<short[N]> {
126 using type = arrow::FixedSizeListArray;
127 using value_type = short;
128};
129template <int N>
130struct arrow_array_for<double[N]> {
131 using type = arrow::FixedSizeListArray;
132 using value_type = double;
133};
134template <int N>
136 using type = arrow::FixedSizeListArray;
138};
139template <int N>
141 using type = arrow::FixedSizeListArray;
143};
144
145#define ARROW_VECTOR_FOR(_type_) \
146 template <> \
147 struct arrow_array_for<std::vector<_type_>> { \
148 using type = arrow::ListArray; \
149 using value_type = _type_; \
150 };
151
156
161
164
165template <typename T>
167template <typename T>
169
170template <class Array>
171using array_element_t = std::decay_t<decltype(std::declval<Array>()[0])>;
172
173template <typename T>
174std::shared_ptr<arrow::DataType> asArrowDataType(int list_size = 1)
175{
176 auto typeGenerator = [](std::shared_ptr<arrow::DataType> const& type, int list_size) -> std::shared_ptr<arrow::DataType> {
177 switch (list_size) {
178 case -1:
179 return arrow::list(type);
180 case 1:
181 return std::move(type);
182 default:
183 return arrow::fixed_size_list(type, list_size);
184 }
185 };
186
187 if constexpr (std::is_arithmetic_v<T>) {
188 if constexpr (std::same_as<T, bool>) {
189 return typeGenerator(arrow::boolean(), list_size);
190 } else if constexpr (std::same_as<T, uint8_t>) {
191 return typeGenerator(arrow::uint8(), list_size);
192 } else if constexpr (std::same_as<T, uint16_t>) {
193 return typeGenerator(arrow::uint16(), list_size);
194 } else if constexpr (std::same_as<T, uint32_t>) {
195 return typeGenerator(arrow::uint32(), list_size);
196 } else if constexpr (std::same_as<T, uint64_t>) {
197 return typeGenerator(arrow::uint64(), list_size);
198 } else if constexpr (std::same_as<T, int8_t>) {
199 return typeGenerator(arrow::int8(), list_size);
200 } else if constexpr (std::same_as<T, int16_t>) {
201 return typeGenerator(arrow::int16(), list_size);
202 } else if constexpr (std::same_as<T, int32_t>) {
203 return typeGenerator(arrow::int32(), list_size);
204 } else if constexpr (std::same_as<T, int64_t>) {
205 return typeGenerator(arrow::int64(), list_size);
206 } else if constexpr (std::same_as<T, float>) {
207 return typeGenerator(arrow::float32(), list_size);
208 } else if constexpr (std::same_as<T, double>) {
209 return typeGenerator(arrow::float64(), list_size);
210 }
211 } else if constexpr (std::is_bounded_array_v<T>) {
212 return asArrowDataType<array_element_t<T>>(std::extent_v<T>);
213 } else if constexpr (o2::framework::is_specialization_v<T, std::vector>) {
214 return asArrowDataType<typename T::value_type>(-1);
215 }
216 return nullptr;
217}
218} // namespace o2::soa
219#endif // O2_FRAMEWORK_ARROWTYPES_H
#define ARROW_VECTOR_FOR(_type_)
Definition ArrowTypes.h:145
unsigned short int uint16_t
Timeframe settings Average distance in cm of events merged into Randomize distance around average distance of merged events Simulate correct bunch interactions instead of placing only the average number of events A value[n] sets TFlen for[n] collisions in Number of bunches per trainsort Number of bunch trains Instantaneous interaction rate Do not simulate border Total number of collisions to be placed in the interior of all time Do not select random but walk over array of events in stride steps h Event generator settings h configStandalone Use CPU for GPU Force usage of the specified GPU device no CPU fallback Continue processing after an error First event to Number of events to Number of iterations to Number of iterations to Number of initial iterations excluded from const char Directory with events to Events directory is and not inside events Show standalone event Run event generator maximum time bin of continuous for triggered for automatic continuous for automatic continuous triggered uint64_t
configStandalone Reconstruction rec_tpc multiplier for the combined cluster track error during track following area in cm for the search of for z only used if decrease removal by when above this local occupancy rowx Shift of TPC revert cluster flag edge bit distance to edge exceeds this ZOffset for secondary multiple of radius max P of electron K BB bands crossing max P of electron p BB bands crossing max P of electron d BB bands crossing max P of electron t BB bands crossing min accepted dEdxMaxTPC of the track Additive extra cluster error for Y2 if edge flag set Additive extra cluster error for Y2 if single set Additive extra cluster error for Y2 if or single set Additive extra cluster error for Z2 if or single set Extra errors to add to track when crossing depending on addErrorsCECrossing Inner portion of y range in sector that is not used in searching for extrapolated track candidates factor2YS for track merging factor2K for track merging Margin in cm of Y position when rejecting edge clusters based on uncorrected track Y Max fraction of shared clusters for tracklet Except pad from cluster finding if total number of charges in a fragment is above this Except pad from cluster finding if number of consecutive charges in a fragment is above this Threshold where a timebin is considered disabling the noisy pad check for that pad Treat signal above this minimum as saturated Smoothing factor for the exponential Highly Ionising Particle tail filter Number of timebins to use for the averaging Add additional custom track errors when crossing relative as fraction of Do N fit passes in final fit of< 0:reject at max -n hits") AddOptionRTC(rejectIFCLowRadiusCluster, uint8_t, 1, "", 0, "Reject clusters that get the IFC mask error during refit") AddOptionRTC(dEdxTruncLow, uint8_t, 2, "", 0, "Low truncation threshold, fraction of 128") AddOptionRTC(dEdxTruncHigh, uint8_t, 77, "", 0, "High truncation threshold, fraction of 128") AddOptionRTC(extrapolationTracking, int8_t, 1, "", 0, "Enable Extrapolation Tracking(prolong tracks to adjacent sectors to find short segments)") AddOptionRTC(disableRefitAttachment, uint8_t, 0, "", 0, "Bitmask to disable certain attachment steps during refit(1:attachment, 2:propagation, 4:loop following, 8:mirroring)") AddOptionRTC(rejectionStrategy, uint8_t, o2::gpu::GPUSettings::RejectionStrategyA, "", 0, "Enable rejection of TPC clusters for compression(0=no, 1=strategy A, 2=strategy B)") AddOptionRTC(mergeLoopersAfterburner, uint8_t, 1, "", 0, "Run afterburner for additional looper merging") AddOptionRTC(compressionTypeMask, uint8_t, o2::gpu::GPUSettings::CompressionFull, "", 0, "TPC Compression mode bits(1=truncate charge/width LSB, 2=differences, 4=track-model)") AddOptionRTC(compressionSortOrder, uint8_t, o2::gpu::GPUSettings::SortTime, "", 0, "Sort order of TPC compression(0=time, 1=pad, 2=Z-time-pad, 3=Z-pad-time, 4=no sorting(use incoming order))") AddOptionRTC(sigBitsCharge, uint8_t, 4, "", 0, "Number of significant bits for TPC cluster charge in compression mode 1") AddOptionRTC(sigBitsWidth, uint8_t, 3, "", 0, "Number of significant bits for TPC cluster width in compression mode 1") AddOptionRTC(dropLoopers, uint8_t, 0, "", 0, "Drop looping tracks starting from second loop") AddOptionRTC(mergerInterpolateErrors, uint8_t, 1, "", 0, "Use interpolation instead of extrapolation for chi2 based cluster rejection") AddOptionRTC(mergerInterpolateRejectAlsoOnCurrentPosition, uint8_t, 1, "", 0, "When using mergerInterpolateErrors, reject based on chi2 twice computed with interpolated and current track position starting from NDF > mergerNonInterpolateRejectMinNDF Merge tracks accross the central electrode Enable PID response Add sub threshold clusters in TPC dEdx computation OR mask of TPC flags that will reject the cluster in alternative dEdx Reject edge clusters based on uncorrected track Y during track fit Increase decrease removal by when on this number of pads from the edge h configStandalone Reconstruction rec_trd Max chi2 for TRD tracklets to be matched to a track Chi2 cut for strict matching mode Number of sigmas for ITS TPC track time error estimate Addition to search road around track prolongation along Y in cm Optimal tracklet rphi resolution in maximum chi2 per attached tracklet for TRD tracks take into account the full likelihood of the slope including exponential tails Tracks with less attached tracklets are discarded after the tracking Post trigger Pile up integration time in BCs fits all iterations but the last one in projections Use material lookup table for TPC refit Apply workarounds to avoid cyclic merge should not be needed tpc dyn configStandalone Processing proc_rtc Replace constant variables static by constexpr expressions Compile RTC in deterministic with NO_FAST_MATH flags and GPUCA_DETERMINISTIC_MODE define Use RTC to optimize GPU code h configStandalone Processing proc_rtctech Use a file lock to serialize access to the cache folder int32_t
configStandalone Reconstruction rec_tpc multiplier for the combined cluster track error during track following area in cm for the search of for z only used if decrease removal by when above this local occupancy rowx Shift of TPC revert cluster flag edge bit distance to edge exceeds this ZOffset for secondary multiple of radius max P of electron K BB bands crossing max P of electron p BB bands crossing max P of electron d BB bands crossing max P of electron t BB bands crossing min accepted dEdxMaxTPC of the track Additive extra cluster error for Y2 if edge flag set Additive extra cluster error for Y2 if single set Additive extra cluster error for Y2 if or single set Additive extra cluster error for Z2 if or single set Extra errors to add to track when crossing depending on addErrorsCECrossing Inner portion of y range in sector that is not used in searching for extrapolated track candidates factor2YS for track merging factor2K for track merging Margin in cm of Y position when rejecting edge clusters based on uncorrected track Y Max fraction of shared clusters for tracklet Except pad from cluster finding if total number of charges in a fragment is above this Except pad from cluster finding if number of consecutive charges in a fragment is above this Threshold where a timebin is considered disabling the noisy pad check for that pad Treat signal above this minimum as saturated Smoothing factor for the exponential Highly Ionising Particle tail filter Number of timebins to use for the averaging Add additional custom track errors when crossing relative as fraction of Do N fit passes in final fit of< 0:reject at max -n hits") AddOptionRTC(rejectIFCLowRadiusCluster, uint8_t, 1, "", 0, "Reject clusters that get the IFC mask error during refit") AddOptionRTC(dEdxTruncLow, uint8_t, 2, "", 0, "Low truncation threshold, fraction of 128") AddOptionRTC(dEdxTruncHigh, uint8_t, 77, "", 0, "High truncation threshold, fraction of 128") AddOptionRTC(extrapolationTracking, int8_t, 1, "", 0, "Enable Extrapolation Tracking(prolong tracks to adjacent sectors to find short segments)") AddOptionRTC(disableRefitAttachment, uint8_t, 0, "", 0, "Bitmask to disable certain attachment steps during refit(1:attachment, 2:propagation, 4:loop following, 8:mirroring)") AddOptionRTC(rejectionStrategy, uint8_t, o2::gpu::GPUSettings::RejectionStrategyA, "", 0, "Enable rejection of TPC clusters for compression(0=no, 1=strategy A, 2=strategy B)") AddOptionRTC(mergeLoopersAfterburner, uint8_t, 1, "", 0, "Run afterburner for additional looper merging") AddOptionRTC(compressionTypeMask, uint8_t, o2::gpu::GPUSettings::CompressionFull, "", 0, "TPC Compression mode bits(1=truncate charge/width LSB, 2=differences, 4=track-model)") AddOptionRTC(compressionSortOrder, uint8_t, o2::gpu::GPUSettings::SortTime, "", 0, "Sort order of TPC compression(0=time, 1=pad, 2=Z-time-pad, 3=Z-pad-time, 4=no sorting(use incoming order))") AddOptionRTC(sigBitsCharge, uint8_t, 4, "", 0, "Number of significant bits for TPC cluster charge in compression mode 1") AddOptionRTC(sigBitsWidth, uint8_t, 3, "", 0, "Number of significant bits for TPC cluster width in compression mode 1") AddOptionRTC(dropLoopers, uint8_t, 0, "", 0, "Drop looping tracks starting from second loop") AddOptionRTC(mergerInterpolateErrors, uint8_t, 1, "", 0, "Use interpolation instead of extrapolation for chi2 based cluster rejection") AddOptionRTC(mergerInterpolateRejectAlsoOnCurrentPosition, uint8_t, 1, "", 0, "When using mergerInterpolateErrors, reject based on chi2 twice computed with interpolated and current track position starting from NDF > mergerNonInterpolateRejectMinNDF Merge tracks accross the central electrode Enable PID response Add sub threshold clusters in TPC dEdx computation OR mask of TPC flags that will reject the cluster in alternative dEdx Reject edge clusters based on uncorrected track Y during track fit Increase decrease removal by when on this number of pads from the edge h configStandalone Reconstruction rec_trd Max chi2 for TRD tracklets to be matched to a track Chi2 cut for strict matching mode Number of sigmas for ITS TPC track time error estimate Addition to search road around track prolongation along Y in cm Optimal tracklet rphi resolution in maximum chi2 per attached tracklet for TRD tracks take into account the full likelihood of the slope including exponential tails Tracks with less attached tracklets are discarded after the tracking Post trigger Pile up integration time in BCs fits all iterations but the last one in projections Use material lookup table for TPC refit Apply workarounds to avoid cyclic merge should not be needed tpc dyn configStandalone Processing proc_rtc bool
configStandalone Reconstruction rec_tpc float
GLsizeiptr size
Definition glcorearb.h:659
GLenum GLint * range
Definition glcorearb.h:1899
GLint GLint GLsizei GLint GLenum GLenum type
Definition glcorearb.h:275
GLintptr offset
Definition glcorearb.h:660
std::decay_t< decltype(std::declval< Array >()[0])> array_element_t
Definition ArrowTypes.h:171
std::shared_ptr< arrow::DataType > asArrowDataType(int list_size=1)
Definition ArrowTypes.h:174
typename arrow_array_for< T >::value_type value_for_t
Definition ArrowTypes.h:168
typename arrow_array_for< T >::type arrow_array_for_t
Definition ArrowTypes.h:166
bool operator!=(ArrowRange const &other) const
Definition ArrowTypes.h:25
std::shared_ptr< arrow::Table > const & operator->() const
Definition ArrowTypes.h:57
ArrowTableRef makeEmpty() const
Definition ArrowTypes.h:47
std::shared_ptr< arrow::Table > tablePtr
Definition ArrowTypes.h:32
ArrowTableRef slice(ArrowRange newRange) const
Definition ArrowTypes.h:52
ArrowTableRef(std::shared_ptr< arrow::Table > table, ArrowRange range_)
Definition ArrowTypes.h:41
ArrowTableRef(std::shared_ptr< arrow::Table > table)
Definition ArrowTypes.h:36
arrow::FixedSizeListArray type
Definition ArrowTypes.h:131
arrow::FixedSizeListArray type
Definition ArrowTypes.h:116
arrow::FixedSizeListArray type
Definition ArrowTypes.h:141
arrow::FixedSizeListArray type
Definition ArrowTypes.h:136
arrow::FixedSizeListArray type
Definition ArrowTypes.h:121
arrow::FixedSizeListArray type
Definition ArrowTypes.h:126
VectorOfTObjectPtrs other