source: trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPNames.cc@ 830

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

import all except CVS
File size: 11.5 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// neutron_hp -- source file
27// J.P. Wellisch, Nov-1996
28// A prototype of the low energy neutron transport model.
29//
30// 24-Jan-07 Enable to use exact data only and add warnig when substitute file is used T. Koi
31// 30-Jan-07 Modified method of searching substitute isotope data by T. Koi
32// 07-06-12 fix memory leaking by T. Koi
33// 07-06-25 Change data selection logic when G4NEUTRONHP_SKIP_MISSING_ISOTOPES is turn on
34// Natural Abundance data are allowed. by T. Koi
35// 07-07-06 Allow _nat_ final state even for isotoped cross sections by T. Koi
36//
37#include "G4NeutronHPNames.hh"
38#include "G4SandiaTable.hh"
39#include "G4HadronicException.hh"
40#include <fstream>
41
42 const G4String G4NeutronHPNames::theString[99] = {"Hydrogen", "Helium",
43 "Lithium", "Berylium", "Boron", "Carbon", "Nitrogen", "Oxygen", "Fluorine",
44 "Neon", "Sodium", "Magnesium", "Aluminum", "Silicon", "Phosphorous",
45 "Sulfur", "Chlorine", "Argon", "Potassium", "Calcium", "Scandium",
46 "Titanium", "Vanadium", "Chromium", "Manganese", "Iron", "Cobalt", "Nickel",
47 "Copper", "Zinc", "Gallium", "Germanium", "Arsenic", "Selenium", "Bromine",
48 "Krypton", "Rubidium", "Strontium", "Yttrium", "Zirconium", "Niobium",
49 "Molybdenum", "Technetium", "Ruthenium", "Rhodium", "Palladium", "Silver",
50 "Cadmium", "Indium", "Tin", "Antimony", "Tellurium", "Iodine", "Xenon",
51 "Cesium", "Barium", "Lanthanum", "Cerium", "Praseodymium", "Neodymium",
52 "Promethium", "Samarium", "Europium", "Gadolinium", "Terbium", "Dysprosium",
53 "Holmium", "Erbium", "Thulium", "Ytterbium", "Lutetium", "Hafnium",
54 "Tantalum", "Tungsten", "Rhenium", "Osmium", "Iridium", "Platinium", "Gold",
55 "Mercury", "Thallium", "Lead", "Bismuth", "Polonium", "Astatine", "Radon",
56 "Francium", "Radium", "Actinium ", "Thorium", "Protactinium", "Uranium",
57 "Neptunium", "Plutonium", "Americium", "Curium", "Berkelium", "Californium",
58 "Einsteinium"};
59
60
61 G4String G4NeutronHPNames::GetName(G4int i) { return theString[i]; }
62
63G4NeutronHPDataUsed G4NeutronHPNames::GetName(G4int A, G4int Z, G4String base, G4String rest, G4bool & aFlag)
64{
65
66 G4NeutronHPDataUsed result;
67 aFlag = true;
68if(getenv("NeutronHPNames")) G4cout << "Names::GetName entered for Z = " << Z << ", A = " << A <<G4endl;
69
70 G4int myA = A;
71 G4int myZ = Z;
72
73 if(Z>92.5&&!getenv("AllowForHeavyElements") )
74 {
75 G4cerr << "Please contact Hans-Peter.Wellisch@cern.ch"<<G4endl;
76 throw G4HadronicException(__FILE__, __LINE__, "G4NeutronHPNames::GetName - data with Z>92 are not provided");
77 }
78
79 G4String * theName = 0;
80 G4String theFileName("");
81
82// G4int inc = 1;
83
84 G4int flip_Z = 1;
85 G4int delta_Z = 0;
86
87 G4int flip_A = 1;
88 G4int delta_A = 0;
89
90 std::ifstream * check = new std::ifstream(".dummy");
91 G4bool first = true;
92if(getenv("NeutronHPNames")) G4cout << "entered GetName!!!"<<G4endl;
93 do
94 {
95 aFlag = true;
96 G4String * biff = new G4String(); // delete here as theName
97 *biff = base+"/"+"CrossSection/"+itoa(myZ)+"_"+itoa(myA)+"_"+theString[myZ-1];
98
99 if(theName!=0) delete theName;
100 theName = biff;
101 result.SetName(*theName);
102 result.SetA(myA);
103 result.SetZ(myZ);
104if(getenv("NeutronHPNames")) G4cout <<"HPWD 1 "<<*theName<<G4endl;
105
106 // T.K. debug for memory leak
107 if ( check != 0 )
108 {
109 check->close();
110 delete check;
111 }
112 check = new std::ifstream(*theName);
113 if ( !(*check) )
114 {
115 check->close();
116 delete check;
117 check = 0;
118 aFlag = false;
119 if ( first )
120 {
121 aFlag = true;
122 first = false;
123 biff = new G4String(); // delete here as theName
124 *biff = base+"/"+"CrossSection/"+itoa(myZ)+"_"+"nat"+"_"+theString[myZ-1];
125 if(theName!=0) delete theName;
126 theName = biff;
127if(getenv("NeutronHPNames")) G4cout <<"HPWD 2 "<<*theName<<G4endl;
128 result.SetName(*theName);
129 G4double natA = myZ/G4SandiaTable::GetZtoA(myZ);
130 result.SetA(natA);
131 result.SetZ(myZ);
132 check = new std::ifstream(*theName);
133 if ( !(*check) )
134 {
135 check->close();
136 delete check;
137 check = 0;
138 aFlag = false;
139 }
140 else
141 {
142 biff = new G4String(); // delete here as theName
143 if(theName!=0) delete theName;
144 *biff = base+"/"+rest+itoa(myZ)+"_"+"nat"+"_"+theString[myZ-1];
145 theName = biff;
146if(getenv("NeutronHPNames")) G4cout <<"HPWD 3 "<<*theName<<G4endl;
147 result.SetName(*theName);
148 G4double natA = myZ/G4SandiaTable::GetZtoA(myZ);
149 result.SetA(natA);
150 result.SetZ(myZ);
151 result.SetNaturalAbundanceFlag();
152 }
153 }
154 }
155 else
156 {
157// 070706 T. Koi Modified
158/*
159 biff = new G4String(); // delete here as theName
160 *biff = base+"/"+rest+itoa(myZ)+"_"+itoa(myA)+"_"+theString[myZ-1];
161 if(theName!=0) delete theName;
162 theName = biff;
163if(getenv("NeutronHPNames")) G4cout <<"HPWD 4 "<<*theName<<G4endl;
164 result.SetName(*theName);
165 result.SetA(myA);
166 result.SetZ(myZ);
167*/
168
169 G4double tmpA = myA;
170 std::ifstream* file = NULL;
171 G4String fileName;
172
173 if ( rest == "/CrossSection/" )
174 {
175
176 fileName = base+"/"+rest+itoa(myZ)+"_"+itoa(myA)+"_"+theString[myZ-1];
177if(getenv("NeutronHPNames")) G4cout <<"HPWD 4a "<<*theName<<G4endl;
178
179 }
180 else
181 {
182
183// For FS
184 fileName = base+"/"+rest+itoa(myZ)+"_"+itoa(myA)+"_"+theString[myZ-1];
185 file = new std::ifstream(fileName);
186
187 if ( *file )
188 {
189
190// isotope FS
191if(getenv("NeutronHPNames")) G4cout <<"HPWD 4b1 "<<*theName<<G4endl;
192 }
193 else
194 {
195
196// _nat_ FS
197 fileName = base+"/"+rest+itoa(myZ)+"_"+"nat"+"_"+theString[myZ-1];
198
199 delete file;
200 file = new std::ifstream(fileName);
201 if ( *file )
202 {
203
204// FS neither isotope nor _nat_
205if(getenv("NeutronHPNames")) G4cout <<"HPWD 4b2a "<<*theName<<G4endl;
206 G4double natA = myZ/G4SandiaTable::GetZtoA(myZ);
207 tmpA = natA;
208 }
209 else
210 {
211if(getenv("NeutronHPNames")) G4cout <<"HPWD 4b2c "<<*theName<<G4endl;
212 }
213 }
214
215 delete file;
216
217 }
218
219 result.SetName(fileName);
220 result.SetA(tmpA);
221 result.SetZ(myZ);
222
223 }
224
225 do
226 {
227// if (std::abs(myZ-Z)>theMaxOffSet||myZ==0||myA==0)
228 if ( delta_Z > theMaxOffSet )
229 {
230 //if ( inc > 0 )
231 //{
232 // inc*= -1;
233 // myZ = Z;
234 // myA = A;
235 //}
236 //else
237 //{
238 G4cout <<"G4NeutronHPNames: Sorry, this material does not come near to any data."<<G4endl;
239 G4cout <<"G4NeutronHPNames: Please make sure G4NEUTRONHPDATA points to the" << G4endl;
240 G4cout <<" directory, the neutron scattering data are located in." << G4endl;
241 G4cout << "G4NeutronHPNames: The material was: A="<<A<<", Z="<<Z<<G4endl;
242 throw G4HadronicException(__FILE__, __LINE__, "In case the data sets are at present not available in the neutron data library, please contact Hans-Peter.Wellisch@cern.ch");
243 delete theName;
244 theFileName = "";
245 return result;
246 //}
247 }
248
249 //if ( std::abs( myA - A ) > theMaxOffSet )
250 if ( delta_A > 2*theMaxOffSet )
251 {
252 delta_A = 0;
253 flip_A = 1;
254
255 first = true;
256
257 if ( flip_Z > 0 )
258 {
259 delta_Z +=1;
260 }
261 myZ = Z + flip_Z * delta_Z;
262 flip_Z *= -1;
263
264 myA = A;
265 if ( myZ > 99 )
266 {
267 myZ = 99;
268 }
269 if ( myZ < 1 )
270 {
271 myZ = 1;
272 }
273
274// myZ += inc;
275 }
276 else
277 {
278 if ( flip_A > 0 )
279 {
280 delta_A += 1;
281 }
282 myA = A + flip_A * delta_A;
283 flip_A *= -1;
284
285 if ( myA < 1 )
286 {
287 myA = 1;
288 }
289
290// myA += inc;
291 }
292
293 }
294 while( myZ == 0 || myA == 0 ); // No meaning
295
296 }
297 while((!check) || (!(*check)));
298
299 if(getenv("NeutronHPNamesLogging") || getenv("NeutronHPNames"))
300 {
301 G4cout << "Names::GetName: last theName proposal = "<< G4endl;
302 G4cout << *theName <<" "<<A<<" "<<Z<<" "<<result.GetName()<<G4endl;
303 }
304
305// administration and anouncement for lacking of exact data in NDL
306 if ( Z != result.GetZ() || A != result.GetA() )
307 {
308 if ( rest == "/CrossSection/" )
309 {
310 G4String reac = base;
311 G4String dir = getenv("G4NEUTRONHPDATA");
312 reac.erase ( 0 , dir.length() );
313 if ( getenv ( "G4NEUTRONHP_SKIP_MISSING_ISOTOPES" ) && !( Z == result.GetZ() && result.IsThisNaturalAbundance() ) )
314 {
315 G4cout << "NeutronHP: " << reac << " file for Z = " << Z << ", A = " << A << " is not found and CrossSection set to 0." << G4endl;
316 G4String new_name = base+"/"+rest+"0_0_Zero";
317 result.SetName( new_name );
318 }
319 else
320 {
321 G4cout << "NeutronHP: " << reac << " file for Z = " << Z << ", A = " << A << " is not found and NeutronHP will use " << result.GetName() << G4endl;
322 }
323 }
324 }
325
326 delete theName;
327 if(aFlag)
328 {
329 check->close();
330 delete check;
331 check = 0;
332 }
333 return result;
334 }
Note: See TracBrowser for help on using the repository browser.