packages/PHField/PHField3DCartesian.h

0001 #ifndef PHFIELD_PHFIELD3DCARTESIAN_H
0002 #define PHFIELD_PHFIELD3DCARTESIAN_H
0003
0004 #include "PHField.h"
0005
0006 #include <limits>
0007 #include <map>
0008 #include <set>
0009 #include <string>
0010 #include <tuple>
0011
0012 class PHField3DCartesian : public PHField
0013 {
0014  public:
0015
0016   //! constructor
0017   explicit PHField3DCartesian(const std::string &fname, const float magfield_rescale = 1.0, const float innerradius = 0, const float outerradius = 1.e10, const float size_z = 1.e10);
0018
0019   //! destructor
0020   ~PHField3DCartesian() override;
0021
0022   //! access field value
0023   //! Follow the convention of G4ElectroMagneticField
0024   //! @param[in]  Point   space time coordinate. x, y, z, t in Geant4/CLHEP units
0025   //! @param[out] Bfield  field value. In the case of magnetic field, the order is Bx, By, Bz in in Geant4/CLHEP units
0026   void GetFieldValue(const double Point[4], double *Bfield) const override;
0027
0028   void GetFieldValue_nocache(const double Point[4], double *Bfield) const override;
0029
0030   private:
0031   std::string filename;
0032   double xmin {1000000};
0033   double xmax {-1000000};
0034   double ymin {1000000};
0035   double ymax {-1000000};
0036   double zmin {1000000};
0037   double zmax {-1000000};
0038   double xstepsize {std::numeric_limits<double>::quiet_NaN()};
0039   double ystepsize {std::numeric_limits<double>::quiet_NaN()};
0040   double zstepsize {std::numeric_limits<double>::quiet_NaN()};
0041
0042   // these are updated in a const method
0043   // to cache previous values
0044   mutable double bf[2][2][2][3]{};
0045   mutable double xkey_save {std::numeric_limits<double>::quiet_NaN()};
0046   mutable double ykey_save {std::numeric_limits<double>::quiet_NaN()};
0047   mutable double zkey_save {std::numeric_limits<double>::quiet_NaN()};
0048   mutable int cache_hits {0};
0049   mutable int cache_misses {0};
0050
0051   typedef std::tuple<float, float, float> trio;
0052   std::map<std::tuple<float, float, float>, std::tuple<float, float, float> > fieldmap;
0053   std::set<float> xvals;
0054   std::set<float> yvals;
0055   std::set<float> zvals;
0056 };
0057
0058 #endif