BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSFieldEInterpolated1D.cc
1/*
2Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway,
3University of London 2001 - 2023.
4
5This file is part of BDSIM.
6
7BDSIM is free software: you can redistribute it and/or modify
8it under the terms of the GNU General Public License as published
9by the Free Software Foundation version 3 of the License.
10
11BDSIM is distributed in the hope that it will be useful, but
12WITHOUT ANY WARRANTY; without even the implied warranty of
13MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14GNU General Public License for more details.
15
16You should have received a copy of the GNU General Public License
17along with BDSIM. If not, see <http://www.gnu.org/licenses/>.
18*/
19#include "BDSDimensionType.hh"
20#include "BDSFieldEInterpolated1D.hh"
21#include "BDSInterpolator1D.hh"
22
23#include "G4ThreeVector.hh"
24
25BDSFieldEInterpolated1D::BDSFieldEInterpolated1D(BDSInterpolator1D* interpolatorIn,
26 const G4Transform3D& offset,
27 G4double eScalingIn):
28 BDSFieldEInterpolated(interpolatorIn, offset, eScalingIn),
29 interpolator(interpolatorIn),
30 dimensionIndex((interpolatorIn->FirstDimension()).underlying()),
31 time((interpolatorIn->FirstDimension()).underlying() > 2)
32{;}
33
34BDSFieldEInterpolated1D::~BDSFieldEInterpolated1D()
35{
36 delete interpolator;
37}
38
39G4ThreeVector BDSFieldEInterpolated1D::GetField(const G4ThreeVector& position,
40 const G4double t) const
41{
42 G4double coordinate = 0;
43 if (time)
44 {coordinate = t;}
45 else
46 {coordinate = position[dimensionIndex];}
47 return interpolator->GetInterpolatedValue(coordinate) * EScaling();
48}
virtual G4ThreeVector GetField(const G4ThreeVector &position, const G4double t=0) const
Return the interpolated field value at a given point.
const G4bool time
Cache of whether to use time coordinate.
const G4int dimensionIndex
Integer index to dimension to use.
BDSInterpolator1D * interpolator
Interpolator the field is based on.
Class to provide scaling and a base class pointer for interpolator fields.
Interface for all 1D interpolators.
G4ThreeVector GetInterpolatedValue(G4double x) const
Public interface to a 1D interpolator. Returns Geant4 type as that's what will be needed.