| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 // 23 //
27 // ------------------------------------------- 24 // -------------------------------------------------------------------
28 // GEANT 4 class header file 25 // GEANT 4 class header file
29 // 26 //
30 // CERN, Geneva, Switzerland 27 // CERN, Geneva, Switzerland
31 // 28 //
32 // File name: G4KM_NucleonEqRhs.hh 29 // File name: G4KM_NucleonEqRhs.hh
33 // 30 //
34 // Author: Alessandro Brunengo (Al 31 // Author: Alessandro Brunengo (Alessandro.Brunengo@ge.infn.it)
35 // 32 //
36 // Creation date: 5 June 2000 33 // Creation date: 5 June 2000
37 // ------------------------------------------- 34 // -------------------------------------------------------------------
38 35
39 #ifndef G4KM_NucleonEqRhs_hh 36 #ifndef G4KM_NucleonEqRhs_hh
40 #define G4KM_NucleonEqRhs_hh 37 #define G4KM_NucleonEqRhs_hh
41 38
42 #include "globals.hh" 39 #include "globals.hh"
43 //#include "G4EquationOfMotion.hh" 40 //#include "G4EquationOfMotion.hh"
44 #include "G4Mag_EqRhs.hh" 41 #include "G4Mag_EqRhs.hh"
45 #include "G4KM_DummyField.hh" // needed by G4M 42 #include "G4KM_DummyField.hh" // needed by G4Mag_EqRhs constructor.
46 #include "G4V3DNucleus.hh" 43 #include "G4V3DNucleus.hh"
47 44
48 //class G4KM_NucleonEqRhs : public G4EquationO 45 //class G4KM_NucleonEqRhs : public G4EquationOfMotion // I'd like
49 class G4KM_NucleonEqRhs : public G4Mag_EqRhs 46 class G4KM_NucleonEqRhs : public G4Mag_EqRhs
50 { 47 {
51 48
52 public: 49 public:
53 G4KM_NucleonEqRhs(G4KM_DummyField *field, G4 50 G4KM_NucleonEqRhs(G4KM_DummyField *field, G4V3DNucleus * nucleus);
54 ~G4KM_NucleonEqRhs(); 51 ~G4KM_NucleonEqRhs();
55 52
56 virtual void EvaluateRhsGivenB(const G4doubl 53 virtual void EvaluateRhsGivenB(const G4double y[], const G4double B[3],
57 G4double dydx[]) const; 54 G4double dydx[]) const;
58 virtual void SetChargeMomentumMass(G4ChargeS << 55 virtual void SetChargeMomentumMass(G4double particleCharge,
59 G4double MomentumXc, 56 G4double MomentumXc,
60 G4double MassXc2); 57 G4double MassXc2);
61 void SetMass(G4double aMass); 58 void SetMass(G4double aMass);
62 59
63 private: 60 private:
64 61
65 // use G4VKM_NuclearDensity as it needs a GetD 62 // use G4VKM_NuclearDensity as it needs a GetDeriv() method.
66 G4V3DNucleus * theNucleus; 63 G4V3DNucleus * theNucleus;
67 G4double factor; 64 G4double factor;
68 G4int A; 65 G4int A;
69 G4double theMass; 66 G4double theMass;
70 }; 67 };
71 68
72 69
73 70
74 inline void G4KM_NucleonEqRhs::SetMass(G4doubl 71 inline void G4KM_NucleonEqRhs::SetMass(G4double aMass)
75 { 72 {
76 theMass = aMass; 73 theMass = aMass;
77 } 74 }
78 75
79 76
80 77
81 inline G4KM_NucleonEqRhs::~G4KM_NucleonEqRhs() 78 inline G4KM_NucleonEqRhs::~G4KM_NucleonEqRhs()
>> 79 { }
>> 80
>> 81
>> 82 // Here by design, but it is unnecessary for nuclear fields
>> 83 inline void
>> 84 G4KM_NucleonEqRhs::SetChargeMomentumMass(G4double ,G4double ,G4double )
82 { } 85 { }
83 86
84 #endif 87 #endif
85 88
86 89
87 90