source: trunk/source/processes/electromagnetic/lowenergy/include/G4QAOLowEnergyLoss.hh@ 1201

Last change on this file since 1201 was 819, checked in by garnier, 17 years ago

import all except CVS
File size: 5.8 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27// ------------------------------------------------------------
28// GEANT 4 class header file
29//
30// History: New Implementation
31//
32// ---------- G4QAOLowEnergyLoss physics process -------
33// by Stephane Chauvie, 21 May 2000
34//
35// Modified:
36// 16/09/2000 S. Chauvie Oscillator for all materials
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 energy loss of low energy antiprotons
42// Further documentation available from http://www.ge.infn.it/geant4/lowE
43
44// ------------------------------------------------------------
45
46
47#ifndef G4QAOLowEnergyLoss_hh
48#define G4QAOLowEnergyLoss_hh 1
49
50#include "G4VLowEnergyModel.hh"
51#include "globals.hh"
52
53class G4QAOLowEnergyLoss : public G4VLowEnergyModel
54{
55public:
56
57 G4QAOLowEnergyLoss(const G4String& name);
58
59 ~G4QAOLowEnergyLoss();
60
61 G4double HighEnergyLimit(const G4ParticleDefinition* aParticle,
62 const G4Material* material) const;
63 // returns the higher limit for model validity
64
65 G4double LowEnergyLimit(const G4ParticleDefinition* aParticle,
66 const G4Material* material) const;
67 // returns the lower limit for model validity
68
69 G4double HighEnergyLimit(const G4ParticleDefinition* aParticle) const;
70 // returns the higher limit for model validity
71
72 G4double LowEnergyLimit(const G4ParticleDefinition* aParticle) const;
73 // returns the lower limit for model validity
74
75 G4bool IsInCharge(const G4DynamicParticle* particle,
76 const G4Material* material) const;
77 // returns true if the model is applicable at that energy for
78 // that particle for that material
79
80 G4bool IsInCharge(const G4ParticleDefinition* aParticle,
81 const G4Material* material) const;
82 // returns true if the model is applicable at that energy for
83 // that particle for that material
84
85 G4double TheValue(const G4DynamicParticle* particle,
86 const G4Material* material);
87 // returns the energy loss via the quantal harmonic oscillator model
88
89 G4double TheValue(const G4ParticleDefinition* aParticle,
90 const G4Material* material,
91 G4double kineticEnergy);
92 // returns the energy loss via the quantal harmonic oscillator model
93
94private:
95
96 G4double EnergyLoss(const G4Material* material,
97 G4double kineticEnergy,
98 G4double zParticle) const;
99 // returns the energy loss via the quantal harmonic oscillator model
100
101 // get number of shell, energy and oscillator strenghts for material
102 G4int GetNumberOfShell(const G4Material* material) const;
103
104 G4double GetShellEnergy(const G4Material* material,G4int nbOfTheShell) const;
105 G4double GetOscillatorEnergy(const G4Material* material,G4int nbOfTheShell) const;
106 G4double GetShellStrength(const G4Material* material,G4int nbOfTheShell) const;
107 G4double GetOccupationNumber(G4int Z, G4int ShellNb) const;
108
109 // calculate stopping number for L's term
110 G4double GetL0(G4double normEnergy) const;
111 // terms in Z^2
112 G4double GetL1(G4double normEnergy) const;
113 // terms in Z^3
114 G4double GetL2(G4double normEnergy) const;
115 // terms in Z^4
116
117 // Z of element at now avaliable for the model
118 static const G4int materialAvailable[6];
119
120 // number, energy and oscillator strenghts
121 // for an harmonic oscillator model of material
122 static const G4int nbofShellForMaterial[6];
123 static G4double alShellEnergy[3];
124 static G4double alShellStrength[3];
125 static G4double siShellEnergy[3];
126 static G4double siShellStrength[3];
127 static G4double cuShellEnergy[4];
128 static G4double cuShellStrength[4];
129 static G4double taShellEnergy[6];
130 static G4double taShellStrength[6];
131 static G4double auShellEnergy[6];
132 static G4double auShellStrength[6];
133 static G4double ptShellEnergy[6];
134 static G4double ptShellStrength[6];
135
136 G4int numberOfMaterials;
137
138 // variable for calculation of stopping number of L's term
139 static const G4double L0[67][2];
140 static const G4double L1[22][2];
141 static const G4double L2[14][2];
142 static const G4int nbOfElectronPerSubShell[1540];
143 static const G4int fNumberOfShells[101];
144
145 G4int sizeL0;
146 G4int sizeL1;
147 G4int sizeL2;
148
149};
150
151#endif
Note: See TracBrowser for help on using the repository browser.