BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSFieldEMGlobalPlacement.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 "BDSFieldEM.hh"
20#include "BDSFieldEMGlobalPlacement.hh"
21
22#include "G4ThreeVector.hh"
23#include "G4Types.hh"
24
25#include <utility>
26
27BDSFieldEMGlobalPlacement::BDSFieldEMGlobalPlacement(BDSFieldEM* fieldIn):
28 field(fieldIn)
29{
30 finiteStrength = field->FiniteStrength();
31}
32
33BDSFieldEMGlobalPlacement::~BDSFieldEMGlobalPlacement()
34{
35 delete field;
36}
37
38std::pair<G4ThreeVector,G4ThreeVector> BDSFieldEMGlobalPlacement::GetFieldTransformed(const G4ThreeVector& position,
39 const G4double t) const
40{
41 if (!finiteStrength)
42 {return std::make_pair(G4ThreeVector(),G4ThreeVector());} // quicker than query
43 else
44 {return GetField(position, t);}
45}
46
47std::pair<G4ThreeVector,G4ThreeVector> BDSFieldEMGlobalPlacement::GetField(const G4ThreeVector& position,
48 const G4double t) const
49{
50 G4bool foundAPlacementVolume = false;
51 G4ThreeVector localPosition = ConvertToLocal(position, foundAPlacementVolume);
52 if (!foundAPlacementVolume)
53 {return std::make_pair(G4ThreeVector(),G4ThreeVector());} // force to 0 outside range of coordinate transforms
54 auto localField = field->GetFieldTransformed(localPosition, t);
55 auto globalField = ConvertAxisToGlobal(localField);
56 return globalField;
57}
virtual std::pair< G4ThreeVector, G4ThreeVector > GetField(const G4ThreeVector &position, const G4double t) const
virtual std::pair< G4ThreeVector, G4ThreeVector > GetFieldTransformed(const G4ThreeVector &position, const G4double t) const
BDSFieldEM * field
The field on which this is based.
Interface for BDSIM electro-magnetic fields that may or may not be local.
Definition: BDSFieldEM.hh:41
virtual std::pair< G4ThreeVector, G4ThreeVector > GetFieldTransformed(const G4ThreeVector &position, const G4double t) const
Get the field value after applying transform for local offset.
Definition: BDSFieldEM.cc:44
G4bool finiteStrength
Flag to cache whether finite nor not.
Definition: BDSFieldEM.hh:88
G4ThreeVector ConvertToLocal(const G4ThreeVector &globalPosition, G4bool &foundAPlacementVolume) const
G4ThreeVector ConvertAxisToGlobal(const G4ThreeVector &localAxis) const