| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // G4CachedMagneticField implementation
27 //
28 // Author: J.Apostolakis, 20 July 2009.
29 // -------------------------------------------
30
31 #include "G4CachedMagneticField.hh"
32
33 G4CachedMagneticField::G4CachedMagneticField(G
34 G
35 : fpMagneticField(pMagField), fDistanceConst
36 fLastLocation(DBL_MAX,DBL_MAX,DBL_MAX), fL
37 {
38 ClearCounts();
39 }
40
41 G4Field* G4CachedMagneticField::Clone() const
42 {
43 // Cannot use copy constructor: need to clon
44
45 auto aF = static_cast<G4MagneticField*>(fpM
46 auto cloned = new G4CachedMagneticField(aF,
47
48 cloned->fLastLocation = fLastLocation;
49 cloned->fLastValue = fLastValue;
50 return cloned;
51 }
52
53 G4CachedMagneticField::~G4CachedMagneticField(
54
55 void
56 G4CachedMagneticField::ReportStatistics()
57 {
58 G4cout << " Cached field: " << G4endl
59 << " Number of calls: " << f
60 << " Number of evaluations : " << f
61 }
62
63 G4CachedMagneticField::
64 G4CachedMagneticField(const G4CachedMagneticFi
65 : G4MagneticField(rightCMF)
66 {
67 fpMagneticField= rightCMF.fpMagneticField;
68 fDistanceConst = rightCMF.fDistanceConst;
69 fLastLocation = rightCMF.fLastLocation;
70 fLastValue = rightCMF.fLastValue;
71 ClearCounts();
72 }
73
74 G4CachedMagneticField&
75 G4CachedMagneticField::operator = (const G4Cac
76 {
77 if (&p == this) { return *this; }
78 G4MagneticField::operator=(p);
79 fpMagneticField= p.fpMagneticField; // NOTE
80 fDistanceConst = p.fDistanceConst;
81 fLastLocation = p.fLastLocation;
82 fLastValue = p.fLastValue;
83 ClearCounts();
84 return *this;
85 }
86
87 void
88 G4CachedMagneticField::GetFieldValue( const G4
89 G4
90 {
91 G4ThreeVector newLocation( Point[0], Point[1
92
93 G4double distSq= (newLocation-fLastLoca
94 ++fCountCalls;
95 if( distSq < fDistanceConst*fDistanceConst )
96 {
97 Bfield[0] = fLastValue.x();
98 Bfield[1] = fLastValue.y();
99 Bfield[2] = fLastValue.z();
100 }
101 else
102 {
103 fpMagneticField->GetFieldValue( Point, Bf
104 ++fCountEvaluations;
105 fLastLocation = G4ThreeVector( Point[0],
106 fLastValue = G4ThreeVector( Bfield[0],
107 }
108 }
109