| 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 // -------------------------------------------
28 // GEANT 4 class header file
29 //
30 // History: New Implementation
31 //
32 // ---------- G4QAOLowEnergyLoss physics
33 // by Stephane Chauvie, 21 Ma
34 //
35 // Modified:
36 // 16/09/2000 S. Chauvie Oscillator for all m
37 // 23/05/2000 MGP Made compliant to design
38 // 01/06/2001 V.Ivanchenko replace names by Z
39 //
40 // Class description:
41 // Quantal Harmonic Oscillator Model for energ
42 // Further documentation available from http:/
43
44 // -------------------------------------------
45
46
47 #ifndef G4QAOLowEnergyLoss_hh
48 #define G4QAOLowEnergyLoss_hh 1
49
50 #include "G4VLowEnergyModel.hh"
51 #include "globals.hh"
52
53 class G4QAOLowEnergyLoss : public G4VLowEnergy
54 {
55 public:
56 explicit G4QAOLowEnergyLoss(const G4String&
57 ~G4QAOLowEnergyLoss();
58
59 G4double HighEnergyLimit(const G4ParticleDef
60 const G4Material* m
61 // returns the higher limit for model validi
62
63 G4double LowEnergyLimit(const G4ParticleDefi
64 const G4Material* ma
65 // returns the lower limit for model validit
66
67 G4double HighEnergyLimit(const G4ParticleDef
68 // returns the higher limit for model validi
69
70 G4double LowEnergyLimit(const G4ParticleDefi
71 // returns the lower limit for model validit
72
73 G4bool IsInCharge(const G4DynamicParticle* p
74 const G4Material* material) const over
75 // returns true if the model is applicable a
76 // that particle for that material
77
78 G4bool IsInCharge(const G4ParticleDefinition
79 const G4Material* material) const over
80 // returns true if the model is applicable a
81 // that particle for that material
82
83 G4double TheValue(const G4DynamicParticle* p
84 const G4Material* material) ove
85 // returns the energy loss via the quantal h
86
87 G4double TheValue(const G4ParticleDefinition
88 const G4Material* material
89 G4double kinet
90 // returns the energy loss via the quantal h
91
92 private:
93 G4double EnergyLoss(const G4Material* materi
94 G4double kineticEn
95 G4double zParticle
96 // returns the energy loss via the quantal h
97
98 // get number of shell, energy and oscillato
99 G4int GetNumberOfShell(const G4Material* mat
100
101 G4double GetShellEnergy(const G4Material* ma
102 G4double GetOscillatorEnergy(const G4Materia
103 G4double GetShellStrength(const G4Material*
104 G4double GetOccupationNumber(G4int Z, G4int
105
106 // calculate stopping number for L's term
107 G4double GetL0(G4double normEnergy) const;
108 // terms in Z^2
109 G4double GetL1(G4double normEnergy) const;
110 // terms in Z^3
111 G4double GetL2(G4double normEnergy) const;
112 // terms in Z^4
113
114 // number, energy and oscillator strengths
115 // for an harmonic oscillator model of mater
116 static const G4int nbofShellForMaterial[6];
117 static const G4double alShellEnergy[3];
118 static const G4double alShellStrength[3];
119 static const G4double siShellEnergy[3];
120 static const G4double siShellStrength[3];
121 static const G4double cuShellEnergy[4];
122 static const G4double cuShellStrength[4];
123 static const G4double taShellEnergy[6];
124 static const G4double taShellStrength[6];
125 static const G4double auShellEnergy[6];
126 static const G4double auShellStrength[6];
127 static const G4double ptShellEnergy[6];
128 static const G4double ptShellStrength[6];
129 // variable for calculation of stopping num
130 static const G4double L0[67][2];
131 static const G4double L1[22][2];
132 static const G4double L2[14][2];
133 static const G4int nbOfElectronPerSubShell[1
134 static const G4int fNumberOfShells[101];
135
136 // Z of element at now avaliable for the mod
137 static const G4int materialAvailable[6];
138
139 G4int numberOfMaterials;
140 G4int sizeL0;
141 G4int sizeL1;
142 G4int sizeL2;
143
144 };
145
146 #endif
147