File indexing completed on 2025-07-18 08:15:07
0001
0002
0003
0004
0005
0006
0007
0008
0009 #pragma once
0010
0011 #include <Acts/Definitions/Algebra.hpp>
0012 #include <Acts/MagneticField/MagneticFieldContext.hpp>
0013 #include <Acts/MagneticField/MagneticFieldProvider.hpp>
0014 #include <Acts/Utilities/Result.hpp>
0015 #include <DD4hep/Detector.h>
0016 #include <gsl/pointers>
0017 #include <memory>
0018 #include <utility>
0019 #include <variant>
0020
0021 namespace eicrecon::BField {
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033 class DD4hepBField final : public Acts::MagneticFieldProvider {
0034 public:
0035 gsl::not_null<const dd4hep::Detector*> m_det;
0036
0037 public:
0038 struct Cache {
0039 Cache(const Acts::MagneticFieldContext& ) {}
0040 };
0041
0042 Acts::MagneticFieldProvider::Cache
0043 makeCache(const Acts::MagneticFieldContext& mctx) const override {
0044 return Acts::MagneticFieldProvider::Cache(std::in_place_type<Cache>, mctx);
0045 }
0046
0047
0048
0049
0050
0051 explicit DD4hepBField(gsl::not_null<const dd4hep::Detector*> det) : m_det(det) {}
0052
0053
0054
0055
0056
0057
0058
0059
0060
0061
0062 Acts::Result<Acts::Vector3> getField(const Acts::Vector3& position,
0063 Acts::MagneticFieldProvider::Cache& cache) const override;
0064
0065 #if Acts_VERSION_MAJOR < 39
0066
0067
0068
0069
0070
0071
0072
0073
0074
0075
0076
0077
0078
0079 Acts::Result<Acts::Vector3>
0080 getFieldGradient(const Acts::Vector3& position, Acts::ActsMatrix<3, 3>& ,
0081 Acts::MagneticFieldProvider::Cache& cache) const override;
0082 #endif
0083 };
0084
0085 using BFieldVariant = std::variant<std::shared_ptr<const DD4hepBField>>;
0086
0087 }