diff --git a/t/ME_data/ME_h_mt.dat b/t/ME_data/ME_h_mt.dat index ba23656..349d686 100644 --- a/t/ME_data/ME_h_mt.dat +++ b/t/ME_data/ME_h_mt.dat @@ -1,807 +1,807 @@ 3.999227056e-08 4.193300057e-05 3.89321526e-05 0.0168082963 4.057962883e-06 4.063245819e-07 6.411897603e-07 7.822040516e-07 6.064514954e-07 -2.202805699e-05 +2.338845036e-05 4.738208766e-05 1.624045121e-06 -0.0005217822994 +0.0005243567971 1.536443326e-05 9.5908372e-06 0.0339113482 8.152325116e-07 608.615459 5.962078524e-06 -0.01340321821 +0.01352350209 1.931630897 0.0005720381031 0.0008001259943 3.300880891e-05 -0.0003564506816 +0.0007712160505 0.0007037298905 0.03792377651 -6.513991407e-07 +1.079417185e-06 2.011798913 8.842231663e-08 0.0463578742 0.1064067277 9.357584012e-07 -0.006517735511 +0.007615558468 0.02272391876 3.606711294e-05 6.725111513e-07 0.005897687516 7.581255699e-05 0.005354681787 0.0006088947447 0.0083356406 0.003352568463 -3.23438533e-05 +3.95296533e-05 0.0001515363281 0.005531947882 -0.03875973775 +0.0388505972 2.312735565e-07 5.479259483e-05 8.480458133 -0.0001854378187 -4.957725793e-06 +0.0001863702054 +7.349046553e-06 5.091269317e-06 -0.0001308303481 -0.1918392997 +0.0001452556797 +0.2136081233 0.0238989834 4.446354691 -0.0002597751966 +0.0002720132705 1.267033388 0.005050109129 -2.130648865e-06 -0.05854510061 +2.848397879e-06 +0.07116846744 0.00280707581 -6.289945688 +7.10418958 7.216599144e-08 4.779195357e-06 -1.655604389e-05 +1.924410023e-05 1.522656761e-05 1.116772322e-06 0.0009539322219 9.279967349e-07 2.28564472 0.001331277201 0.001585734932 26972.25526 -8.555976389e-07 +1.2516302e-06 0.0006894623725 1.065245235e-07 0.01968534835 1.951627129e-06 4.401131811e-05 0.02561966521 49705.99132 5.383396351e-06 1.07104069e-06 2.808652999e-06 0.0001066898908 0.0002239236027 7.851273644e-08 0.002704063747 -0.0005663491795 +0.0005842786082 0.07550986067 -0.00128549481 -3.590772606e-05 +0.001335627914 +6.766626128e-05 0.1149401044 0.183403668 0.000113873958 -1.730096598e-05 -0.02430322741 +1.844065382e-05 +0.02689433062 0.001076022562 5.066109523e-06 -0.1068111238 +0.1097193795 9.253312815e-05 -91.21436585 +92.12237679 0.01945125628 8.958261209e-05 -1.678018826e-06 -0.0002943803381 +3.849412514e-06 +0.0003305165961 0.0001299537598 1.170328977e-06 0.0001328201344 3034.923139 1.378838481 -0.08550576557 -36.1778481 +0.08677679017 +36.20174579 1.785612005e-06 -9.054761711e-06 +1.032516793e-05 0.001615788007 1.23208212 0.0003642605012 2.708222463e-06 4.719890307e-06 0.001710230342 0.3182780512 -0.0008830648263 -0.08778053115 +0.0009491003377 +0.09165888221 0.2629019547 0.02748771871 -1.833687248e-06 -0.002690141449 +3.179028152e-06 +0.002811209731 1.89343446e-06 1.77338727e-06 2.291261955e-05 1.31623411e-05 0.02865008699 -0.0001316279098 +0.0001547359261 5.992647136 427898745.7 0.0002908564713 180.1583563 1.149638817 -4.549308837e-06 +4.915348327e-06 0.02938178265 12.41064482 0.0005867842208 1.235864829e-06 -4.615116294e-05 +6.021062978e-05 0.7228122449 -0.0003119308893 +0.0003703954772 1.100546805e-07 0.0007569561799 3.456980613e-06 0.001902875586 1.766220045 0.0006636538355 435.7550934 1.468265667e-05 -0.003868037348 +0.003895447939 0.005375424637 9.895862027e-06 0.002634610055 0.04246482577 1.811370936e-05 0.0008533026641 -1.935841576 -1.965419689e-06 +1.973172924 +2.216040365e-06 0.0001940609469 6.432993714e-05 436.2325667 8.585463564e-05 0.003171908449 3.1453383e-07 0.10708771 1.134414955e-07 5.434578249e-06 6.681002905e-06 -0.00246782063 -0.00431428184 +0.002957334168 +0.005022320071 2.421882039e-05 3.783451257e-08 5.144794528e-06 -8.550628039e-05 +0.0001067748459 13660.98467 117.1131418 387.5973894 -0.2157832996 +0.2286644594 2.335478847e-06 0.0001022288321 0.004686932429 0.4278874185 -0.000106415896 +0.0001089502807 1.985873409e-06 2.25291051e-07 2.945623146e-05 0.3581923626 -6.785775177e-06 +8.204917967e-06 0.004893941209 0.0005325322641 1.060152245e-07 1.118231806e-05 2.203633348e-06 -3.969211813e-06 +4.953456305e-06 0.007400958058 5.769269726e-05 1.899358848e-06 0.002291707993 0.007504512648 0.0946266303 -0.07641638751 -0.001194744559 +0.07677680285 +0.001406412108 0.000131740617 3.126413693e-07 6.676941465e-07 2.247981499e-07 1.394670123e-07 0.0003117083998 0.003124792631 3.196614452e-07 0.0005145950558 1.196095567 1.715539058e-05 2.527890627e-05 -2.665407449e-06 +4.198681685e-06 4.483788358e-07 3.597134253e-05 6.514469651e-07 0.002675227421 0.1700338633 0.0151159864 -1.877463443e-05 -0.001018686854 -0.0001769414038 +2.549964327e-05 +0.00103462801 +0.0002092909388 3.068873861e-06 1176.630703 9.378378651e-06 7.10431401e-07 4.541139018e-05 3.954117424e-06 0.003271203855 0.01288628186 0.01758893926 0.0001040458517 0.08587993013 7.493297002e-07 0.0303353991 2.525780749e-07 1.075360445e-06 0.0005909996117 3.980399911e-06 67.71892274 -8.802422373e-06 +8.935469327e-06 3.997018327e-05 -2.690317449e-06 +3.48516788e-06 0.1177580565 1.362864165e-07 0.00719011227 3.379854611e-06 1.335802906e-08 0.0008096413047 8.017608649e-08 5.475646202e-07 9.89836826e-07 0.0001097524935 0.005006547178 0.01328382633 3.708077203e-10 1.7738466e-09 6.195338978e-10 1.375292134e-08 -7.712749689e-06 +1.198886719e-05 1.087672399e-10 7.612820481e-06 1.631362221e-06 7.195255533e-08 8.307139358e-09 2.261143452e-08 -2.267767047e-07 +3.176205908e-07 7.680504754e-08 1.175875499e-07 0.0006018013011 1.467361257e-05 1.350450721e-07 -2.569869939e-07 -7.718323457e-07 +4.010394331e-07 +7.773787384e-07 0.0170653268 0.0001792376906 3.072166197e-06 22266.55115 0.06729772649 -1.065337157 +1.156967792 7.103072447e-07 0.0008701526857 0.3457109502 0.001639388444 6.829709377e-07 103.0220399 1.724918219e-06 0.006442680294 0.0002924073612 0.0061407479 3.014208435e-07 -1.065501509e-05 +1.067129219e-05 5.819752501e-07 0.002933183567 1.390851143e-08 1.359009867e-06 0.02819012119 1.30006471e-09 78.00480652 3.27001628e-06 2.84175006e-07 -3.993660531e-07 +4.457488422e-07 0.0003605232041 1.871558115e-06 2.854501301e-07 0.02600693272 0.0004600409654 -0.001348216483 +0.001476738524 2.488693144e-08 2.395438283e-11 4.248528123e-07 1.389769895e-05 8.049135616e-07 2.12824677e-07 8.405933433e-07 2.36621812e-06 0.06259159217 1.831343487e-07 2.1663617e-08 6.323481085e-10 1.339773875e-07 4.883049792e-06 1.705369167e-05 0.0001013417062 0.01536206656 1.67093181 5.039245627e-10 1.141488157e-05 3.188719492e-07 3.24135522e-10 0.0002434524475 1.383868594e-06 0.498790572 9.03612161e-05 0.02764416817 0.01341805272 7.529041731e-07 0.0002257161473 0.0003717671499 1.858022341e-08 5.715357136e-07 2.477427879e-09 0.0001517676839 2.326424939e-11 1.729565889e-05 7.354412281e-09 -6.877321994e-06 +9.505701399e-06 0.0005289104709 0.02479967619 6.212490017e-06 1.714682348e-08 7.916835703e-09 0.00934088881 4.083237353e-07 0.004696390685 6.761095529e-10 0.003004586166 6.030026635e-05 8.803060787e-08 6.957661086e-08 9.612925872e-08 1.628346787e-05 5.143053469e-07 8.784807581e-08 1.024493405e-08 3.137825282e-07 42.00569276 0.004640201019 1.146513701e-05 -4.057040434e-08 +5.876691844e-08 3.381428808e-06 2.603443408e-07 -0.002272672712 +0.002350385674 1.276808422e-05 0.06226820245 1.778585408e-08 2.378758757e-07 9.7714181e-05 0.387678215 4.398556239e-07 3.0060972e-08 3.874778496e-05 -9.119354143e-06 +1.0694862e-05 1.402325829e-06 4.695875643e-06 2.383017224e-07 2.417306793e-06 0.0002708928034 0.01384277314 1.710909166e-07 1.217157554 177.0639344 0.03861638077 7.433640647e-05 0.0002263704383 0.002948183736 0.008773163623 5.353914571e-08 360.0073872 6.478760385e-07 2.750359008e-07 8.589493477e-08 4.263303371e-07 2.005861581e-05 5.503740271e-08 0.00180575553 16.03887797 7.922076143 5.036869078e-12 1.277093114e-06 0.003373959327 1.743675895e-05 1.778007998e-06 0.0002212911594 9.259789313e-05 8.539690283e-05 2.423945599e-06 0.0001233083611 -6.913916933e-09 +9.174270377e-09 2.570708518e-07 0.0008554774811 -0.2986161102 +0.3056275035 0.0006441353055 1.203219038e-07 4.75874322e-07 258.1003097 2.768523645e-09 1.270156521e-11 1.46097374e-09 34.87523283 8.506779165e-05 4.671813569e-06 1.907838682e-05 3.716521365e-07 1.888947852e-08 -0.02669197841 +0.02682389012 15.26083816 2.873162527e-10 1.360352114e-06 9.527755241e-07 0.0003084354836 0.0009838761157 1.302357534e-05 0.004623681617 0.01694286548 9.642283586e-07 8.499293031e-06 8.770386824e-06 0.001276821989 6.668543971e-07 5.652013261e-11 0.238779506 1.624068992e-07 0.0002438443448 0.002313588467 7.53704251e-09 7.248545824e-12 0.02160343338 6.465648118e-09 7.430463219e-08 0.000611761207 2.262200983e-07 8.810901655e-07 1.688360234e-07 -2.91786064e-06 +3.452136087e-06 5.480913194e-08 0.0004601409972 3.364979337e-10 8.967694023e-09 0.1345541103 2.492529354e-06 9.220762636e-06 7.386130159e-06 3.598519539 8.045779179e-06 4.730716431e-07 2.192459188e-07 2.194928052e-06 3.130712514e-07 1.43402212e-05 0.0002608840076 23.2295584 -1.440575153e-08 +1.820768009e-08 4.161369951e-06 0.001762466634 4.255417524e-09 3320585.541 0.0009563172933 0.03635196251 0.5161165209 0.004796521016 2.31018301e-06 0.06003700358 0.0004342970655 2.009638611e-08 0.0001808310494 0.0003251077547 -4.320349927e-07 +6.445298688e-07 2.882854303e-08 0.0422146142 8.524423251e-07 19.01825412 0.000339795268 1.575704657e-07 0.0004837774983 4.132343892e-07 5.931203648e-09 0.01135109493 5.877226405e-05 1.480445352e-08 0.3379228367 1.786471685e-06 1.021098483e-07 1.079160999e-08 -9.622764911e-09 +1.26575267e-08 0.002978083967 4.588139504e-06 6.6721226e-11 3.261090243e-06 0.06059954782 4.847476637e-09 1.104977758e-06 -0.0003141612317 -0.0006548050265 +0.0003977585716 +0.0006607921183 0.000225642662 4.294099642e-06 1.185934783e-07 0.001042912926 8.03811263e-06 1.755333168e-10 6.77273477e-08 598.9038787 1.572919365e-07 0.04266693618 2.739698989e-08 0.0005135423693 1.264426626e-10 2.608778512e-07 1.857091679e-06 1.075074651e-10 1.188966384e-07 0.03379476153 0.0006473640268 3.624390296e-09 1.237883305e-05 1.016117525e-05 9.687850834e-06 0.0002738964397 9.038889191e-09 7.786947927e-07 3.198680693e-05 1.608278345e-11 0.003353487144 6.107229663e-08 0.0003768704931 0.0002005945404 3.670401818e-11 9.030155872e-09 2.042723216e-08 2.386955262e-08 4.835084799e-09 3.685817471e-11 1.066443786e-07 2.652667148e-07 3.649486136e-11 0.02477462094 3.517378902e-05 3.080027841e-09 1.018709237e-09 2.870228082e-10 8.744977958e-06 5.273214292e-09 3.168385461e-11 0.001556612382 -3.116051041e-07 +7.188695072e-07 2.492191426e-09 7.615392241e-10 0.01004217426 1.101368692e-06 3.918013244e-06 3.713827343e-09 6.857767402e-09 1.290966146e-09 6.224822208e-12 3.376319521e-06 -3.94050205e-09 +4.097980204e-09 7.633357106e-05 6.164135969e-05 1.630673977e-06 7.08423096e-09 4.080955534e-10 232770.6483 0.008068531153 8.666549381e-10 2.246457597e-10 3.222301462e-08 1.395629818e-06 2.196466632e-08 1.592802071e-08 4.326999161e-05 0.0005042839931 0.00197308444 3.032370137e-05 0.0005857579152 7.171987925e-06 6.864328965e-06 4.00109422e-06 0.005161817275 2.375668536e-09 9.778482285e-06 1.181506681e-08 1.503749258e-08 3.203838424e-11 0.3676119133 1.281881299e-09 0.0004415836065 -1.190288394e-07 +1.205682724e-07 8.77428226e-07 7.210177682e-11 1.45886715e-06 0.02121994342 0.01168544154 3.982811877e-09 8.600471911e-09 1.434258247e-08 8.098070419e-10 0.0007724068717 2.490831815e-12 2.21358682e-06 -0.2873305685 -5.652944325e-06 +0.3039319263 +1.06879986e-05 2.722304477e-09 3.122833335e-06 1.198571403e-06 4.489385847e-09 2.67381748e-08 3.471377557e-09 5.756644016e-06 6.839839041e-05 2.209846587e-07 8.291924237e-08 5.242886271e-07 6.844201737e-11 1.726048304e-06 9.225948144e-07 0.000185717301 0.002602315394 6.138200606e-11 0.0001818846094 7.733889705e-14 1.358383639e-08 1.07000594e-06 4.579302999e-07 7.849169727e-07 2.448333585 7.535064872e-11 3.779457663e-06 1.495153092e-07 3.954406596e-09 8.122956595e-09 0.0001228391098 2.269361403e-11 0.001844497418 0.01115645921 -2.08005184e-06 +2.597690032e-06 4.660852875e-11 4.199474743e-10 4.628821128e-05 0.0004980806998 7.598116953e-09 4.143456828e-13 0.002644658829 2.564668699e-08 6.50580785e-08 3.820040563e-11 0.2398554094 9.400695396e-08 0.002358652784 1.171877918e-11 1.637975849e-09 8.942959658e-08 3.408065038e-10 6.371021504e-07 8.862071381e-07 8.07800753e-13 1.275141633e-07 5.065134208e-08 3.723768129e-05 1.586955855e-06 4.630822765e-08 3.869149299e-05 0.0002373420357 8.091381421e-08 2.031825403e-10 1.493296083e-05 2.814058365e-06 4.293942142e-10 1.227720708e-09 9.668701181e-10 5.082776865e-06 0.001170366525 -2.950213679e-07 +3.153572708e-07 2.929772297e-06 3.908673483e-06 -1.641984659e-08 +1.837314879e-08 9.22544372e-06 0.004601199051 1.245696104e-11 1.623356072e-06 0.0004075079271 0.06223654452 2.229639676e-12 1.548827536e-10 0.2098995653 0.02739866019 1.922659779e-12 0.001586041269 1.845383825 0.8609993903 0.09603688052 5.904130818e-11 9.958673683e-12 0.05389921105 0.0001869127049 6.618563858e-08 3.845467031e-07 9.065990802e-07 0.0003895671306 2.782863325e-11 0.0002689153495 5.924260979e-06 0.0005747447785 4.117458599e-09 5.117012888e-08 3.396356612e-09 8.109913215e-12 4.329382871e-08 4.881990436e-07 1.276280038e-11 9.120675139e-12 7.192434752e-07 3.913876447e-11 1.9517939e-08 0.11670582 1.114874302e-06 -8.08838007e-07 +1.051244942e-06 2.596368516e-07 0.0001277437573 0.0004507490086 5.934514686e-10 9.045316154e-06 1.1680623e-08 4.405302098e-07 1.103299007e-15 7.105364025e-08 6.778113877e-05 2.587602347e-06 2.630707294e-07 2.256922483e-10 1.585065727e-05 -9.078985718e-05 +0.0001137727537 1.439932327e-09 0.003315113132 4.417287152e-05 3.208255435e-09 7.222413113e-10 0.000300699488 4.624336745e-05 8.037199117e-10 6.718759329e-12 3.83547078e-05 7.94104255e-05 4.496517789e-05 1.351102043e-06 9.177810983e-10 8.693619098e-07 1.929511565e-06 4.773961763e-09 1.896551164e-08 1.552346596e-08 9.245502051e-06 0.0001037407773 6.368845512e-07 4.041904863e-06 1.506851808e-08 9.717844454e-12 2.297581007e-05 0.0234227485 0.00306449787 5.437659639e-08 0.001404494601 1.582456182e-12 10.74537802 0.0006876957837 3.080472252e-07 4.979458387e-09 2.191738006e-08 diff --git a/t/ME_data/ME_h_mtmb.dat b/t/ME_data/ME_h_mtmb.dat index 463c576..8370f81 100644 --- a/t/ME_data/ME_h_mtmb.dat +++ b/t/ME_data/ME_h_mtmb.dat @@ -1,807 +1,807 @@ 4.043634525e-08 3.947735809e-05 3.978333671e-05 0.01707962353 4.133534858e-06 4.12758231e-07 6.527379321e-07 7.807924441e-07 6.213247187e-07 -2.158745738e-05 +2.292064051e-05 4.929604661e-05 1.648537835e-06 -0.0005294944592 +0.0005321070091 1.551221705e-05 9.732109994e-06 0.03508523471 8.385697667e-07 608.0814044 6.0749939e-06 -0.01370931025 +0.01383234108 1.96358809 0.0005777778828 0.000827073758 3.396653466e-05 -0.0003663185905 +0.0007925662403 0.000712177849 0.03945983226 -6.694529237e-07 +1.109333656e-06 2.002990871 9.107975513e-08 0.04804489153 0.1078671772 9.503786206e-07 -0.00665553606 +0.007776569624 0.02341176253 3.684250185e-05 6.861318278e-07 0.005400516499 7.647108683e-05 0.005354950581 0.0006200616798 0.008448060377 0.003448403835 -3.215759836e-05 +3.93020183e-05 0.0001546133101 0.005591549046 -0.03949683081 +0.03958941813 2.350530066e-07 5.66847762e-05 8.602097332 -0.0001911650059 -5.061383546e-06 +0.0001921261891 +7.502702823e-06 5.119802912e-06 -0.0001348170624 -0.1914012121 +0.0001496819684 +0.213120324 0.02428150541 4.539435975 -0.0002614830785 +0.0002738016111 1.314349378 0.005127584311 -2.106730988e-06 -0.06042114066 +2.816422817e-06 +0.07344901515 0.002877526755 -6.46292188 +7.299557828 7.320533172e-08 4.686272323e-06 -1.674918088e-05 +1.946859515e-05 1.561186199e-05 1.127468541e-06 0.0009956223126 9.322344484e-07 2.355779109 0.001354263217 0.001607437479 27300.17114 -8.763842677e-07 +1.282038387e-06 0.0006972261753 1.071457504e-07 0.02005686587 1.961498186e-06 4.460773766e-05 0.02601668696 50724.57307 5.480882962e-06 1.076584976e-06 2.869362531e-06 0.0001072518549 0.0002293795185 8.022114184e-08 0.002698012227 -0.0005830909647 +0.0006015504033 0.07725070316 -0.001325583365 -3.695641615e-05 +0.001377279885 +6.964246376e-05 0.1179283377 0.187466976 0.0001155356229 -1.778799047e-05 -0.02462347248 +1.895976068e-05 +0.02724871881 0.001087874101 5.154078743e-06 -0.1031621088 +0.1059710091 9.484855944e-05 -93.85073623 +94.78499143 0.02033391802 9.092704303e-05 -1.705551286e-06 -0.0002974887579 +3.912572591e-06 +0.0003340065858 0.0001322670062 1.193514887e-06 0.0001344144727 3095.568362 1.406347883 -0.08812981858 -35.91258076 +0.08943984916 +35.93630323 1.82576512e-06 -9.27160445e-06 +1.057243426e-05 0.001666287395 1.269161514 0.0003682229129 2.790042514e-06 4.841321328e-06 0.001729869633 0.3293946648 -0.0008921922127 -0.08111101834 +0.0009589102693 +0.08469469459 0.2664463909 0.02792004132 -1.869728836e-06 -0.002768079641 +3.241512757e-06 +0.002892655486 1.908915912e-06 1.7932868e-06 2.342012655e-05 1.3539045e-05 0.02918304927 -0.0001341037055 +0.0001576463616 6.099926061 438766540.3 0.0003006556126 182.0743425 1.166556788 -4.685195193e-06 +5.062168162e-06 0.02991092642 12.79082303 0.0005990054968 1.245221356e-06 -4.607261468e-05 +6.010815261e-05 0.729519575 -0.0003077751881 +0.0003654608813 1.113634813e-07 0.0007837766949 3.49282179e-06 0.001938345149 1.741954359 0.000679566605 447.8606334 1.502685811e-05 -0.003972927121 +0.004001081006 0.005572294955 1.00407234e-05 0.002697149585 0.04268640856 1.838639564e-05 0.0008683709605 -1.984375256 -2.023641764e-06 +2.022642541 +2.281686632e-06 0.000198987372 6.590604109e-05 451.0223058 8.870879716e-05 0.003243511862 3.187920803e-07 0.1107987679 1.145759171e-07 5.532565231e-06 6.903234805e-06 -0.002276208464 -0.004410664221 +0.002727714074 +0.005134520243 2.467146722e-05 3.863236208e-08 5.225855864e-06 -8.615548536e-05 +0.0001075855321 13702.18793 118.8790931 398.6816605 -0.2170834277 +0.2300421985 2.355356756e-06 0.0001050726086 0.004854583584 0.4416609149 -0.0001091431342 +0.0001117424705 2.015137689e-06 2.270813172e-07 3.065037179e-05 0.3695610208 -6.626431375e-06 +8.012249806e-06 0.005034520677 0.0005496323042 1.077013565e-07 1.134883954e-05 2.217217266e-06 -4.053576036e-06 +5.058740304e-06 0.007630107217 5.870346284e-05 1.937770986e-06 0.002308560802 0.00756714829 0.09713373161 -0.07865638607 -0.001166988788 +0.07902736629 +0.001373738971 0.0001360267898 3.146930172e-07 6.721239549e-07 2.308081382e-07 1.402445764e-07 0.000317796034 0.003250721897 3.220983676e-07 0.0005361385069 1.1674911 1.701754998e-05 2.555659626e-05 -2.72704331e-06 +4.295773543e-06 4.575225674e-07 3.723503521e-05 6.625793946e-07 0.002713331741 0.1757705927 0.01528287892 -1.914346383e-05 -0.001027883901 -0.0001819228709 +2.600058607e-05 +0.00104396898 +0.000215183149 3.169381656e-06 1197.183436 9.459688271e-06 7.151860571e-07 4.709453403e-05 4.024117529e-06 0.003399032872 0.01326938543 0.0181234889 0.0001070400499 0.08816817298 7.557106621e-07 0.03150507468 2.546159018e-07 1.083218854e-06 0.0006029839411 4.109536026e-06 69.62539127 -9.053551382e-06 +9.190394104e-06 3.976334138e-05 -2.74610264e-06 +3.557434726e-06 0.1217200425 1.39072939e-07 0.007475099982 3.40898863e-06 1.355970335e-08 0.0008211141398 8.085185969e-08 5.540708628e-07 9.793920273e-07 0.0001116678867 0.005087064408 0.01371958748 3.747341913e-10 1.784590611e-09 6.285401932e-10 1.403711889e-08 -7.942905268e-06 +1.234662607e-05 1.099705834e-10 7.694717507e-06 1.645454731e-06 7.281682073e-08 8.437713425e-09 2.293298229e-08 -2.337785459e-07 +3.274272811e-07 7.748072516e-08 1.184163566e-07 0.0006048814428 1.496340286e-05 1.358047852e-07 -2.641795688e-07 -7.894816929e-07 +4.122637606e-07 +7.951549138e-07 0.01733476065 0.0001847336643 3.09120607e-06 22646.44954 0.06825326182 -1.097908477 +1.192340601 7.194716543e-07 0.0008831460117 0.3510995791 0.001703720233 6.923341702e-07 104.5021475 1.817413974e-06 0.006274458892 0.000298133205 0.006212019474 2.971106209e-07 -1.073536231e-05 +1.075176214e-05 6.001267547e-07 0.002926676713 1.394485679e-08 1.378850246e-06 0.02873430354 1.324231029e-09 78.68865757 3.323528263e-06 2.901477846e-07 -4.025907172e-07 +4.493480222e-07 0.0003762267637 1.888509349e-06 2.984779601e-07 0.02649146415 0.0004735978501 -0.001390723358 +0.001523297471 2.511800946e-08 2.406174629e-11 4.343797764e-07 1.413682848e-05 8.187236684e-07 2.137756399e-07 8.487602643e-07 2.382894999e-06 0.06385257472 1.843038529e-07 2.200340736e-08 6.500296312e-10 1.369543749e-07 5.034654856e-06 1.729239287e-05 0.0001029880478 0.01584013323 1.717937137 5.051793703e-10 1.169719308e-05 3.274892339e-07 3.257575264e-10 0.0002494213231 1.395452311e-06 0.503263349 9.183628731e-05 0.02652180806 0.01351083505 7.697372696e-07 0.0002345916253 0.0003785854477 1.911331681e-08 5.741981368e-07 2.496303102e-09 0.0001538309849 2.342705363e-11 1.733704579e-05 7.440708121e-09 -6.97019919e-06 +9.634074462e-06 0.0005366276141 0.02581853732 6.351131637e-06 1.725302693e-08 7.963798782e-09 0.009484858796 4.13192111e-07 0.004760868685 6.868174526e-10 0.003035550995 6.122424548e-05 8.912863951e-08 7.049197773e-08 9.784536423e-08 1.644442743e-05 5.17689779e-07 8.945719428e-08 1.043648257e-08 3.220171418e-07 42.8640364 0.004721640337 1.17318728e-05 -4.163036826e-08 +6.030229416e-08 3.442022366e-06 2.613636766e-07 -0.002337930201 +0.002417874611 1.311028403e-05 0.06454239294 1.788371853e-08 2.431129133e-07 9.887586539e-05 0.393826783 4.435492073e-07 3.017594495e-08 3.944030937e-05 -8.95155309e-06 +1.049807075e-05 1.480534305e-06 4.80208438e-06 2.278546025e-07 2.431285271e-06 0.0002795248679 0.01414217835 1.721740596e-07 1.225137873 184.4942235 0.03969271858 7.5889168e-05 0.0002307099487 0.002991500607 0.008823843341 5.377454625e-08 363.5972586 6.577363411e-07 2.768216121e-07 8.808374953e-08 4.320144869e-07 2.087813482e-05 5.563234756e-08 0.001828101352 16.31099206 8.047143214 5.054980175e-12 1.31165661e-06 0.003453525346 1.775158542e-05 1.735750599e-06 0.0002258034827 9.302350911e-05 8.722704625e-05 2.463340497e-06 0.0001266405389 -7.186989683e-09 +9.536618271e-09 2.616936561e-07 0.0008671500774 -0.3034312162 +0.3105556664 0.0006552874753 1.250994614e-07 4.836707652e-07 267.0859855 2.784590138e-09 1.277933942e-11 1.625022487e-09 35.42273613 8.653074513e-05 4.787553737e-06 1.988624796e-05 3.776413822e-07 1.896902599e-08 -0.02757075064 +0.02770700524 15.7257863 2.903221186e-10 1.383656135e-06 9.796122701e-07 0.000313466982 0.001008582166 1.326553759e-05 0.004678660451 0.01717544418 9.763673978e-07 8.769786833e-06 8.887711201e-06 0.001304503338 6.699929944e-07 5.668340919e-11 0.2432408562 1.65579253e-07 0.000246842309 0.002380529019 7.590644517e-09 7.382339234e-12 0.02191318318 6.595093824e-09 7.59955463e-08 0.0006241836255 2.306006909e-07 8.949234518e-07 1.745865901e-07 -2.984765317e-06 +3.531291359e-06 5.515085903e-08 0.0004672383281 3.41289611e-10 9.072028681e-09 0.1360437848 2.516296619e-06 9.347242718e-06 7.560996071e-06 3.650323724 8.084973509e-06 4.821049366e-07 2.262423437e-07 2.237305224e-06 3.165244113e-07 1.440989373e-05 0.0002661269589 24.06039199 -1.471908477e-08 +1.860370743e-08 4.189979743e-06 0.001821016076 4.36609224e-09 3372180.49 0.0009765702499 0.03671795044 0.5266849257 0.004945279364 2.349850855e-06 0.06096774779 0.000440582517 2.030747703e-08 0.0001828061992 0.0003290502891 -4.360219646e-07 +6.504778186e-07 2.945465959e-08 0.04314824277 8.695492983e-07 19.3255897 0.0003433917438 1.623604604e-07 0.0004985403516 4.206613628e-07 6.048976857e-09 0.01167821222 5.984173029e-05 1.495919742e-08 0.349858651 1.798434603e-06 1.039820674e-07 1.118772128e-08 -9.889699357e-09 +1.30086451e-08 0.003069395217 4.6681898e-06 6.697088362e-11 3.350850437e-06 0.06185962328 4.883448988e-09 1.110833085e-06 -0.0003200423438 -0.0006747162613 +0.0004052046295 +0.0006808854079 0.0002346905421 4.330939823e-06 1.219218909e-07 0.00107853716 8.281691391e-06 1.767444691e-10 6.900019008e-08 605.6626124 1.58440332e-07 0.04329139184 2.752138714e-08 0.0005185712214 1.286238525e-10 2.639855477e-07 1.869433683e-06 1.079551089e-10 1.229225358e-07 0.03468655392 0.0006584274735 3.689112319e-09 1.260350348e-05 1.040877611e-05 9.878538403e-06 0.00027695052 9.100037836e-09 8.033944341e-07 3.315013452e-05 1.633327055e-11 0.003414227642 6.245172231e-08 0.0003805246583 0.0002060494734 3.729836801e-11 9.13272934e-09 2.051216421e-08 2.436679394e-08 4.8855578e-09 3.735846979e-11 1.088675628e-07 2.683031236e-07 3.706588258e-11 0.02523141233 3.549809347e-05 3.149827076e-09 1.040176277e-09 2.910870047e-10 8.86122235e-06 5.308664537e-09 3.218598972e-11 0.001577139281 -3.19884756e-07 +7.379705718e-07 2.55834443e-09 7.691641781e-10 0.01021181751 1.105572103e-06 3.947538433e-06 3.734235002e-09 6.899637353e-09 1.304337436e-09 6.538762872e-12 3.449997548e-06 -4.050730086e-09 +4.212613391e-09 7.962318534e-05 6.295652823e-05 1.652753715e-06 7.186541052e-09 4.161402504e-10 240986.4156 0.00837085449 8.760965452e-10 2.273546864e-10 3.315418794e-08 1.416471539e-06 2.236464809e-08 1.633966009e-08 4.518163364e-05 0.0005095243003 0.001992318371 3.067811242e-05 0.000600451339 7.236523643e-06 6.932191963e-06 4.091606026e-06 0.00532740987 2.395091742e-09 9.944229055e-06 1.205282983e-08 1.590006487e-08 3.215700583e-11 0.3737229127 1.296635624e-09 0.0004534080133 -1.217859786e-07 +1.233610704e-07 8.898892246e-07 7.255073266e-11 1.472337758e-06 0.02150215784 0.01192727543 4.01006022e-09 8.731812041e-09 1.44194837e-08 8.229024075e-10 0.0007855372767 2.506773438e-12 2.241537936e-06 -0.2916665128 -5.828120261e-06 +0.3085183925 +1.101920302e-05 2.752322662e-09 3.158084387e-06 1.215204067e-06 4.56919432e-09 2.790752035e-08 3.511291959e-09 5.863328277e-06 6.94327717e-05 2.22257262e-07 8.350918192e-08 5.331043956e-07 6.969438883e-11 1.760059719e-06 9.217915003e-07 0.000186591209 0.002628593815 6.211758257e-11 0.0001834032575 7.768284292e-14 1.377546165e-08 1.091734431e-06 4.677996963e-07 7.943499964e-07 2.545817678 7.562150253e-11 3.799329956e-06 1.510534633e-07 3.983823634e-09 8.184966815e-09 0.000124713148 2.309165194e-11 0.001872087514 0.01153490791 -2.14659923e-06 +2.680798294e-06 4.693787525e-11 4.26518767e-10 4.684595348e-05 0.0005076579094 7.644830233e-09 4.193877203e-13 0.002676923954 2.654050699e-08 6.800195443e-08 3.888064789e-11 0.2446884347 9.628489292e-08 0.002537052711 1.189694289e-11 1.656137485e-09 9.040798991e-08 3.488787289e-10 6.557255554e-07 9.083676565e-07 8.150458329e-13 1.315321346e-07 5.098465051e-08 3.786448639e-05 1.600711068e-06 4.658140964e-08 3.906381479e-05 0.0002389891749 8.287668941e-08 2.051670765e-10 1.60134009e-05 2.847282071e-06 4.308493985e-10 1.238922307e-09 9.741854048e-10 5.146899686e-06 0.001209019717 -3.025450861e-07 +3.233996009e-07 2.974261102e-06 3.899610282e-06 -1.695222355e-08 +1.896885723e-08 9.332118489e-06 0.004694117038 1.266652132e-11 1.65446818e-06 0.0004100868952 0.06302457519 2.236229887e-12 1.584413907e-10 0.2179964362 0.0278380659 1.936978481e-12 0.001609711018 1.892891439 0.8712003892 0.09866862844 5.934210083e-11 1.002311696e-11 0.05494208606 0.000190347461 6.722151007e-08 4.030694324e-07 9.278979899e-07 0.0003978128465 2.803383506e-11 0.0002738201279 5.960007001e-06 0.0005814501443 4.157124101e-09 5.198894043e-08 3.465884394e-09 8.167365756e-12 4.450019241e-08 4.893968894e-07 1.288321375e-11 9.163290302e-12 7.451831479e-07 3.940383302e-11 1.996508632e-08 0.1236118645 1.182043812e-06 -8.323870133e-07 +1.081851533e-06 2.618934064e-07 0.0001307650885 0.0004571929478 6.07226713e-10 9.141057675e-06 1.18356058e-08 4.465375632e-07 1.107633947e-15 7.416788486e-08 6.870175273e-05 2.600618361e-06 2.647425108e-07 2.279494937e-10 1.611105573e-05 -9.315734048e-05 +0.0001167395509 1.474187502e-09 0.003419123567 4.527389745e-05 3.271142207e-09 7.434473088e-10 0.0003102710133 4.710559791e-05 8.139484815e-10 6.749221433e-12 3.882836934e-05 8.016924727e-05 4.578651435e-05 1.374899461e-06 9.31258077e-10 9.107579154e-07 2.005597861e-06 4.853545363e-09 1.906194275e-08 1.569062034e-08 9.37420375e-06 0.0001051157618 6.464238034e-07 4.092696924e-06 1.538436694e-08 9.760057524e-12 2.626691189e-05 0.02379715614 0.00314703493 5.484543738e-08 0.00141690332 1.612112513e-12 11.09088062 0.0007046654234 3.111374555e-07 5.038033967e-09 2.21156295e-08 diff --git a/t/test_ME_generic.cc b/t/test_ME_generic.cc index c24d455..7975cdd 100644 --- a/t/test_ME_generic.cc +++ b/t/test_ME_generic.cc @@ -1,104 +1,104 @@ // Generic tester for the ME for a given set of PSP // reference weights and PSP (as LHE file) have to be given as _individual_ files #include #include "LHEF/LHEF.h" #include "RHEJ/MatrixElement.hh" #include "RHEJ/Event.hh" #include "RHEJ/YAMLreader.hh" #include "RHEJ/stream.hh" constexpr double alpha_s = 0.118; -constexpr double mu = 1234.; // coupling fixed mu doesn't matter +constexpr double mu = 1234.; // coupling fixed, mu doesn't matter constexpr double ep = 1e-6; void dump(RHEJ::UnclusteredEvent const & ev, RHEJ::Config config){ { LHEF::Writer writer{std::cout}; std::cout << std::setprecision(6); RHEJ::Event out_ev{ev, config.resummation_jets.def, config.resummation_jets.min_pt}; writer.hepeup = to_HEPEUP(std::move(out_ev), nullptr); writer.writeEvent(); } std::cout << "Rapidity ordering:\n"; for(const auto & part: ev.outgoing){ std::cout << std::setw(2) << part.type << ": "<< std::setw(7) << part.rapidity() << std::endl; } } int main(int argn, char** argv){ if(argn != 4 && argn != 5){ std::cerr << "\n# Usage:\n."<< argv[0] <<" config.yml ME_weights input_file.lhe\n\n"; return EXIT_FAILURE; } bool OUTPUT_MODE = false; if(argn == 5 && std::string("OUTPUT")==std::string(argv[4])) OUTPUT_MODE = true; const RHEJ::Config config = RHEJ::load_config(argv[1]); std::fstream wgt_file; if ( OUTPUT_MODE ) { std::cout << "_______________________USING OUTPUT MODE!_______________________" << std::endl; wgt_file.open(argv[2], std::fstream::out); wgt_file.precision(10); } else { wgt_file.open(argv[2], std::fstream::in); } RHEJ::istream in{argv[3]}; LHEF::Reader reader{in}; RHEJ::MatrixElement ME{ [](double){ return alpha_s; }, RHEJ::to_MatrixElementConfig(config) }; double max_ratio = 0.; size_t idx_max_ratio = 0; RHEJ::UnclusteredEvent ev_max_ratio; double av_ratio = 0; size_t i = 0; while(reader.readEvent()){ ++i; RHEJ::UnclusteredEvent event{reader.hepeup}; const double our_ME = ME.tree( mu, event.incoming, event.outgoing, true ); if ( OUTPUT_MODE ) { wgt_file << our_ME << std::endl; } else { std::string line; if(!std::getline(wgt_file,line)) break; const double ref_ME = std::stod(line); const double diff = std::abs(our_ME/ref_ME-1.); av_ratio+=diff; if( diff > max_ratio ) { max_ratio = diff; idx_max_ratio = i; ev_max_ratio = event; } if( diff > ep ){ size_t precision(std::cout.precision()); std::cout.precision(16); std::cout<< "Large difference in PSP " << i << "\nis: "< difference: " << diff << std::endl; std::cout.precision(precision); dump(event, config); return EXIT_FAILURE; } } } wgt_file.close(); if ( !OUTPUT_MODE ) { size_t precision(std::cout.precision()); std::cout.precision(16); std::cout << "Avg ratio after " << i << " PSP: " << av_ratio/i << std::endl; std::cout << "maximal ratio at " << idx_max_ratio << ": " << max_ratio << std::endl; std::cout.precision(precision); } return EXIT_SUCCESS; }