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