BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSFieldMagMultipoleOuterDualOld.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 "BDSFieldMagMultipoleOuterDualOld.hh"
20#include "BDSFieldMagMultipoleOuterOld.hh"
21
22#include "G4ThreeVector.hh"
23
24BDSFieldMagMultipoleOuterDualOld::BDSFieldMagMultipoleOuterDualOld(G4int orderIn,
25 G4double poleTipRadiusIn,
26 const BDSFieldMag* innerFieldIn,
27 G4bool kPositive,
28 G4double separation,
29 G4bool secondFieldOnLeft,
30 G4double arbitraryScaling):
31 fieldBase(nullptr)
32{
33 fieldBase = new BDSFieldMagMultipoleOuterOld(orderIn, poleTipRadiusIn, innerFieldIn, kPositive, arbitraryScaling);
34 G4double offsetX = secondFieldOnLeft ? -separation : separation;
35 offset = G4ThreeVector(offsetX,0,0);
36}
37
38BDSFieldMagMultipoleOuterDualOld::~BDSFieldMagMultipoleOuterDualOld()
39{
40 delete fieldBase;
41}
42
43G4ThreeVector BDSFieldMagMultipoleOuterDualOld::GetField(const G4ThreeVector& position,
44 const G4double t) const
45{
46 G4ThreeVector aSide = fieldBase->GetField(position, t);
47 G4ThreeVector shiftedPos = position + offset;
48 G4ThreeVector bSide = fieldBase->GetField(shiftedPos,t);
49 bSide *= -1.0; // opposite to whatever 'a' side is
50 G4ThreeVector result = aSide + bSide;
51 return result;
52}
virtual G4ThreeVector GetField(const G4ThreeVector &position, const double t=0) const
Access the field value.
A simple paramaterisation of N-Pole outer yoke magnetic field.
virtual G4ThreeVector GetField(const G4ThreeVector &position, const double t=0) const
Access the field value.
Interface for static magnetic fields that may or may not be local.
Definition: BDSFieldMag.hh:39