| 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 // Author: Maria Grazia Pia (Maria.Grazia.Pia@
29 //
30 // History:
31 // -----------
32 // 1 Aug 2001 MGP Created
33 // 31 Jul 2008 MGP Revised and rename
34 //
35 //
36 // -------------------------------------------
37
38 #include "G4CompositeDataSet.hh"
39 #include "G4DataSet.hh"
40 #include "G4IInterpolator.hh"
41 #include <fstream>
42 #include <sstream>
43
44 G4CompositeDataSet::G4CompositeDataSet(G4IInte
45 G4double eUnit,
46 G4double dataUnit,
47 G4int zMin,
48 G4int zMax)
49 :
50 algorithm(algo),
51 unitEnergies(eUnit),
52 unitData(dataUnit),
53 minZ(zMin),
54 maxZ(zMax)
55 {
56 if (algorithm == 0)
57 G4Exception("G4CompositeDataSet::G4Composi
58 "pii00000001",
59 FatalException,
60 "Interpolation == 0");
61 }
62
63
64
65 G4CompositeDataSet::~G4CompositeDataSet()
66 {
67 CleanUpComponents();
68 if (algorithm) delete algorithm;
69 }
70
71
72 G4double G4CompositeDataSet::FindValue(G4doubl
73 {
74 const G4IDataSet* component(GetComponent(com
75
76 if (component) return component->FindValue(e
77
78 std::ostringstream message;
79 message << "G4CompositeDataSet::FindValue -
80
81 G4Exception("G4CompositeDataSet::FindValue"
82 "pii00000010",
83 FatalException,
84 message.str().c_str());
85
86 return 0.;
87 }
88
89 void G4CompositeDataSet::PrintData(void) const
90 {
91 const G4int n = (G4int)NumberOfComponents();
92
93 G4cout << "The data set has " << n << " comp
94 G4cout << G4endl;
95
96 G4int i(0);
97
98 while (i<n)
99 {
100 G4cout << "--- Component " << i << " ---
101 GetComponent(i)->PrintData();
102 ++i;
103 }
104 }
105
106 void G4CompositeDataSet::SetEnergiesData(G4Dat
107 {
108 G4IDataSet * component(components[componentI
109
110 if (component)
111 {
112 component->SetEnergiesData(energies, dat
113 return;
114 }
115
116 std::ostringstream message;
117 message << "G4CompositeDataSet::SetEnergiesD
118
119 G4Exception("G4CompositeDataSet::SetEnergies
120 "pii00000020",
121 FatalException,
122 message.str().c_str());
123
124 }
125
126 G4bool G4CompositeDataSet::LoadData(const G4St
127 {
128 CleanUpComponents();
129
130 for (G4int z(minZ); z<maxZ; ++z)
131 {
132 G4IDataSet* component = new G4DataSet(z,
133 if (!component->LoadData(argFileName))
134 {
135 delete component;
136 return false;
137 }
138 AddComponent(component);
139 }
140 return true;
141 }
142
143
144
145 G4bool G4CompositeDataSet::SaveData(const G4St
146 {
147 for (G4int z=minZ; z<maxZ; ++z)
148 {
149 const G4IDataSet* component(GetComponent
150
151 if (!component)
152 {
153 std::ostringstream message;
154 message << "G4CompositeDataSet::SaveData -
155 G4Exception("G4CompositeDataSet::SaveData"
156 "pii00000030",
157 FatalException,
158 message.str().c_str());
159 }
160
161 if (!component->SaveData(argFileName))
162 return false;
163 }
164
165 return true;
166 }
167
168 void G4CompositeDataSet::CleanUpComponents(voi
169 {
170 while (!components.empty())
171 {
172 if (components.back())
173 delete components.back();
174 components.pop_back();
175 }
176 }
177
178
179 G4double G4CompositeDataSet::RandomSelect(G4in
180 {
181 G4double value = 0.;
182 if (componentId >= 0 && componentId < (G4int
183 {
184 const G4IDataSet* dataSet = GetComponent
185 value = dataSet->RandomSelect();
186 }
187 return value;
188 }
189