| 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 // -------------------------------------------
27 // GEANT4 Class file
28 //
29 //
30 // File name: G4LightIonQMDNucleus.hh
31 //
32 // Author: Koi, Tatsumi (tkoi@slac.stanfo
33 //
34 // Creation date: 3 April 2007
35 // -------------------------------------------
36 // 230309 Skyrme-QMD parameters added by Y-H.
37 // 230309 Total energy evaluated by Lorentz co
38 #ifndef G4LightIonQMDNucleus_hh
39 #define G4LightIonQMDNucleus_hh
40
41 #include "G4QMDSystem.hh"
42 #include "G4LightIonQMDParameters.hh"
43
44 class G4LightIonQMDNucleus : public G4QMDSyste
45 {
46 public:
47 G4LightIonQMDNucleus();
48 //virtual ~G4LightIonQMDNucleus();
49
50 G4LorentzVector Get4Momentum();
51
52 // Number of Nucleons (Proton or Neutron
53 G4int GetMassNumber();
54
55 // Number of Protons
56 G4int GetAtomicNumber();
57
58 void CalEnergyAndAngularMomentumInCM();
59
60 // rest mass from G4NucleiPropertiesTabl
61 G4double GetNuclearMass();
62
63 void SetTotalPotential( G4double x ){ po
64 G4double GetExcitationEnergy(){ return e
65
66 G4int GetAngularMomentum(){ return jj; }
67
68 private:
69
70 G4double hbc;
71
72 std::vector < G4ThreeVector > rcm, pcm;
73 std::vector < G4double > es;
74 G4int jj;
75
76 G4double potentialEnergy;
77 G4double excitationEnergy;
78 //G4double bindingEnergy;
79
80 //G4double kineticEnergyPerNucleon;
81 //G4double bindingEnergyPerNucleon;
82 //G4double potentialEnergyPerNucleon;
83
84 G4double GetNuclPotential( G4int );
85
86 //G4double rclds;
87
88 G4double rho0;
89 G4double epsx , epscl;
90
91 G4double cpc;
92
93 //G4int icoul, irelcr;
94 G4int irelcr;
95 G4double gamm, c0, c3, cs, cl, wl;
96 //G4double c0w, c3w, clw, c0sw;
97 G4double c0w, clw, c0sw;
98
99 G4double c0g,c3g,csg,pag;
100
101 G4double cpw,cph;
102
103 // Skyrme-QMD
104 G4double gtau0;
105 G4double g0;
106 G4double g0iso;
107 G4double eta;
108 G4double kappas;
109
110 G4double pag_tau;
111 G4double cg0;
112 G4double cgtau0;
113 };
114
115 #endif
116