OPALX (Object Oriented Parallel Accelerator Library for Exascal) master (dc2a29eed580)
OPALX
Loading...
Searching...
No Matches
TestLinearBreitWheelerSpectrum.cpp
Go to the documentation of this file.
1
17#include "LinearBreitWheelerBenchmarkCommon.h"
18#include "Utilities/Options.h"
19#include "gtest/gtest.h"
20
21#include <filesystem>
22
23namespace {
24 std::filesystem::path referenceElectronEnergyPath() {
25 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
26 / "cain_linear_breit_wheeler_head_on_electron_energy.csv";
27 }
28
29 std::filesystem::path referencePositronEnergyPath() {
30 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
31 / "cain_linear_breit_wheeler_head_on_positron_energy.csv";
32 }
33
34 std::filesystem::path referenceElectronThetaPath() {
35 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
36 / "cain_linear_breit_wheeler_head_on_electron_theta.csv";
37 }
38
39 std::filesystem::path referencePositronThetaPath() {
40 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
41 / "cain_linear_breit_wheeler_head_on_positron_theta.csv";
42 }
43
44 std::filesystem::path referenceElectronJointPath() {
45 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
46 / "cain_linear_breit_wheeler_head_on_electron_joint.csv";
47 }
48
49 std::filesystem::path referencePositronJointPath() {
50 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
51 / "cain_linear_breit_wheeler_head_on_positron_joint.csv";
52 }
53
54 std::filesystem::path referenceFinitePhotonBeamElectronEnergyPath() {
55 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
56 / "cain_linear_breit_wheeler_finite_photon_beam_electron_energy.csv";
57 }
58
59 std::filesystem::path referenceFinitePhotonBeamPositronEnergyPath() {
60 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
61 / "cain_linear_breit_wheeler_finite_photon_beam_positron_energy.csv";
62 }
63
64 std::filesystem::path referenceFinitePhotonBeamElectronThetaPath() {
65 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
66 / "cain_linear_breit_wheeler_finite_photon_beam_electron_theta.csv";
67 }
68
69 std::filesystem::path referenceFinitePhotonBeamPositronThetaPath() {
70 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
71 / "cain_linear_breit_wheeler_finite_photon_beam_positron_theta.csv";
72 }
73
74 std::filesystem::path referenceFinitePhotonBeamElectronJointPath() {
75 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
76 / "cain_linear_breit_wheeler_finite_photon_beam_electron_joint.csv";
77 }
78
79 std::filesystem::path referenceFinitePhotonBeamPositronJointPath() {
80 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
81 / "cain_linear_breit_wheeler_finite_photon_beam_positron_joint.csv";
82 }
83
84 std::filesystem::path referenceFinitePhotonBeamEnergySpreadElectronEnergyPath() {
85 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
86 / "cain_linear_breit_wheeler_finite_photon_beam_energy_spread_electron_energy.csv";
87 }
88
89 std::filesystem::path referenceFinitePhotonBeamEnergySpreadPositronEnergyPath() {
90 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
91 / "cain_linear_breit_wheeler_finite_photon_beam_energy_spread_positron_energy.csv";
92 }
93
94 std::filesystem::path referenceFinitePhotonBeamEnergySpreadElectronThetaPath() {
95 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
96 / "cain_linear_breit_wheeler_finite_photon_beam_energy_spread_electron_theta.csv";
97 }
98
99 std::filesystem::path referenceFinitePhotonBeamEnergySpreadPositronThetaPath() {
100 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
101 / "cain_linear_breit_wheeler_finite_photon_beam_energy_spread_positron_theta.csv";
102 }
103
104 std::filesystem::path referenceFinitePhotonBeamEnergySpreadElectronJointPath() {
105 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
106 / "cain_linear_breit_wheeler_finite_photon_beam_energy_spread_electron_joint.csv";
107 }
108
109 std::filesystem::path referenceFinitePhotonBeamEnergySpreadPositronJointPath() {
110 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
111 / "cain_linear_breit_wheeler_finite_photon_beam_energy_spread_positron_joint.csv";
112 }
113
114 std::filesystem::path referenceOverlapElectronEnergyPath() {
115 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
116 / "cain_linear_breit_wheeler_overlap_electron_energy.csv";
117 }
118
119 std::filesystem::path referenceOverlapPositronEnergyPath() {
120 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
121 / "cain_linear_breit_wheeler_overlap_positron_energy.csv";
122 }
123
124 std::filesystem::path referenceOverlapElectronThetaPath() {
125 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
126 / "cain_linear_breit_wheeler_overlap_electron_theta.csv";
127 }
128
129 std::filesystem::path referenceOverlapPositronThetaPath() {
130 return std::filesystem::path(OPALX_TEST_SOURCE_DIR) / "data"
131 / "cain_linear_breit_wheeler_overlap_positron_theta.csv";
132 }
133} // namespace
134
135TEST(TestLinearBreitWheelerSpectrum, ElectronEnergySpectrumMatchesCainReference) {
136 const int previousSeed = Options::seed;
137 Options::seed = 13579;
138
139 LinearBreitWheelerBenchmark::HistogramConfig config;
140 config.minValue = 0.0;
141 config.maxValue = 0.5;
142 const auto opalx = LinearBreitWheelerBenchmark::sampleHistogram(
143 config, LinearBreitWheelerBenchmark::FinalState::Electron,
144 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
145 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(referenceElectronEnergyPath());
146
147 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
148 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 1.5e-2);
149 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
150 EXPECT_NEAR(
151 LinearBreitWheelerBenchmark::histogramMean(opalx),
152 LinearBreitWheelerBenchmark::histogramMean(cain),
153 LinearBreitWheelerBenchmark::histogramMean(cain) * 5.0e-2);
154 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.12);
155
156 Options::seed = previousSeed;
157}
158
159TEST(TestLinearBreitWheelerSpectrum, PositronEnergySpectrumMatchesCainReference) {
160 const int previousSeed = Options::seed;
161 Options::seed = 13579;
162
163 LinearBreitWheelerBenchmark::HistogramConfig config;
164 config.minValue = 0.0;
165 config.maxValue = 0.5;
166 const auto opalx = LinearBreitWheelerBenchmark::sampleHistogram(
167 config, LinearBreitWheelerBenchmark::FinalState::Positron,
168 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
169 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(referencePositronEnergyPath());
170
171 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
172 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 1.5e-2);
173 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
174 EXPECT_NEAR(
175 LinearBreitWheelerBenchmark::histogramMean(opalx),
176 LinearBreitWheelerBenchmark::histogramMean(cain),
177 LinearBreitWheelerBenchmark::histogramMean(cain) * 5.0e-2);
178 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.12);
179
180 Options::seed = previousSeed;
181}
182
183TEST(TestLinearBreitWheelerSpectrum, ElectronAngularSpectrumMatchesCainReference) {
184 const int previousSeed = Options::seed;
185 Options::seed = 13579;
186
187 LinearBreitWheelerBenchmark::HistogramConfig config;
188 config.minValue = 0.0;
189 config.maxValue = 0.0045;
190 const auto opalx = LinearBreitWheelerBenchmark::sampleHistogram(
191 config, LinearBreitWheelerBenchmark::FinalState::Electron,
192 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
193 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(referenceElectronThetaPath());
194
195 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
196 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
197 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
198 EXPECT_NEAR(
199 LinearBreitWheelerBenchmark::histogramMean(opalx),
200 LinearBreitWheelerBenchmark::histogramMean(cain),
201 LinearBreitWheelerBenchmark::histogramMean(cain) * 7.0e-2);
202 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.15);
203
204 Options::seed = previousSeed;
205}
206
207TEST(TestLinearBreitWheelerSpectrum, PositronAngularSpectrumMatchesCainReference) {
208 const int previousSeed = Options::seed;
209 Options::seed = 13579;
210
211 LinearBreitWheelerBenchmark::HistogramConfig config;
212 config.minValue = 0.0;
213 config.maxValue = 0.0045;
214 const auto opalx = LinearBreitWheelerBenchmark::sampleHistogram(
215 config, LinearBreitWheelerBenchmark::FinalState::Positron,
216 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
217 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(referencePositronThetaPath());
218
219 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
220 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
221 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
222 EXPECT_NEAR(
223 LinearBreitWheelerBenchmark::histogramMean(opalx),
224 LinearBreitWheelerBenchmark::histogramMean(cain),
225 LinearBreitWheelerBenchmark::histogramMean(cain) * 7.0e-2);
226 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.15);
227
228 Options::seed = previousSeed;
229}
230
231TEST(TestLinearBreitWheelerSpectrum, ElectronJointSpectrumMatchesCainReference) {
232 const int previousSeed = Options::seed;
233 Options::seed = 13579;
234
235 LinearBreitWheelerBenchmark::JointHistogramConfig config;
236 const auto opalx = LinearBreitWheelerBenchmark::sampleJointHistogram(
237 config, LinearBreitWheelerBenchmark::FinalState::Electron, 250000);
238 const auto cain =
239 LinearBreitWheelerBenchmark::readJointHistogramCSV(referenceElectronJointPath());
240
241 ASSERT_EQ(opalx.energyCentersGeV.size(), cain.energyCentersGeV.size());
242 ASSERT_EQ(opalx.thetaCentersRad.size(), cain.thetaCentersRad.size());
243 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(opalx), 1.0, 2.0e-2);
244 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(cain), 1.0, 2.0e-3);
245 EXPECT_NEAR(
246 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(opalx),
247 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain),
248 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain) * 5.0e-2);
249 EXPECT_NEAR(
250 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(opalx),
251 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain),
252 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain) * 7.0e-2);
253 EXPECT_LT(LinearBreitWheelerBenchmark::jointHistogramL1Distance(opalx, cain), 0.16);
254
255 Options::seed = previousSeed;
256}
257
258TEST(TestLinearBreitWheelerSpectrum, PositronJointSpectrumMatchesCainReference) {
259 const int previousSeed = Options::seed;
260 Options::seed = 13579;
261
262 LinearBreitWheelerBenchmark::JointHistogramConfig config;
263 const auto opalx = LinearBreitWheelerBenchmark::sampleJointHistogram(
264 config, LinearBreitWheelerBenchmark::FinalState::Positron, 250000);
265 const auto cain =
266 LinearBreitWheelerBenchmark::readJointHistogramCSV(referencePositronJointPath());
267
268 ASSERT_EQ(opalx.energyCentersGeV.size(), cain.energyCentersGeV.size());
269 ASSERT_EQ(opalx.thetaCentersRad.size(), cain.thetaCentersRad.size());
270 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(opalx), 1.0, 2.0e-2);
271 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(cain), 1.0, 2.0e-3);
272 EXPECT_NEAR(
273 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(opalx),
274 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain),
275 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain) * 5.0e-2);
276 EXPECT_NEAR(
277 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(opalx),
278 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain),
279 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain) * 7.0e-2);
280 EXPECT_LT(LinearBreitWheelerBenchmark::jointHistogramL1Distance(opalx, cain), 0.16);
281
282 Options::seed = previousSeed;
283}
284
285TEST(TestLinearBreitWheelerSpectrum, FinitePhotonBeamElectronEnergySpectrumMatchesCainReference) {
286 const int previousSeed = Options::seed;
287 Options::seed = 13579;
288
289 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
290 config.minValue = 0.0;
291 config.maxValue = 0.5;
292 config.sigmaThetaXRad = 1.0e-3;
293 config.sigmaThetaYRad = 1.0e-3;
294 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
295 config, LinearBreitWheelerBenchmark::FinalState::Electron,
296 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
297 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
298 referenceFinitePhotonBeamElectronEnergyPath());
299
300 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
301 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
302 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
303 EXPECT_NEAR(
304 LinearBreitWheelerBenchmark::histogramMean(opalx),
305 LinearBreitWheelerBenchmark::histogramMean(cain),
306 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
307 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.18);
308
309 Options::seed = previousSeed;
310}
311
312TEST(TestLinearBreitWheelerSpectrum, FinitePhotonBeamPositronEnergySpectrumMatchesCainReference) {
313 const int previousSeed = Options::seed;
314 Options::seed = 13579;
315
316 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
317 config.minValue = 0.0;
318 config.maxValue = 0.5;
319 config.sigmaThetaXRad = 1.0e-3;
320 config.sigmaThetaYRad = 1.0e-3;
321 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
322 config, LinearBreitWheelerBenchmark::FinalState::Positron,
323 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
324 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
325 referenceFinitePhotonBeamPositronEnergyPath());
326
327 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
328 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
329 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
330 EXPECT_NEAR(
331 LinearBreitWheelerBenchmark::histogramMean(opalx),
332 LinearBreitWheelerBenchmark::histogramMean(cain),
333 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
334 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.18);
335
336 Options::seed = previousSeed;
337}
338
339TEST(TestLinearBreitWheelerSpectrum, FinitePhotonBeamElectronAngularSpectrumMatchesCainReference) {
340 const int previousSeed = Options::seed;
341 Options::seed = 13579;
342
343 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
344 config.minValue = 0.0;
345 config.maxValue = 0.0060;
346 config.sigmaThetaXRad = 1.0e-3;
347 config.sigmaThetaYRad = 1.0e-3;
348 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
349 config, LinearBreitWheelerBenchmark::FinalState::Electron,
350 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
351 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
352 referenceFinitePhotonBeamElectronThetaPath());
353
354 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
355 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.5e-2);
356 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 4.0e-3);
357 EXPECT_NEAR(
358 LinearBreitWheelerBenchmark::histogramMean(opalx),
359 LinearBreitWheelerBenchmark::histogramMean(cain),
360 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
361 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.20);
362
363 Options::seed = previousSeed;
364}
365
366TEST(TestLinearBreitWheelerSpectrum, FinitePhotonBeamPositronAngularSpectrumMatchesCainReference) {
367 const int previousSeed = Options::seed;
368 Options::seed = 13579;
369
370 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
371 config.minValue = 0.0;
372 config.maxValue = 0.0060;
373 config.sigmaThetaXRad = 1.0e-3;
374 config.sigmaThetaYRad = 1.0e-3;
375 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
376 config, LinearBreitWheelerBenchmark::FinalState::Positron,
377 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
378 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
379 referenceFinitePhotonBeamPositronThetaPath());
380
381 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
382 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.5e-2);
383 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 4.0e-3);
384 EXPECT_NEAR(
385 LinearBreitWheelerBenchmark::histogramMean(opalx),
386 LinearBreitWheelerBenchmark::histogramMean(cain),
387 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
388 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.20);
389
390 Options::seed = previousSeed;
391}
392
393TEST(TestLinearBreitWheelerSpectrum, FinitePhotonBeamElectronJointSpectrumMatchesCainReference) {
394 const int previousSeed = Options::seed;
395 Options::seed = 13579;
396
397 LinearBreitWheelerBenchmark::FinitePhotonBeamJointConfig config;
398 config.energyMinGeV = 0.0;
399 config.energyMaxGeV = 0.5;
400 config.thetaMinRad = 0.0;
401 config.thetaMaxRad = 0.0060;
402 config.sigmaThetaXRad = 1.0e-3;
403 config.sigmaThetaYRad = 1.0e-3;
404 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamJointHistogram(
405 config, LinearBreitWheelerBenchmark::FinalState::Electron, 250000);
406 const auto cain = LinearBreitWheelerBenchmark::readJointHistogramCSV(
407 referenceFinitePhotonBeamElectronJointPath());
408
409 ASSERT_EQ(opalx.energyCentersGeV.size(), cain.energyCentersGeV.size());
410 ASSERT_EQ(opalx.thetaCentersRad.size(), cain.thetaCentersRad.size());
411 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(opalx), 1.0, 2.5e-2);
412 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(cain), 1.0, 4.0e-3);
413 EXPECT_NEAR(
414 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(opalx),
415 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain),
416 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain) * 8.0e-2);
417 EXPECT_NEAR(
418 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(opalx),
419 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain),
420 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain) * 8.5e-2);
421 EXPECT_LT(LinearBreitWheelerBenchmark::jointHistogramL1Distance(opalx, cain), 0.46);
422
423 Options::seed = previousSeed;
424}
425
426TEST(TestLinearBreitWheelerSpectrum, FinitePhotonBeamPositronJointSpectrumMatchesCainReference) {
427 const int previousSeed = Options::seed;
428 Options::seed = 13579;
429
430 LinearBreitWheelerBenchmark::FinitePhotonBeamJointConfig config;
431 config.energyMinGeV = 0.0;
432 config.energyMaxGeV = 0.5;
433 config.thetaMinRad = 0.0;
434 config.thetaMaxRad = 0.0060;
435 config.sigmaThetaXRad = 1.0e-3;
436 config.sigmaThetaYRad = 1.0e-3;
437 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamJointHistogram(
438 config, LinearBreitWheelerBenchmark::FinalState::Positron, 250000);
439 const auto cain = LinearBreitWheelerBenchmark::readJointHistogramCSV(
440 referenceFinitePhotonBeamPositronJointPath());
441
442 ASSERT_EQ(opalx.energyCentersGeV.size(), cain.energyCentersGeV.size());
443 ASSERT_EQ(opalx.thetaCentersRad.size(), cain.thetaCentersRad.size());
444 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(opalx), 1.0, 2.5e-2);
445 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(cain), 1.0, 4.0e-3);
446 EXPECT_NEAR(
447 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(opalx),
448 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain),
449 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain) * 8.0e-2);
450 EXPECT_NEAR(
451 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(opalx),
452 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain),
453 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain) * 8.5e-2);
454 EXPECT_LT(LinearBreitWheelerBenchmark::jointHistogramL1Distance(opalx, cain), 0.46);
455
456 Options::seed = previousSeed;
457}
458
459TEST(TestLinearBreitWheelerSpectrum,
460 FinitePhotonBeamEnergySpreadElectronEnergySpectrumMatchesCainReference) {
461 const int previousSeed = Options::seed;
462 Options::seed = 13579;
463
464 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
465 config.minValue = 0.0;
466 config.maxValue = 0.5;
467 config.sigmaThetaXRad = 1.0e-3;
468 config.sigmaThetaYRad = 1.0e-3;
469 config.relativeEnergySpread = 1.0e-3;
470 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
471 config, LinearBreitWheelerBenchmark::FinalState::Electron,
472 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
473 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
474 referenceFinitePhotonBeamEnergySpreadElectronEnergyPath());
475
476 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
477 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
478 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
479 EXPECT_NEAR(
480 LinearBreitWheelerBenchmark::histogramMean(opalx),
481 LinearBreitWheelerBenchmark::histogramMean(cain),
482 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
483 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.18);
484
485 Options::seed = previousSeed;
486}
487
488TEST(TestLinearBreitWheelerSpectrum,
489 FinitePhotonBeamEnergySpreadPositronEnergySpectrumMatchesCainReference) {
490 const int previousSeed = Options::seed;
491 Options::seed = 13579;
492
493 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
494 config.minValue = 0.0;
495 config.maxValue = 0.5;
496 config.sigmaThetaXRad = 1.0e-3;
497 config.sigmaThetaYRad = 1.0e-3;
498 config.relativeEnergySpread = 1.0e-3;
499 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
500 config, LinearBreitWheelerBenchmark::FinalState::Positron,
501 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
502 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
503 referenceFinitePhotonBeamEnergySpreadPositronEnergyPath());
504
505 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
506 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
507 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 5.0e-4);
508 EXPECT_NEAR(
509 LinearBreitWheelerBenchmark::histogramMean(opalx),
510 LinearBreitWheelerBenchmark::histogramMean(cain),
511 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
512 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.18);
513
514 Options::seed = previousSeed;
515}
516
517TEST(TestLinearBreitWheelerSpectrum,
518 FinitePhotonBeamEnergySpreadElectronAngularSpectrumMatchesCainReference) {
519 const int previousSeed = Options::seed;
520 Options::seed = 13579;
521
522 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
523 config.minValue = 0.0;
524 config.maxValue = 0.0060;
525 config.sigmaThetaXRad = 1.0e-3;
526 config.sigmaThetaYRad = 1.0e-3;
527 config.relativeEnergySpread = 1.0e-3;
528 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
529 config, LinearBreitWheelerBenchmark::FinalState::Electron,
530 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
531 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
532 referenceFinitePhotonBeamEnergySpreadElectronThetaPath());
533
534 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
535 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.5e-2);
536 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 4.0e-3);
537 EXPECT_NEAR(
538 LinearBreitWheelerBenchmark::histogramMean(opalx),
539 LinearBreitWheelerBenchmark::histogramMean(cain),
540 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
541 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.20);
542
543 Options::seed = previousSeed;
544}
545
546TEST(TestLinearBreitWheelerSpectrum,
547 FinitePhotonBeamEnergySpreadPositronAngularSpectrumMatchesCainReference) {
548 const int previousSeed = Options::seed;
549 Options::seed = 13579;
550
551 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
552 config.minValue = 0.0;
553 config.maxValue = 0.0060;
554 config.sigmaThetaXRad = 1.0e-3;
555 config.sigmaThetaYRad = 1.0e-3;
556 config.relativeEnergySpread = 1.0e-3;
557 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
558 config, LinearBreitWheelerBenchmark::FinalState::Positron,
559 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
560 const auto cain = LinearBreitWheelerBenchmark::readHistogramCSV(
561 referenceFinitePhotonBeamEnergySpreadPositronThetaPath());
562
563 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
564 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.5e-2);
565 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(cain), 1.0, 4.0e-3);
566 EXPECT_NEAR(
567 LinearBreitWheelerBenchmark::histogramMean(opalx),
568 LinearBreitWheelerBenchmark::histogramMean(cain),
569 LinearBreitWheelerBenchmark::histogramMean(cain) * 8.0e-2);
570 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.20);
571
572 Options::seed = previousSeed;
573}
574
575TEST(TestLinearBreitWheelerSpectrum,
576 FinitePhotonBeamEnergySpreadElectronJointSpectrumMatchesCainReference) {
577 const int previousSeed = Options::seed;
578 Options::seed = 13579;
579
580 LinearBreitWheelerBenchmark::FinitePhotonBeamJointConfig config;
581 config.energyMinGeV = 0.0;
582 config.energyMaxGeV = 0.5;
583 config.thetaMinRad = 0.0;
584 config.thetaMaxRad = 0.0060;
585 config.sigmaThetaXRad = 1.0e-3;
586 config.sigmaThetaYRad = 1.0e-3;
587 config.relativeEnergySpread = 1.0e-3;
588 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamJointHistogram(
589 config, LinearBreitWheelerBenchmark::FinalState::Electron, 250000);
590 const auto cain = LinearBreitWheelerBenchmark::readJointHistogramCSV(
591 referenceFinitePhotonBeamEnergySpreadElectronJointPath());
592
593 ASSERT_EQ(opalx.energyCentersGeV.size(), cain.energyCentersGeV.size());
594 ASSERT_EQ(opalx.thetaCentersRad.size(), cain.thetaCentersRad.size());
595 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(opalx), 1.0, 2.5e-2);
596 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(cain), 1.0, 4.0e-3);
597 EXPECT_NEAR(
598 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(opalx),
599 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain),
600 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain) * 8.5e-2);
601 EXPECT_NEAR(
602 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(opalx),
603 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain),
604 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain) * 8.5e-2);
605 EXPECT_LT(LinearBreitWheelerBenchmark::jointHistogramL1Distance(opalx, cain), 0.46);
606
607 Options::seed = previousSeed;
608}
609
610TEST(TestLinearBreitWheelerSpectrum,
611 FinitePhotonBeamEnergySpreadPositronJointSpectrumMatchesCainReference) {
612 const int previousSeed = Options::seed;
613 Options::seed = 13579;
614
615 LinearBreitWheelerBenchmark::FinitePhotonBeamJointConfig config;
616 config.energyMinGeV = 0.0;
617 config.energyMaxGeV = 0.5;
618 config.thetaMinRad = 0.0;
619 config.thetaMaxRad = 0.0060;
620 config.sigmaThetaXRad = 1.0e-3;
621 config.sigmaThetaYRad = 1.0e-3;
622 config.relativeEnergySpread = 1.0e-3;
623 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamJointHistogram(
624 config, LinearBreitWheelerBenchmark::FinalState::Positron, 250000);
625 const auto cain = LinearBreitWheelerBenchmark::readJointHistogramCSV(
626 referenceFinitePhotonBeamEnergySpreadPositronJointPath());
627
628 ASSERT_EQ(opalx.energyCentersGeV.size(), cain.energyCentersGeV.size());
629 ASSERT_EQ(opalx.thetaCentersRad.size(), cain.thetaCentersRad.size());
630 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(opalx), 1.0, 2.5e-2);
631 EXPECT_NEAR(LinearBreitWheelerBenchmark::jointHistogramArea(cain), 1.0, 4.0e-3);
632 EXPECT_NEAR(
633 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(opalx),
634 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain),
635 LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV(cain) * 8.5e-2);
636 EXPECT_NEAR(
637 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(opalx),
638 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain),
639 LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad(cain) * 8.5e-2);
640 EXPECT_LT(LinearBreitWheelerBenchmark::jointHistogramL1Distance(opalx, cain), 0.46);
641
642 Options::seed = previousSeed;
643}
644
645TEST(TestLinearBreitWheelerSpectrum, OverlapWeightedElectronEnergySpectrumMatchesCainReference) {
646 const int previousSeed = Options::seed;
647 Options::seed = 13579;
648
649 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
650 config.minValue = 0.0;
651 config.maxValue = 0.5;
652 config.sigmaThetaXRad = 1.0e-3;
653 config.sigmaThetaYRad = 1.0e-3;
654 config.sigmaX_m = 1.0e-6;
655 config.sigmaY_m = 1.0e-6;
656 config.sigmaS_m = 5.0e-12 * Physics::c;
657 config.laserRayleighX_m = 1.0e-6;
658 config.laserRayleighY_m = 1.0e-6;
659 config.laserSigmaT_m = 5.0e-12 * Physics::c;
660 config.overlapWeighting = true;
661 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
662 config, LinearBreitWheelerBenchmark::FinalState::Electron,
663 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
664 const auto cain =
665 LinearBreitWheelerBenchmark::readHistogramCSV(referenceOverlapElectronEnergyPath());
666
667 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
668 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
669 EXPECT_NEAR(
670 LinearBreitWheelerBenchmark::histogramMean(opalx),
671 LinearBreitWheelerBenchmark::histogramMean(cain),
672 LinearBreitWheelerBenchmark::histogramMean(cain) * 1.2e-1);
673 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.28);
674
675 Options::seed = previousSeed;
676}
677
678TEST(TestLinearBreitWheelerSpectrum, OverlapWeightedPositronEnergySpectrumMatchesCainReference) {
679 const int previousSeed = Options::seed;
680 Options::seed = 13579;
681
682 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
683 config.minValue = 0.0;
684 config.maxValue = 0.5;
685 config.sigmaThetaXRad = 1.0e-3;
686 config.sigmaThetaYRad = 1.0e-3;
687 config.sigmaX_m = 1.0e-6;
688 config.sigmaY_m = 1.0e-6;
689 config.sigmaS_m = 5.0e-12 * Physics::c;
690 config.laserRayleighX_m = 1.0e-6;
691 config.laserRayleighY_m = 1.0e-6;
692 config.laserSigmaT_m = 5.0e-12 * Physics::c;
693 config.overlapWeighting = true;
694 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
695 config, LinearBreitWheelerBenchmark::FinalState::Positron,
696 LinearBreitWheelerBenchmark::Observable::Energy, 250000);
697 const auto cain =
698 LinearBreitWheelerBenchmark::readHistogramCSV(referenceOverlapPositronEnergyPath());
699
700 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
701 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.0e-2);
702 EXPECT_NEAR(
703 LinearBreitWheelerBenchmark::histogramMean(opalx),
704 LinearBreitWheelerBenchmark::histogramMean(cain),
705 LinearBreitWheelerBenchmark::histogramMean(cain) * 1.2e-1);
706 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.28);
707
708 Options::seed = previousSeed;
709}
710
711TEST(TestLinearBreitWheelerSpectrum, OverlapWeightedElectronAngularSpectrumMatchesCainReference) {
712 const int previousSeed = Options::seed;
713 Options::seed = 13579;
714
715 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
716 config.minValue = 0.0;
717 config.maxValue = 0.0060;
718 config.sigmaThetaXRad = 1.0e-3;
719 config.sigmaThetaYRad = 1.0e-3;
720 config.sigmaX_m = 1.0e-6;
721 config.sigmaY_m = 1.0e-6;
722 config.sigmaS_m = 5.0e-12 * Physics::c;
723 config.laserRayleighX_m = 1.0e-6;
724 config.laserRayleighY_m = 1.0e-6;
725 config.laserSigmaT_m = 5.0e-12 * Physics::c;
726 config.overlapWeighting = true;
727 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
728 config, LinearBreitWheelerBenchmark::FinalState::Electron,
729 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
730 const auto cain =
731 LinearBreitWheelerBenchmark::readHistogramCSV(referenceOverlapElectronThetaPath());
732
733 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
734 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.5e-2);
735 EXPECT_NEAR(
736 LinearBreitWheelerBenchmark::histogramMean(opalx),
737 LinearBreitWheelerBenchmark::histogramMean(cain),
738 LinearBreitWheelerBenchmark::histogramMean(cain) * 1.2e-1);
739 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.32);
740
741 Options::seed = previousSeed;
742}
743
744TEST(TestLinearBreitWheelerSpectrum, OverlapWeightedPositronAngularSpectrumMatchesCainReference) {
745 const int previousSeed = Options::seed;
746 Options::seed = 13579;
747
748 LinearBreitWheelerBenchmark::FinitePhotonBeamConfig config;
749 config.minValue = 0.0;
750 config.maxValue = 0.0060;
751 config.sigmaThetaXRad = 1.0e-3;
752 config.sigmaThetaYRad = 1.0e-3;
753 config.sigmaX_m = 1.0e-6;
754 config.sigmaY_m = 1.0e-6;
755 config.sigmaS_m = 5.0e-12 * Physics::c;
756 config.laserRayleighX_m = 1.0e-6;
757 config.laserRayleighY_m = 1.0e-6;
758 config.laserSigmaT_m = 5.0e-12 * Physics::c;
759 config.overlapWeighting = true;
760 const auto opalx = LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram(
761 config, LinearBreitWheelerBenchmark::FinalState::Positron,
762 LinearBreitWheelerBenchmark::Observable::Theta, 250000);
763 const auto cain =
764 LinearBreitWheelerBenchmark::readHistogramCSV(referenceOverlapPositronThetaPath());
765
766 ASSERT_EQ(opalx.centers.size(), cain.centers.size());
767 EXPECT_NEAR(LinearBreitWheelerBenchmark::histogramArea(opalx), 1.0, 2.5e-2);
768 EXPECT_NEAR(
769 LinearBreitWheelerBenchmark::histogramMean(opalx),
770 LinearBreitWheelerBenchmark::histogramMean(cain),
771 LinearBreitWheelerBenchmark::histogramMean(cain) * 1.2e-1);
772 EXPECT_LT(LinearBreitWheelerBenchmark::histogramL1Distance(opalx, cain), 0.32);
773
774 Options::seed = previousSeed;
775}
TEST
TEST(TestLinearBreitWheelerSpectrum, ElectronEnergySpectrumMatchesCainReference)
Definition TestLinearBreitWheelerSpectrum.cpp:135
LinearBreitWheelerBenchmark::histogramMean
double histogramMean(const Histogram &histogram)
Definition LinearBreitWheelerBenchmarkCommon.h:797
LinearBreitWheelerBenchmark::histogramL1Distance
double histogramL1Distance(const Histogram &lhs, const Histogram &rhs)
Definition LinearBreitWheelerBenchmarkCommon.h:805
LinearBreitWheelerBenchmark::sampleFinitePhotonBeamJointHistogram
JointHistogram sampleFinitePhotonBeamJointHistogram(const FinitePhotonBeamJointConfig &config, FinalState state, std::size_t sampleCount, std::uint64_t streamIndex=0)
Sample the joint Breit-Wheeler laboratory distribution in energy and polar angle.
Definition LinearBreitWheelerBenchmarkCommon.h:465
LinearBreitWheelerBenchmark::sampleFinitePhotonBeamHistogram
Histogram sampleFinitePhotonBeamHistogram(const FinitePhotonBeamConfig &config, FinalState state, Observable observable, std::size_t sampleCount, std::uint64_t streamIndex=0)
Sample a one-dimensional Breit-Wheeler histogram for a finite incoming photon beam.
Definition LinearBreitWheelerBenchmarkCommon.h:327
LinearBreitWheelerBenchmark::sampleHistogram
Histogram sampleHistogram(const HistogramConfig &config, FinalState state, Observable observable, std::size_t sampleCount, std::uint64_t streamIndex=0)
Sample a one-dimensional Breit-Wheeler benchmark histogram.
Definition LinearBreitWheelerBenchmarkCommon.h:403
LinearBreitWheelerBenchmark::sampleJointHistogram
JointHistogram sampleJointHistogram(const JointHistogramConfig &config, FinalState state, std::size_t sampleCount, std::uint64_t streamIndex=0)
Definition LinearBreitWheelerBenchmarkCommon.h:554
LinearBreitWheelerBenchmark::jointHistogramMeanEnergyGeV
double jointHistogramMeanEnergyGeV(const JointHistogram &histogram)
Definition LinearBreitWheelerBenchmarkCommon.h:826
LinearBreitWheelerBenchmark::histogramArea
double histogramArea(const Histogram &histogram)
Definition LinearBreitWheelerBenchmarkCommon.h:789
LinearBreitWheelerBenchmark::readJointHistogramCSV
JointHistogram readJointHistogramCSV(const std::filesystem::path &inputPath)
Definition LinearBreitWheelerBenchmarkCommon.h:701
LinearBreitWheelerBenchmark::jointHistogramL1Distance
double jointHistogramL1Distance(const JointHistogram &lhs, const JointHistogram &rhs)
Definition LinearBreitWheelerBenchmarkCommon.h:848
LinearBreitWheelerBenchmark::jointHistogramArea
double jointHistogramArea(const JointHistogram &histogram)
Definition LinearBreitWheelerBenchmarkCommon.h:816
LinearBreitWheelerBenchmark::readHistogramCSV
Histogram readHistogramCSV(const std::filesystem::path &inputPath)
Definition LinearBreitWheelerBenchmarkCommon.h:664
LinearBreitWheelerBenchmark::jointHistogramMeanThetaRad
double jointHistogramMeanThetaRad(const JointHistogram &histogram)
Definition LinearBreitWheelerBenchmarkCommon.h:837
LinearBreitWheelerBenchmark::FinitePhotonBeamConfig
Configuration for folding the linear Breit-Wheeler kernel over a finite incoming photon beam.
Definition LinearBreitWheelerBenchmarkCommon.h:41
Options::seed
int seed
The current random seed.
Definition Options.cpp:37
Physics::c
constexpr double c
The velocity of light in m/s.
Definition Physics.h:60