1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
|
import static org.junit.jupiter.api.Assertions.assertEquals;
import static org.junit.jupiter.api.Assertions.assertArrayEquals;
import static org.junit.jupiter.api.Assertions.assertThrows;
import java.util.stream.Stream;
import org.junit.jupiter.api.Test;
import org.junit.jupiter.api.BeforeAll;
import org.junit.jupiter.api.DisplayName;
import org.junit.jupiter.params.ParameterizedTest;
import org.junit.jupiter.params.provider.ValueSource;
import org.junit.jupiter.params.provider.MethodSource;
import org.junit.jupiter.params.provider.NullAndEmptySource;
import org.junit.jupiter.params.provider.Arguments;
import static org.junit.jupiter.params.provider.Arguments.arguments;
import com.github.tschoonj.xraylib.Xraylib;
public class TestCrossSectionsFluorLine {
static Stream<Arguments> goodValuesProvider() {
return Stream.of(
arguments(29, Xraylib.L3M5_LINE, 10.0, 0.13198670698075143),
arguments(29, Xraylib.L1M5_LINE, 10.0, 7.723944209880828e-06),
arguments(29, Xraylib.L1M2_LINE, 10.0, 0.0018343168459245755),
arguments(29, Xraylib.KL3_LINE, 10.0, 49.21901698835919),
arguments(26, Xraylib.L2M4_LINE, 10.0, 0.0200667),
arguments(26, Xraylib.L1M2_LINE, 10.0, 0.000830915),
arguments(29, Xraylib.KL3_LINE, 999.0, 1.1827911388054846E-4)
);
}
@ParameterizedTest(name="test_good_values {index} -> {0} {1} {2}")
@MethodSource("goodValuesProvider")
public void test_good_values(int Z, int line, double energy, double expected) {
double cs = Xraylib.CS_FluorLine(Z, line, energy);
assertEquals(cs, expected, 1E-6);
}
static Stream<Arguments> badValuesProvider() {
return Stream.of(
arguments(26, Xraylib.L3M5_LINE, 10.0, Xraylib.UNAVAILABLE_CK),
arguments(0, Xraylib.KL3_LINE, 10.0, Xraylib.Z_OUT_OF_RANGE),
arguments(Xraylib.ZMAX, Xraylib.KL3_LINE, 10.0, Xraylib.INVALID_LINE),
arguments(Xraylib.ZMAX + 1, Xraylib.KL3_LINE, 10.0, Xraylib.Z_OUT_OF_RANGE),
arguments(1, Xraylib.KL3_LINE, 10.0, Xraylib.INVALID_LINE),
arguments(92, Xraylib.M5N7_LINE, 10.0, Xraylib.INVALID_LINE),
arguments(26, Xraylib.KL3_LINE, 0.0, Xraylib.NEGATIVE_ENERGY),
arguments(92, Xraylib.L3M5_LINE, 10.0, Xraylib.TOO_LOW_EXCITATION_ENERGY),
arguments(26, Xraylib.KL3_LINE, 1001.0, Xraylib.SPLINT_X_TOO_HIGH)
);
}
@ParameterizedTest(name="test_bad_values {index} -> {0} {1} {2}")
@MethodSource("badValuesProvider")
public void test_bad_values(int Z, int line, double energy, String message) {
IllegalArgumentException exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = Xraylib.CS_FluorLine(Z, line, energy);
}, message);
assertEquals(exc.getMessage(), message);
}
@Test
public void test_deprecated_siegbahn_macros() {
// KA_LINE
double cs = Xraylib.CS_FluorLine(26, Xraylib.KL3_LINE, 10.0);
cs += Xraylib.CS_FluorLine(26, Xraylib.KL2_LINE, 10.0);
assertEquals(cs, Xraylib.CS_FluorLine(26, Xraylib.KA_LINE, 10.0), 1E-6);
// LA_LINE
cs = Xraylib.CS_FluorLine(92, Xraylib.L3M5_LINE, 30.0);
cs += Xraylib.CS_FluorLine(92, Xraylib.L3M4_LINE, 30.0);
assertEquals(cs, Xraylib.CS_FluorLine(92, Xraylib.LA_LINE, 30.0), 1E-6);
// LB_LINE
final int[] lb_line_macros = new int[]{
Xraylib.LB1_LINE, Xraylib.LB2_LINE, Xraylib.LB3_LINE, Xraylib.LB4_LINE,
Xraylib.LB5_LINE, Xraylib.LB6_LINE, Xraylib.LB7_LINE, Xraylib.LB9_LINE,
Xraylib.LB10_LINE, Xraylib.LB15_LINE, Xraylib.LB17_LINE, Xraylib.L3N6_LINE,
Xraylib.L3N7_LINE
};
cs = 0.0;
for (int line: lb_line_macros) {
cs += Xraylib.CS_FluorLine(92, line, 30.0);
}
assertEquals(cs, Xraylib.CS_FluorLine(92, Xraylib.LB_LINE, 30.0), 1E-6);
}
private static class LineMapping {
public final int line_lower;
public final int line_upper;
public final int shell;
public LineMapping(final int line_lower, final int line_upper, final int shell) {
this.line_lower = line_lower;
this.line_upper = line_upper;
this.shell = shell;
}
}
private static final LineMapping[] line_mappings = new LineMapping[]{
new LineMapping(Xraylib.KN5_LINE, Xraylib.KB_LINE, Xraylib.K_SHELL),
new LineMapping(Xraylib.L1P5_LINE, Xraylib.L1L2_LINE, Xraylib.L1_SHELL),
new LineMapping(Xraylib.L2Q1_LINE, Xraylib.L2L3_LINE, Xraylib.L2_SHELL),
new LineMapping(Xraylib.L3Q1_LINE, Xraylib.L3M1_LINE, Xraylib.L3_SHELL),
};
@Test
public void test_fluor_shell_all() {
for (LineMapping mapping: line_mappings) {
double cs = Xraylib.CS_FluorShell(92, mapping.shell, 120.0);
double cs2 = 0;
double rr = 0;
for (int j = mapping.line_lower ; j <= mapping.line_upper ; j++) {
try {
rr += Xraylib.RadRate(92, j);
cs2 += Xraylib.CS_FluorLine(92, j, 120.0);
} catch (IllegalArgumentException e) {
continue;
}
}
assertEquals(cs2, rr * cs, 1E-6);
}
}
static Stream<Arguments> badShellValuesProvider() {
return Stream.of(
arguments(0, Xraylib.K_SHELL, 10.0, Xraylib.Z_OUT_OF_RANGE),
arguments(Xraylib.ZMAX, Xraylib.K_SHELL, 10.0, Xraylib.INVALID_SHELL),
arguments(Xraylib.ZMAX + 1, Xraylib.K_SHELL, 10.0, Xraylib.Z_OUT_OF_RANGE),
arguments(1, Xraylib.K_SHELL, 10.0, Xraylib.INVALID_SHELL),
arguments(92, Xraylib.M1_SHELL, 10.0, Xraylib.INVALID_SHELL),
arguments(92, Xraylib.KL3_LINE, 10.0, Xraylib.INVALID_SHELL),
arguments(26, Xraylib.K_SHELL, 0.0, Xraylib.NEGATIVE_ENERGY),
arguments(26, Xraylib.K_SHELL, 1001.0, Xraylib.SPLINT_X_TOO_HIGH),
arguments(26, Xraylib.K_SHELL, 5, Xraylib.TOO_LOW_EXCITATION_ENERGY)
);
}
@ParameterizedTest(name="test_bad_shell_values {index} -> {0} {1} {2}")
@MethodSource("badShellValuesProvider")
public void test_bad_shell_values(int Z, int shell, double energy, String message) {
IllegalArgumentException exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = Xraylib.CS_FluorShell(Z, shell, energy);
}, message);
assertEquals(exc.getMessage(), message);
}
}
|
