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 | // $Id: G4DNACrossSectionDataSetTest.cc,v 1.2 2006/06/29 19:43:55 gunter Exp $ |
---|
28 | // GEANT4 tag $Name: geant4-09-04-ref-00 $ |
---|
29 | // |
---|
30 | // Author: Riccardo Capra <capra@ge.infn.it> |
---|
31 | // |
---|
32 | // History: |
---|
33 | // ----------- |
---|
34 | // 30 Jun 2005 RC Created |
---|
35 | |
---|
36 | |
---|
37 | #include "G4DNACrossSectionDataSet.hh" |
---|
38 | #include "G4LogLogInterpolation.hh" |
---|
39 | #include <fstream> |
---|
40 | #include <unistd.h> |
---|
41 | |
---|
42 | #define EnergyUMIS 2 //eV |
---|
43 | #define DataUMIS 3 //barn |
---|
44 | |
---|
45 | int main(int argc, char **argv) |
---|
46 | { |
---|
47 | const char *inFile; |
---|
48 | |
---|
49 | if (argc>2) |
---|
50 | { |
---|
51 | G4cout << "Syntax: " << argv[0] << " [filename]" << G4endl; |
---|
52 | G4cout << G4endl; |
---|
53 | G4cout << "(!) NOTE: some files will be written in this directory (pattern: G4DNACrossSectionDataSetTestFile*.dat)" << G4endl; |
---|
54 | return 0; |
---|
55 | } |
---|
56 | else if (argc==2) |
---|
57 | inFile=argv[1]; |
---|
58 | else |
---|
59 | { |
---|
60 | inFile="G4DNACrossSectionDataSetTestFile"; |
---|
61 | std::ofstream out("./G4DNACrossSectionDataSetTestFile.dat"); |
---|
62 | out << "1.00000000e+05 8.14416702e+02 6.76315789e+02 4.88340749e+02 8.39880339e+01 7.89635115e-02 2.06314024e+03\n" |
---|
63 | "2.00000000e+05 5.28359350e+02 4.63221618e+02 3.44771389e+02 7.30146866e+01 2.52677590e-01 1.40961972e+03\n\n" |
---|
64 | " 3.00000000e+05 3.95840046e+02 3.53745700e+02 2.65748336e+02\t5.98188864e+01 4.09845238e-01 1.07556281e+03\n" |
---|
65 | "4.00000000e+05 3.19343880e+02 2.88487842e+02 2.17686875e+02 5.03590031e+01 5.32869987e-01 8.76410470e+02\r\n" |
---|
66 | "5.00000000e+05 2.69156735e+02 2.44940828e+02 1.85347963e+02 4.34715722e+01 6.24418579e-01 7.43541517e+02\n" |
---|
67 | "1.00000000e+06 1.55452794e+02 1.43970277e+02 1.09881243e+02 2.64352725e+01 8.10553505e-01 4.36550140e+02 \n" |
---|
68 | "1.50000000e+06 1.11668214e+02 1.04282494e+02 7.98560948e+01\t1.93490562e+01 8.14261386e-01 3.15970120e+02\n" |
---|
69 | "2.00000000e+06 8.80312507e+01 8.26279732e+01 6.33649352e+01 1.54081437e+01 7.71586134e-01 2.50203889e+02\n" |
---|
70 | " 2.50000000e+06 7.30696756e+01 6.88018974e+01 5.28437021e+01 1.28756775e+01 7.19674081e-01 2.08310627e+02\n\r" |
---|
71 | "3.00000000e+06 6.26612557e+01 5.91740007e+01 4.55044481e+01 1.11003552e+01 6.69053959e-01 1.79109114e+02\n" |
---|
72 | " 3.50000000e+06 5.50134904e+01 5.20584986e+01 4.00630585e+01 9.78183471e+00 6.23068350e-01 1.57539951e+02\n" |
---|
73 | "4.00000000e+06 4.91088724e+01 4.65679165e+01 3.58555442e+01 8.76079088e+00 5.82043180e-01 1.40875167e+02\n\n # comment\n" |
---|
74 | "4.50000000e+06 4.44233138e+01 4.21789727e+01 3.25041581e+01 7.94478293e+00 5.45740594e-01 1.27596968e+02\n#comment\r\n" |
---|
75 | "5.00000000e+06 4.05969944e+01 3.85988273e+01 2.97569104e+01 7.27564958e+00 5.13600786e-01 1.16741982e+02\n" |
---|
76 | "6.00000000e+06 3.47166562e+01 3.30808108e+01 2.55256187e+01 6.24407052e+00 4.59695847e-01 1.00026852e+02\n" |
---|
77 | "7.00000000e+06 3.03909338e+01 2.90156224e+01 2.24051103e+01 5.48359988e+00 4.16287624e-01 8.77115540e+01\n" |
---|
78 | "8.00000000e+06 2.70812879e+01 2.58941898e+01 2.00065525e+01\t4.89722388e+00 3.80739730e-01 7.82599938e+01\n" |
---|
79 | "9.00000000e+06 2.44599065e+01 2.34131260e+01 1.80955062e+01 4.43069001e+00 3.51212438e-01 7.07504411e+01\n" |
---|
80 | "1.00000000e+07 2.23229847e+01 2.13902553e+01 1.65394526e+01 4.05000347e+00 3.26180582e-01 6.46288766e+01 #comment \n" |
---|
81 | "1.10000000e+07 2.05485697e+01 1.97086000e+01 1.52450755e+01 3.73311927e+00 3.04729195e-01 5.95400936e+01#comment\n" |
---|
82 | "1.20000000e+07 1.90475083e+01 1.82869443e+01 1.41492895e+01 3.46491210e+00 2.86140746e-01 5.52347951e+01\n" |
---|
83 | "1.30000000e+07 1.77638183e+01 1.70693653e+01 1.32099385e+01 3.23474125e+00 \t 2.69857202e-01 5.15477205e+01 \n" |
---|
84 | " 1.40000000e+07 1.66514478e+01 1.60114678e+01 1.23958510e+01 3.03496833e+00 2.55467589e-01 4.83492025e+01\n" |
---|
85 | "1.50000000e+07 1.56760259e+01 1.50837486e+01 1.16792471e+01 2.85979032e+00 2.42653835e-01 4.55414657e+01\n" |
---|
86 | "1.60000000e+07 1.48138030e+01 1.42630410e+01 1.10462104e+01 2.70486989e+00 2.31166719e-01 4.30590910e+01 \n\n"; |
---|
87 | |
---|
88 | if (!out.good()) |
---|
89 | { |
---|
90 | G4cout << argv[0] << ": Failed to create \"./G4DNACrossSectionDataSetTestFile.dat\" file" << G4endl; |
---|
91 | return 0; |
---|
92 | } |
---|
93 | |
---|
94 | out.close(); |
---|
95 | setenv("G4LEDATA", ".", 1); |
---|
96 | } |
---|
97 | |
---|
98 | G4LogLogInterpolation *algorithm(new G4LogLogInterpolation); |
---|
99 | |
---|
100 | G4DNACrossSectionDataSet dataset(algorithm, EnergyUMIS, DataUMIS); |
---|
101 | |
---|
102 | G4cout << argv[0] << ": G4DNACrossSectionDataSet::Load: " << G4endl; |
---|
103 | |
---|
104 | try |
---|
105 | { |
---|
106 | if (!dataset.LoadData(inFile)) |
---|
107 | { |
---|
108 | G4cout << " - Failed with false" << G4endl; |
---|
109 | return 0; |
---|
110 | } |
---|
111 | } |
---|
112 | catch (...) |
---|
113 | { |
---|
114 | G4cout << " - Failed with exception" << G4endl; |
---|
115 | return 0; |
---|
116 | } |
---|
117 | |
---|
118 | if (argc==1) |
---|
119 | unlink("./G4DNACrossSectionDataSetTestFile.dat"); |
---|
120 | |
---|
121 | G4cout << " - Passed" << G4endl; |
---|
122 | |
---|
123 | G4cout << argv[0] << ": G4DNACrossSectionDataSet::PrintData: " << G4endl; |
---|
124 | dataset.PrintData(); |
---|
125 | |
---|
126 | setenv("G4LEDATA", ".", 1); |
---|
127 | G4cout << argv[0] << ": G4DNACrossSectionDataSet::Save: " << G4endl; |
---|
128 | try |
---|
129 | { |
---|
130 | if (!dataset.SaveData("G4DNACrossSectionDataSetTestFile-2")) |
---|
131 | { |
---|
132 | G4cout << " - Failed with false" << G4endl; |
---|
133 | return 0; |
---|
134 | } |
---|
135 | } |
---|
136 | catch (...) |
---|
137 | { |
---|
138 | G4cout << " - Failed with exception" << G4endl; |
---|
139 | return 0; |
---|
140 | } |
---|
141 | |
---|
142 | G4DNACrossSectionDataSet dataset2(algorithm->Clone(), EnergyUMIS, DataUMIS); |
---|
143 | try |
---|
144 | { |
---|
145 | if (!dataset2.LoadData("G4DNACrossSectionDataSetTestFile-2")) |
---|
146 | { |
---|
147 | G4cout << " - Failed with file missing" << G4endl; |
---|
148 | return 0; |
---|
149 | } |
---|
150 | } |
---|
151 | catch (...) |
---|
152 | { |
---|
153 | G4cout << " - Failed with file missing" << G4endl; |
---|
154 | return 0; |
---|
155 | } |
---|
156 | |
---|
157 | unlink("./G4DNACrossSectionDataSetTestFile-2.dat"); |
---|
158 | |
---|
159 | G4cout << " - Passed" << G4endl; |
---|
160 | |
---|
161 | G4cout << argv[0] << ": Comparing data:" << G4endl; |
---|
162 | |
---|
163 | size_t n(dataset.NumberOfComponents()); |
---|
164 | |
---|
165 | if (n!=dataset2.NumberOfComponents()) |
---|
166 | { |
---|
167 | G4cout << " - Failed different number of components" << G4endl; |
---|
168 | return 0; |
---|
169 | } |
---|
170 | |
---|
171 | while (n>0) |
---|
172 | { |
---|
173 | n--; |
---|
174 | |
---|
175 | const G4DataVector &e(dataset.GetComponent(n)->GetEnergies(0)); |
---|
176 | const G4DataVector &e2(dataset2.GetComponent(n)->GetEnergies(0)); |
---|
177 | const G4DataVector &d(dataset.GetComponent(n)->GetData(0)); |
---|
178 | const G4DataVector &d2(dataset2.GetComponent(n)->GetData(0)); |
---|
179 | |
---|
180 | if (e.size()!=e2.size() || d.size()!=d2.size() || e.size()!=d.size()) |
---|
181 | { |
---|
182 | G4cout << " - Failed different number of rows" << G4endl; |
---|
183 | return 0; |
---|
184 | } |
---|
185 | |
---|
186 | size_t m(e.size()); |
---|
187 | |
---|
188 | while (m>0) |
---|
189 | { |
---|
190 | m--; |
---|
191 | |
---|
192 | if (e[m]!=e2[m]) |
---|
193 | { |
---|
194 | G4cout << " - Failed different energies (line " << (m+1) << " - " << e[m] << ", " << e2[m] << ")" << G4endl; |
---|
195 | return 0; |
---|
196 | } |
---|
197 | |
---|
198 | if (d[m]!=d2[m]) |
---|
199 | { |
---|
200 | G4cout << " - Failed different data (line " << (m+1) << ", column " << (n+2) << " - " << d[m] << ", " << d2[m] << ")" << G4endl; |
---|
201 | return 0; |
---|
202 | } |
---|
203 | } |
---|
204 | } |
---|
205 | |
---|
206 | G4cout << " - Passed" << G4endl; |
---|
207 | G4cout << G4endl; |
---|
208 | G4cout << "(!) NOTE: Files G4DNACrossSectionDataSetTestFile*.dat were written and removed from this directory" << G4endl; |
---|
209 | |
---|
210 | return 0; |
---|
211 | } |
---|