From 49420614bc7ed4ed8d4f86da1300342af34c1414 Mon Sep 17 00:00:00 2001 From: "C. Weaver" Date: Mon, 27 Dec 2021 11:17:32 -0500 Subject: [PATCH] Check fewer decimal places at the glashow resonance --- test/glashow_resonance.output.txt | 400 +++++++++++++++--------------- test/glashow_resonance.test.cpp | 3 + 2 files changed, 203 insertions(+), 200 deletions(-) diff --git a/test/glashow_resonance.output.txt b/test/glashow_resonance.output.txt index dcadf3fd..e533d849 100644 --- a/test/glashow_resonance.output.txt +++ b/test/glashow_resonance.output.txt @@ -3,203 +3,203 @@ Checking single-differential cross-section Evolving a nue_bar line spectrum at the Glashow resonance Propagating. . . Propagation done -0 10000 0 0 0 0.998632 0.000245023 0.000245023 -1 10351.4 0 0 0 0.998604 0.000245176 0.000245176 -2 10715.2 0 0 0 0.998554 0.000253851 0.000253851 -3 11091.7 0 0 0 0.998503 0.000262835 0.000262835 -4 11481.5 0 0 0 0.99845 0.000272137 0.000272137 -5 11885 0 0 0 0.998396 0.000281771 0.000281771 -6 12302.7 0 0 0 0.998339 0.000291747 0.000291747 -7 12735 0 0 0 0.99828 0.000302078 0.000302078 -8 13182.6 0 0 0 0.99822 0.000312777 0.000312777 -9 13645.8 0 0 0 0.998157 0.000323857 0.000323857 -10 14125.4 0 0 0 0.998092 0.000335332 0.000335332 -11 14621.8 0 0 0 0.998024 0.000347216 0.000347216 -12 15135.6 0 0 0 0.997955 0.000359524 0.000359524 -13 15667.5 0 0 0 0.997882 0.000372272 0.000372272 -14 16218.1 0 0 0 0.997808 0.000385474 0.000385474 -15 16788 0 0 0 0.99773 0.000399149 0.000399149 -16 17378 0 0 0 0.99765 0.000413312 0.000413312 -17 17988.7 0 0 0 0.997567 0.000427982 0.000427982 -18 18620.9 0 0 0 0.997481 0.000443178 0.000443178 -19 19275.2 0 0 0 0.997392 0.000458917 0.000458917 -20 19952.6 0 0 0 0.9973 0.000475221 0.000475221 -21 20653.8 0 0 0 0.997205 0.00049211 0.00049211 -22 21379.6 0 0 0 0.997106 0.000509604 0.000509604 -23 22130.9 0 0 0 0.997004 0.000527727 0.000527727 -24 22908.7 0 0 0 0.996898 0.000546502 0.000546502 -25 23713.7 0 0 0 0.996788 0.000565952 0.000565952 -26 24547.1 0 0 0 0.996675 0.000586102 0.000586102 -27 25409.7 0 0 0 0.996557 0.000606979 0.000606979 -28 26302.7 0 0 0 0.996435 0.000628608 0.000628608 -29 27227 0 0 0 0.996309 0.000651018 0.000651018 -30 28183.8 0 0 0 0.996179 0.000674238 0.000674238 -31 29174.3 0 0 0 0.996044 0.000698298 0.000698298 -32 30199.5 0 0 0 0.995904 0.000723229 0.000723229 -33 31260.8 0 0 0 0.995759 0.000749063 0.000749063 -34 32359.4 0 0 0 0.995609 0.000775834 0.000775834 -35 33496.5 0 0 0 0.995453 0.000803578 0.000803578 -36 34673.7 0 0 0 0.995292 0.00083233 0.00083233 -37 35892.2 0 0 0 0.995126 0.000862129 0.000862129 -38 37153.5 0 0 0 0.994953 0.000893014 0.000893014 -39 38459.2 0 0 0 0.994774 0.000925025 0.000925025 -40 39810.7 0 0 0 0.994589 0.000958206 0.000958206 -41 41209.8 0 0 0 0.994397 0.000992601 0.000992601 -42 42658 0 0 0 0.994198 0.00102826 0.00102826 -43 44157 0 0 0 0.993993 0.00106522 0.00106522 -44 45708.8 0 0 0 0.993779 0.00110354 0.00110354 -45 47315.1 0 0 0 0.993559 0.00114327 0.00114327 -46 48977.9 0 0 0 0.99333 0.00118446 0.00118446 -47 50699.1 0 0 0 0.993093 0.00122718 0.00122718 -48 52480.7 0 0 0 0.992848 0.00127147 0.00127147 -49 54325 0 0 0 0.992594 0.00131741 0.00131741 -50 56234.1 0 0 0 0.99233 0.00136504 0.00136504 -51 58210.3 0 0 0 0.992058 0.00141445 0.00141445 -52 60256 0 0 0 0.991775 0.0014657 0.0014657 -53 62373.5 0 0 0 0.991483 0.00151885 0.00151885 -54 64565.4 0 0 0 0.991179 0.001574 0.001574 -55 66834.4 0 0 0 0.990865 0.00163121 0.00163121 -56 69183.1 0 0 0 0.99054 0.00169056 0.00169056 -57 71614.3 0 0 0 0.990202 0.00175215 0.00175215 -58 74131 0 0 0 0.989853 0.00181605 0.00181605 -59 76736.1 0 0 0 0.989491 0.00188237 0.00188237 -60 79432.8 0 0 0 0.989116 0.0019512 0.0019512 -61 82224.3 0 0 0 0.988727 0.00202264 0.00202264 -62 85113.8 0 0 0 0.988324 0.0020968 0.0020968 -63 88104.9 0 0 0 0.987906 0.00217379 0.00217379 -64 91201.1 0 0 0 0.987473 0.00225372 0.00225372 -65 94406.1 0 0 0 0.987025 0.00233671 0.00233671 -66 97723.7 0 0 0 0.98656 0.00242289 0.00242289 -67 101158 0 0 0 0.986078 0.00251239 0.00251239 -68 104713 0 0 0 0.985579 0.00260534 0.00260534 -69 108393 0 0 0 0.985061 0.0027019 0.0027019 -70 112202 0 0 0 0.984524 0.00280222 0.00280222 -71 116145 0 0 0 0.983967 0.00290644 0.00290644 -72 120226 0 0 0 0.98339 0.00301474 0.00301474 -73 124451 0 0 0 0.982792 0.00312728 0.00312728 -74 128825 0 0 0 0.982172 0.00324425 0.00324425 -75 133352 0 0 0 0.981528 0.00336584 0.00336584 -76 138038 0 0 0 0.980861 0.00349225 0.00349225 -77 142889 0 0 0 0.980169 0.00362368 0.00362368 -78 147911 0 0 0 0.979452 0.00376035 0.00376035 -79 153109 0 0 0 0.978707 0.00390249 0.00390249 -80 158489 0 0 0 0.977935 0.00405035 0.00405035 -81 164059 0 0 0 0.977134 0.00420416 0.00420416 -82 169824 0 0 0 0.976303 0.00436419 0.00436419 -83 175792 0 0 0 0.975441 0.00453073 0.00453073 -84 181970 0 0 0 0.974546 0.00470405 0.00470405 -85 188365 0 0 0 0.973617 0.00488447 0.00488447 -86 194984 0 0 0 0.972653 0.0050723 0.0050723 -87 201837 0 0 0 0.971653 0.00526787 0.00526787 -88 208930 0 0 0 0.970614 0.00547155 0.00547155 -89 216272 0 0 0 0.969536 0.0056837 0.0056837 -90 223872 0 0 0 0.968416 0.0059047 0.0059047 -91 231739 0 0 0 0.967253 0.00613497 0.00613497 -92 239883 0 0 0 0.966046 0.00637494 0.00637494 -93 248313 0 0 0 0.964792 0.00662504 0.00662504 -94 257040 0 0 0 0.963489 0.00688576 0.00688576 -95 266073 0 0 0 0.962135 0.0071576 0.0071576 -96 275423 0 0 0 0.960728 0.00744107 0.00744107 -97 285102 0 0 0 0.959266 0.00773672 0.00773672 -98 295121 0 0 0 0.957746 0.00804512 0.00804512 -99 305492 0 0 0 0.956167 0.00836689 0.00836689 -100 316228 0 0 0 0.954524 0.00870266 0.00870266 -101 327341 0 0 0 0.952815 0.00905309 0.00905309 -102 338844 0 0 0 0.951038 0.0094189 0.0094189 -103 350752 0 0 0 0.94919 0.0098008 0.0098008 -104 363078 0 0 0 0.947266 0.0101996 0.0101996 -105 375837 0 0 0 0.945264 0.010616 0.010616 -106 389045 0 0 0 0.94318 0.0110511 0.0110511 -107 402717 0 0 0 0.941009 0.0115055 0.0115055 -108 416869 0 0 0 0.938749 0.0119803 0.0119803 -109 431519 0 0 0 0.936395 0.0124765 0.0124765 -110 446684 0 0 0 0.933941 0.012995 0.012995 -111 462381 0 0 0 0.931384 0.013537 0.013537 -112 478630 0 0 0 0.928718 0.0141036 0.0141036 -113 495450 0 0 0 0.925938 0.0146959 0.0146959 -114 512861 0 0 0 0.923037 0.0153152 0.0153152 -115 530884 0 0 0 0.92001 0.0159628 0.0159628 -116 549541 0 0 0 0.916851 0.01664 0.01664 -117 568853 0 0 0 0.913552 0.0173482 0.0173482 -118 588844 0 0 0 0.910106 0.0180889 0.0180889 -119 609537 0 0 0 0.906505 0.0188635 0.0188635 -120 630957 0 0 0 0.90274 0.0196737 0.0196737 -121 653131 0 0 0 0.898804 0.020521 0.020521 -122 676083 0 0 0 0.894685 0.0214071 0.0214071 -123 699842 0 0 0 0.890374 0.0223338 0.0223338 -124 724436 0 0 0 0.88586 0.0233028 0.0233028 -125 749894 0 0 0 0.881132 0.024316 0.024316 -126 776247 0 0 0 0.876175 0.0253752 0.0253752 -127 803526 0 0 0 0.870979 0.0264823 0.0264823 -128 831764 0 0 0 0.865526 0.0276392 0.0276392 -129 860994 0 0 0 0.859803 0.0288478 0.0288478 -130 891251 0 0 0 0.853793 0.0301103 0.0301103 -131 922571 0 0 0 0.847477 0.0314283 0.0314283 -132 954993 0 0 0 0.840837 0.032804 0.032804 -133 988553 0 0 0 0.833852 0.0342392 0.0342392 -134 1.02329e+06 0 0 0 0.8265 0.0357358 0.0357358 -135 1.05925e+06 0 0 0 0.818757 0.0372954 0.0372954 -136 1.09648e+06 0 0 0 0.810597 0.0389197 0.0389197 -137 1.13501e+06 0 0 0 0.801993 0.0406103 0.0406103 -138 1.1749e+06 0 0 0 0.792916 0.0423684 0.0423684 -139 1.21619e+06 0 0 0 0.783333 0.0441951 0.0441951 -140 1.25893e+06 0 0 0 0.773212 0.0460913 0.0460913 -141 1.30317e+06 0 0 0 0.762516 0.0480575 0.0480575 -142 1.34896e+06 0 0 0 0.751205 0.0500938 0.0500938 -143 1.39637e+06 0 0 0 0.73924 0.0521999 0.0521999 -144 1.44544e+06 0 0 0 0.726576 0.054375 0.054375 -145 1.49624e+06 0 0 0 0.713166 0.0566176 0.0566176 -146 1.54882e+06 0 0 0 0.698962 0.0589256 0.0589256 -147 1.60325e+06 0 0 0 0.683913 0.061296 0.061296 -148 1.65959e+06 0 0 0 0.667965 0.0637251 0.0637251 -149 1.71791e+06 0 0 0 0.651063 0.0662081 0.0662081 -150 1.77828e+06 0 0 0 0.633151 0.0687389 0.0687389 -151 1.84077e+06 0 0 0 0.614172 0.0713105 0.0713105 -152 1.90546e+06 0 0 0 0.59407 0.0739142 0.0739142 -153 1.97242e+06 0 0 0 0.572793 0.0765402 0.0765402 -154 2.04174e+06 0 0 0 0.550289 0.0791766 0.0791766 -155 2.11349e+06 0 0 0 0.526518 0.0818104 0.0818104 -156 2.18776e+06 0 0 0 0.501446 0.0844263 0.0844263 -157 2.26464e+06 0 0 0 0.475054 0.0870076 0.0870076 -158 2.34423e+06 0 0 0 0.447343 0.0895357 0.0895357 -159 2.42661e+06 0 0 0 0.418338 0.0919901 0.0919901 -160 2.51189e+06 0 0 0 0.388099 0.0943491 0.0943491 -161 2.60016e+06 0 0 0 0.356724 0.0965895 0.0965895 -162 2.69153e+06 0 0 0 0.324366 0.0986873 0.0986873 -163 2.78612e+06 0 0 0 0.291238 0.100618 0.100618 -164 2.88403e+06 0 0 0 0.257627 0.102358 0.102358 -165 2.98538e+06 0 0 0 0.223905 0.103885 0.103885 -166 3.0903e+06 0 0 0 0.190535 0.105178 0.105178 -167 3.1989e+06 0 0 0 0.158074 0.106221 0.106221 -168 3.31131e+06 0 0 0 0.127163 0.107003 0.107003 -169 3.42768e+06 0 0 0 0.0985062 0.107517 0.107517 -170 3.54813e+06 0 0 0 0.0728213 0.107766 0.107766 -171 3.67282e+06 0 0 0 0.0507711 0.107757 0.107757 -172 3.80189e+06 0 0 0 0.0328623 0.107503 0.107503 -173 3.9355e+06 0 0 0 0.0193301 0.107022 0.107022 -174 4.0738e+06 0 0 0 0.0100337 0.106328 0.106328 -175 4.21697e+06 0 0 0 0.00441127 0.10543 0.10543 -176 4.36516e+06 0 0 0 0.00155155 0.104325 0.104325 -177 4.51856e+06 0 0 0 0.000406273 0.102995 0.102995 -178 4.67735e+06 0 0 0 7.73392e-05 0.101404 0.101404 -179 4.84172e+06 0 0 0 1.66218e-05 0.0994966 0.0994966 -180 5.01187e+06 0 0 0 8.46021e-06 0.0971968 0.0971968 -181 5.188e+06 0 0 0 5.88048e-06 0.0943984 0.0943984 -182 5.37032e+06 0 0 0 4.0913e-06 0.0909445 0.0909445 -183 5.55904e+06 0 0 0 2.72291e-06 0.0865866 0.0865866 -184 5.7544e+06 0 0 0 1.66087e-06 0.0808876 0.0808876 -185 5.95662e+06 0 0 0 8.62428e-07 0.0729707 0.0729707 -186 6.16595e+06 0 0 0 3.32173e-07 0.0615168 0.0615168 -187 6.38263e+06 0 0 0 1.25725e-07 0.0564797 0.0564797 -188 6.60693e+06 0 0 0 2.23414e-07 0.0571122 0.0571122 -189 6.83912e+06 0 0 0 8.0776e-07 0.0565696 0.0565696 -190 7.07946e+06 0 0 0 2.10193e-06 0.0548885 0.0548885 -191 7.32825e+06 0 0 0 4.40761e-06 0.05236 0.05236 -192 7.58578e+06 0 0 0 8.22327e-06 0.049135 0.049135 -193 7.85236e+06 0 0 0 1.43941e-05 0.0452795 0.0452795 -194 8.12831e+06 0 0 0 2.44653e-05 0.0408147 0.0408147 -195 8.41395e+06 0 0 0 4.23052e-05 0.0357352 0.0357352 -196 8.70964e+06 0 0 0 8.80122e-05 0.0300185 0.0300185 -197 9.01571e+06 0 0 0 0.000261816 0.0236299 0.0236299 -198 9.33254e+06 0 0 0 0.000864497 0.0165277 0.0165277 -199 9.66051e+06 0 0 0 0.0024616 0.00866633 0.00866633 +0 1e+04 0 0 0 0.999 0.000245 0.000245 +1 1.04e+04 0 0 0 0.999 0.000245 0.000245 +2 1.07e+04 0 0 0 0.999 0.000254 0.000254 +3 1.11e+04 0 0 0 0.999 0.000263 0.000263 +4 1.15e+04 0 0 0 0.998 0.000272 0.000272 +5 1.19e+04 0 0 0 0.998 0.000282 0.000282 +6 1.23e+04 0 0 0 0.998 0.000292 0.000292 +7 1.27e+04 0 0 0 0.998 0.000302 0.000302 +8 1.32e+04 0 0 0 0.998 0.000313 0.000313 +9 1.36e+04 0 0 0 0.998 0.000324 0.000324 +10 1.41e+04 0 0 0 0.998 0.000335 0.000335 +11 1.46e+04 0 0 0 0.998 0.000347 0.000347 +12 1.51e+04 0 0 0 0.998 0.00036 0.00036 +13 1.57e+04 0 0 0 0.998 0.000372 0.000372 +14 1.62e+04 0 0 0 0.998 0.000385 0.000385 +15 1.68e+04 0 0 0 0.998 0.000399 0.000399 +16 1.74e+04 0 0 0 0.998 0.000413 0.000413 +17 1.8e+04 0 0 0 0.998 0.000428 0.000428 +18 1.86e+04 0 0 0 0.997 0.000443 0.000443 +19 1.93e+04 0 0 0 0.997 0.000459 0.000459 +20 2e+04 0 0 0 0.997 0.000475 0.000475 +21 2.07e+04 0 0 0 0.997 0.000492 0.000492 +22 2.14e+04 0 0 0 0.997 0.00051 0.00051 +23 2.21e+04 0 0 0 0.997 0.000528 0.000528 +24 2.29e+04 0 0 0 0.997 0.000547 0.000547 +25 2.37e+04 0 0 0 0.997 0.000566 0.000566 +26 2.45e+04 0 0 0 0.997 0.000586 0.000586 +27 2.54e+04 0 0 0 0.997 0.000607 0.000607 +28 2.63e+04 0 0 0 0.996 0.000629 0.000629 +29 2.72e+04 0 0 0 0.996 0.000651 0.000651 +30 2.82e+04 0 0 0 0.996 0.000674 0.000674 +31 2.92e+04 0 0 0 0.996 0.000698 0.000698 +32 3.02e+04 0 0 0 0.996 0.000723 0.000723 +33 3.13e+04 0 0 0 0.996 0.000749 0.000749 +34 3.24e+04 0 0 0 0.996 0.000776 0.000776 +35 3.35e+04 0 0 0 0.995 0.000804 0.000804 +36 3.47e+04 0 0 0 0.995 0.000832 0.000832 +37 3.59e+04 0 0 0 0.995 0.000862 0.000862 +38 3.72e+04 0 0 0 0.995 0.000893 0.000893 +39 3.85e+04 0 0 0 0.995 0.000925 0.000925 +40 3.98e+04 0 0 0 0.995 0.000958 0.000958 +41 4.12e+04 0 0 0 0.994 0.000993 0.000993 +42 4.27e+04 0 0 0 0.994 0.00103 0.00103 +43 4.42e+04 0 0 0 0.994 0.00107 0.00107 +44 4.57e+04 0 0 0 0.994 0.0011 0.0011 +45 4.73e+04 0 0 0 0.994 0.00114 0.00114 +46 4.9e+04 0 0 0 0.993 0.00118 0.00118 +47 5.07e+04 0 0 0 0.993 0.00123 0.00123 +48 5.25e+04 0 0 0 0.993 0.00127 0.00127 +49 5.43e+04 0 0 0 0.993 0.00132 0.00132 +50 5.62e+04 0 0 0 0.992 0.00137 0.00137 +51 5.82e+04 0 0 0 0.992 0.00141 0.00141 +52 6.03e+04 0 0 0 0.992 0.00147 0.00147 +53 6.24e+04 0 0 0 0.991 0.00152 0.00152 +54 6.46e+04 0 0 0 0.991 0.00157 0.00157 +55 6.68e+04 0 0 0 0.991 0.00163 0.00163 +56 6.92e+04 0 0 0 0.991 0.00169 0.00169 +57 7.16e+04 0 0 0 0.99 0.00175 0.00175 +58 7.41e+04 0 0 0 0.99 0.00182 0.00182 +59 7.67e+04 0 0 0 0.989 0.00188 0.00188 +60 7.94e+04 0 0 0 0.989 0.00195 0.00195 +61 8.22e+04 0 0 0 0.989 0.00202 0.00202 +62 8.51e+04 0 0 0 0.988 0.0021 0.0021 +63 8.81e+04 0 0 0 0.988 0.00217 0.00217 +64 9.12e+04 0 0 0 0.987 0.00225 0.00225 +65 9.44e+04 0 0 0 0.987 0.00234 0.00234 +66 9.77e+04 0 0 0 0.987 0.00242 0.00242 +67 1.01e+05 0 0 0 0.986 0.00251 0.00251 +68 1.05e+05 0 0 0 0.986 0.00261 0.00261 +69 1.08e+05 0 0 0 0.985 0.0027 0.0027 +70 1.12e+05 0 0 0 0.985 0.0028 0.0028 +71 1.16e+05 0 0 0 0.984 0.00291 0.00291 +72 1.2e+05 0 0 0 0.983 0.00301 0.00301 +73 1.24e+05 0 0 0 0.983 0.00313 0.00313 +74 1.29e+05 0 0 0 0.982 0.00324 0.00324 +75 1.33e+05 0 0 0 0.982 0.00337 0.00337 +76 1.38e+05 0 0 0 0.981 0.00349 0.00349 +77 1.43e+05 0 0 0 0.98 0.00362 0.00362 +78 1.48e+05 0 0 0 0.979 0.00376 0.00376 +79 1.53e+05 0 0 0 0.979 0.0039 0.0039 +80 1.58e+05 0 0 0 0.978 0.00405 0.00405 +81 1.64e+05 0 0 0 0.977 0.0042 0.0042 +82 1.7e+05 0 0 0 0.976 0.00436 0.00436 +83 1.76e+05 0 0 0 0.975 0.00453 0.00453 +84 1.82e+05 0 0 0 0.975 0.0047 0.0047 +85 1.88e+05 0 0 0 0.974 0.00488 0.00488 +86 1.95e+05 0 0 0 0.973 0.00507 0.00507 +87 2.02e+05 0 0 0 0.972 0.00527 0.00527 +88 2.09e+05 0 0 0 0.971 0.00547 0.00547 +89 2.16e+05 0 0 0 0.97 0.00568 0.00568 +90 2.24e+05 0 0 0 0.968 0.0059 0.0059 +91 2.32e+05 0 0 0 0.967 0.00613 0.00613 +92 2.4e+05 0 0 0 0.966 0.00637 0.00637 +93 2.48e+05 0 0 0 0.965 0.00663 0.00663 +94 2.57e+05 0 0 0 0.963 0.00689 0.00689 +95 2.66e+05 0 0 0 0.962 0.00716 0.00716 +96 2.75e+05 0 0 0 0.961 0.00744 0.00744 +97 2.85e+05 0 0 0 0.959 0.00774 0.00774 +98 2.95e+05 0 0 0 0.958 0.00805 0.00805 +99 3.05e+05 0 0 0 0.956 0.00837 0.00837 +100 3.16e+05 0 0 0 0.955 0.0087 0.0087 +101 3.27e+05 0 0 0 0.953 0.00905 0.00905 +102 3.39e+05 0 0 0 0.951 0.00942 0.00942 +103 3.51e+05 0 0 0 0.949 0.0098 0.0098 +104 3.63e+05 0 0 0 0.947 0.0102 0.0102 +105 3.76e+05 0 0 0 0.945 0.0106 0.0106 +106 3.89e+05 0 0 0 0.943 0.0111 0.0111 +107 4.03e+05 0 0 0 0.941 0.0115 0.0115 +108 4.17e+05 0 0 0 0.939 0.012 0.012 +109 4.32e+05 0 0 0 0.936 0.0125 0.0125 +110 4.47e+05 0 0 0 0.934 0.013 0.013 +111 4.62e+05 0 0 0 0.931 0.0135 0.0135 +112 4.79e+05 0 0 0 0.929 0.0141 0.0141 +113 4.95e+05 0 0 0 0.926 0.0147 0.0147 +114 5.13e+05 0 0 0 0.923 0.0153 0.0153 +115 5.31e+05 0 0 0 0.92 0.016 0.016 +116 5.5e+05 0 0 0 0.917 0.0166 0.0166 +117 5.69e+05 0 0 0 0.914 0.0173 0.0173 +118 5.89e+05 0 0 0 0.91 0.0181 0.0181 +119 6.1e+05 0 0 0 0.907 0.0189 0.0189 +120 6.31e+05 0 0 0 0.903 0.0197 0.0197 +121 6.53e+05 0 0 0 0.899 0.0205 0.0205 +122 6.76e+05 0 0 0 0.895 0.0214 0.0214 +123 7e+05 0 0 0 0.89 0.0223 0.0223 +124 7.24e+05 0 0 0 0.886 0.0233 0.0233 +125 7.5e+05 0 0 0 0.881 0.0243 0.0243 +126 7.76e+05 0 0 0 0.876 0.0254 0.0254 +127 8.04e+05 0 0 0 0.871 0.0265 0.0265 +128 8.32e+05 0 0 0 0.866 0.0276 0.0276 +129 8.61e+05 0 0 0 0.86 0.0288 0.0288 +130 8.91e+05 0 0 0 0.854 0.0301 0.0301 +131 9.23e+05 0 0 0 0.847 0.0314 0.0314 +132 9.55e+05 0 0 0 0.841 0.0328 0.0328 +133 9.89e+05 0 0 0 0.834 0.0342 0.0342 +134 1.02e+06 0 0 0 0.827 0.0357 0.0357 +135 1.06e+06 0 0 0 0.819 0.0373 0.0373 +136 1.1e+06 0 0 0 0.811 0.0389 0.0389 +137 1.14e+06 0 0 0 0.802 0.0406 0.0406 +138 1.17e+06 0 0 0 0.793 0.0424 0.0424 +139 1.22e+06 0 0 0 0.783 0.0442 0.0442 +140 1.26e+06 0 0 0 0.773 0.0461 0.0461 +141 1.3e+06 0 0 0 0.763 0.0481 0.0481 +142 1.35e+06 0 0 0 0.751 0.0501 0.0501 +143 1.4e+06 0 0 0 0.739 0.0522 0.0522 +144 1.45e+06 0 0 0 0.727 0.0544 0.0544 +145 1.5e+06 0 0 0 0.713 0.0566 0.0566 +146 1.55e+06 0 0 0 0.699 0.0589 0.0589 +147 1.6e+06 0 0 0 0.684 0.0613 0.0613 +148 1.66e+06 0 0 0 0.668 0.0637 0.0637 +149 1.72e+06 0 0 0 0.651 0.0662 0.0662 +150 1.78e+06 0 0 0 0.633 0.0687 0.0687 +151 1.84e+06 0 0 0 0.614 0.0713 0.0713 +152 1.91e+06 0 0 0 0.594 0.0739 0.0739 +153 1.97e+06 0 0 0 0.573 0.0765 0.0765 +154 2.04e+06 0 0 0 0.55 0.0792 0.0792 +155 2.11e+06 0 0 0 0.527 0.0818 0.0818 +156 2.19e+06 0 0 0 0.501 0.0844 0.0844 +157 2.26e+06 0 0 0 0.475 0.087 0.087 +158 2.34e+06 0 0 0 0.447 0.0895 0.0895 +159 2.43e+06 0 0 0 0.418 0.092 0.092 +160 2.51e+06 0 0 0 0.388 0.0943 0.0943 +161 2.6e+06 0 0 0 0.357 0.0966 0.0966 +162 2.69e+06 0 0 0 0.324 0.0987 0.0987 +163 2.79e+06 0 0 0 0.291 0.101 0.101 +164 2.88e+06 0 0 0 0.258 0.102 0.102 +165 2.99e+06 0 0 0 0.224 0.104 0.104 +166 3.09e+06 0 0 0 0.191 0.105 0.105 +167 3.2e+06 0 0 0 0.158 0.106 0.106 +168 3.31e+06 0 0 0 0.127 0.107 0.107 +169 3.43e+06 0 0 0 0.0985 0.108 0.108 +170 3.55e+06 0 0 0 0.0728 0.108 0.108 +171 3.67e+06 0 0 0 0.0508 0.108 0.108 +172 3.8e+06 0 0 0 0.0329 0.108 0.108 +173 3.94e+06 0 0 0 0.0193 0.107 0.107 +174 4.07e+06 0 0 0 0.01 0.106 0.106 +175 4.22e+06 0 0 0 0.00441 0.105 0.105 +176 4.37e+06 0 0 0 0.00155 0.104 0.104 +177 4.52e+06 0 0 0 0.000406 0.103 0.103 +178 4.68e+06 0 0 0 7.73e-05 0.101 0.101 +179 4.84e+06 0 0 0 1.66e-05 0.0995 0.0995 +180 5.01e+06 0 0 0 8.46e-06 0.0972 0.0972 +181 5.19e+06 0 0 0 5.88e-06 0.0944 0.0944 +182 5.37e+06 0 0 0 4.09e-06 0.0909 0.0909 +183 5.56e+06 0 0 0 2.72e-06 0.0866 0.0866 +184 5.75e+06 0 0 0 1.66e-06 0.0809 0.0809 +185 5.96e+06 0 0 0 8.62e-07 0.073 0.073 +186 6.17e+06 0 0 0 3.32e-07 0.0615 0.0615 +187 6.38e+06 0 0 0 1.26e-07 0.0565 0.0565 +188 6.61e+06 0 0 0 2.23e-07 0.0571 0.0571 +189 6.84e+06 0 0 0 8.08e-07 0.0566 0.0566 +190 7.08e+06 0 0 0 2.1e-06 0.0549 0.0549 +191 7.33e+06 0 0 0 4.41e-06 0.0524 0.0524 +192 7.59e+06 0 0 0 8.22e-06 0.0491 0.0491 +193 7.85e+06 0 0 0 1.44e-05 0.0453 0.0453 +194 8.13e+06 0 0 0 2.45e-05 0.0408 0.0408 +195 8.41e+06 0 0 0 4.23e-05 0.0357 0.0357 +196 8.71e+06 0 0 0 8.8e-05 0.03 0.03 +197 9.02e+06 0 0 0 0.000262 0.0236 0.0236 +198 9.33e+06 0 0 0 0.000864 0.0165 0.0165 +199 9.66e+06 0 0 0 0.00246 0.00867 0.00867 diff --git a/test/glashow_resonance.test.cpp b/test/glashow_resonance.test.cpp index 71d01e8a..6fb67091 100644 --- a/test/glashow_resonance.test.cpp +++ b/test/glashow_resonance.test.cpp @@ -90,6 +90,9 @@ int main (int argc, char const *argv[]) squid.EvolveState(); std::cout << "Propagation done" << std::endl; std::ostream &output = std::cout; + // One part in a thousand is a somewhat loose tolerance, but this is necessary because the key area + // around the resonance has such a suppression of the (nuebar) flux that it has much higher relative error. + output.precision(3); for (auto i=0; i < num_steps-1; i++) { double ee = E_range[i]; double de = (E_range[i+1]-E_range[i])/units.GeV;