Version: SMASH-3.1
parametrizations_data.h
Go to the documentation of this file.
1 /*
2  *
3  * Copyright (c) 2013-2020,2023-2024
4  * SMASH Team
5  *
6  * GNU General Public License (GPLv3 or later)
7  *
8  */
9 #ifndef SRC_INCLUDE_SMASH_PARAMETRIZATIONS_DATA_H_
10 #define SRC_INCLUDE_SMASH_PARAMETRIZATIONS_DATA_H_
11 
12 #include <initializer_list>
13 #include <memory>
14 
15 #include "interpolation.h"
16 
17 namespace smash {
18 
20 const std::initializer_list<double> KMINUSN_TOT_PLAB = {
21  0.627, 0.728333, 0.846, 0.924, 1.01033, 1.09367, 1.259, 1.48,
22  1.72, 2.05333, 2.46667, 2.81667, 3.31667, 4.32667, 5.54333, 6.88333,
23  8.21667, 9.33333, 10.6667, 11.7667, 13.1, 14.1, 15., 16.3333,
24  18., 19.3333, 21.6667, 23.3333, 26.6667, 28.3333, 31.6667, 33.3333,
25  35., 36.6667, 38.3333, 41.6667, 43.3333, 46.6667, 48.3333, 50.,
26  50., 51.6667, 53.3333, 60., 75., 90., 106.667, 123.333,
27  140., 156.667, 173.333, 190., 213.333, 240., 276.667, 280.,
28  310.};
29 
31 static std::unique_ptr<InterpolateDataLinear<double>>
33 
35 const std::initializer_list<double> KMINUSN_TOT_SIG = {
36  26.2, 29.1333, 30.8, 33.9667, 36.1667, 35.5667, 30.8, 26.7,
37  24.4667, 23.2, 22.5667, 22.2667, 21.5667, 21.4, 21.4667, 21.2,
38  20.9333, 20.4667, 20.6333, 20.5667, 20.4, 20.3033, 20.1033, 20.17,
39  19.9, 19.75, 19.55, 19.6667, 19.5267, 19.66, 19.81, 19.8233,
40  19.7367, 19.57, 19.5333, 19.4067, 19.6733, 19.8033, 19.7867, 19.6867,
41  19.7233, 19.7633, 19.93, 19.92, 20.0267, 19.88, 19.99, 20.08,
42  20.0967, 20.1267, 20.26, 20.3567, 20.4567, 20.5833, 20.84, 20.9,
43  21.07};
44 
46 const std::initializer_list<double> KMINUSP_ELASTIC_P_LAB = {
47  0.03000, 0.05000, 0.06287, 0.07000, 0.07044, 0.07300, 0.08730,
48  0.09000, 0.11000, 0.11000, 0.11210, 0.11300, 0.12262, 0.13000,
49  0.13000, 0.13722, 0.14552, 0.15000, 0.15000, 0.15300, 0.15909,
50  0.16269, 0.17000, 0.17000, 0.17470, 0.18768, 0.18916, 0.19000,
51  0.19000, 0.19300, 0.20534, 0.21000, 0.21000, 0.21554, 0.22049,
52  0.22500, 0.23000, 0.23000, 0.23300, 0.23500, 0.23944, 0.24500,
53  0.24623, 0.25000, 0.25500, 0.26154, 0.26500, 0.27000, 0.27500,
54  0.27618, 0.28227, 0.28500, 0.29000, 0.29290, 0.29300, 0.29500,
55  0.30000, 0.30191, 0.30500, 0.31500, 0.32500, 0.33500, 0.34500,
56  0.35000, 0.35012, 0.35500, 0.36500, 0.37500, 0.38500, 0.38700,
57  0.38986, 0.39200, 0.39500, 0.40000, 0.40500, 0.41500, 0.42000,
58  0.42500, 0.43400, 0.43407, 0.43500, 0.43600, 0.44500, 0.45500,
59  0.45500, 0.46500, 0.47500, 0.49500, 0.51300, 0.51300, 0.51400,
60  0.52228, 0.53400, 0.55400, 0.57300, 0.59700, 0.61000, 0.61700,
61  0.62004, 0.63700, 0.64200, 0.65800, 0.66000, 0.67200, 0.67700,
62  0.69900, 0.70000, 0.70300, 0.70800, 0.71900, 0.72500, 0.73000,
63  0.74000, 0.74000, 0.74100, 0.75800, 0.75996, 0.76100, 0.76800,
64  0.77300, 0.77700, 0.77700, 0.78000, 0.78500, 0.79300, 0.80200,
65  0.80600, 0.80600, 0.81000, 0.82000, 0.82000, 0.83300, 0.83800,
66  0.83800, 0.84999, 0.85300, 0.85300, 0.85600, 0.86000, 0.87000,
67  0.87400, 0.87400, 0.87600, 0.89400, 0.89400, 0.89900, 0.90000,
68  0.90400, 0.90400, 0.90500, 0.91600, 0.91600, 0.92200, 0.92500,
69  0.93500, 0.93500, 0.93500, 0.94000, 0.94300, 0.94500, 0.95400,
70  0.95400, 0.95500, 0.96000, 0.96500, 0.97000, 0.97000, 0.98000,
71  0.99000, 0.99100, 0.99100, 1.00500, 1.02000, 1.02200, 1.02200,
72  1.04000, 1.04400, 1.04400, 1.04500, 1.06000, 1.06100, 1.06100,
73  1.08000, 1.08000, 1.08500, 1.10000, 1.10200, 1.10200, 1.11700,
74  1.11700, 1.12500, 1.13400, 1.13400, 1.13800, 1.14000, 1.15000,
75  1.15000, 1.15300, 1.15300, 1.16100, 1.16500, 1.17400, 1.17400,
76  1.17900, 1.18000, 1.18300, 1.18300, 1.20100, 1.20500, 1.22000,
77  1.22600, 1.22600, 1.23300, 1.24500, 1.25300, 1.26000, 1.26000,
78  1.26300, 1.26300, 1.27600, 1.28500, 1.29600, 1.30000, 1.31600,
79  1.31600, 1.32000, 1.32800, 1.34000, 1.35500, 1.36800, 1.36800,
80  1.38000, 1.38300, 1.41500, 1.41500, 1.42300, 1.43300, 1.46200,
81  1.46500, 1.48300, 1.51300, 1.51400, 1.53000, 1.53400, 1.54500,
82  1.54600, 1.58400, 1.60600, 1.60600, 1.63400, 1.65200, 1.65300,
83  1.68000, 1.68400, 1.70500, 1.70500, 1.73400, 1.73900, 1.74100,
84  1.78400, 1.80000, 1.80000, 1.81500, 1.84300, 1.84300, 1.88400,
85  1.93400, 1.93400, 1.98400, 2.00000, 2.03100, 2.03400, 2.08400,
86  2.13400, 2.13500, 2.17500, 2.23400, 2.24000, 2.28400, 2.32500,
87  2.33100, 2.37400, 2.41200, 2.51600, 2.66000, 2.66000, 2.98500,
88  3.00000, 3.00000, 3.46000, 3.59000, 3.65000, 4.10000, 4.20000,
89  4.60000, 5.00000, 5.50000, 6.00000, 7.20020, 9.00010, 10.12000,
90  14.30000, 14.30000, 25.20000, 32.10000, 40.10000, 50.00000, 70.00000,
91  100.00000, 140.00000, 147.00000, 175.00000};
92 
94 const std::initializer_list<double> KMINUSP_ELASTIC_SIG = {
95  313.50, 103.60, 113.00, 44.800, 58.500, 187.00, 92.000, 71.500, 92.800,
96  87.290, 82.000, 105.00, 59.400, 40.400, 79.220, 82.000, 49.000, 41.400,
97  69.610, 108.00, 53.900, 98.000, 75.760, 32.800, 45.000, 73.000, 66.000,
98  59.090, 53.300, 68.000, 37.000, 53.300, 60.490, 48.000, 41.000, 65.410,
99  62.900, 55.690, 50.000, 55.740, 41.200, 53.240, 37.000, 51.500, 49.220,
100  43.600, 47.710, 58.060, 48.750, 30.000, 44.900, 39.420, 38.270, 47.800,
101  48.200, 41.220, 44.500, 44.200, 40.360, 37.020, 40.280, 37.840, 37.260,
102  34.000, 33.500, 34.770, 34.210, 36.670, 33.890, 31.900, 34.700, 34.000,
103  38.470, 38.900, 32.060, 32.590, 48.400, 31.190, 30.600, 32.800, 26.830,
104  25.800, 28.830, 23.800, 30.320, 31.990, 23.100, 21.500, 26.500, 27.600,
105  21.700, 35.000, 19.300, 19.100, 17.500, 17.700, 17.660, 18.600, 16.000,
106  16.000, 17.230, 16.400, 12.100, 16.220, 15.600, 15.200, 14.200, 15.220,
107  13.500, 14.200, 11.500, 14.070, 15.900, 14.000, 11.600, 16.930, 16.700,
108  17.300, 15.200, 18.600, 18.300, 18.700, 17.900, 19.370, 20.500, 19.300,
109  20.000, 20.450, 20.670, 18.700, 19.300, 19.630, 19.800, 19.530, 22.400,
110  19.100, 19.390, 19.780, 19.500, 21.360, 20.100, 20.310, 21.070, 21.600,
111  21.660, 22.020, 21.500, 20.900, 20.710, 21.340, 20.850, 20.400, 22.340,
112  20.940, 20.120, 20.100, 19.980, 21.800, 21.010, 19.330, 20.640, 20.700,
113  20.610, 21.270, 20.340, 20.400, 20.720, 22.400, 21.220, 21.400, 22.150,
114  21.750, 20.800, 22.100, 22.150, 23.300, 21.500, 21.460, 22.220, 21.200,
115  20.600, 20.560, 18.700, 18.740, 19.830, 18.300, 18.600, 18.770, 17.600,
116  17.820, 17.890, 17.000, 17.750, 15.700, 17.000, 18.300, 17.300, 17.200,
117  17.230, 15.300, 15.390, 16.500, 16.460, 14.660, 16.800, 15.900, 15.890,
118  12.680, 13.890, 15.700, 11.300, 11.810, 13.110, 12.320, 11.870, 15.200,
119  14.000, 10.870, 10.870, 11.440, 10.960, 11.570, 12.000, 10.200, 11.200,
120  10.260, 9.7400, 14.400, 9.5300, 10.300, 10.300, 16.600, 10.500, 8.8300,
121  8.8300, 8.4200, 8.6000, 9.1100, 9.1100, 8.3000, 7.7000, 7.7000, 8.1500,
122  8.8000, 8.0600, 8.0600, 8.7000, 8.8600, 8.8600, 9.2000, 8.4000, 8.4000,
123  8.1900, 8.9000, 9.0800, 9.0800, 9.5000, 8.5100, 8.5100, 9.0000, 8.1300,
124  8.1300, 6.9500, 7.8600, 7.8600, 9.0000, 7.9000, 7.4800, 7.8000, 7.4600,
125  7.3100, 7.1000, 7.4000, 7.1000, 7.2600, 7.9000, 6.8700, 6.2000, 7.3000,
126  6.5000, 6.7600, 6.6000, 6.3200, 6.2100, 5.7000, 6.0400, 4.9500, 5.0600,
127  4.9500, 4.9400, 4.4000, 4.6000, 4.3000, 4.5000, 4.2000, 3.8400, 4.1000,
128  3.6200, 4.2300, 3.9500, 3.2400, 2.9600, 3.0100, 2.4600, 2.5600, 2.3300,
129  2.5400, 2.5300, 2.5100, 2.5200, 2.7400, 2.5900};
130 
132 static std::unique_ptr<InterpolateDataLinear<double>>
134 
136 const std::initializer_list<double> KMINUSP_TOT_PLAB = {
137  0.245, 0.255, 0.265, 0.275, 0.285, 0.293, 0.293, 0.295,
138  0.300, 0.305, 0.315, 0.325, 0.335, 0.345, 0.350, 0.350,
139  0.355, 0.360, 0.365, 0.375, 0.385, 0.385, 0.387, 0.390,
140  0.392, 0.395, 0.400, 0.405, 0.405, 0.415, 0.419, 0.425,
141  0.434, 0.434, 0.435, 0.436, 0.440, 0.445, 0.446, 0.455,
142  0.455, 0.464, 0.475, 0.475, 0.483, 0.494, 0.495, 0.506,
143  0.506, 0.513, 0.513, 0.514, 0.520, 0.534, 0.534, 0.536,
144  0.547, 0.554, 0.561, 0.566, 0.573, 0.574, 0.586, 0.590,
145  0.596, 0.597, 0.599, 0.604, 0.612, 0.617, 0.620, 0.623,
146  0.625, 0.627, 0.630, 0.637, 0.639, 0.653, 0.655, 0.657,
147  0.658, 0.667, 0.677, 0.680, 0.686, 0.687, 0.693, 0.699,
148  0.707, 0.708, 0.712, 0.719, 0.720, 0.725, 0.732, 0.736,
149  0.740, 0.741, 0.744, 0.752, 0.756, 0.760, 0.761, 0.767,
150  0.768, 0.768, 0.773, 0.777, 0.782, 0.793, 0.794, 0.795,
151  0.802, 0.806, 0.819, 0.820, 0.822, 0.838, 0.839, 0.840,
152  0.848, 0.850, 0.853, 0.860, 0.874, 0.874, 0.881, 0.894,
153  0.894, 0.899, 0.900, 0.900, 0.904, 0.916, 0.925, 0.931,
154  0.935, 0.936, 0.952, 0.954, 0.965, 0.969, 0.970, 0.975,
155  0.978, 0.980, 0.987, 0.991, 0.991, 1.020, 1.022, 1.022,
156  1.023, 1.030, 1.043, 1.044, 1.061, 1.067, 1.068, 1.070,
157  1.080, 1.080, 1.093, 1.100, 1.101, 1.102, 1.115, 1.117,
158  1.119, 1.130, 1.134, 1.144, 1.150, 1.153, 1.169, 1.169,
159  1.174, 1.179, 1.183, 1.194, 1.219, 1.220, 1.226, 1.230,
160  1.230, 1.244, 1.263, 1.270, 1.274, 1.283, 1.306, 1.316,
161  1.320, 1.325, 1.328, 1.336, 1.350, 1.368, 1.370, 1.375,
162  1.383, 1.384, 1.406, 1.415, 1.433, 1.434, 1.442, 1.462,
163  1.480, 1.483, 1.485, 1.488, 1.509, 1.514, 1.521, 1.534,
164  1.546, 1.549, 1.562, 1.583, 1.584, 1.600, 1.606, 1.614,
165  1.634, 1.637, 1.653, 1.663, 1.684, 1.687, 1.705, 1.710,
166  1.730, 1.734, 1.735, 1.741, 1.763, 1.784, 1.785, 1.800,
167  1.805, 1.835, 1.840, 1.843, 1.855, 1.879, 1.884, 1.885,
168  1.903, 1.924, 1.929, 1.934, 1.950, 1.953, 1.984, 1.985,
169  2.005, 2.031, 2.034, 2.035, 2.052, 2.080, 2.084, 2.104,
170  2.124, 2.134, 2.135, 2.159, 2.164, 2.175, 2.206, 2.207,
171  2.234, 2.240, 2.250, 2.257, 2.284, 2.301, 2.308, 2.325,
172  2.331, 2.350, 2.351, 2.356, 2.374, 2.392, 2.400, 2.404,
173  2.423, 2.450, 2.467, 2.480, 2.500, 2.514, 2.550, 2.558,
174  2.600, 2.608, 2.630, 2.650, 2.657, 2.660, 2.660, 2.670,
175  2.700, 2.700, 2.750, 2.800, 2.850, 2.880, 2.900, 2.950,
176  2.970, 2.980, 3.000, 3.000, 3.050, 3.090, 3.100, 3.150,
177  3.190, 3.200, 3.250, 3.300, 3.300, 3.500, 3.710, 3.920,
178  3.980, 4.000, 4.130, 4.200, 4.250, 4.340, 4.760, 5.180,
179  5.500, 6.000, 6.650, 7.000, 7.000, 8.000, 8.500, 10.000,
180  10.000, 10.000, 10.000, 10.000, 10.100, 11.500, 12.000, 12.000,
181  12.600, 13.300, 13.400, 14.000, 14.500, 15.000, 15.800, 16.000,
182  17.500, 18.000, 19.100, 20.000, 20.000, 20.000, 21.000, 21.000,
183  23.000, 25.000, 25.000, 25.000, 27.500, 30.000, 30.000, 30.000,
184  31.500, 32.000, 33.500, 35.000, 35.000, 35.500, 37.500, 40.000,
185  40.000, 40.000, 40.000, 42.500, 45.000, 45.000, 45.000, 47.500,
186  50.000, 50.000, 50.000, 50.000, 50.000, 52.500, 55.000, 55.000,
187  55.000, 70.000, 100.000, 100.000, 100.000, 120.000, 147.000, 150.000,
188  150.000, 170.000, 175.000, 200.000, 200.000, 240.000, 280.000, 310.000};
189 
191 static std::unique_ptr<InterpolateDataLinear<double>>
193 
195 const std::initializer_list<double> KMINUSP_TOT_SIG = {
196  113.80, 98.00, 94.00, 96.70, 75.10, 89.30, 90.70, 82.50, 79.40, 78.60,
197  70.90, 76.00, 71.50, 71.30, 70.20, 69.30, 68.80, 69.32, 70.00, 81.80,
198  100.19, 88.30, 69.60, 77.80, 80.60, 91.60, 76.20, 69.80, 76.60, 70.70,
199  58.61, 66.70, 61.50, 62.80, 54.90, 56.70, 58.41, 55.70, 50.27, 57.20,
200  48.60, 46.82, 45.30, 48.45, 44.19, 43.19, 42.30, 42.70, 44.71, 50.10,
201  49.50, 40.70, 40.11, 40.50, 39.55, 42.15, 38.34, 39.30, 37.78, 39.10,
202  35.60, 36.74, 36.32, 32.99, 37.80, 36.50, 36.03, 38.59, 35.07, 36.30,
203  32.40, 34.40, 34.75, 36.36, 35.40, 32.70, 34.41, 34.31, 34.44, 35.40,
204  34.90, 34.38, 33.70, 34.34, 35.33, 34.16, 34.81, 33.10, 34.63, 32.60,
205  35.37, 35.40, 34.69, 33.30, 35.26, 35.64, 36.80, 31.90, 37.65, 36.40,
206  38.84, 40.10, 38.10, 38.57, 38.50, 39.87, 40.60, 40.70, 40.31, 40.90,
207  40.96, 40.53, 40.20, 41.30, 40.67, 39.70, 40.29, 37.30, 39.86, 36.10,
208  40.58, 40.60, 38.10, 41.13, 41.75, 41.30, 42.36, 43.70, 42.70, 43.41,
209  43.41, 52.00, 42.10, 43.20, 44.86, 44.10, 42.40, 45.46, 46.24, 42.90,
210  47.33, 47.16, 43.20, 47.19, 47.96, 47.10, 45.30, 45.80, 49.79, 51.70,
211  51.27, 47.70, 52.19, 48.40, 52.17, 51.90, 50.30, 49.40, 51.46, 49.30,
212  50.51, 46.90, 46.13, 43.60, 44.50, 43.50, 44.18, 40.90, 42.23, 43.74,
213  39.30, 39.80, 45.00, 40.50, 39.18, 38.67, 42.00, 39.01, 40.70, 36.92,
214  35.21, 35.51, 33.40, 35.41, 33.80, 33.39, 32.36, 32.52, 32.11, 32.50,
215  31.17, 29.68, 31.03, 30.60, 30.50, 31.12, 31.20, 32.61, 31.27, 30.90,
216  31.20, 31.09, 31.63, 29.78, 31.90, 31.83, 32.28, 31.68, 32.50, 32.90,
217  33.35, 32.83, 33.47, 29.87, 33.53, 33.80, 32.16, 34.00, 33.97, 34.17,
218  34.20, 32.50, 33.19, 34.42, 34.20, 34.07, 32.79, 34.19, 33.80, 33.84,
219  33.20, 33.66, 32.50, 33.20, 33.19, 33.42, 32.77, 32.30, 32.26, 32.94,
220  31.92, 31.53, 31.19, 31.32, 30.72, 30.85, 30.40, 30.81, 30.18, 30.40,
221  30.49, 30.14, 30.50, 30.46, 30.10, 30.17, 30.37, 30.71, 29.90, 29.99,
222  29.86, 29.89, 29.80, 29.70, 29.81, 29.60, 31.77, 29.83, 29.79, 29.70,
223  29.87, 30.05, 30.75, 29.79, 29.79, 29.46, 29.50, 29.63, 29.37, 29.30,
224  31.89, 29.44, 29.38, 28.81, 28.90, 29.07, 29.10, 28.86, 28.88, 28.87,
225  28.50, 26.90, 28.46, 28.37, 28.30, 28.13, 28.21, 27.88, 30.30, 28.10,
226  27.96, 30.40, 30.20, 27.41, 28.02, 30.90, 27.83, 27.72, 27.73, 26.65,
227  27.58, 27.51, 25.30, 27.76, 24.90, 27.42, 27.20, 26.25, 27.00, 26.98,
228  26.75, 26.78, 26.70, 26.73, 27.00, 27.15, 26.90, 25.95, 25.40, 28.20,
229  25.14, 25.40, 25.80, 25.84, 24.74, 24.03, 24.30, 24.00, 23.80, 25.00,
230  23.50, 23.60, 24.60, 23.20, 22.90, 22.60, 22.67, 22.50, 22.50, 23.30,
231  22.90, 21.60, 21.60, 22.10, 21.50, 21.50, 21.60, 21.60, 21.90, 21.30,
232  21.90, 21.00, 21.60, 21.20, 21.02, 22.40, 21.16, 21.09, 20.89, 20.72,
233  20.70, 20.81, 20.72, 21.30, 20.84, 20.57, 20.50, 19.90, 20.55, 20.80,
234  20.54, 20.43, 20.39, 20.90, 20.59, 20.62, 20.49, 20.35, 20.60, 20.57,
235  20.45, 20.45, 21.00, 20.48, 20.30, 20.25, 20.40, 20.32, 21.50, 20.49,
236  20.23, 20.38, 20.80, 20.45, 20.41, 20.59, 21.50, 20.60, 20.57, 20.65,
237  21.40, 20.79, 20.84, 21.30, 21.32, 21.45};
238 
240 const std::initializer_list<double> KMINUSP_RES_SQRTS = {
241  1.4355, 1.4530, 1.4730, 1.4929, 1.5129, 1.5329, 1.5529, 1.5729,
242  1.5929, 1.6128, 1.6328, 1.6528, 1.6728, 1.6928, 1.7127, 1.7327,
243  1.7527, 1.7727, 1.7927, 1.8127, 1.8326, 1.8526, 1.8726, 1.8926,
244  1.9126, 1.9325, 1.9525, 1.9725, 1.9925, 2.0125, 2.0325, 2.0524,
245  2.0724, 2.0924, 2.1124, 2.1324, 2.1523, 2.1723, 2.1923, 2.2123,
246  2.2323, 2.2523, 2.2722, 2.2922, 2.3122, 2.3322, 2.3522, 2.3721,
247  2.3921, 2.4121, 2.4321, 2.4521, 2.4721, 2.4920, 2.5120, 2.5320};
248 
256 const std::initializer_list<double> KMINUSP_RES_SIG = {
257  0.46856081097, 0.68994120313, 1.00475205957, 1.66469547655,
258  5.63530052434, 4.17372028288, 4.50737036469, 8.00913400697,
259  0.29205365102, 2.72859364291, 3.30822314603, 4.44740017628,
260  4.95697831919, 5.05350905117, 4.87562017799, 7.48383422000,
261  8.29845755438, 9.71940157530, 11.10200040600, 12.00610574630,
262  10.06137989140, 7.47886042856, 6.11390219499, 5.19531126779,
263  4.38090191191, 3.95316327084, 3.53446044755, 3.46497827089,
264  3.63741875589, 3.77762079044, 0.87409952036, 4.19070149234,
265  4.38097308237, 4.27752586136, 4.12637945445, 3.70027602474,
266  3.31806303484, 2.88526838044, 2.58141493751, 2.36391939397,
267  2.18133708906, 1.39193162095, 2.03247269918, 2.00726146262,
268  2.13817978212, 2.16907178433, 2.08118209913, 1.83166338166,
269  1.56038155638, 1.27216056674, 1.03167072054, 0.85006416230,
270  0.39627220898, 0.57172926654, 0.51129452389, 0.44626386026};
271 
273 static std::unique_ptr<InterpolateDataSpline>
275 
281 // clang-format off
282 const std::initializer_list<double> KPLUSN_TOT_PLAB = {
283  0.770, 0.888, 0.939, 0.970, 0.989, 1.040, 1.091, 1.141,
284  /*1.170,*/ 1.191, 1.242, 1.292, 1.300, 1.342, 1.392, 1.440,
285  1.442, 1.492, 1.550, 1.593, 1.600, 1.643, 1.690, 1.693,
286  1.700, 1.743, 1.750, 1.793, 1.800, 1.850, 1.893, 1.900,
287  1.950, 1.970, 1.993, 2.000, 2.050, 2.093, 2.100, 2.150,
288  2.193, 2.200, 2.260, 2.300, 2.350, 2.393, 2.400, 2.450,
289  2.500, 2.550, 2.550, 2.600, 2.650, 2.700, 2.750, 2.800,
290  2.830, 2.850, 2.900, 2.950, 3.000, 3.050, 3.100, 3.150,
291  3.200, 3.250, 3.300, 6.000, 8.000, 10.000, 12.000, 14.000,
292  15.000, 16.000, 18.000, 20.000, 20.000, 25.000, 30.000, 35.000,
293  35.000, 40.000, 45.000, 50.000, 50.000, 50.000, 55.000, 70.000,
294  100.000, 100.000, 120.000, 150.000, 150.000, 170.000, 200.000, 200.000,
295  240.000, 280.000, 310.000};
296 // clang-format on
297 
303 const std::initializer_list<double> KPLUSN_TOT_SIG = {
304  15.50, 16.85, 17.60, 17.80, 18.53, 18.91, 20.61, 21.25, /*18.20,*/ 20.87,
305  20.26, 19.68, 18.50, 19.32, 19.22, 18.10, 19.07, 18.95, 18.91,
306  18.79, 18.89, 18.67, 18.50, 18.69, 18.83, 18.88, 18.86, 18.73,
307  18.53, 18.66, 18.50, 18.69, 18.70, 18.60, 18.55, 18.79, 18.54,
308  18.67, 18.49, 18.43, 18.40, 18.40, 17.70, 18.27, 18.26, 18.63,
309  18.09, 18.25, 18.11, 17.10, 18.17, 18.09, 18.02, 18.11, 18.06,
310  18.01, 17.50, 17.95, 17.85, 17.81, 17.81, 17.83, 17.85, 17.61,
311  17.61, 17.66, 17.55, 17.50, 17.60, 17.50, 17.60, 17.50, 17.87,
312  17.40, 17.60, 17.94, 17.70, 17.78, 17.69, 18.29, 18.12, 18.15,
313  18.30, 18.66, 18.56, 18.02, 18.43, 18.60, 19.04, 18.99, 19.23,
314  19.63, 19.55, 19.74, 19.72, 19.82, 20.37, 20.61, 20.80};
315 
317 static std::unique_ptr<InterpolateDataLinear<double>>
319 
321 const std::initializer_list<double> KPLUSP_TOT_PLAB = {
322  0.178, 0.265, 0.321, 0.351, 0.366, 0.405, 0.440, 0.451,
323  0.475, 0.475, 0.506, 0.522, 0.536, 0.566, 0.569, 0.588,
324  0.589, 0.592, 0.593, 0.596, 0.618, 0.620, 0.627, 0.643,
325  0.644, 0.657, 0.668, 0.686, 0.698, 0.713, 0.717, 0.727,
326  0.757, 0.768, 0.770, 0.786, 0.817, 0.823, 0.846, 0.864,
327  0.864, 0.881, 0.891, 0.900, 0.904, 0.916, 0.938, 0.942,
328  0.951, 0.969, 0.969, 0.970, 0.970, 0.985, 0.992, 1.020,
329  1.029, 1.043, 1.055, 1.060, 1.084, 1.090, 1.094, 1.125,
330  1.130, 1.140, 1.144, 1.160, 1.170, 1.189, 1.194, 1.207,
331  1.210, 1.238, 1.245, 1.250, 1.293, 1.295, 1.300, 1.320,
332  1.345, 1.347, 1.380, 1.395, 1.408, 1.440, 1.445, 1.455,
333  1.468, 1.480, 1.495, 1.550, 1.563, 1.596, 1.600, 1.646,
334  1.690, 1.696, 1.700, 1.746, 1.750, 1.796, 1.800, 1.850,
335  1.896, 1.900, 1.945, 1.950, 1.960, 1.970, 1.996, 2.000,
336  2.050, 2.096, 2.100, 2.150, 2.196, 2.200, 2.260, 2.300,
337  2.350, 2.396, 2.400, 2.450, 2.473, 2.500, 2.530, 2.550,
338  2.550, 2.600, 2.650, 2.700, 2.750, 2.760, 2.800, 2.830,
339  2.850, 2.900, 2.950, 3.000, 3.050, 3.100, 3.150, 3.200,
340  3.200, 3.250, 3.250, 3.300, 3.700, 4.000, 4.200, 4.750,
341  5.000, 5.500, 6.000, 7.000, 7.000, 8.000, 8.200, 8.500,
342  10.000, 10.000, 10.000, 10.000, 10.900, 11.500, 12.000, 12.500,
343  13.400, 14.000, 15.000, 15.000, 16.000, 16.000, 16.900, 18.000,
344  19.000, 20.000, 20.000, 25.000, 30.000, 32.000, 35.000, 35.000,
345  40.000, 42.500, 45.000, 50.000, 50.000, 50.000, 52.200, 55.000,
346  70.000, 100.000, 100.000, 100.000, 100.000, 120.000, 147.000, 150.000,
347  150.000, 170.000, 175.000, 200.000, 200.000, 240.000, 280.000, 310.000};
348 
350 const std::initializer_list<double> KPLUSP_TOT_SIG = {
351  11.40, 13.00, 14.00, 12.20, 13.20, 13.69, 12.81, 16.30, 12.70, 13.58, 13.02,
352  15.20, 12.09, 13.20, 12.70, 12.60, 16.30, 14.36, 13.05, 13.04, 12.65, 12.91,
353  12.18, 12.50, 12.88, 12.43, 13.10, 11.25, 12.60, 11.14, 11.10, 12.45, 12.65,
354  11.65, 13.00, 12.80, 13.20, 12.97, 13.45, 14.07, 13.21, 13.90, 14.39, 13.10,
355  14.23, 14.20, 14.59, 15.57, 14.95, 15.28, 15.63, 15.40, 15.25, 16.20, 15.97,
356  16.10, 15.69, 17.39, 16.95, 16.40, 17.04, 17.60, 17.12, 17.55, 18.08, 18.02,
357  18.09, 17.95, 18.10, 18.06, 18.47, 19.85, 18.58, 18.11, 18.54, 20.71, 18.44,
358  18.61, 17.90, 19.33, 18.44, 18.27, 18.64, 18.27, 17.97, 18.10, 18.04, 18.20,
359  17.94, 18.04, 17.93, 17.70, 17.66, 17.75, 17.71, 17.86, 17.50, 17.85, 17.73,
360  17.80, 17.83, 17.80, 17.98, 17.77, 17.81, 17.79, 17.41, 17.75, 19.40, 16.90,
361  17.60, 17.63, 17.72, 17.51, 17.56, 17.57, 17.54, 17.60, 17.10, 17.44, 17.52,
362  17.55, 17.56, 17.48, 17.25, 17.49, 17.47, 17.10, 17.44, 17.50, 17.47, 17.41,
363  17.41, 17.41, 17.40, 16.70, 17.30, 17.34, 17.30, 17.19, 17.14, 17.08, 17.15,
364  17.13, 17.13, 17.13, 17.50, 17.14, 21.00, 17.60, 17.10, 21.30, 17.20, 17.90,
365  17.00, 17.20, 18.40, 17.30, 17.20, 18.70, 17.20, 17.30, 18.80, 17.30, 18.10,
366  19.00, 17.30, 18.30, 17.50, 17.40, 18.50, 17.31, 17.10, 17.00, 18.80, 17.10,
367  17.30, 17.50, 17.42, 17.68, 17.72, 18.40, 17.82, 17.80, 18.05, 17.91, 17.88,
368  18.06, 18.03, 18.37, 18.28, 18.17, 18.52, 18.40, 18.88, 18.70, 18.85, 19.14,
369  19.52, 19.36, 19.33, 19.64, 18.20, 19.91, 19.84, 20.22, 20.45, 20.67};
370 
372 static std::unique_ptr<InterpolateDataLinear<double>>
374 
376 const std::initializer_list<double> PIMINUSP_ELASTIC_P_LAB = {
377  0.09875, 0.14956, 0.21648, 0.21885, 0.22828, 0.24684, 0.25599, 0.26733,
378  0.27071, 0.2752, 0.29303, 0.29303, 0.33138, 0.33571, 0.33788, 0.35052,
379  0.378, 0.38261, 0.404, 0.40626, 0.408, 0.42188, 0.427, 0.44888,
380  0.452, 0.471, 0.49008, 0.49008, 0.49008, 0.509, 0.523, 0.52845,
381  0.53155, 0.547, 0.54911, 0.54911, 0.556, 0.565, 0.57281, 0.582,
382  0.586, 0.609, 0.6139, 0.61698, 0.625, 0.64054, 0.657, 0.65793,
383  0.658, 0.6753, 0.683, 0.687, 0.69061, 0.699, 0.70692, 0.71399,
384  0.72628, 0.731, 0.73257, 0.74257, 0.75, 0.76189, 0.775, 0.77714,
385  0.77714, 0.77827, 0.798, 0.82586, 0.83803, 0.848, 0.84954, 0.854,
386  0.87466, 0.90386, 0.91903, 0.924, 0.95947, 1.03016, 1.106, 1.12339,
387  1.164, 1.165, 1.174, 1.214, 1.21659, 1.23, 1.25, 1.26,
388  1.279, 1.323, 1.339, 1.347, 1.355, 1.365, 1.375, 1.385,
389  1.395, 1.405, 1.415, 1.425, 1.435, 1.445, 1.455, 1.465,
390  1.475, 1.485, 1.495, 1.497, 1.499, 1.503, 1.5031, 1.509,
391  1.567, 1.59, 1.6, 1.603, 1.71, 1.85, 2.1, 2.14,
392  2.26, 2.29, 2.7, 2.75, 2.77, 2.7999, 3., 3.15};
393 
395 const std::initializer_list<double> PIMINUSP_ELASTIC_SIG = {
396  1.847, 2.9, 9.6, 11.3, 12.8, 17., 20.1, 21.4, 22.5,
397  21.2, 22.5, 18.2, 19.5, 16., 17.4, 15.1, 12.29, 12.4,
398  10.1, 13.8, 10.41, 11.4, 9., 10.3, 8.9, 9.2, 10.42,
399  10.9, 10.8, 9.15, 9.8, 11.4, 11.4, 9.99, 13., 13.,
400  10.24, 10.8, 12.19, 16.2, 11.34, 12.86, 13.71, 13.9, 12.19,
401  14.8, 13.92, 16.2, 15.32, 16.98, 18.9, 17.07, 18.86, 19.07,
402  19.95, 20.5, 19.87, 18.9, 16.6, 19.4, 19.91, 18.94, 17.56,
403  17.19, 17.82, 16.1, 14.91, 15.75, 14.9, 13.2, 14.1, 14.47,
404  14.4, 14.8, 14.1, 18.8, 18., 18.6, 17.95, 17.7, 13.66,
405  15.01, 15.73, 12.45, 14.1, 14.6, 13.31, 13.8, 12.8, 13.09,
406  12.627, 12.987, 12.763, 12.367, 12.852, 12.67, 12.126, 12.972, 12.478,
407  12.594, 12.532, 11.801, 11.568, 11.413, 11.119, 11.643, 11.368, 11.523,
408  11.163, 11.69, 10., 10.39, 10.21, 9.65, 9., 9.82, 10.4,
409  11.1, 9.69, 9.3, 8.91, 8.5, 7.7, 7.2, 7.2, 7.8,
410  7.57, 6.1};
411 
413 static std::unique_ptr<InterpolateDataLinear<double>>
415 
417 const std::initializer_list<double> PIMINUSP_LAMBDAK0_P_LAB = {
418  0.904, 0.91, 0.919, 0.922, 0.926, 0.93, 0.931, 0.942, 0.945, 0.958, 0.964,
419  0.968, 0.98, 0.98, 0.983, 0.993, 0.997, 1.001, 1.002, 1.007, 1.012, 1.017,
420  1.02, 1.02, 1.021, 1.022, 1.023, 1.027, 1.031, 1.035, 1.04, 1.04, 1.04,
421  1.043, 1.048, 1.054, 1.059, 1.065, 1.078, 1.081, 1.091, 1.091, 1.094, 1.095,
422  1.097, 1.116, 1.129, 1.13, 1.135, 1.144, 1.159, 1.194, 1.223, 1.235, 1.277,
423  1.287, 1.326, 1.331, 1.332, 1.395, 1.433, 1.455, 1.5, 1.508, 1.515, 1.575,
424  1.59, 1.6, 1.615, 1.645, 1.69, 1.69, 1.705, 1.775, 1.85, 1.875, 1.94,
425  1.95, 1.98, 2., 2.05, 2.05, 2.055, 2.115, 2.14, 2.15, 2.235, 2.25,
426  2.35, 2.375, 2.494, 2.605, 2.7, 2.75, 2.75, 2.86, 3.01, 3.125, 3.21,
427  3.885, 3.9, 3.95, 4.16, 4.5, 6., 8., 10.};
428 
430 const std::initializer_list<double> PIMINUSP_LAMBDAK0_SIG = {
431  0.056, 0.122, 0.18, 0.14, 0.227, 0.212, 0.13, 0.3, 0.336, 0.43,
432  0.427, 0.52, 0.467, 0.45, 0.576, 0.59, 0.652, 0.56, 0.588, 0.634,
433  0.686, 0.665, 0.67, 0.69, 0.809, 0.675, 0.94, 0.737, 0.734, 0.73,
434  0.926, 0.59, 0.92, 0.57, 0.568, 0.651, 0.899, 0.64, 0.794, 0.58,
435  0.82, 0.58, 0.7, 0.68, 0.729, 0.575, 0.592, 0.462, 0.541, 0.64,
436  0.48, 0.58, 0.46, 0.485, 0.447, 0.25, 0.367, 0.32, 0.29, 0.25,
437  0.32, 0.29, 0.334, 0.214, 0.22, 0.21, 0.214, 0.238, 0.208, 0.16,
438  0.199, 0.174, 0.14, 0.13, 0.181, 0.16, 0.185, 0.182, 0.184, 0.15,
439  0.182, 0.179, 0.11, 0.16, 0.162, 0.192, 0.15, 0.172, 0.174, 0.12,
440  0.16, 0.106, 0.12, 0.09, 0.09, 0.109, 0.084, 0.094, 0.087, 0.067,
441  0.058, 0.0644, 0.049, 0.054, 0.038, 0.0221, 0.0157};
442 
444 static std::unique_ptr<InterpolateDataLinear<double>>
446 
448 const std::initializer_list<double> PIMINUSP_SIGMAMINUSKPLUS_P_LAB = {
449  1.091, 1.128, 1.17, 1.22, 1.235, 1.284, 1.326, 1.5, 1.59,
450  1.615, 1.69, 1.69, 1.85, 1.94, 1.95, 1.98, 2.05, 2.14,
451  2.15, 2.25, 2.35, 2.494, 2.61, 2.7, 2.75, 2.86, 3.,
452  3.01, 3.13, 3.21, 3.89, 3.95, 4., 4.16};
453 
455 const std::initializer_list<double> PIMINUSP_SIGMAMINUSKPLUS_SIG = {
456  0.25, 0.218, 0.231, 0.27, 0.235, 0.209, 0.245, 0.242, 0.262,
457  0.18, 0.153, 0.19, 0.099, 0.098, 0.099, 0.09, 0.087, 0.069,
458  0.065, 0.057, 0.053, 0.051, 0.03, 0.031, 0.032, 0.022, 0.015,
459  0.022, 0.0155, 0.0145, 0.0085, 0.0096, 0.005, 0.0045};
460 
465 static std::unique_ptr<InterpolateDataLinear<double>>
467 
469 const std::initializer_list<double> PIMINUSP_SIGMA0K0_RES_SQRTS = {
470  1.5, 1.516, 1.532, 1.548, 1.564, 1.58, 1.596, 1.612, 1.628, 1.644, 1.66,
471  1.676, 1.692, 1.708, 1.724, 1.74, 1.756, 1.772, 1.788, 1.804, 1.82, 1.836,
472  1.852, 1.868, 1.884, 1.9, 1.916, 1.932, 1.948, 1.964, 1.98, 1.996, 2.012,
473  2.028, 2.044, 2.06, 2.076, 2.092, 2.108, 2.124, 2.14, 2.156, 2.172, 2.188,
474  2.204, 2.22, 2.236, 2.252, 2.268, 2.284, 2.3, 2.316, 2.332, 2.348, 2.364,
475  2.38, 2.396, 2.412, 2.428, 2.444, 2.46, 2.476, 2.492, 2.508, 2.524, 2.54,
476  2.556, 2.572, 2.588, 2.604, 2.62, 2.636, 2.652, 2.668, 2.684, 2.7, 2.716,
477  2.732, 2.748, 2.764, 2.78, 2.796, 2.812, 2.828, 2.844, 2.86, 2.876, 2.892,
478  2.908, 2.924, 2.94, 2.956, 2.972, 2.988, 3.004, 3.02, 3.036, 3.052, 3.068,
479  3.084, 3.1, 3.116, 3.132, 3.148, 3.164, 3.18};
480 
487 const std::initializer_list<double> PIMINUSP_SIGMA0K0_RES_SIG = {
488  0., 0., 0., 0., 0., 0.,
489  0., 0., 0., 0., 0., 0.,
490  0.0386981, 0.09589789, 0.11956695, 0.11685363, 0.12053117, 0.13208736,
491  0.14949223, 0.16688579, 0.18830654, 0.20611132, 0.22231072, 0.23099061,
492  0.23734563, 0.23334048, 0.22794051, 0.21559531, 0.20134617, 0.18763246,
493  0.1723282, 0.15814744, 0.14757816, 0.13750278, 0.12698656, 0.11719809,
494  0.11024985, 0.1044732, 0.09623321, 0.09092108, 0.08670191, 0.08147493,
495  0.0772165, 0.07346243, 0.0719974, 0.06805902, 0.06496733, 0.06264939,
496  0.05904799, 0.05762721, 0.05588871, 0.05393479, 0.0517673, 0.05165839,
497  0.05087591, 0.04885535, 0.04730724, 0.04651682, 0.04604065, 0.04529776,
498  0.04406593, 0.04367817, 0.04230014, 0.04144308, 0.04171145, 0.04073006,
499  0.03996921, 0.03902337, 0.03939531, 0.03895125, 0.03904553, 0.03816119,
500  0.03772662, 0.03710955, 0.0361001, 0.03632378, 0.03549849, 0.03549987,
501  0.03527251, 0.034535, 0.03314715, 0.0335742, 0.03326698, 0.0330181,
502  0.0324203, 0.03227253, 0.0315376, 0.03065083, 0.03041305, 0.03023753,
503  0.03008669, 0.02900321, 0.02827017, 0.02805024, 0.02785525, 0.02753706,
504  0.02692862, 0.02603758, 0.02591122, 0.02537291, 0.02467199, 0.02466657,
505  0.02370074, 0.02353027, 0.02362089, 0.0230085};
506 
511 static std::unique_ptr<InterpolateDataLinear<double>>
513 
515 const std::initializer_list<double> PIMINUSP_RES_SQRTS = {
516  1.1438620, 1.1482410, 1.1514750, 1.1566800, 1.1572040, 1.1579910, 1.1665900,
517  1.1679230, 1.1701290, 1.1720720, 1.1764440, 1.1765450, 1.1773210, 1.1786720,
518  1.1823280, 1.1845210, 1.1869010, 1.1881750, 1.1884820, 1.1895390, 1.1931590,
519  1.1955810, 1.1984320, 1.1998050, 1.2002920, 1.2018580, 1.2047590, 1.2057700,
520  1.2077800, 1.2087440, 1.2101220, 1.2117760, 1.2119900, 1.2150870, 1.2151980,
521  1.2156390, 1.2158600, 1.2178480, 1.2190910, 1.2203690, 1.2204730, 1.2204940,
522  1.2228010, 1.2255380, 1.2260010, 1.2271080, 1.2274050, 1.2291260, 1.2294580,
523  1.2308410, 1.2338980, 1.2342780, 1.2350390, 1.2363740, 1.2374810, 1.2398320,
524  1.2411660, 1.2412150, 1.2447480, 1.2455710, 1.2456470, 1.2464000, 1.2483350,
525  1.2501600, 1.2508920, 1.2515830, 1.2516660, 1.2519140, 1.2539170, 1.2554980,
526  1.2569070, 1.2579900, 1.2584040, 1.2603780, 1.2605570, 1.2613850, 1.2619360,
527  1.2628740, 1.2633150, 1.2653620, 1.2658450, 1.2666570, 1.2677530, 1.2688130,
528  1.2701420, 1.2708920, 1.2722680, 1.2726330, 1.2750260, 1.2774170, 1.2799030,
529  1.2801710, 1.2812000, 1.2820710, 1.2823940, 1.2825720, 1.2860820, 1.2885750,
530  1.2897730, 1.2914760, 1.2934590, 1.2971470, 1.3008370, 1.3017160, 1.3045220,
531  1.3052980, 1.3053930, 1.3082150, 1.3112370, 1.3119160, 1.3119980, 1.3156860,
532  1.3173860, 1.3194780, 1.3196130, 1.3213980, 1.3231680, 1.3247750, 1.3268590,
533  1.3302420, 1.3302480, 1.3315140, 1.3344390, 1.3382320, 1.3407860, 1.3415830,
534  1.3420180, 1.3477680, 1.3509370, 1.3516040, 1.3530900, 1.3547200, 1.3609120,
535  1.3616290, 1.3636960, 1.3648880, 1.3675350, 1.3675540, 1.3685060, 1.3753530,
536  1.3767680, 1.3780830, 1.3889330, 1.3911860, 1.3920020, 1.3956730, 1.3973910,
537  1.3990370, 1.3997060, 1.4003760, 1.4009540, 1.4041980, 1.4090600, 1.4106690,
538  1.4157070, 1.4171170, 1.4177600, 1.4179530, 1.4190210, 1.4223220, 1.4252400,
539  1.4275900, 1.4289080, 1.4299430, 1.4387270, 1.4420430, 1.4442120, 1.4471120,
540  1.4484790, 1.4534290, 1.4564910, 1.4577870, 1.4581640, 1.4600930, 1.4647410,
541  1.4653680, 1.4677830, 1.4698700, 1.4734890, 1.4759800, 1.4773430, 1.4790890,
542  1.4805180, 1.4846720, 1.4852910, 1.4864480, 1.4868440, 1.4877650, 1.4889450,
543  1.4896180, 1.4912540, 1.4914700, 1.4918830, 1.4931410, 1.4952340, 1.4962800,
544  1.4970120, 1.4988550, 1.5025360, 1.5039210, 1.5049860, 1.5056590, 1.5062770,
545  1.5068210, 1.5070170, 1.5098750, 1.5120090, 1.5143540, 1.5147500, 1.5180690,
546  1.5194050, 1.5201020, 1.5202180, 1.5208240, 1.5217760, 1.5220360, 1.5303840,
547  1.5317880, 1.5323060, 1.5353150, 1.5359820, 1.5383300, 1.5407180, 1.5413170,
548  1.5425990, 1.5431130, 1.5439750, 1.5493680, 1.5502810, 1.5508760, 1.5516930,
549  1.5560580, 1.5586050, 1.5607390, 1.5633430, 1.5651170, 1.5668890, 1.5680450,
550  1.5704270, 1.5757210, 1.5798270, 1.5817360, 1.5827540, 1.5840110, 1.5862590,
551  1.5912250, 1.5944120, 1.5970180, 1.5973140, 1.5990530, 1.5994180, 1.6016820,
552  1.6048380, 1.6088750, 1.6111040, 1.6133580, 1.6152010, 1.6168420, 1.6186420,
553  1.6232190, 1.6249900, 1.6261880, 1.6295560, 1.6323380, 1.6368790, 1.6374460,
554  1.6376790, 1.6398760, 1.6404590, 1.6421720, 1.6456730, 1.6495730, 1.6515580,
555  1.6518560, 1.6518680, 1.6541320, 1.6549290, 1.6603630, 1.6631810, 1.6644430,
556  1.6655580, 1.6676540, 1.6688130, 1.6689020, 1.6693530, 1.6699370, 1.6716500,
557  1.6722390, 1.6794480, 1.6800020, 1.6800180, 1.6800570, 1.6822520, 1.6828430,
558  1.6860770, 1.6932300, 1.6933350, 1.6955350, 1.6961330, 1.6977940, 1.6987140,
559  1.7044260, 1.7069070, 1.7093030, 1.7115550, 1.7123480, 1.7126470, 1.7159200,
560  1.7199390, 1.7204800, 1.7235710, 1.7252330, 1.7264200, 1.7268130, 1.7273740,
561  1.7339520, 1.7344240, 1.7403830, 1.7473240, 1.7478990, 1.7479100, 1.7500400,
562  1.7531770, 1.7605990, 1.7611280, 1.7664090, 1.7691810, 1.7737780, 1.7764030,
563  1.7766390, 1.7795470, 1.7797780, 1.7807510, 1.7829640, 1.7868640, 1.7925930,
564  1.7998580, 1.8007640, 1.8055470, 1.8127620, 1.8178990, 1.8193350, 1.8209750,
565  1.8214510, 1.8214820, 1.8255780, 1.8286820, 1.8367850, 1.8383080, 1.8414410,
566  1.8469150, 1.8509520, 1.8554840, 1.8584990, 1.8635130, 1.8710100, 1.8752210,
567  1.8759920, 1.8774990, 1.8787820, 1.8804650, 1.8854230, 1.8883910, 1.8911140,
568  1.8972700, 1.9007120, 1.9061080, 1.9129540, 1.9183170, 1.9251210, 1.9270610,
569  1.9280640, 1.9303930, 1.9372130, 1.9400900, 1.9449130, 1.9492320, 1.9497830,
570  1.9611780, 1.9645110, 1.9688770, 1.9730540, 1.9756480, 1.9796740, 1.9878660,
571  1.9881530, 2.0115230, 2.0159570, 2.0332500, 2.0519420, 2.0579350, 2.0633810,
572  2.0796300, 2.0796340, 2.0930840, 2.1011130, 2.1056180, 2.1219370, 2.1238140,
573  2.1373130, 2.1443060, 2.1587020, 2.1655840, 2.1802270, 2.1875130, 2.2024370,
574  2.2099850, 2.2303110, 2.2507840, 2.2686000, 2.2710750, 2.3119360, 2.3312970,
575  2.3536850, 2.3734980, 2.3742870, 2.3923640, 2.4103070, 2.4126380, 2.4300470,
576  2.4511570, 2.4625960, 2.4890840, 2.5264450, 2.5632640, 2.5923800, 2.5995640,
577  2.6103560, 2.6353660, 2.6706900, 2.6986180, 2.7055550, 2.7399780, 2.7739760,
578  2.8075630, 2.8407540, 2.8670310, 2.8768230, 2.9092270, 2.9316960, 2.9412740,
579  2.9729770, 3.0043460, 3.0353910, 3.0569360, 3.0661230, 3.0965510, 3.0988120,
580  3.1266830, 3.1565270, 3.1772520, 3.1860930, 3.1919740, 3.2153880, 3.2444180,
581  3.2731920, 3.3017150, 3.3328080, 3.3608260, 3.3886130, 3.4161730, 3.4189170,
582  3.4435140, 3.4706390, 3.4868130, 3.4975540, 3.5507730, 3.6032070, 3.6548900,
583  3.6574550, 3.7058520, 3.7310720, 3.7461230, 3.7955240, 3.9886240, 4.0003630,
584  4.2614590, 4.4340360, 4.6809920, 4.7288360, 4.8386340};
590 const std::initializer_list<double> PIMINUSP_RES_SIG = {
591  2.727771, 3.571329, 4.034613, 4.882370, 4.975210, 5.117328, 6.891753,
592  7.205112, 7.746671, 8.247897, 9.456023, 9.485245, 9.711213, 10.112155,
593  11.241098, 11.944095, 12.723382, 13.144982, 13.246803, 13.598222, 14.798669,
594  15.586762, 16.480923, 16.893047, 17.035961, 17.481312, 18.242040, 18.484793,
595  18.928599, 19.122060, 19.375409, 19.642156, 19.673564, 20.048232, 20.058785,
596  20.099053, 20.118021, 20.253818, 20.307050, 20.336885, 20.338216, 20.338463,
597  20.326376, 20.215612, 20.187220, 20.108831, 20.085299, 19.929646, 19.896012,
598  19.744140, 19.348620, 19.294422, 19.183131, 18.979071, 18.802318, 18.407836,
599  18.174742, 18.166186, 17.527526, 17.376048, 17.362015, 17.222750, 16.864140,
600  16.526686, 16.391817, 16.265062, 16.249887, 16.204412, 15.841202, 15.559769,
601  15.313297, 15.126607, 15.055993, 14.725313, 14.695755, 14.560464, 14.471319,
602  14.321746, 14.252233, 13.937248, 13.864858, 13.744437, 13.585398, 13.434918,
603  13.251312, 13.150076, 12.968950, 12.921955, 12.624439, 12.345808, 12.074788,
604  12.046720, 11.940788, 11.853593, 11.821900, 11.804500, 11.481213, 11.273034,
605  11.179088, 11.052115, 10.913921, 10.686029, 10.491401, 10.449671, 10.328092,
606  10.297759, 10.294122, 10.194918, 10.106365, 10.088903, 10.086866, 10.009062,
607  9.982275, 9.955883, 9.954419, 9.937770, 9.926040, 9.919670, 9.917476,
608  9.925773, 9.925804, 9.932451, 9.954887, 9.999073, 10.035233, 10.047444,
609  10.054284, 10.154801, 10.215371, 10.228372, 10.257565, 10.289514, 10.404053,
610  10.416507, 10.451255, 10.469744, 10.506548, 10.506807, 10.518891, 10.583211,
611  10.588553, 10.591846, 10.551600, 10.528871, 10.519808, 10.468089, 10.437574,
612  10.407237, 10.394638, 10.381921, 10.370850, 10.307754, 10.204285, 10.171181,
613  10.069177, 10.040061, 10.027313, 10.023558, 10.003399, 9.948213, 9.908627,
614  9.883310, 9.872680, 9.866237, 9.877111, 9.919510, 9.961662, 10.035629,
615  10.078320, 10.277399, 10.432873, 10.507124, 10.529744, 10.652615, 10.997871,
616  11.049584, 11.260341, 11.456495, 11.825114, 12.094558, 12.246925, 12.447007,
617  12.613749, 13.107598, 13.181143, 13.361513, 13.418153, 13.534148, 13.667426,
618  13.738379, 13.898388, 13.918290, 13.955613, 14.063025, 14.218190, 14.284287,
619  14.326422, 14.417247, 14.531076, 14.550199, 14.554711, 14.553183, 14.549119,
620  14.543450, 14.540937, 14.476394, 14.395865, 14.277562, 14.253650, 14.025681,
621  13.921779, 13.865112, 13.855605, 13.804923, 13.723295, 13.700529, 12.899844,
622  12.759289, 12.707480, 12.406185, 12.339343, 12.108105, 11.880783, 11.825126,
623  11.708153, 11.661998, 11.585849, 11.139385, 11.070268, 11.026165, 10.967131,
624  10.676650, 10.527975, 10.416330, 10.296090, 10.224045, 10.160037, 10.122722,
625  10.056953, 9.963186, 9.940076, 9.944406, 9.950619, 9.962028, 9.992798,
626  10.107907, 10.216027, 10.324590, 10.338056, 10.421938, 10.440565, 10.563991,
627  10.758763, 11.046275, 11.219098, 11.406161, 11.568358, 11.719190, 11.891153,
628  12.356378, 12.545937, 12.676903, 13.054393, 13.374187, 13.902515, 13.968362,
629  13.995323, 14.248969, 14.315704, 14.510032, 14.894839, 15.294828, 15.483193,
630  15.510420, 15.511442, 15.708221, 15.771604, 16.129790, 16.258504, 16.302184,
631  16.332919, 16.371389, 16.381668, 16.382132, 16.383765, 16.384106, 16.373504,
632  16.365846, 16.094248, 16.061361, 16.060351, 16.057988, 15.912247, 15.869115,
633  15.604058, 14.881278, 14.869656, 14.615127, 14.544006, 14.343831, 14.231480,
634  13.518719, 13.205142, 12.904069, 12.624201, 12.524907, 12.487428, 12.080104,
635  11.596000, 11.532586, 11.165458, 10.970463, 10.834819, 10.790496, 10.727909,
636  10.033831, 9.986213, 9.422328, 8.843849, 8.799297, 8.798477, 8.637944,
637  8.413399, 7.933449, 7.901865, 7.604069, 7.459614, 7.237359, 7.119512,
638  7.109214, 6.986207, 6.976764, 6.937338, 6.850621, 6.707226, 6.516829,
639  6.307029, 6.283124, 6.165071, 6.009385, 5.914360, 5.889920, 5.863003,
640  5.855379, 5.854890, 5.792956, 5.750120, 5.653524, 5.637576, 5.606512,
641  5.557455, 5.524890, 5.491104, 5.470122, 5.436638, 5.387926, 5.360343,
642  5.355253, 5.345255, 5.336686, 5.325224, 5.289797, 5.267817, 5.247120,
643  5.194926, 5.164390, 5.113899, 5.046572, 4.990838, 4.916877, 4.895353,
644  4.884239, 4.858425, 4.779720, 4.745950, 4.688505, 4.637428, 4.630918,
645  4.491197, 4.449922, 4.395158, 4.342866, 4.310165, 4.260304, 4.156622,
646  4.153074, 3.863210, 3.807705, 3.599068, 3.381739, 3.313336, 3.252640,
647  3.078632, 3.078586, 2.942304, 2.865732, 2.822595, 2.680616, 2.664734,
648  2.559076, 2.508390, 2.409106, 2.363398, 2.270298, 2.225212, 2.133668,
649  2.087464, 1.964694, 1.842465, 1.737800, 1.723454, 1.497293, 1.398065,
650  1.291078, 1.202957, 1.199580, 1.125166, 1.056454, 1.047924, 0.986448,
651  0.917631, 0.882818, 0.808313, 0.716179, 0.638085, 0.583949, 0.571516,
652  0.553453, 0.514354, 0.465001, 0.430238, 0.422140, 0.384707, 0.351896,
653  0.323011, 0.297502, 0.279222, 0.272812, 0.253013, 0.240484, 0.235418,
654  0.219741, 0.205756, 0.193259, 0.185301, 0.182074, 0.172035, 0.171327,
655  0.163000, 0.154840, 0.149582, 0.147433, 0.146033, 0.140669, 0.134455,
656  0.128697, 0.123321, 0.117770, 0.112988, 0.108418, 0.104023, 0.103592,
657  0.099782, 0.095677, 0.093278, 0.091704, 0.084136, 0.077092, 0.070601,
658  0.070291, 0.064685, 0.061942, 0.060365, 0.055497, 0.040625, 0.039905,
659  0.027723, 0.022456, 0.017122, 0.016299, 0.014606};
660 
662 static std::unique_ptr<InterpolateDataSpline>
664 
666 const std::initializer_list<double> PIPLUSP_ELASTIC_P_LAB = {
667  0.09875, 0.13984, 0.14956, 0.33138, 0.378, 0.408, 0.4093, 0.427,
668  0.42736, 0.453, 0.471, 0.49008, 0.509, 0.5214, 0.53155, 0.547,
669  0.57281, 0.574, 0.586, 0.5891, 0.59, 0.591, 0.6139, 0.625,
670  0.625, 0.635, 0.645, 0.645, 0.657, 0.65793, 0.6753, 0.687,
671  0.698, 0.7, 0.70692, 0.725, 0.72628, 0.72628, 0.752, 0.752,
672  0.77714, 0.807, 0.809, 0.82, 0.82586, 0.85, 0.895, 0.895,
673  0.9, 0.91, 0.93926, 0.945, 0.945, 0.995, 0.995, 1.0029,
674  1.02, 1.04, 1.04, 1.0402, 1.05, 1.12091, 1.12091, 1.195,
675  1.207, 1.2217, 1.232, 1.375, 1.3925, 1.44384, 1.46, 1.493,
676  1.5, 1.53, 1.585, 1.585, 1.6, 1.68881, 1.69, 1.77,
677  1.8, 1.869, 1.9, 1.99, 2., 2.0199, 2.07, 2.077,
678  2.11, 2.19, 2.3, 2.35, 2.5, 2.77, 2.9, 3.,
679  3.05, 3.56, 3.63, 3.65, 3.67, 3.9, 4., 5.,
680  5., 6., 6.8001, 8., 8.04, 8.8, 10.8, 11.7,
681  12.8, 14.8, 16., 16.2, 16.7, 29., 32.1, 43.,
682  50., 60., 70., 100., 140., 147., 175., 200.,
683  250.};
684 
686 const std::initializer_list<double> PIPLUSP_ELASTIC_SIG = {
687  6.15, 15.8, 20.4, 140.9, 91.6, 71.6, 67.5, 57.3, 60.19, 46.6, 40.8,
688  38.74, 30.6, 29.6, 30.59, 24.74, 24.31, 28.16, 19.83, 20.4, 20.64, 20.63,
689  19.55, 15.75, 18.5, 17.2, 14.85, 16.16, 14.71, 15.32, 14.38, 12.2, 12.93,
690  12.96, 12.17, 11.6, 11.06, 11.5, 10.62, 10.55, 8.82, 9.36, 8.97, 9.1,
691  8.02, 8.38, 8.37, 8.14, 11.1, 9.3, 11., 10.32, 9.87, 11.7, 11.15,
692  12.05, 11.1, 12.37, 11.83, 10.3, 12.8, 14.54, 15.3, 14.3, 13.59, 12.6,
693  12.3, 17.87, 16.5, 19.31, 18.73, 16.68, 13.8, 15.86, 16.7, 16.7, 15.05,
694  13.04, 13.57, 12.46, 12.3, 11.81, 10.9, 9.84, 10.6, 9.1, 9.52, 9.46,
695  9.44, 9.15, 8.45, 10.2, 6.9, 7.7, 8.3, 7.84, 7.02, 6.93, 7.07,
696  6.88, 7.15, 6.5, 6.4, 5.85, 5.79, 5.33, 5.47, 4.9, 4.9, 5.02,
697  4.75, 4.2, 4.54, 4.46, 4.21, 4.21, 3.98, 3.19, 3.37, 3.16, 3.29,
698  3.1, 3.35, 3.3, 3.39, 3.24, 3.37, 3.17, 3.3};
699 
701 static std::unique_ptr<InterpolateDataLinear<double>>
703 
705 const std::initializer_list<double> PIPLUSP_SIGMAPLUSKPLUS_P_LAB = {
706  1.041, 1.105, 1.111, 1.15, 1.157, 1.17, 1.195, 1.206, 1.218, 1.222, 1.265,
707  1.28, 1.282, 1.328, 1.34, 1.377, 1.39, 1.41, 1.419, 1.43, 1.456, 1.49,
708  1.508, 1.518, 1.549, 1.55, 1.58, 1.582, 1.614, 1.63, 1.68, 1.687, 1.7,
709  1.712, 1.76, 1.77, 1.775, 1.808, 1.84, 1.879, 1.906, 1.95, 1.971, 1.997,
710  2.067, 2.08, 2.099, 2.152, 2.197, 2.241, 2.291, 2.344, 2.379, 2.437, 2.473,
711  2.77, 3.23, 3.71, 4., 5., 5.5, 7., 10.3, 12., 16.};
712 
714 const std::initializer_list<double> PIPLUSP_SIGMAPLUSKPLUS_SIG = {
715  0.034, 0.146, 0.144, 0.214, 0.248, 0.205, 0.24, 0.214, 0.242, 0.25,
716  0.278, 0.34, 0.369, 0.412, 0.4, 0.467, 0.44, 0.49, 0.523, 0.51,
717  0.662, 0.529, 0.692, 0.545, 0.604, 0.53, 0.53, 0.465, 0.494, 0.47,
718  0.505, 0.434, 0.47, 0.436, 0.38, 0.415, 0.418, 0.393, 0.405, 0.332,
719  0.37, 0.31, 0.338, 0.298, 0.3, 0.29, 0.28, 0.273, 0.26, 0.25,
720  0.23, 0.242, 0.22, 0.217, 0.234, 0.165, 0.168, 0.104, 0.059, 0.059,
721  0.0297, 0.0371, 0.02, 0.0202, 0.0143};
722 
727 static std::unique_ptr<InterpolateDataLinear<double>>
729 
731 const std::initializer_list<double> PIPLUSP_RES_SQRTS = {
732  1.1173610, 1.1241380, 1.1358180, 1.1371030, 1.1380990, 1.1424360, 1.1457360,
733  1.1465540, 1.1482410, 1.1566800, 1.1579910, 1.1583200, 1.1604390, 1.1665900,
734  1.1679230, 1.1685180, 1.1749480, 1.1757030, 1.1768490, 1.1773210, 1.1786720,
735  1.1829250, 1.1881750, 1.1884820, 1.1892660, 1.1892730, 1.1895390, 1.1913820,
736  1.1924340, 1.1940070, 1.1947940, 1.1955810, 1.1984320, 1.1998050, 1.2002920,
737  1.2010810, 1.2018580, 1.2048200, 1.2049860, 1.2087440, 1.2090200, 1.2096600,
738  1.2101220, 1.2112110, 1.2119900, 1.2127620, 1.2143280, 1.2156390, 1.2158600,
739  1.2166330, 1.2189530, 1.2190910, 1.2204730, 1.2212600, 1.2231810, 1.2251030,
740  1.2255180, 1.2260010, 1.2262570, 1.2266370, 1.2274050, 1.2294580, 1.2300810,
741  1.2308410, 1.2342780, 1.2350390, 1.2358000, 1.2363740, 1.2398320, 1.2412150,
742  1.2426250, 1.2428740, 1.2456470, 1.2508920, 1.2515830, 1.2531650, 1.2602810,
743  1.2605570, 1.2606400, 1.2619360, 1.2621290, 1.2688130, 1.2708920, 1.2722680,
744  1.2743520, 1.2792100, 1.2812000, 1.2825720, 1.2908540, 1.2914760, 1.2921600,
745  1.3000800, 1.3017160, 1.3112370, 1.3119160, 1.3174950, 1.3196130, 1.3213980,
746  1.3247750, 1.3252880, 1.3315140, 1.3382320, 1.3415830, 1.3442600, 1.3463790,
747  1.3516040, 1.3602490, 1.3609120, 1.3616290, 1.3648880, 1.3657620, 1.3698540,
748  1.3761100, 1.3780830, 1.3808030, 1.3889330, 1.3911860, 1.3943310, 1.3944480,
749  1.4041980, 1.4077270, 1.4106690, 1.4157070, 1.4171170, 1.4210000, 1.4282490,
750  1.4293040, 1.4299430, 1.4328100, 1.4341490, 1.4419790, 1.4420430, 1.4465240,
751  1.4490100, 1.4534290, 1.4536500, 1.4553190, 1.4647410, 1.4672450, 1.4759800,
752  1.4824160, 1.4855750, 1.4865280, 1.4877650, 1.4912540, 1.4956530, 1.4988550,
753  1.4994080, 1.5025360, 1.5093620, 1.5098750, 1.5118750, 1.5141410, 1.5208240,
754  1.5279280, 1.5321730, 1.5323060, 1.5347140, 1.5425990, 1.5431130, 1.5462770,
755  1.5508760, 1.5510910, 1.5526630, 1.5564860, 1.5633430, 1.5651170, 1.5668890,
756  1.5702910, 1.5715220, 1.5757210, 1.5786540, 1.5862590, 1.5961540, 1.5973140,
757  1.6058020, 1.6129040, 1.6139380, 1.6186080, 1.6284110, 1.6323380, 1.6404590,
758  1.6416010, 1.6467170, 1.6498710, 1.6518170, 1.6529970, 1.6581190, 1.6626220,
759  1.6666740, 1.6671250, 1.6688130, 1.6722390, 1.6822130, 1.6850290, 1.6961330,
760  1.6981660, 1.7011210, 1.7175570, 1.7204800, 1.7235710, 1.7392590, 1.7425210,
761  1.7500400, 1.7607570, 1.7632420, 1.7713600, 1.7766390, 1.7790240, 1.7903550,
762  1.7924370, 1.7946710, 1.8028500, 1.8096730, 1.8194380, 1.8235340, 1.8286820,
763  1.8301710, 1.8434770, 1.8494390, 1.8536470, 1.8642640, 1.8700130, 1.8737520,
764  1.8787820, 1.8824500, 1.8908620, 1.8912220, 1.8936740, 1.8977620, 1.9035830,
765  1.9158620, 1.9182680, 1.9183170, 1.9280640, 1.9328680, 1.9381770, 1.9421740,
766  1.9517610, 1.9568860, 1.9613690, 1.9706840, 1.9756480, 1.9777850, 1.9801460,
767  1.9893750, 1.9983050, 2.0031410, 2.0085470, 2.0341700, 2.0359620, 2.0383030,
768  2.0561620, 2.0633810, 2.0810910, 2.0815680, 2.0930840, 2.1104420, 2.1256810,
769  2.1276610, 2.1386260, 2.1520030, 2.1603920, 2.1736890, 2.1767910, 2.1802270,
770  2.1815150, 2.1917870, 2.2020120, 2.2058340, 2.2164220, 2.2478710, 2.2625000,
771  2.2932290, 2.3035850, 2.3337060, 2.3479530, 2.3532870, 2.3536850, 2.3734980,
772  2.3875670, 2.4126380, 2.4265400, 2.4511570, 2.4686940, 2.4890840, 2.5064020,
773  2.5264450, 2.5435540, 2.5449200, 2.5632640, 2.5801700, 2.5995640, 2.6162720,
774  2.6353660, 2.6554170, 2.6706900, 2.7055550, 2.7251500, 2.7399780, 2.7604270,
775  2.7739760, 2.7965000, 2.8075630, 2.8407540, 2.8660730, 2.8768230, 2.8995430,
776  2.9092270, 2.9339980, 2.9412740, 2.9729770, 3.0003830, 3.0043460, 3.0353910,
777  3.0569360, 3.0653340, 3.0661230, 3.0965510, 3.1266830, 3.1565270, 3.1860930,
778  3.2153880, 3.2444180, 3.2731920, 3.3017150, 3.3328080, 3.3608260, 3.3886130,
779  3.4161730, 3.4189170, 3.4435140, 3.4706390, 3.4868130, 3.4975540, 3.5242640,
780  3.5507730, 3.6032070, 3.6548900, 3.7058520, 3.7461230, 3.9246250, 3.9886240,
781  4.4000620, 4.4340360, 4.8192100, 4.8386340};
787 const std::initializer_list<double> PIPLUSP_RES_SIG = {
788  6.022200, 10.038871, 17.441789, 18.453186, 19.268833, 23.150742,
789  26.498913, 27.386442, 29.291426, 40.496896, 42.513388, 43.030735,
790  46.492651, 57.799559, 60.507365, 61.744072, 76.313703, 78.161300,
791  81.021590, 82.218538, 85.698091, 97.146898, 112.028626, 112.912122,
792  115.173370, 115.193049, 115.960670, 121.274197, 124.294133, 128.769360,
793  130.983774, 133.177973, 140.888488, 144.431558, 145.658198, 147.606924,
794  149.472691, 156.080287, 156.422132, 163.346112, 163.784724, 164.765987,
795  165.439046, 166.911553, 167.863958, 168.723549, 170.204028, 171.170636,
796  171.308879, 171.736988, 172.504248, 172.525712, 172.593592, 172.514719,
797  171.976396, 170.970613, 170.695374, 170.349193, 170.155496, 169.854117,
798  169.198023, 167.145300, 166.443383, 165.538602, 160.877328, 159.734532,
799  158.557081, 157.647233, 151.838017, 149.388278, 146.835348, 146.380094,
800  141.238822, 131.363425, 130.066638, 127.112163, 114.248942, 113.768689,
801  113.624927, 111.391996, 111.062599, 100.192520, 97.031881, 94.998359,
802  92.007883, 85.447008, 82.920313, 81.231779, 71.911700, 71.268389,
803  70.570322, 63.112285, 61.706466, 54.325885, 53.847853, 50.135114,
804  48.820524, 47.750605, 45.816616, 45.532806, 42.281080, 39.137415,
805  37.695516, 36.598844, 35.763674, 33.820390, 30.931186, 30.724707,
806  30.504071, 29.529065, 29.275637, 28.130677, 26.504830, 26.021034,
807  25.375549, 23.583128, 23.119865, 22.496127, 22.473352, 20.695879,
808  20.104852, 19.631763, 18.860276, 18.652732, 18.099067, 17.131947,
809  16.997986, 16.917638, 16.564509, 16.403533, 15.511183, 15.504274,
810  15.029072, 14.775835, 14.342938, 14.321885, 14.164381, 13.330021,
811  13.123245, 12.447103, 11.992319, 11.781862, 11.720067, 11.640951,
812  11.424298, 11.165000, 10.986072, 10.955990, 10.790078, 10.454817,
813  10.431082, 10.340373, 10.241359, 9.973671, 9.727943, 9.600744,
814  9.597007, 9.531720, 9.352521, 9.342579, 9.286710, 9.220119,
815  9.217412, 9.198613, 9.160760, 9.116223, 9.108690, 9.102419,
816  9.092648, 9.089575, 9.079161, 9.070168, 9.030578, 8.920964,
817  8.903042, 8.745744, 8.587201, 8.563189, 8.451960, 8.222839,
818  8.135503, 7.964442, 7.941527, 7.840761, 7.779875, 7.743029,
819  7.721027, 7.623950, 7.539104, 7.460968, 7.452091, 7.418960,
820  7.352368, 7.151795, 7.095082, 6.872018, 6.831788, 6.775597,
821  6.499647, 6.458988, 6.421107, 6.281186, 6.267541, 6.253040,
822  6.283141, 6.299494, 6.373666, 6.441742, 6.477662, 6.688097,
823  6.733830, 6.785977, 6.995237, 7.190716, 7.500324, 7.636831,
824  7.812619, 7.864044, 8.315853, 8.503866, 8.626044, 8.879732,
825  8.974573, 9.017728, 9.051747, 9.056952, 9.008060, 9.004429,
826  8.972971, 8.906696, 8.781147, 8.419200, 8.335117, 8.333399,
827  7.962610, 7.768317, 7.542001, 7.372803, 6.958122, 6.731881,
828  6.543614, 6.150650, 5.943860, 5.858059, 5.766199, 5.407669,
829  5.077684, 4.907944, 4.721214, 3.920965, 3.868515, 3.802953,
830  3.320499, 3.139608, 2.718237, 2.708015, 2.448100, 2.100843,
831  1.843201, 1.814261, 1.669739, 1.527892, 1.453866, 1.356321,
832  1.336277, 1.315603, 1.308211, 1.252123, 1.203159, 1.186153,
833  1.142877, 1.035775, 0.993647, 0.915748, 0.892059, 0.828894,
834  0.801562, 0.791697, 0.790969, 0.756008, 0.732619, 0.693611,
835  0.673343, 0.639635, 0.617218, 0.592709, 0.573136, 0.551832,
836  0.534734, 0.533411, 0.516219, 0.501298, 0.485222, 0.472219,
837  0.458270, 0.444610, 0.434844, 0.414467, 0.404081, 0.396683,
838  0.387082, 0.381073, 0.371641, 0.367241, 0.354831, 0.346043,
839  0.342460, 0.335135, 0.332105, 0.324567, 0.322406, 0.313236,
840  0.305601, 0.304518, 0.296208, 0.290622, 0.288484, 0.288284,
841  0.280729, 0.273519, 0.266620, 0.259982, 0.253547, 0.247248,
842  0.241018, 0.234792, 0.227884, 0.221501, 0.214989, 0.208334,
843  0.207660, 0.201536, 0.194611, 0.190406, 0.187587, 0.180499,
844  0.173394, 0.159321, 0.145738, 0.132952, 0.123434, 0.088815,
845  0.079356, 0.042881, 0.041067, 0.026625, 0.026107};
846 
848 static std::unique_ptr<InterpolateDataSpline>
850 
852 const std::initializer_list<double> PIPLUSP_TOT_SQRTS = {
853  1.0825000, 1.0925000, 1.1050000, 1.1175000, 1.1300000, 1.1425000, 1.1550000,
854  1.1675000, 1.1800000, 1.1925000, 1.2050000, 1.2175000, 1.2300000, 1.2425000,
855  1.2550000, 1.2675000, 1.2800000, 1.2925000, 1.3050000, 1.3175000, 1.3300000,
856  1.3425000, 1.3550000, 1.3675000, 1.3800000, 1.3925000, 1.4050000, 1.4175000,
857  1.4300000, 1.4425000, 1.4550000, 1.4675000, 1.4800000, 1.4925000, 1.5050000,
858  1.5175000, 1.5300000, 1.5425000, 1.5550000, 1.5675000, 1.5800000, 1.5925000,
859  1.6050000, 1.6175000, 1.6300000, 1.6425000, 1.6550000, 1.6675000, 1.6800000,
860  1.6925000, 1.7050000, 1.7175000, 1.7300000, 1.7425000, 1.7550000, 1.7675000,
861  1.7800000, 1.7925000, 1.8050000, 1.8175000, 1.8300000, 1.8425000, 1.8550000,
862  1.8675000, 1.8800000, 1.8925000, 1.9050000, 1.9175000, 1.9300000, 1.9425000,
863  1.9550000, 1.9675000, 1.9800000, 1.9925000, 2.0050000, 2.0175000, 2.0300000,
864  2.0425000, 2.0550000, 2.0675000, 2.0800000, 2.0925000, 2.1050000, 2.1175000,
865  2.1300000, 2.1425000, 2.1550000, 2.1675000, 2.1800000, 2.1925000, 2.2050000,
866  2.2175000, 2.2300000, 2.2425000, 2.2550000, 2.2675000, 2.2800000, 2.2925000,
867  2.3050000, 2.3175000, 2.3300000, 2.3425000, 2.3550000, 2.3675000, 2.3800000,
868  2.3925000, 2.4050000, 2.4175000, 2.4300000, 2.4425000, 2.4550000, 2.4675000,
869  2.4800000, 2.4925000, 2.5050000, 2.5175000, 2.5300000, 2.5425000, 2.5550000,
870  2.5675000, 2.5800000, 2.5925000, 2.6050000, 2.6175000, 2.6300000, 2.6425000,
871  2.6550000, 2.6675000, 2.6800000, 2.6925000, 2.7050000, 2.7175000, 2.7300000,
872  2.7425000, 2.7550000, 2.7675000, 2.7800000, 2.7925000, 2.8050000, 2.8175000,
873  2.8300000, 2.8425000, 2.8550000, 2.8675000, 2.8800000, 2.8925000, 2.9050000,
874  2.9175000, 2.9300000, 2.9425000, 2.9550000, 2.9675000, 2.9800000, 2.9925000,
875  3.0050000, 3.0175000, 3.0300000, 3.0425000, 3.0550000, 3.0675000, 3.0800000,
876  3.0925000, 3.1050000, 3.1175000, 3.1300000, 3.1425000, 3.1550000, 3.1675000,
877  3.1800000, 3.1925000, 3.2050000, 3.2175000, 3.2300000, 3.2425000, 3.2550000,
878  3.2675000, 3.2800000, 3.2925000, 3.3050000, 3.3175000, 3.3300000, 3.3425000,
879  3.3550000, 3.3675000, 3.3800000, 3.3925000, 3.4050000, 3.4175000, 3.4300000,
880  3.4425000, 3.4550000, 3.4675000, 3.4800000, 3.4925000, 3.5050000, 3.5175000,
881  3.5300000, 3.5425000, 3.5550000, 3.5675000, 3.5800000, 3.5925000, 3.6050000,
882  3.6175000, 3.6300000, 3.6425000, 3.6550000, 3.6675000, 3.6800000, 3.6925000,
883  3.7050000, 3.7175000, 3.7300000, 3.7425000, 3.7550000, 3.7675000, 3.7800000,
884  3.7925000, 3.8050000, 3.8175000, 3.8300000, 3.8425000, 3.8550000, 3.8675000,
885  3.8800000, 3.8925000, 3.9050000, 3.9175000, 3.9300000, 3.9425000, 3.9550000,
886  3.9675000, 3.9800000, 3.9925000, 4.0050000, 4.0175000, 4.0300000, 4.0425000,
887  4.0550000, 4.0675000, 4.0800000, 4.0925000, 4.1050000, 4.1175000, 4.1300000,
888  4.1425000, 4.1550000, 4.1675000, 4.1800000, 4.1925000, 4.2050000, 4.2175000,
889  4.2300000, 4.2425000, 4.2550000, 4.2675000, 4.2800000, 4.2925000, 4.3050000,
890  4.3175000, 4.3300000, 4.3425000, 4.3550000, 4.3675000, 4.3800000, 4.3925000,
891  4.4050000, 4.4175000, 4.4300000, 4.4425000, 4.4550000, 4.4675000, 4.4800000,
892  4.4925000, 4.5050000, 4.5175000, 4.5300000, 4.5425000, 4.5550000, 4.5675000,
893  4.5800000, 4.5925000, 4.6050000, 4.6175000, 4.6300000, 4.6425000, 4.6550000,
894  4.6675000, 4.6800000, 4.6925000, 4.7050000, 4.7175000, 4.7300000, 4.7425000,
895  4.7550000, 4.7675000, 4.7800000, 4.7925000, 4.8050000, 4.8175000, 4.8300000,
896  4.8425000, 4.8550000, 4.8675000, 4.8800000, 4.8925000, 4.9050000, 4.9175000,
897  4.9300000, 4.9425000, 4.9550000, 4.9675000, 4.9800000, 4.9925000, 5.0050000,
898  5.0175000, 5.0300000, 5.0425000, 5.0550000, 5.0675000, 5.0800000, 5.0925000,
899  5.1050000, 5.1175000, 5.1300000, 5.1425000, 5.1550000, 5.1675000, 5.1800000,
900  5.1925000, 5.2050000, 5.2175000, 5.2300000, 5.2425000, 5.2550000, 5.2675000,
901  5.2800000, 5.2925000, 5.3050000, 5.3175000, 5.3300000, 5.3425000, 5.3550000,
902  5.3675000, 5.3800000, 5.3925000, 5.4050000, 5.4175000, 5.4300000, 5.4425000,
903  5.4550000, 5.4675000, 5.4800000, 5.4925000, 5.5050000, 5.5175000, 5.5300000,
904  5.5425000, 5.5550000, 5.5675000, 5.5800000, 5.5925000, 5.6050000, 5.6175000,
905  5.6300000, 5.6425000, 5.6550000, 5.6675000, 5.6800000, 5.6925000, 5.7050000,
906  5.7175000, 5.7300000, 5.7425000, 5.7550000, 5.7675000, 5.7800000, 5.7925000,
907  5.8050000, 5.8175000, 5.8300000, 5.8425000, 5.8550000, 5.8675000, 5.8800000,
908  5.8925000, 5.9050000, 5.9175000, 5.9300000, 5.9425000, 5.9550000, 5.9675000,
909  5.9800000, 5.9925000, 6.0050000, 6.0175000, 6.0300000, 6.0425000, 6.0550000,
910  6.0675000, 6.0800000};
911 
914 const std::initializer_list<double> PIPLUSP_TOT_SIG = {
915  0.091281, 0.962113, 3.137700, 7.044444, 13.412202, 23.289333,
916  38.047868, 59.689265, 89.229993, 124.462341, 156.483158, 172.253319,
917  166.612043, 147.129923, 123.805319, 102.380895, 84.566587, 70.376965,
918  59.179948, 50.333752, 43.276005, 37.513856, 32.877710, 29.103911,
919  25.939398, 23.253257, 21.072625, 19.157278, 17.610095, 16.260927,
920  15.138353, 14.257853, 13.538350, 13.031891, 14.313882, 14.073960,
921  13.809499, 13.576393, 14.170392, 15.239011, 16.399310, 17.349004,
922  18.366929, 19.247360, 20.163402, 21.016895, 21.832534, 22.535146,
923  22.831721, 22.892679, 23.162585, 23.586223, 24.278299, 24.982553,
924  25.949646, 27.128727, 28.602172, 30.467640, 32.685130, 35.073372,
925  37.416683, 39.678407, 41.743503, 43.400009, 43.966275, 43.854866,
926  43.074254, 41.741346, 40.042463, 38.062695, 35.889092, 33.816300,
927  31.944858, 30.317878, 28.952918, 27.868056, 27.147320, 26.684076,
928  26.446123, 26.415282, 26.518282, 26.827488, 27.222833, 27.727514,
929  28.203844, 28.749745, 29.174557, 29.547540, 29.783788, 29.897626,
930  29.864765, 29.807407, 29.752597, 29.669430, 29.646424, 29.551999,
931  29.538956, 29.454628, 29.350158, 29.310488, 29.255296, 29.241534,
932  29.154373, 29.149326, 29.060770, 29.018977, 28.952918, 28.921213,
933  28.837072, 28.802318, 28.770130, 28.702180, 28.652110, 28.635464,
934  28.572797, 28.524768, 28.474279, 28.462450, 28.416912, 28.352874,
935  28.306302, 28.270080, 28.248553, 28.182611, 28.117140, 28.114662,
936  28.055208, 28.008394, 27.956379, 27.933125, 27.906964, 27.847823,
937  27.798540, 27.783755, 27.738047, 27.692756, 27.658220, 27.611654,
938  27.603484, 27.529951, 27.476844, 27.464588, 27.427920, 27.401728,
939  27.405919, 27.342361, 27.317915, 27.261654, 27.284158, 27.208205,
940  27.171635, 27.147320, 27.161214, 27.102165, 27.060074, 27.025747,
941  26.972840, 26.969788, 26.928750, 26.897369, 26.868160, 26.851232,
942  26.802864, 26.791929, 26.744823, 26.735150, 26.709613, 26.677266,
943  26.662795, 26.608529, 26.590440, 26.537961, 26.545341, 26.501630,
944  26.479427, 26.446123, 26.446123, 26.423981, 26.374161, 26.366965,
945  26.315283, 26.309774, 26.291158, 26.268427, 26.191790, 26.244786,
946  26.179172, 26.168236, 26.140775, 26.126859, 26.110955, 26.063242,
947  26.056426, 26.028541, 26.015529, 26.015529, 25.977463, 25.916705,
948  25.944156, 25.905724, 25.893175, 25.866004, 25.810503, 25.872783,
949  25.785789, 25.754155, 25.764323, 25.748223, 25.730429, 25.704618,
950  25.703812, 25.646923, 25.623959, 25.632149, 25.618044, 25.588469,
951  25.581018, 25.564874, 25.534019, 25.524119, 25.512966, 25.484652,
952  25.456338, 25.464596, 25.446900, 25.415528, 25.402783, 25.360023,
953  25.342877, 25.364548, 25.314549, 25.335977, 25.305725, 25.279328,
954  25.277569, 25.294896, 25.262407, 25.227513, 25.239039, 25.191339,
955  25.200138, 25.164943, 25.153243, 25.138619, 25.125945, 25.131249,
956  25.087728, 25.098877, 25.058440, 25.044608, 25.066670, 25.024553,
957  25.014553, 24.987220, 25.002448, 24.962332, 24.947554, 24.927163,
958  24.948191, 24.884555, 24.891924, 24.899292, 24.884555, 24.867105,
959  24.827071, 24.820123, 24.844668, 24.824725, 24.814753, 24.786011,
960  24.783025, 24.744951, 24.755690, 24.744951, 24.731693, 24.718849,
961  24.707828, 24.661424, 24.655458, 24.654873, 24.652143, 24.633581,
962  24.640542, 24.625626, 24.605739, 24.605739, 24.605739, 24.579710,
963  24.580499, 24.536329, 24.545006, 24.520905, 24.524081, 24.518064,
964  24.522448, 24.526414, 24.490056, 24.466920, 24.466920, 24.466920,
965  24.450948, 24.411549, 24.449617, 24.392383, 24.377932, 24.383864,
966  24.383864, 24.393093, 24.359255, 24.328494, 24.344779, 24.328494,
967  24.328494, 24.328494, 24.278299, 24.300887, 24.293985, 24.297819,
968  24.259477, 24.259477, 24.236471, 24.231870, 24.270374, 24.205797,
969  24.236471, 24.198128, 24.190460, 24.190460, 24.183216, 24.190460,
970  24.190460, 24.184475, 24.135403, 24.113993, 24.180628, 24.121639,
971  24.160965, 24.098699, 24.095832, 24.087229, 24.098699, 24.087229,
972  24.073994, 24.052819, 24.052819, 24.052819, 24.052819, 24.039094,
973  24.029944, 24.027865, 24.010588, 23.984195, 24.007069, 23.970470,
974  23.995632, 23.952170, 23.955938, 23.970470, 23.957801, 23.972757,
975  23.944981, 23.925374, 23.915570, 23.935177, 23.933870, 23.915570,
976  23.905043, 23.897323, 23.889910, 23.881357, 23.892761, 23.888199,
977  23.858547, 23.888199, 23.853985, 23.866693, 23.851168, 23.847143,
978  23.847143, 23.828480, 23.824333, 23.827017, 23.806086, 23.810297,
979  23.794815, 23.799770, 23.790119, 23.778715, 23.778715, 23.747223,
980  23.766309, 23.759220, 23.741498, 23.778715, 23.747223, 23.751422,
981  23.757168, 23.726229, 23.700736, 23.714497, 23.733227};
982 
984 static std::unique_ptr<InterpolateDataLinear<double>>
986 
988 const std::initializer_list<double> PIMINUSP_TOT_SQRTS = {
989  1.0825000, 1.0883300, 1.0966700, 1.1050000, 1.1133300, 1.1216700, 1.1300000,
990  1.1383300, 1.1466700, 1.1550000, 1.1633300, 1.1716700, 1.1800000, 1.1883300,
991  1.1966700, 1.2050000, 1.2133300, 1.2216700, 1.2300000, 1.2383300, 1.2466700,
992  1.2550000, 1.2633300, 1.2716700, 1.2800000, 1.2883300, 1.2966700, 1.3050000,
993  1.3133300, 1.3216700, 1.3300000, 1.3383300, 1.3466700, 1.3550000, 1.3633300,
994  1.3716700, 1.3800000, 1.3883300, 1.3966700, 1.4050000, 1.4133300, 1.4216700,
995  1.4300000, 1.4383300, 1.4466700, 1.4550000, 1.4633300, 1.4716700, 1.4800000,
996  1.4883300, 1.4966700, 1.5050000, 1.5133300, 1.5216700, 1.5300000, 1.5383300,
997  1.5466700, 1.5550000, 1.5633300, 1.5716700, 1.5800000, 1.5883300, 1.5966700,
998  1.6050000, 1.6133300, 1.6216700, 1.6300000, 1.6383300, 1.6466700, 1.6550000,
999  1.6633300, 1.6716700, 1.6800000, 1.6883300, 1.6966700, 1.7050000, 1.7133300,
1000  1.7216700, 1.7300000, 1.7383300, 1.7466700, 1.7550000, 1.7633300, 1.7716700,
1001  1.7800000, 1.7883300, 1.7966700, 1.8050000, 1.8133300, 1.8216700, 1.8300000,
1002  1.8383300, 1.8466700, 1.8550000, 1.8633300, 1.8716700, 1.8800000, 1.8883300,
1003  1.8966700, 1.9050000, 1.9133300, 1.9216700, 1.9300000, 1.9383300, 1.9466700,
1004  1.9550000, 1.9633300, 1.9716700, 1.9800000, 1.9883300, 1.9966700, 2.0050000,
1005  2.0133300, 2.0216700, 2.0300000, 2.0383300, 2.0466700, 2.0550000, 2.0633300,
1006  2.0716700, 2.0800000, 2.0883300, 2.0966700, 2.1050000, 2.1133300, 2.1216700,
1007  2.1300000, 2.1383300, 2.1466700, 2.1550000, 2.1633300, 2.1716700, 2.1800000,
1008  2.1883300, 2.1966700, 2.2050000, 2.2133300, 2.2216700, 2.2300000, 2.2383300,
1009  2.2466700, 2.2550000, 2.2633300, 2.2716700, 2.2800000, 2.2883300, 2.2966700,
1010  2.3050000, 2.3133300, 2.3216700, 2.3300000, 2.3383300, 2.3466700, 2.3550000,
1011  2.3633300, 2.3716700, 2.3800000, 2.3883300, 2.3966700, 2.4050000, 2.4133300,
1012  2.4216700, 2.4300000, 2.4383300, 2.4466700, 2.4550000, 2.4633300, 2.4716700,
1013  2.4800000, 2.4883300, 2.4966700, 2.5050000, 2.5133300, 2.5216700, 2.5300000,
1014  2.5383300, 2.5466700, 2.5550000, 2.5633300, 2.5716700, 2.5800000, 2.5883300,
1015  2.5966700, 2.6050000, 2.6133300, 2.6216700, 2.6300000, 2.6383300, 2.6466700,
1016  2.6550000, 2.6633300, 2.6716700, 2.6800000, 2.6883300, 2.6966700, 2.7050000,
1017  2.7133300, 2.7216700, 2.7300000, 2.7383300, 2.7466700, 2.7550000, 2.7633300,
1018  2.7716700, 2.7800000, 2.7883300, 2.7966700, 2.8050000, 2.8133300, 2.8216700,
1019  2.8300000, 2.8383300, 2.8466700, 2.8550000, 2.8633300, 2.8716700, 2.8800000,
1020  2.8883300, 2.8966700, 2.9050000, 2.9133300, 2.9216700, 2.9300000, 2.9383300,
1021  2.9466700, 2.9550000, 2.9633300, 2.9716700, 2.9800000, 2.9883300, 2.9966700,
1022  3.0050000, 3.0133300, 3.0216700, 3.0300000, 3.0383300, 3.0466700, 3.0550000,
1023  3.0633300, 3.0716700, 3.0800000, 3.0883300, 3.0966700, 3.1050000, 3.1133300,
1024  3.1216700, 3.1300000, 3.1383300, 3.1466700, 3.1550000, 3.1633300, 3.1716700,
1025  3.1800000, 3.1883300, 3.1966700, 3.2050000, 3.2133300, 3.2216700, 3.2300000,
1026  3.2383300, 3.2466700, 3.2550000, 3.2633300, 3.2716700, 3.2800000, 3.2883300,
1027  3.2966700, 3.3050000, 3.3133300, 3.3216700, 3.3300000, 3.3383300, 3.3466700,
1028  3.3550000, 3.3633300, 3.3716700, 3.3800000, 3.3883300, 3.3966700, 3.4050000,
1029  3.4133300, 3.4216700, 3.4300000, 3.4383300, 3.4466700, 3.4550000, 3.4633300,
1030  3.4716700, 3.4800000, 3.4883300, 3.4966700, 3.5050000, 3.5133300, 3.5216700,
1031  3.5300000, 3.5383300, 3.5466700, 3.5550000, 3.5633300, 3.5716700, 3.5800000,
1032  3.5883300, 3.5966700, 3.6050000, 3.6133300, 3.6216700, 3.6300000, 3.6383300,
1033  3.6466700, 3.6550000, 3.6633300, 3.6716700, 3.6800000, 3.6883300, 3.6966700,
1034  3.7050000, 3.7133300, 3.7216700, 3.7300000, 3.7383300, 3.7466700, 3.7550000,
1035  3.7633300, 3.7716700, 3.7800000, 3.7883300, 3.7966700, 3.8050000, 3.8133300,
1036  3.8216700, 3.8300000, 3.8383300, 3.8466700, 3.8550000, 3.8633300, 3.8716700,
1037  3.8800000, 3.8883300, 3.8966700, 3.9050000, 3.9133300, 3.9216700, 3.9300000,
1038  3.9383300, 3.9466700, 3.9550000, 3.9633300, 3.9716700, 3.9800000, 3.9883300,
1039  3.9966700, 4.0050000, 4.0133300, 4.0216700, 4.0300000, 4.0383300, 4.0466700,
1040  4.0550000, 4.0633300, 4.0716700, 4.0800000, 4.0883300, 4.0966700, 4.1050000,
1041  4.1133300, 4.1216700, 4.1300000, 4.1383300, 4.1466700, 4.1550000, 4.1633300,
1042  4.1716700, 4.1800000, 4.1883300, 4.1966700, 4.2050000, 4.2133300, 4.2216700,
1043  4.2300000, 4.2383300, 4.2466700, 4.2550000, 4.2633300, 4.2716700, 4.2800000,
1044  4.2883300, 4.2966700, 4.3050000, 4.3133300, 4.3216700, 4.3300000, 4.3383300,
1045  4.3466700, 4.3550000, 4.3633300, 4.3716700, 4.3800000, 4.3883300, 4.3966700,
1046  4.4050000, 4.4133300, 4.4216700, 4.4300000, 4.4383300, 4.4466700, 4.4550000,
1047  4.4633300, 4.4716700, 4.4800000, 4.4883300, 4.4966700, 4.5050000, 4.5133300,
1048  4.5216700, 4.5300000, 4.5383300, 4.5466700, 4.5550000, 4.5633300, 4.5716700,
1049  4.5800000, 4.5883300, 4.5966700, 4.6050000, 4.6133300, 4.6216700, 4.6300000,
1050  4.6383300, 4.6466700, 4.6550000, 4.6633300, 4.6716700, 4.6800000, 4.6883300,
1051  4.6966700, 4.7050000, 4.7133300, 4.7216700, 4.7300000, 4.7383300, 4.7466700,
1052  4.7550000, 4.7633300, 4.7716700, 4.7800000, 4.7883300, 4.7966700, 4.8050000,
1053  4.8133300, 4.8216700, 4.8300000, 4.8383300, 4.8466700, 4.8550000, 4.8633300,
1054  4.8716700, 4.8800000, 4.8883300, 4.8966700, 4.9050000, 4.9133300, 4.9216700,
1055  4.9300000, 4.9383300, 4.9466700, 4.9550000, 4.9633300, 4.9716700, 4.9800000,
1056  4.9883300, 4.9966700, 5.0050000, 5.0133300, 5.0216700, 5.0300000, 5.0383300,
1057  5.0466700, 5.0550000, 5.0633300, 5.0716700, 5.0800000, 5.0883300, 5.0966700,
1058  5.1050000, 5.1133300, 5.1216700, 5.1300000, 5.1383300, 5.1466700, 5.1550000,
1059  5.1633300, 5.1716700, 5.1800000, 5.1883300, 5.1966700, 5.2050000, 5.2133300,
1060  5.2216700, 5.2300000, 5.2383300, 5.2466700, 5.2550000, 5.2633300, 5.2716700,
1061  5.2800000, 5.2883300, 5.2966700, 5.3050000, 5.3133300, 5.3216700, 5.3300000,
1062  5.3383300, 5.3466700, 5.3550000, 5.3633300, 5.3716700, 5.3800000, 5.3883300,
1063  5.3966700, 5.4050000, 5.4133300, 5.4216700, 5.4300000, 5.4383300, 5.4466700,
1064  5.4550000, 5.4633300, 5.4716700, 5.4800000, 5.4883300, 5.4966700, 5.5050000,
1065  5.5133300, 5.5216700, 5.5300000, 5.5383300, 5.5466700, 5.5550000, 5.5633300,
1066  5.5716700, 5.5800000, 5.5883300, 5.5966700, 5.6050000, 5.6133300, 5.6216700,
1067  5.6300000, 5.6383300, 5.6466700, 5.6550000, 5.6633300, 5.6716700, 5.6800000,
1068  5.6883300, 5.6966700, 5.7050000, 5.7133300, 5.7216700, 5.7300000, 5.7383300,
1069  5.7466700, 5.7550000, 5.7633300, 5.7716700, 5.7800000, 5.7883300, 5.7966700,
1070  5.8050000, 5.8133300, 5.8216700, 5.8300000, 5.8383300, 5.8466700, 5.8550000,
1071  5.8633300, 5.8716700, 5.8800000, 5.8883300, 5.8966700, 5.9050000, 5.9133300,
1072  5.9216700, 5.9300000, 5.9383300, 5.9466700, 5.9550000, 5.9633300, 5.9716700,
1073  5.9800000, 5.9883300, 5.9966700, 6.0050000, 6.0133300, 6.0216700, 6.0300000,
1074  6.0383300, 6.0466700, 6.0550000, 6.0633300, 6.0716700, 6.0800000};
1075 
1079 const std::initializer_list<double> PIMINUSP_TOT_SIG = {
1080  0.152858, 0.327257, 0.691281, 1.295780, 2.171295, 3.386768, 5.008652,
1081  7.163028, 10.009071, 13.673280, 18.337398, 24.205797, 31.337485, 39.349561,
1082  47.474009, 54.547896, 59.146810, 60.567825, 58.976519, 55.362913, 50.786356,
1083  46.100703, 41.835492, 38.112119, 35.041374, 32.559396, 30.666500, 30.553421,
1084  29.535834, 28.802318, 28.313655, 28.032401, 28.055208, 28.090181, 28.137784,
1085  28.269045, 28.295555, 28.379242, 28.427585, 28.479308, 28.552234, 28.732221,
1086  29.075600, 29.568220, 30.403254, 31.644870, 33.252611, 35.283816, 37.447700,
1087  40.095296, 41.869358, 42.524230, 42.374609, 41.784164, 40.663023, 39.329058,
1088  38.163191, 37.325262, 36.813379, 36.643537, 36.891949, 37.565398, 38.555171,
1089  39.915579, 41.583730, 43.604668, 46.118370, 48.847730, 51.552357, 53.970086,
1090  55.750872, 56.410438, 55.732488, 54.095957, 52.248099, 50.454878, 48.464654,
1091  46.396081, 44.321021, 42.528800, 41.007405, 39.883259, 38.905444, 38.088952,
1092  37.510993, 37.135283, 36.952662, 36.917412, 36.883476, 36.847426, 36.894017,
1093  36.917412, 36.918139, 36.962335, 37.068928, 37.154242, 37.180553, 37.154242,
1094  37.068928, 36.859807, 36.643537, 36.401634, 36.075176, 35.702515, 35.319773,
1095  34.830399, 34.297869, 33.816300, 33.360011, 33.017105, 32.772630, 32.553345,
1096  32.440594, 32.350997, 32.290537, 32.350997, 32.365473, 32.365473, 32.551776,
1097  32.743464, 32.874786, 33.109988, 33.345355, 33.610224, 33.816300, 34.091913,
1098  34.317933, 34.559182, 34.746146, 34.931242, 35.068813, 35.165889, 35.221325,
1099  35.298935, 35.275177, 35.210237, 35.144506, 35.043288, 35.022017, 34.947613,
1100  34.897018, 34.813921, 34.751650, 34.712273, 34.656700, 34.610715, 34.507930,
1101  34.471322, 34.400888, 34.395470, 34.284620, 34.236712, 34.197676, 34.130444,
1102  34.098418, 34.006726, 33.952272, 33.940617, 33.843494, 33.816300, 33.712716,
1103  33.716827, 33.653526, 33.599399, 33.539859, 33.491145, 33.423650, 33.365154,
1104  33.329156, 33.270394, 33.236816, 33.197860, 33.151440, 33.080759, 33.029387,
1105  32.996307, 32.970622, 32.880007, 32.855007, 32.822631, 32.746829, 32.714297,
1106  32.685130, 32.667349, 32.611272, 32.578447, 32.505949, 32.489631, 32.466293,
1107  32.378776, 32.323932, 32.297227, 32.254005, 32.218483, 32.187546, 32.122138,
1108  32.092771, 32.057976, 32.027580, 31.941751, 31.920623, 31.900251, 31.867858,
1109  31.833146, 31.794096, 31.746677, 31.730871, 31.644153, 31.602765, 31.604736,
1110  31.541747, 31.501713, 31.473118, 31.478837, 31.415927, 31.396316, 31.332255,
1111  31.320053, 31.259043, 31.247865, 31.180798, 31.147264, 31.126628, 31.102553,
1112  31.085208, 31.040114, 30.995749, 30.968754, 30.918977, 30.881578, 30.858871,
1113  30.821891, 30.790750, 30.790750, 30.702000, 30.738979, 30.668718, 30.645144,
1114  30.608098, 30.563141, 30.557977, 30.540764, 30.523551, 30.471427, 30.439310,
1115  30.414960, 30.349764, 30.366302, 30.306723, 30.306723, 30.229656, 30.213892,
1116  30.223234, 30.185868, 30.164535, 30.120609, 30.094985, 30.069361, 30.026206,
1117  30.018114, 29.996207, 29.969042, 29.927909, 29.898843, 29.873285, 29.864765,
1118  29.818878, 29.769167, 29.826526, 29.750048, 29.739619, 29.711809, 29.629659,
1119  29.614723, 29.610100, 29.581040, 29.559245, 29.546564, 29.489013, 29.478087,
1120  29.428813, 29.454628, 29.427364, 29.407074, 29.379478, 29.336245, 29.350158,
1121  29.296690, 29.296690, 29.272160, 29.255296, 29.230065, 29.212043, 29.201230,
1122  29.129142, 29.103911, 29.103911, 29.095029, 29.022607, 29.058613, 29.028414,
1123  28.987762, 28.978083, 28.952918, 28.911085, 28.905855, 28.877618, 28.839968,
1124  28.843391, 28.814868, 28.802318, 28.779209, 28.784646, 28.727214, 28.693076,
1125  28.670886, 28.670886, 28.662839, 28.652110, 28.622147, 28.625673, 28.536868,
1126  28.582197, 28.536868, 28.530386, 28.502296, 28.472411, 28.446262, 28.449559,
1127  28.408907, 28.388032, 28.379242, 28.372797, 28.352874, 28.331584, 28.319756,
1128  28.341410, 28.267714, 28.220403, 28.253521, 28.213159, 28.203844, 28.187329,
1129  28.161377, 28.135243, 28.137784, 28.102146, 28.120236, 28.075476, 28.055208,
1130  28.045325, 28.018147, 28.034030, 28.001301, 27.995910, 27.959908, 27.906964,
1131  27.906964, 27.906964, 27.897107, 27.898750, 27.894643, 27.823798, 27.836559,
1132  27.839759, 27.798540, 27.771969, 27.773898, 27.750439, 27.716982, 27.718897,
1133  27.722248, 27.707502, 27.665981, 27.666951, 27.630086, 27.631315, 27.611654,
1134  27.603484, 27.603484, 27.574888, 27.582241, 27.577716, 27.521781, 27.523415,
1135  27.489099, 27.464588, 27.464588, 27.464588, 27.415697, 27.387392, 27.415697,
1136  27.396893, 27.376585, 27.336249, 27.342361, 27.317915, 27.317915, 27.305725,
1137  27.317915, 27.269155, 27.291319, 27.202981, 27.220395, 27.249651, 27.230147,
1138  27.192532, 27.171635, 27.171635, 27.171635, 27.123006, 27.135163, 27.150794,
1139  27.089573, 27.114901, 27.062219, 27.060074, 27.050062, 27.036168, 27.036969,
1140  27.025747, 27.017189, 27.012520, 27.004962, 26.996648, 26.962093, 26.941513,
1141  26.947404, 26.909351, 26.891444, 26.912584, 26.880252, 26.880252, 26.863181,
1142  26.849159, 26.864130, 26.798238, 26.807701, 26.793191, 26.799640, 26.771425,
1143  26.748969, 26.751272, 26.748341, 26.735150, 26.711032, 26.704685, 26.712887,
1144  26.691108, 26.669373, 26.646098, 26.666603, 26.635662, 26.634966, 26.611113,
1145  26.615977, 26.590440, 26.598480, 26.582845, 26.542335, 26.558370, 26.503850,
1146  26.560058, 26.531016, 26.534934, 26.532713, 26.501630, 26.472363, 26.506888,
1147  26.455744, 26.453719, 26.455744, 26.422136, 26.426497, 26.439581, 26.428133,
1148  26.396303, 26.353601, 26.386155, 26.392152, 26.350174, 26.361462, 26.330984,
1149  26.338180, 26.315906, 26.335412, 26.302199, 26.295023, 26.302199, 26.287091,
1150  26.270303, 26.271982, 26.236957, 26.239404, 26.243098, 26.230433, 26.201727,
1151  26.223909, 26.210861, 26.180749, 26.167111, 26.194550, 26.170629, 26.152161,
1152  26.158668, 26.130040, 26.144354, 26.127995, 26.131829, 26.128533, 26.097322,
1153  26.095050, 26.055290, 26.075170, 26.063242, 26.063242, 26.015529, 26.029843,
1154  26.029843, 26.029843, 26.015529, 26.015529, 25.971607, 26.005333, 25.979842,
1155  25.974744, 25.979842, 25.947554, 25.954352, 25.944156, 25.913567, 25.937667,
1156  25.920365, 25.915606, 25.929881, 25.894744, 25.872783, 25.872783, 25.872783,
1157  25.844312, 25.849057, 25.854988, 25.830077, 25.837194, 25.832110, 25.830077,
1158  25.789045, 25.808724, 25.773135, 25.809514, 25.789045, 25.780672, 25.742808,
1159  25.744665, 25.782875, 25.740597, 25.766018, 25.750765, 25.730429, 25.702037,
1160  25.717524, 25.704618, 25.680325, 25.689869, 25.697669, 25.710149, 25.685260,
1161  25.673645, 25.640091, 25.667336, 25.625827, 25.648241, 25.615509, 25.613520,
1162  25.632149, 25.618044, 25.616861, 25.620015, 25.588469, 25.588469, 25.573567,
1163  25.588469, 25.571813, 25.544909, 25.535381, 25.568245, 25.494090, 25.536989,
1164  25.499989, 25.524119, 25.505887, 25.517685, 25.531841, 25.464596, 25.496449,
1165  25.494090, 25.459770, 25.482292, 25.458698, 25.461057, 25.469253};
1166 
1168 static std::unique_ptr<InterpolateDataLinear<double>>
1170 
1172 const std::initializer_list<double> PIPLUSPIMINUS_TOT_SQRTS = {
1173  0.2825000, 0.2882500, 0.2965000, 0.3047500, 0.3130000, 0.3212500, 0.3295000,
1174  0.3377500, 0.3460000, 0.3542500, 0.3625000, 0.3707500, 0.3790000, 0.3872500,
1175  0.3955000, 0.4037500, 0.4120000, 0.4202500, 0.4285000, 0.4367500, 0.4450000,
1176  0.4532500, 0.4615000, 0.4697500, 0.4780000, 0.4862500, 0.4945000, 0.5027500,
1177  0.5110000, 0.5192500, 0.5275000, 0.5357500, 0.5440000, 0.5522500, 0.5605000,
1178  0.5687500, 0.5770000, 0.5852500, 0.5935000, 0.6017500, 0.6100000, 0.6182500,
1179  0.6265000, 0.6347500, 0.6430000, 0.6512500, 0.6595000, 0.6677500, 0.6760000,
1180  0.6842500, 0.6925000, 0.7007500, 0.7090000, 0.7172500, 0.7255000, 0.7337500,
1181  0.7420000, 0.7502500, 0.7585000, 0.7667500, 0.7750000, 0.7832500, 0.7915000,
1182  0.7997500, 0.8080000, 0.8162500, 0.8245000, 0.8327500, 0.8410000, 0.8492500,
1183  0.8575000, 0.8657500, 0.8740000, 0.8822500, 0.8905000, 0.8987500, 0.9070000,
1184  0.9152500, 0.9235000, 0.9317500, 0.9400000, 0.9482500, 0.9565000, 0.9647500,
1185  0.9730000, 0.9812500, 0.9895000, 0.9977500, 1.0060000, 1.0142500, 1.0225000,
1186  1.0307500, 1.0390000, 1.0472500, 1.0555000, 1.0637500, 1.0720000, 1.0802500,
1187  1.0885000, 1.0967500, 1.1050000, 1.1132500, 1.1215000, 1.1297500, 1.1380000,
1188  1.1462500, 1.1545000, 1.1627500, 1.1710000, 1.1792500, 1.1875000, 1.1957500,
1189  1.2040000, 1.2122500, 1.2205000, 1.2287500, 1.2370000, 1.2452500, 1.2535000,
1190  1.2617500, 1.2700000, 1.2782500, 1.2865000, 1.2947500, 1.3030000, 1.3112500,
1191  1.3195000, 1.3277500, 1.3360000, 1.3442500, 1.3525000, 1.3607500, 1.3690000,
1192  1.3772500, 1.3855000, 1.3937500, 1.4020000, 1.4102500, 1.4185000, 1.4267500,
1193  1.4350000, 1.4432500, 1.4515000, 1.4597500, 1.4680000, 1.4762500, 1.4845000,
1194  1.4927500, 1.5010000, 1.5092500, 1.5175000, 1.5257500, 1.5340000, 1.5422500,
1195  1.5505000, 1.5587500, 1.5670000, 1.5752500, 1.5835000, 1.5917500, 1.6000000,
1196  1.6082500, 1.6165000, 1.6247500, 1.6330000, 1.6412500, 1.6495000, 1.6577500,
1197  1.6660000, 1.6742500, 1.6825000, 1.6907500, 1.6990000, 1.7072500, 1.7155000,
1198  1.7237500, 1.7320000, 1.7402500, 1.7485000, 1.7567500, 1.7650000, 1.7732500,
1199  1.7815000, 1.7897500, 1.7980000, 1.8062500, 1.8145000, 1.8227500, 1.8310000,
1200  1.8392500, 1.8475000, 1.8557500, 1.8640000, 1.8722500, 1.8805000, 1.8887500,
1201  1.8970000, 1.9052500, 1.9135000, 1.9217500, 1.9300000, 1.9382500, 1.9465000,
1202  1.9547500, 1.9630000, 1.9712500, 1.9795000, 1.9877500, 1.9960000, 2.0042500,
1203  2.0125000, 2.0207500, 2.0290000, 2.0372500, 2.0455000, 2.0537500, 2.0620000,
1204  2.0702500, 2.0785000, 2.0867500, 2.0950000, 2.1032500, 2.1115000, 2.1197500,
1205  2.1280000, 2.1362500, 2.1445000, 2.1527500, 2.1610000, 2.1692500, 2.1775000,
1206  2.1857500, 2.1940000, 2.2022500, 2.2105000, 2.2187500, 2.2270000, 2.2352500,
1207  2.2435000, 2.2517500, 2.2600000, 2.2682500, 2.2765000, 2.2847500, 2.2930000,
1208  2.3012500, 2.3095000, 2.3177500, 2.3260000, 2.3342500, 2.3425000, 2.3507500,
1209  2.3590000, 2.3672500, 2.3755000, 2.3837500, 2.3920000, 2.4002500, 2.4085000,
1210  2.4167500, 2.4250000, 2.4332500, 2.4415000, 2.4497500, 2.4580000, 2.4662500,
1211  2.4745000, 2.4827500, 2.4910000, 2.4992500, 2.5075000, 2.5157500, 2.5240000,
1212  2.5322500, 2.5405000, 2.5487500, 2.5570000, 2.5652500, 2.5735000, 2.5817500,
1213  2.5900000, 2.5982500, 2.6065000, 2.6147500, 2.6230000, 2.6312500, 2.6395000,
1214  2.6477500, 2.6560000, 2.6642500, 2.6725000, 2.6807500, 2.6890000, 2.6972500,
1215  2.7055000, 2.7137500, 2.7220000, 2.7302500, 2.7385000, 2.7467500, 2.7550000,
1216  2.7632500, 2.7715000, 2.7797500, 2.7880000, 2.7962500, 2.8045000, 2.8127500,
1217  2.8210000, 2.8292500, 2.8375000, 2.8457500, 2.8540000, 2.8622500, 2.8705000,
1218  2.8787500, 2.8870000, 2.8952500, 2.9035000, 2.9117500, 2.9200000, 2.9282500,
1219  2.9365000, 2.9447500, 2.9530000, 2.9612500, 2.9695000, 2.9777500, 2.9860000,
1220  2.9942500, 3.0025000, 3.0107500, 3.0190000, 3.0272500, 3.0355000, 3.0437500,
1221  3.0520000, 3.0602500, 3.0685000, 3.0767500, 3.0850000, 3.0932500, 3.1015000,
1222  3.1097500, 3.1180000, 3.1262500, 3.1345000, 3.1427500, 3.1510000, 3.1592500,
1223  3.1675000, 3.1757500, 3.1840000, 3.1922500, 3.2005000, 3.2087500, 3.2170000,
1224  3.2252500, 3.2335000, 3.2417500, 3.2500000, 3.2582500, 3.2665000, 3.2747500,
1225  3.2830000, 3.2912500, 3.2995000, 3.3077500, 3.3160000, 3.3242500, 3.3325000,
1226  3.3407500, 3.3490000, 3.3572500, 3.3655000, 3.3737500, 3.3820000, 3.3902500,
1227  3.3985000, 3.4067500, 3.4150000, 3.4232500, 3.4315000, 3.4397500, 3.4480000,
1228  3.4562500, 3.4645000, 3.4727500, 3.4810000, 3.4892500, 3.4975000, 3.5057500,
1229  3.5140000, 3.5222500, 3.5305000, 3.5387500, 3.5470000, 3.5552500, 3.5635000,
1230  3.5717500, 3.5800000};
1231 
1236 const std::initializer_list<double> PIPLUSPIMINUS_TOT_SIG = {
1237  7.542964, 7.673451, 7.808281, 7.968326, 8.162403, 8.344765,
1238  8.547502, 8.800306, 9.065026, 9.321825, 9.611654, 9.940196,
1239  10.288609, 10.613969, 11.005539, 11.404178, 11.848546, 12.323925,
1240  12.794033, 13.314148, 13.840777, 14.420107, 15.045962, 15.647782,
1241  16.342663, 17.048769, 17.839317, 18.666885, 19.565381, 20.496500,
1242  21.491759, 22.596055, 23.730551, 24.971220, 26.302199, 27.759113,
1243  29.335179, 31.024504, 32.865649, 34.876427, 37.040490, 39.458460,
1244  42.161261, 45.072788, 48.319055, 51.951732, 55.990250, 60.480567,
1245  65.564842, 71.157859, 77.416582, 84.277474, 91.768578, 99.818020,
1246  108.102989, 116.396756, 123.905819, 130.079123, 134.203733, 135.801247,
1247  134.727977, 131.660758, 124.644799, 116.878612, 108.472427, 99.818020,
1248  91.346028, 83.468975, 76.158023, 69.564613, 63.570324, 58.305449,
1249  53.511879, 49.283931, 45.496384, 42.190536, 39.169784, 36.643537,
1250  34.352628, 32.365473, 30.730350, 29.535218, 28.812358, 28.802318,
1251  29.711809, 31.259043, 28.587904, 24.825775, 22.661647, 21.015246,
1252  19.699897, 18.614427, 17.671459, 16.856412, 16.250363, 15.718505,
1253  15.299725, 14.957123, 14.741138, 14.633735, 14.548127, 14.633735,
1254  14.813002, 15.042436, 15.406763, 15.851665, 16.456153, 17.135101,
1255  17.971219, 18.960641, 20.033004, 21.269245, 22.586906, 23.973359,
1256  25.402783, 26.802864, 28.055208, 29.003249, 29.830002, 30.141108,
1257  30.155673, 29.781334, 29.124095, 28.278359, 27.220395, 26.187374,
1258  25.024553, 23.976570, 22.907156, 21.927853, 21.035650, 20.223063,
1259  19.491043, 18.844883, 18.259722, 17.750129, 20.987447, 20.590816,
1260  20.208454, 19.880666, 19.651069, 19.473283, 19.318242, 19.236168,
1261  19.185035, 19.229350, 19.236168, 19.256687, 19.243410, 19.163981,
1262  18.991124, 18.780798, 18.601081, 18.385386, 18.253458, 18.195077,
1263  18.214169, 18.329400, 18.531624, 18.844883, 19.148512, 19.606680,
1264  20.093646, 20.643847, 21.245482, 21.872787, 22.564686, 23.256580,
1265  23.936686, 24.501625, 25.015219, 25.316585, 25.376313, 25.235334,
1266  24.818859, 24.224968, 23.482239, 22.574209, 21.698391, 20.792680,
1267  19.947372, 19.169231, 18.410461, 17.782089, 17.203361, 16.653492,
1268  16.250735, 15.880216, 15.600683, 15.382828, 15.264160, 15.130098,
1269  15.055833, 15.024986, 14.997841, 15.047607, 15.089920, 15.195184,
1270  15.289821, 15.433144, 15.552319, 15.703251, 15.850734, 16.013400,
1271  16.180642, 16.353993, 16.512996, 16.671103, 16.825781, 16.971669,
1272  17.103896, 17.198316, 17.286529, 17.371720, 17.417153, 17.436625,
1273  17.491787, 17.475698, 17.541506, 17.553845, 17.595356, 17.636175,
1274  17.656757, 17.687193, 17.776353, 17.807680, 17.861820, 17.947461,
1275  18.026655, 18.102499, 18.166941, 18.254545, 18.332066, 18.397422,
1276  18.475657, 18.555148, 18.614427, 18.694648, 18.747651, 18.819145,
1277  18.861595, 18.892053, 18.977576, 19.028763, 19.072674, 19.100573,
1278  19.128790, 19.191990, 19.201106, 19.200075, 19.198409, 19.236168,
1279  19.236168, 19.236168, 19.214512, 19.195262, 19.205489, 19.184497,
1280  19.158662, 19.148512, 19.137993, 19.093062, 19.088985, 19.018307,
1281  19.052287, 19.000534, 19.004716, 18.991124, 18.991124, 18.954014,
1282  18.960641, 18.913530, 18.880276, 18.869191, 18.860509, 18.830810,
1283  18.801668, 18.785048, 18.779635, 18.761155, 18.747651, 18.740524,
1284  18.714610, 18.682417, 18.663779, 18.646694, 18.649218, 18.626503,
1285  18.626503, 18.607925, 18.541974, 18.541264, 18.533614, 18.515037,
1286  18.505748, 18.505748, 18.465627, 18.471359, 18.433531, 18.436969,
1287  18.415476, 18.397422, 18.371272, 18.370389, 18.349395, 18.330854,
1288  18.343396, 18.299693, 18.309041, 18.272914, 18.265416, 18.231251,
1289  18.215592, 18.238843, 18.232804, 18.193670, 18.195663, 18.153313,
1290  18.145839, 18.128813, 18.119354, 18.108594, 18.094760, 18.062410,
1291  18.062928, 18.055032, 18.026655, 18.021490, 18.026655, 17.993657,
1292  17.973859, 17.957360, 17.974387, 17.931622, 17.915288, 17.928827,
1293  17.907864, 17.880541, 17.899406, 17.860504, 17.848664, 17.826854,
1294  17.838217, 17.848664, 17.789465, 17.795696, 17.789465, 17.777043,
1295  17.762233, 17.763241, 17.745212, 17.745212, 17.706166, 17.704238,
1296  17.678834, 17.671459, 17.671459, 17.664540, 17.657345, 17.635663,
1297  17.583249, 17.612652, 17.606771, 17.605301, 17.612652, 17.579982,
1298  17.564537, 17.553845, 17.524540, 17.530401, 17.514772, 17.498026,
1299  17.503608, 17.500563, 17.471101, 17.465930, 17.445642, 17.461743,
1300  17.436625, 17.407418, 17.394142, 17.417153, 17.407418, 17.353176,
1301  17.363607, 17.389893, 17.358739, 17.354845, 17.334400, 17.319797,
1302  17.328783, 17.300391, 17.290688, 17.283970, 17.266057, 17.280985,
1303  17.266057, 17.232470, 17.268048, 17.238292, 17.203361};
1304 
1306 static std::unique_ptr<InterpolateDataLinear<double>>
1308 
1310 const std::initializer_list<double> PIZEROPIZERO_TOT_SQRTS = {
1311  0.2825000, 0.2882500, 0.2965000, 0.3047500, 0.3130000, 0.3212500, 0.3295000,
1312  0.3377500, 0.3460000, 0.3542500, 0.3625000, 0.3707500, 0.3790000, 0.3872500,
1313  0.3955000, 0.4037500, 0.4120000, 0.4202500, 0.4285000, 0.4367500, 0.4450000,
1314  0.4532500, 0.4615000, 0.4697500, 0.4780000, 0.4862500, 0.4945000, 0.5027500,
1315  0.5110000, 0.5192500, 0.5275000, 0.5357500, 0.5440000, 0.5522500, 0.5605000,
1316  0.5687500, 0.5770000, 0.5852500, 0.5935000, 0.6017500, 0.6100000, 0.6182500,
1317  0.6265000, 0.6347500, 0.6430000, 0.6512500, 0.6595000, 0.6677500, 0.6760000,
1318  0.6842500, 0.6925000, 0.7007500, 0.7090000, 0.7172500, 0.7255000, 0.7337500,
1319  0.7420000, 0.7502500, 0.7585000, 0.7667500, 0.7750000, 0.7832500, 0.7915000,
1320  0.7997500, 0.8080000, 0.8162500, 0.8245000, 0.8327500, 0.8410000, 0.8492500,
1321  0.8575000, 0.8657500, 0.8740000, 0.8822500, 0.8905000, 0.8987500, 0.9070000,
1322  0.9152500, 0.9235000, 0.9317500, 0.9400000, 0.9482500, 0.9565000, 0.9647500,
1323  0.9730000, 0.9812500, 0.9895000, 0.9977500, 1.0060000, 1.0142500, 1.0225000,
1324  1.0307500, 1.0390000, 1.0472500, 1.0555000, 1.0637500, 1.0720000, 1.0802500,
1325  1.0885000, 1.0967500, 1.1050000, 1.1132500, 1.1215000, 1.1297500, 1.1380000,
1326  1.1462500, 1.1545000, 1.1627500, 1.1710000, 1.1792500, 1.1875000, 1.1957500,
1327  1.2040000, 1.2122500, 1.2205000, 1.2287500, 1.2370000, 1.2452500, 1.2535000,
1328  1.2617500, 1.2700000, 1.2782500, 1.2865000, 1.2947500, 1.3030000, 1.3112500,
1329  1.3195000, 1.3277500, 1.3360000, 1.3442500, 1.3525000, 1.3607500, 1.3690000,
1330  1.3772500, 1.3855000, 1.3937500, 1.4020000, 1.4102500, 1.4185000, 1.4267500,
1331  1.4350000, 1.4432500, 1.4515000, 1.4597500, 1.4680000, 1.4762500, 1.4845000,
1332  1.4927500, 1.5010000, 1.5092500, 1.5175000, 1.5257500, 1.5340000, 1.5422500,
1333  1.5505000, 1.5587500, 1.5670000, 1.5752500, 1.5835000, 1.5917500, 1.6000000,
1334  1.6082500, 1.6165000, 1.6247500, 1.6330000, 1.6412500, 1.6495000, 1.6577500,
1335  1.6660000, 1.6742500, 1.6825000, 1.6907500, 1.6990000, 1.7072500, 1.7155000,
1336  1.7237500, 1.7320000, 1.7402500, 1.7485000, 1.7567500, 1.7650000, 1.7732500,
1337  1.7815000, 1.7897500, 1.7980000, 1.8062500, 1.8145000, 1.8227500, 1.8310000,
1338  1.8392500, 1.8475000, 1.8557500, 1.8640000, 1.8722500, 1.8805000, 1.8887500,
1339  1.8970000, 1.9052500, 1.9135000, 1.9217500, 1.9300000, 1.9382500, 1.9465000,
1340  1.9547500, 1.9630000, 1.9712500, 1.9795000, 1.9877500, 1.9960000, 2.0042500,
1341  2.0125000, 2.0207500, 2.0290000, 2.0372500, 2.0455000, 2.0537500, 2.0620000,
1342  2.0702500, 2.0785000, 2.0867500, 2.0950000, 2.1032500, 2.1115000, 2.1197500,
1343  2.1280000, 2.1362500, 2.1445000, 2.1527500, 2.1610000, 2.1692500, 2.1775000,
1344  2.1857500, 2.1940000, 2.2022500, 2.2105000, 2.2187500, 2.2270000, 2.2352500,
1345  2.2435000, 2.2517500, 2.2600000, 2.2682500, 2.2765000, 2.2847500, 2.2930000,
1346  2.3012500, 2.3095000, 2.3177500, 2.3260000, 2.3342500, 2.3425000, 2.3507500,
1347  2.3590000, 2.3672500, 2.3755000, 2.3837500, 2.3920000, 2.4002500, 2.4085000,
1348  2.4167500, 2.4250000, 2.4332500, 2.4415000, 2.4497500, 2.4580000, 2.4662500,
1349  2.4745000, 2.4827500, 2.4910000, 2.4992500, 2.5075000, 2.5157500, 2.5240000,
1350  2.5322500, 2.5405000, 2.5487500, 2.5570000, 2.5652500, 2.5735000, 2.5817500,
1351  2.5900000, 2.5982500, 2.6065000, 2.6147500, 2.6230000, 2.6312500, 2.6395000,
1352  2.6477500, 2.6560000, 2.6642500, 2.6725000, 2.6807500, 2.6890000, 2.6972500,
1353  2.7055000, 2.7137500, 2.7220000, 2.7302500, 2.7385000, 2.7467500, 2.7550000,
1354  2.7632500, 2.7715000, 2.7797500, 2.7880000, 2.7962500, 2.8045000, 2.8127500,
1355  2.8210000, 2.8292500, 2.8375000, 2.8457500, 2.8540000, 2.8622500, 2.8705000,
1356  2.8787500, 2.8870000, 2.8952500, 2.9035000, 2.9117500, 2.9200000, 2.9282500,
1357  2.9365000, 2.9447500, 2.9530000, 2.9612500, 2.9695000, 2.9777500, 2.9860000,
1358  2.9942500, 3.0025000, 3.0107500, 3.0190000, 3.0272500, 3.0355000, 3.0437500,
1359  3.0520000, 3.0602500, 3.0685000, 3.0767500, 3.0850000, 3.0932500, 3.1015000,
1360  3.1097500, 3.1180000, 3.1262500, 3.1345000, 3.1427500, 3.1510000, 3.1592500,
1361  3.1675000, 3.1757500, 3.1840000, 3.1922500, 3.2005000, 3.2087500, 3.2170000,
1362  3.2252500, 3.2335000, 3.2417500, 3.2500000, 3.2582500, 3.2665000, 3.2747500,
1363  3.2830000, 3.2912500, 3.2995000, 3.3077500, 3.3160000, 3.3242500, 3.3325000,
1364  3.3407500, 3.3490000, 3.3572500, 3.3655000, 3.3737500, 3.3820000, 3.3902500,
1365  3.3985000, 3.4067500, 3.4150000, 3.4232500, 3.4315000, 3.4397500, 3.4480000,
1366  3.4562500, 3.4645000, 3.4727500, 3.4810000, 3.4892500, 3.4975000, 3.5057500,
1367  3.5140000, 3.5222500, 3.5305000, 3.5387500, 3.5470000, 3.5552500, 3.5635000,
1368  3.5717500, 3.5800000};
1369 
1374 const std::initializer_list<double> PIZEROPIZERO_TOT_SIG = {
1375  12.539807, 12.634504, 12.767629, 12.867964, 12.993971, 13.127818, 13.273229,
1376  13.375527, 13.564121, 13.684778, 13.860326, 14.039581, 14.208049, 14.386911,
1377  14.548127, 14.799936, 14.993316, 15.229361, 15.439161, 15.685097, 15.926147,
1378  16.203907, 16.502661, 16.779835, 17.076805, 17.404762, 17.734395, 18.132596,
1379  18.445567, 18.794397, 19.222533, 19.606680, 20.049162, 20.484951, 20.931352,
1380  21.382465, 21.903969, 22.412823, 22.921125, 23.451911, 23.974391, 24.567879,
1381  25.040152, 25.659449, 26.158668, 26.711032, 27.244775, 27.744367, 28.259730,
1382  28.677145, 29.117673, 29.483160, 29.788287, 30.046067, 30.242995, 30.373537,
1383  30.444857, 30.377499, 30.254851, 30.137691, 30.018114, 30.501173, 29.255296,
1384  28.602172, 28.042854, 27.433158, 26.807701, 26.115726, 25.446900, 24.744951,
1385  24.052819, 23.370504, 22.689119, 22.035327, 21.417861, 20.803332, 20.260109,
1386  19.780842, 19.383980, 19.113450, 18.887950, 18.991124, 19.400447, 20.358306,
1387  22.200114, 24.667611, 22.661647, 19.555056, 18.035168, 16.925566, 16.071858,
1388  15.467239, 14.930075, 14.491185, 14.168509, 13.941210, 13.760948, 13.647822,
1389  13.600115, 13.643387, 13.731991, 13.927793, 14.175606, 14.484077, 14.947287,
1390  15.430521, 16.013892, 16.716109, 17.571940, 18.505748, 19.568558, 20.775941,
1391  22.035327, 23.298352, 24.683079, 25.948914, 27.025747, 27.796075, 28.390229,
1392  28.577203, 28.402681, 27.894643, 27.052272, 26.039385, 24.861288, 23.661747,
1393  22.412823, 21.182706, 19.994666, 18.893577, 17.884184, 16.948617, 16.102798,
1394  15.304001, 14.633735, 13.997562, 13.439709, 12.968691, 12.468981, 12.133207,
1395  11.834085, 11.610181, 11.411473, 11.302499, 11.215682, 11.215682, 11.253303,
1396  11.253303, 11.122023, 11.050702, 10.806568, 10.560053, 10.255979, 9.947575,
1397  9.684172, 9.481768, 9.311605, 9.177888, 9.160884, 9.143959, 9.160884,
1398  9.197321, 9.267257, 9.340849, 9.438670, 9.503318, 9.605723, 9.735141,
1399  9.852035, 9.978147, 10.142061, 10.278820, 10.386891, 10.577448, 10.725117,
1400  10.843403, 10.935884, 11.028757, 11.106479, 11.186864, 11.247033, 11.309734,
1401  11.442112, 11.528316, 11.618150, 11.785881, 11.946822, 12.076282, 12.222857,
1402  12.409723, 12.612382, 12.784351, 12.988915, 13.171324, 13.392642, 13.596084,
1403  13.801679, 14.032578, 14.260965, 14.467492, 14.741138, 14.910756, 15.154864,
1404  15.393804, 15.614501, 15.836769, 16.060607, 16.243678, 16.409645, 16.619470,
1405  16.776007, 16.899633, 17.029494, 17.135670, 17.203361, 17.284866, 17.319797,
1406  17.362280, 17.419935, 17.428836, 17.470116, 17.500563, 17.564537, 17.612652,
1407  17.671459, 17.750129, 17.842086, 17.899406, 17.973859, 18.026655, 18.135004,
1408  18.209145, 18.285411, 18.374479, 18.445567, 18.556062, 18.586251, 18.662847,
1409  18.738997, 18.828677, 18.869191, 18.896287, 18.945399, 19.072674, 19.098159,
1410  19.113450, 19.113450, 19.151209, 19.167991, 19.222533, 19.222533, 19.236168,
1411  19.236168, 19.236168, 19.222533, 19.227403, 19.195262, 19.195262, 19.192340,
1412  19.144129, 19.132329, 19.104712, 19.072674, 19.046727, 19.072674, 19.027822,
1413  19.013365, 18.991124, 18.991124, 18.926093, 18.936931, 18.909835, 18.869191,
1414  18.869191, 18.844883, 18.838806, 18.828677, 18.774659, 18.780798, 18.747651,
1415  18.747651, 18.720729, 18.747651, 18.683038, 18.670557, 18.650732, 18.661116,
1416  18.614427, 18.626503, 18.586251, 18.613086, 18.596314, 18.546000, 18.541264,
1417  18.505748, 18.488553, 18.454164, 18.471359, 18.415476, 18.445567, 18.385386,
1418  18.385386, 18.385386, 18.316832, 18.352667, 18.305406, 18.337398, 18.282555,
1419  18.265416, 18.278746, 18.265416, 18.235522, 18.245487, 18.188545, 18.189322,
1420  18.198984, 18.145839, 18.183207, 18.130941, 18.145839, 18.089752, 18.125975,
1421  18.077734, 18.056451, 18.043681, 18.036587, 18.048325, 18.026655, 18.005057,
1422  17.992715, 18.015856, 17.947461, 17.987058, 17.951060, 17.947461, 17.919743,
1423  17.907864, 17.907864, 17.890949, 17.875573, 17.880541, 17.884184, 17.789465,
1424  17.816788, 17.821755, 17.809198, 17.789465, 17.777664, 17.718661, 17.735000,
1425  17.710794, 17.765864, 17.716846, 17.671459, 17.705175, 17.684570, 17.659697,
1426  17.671459, 17.629454, 17.647936, 17.660767, 17.600891, 17.612652, 17.593050,
1427  17.597950, 17.570647, 17.570647, 17.553845, 17.527796, 17.534309, 17.539193,
1428  17.522586, 17.485467, 17.475698, 17.527796, 17.461743, 17.456161, 17.436625,
1429  17.436625, 17.413259, 17.413259, 17.417153, 17.394142, 17.382704, 17.383521,
1430  17.358739, 17.358739, 17.319797, 17.319797, 17.337770, 17.319797, 17.301884,
1431  17.319797, 17.273223, 17.306859, 17.290688, 17.242173, 17.290688, 17.244945,
1432  17.229236, 17.219995};
1433 
1435 static std::unique_ptr<InterpolateDataLinear<double>>
1437 
1438 } // namespace smash
1439 
1440 #endif // SRC_INCLUDE_SMASH_PARAMETRIZATIONS_DATA_H_
Definition: action.h:24
const std::initializer_list< double > KMINUSN_TOT_PLAB
PDG data on K- n total cross section: momentum in lab frame.
const std::initializer_list< double > PIPLUSP_RES_SQRTS
Center-of-mass energy.
const std::initializer_list< double > KMINUSP_RES_SQRTS
Center-of-mass energy list for K̅⁻ N⁺
const std::initializer_list< double > PIMINUSP_ELASTIC_P_LAB
PDG data on pi- p elastic cross section: momentum in lab frame.
const std::initializer_list< double > PIPLUSPIMINUS_TOT_SQRTS
Center-of-mass energy.
static std::unique_ptr< InterpolateDataLinear< double > > piplusp_sigmapluskplus_interpolation
An interpolation that gets lazily filled using the PIPLUSP_SIGMAPLUSKPLUS_SIG data.
static std::unique_ptr< InterpolateDataSpline > piplusp_elastic_res_interpolation
A null interpolation that gets filled using the PIPLUSP_RES data.
const std::initializer_list< double > KPLUSP_TOT_SIG
PDG data on K+ p total cross section: cross section.
const std::initializer_list< double > KMINUSP_TOT_PLAB
PDG smoothed data on K- p total cross section: momentum in lab frame.
const std::initializer_list< double > KMINUSP_RES_SIG
Elastic K̅⁻ N⁺ cross section contributions from decays.
const std::initializer_list< double > PIPLUSP_SIGMAPLUSKPLUS_SIG
PDG data on pi+ p to Sigma+ K+ section: cross section.
const std::initializer_list< double > PIMINUSP_RES_SIG
Elastic π⁻N⁺ cross section contributions from decays.
const std::initializer_list< double > PIPLUSP_SIGMAPLUSKPLUS_P_LAB
PDG data on pi+ p to Sigma+ K+ cross section: momentum in lab frame.
static std::unique_ptr< InterpolateDataSpline > kminusp_elastic_res_interpolation
An interpolation that gets lazily filled using the KMINUSP_RES data.
const std::initializer_list< double > KPLUSN_TOT_SIG
PDG data on K+ n total cross section: cross section.
const std::initializer_list< double > PIPLUSP_RES_SIG
Elastic π⁺N⁺ cross section contributions from decays.
const std::initializer_list< double > PIMINUSP_ELASTIC_SIG
PDG data on pi- p elastic cross section: cross section.
const std::initializer_list< double > KPLUSN_TOT_PLAB
PDG data on K+ n total cross section: momentum in lab frame.
static std::unique_ptr< InterpolateDataLinear< double > > kminusp_total_interpolation
An interpolation that gets lazily filled using the KMINUSP_TOT data.
static std::unique_ptr< InterpolateDataLinear< double > > piplusp_elastic_interpolation
An interpolation that gets lazily filled using the PIPLUSP_ELASTIC_SIG data.
const std::initializer_list< double > PIMINUSP_TOT_SQRTS
Center-of-mass energy.
static std::unique_ptr< InterpolateDataLinear< double > > piminusp_elastic_interpolation
An interpolation that gets lazily filled using the PIMINUSP_ELASTIC data.
const std::initializer_list< double > KMINUSP_ELASTIC_P_LAB
PDG data on K- p elastic cross section: momentum in lab frame.
static std::unique_ptr< InterpolateDataLinear< double > > kplusn_total_interpolation
An interpolation that gets lazily filled using the KPLUSN_TOT data.
const std::initializer_list< double > PIMINUSP_LAMBDAK0_P_LAB
PDG data on pi- p to Lambda K0 cross section: momentum in lab frame.
const std::initializer_list< double > KPLUSP_TOT_PLAB
PDG data on K+ p total cross section: momentum in lab frame.
static std::unique_ptr< InterpolateDataLinear< double > > kplusp_total_interpolation
An interpolation that gets lazily filled using the KPLUSP_TOT data.
const std::initializer_list< double > PIPLUSPIMINUS_TOT_SIG
Total π⁺ π⁻ cross section parametrized from bottom-up SMASH-3.0, using the hadronic list from PDG2018...
static std::unique_ptr< InterpolateDataLinear< double > > kminusp_elastic_interpolation
An interpolation that gets lazily filled using the KMINUSP_ELASTIC data.
const std::initializer_list< double > PIMINUSP_SIGMAMINUSKPLUS_P_LAB
PDG data on pi- p to Sigma- K+ cross section: momentum in lab frame.
const std::initializer_list< double > KMINUSN_TOT_SIG
PDG data on K- n total cross section: cross section.
const std::initializer_list< double > PIMINUSP_SIGMA0K0_RES_SIG
pi- p to Sigma0 K0 cross section: cross section
const std::initializer_list< double > PIPLUSP_ELASTIC_P_LAB
PDG data on pi+ p elastic cross section: momentum in lab frame.
const std::initializer_list< double > PIZEROPIZERO_TOT_SIG
Total π⁰ π⁰ cross section parametrized from bottom-up SMASH-3.0 using the hadronic list from PDG2018.
const std::initializer_list< double > PIMINUSP_TOT_SIG
Total p π⁻ cross section parametrized from bottom-up SMASH-3.0, using the hadronic list from PDG2018.
static std::unique_ptr< InterpolateDataLinear< double > > piplusp_total_interpolation
An interpolation that gets lazily filled using the PIPLUSP_TOT data.
const std::initializer_list< double > PIPLUSP_ELASTIC_SIG
PDG data on pi+ p elastic cross section: cross section.
const std::initializer_list< double > PIMINUSP_SIGMA0K0_RES_SQRTS
pi- p to Sigma0 K0 cross section: square root s
const std::initializer_list< double > PIMINUSP_SIGMAMINUSKPLUS_SIG
PDG data on pi- p to Sigma- K+ cross section: cross section.
const std::initializer_list< double > PIPLUSP_TOT_SQRTS
Center-of-mass energy.
const std::initializer_list< double > KMINUSP_ELASTIC_SIG
PDG data on K- p elastic cross section: cross section.
static std::unique_ptr< InterpolateDataLinear< double > > piminusp_sigma0k0_interpolation
An interpolation that gets lazily filled using the PIMINUSP_SIGMA0K0_RES data.
const std::initializer_list< double > PIPLUSP_TOT_SIG
Total p π⁺ cross section parametrized from bottom-up SMASH-3.0, using the hadronic list from PDG2018.
static std::unique_ptr< InterpolateDataLinear< double > > piminusp_total_interpolation
An interpolation that gets lazily filled using the PIMINUSP_TOT data.
const std::initializer_list< double > PIMINUSP_RES_SQRTS
Center-of-mass energy.
static std::unique_ptr< InterpolateDataLinear< double > > pizeropizero_total_interpolation
An interpolation that gets lazily filled using the PIZEROPIZERO_TOT data.
static std::unique_ptr< InterpolateDataLinear< double > > pipluspiminus_total_interpolation
An interpolation that gets lazily filled using the PIPLUSPIMINUS_TOT data.
static std::unique_ptr< InterpolateDataLinear< double > > kminusn_total_interpolation
An interpolation that gets lazily filled using the KMINUSN_TOT data.
const std::initializer_list< double > KMINUSP_TOT_SIG
PDG smoothed data on K- p total cross section: cross section.
const std::initializer_list< double > PIZEROPIZERO_TOT_SQRTS
Center-of-mass energy.
static std::unique_ptr< InterpolateDataLinear< double > > piminusp_lambdak0_interpolation
An interpolation that gets lazily filled using the PIMINUSP_LAMBDAK0 data.
static std::unique_ptr< InterpolateDataLinear< double > > piminusp_sigmaminuskplus_interpolation
An interpolation that gets lazily filled using the PIMINUSP_SIGMAMINUSKPLUS data.
static std::unique_ptr< InterpolateDataSpline > piminusp_elastic_res_interpolation
An interpolation that gets lazily filled using the PIMINUSP_RES data.
const std::initializer_list< double > PIMINUSP_LAMBDAK0_SIG
PDG data on pi- p to Lambda K0 cross section: cross section.