BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSFieldEMInterpolated4D.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 "BDSFieldEMInterpolated4D.hh"
20#include "BDSInterpolator4D.hh"
21
22#include "G4ThreeVector.hh"
23#include "G4Types.hh"
24
25#include <utility>
26
27BDSFieldEMInterpolated4D::BDSFieldEMInterpolated4D(BDSInterpolator4D* eInterpolatorIn,
28 BDSInterpolator4D* bInterpolatorIn,
29 const G4Transform3D& offset,
30 G4double eScalingIn,
31 G4double bScalingIn):
32 BDSFieldEMInterpolated(eInterpolatorIn, bInterpolatorIn, offset, eScalingIn, bScalingIn),
33 eInterpolator(eInterpolatorIn),
34 bInterpolator(bInterpolatorIn)
35{;}
36
37BDSFieldEMInterpolated4D::~BDSFieldEMInterpolated4D()
38{
39 delete eInterpolator;
40 delete bInterpolator;
41}
42
43std::pair<G4ThreeVector,G4ThreeVector> BDSFieldEMInterpolated4D::GetField(const G4ThreeVector& position,
44 const G4double t) const
45{
46 G4ThreeVector e = eInterpolator->GetInterpolatedValue(position[0],position[1],position[2],t) * EScaling();
47 G4ThreeVector b = bInterpolator->GetInterpolatedValue(position[0],position[1],position[2],t) * BScaling();
48 return std::make_pair(b,e);
49}
BDSInterpolator4D * eInterpolator
E Interpolator the field is based on.
virtual std::pair< G4ThreeVector, G4ThreeVector > GetField(const G4ThreeVector &position, const G4double t=0) const
Return the interpolated field value at a given point.
BDSInterpolator4D * bInterpolator
B Interpolator the field is based on.
Class to provide scaling and a base class pointer for interpolator fields.
G4double EScaling() const
Accessor.
G4double BScaling() const
Accessor.
Interface for all 4D interpolators.
G4ThreeVector GetInterpolatedValue(G4double x, G4double y, G4double z, G4double t) const
Public interface to a 4D interpolator. Returns Geant4 type as that's what will be needed.