BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSFieldMagGlobalPlacement.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 "BDSFieldMag.hh"
20#include "BDSFieldMagGlobalPlacement.hh"
21
22#include "G4ThreeVector.hh"
23#include "G4Types.hh"
24
25BDSFieldMagGlobalPlacement::BDSFieldMagGlobalPlacement(BDSFieldMag* fieldIn):
26 field(fieldIn)
27{
28 finiteStrength = field->FiniteStrength();
29}
30
31BDSFieldMagGlobalPlacement::~BDSFieldMagGlobalPlacement()
32{
33 delete field;
34}
35
36G4ThreeVector BDSFieldMagGlobalPlacement::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 BDSFieldMagGlobalPlacement::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}
virtual G4ThreeVector GetFieldTransformed(const G4ThreeVector &position, const G4double t) const
virtual G4ThreeVector GetField(const G4ThreeVector &position, const G4double t=0) const
BDSFieldMag * field
The field on which this is based.
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
G4ThreeVector ConvertToLocal(const G4ThreeVector &globalPosition, G4bool &foundAPlacementVolume) const
G4ThreeVector ConvertAxisToGlobal(const G4ThreeVector &localAxis) const