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