source: trunk/source/materials/include/G4SimpleMaterialStoppingICRU73.hh @ 1041

Last change on this file since 1041 was 986, checked in by garnier, 15 years ago

fichiers manquants
File size: 6.6 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// $Id: G4SimpleMaterialStoppingICRU73.hh,v 1.5 2009/02/19 12:19:53 vnivanch Exp $
27// GEANT4 tag $Name: geant4-09-02-ref-02 $
28
29#ifndef G4SimpleMaterialStoppingICRU73_h
30#define G4SimpleMaterialStoppingICRU73_h 1
31
32//---------------------------------------------------------------------------
33//
34// ClassName:   G4SimpleMaterialStoppingICRU73
35//
36// Description: Data on stopping powers for light ions in compounds
37//
38// Author:      Ivantchenko 10.07.2008
39//
40// Modifications:
41//
42//----------------------------------------------------------------------------
43//
44// Class Description:
45//
46// Data on Stopping Powers from the ICRU73 report
47//
48
49//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
50
51
52#include "globals.hh"
53#include "G4LPhysicsFreeVector.hh"
54#include <vector>
55
56class G4SimpleMaterialStoppingICRU73
57{
58public:
59
60  G4SimpleMaterialStoppingICRU73(G4bool splineFlag = true);
61
62  ~G4SimpleMaterialStoppingICRU73();
63
64  G4double GetDEDX(G4int ionZ, G4int idxMaterial, G4double kinEnergy);
65
66  inline G4double GetDEDX(G4int ionZ, const G4String& NameMaterial, 
67                          G4double kinEnergy);
68
69  inline G4int GetMaterialIndex(const G4String& NameMaterial);
70
71  // Function returns an unique index (>=0) for each ion-material couple (the
72  // return value is -1 if the couple is not found):
73  inline G4int GetIonMaterialCoupleIndex(
74                        G4int atomicNumber,            // Atomic number of ion
75                        const G4String& materialName); // Material name
76
77  inline G4double GetDensity(G4int idxMaterial);
78
79  inline G4String GetMaterialName(G4int idxMaterial);
80
81  inline G4PhysicsVector* GetPhysicsVector(G4int ionZ, G4int idxMaterial);
82
83  inline G4PhysicsVector* GetPhysicsVector(G4int ionZ, 
84                                           const G4String& NameMaterial);
85
86  inline G4double GetLowerEnergyBoundary();
87
88  inline G4double GetUpperEnergyBoundary();
89
90private:
91
92  void Initialise();
93
94  void AddData(G4double* energy, G4double* stoppower, G4double factor);
95
96  // hide assignment operator
97  G4SimpleMaterialStoppingICRU73 & operator=(const G4SimpleMaterialStoppingICRU73 &right);
98  G4SimpleMaterialStoppingICRU73(const G4SimpleMaterialStoppingICRU73&);
99
100  G4bool   spline;
101  G4int    Z[16];
102  G4double A[16];
103  G4String MatName[25];
104  G4double Density[25];
105
106  // Lower and upper energy boundaries for dE/dx vectors:
107  G4double lowerEnergyBoundary;
108  G4double upperEnergyBoundary;
109
110  std::vector<G4LPhysicsFreeVector*>  dedx;
111};
112
113//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
114
115inline G4double
116G4SimpleMaterialStoppingICRU73::GetDEDX(G4int ionZ, const G4String& NameMaterial, 
117                                        G4double kinEnergy)
118{
119  return GetDEDX(ionZ, GetMaterialIndex(NameMaterial), kinEnergy);
120}
121
122//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
123
124inline G4int
125G4SimpleMaterialStoppingICRU73::GetMaterialIndex(const G4String& NameMaterial)
126{
127  G4int idx = -1;
128  for (G4int i=0; i<25; i++){
129    if(MatName[i] == NameMaterial) {
130      idx = i;
131      break;
132    }
133  }
134  return idx;
135}
136
137//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
138
139inline G4int
140G4SimpleMaterialStoppingICRU73::GetIonMaterialCoupleIndex(
141                       G4int atomicNumber,             
142                       const G4String& materialName) 
143{ 
144  G4int idx = -1;
145  if(atomicNumber >= 3 && atomicNumber <= 18) {
146
147    G4int materialIndex = GetMaterialIndex(materialName);
148    if(materialIndex >= 0) idx = materialIndex * 16 + atomicNumber - 3;
149  }
150
151  return idx;
152}
153
154//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
155
156inline G4double
157G4SimpleMaterialStoppingICRU73::GetDensity(G4int idxMaterial)
158{
159  G4double d = 0.0;
160  if( idxMaterial >= 0 && idxMaterial <= 24) d = Density[idxMaterial];
161  return d;
162}
163
164//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
165
166inline 
167G4String G4SimpleMaterialStoppingICRU73::GetMaterialName(G4int idxMaterial)
168{
169  G4String s = "";
170  if( idxMaterial >= 0 && idxMaterial <= 24) s = MatName[idxMaterial];
171  return s;
172}
173
174//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
175
176inline G4PhysicsVector* 
177G4SimpleMaterialStoppingICRU73::GetPhysicsVector(G4int ionZ, G4int idxMaterial)
178{
179  G4PhysicsVector* v = 0;
180  if(ionZ >= 3 && ionZ <= 18 && idxMaterial >= 0 && idxMaterial <= 24) { 
181    v = dedx[idxMaterial*16 + ionZ - 3];
182  }
183  return v;
184}
185
186//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
187
188inline G4PhysicsVector* 
189G4SimpleMaterialStoppingICRU73::GetPhysicsVector(G4int ionZ, 
190                                                 const G4String& NameMaterial)
191{
192  return GetPhysicsVector(ionZ, GetMaterialIndex(NameMaterial));
193}
194
195//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
196
197inline G4double
198G4SimpleMaterialStoppingICRU73::GetLowerEnergyBoundary() {
199
200  return lowerEnergyBoundary;
201}
202
203//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
204
205inline G4double
206G4SimpleMaterialStoppingICRU73::GetUpperEnergyBoundary() {
207
208  return upperEnergyBoundary;
209}
210
211
212#endif
Note: See TracBrowser for help on using the repository browser.