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
|
import static org.junit.jupiter.api.Assertions.assertEquals;
import static org.junit.jupiter.api.Assertions.assertThrows;
import org.junit.jupiter.api.Test;
import com.github.tschoonj.xraylib.Xraylib;
public class TestComptonProfiles {
@Test
public void test_pz_0() {
double profile = Xraylib.ComptonProfile(26, 0.0);
assertEquals(profile, 7.060, 1E-6);
profile = Xraylib.ComptonProfile_Partial(26, Xraylib.N1_SHELL, 0.0);
assertEquals(profile, 1.550, 1E-6);
double profile1 = Xraylib.ComptonProfile_Partial(26, Xraylib.L2_SHELL, 0.0);
double profile2 = Xraylib.ComptonProfile_Partial(26, Xraylib.L3_SHELL, 0.0);
assertEquals(profile1, profile2, 1E-6);
assertEquals(profile1, 0.065, 1E-6);
}
@Test
public void test_pz_100() {
double profile = Xraylib.ComptonProfile(26, 100.0);
assertEquals(profile, 1.8E-5, 1E-8);
profile = Xraylib.ComptonProfile_Partial(26, Xraylib.N1_SHELL, 100.0);
assertEquals(profile, 5.1E-9, 1E-12);
double profile1 = Xraylib.ComptonProfile_Partial(26, Xraylib.L2_SHELL, 100.0);
double profile2 = Xraylib.ComptonProfile_Partial(26, Xraylib.L3_SHELL, 100.0);
assertEquals(profile1, profile2, 1E-10);
assertEquals(profile1, 1.1E-8, 1E-10);
}
@Test
public void test_pz_50() {
double profile = Xraylib.ComptonProfile(26, 50.0);
assertEquals(profile, 0.0006843950273082384, 1E-8);
profile = Xraylib.ComptonProfile_Partial(26, Xraylib.N1_SHELL, 50.0);
assertEquals(profile, 2.4322755767709126e-07, 1E-12);
double profile1 = Xraylib.ComptonProfile_Partial(26, Xraylib.L2_SHELL, 50.0);
double profile2 = Xraylib.ComptonProfile_Partial(26, Xraylib.L3_SHELL, 50.0);
assertEquals(profile1, profile2, 1E-10);
assertEquals(profile1, 2.026953933016568e-06, 1E-10);
}
@Test
public void test_bad_input() {
IllegalArgumentException exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile(0, 0.0);
});
assertEquals(exc.getMessage(), Xraylib.Z_OUT_OF_RANGE);
Xraylib.ComptonProfile(102, 0.0);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile(103, 0.0);
});
assertEquals(exc.getMessage(), Xraylib.Z_OUT_OF_RANGE);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile(26, -1.0);
});
assertEquals(exc.getMessage(), Xraylib.NEGATIVE_PZ);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile(26, 101);
});
assertEquals(exc.getMessage(), Xraylib.SPLINT_X_TOO_HIGH);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile_Partial(0, Xraylib.K_SHELL, 0.0);
});
assertEquals(exc.getMessage(), Xraylib.Z_OUT_OF_RANGE);
Xraylib.ComptonProfile_Partial(102, Xraylib.K_SHELL, 0.0);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile_Partial(103, Xraylib.K_SHELL, 0.0);
});
assertEquals(exc.getMessage(), Xraylib.Z_OUT_OF_RANGE);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile_Partial(26, Xraylib.K_SHELL, -1.0);
});
assertEquals(exc.getMessage(), Xraylib.NEGATIVE_PZ);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile_Partial(26, -1, 0.0);
});
assertEquals(exc.getMessage(), Xraylib.INVALID_SHELL);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile_Partial(26, Xraylib.N2_SHELL, 0.0);
});
assertEquals(exc.getMessage(), Xraylib.INVALID_SHELL);
exc = assertThrows(IllegalArgumentException.class, () -> {
double profile = Xraylib.ComptonProfile_Partial(26, Xraylib.K_SHELL, 101);
});
assertEquals(exc.getMessage(), Xraylib.SPLINT_X_TOO_HIGH);
}
}
|
