| 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 // G4EqMagElectricField implementation
27 //
28 // This is the standard right-hand side for eq
29 //
30 // The only case another is required is when u
31 // frame ... or extending the class to include
32 // e.g., an electric field
33 //
34 // Created: V.Grichine, 10.11.1998
35 // -------------------------------------------
36
37 #include "G4EqMagElectricField.hh"
38 #include "globals.hh"
39 #include "G4PhysicalConstants.hh"
40 #include "G4SystemOfUnits.hh"
41
42 G4EqMagElectricField::G4EqMagElectricField(G4E
43 : G4EquationOfMotion( emField )
44 {
45 }
46
47 G4EqMagElectricField::~G4EqMagElectricField()
48
49 void
50 G4EqMagElectricField::SetChargeMomentumMass(G4
51 G4
52 G4
53 {
54 G4double pcharge = particleCharge.GetCharge
55 fElectroMagCof = eplus*pcharge*c_light ;
56 fMassCof = particleMass*particleMass ;
57 }
58
59 void
60 G4EqMagElectricField::EvaluateRhsGivenB(const
61 const
62
63 {
64 // Components of y:
65 // 0-2 dr/ds,
66 // 3-5 dp/ds - momentum derivatives
67
68 G4double pSquared = y[3]*y[3] + y[4]*y[4] +
69
70 G4double Energy = std::sqrt( pSquared + f
71 G4double cof2 = Energy/c_light ;
72
73 G4double pModuleInverse = 1.0/std::sqrt(pS
74
75 G4double inverse_velocity = Energy * pModul
76
77 G4double cof1 = fElectroMagCof*pModuleI
78
79 dydx[0] = y[3]*pModuleInverse ;
80 dydx[1] = y[4]*pModuleInverse ;
81 dydx[2] = y[5]*pModuleInverse ;
82
83 dydx[3] = cof1*(cof2*Field[3] + (y[4]*Field
84
85 dydx[4] = cof1*(cof2*Field[4] + (y[5]*Field
86
87 dydx[5] = cof1*(cof2*Field[5] + (y[3]*Field
88
89 dydx[6] = 0.;//not used
90
91 // Lab Time of flight
92 //
93 dydx[7] = inverse_velocity;
94
95 return;
96 }
97