19#include "BDSFieldEM.hh"
22#include "G4ThreeVector.hh"
23#include "G4Transform3D.hh"
29 transform(G4Transform3D::Identity),
30 inverseTransform(G4Transform3D::Identity)
35 transform(transformIn),
36 inverseTransform(transformIn.inverse())
40 const G4double t)
const
43 {
return std::make_pair(G4ThreeVector(), G4ThreeVector());}
44 else if (
transform != G4Transform3D::Identity)
46 G4ThreeVector transformedPosition =
transform * (HepGeom::Point3D<G4double>)position;
47 auto field =
GetField(transformedPosition, t);
48 G4ThreeVector transformedBField =
transform * (HepGeom::Vector3D<G4double>)field.first;
49 G4ThreeVector transformedEField =
transform * (HepGeom::Vector3D<G4double>)field.second;
50 return std::make_pair(transformedBField, transformedEField);
57 G4double* field)
const
59 auto fieldValue =
GetFieldTransformed(G4ThreeVector(point[0], point[1], point[2]), point[3]);
60 field[0] = fieldValue.first[0];
61 field[1] = fieldValue.first[1];
62 field[2] = fieldValue.first[2];
63 field[3] = fieldValue.second[0];
64 field[4] = fieldValue.second[1];
65 field[5] = fieldValue.second[2];
G4Transform3D inverseTransform
The complimentary transform used to initially rotate the point of query.
virtual void SetTransform(const G4Transform3D &transformIn)
virtual std::pair< G4ThreeVector, G4ThreeVector > GetField(const G4ThreeVector &position, const G4double t=0) const =0
G4Transform3D transform
Transform to apply for the field relative to the local coordinates of the geometry.
virtual void GetFieldValue(const G4double point[4], G4double *field) const
virtual std::pair< G4ThreeVector, G4ThreeVector > GetFieldTransformed(const G4ThreeVector &position, const G4double t) const
Get the field value after applying transform for local offset.
G4bool finiteStrength
Flag to cache whether finite nor not.