APFEL - Vrap agreement.

[scarlehoff]

Just for the NLO grids.

:heavy_check_mark: E605

Good agreement (per mile) up to a global factor (i.e., for all masses and rapidities) of ~54.37. I would call this a success. I will investigate what the global factor is, but for sure not due to the way I'm filling the grid.

Numbers below

d validphys pineappl*54.35 ratio ``` 0 393.5439112199864 393.3844725423387 1.0004052998752526 1 336.4336411779308 336.70153333003816 0.9992043631359268 2 288.3770877707228 288.4759960780084 0.9996571350523777 3 248.47075046525396 248.1133030888043 1.0014406618750376 4 213.13876670878267 212.94265574640602 1.0009209566852129 5 183.30722013944538 183.43137517219685 0.9993231526905638 6 157.75121470078332 157.65752926757995 1.000594233803096 7 136.06855846579222 135.97662785526327 1.000676076558001 8 117.05920552147818 116.97340286072922 1.0007335228235696 9 100.6331879246386 100.67748978897846 0.9995599625652892 10 31.241466846486865 31.27080935092567 0.9990616646947151 11 27.048219881547837 27.057840189763425 0.9996444539494609 12 23.4317824073186 23.429329780164153 1.0001046819169588 17 410.8828300457184 410.81949628111386 1.000154164457086 18 352.52233283798785 352.3924783716206 1.000368493865043 19 302.3811319358382 302.5667512419804 0.9993865178332376 20 260.7771691633508 260.757323864068 1.0000761063926746 21 224.14984732196297 224.06702991992194 1.000369609942483 22 193.33300621375426 193.20240853992468 1.0006759629696986 23 166.2242367038938 166.1822728097973 1.0002525172714694 24 143.5738798382832 143.41156817593867 1.001131789188341 25 123.6486426492435 123.44594803951385 1.0016419705381077 26 106.24863121040767 106.29088875272711 0.999602434951713 27 32.62520724475262 32.67171944507492 0.9985763773345174 28 28.24466884907289 28.170676363391152 1.0026265782449546 29 24.296102143112716 24.295250597270076 1.000035049889246 35 425.3125689163745 425.19022647668555 1.0002877357758262 36 365.5820521246038 365.17444280292887 1.001116204405068 37 314.01384898564345 313.9061552935672 1.0003430760762735 38 270.62623055121037 270.8244982854522 0.9992679106377118 39 233.5108147185213 233.13434272108202 1.001614828570708 40 201.63788653400906 201.38212839196174 1.0012700140975246 41 173.26694377641405 173.54052717532085 0.9984235186825818 42 150.22365347979905 150.04748557375328 1.0011740810275636 43 129.4496490450385 129.3159545849007 1.0010338589740682 44 111.25261883992381 111.4598317884077 0.9981409181661314 45 34.37354572827482 34.419313442145395 0.9986702897503315 46 29.760844465379044 29.683883002299755 1.0025927020084713 47 25.55764250535603 25.60205904395414 0.9982651184999669 53 434.9771201596417 435.0411816534994 0.9998527461386202 54 374.314340795701 374.21041966781377 1.0002777077345402 55 322.09924323821144 322.2148015341709 0.9996413625463224 56 278.43833957120876 278.46189433708435 0.9999154111698778 57 240.00606838543223 239.9839806675437 1.0000920383011696 58 207.70444017195518 207.5346975496277 1.0008178999672424 59 179.11495935747627 179.05341575235195 1.000343716453918 60 155.1168704183285 154.9865325757307 1.0008409623754508 61 134.05206759922038 133.82554439102978 1.0016926754097761 62 115.5726114942246 115.59076970299151 0.9998429095263094 63 36.330573229961665 36.37010772536869 0.9989129948224088 64 31.517630222010787 31.437231807481414 1.0025574266532666 65 27.162948885401903 27.174879862927938 0.9995609556477815 71 440.62011431381234 440.44605301246764 1.000395193236843 72 378.7284451964183 378.9812384812921 0.9993329662283895 73 326.3087142425274 326.44854448078866 0.9995716622401131 74 282.6253970020186 282.27477921035523 1.0012421151923108 75 243.7617557552672 243.59721567001924 1.0006754596303384 76 210.8338177748659 210.98131042595432 0.999300920774496 77 182.44913238296772 182.35069609652788 1.0005398185394792 78 158.1705695864468 158.13505759101855 1.0002245675055816 79 136.85363773274585 136.7776208598398 1.0005557698140104 80 118.36248901524112 118.36018803640785 1.000019440479704 81 38.09677624880653 38.14001519644653 0.9988663101622459 82 33.112730708328485 33.09805659166506 1.0004433528180963 83 28.730775062180292 28.729233744184494 1.000053649812227 89 439.75738322970597 440.52081341649 0.9982669827088005 90 378.84007990758295 379.07711663954655 0.9993747004987669 91 326.55641591499756 326.57464840937524 0.9999441705151747 92 282.2344642438793 282.4321243935085 0.9993001499030832 93 244.04576918075077 243.82588148381754 1.0009018226268478 94 211.31991951251388 211.2884724305839 1.0001488348207936 95 182.53657719083876 182.73698030516033 0.9989033248005581 96 158.82052890281926 158.60878192915655 1.0013350267941488 97 137.43868085249338 137.38474998563908 1.0003925535174751 98 119.14350825258 119.08048791654268 1.000529224704567 99 39.11719987416068 39.06362953569447 1.001371361010304 100 34.07798149132137 34.00003420666872 1.0022925648891658 101 29.619835961654623 29.606044487534096 1.0004658330540013 107 179.86997317177102 179.85072797688966 1.0001070064886468 108 155.78559481554112 156.10155300720822 0.9979759445977292 109 135.35532020804663 135.2241923100314 1.0009697073857509 110 117.26080406897252 117.22786933640178 1.0002809462694935 111 38.64666109582427 38.59324200842349 1.0013841565160326 112 33.524710080321825 33.58618444036067 0.9981696533541045 113 29.301647840704575 29.234273700528245 1.0023046284941608 ```

Numbers without factor applied

``` data vp pine ratio 0 393.5439112199864 7.237984775388016 54.37202804821998 1 336.4336411779308 6.195060410856268 54.30675713643762 2 288.3770877707228 5.307746017994634 54.331365290096734 3 248.47075046525396 4.565102172747089 54.42829997290829 4 213.13876670878267 3.9179881462080224 54.400053995841326 5 183.30722013944538 3.375002303076299 54.31321334873215 6 157.75121470078332 2.900782507223182 54.38229660719826 7 136.06855846579222 2.5018698777417345 54.38674476092735 8 117.05920552147818 2.1522245236564714 54.38986696546101 9 100.6331879246386 1.8523917164485457 54.32608396542347 10 31.241466846486865 0.5753598776619259 54.29900147615777 11 27.048219881547837 0.49784434571781827 54.3306760721532 12 23.4317824073186 0.4310824246580341 54.35568946218671 17 410.8828300457184 7.558776380517274 54.35837883824263 18 352.52233283798785 6.483762251547757 54.37002764156509 19 302.3811319358382 5.567005542630734 54.31665724423646 20 260.7771691633508 4.797742849384876 54.354136382441865 21 224.14984732196297 4.122668443788812 54.37008830037395 22 193.33300621375426 3.55478212584958 54.38673858740312 23 166.2242367038938 3.0576315144396924 54.36372431370437 24 143.5738798382832 2.638667307744962 54.41151274238632 25 123.6486426492435 2.271314591343401 54.43924109874615 26 106.24863121040767 1.9556741260851354 54.3283923396256 27 32.62520724475262 0.6011355923656839 54.27262610813102 28 28.24466884907289 0.5183197122979053 54.49275452761328 29 24.296102143112716 0.4470147304005534 54.35190496148052 35 425.3125689163745 7.823187239681427 54.36563843941615 36 365.5820521246038 6.718940989934294 54.41066570941546 37 314.01384898564345 5.775642231712369 54.368646184745465 38 270.62623055121037 4.982971449594337 54.31021094315963 39 233.5108147185213 4.289500326054867 54.43776593281799 40 201.63788653400906 3.7052829510940524 54.419025266200464 41 173.26694377641405 3.193017979306731 54.264318240398325 42 150.22365347979905 2.7607633040248993 54.413811303848085 43 129.4496490450385 2.379318391626508 54.40619023524061 44 111.25261883992381 2.0507788737517516 54.248958902329235 45 34.37354572827482 0.6332900357340459 54.27773024793052 46 29.760844465379044 0.5461616007782844 54.49091335416041 47 25.55764250535603 0.4710590440469943 54.255709190473205 53 434.9771201596417 8.004437564921792 54.341996752634 54 374.314340795701 6.885196314035212 54.36509341537226 55 322.09924323821144 5.928515207620439 54.33050805439263 56 278.43833957120876 5.12349391604571 54.34540259708286 57 240.00606838543223 4.415528623137878 54.35500228166856 58 207.70444017195518 3.8184856954853306 54.39445286321961 59 179.11495935747627 3.294451071800404 54.368680989270445 60 155.1168704183285 2.8516381338680903 54.395706305105755 61 134.05206759922038 2.4622915251339426 54.44199690852133 62 115.5726114942246 2.1267850911313984 54.34146213275492 63 36.330573229961665 0.6691832148181911 54.29092126859793 64 31.517630222010787 0.5784219283805228 54.488996138605046 65 27.162948885401903 0.4999977895662914 54.32613793945693 71 440.62011431381234 8.1038832201006 54.37147875242241 72 378.7284451964183 6.972975869021013 54.31374671451297 73 326.3087142425274 6.006412961928034 54.326719842750144 74 282.6253970020186 5.1936481915428745 54.417508960702094 75 243.7617557552672 4.482009487948836 54.3867112309089 76 210.8338177748659 3.881900835804127 54.31200504409386 77 182.44913238296772 3.3551186034319755 54.3793391376207 78 158.1705695864468 2.9095686769276643 54.362205243928365 79 136.85363773274585 2.5166075595186714 54.380206089391464 80 118.36248901524112 2.177740350255894 54.35105659007191 81 38.09677624880653 0.7017482096862286 54.28838395731806 82 33.112730708328485 0.6089798820913535 54.374096225663536 83 28.730775062180292 0.5285967570227138 54.35291586729453 89 439.75738322970597 8.105258756513155 54.25581051022331 90 378.84007990758295 6.974739956569393 54.316014972107986 91 326.55641591499756 6.0087331814052485 54.346965667499745 92 282.2344642438793 5.196543227111472 54.311963147232575 93 244.04576918075077 4.486216770631417 54.399014059769186 94 211.31991951251388 3.8875523906271185 54.35808917251014 95 182.53657719083876 3.3622259485770067 54.29039570291034 96 158.82052890281926 2.918284856102236 54.42255870626199 97 137.43868085249338 2.527778288604215 54.371335283674775 98 119.14350825258 2.190993337930868 54.37876336269322 99 39.11719987416068 0.7187420337754272 54.42453347091003 100 34.07798149132137 0.6255756063784493 54.47460090172616 101 29.619835961654623 0.5447294293934516 54.375318026484976 107 179.86997317177102 3.3091210299335723 54.35581580265795 108 155.78559481554112 2.872153689185064 54.23999258888658 109 135.35532020804663 2.4880256174798787 54.402703596415556 110 117.26080406897252 2.1569065195290116 54.365269429746974 111 38.64666109582427 0.7100872494650136 54.425228906646375 112 33.524710080321825 0.617961075259626 54.25052065979557 113 29.301647840704575 0.5378891205248987 54.47525655865764 ```

:heavy_check_mark: E886

For this one I can reproduce the results of the NNPDF datasets DYE886P and DYE886R (which is a ratio of the deut over proton). For the first one (DYE88P) it is necessary to multiply the vrap cross section by the same factor of 54 as before and (mass/sqrtS)**3. This factor is found in the runEXXX.cc scripts of externals but it is not necessary for E605. Not sure about what the issue is. In any case, for DYE886P and DYE886R the agreement is more than excelent, numbers bellow.

DYE886P

``` data vp pine*54.35 ratio 0 5.465659737648973 5.464610076069676 1.0001920835274034 1 4.783060877086962 4.783517182758256 0.9999046087525433 2 4.161363924387236 4.162449851821882 0.9997391133891568 3 3.603105659202399 3.5993189585132264 1.0010520603294177 4 3.0973978929542634 3.0931228658390477 1.001382107113309 5 2.645130778807017 2.6391149278011814 1.0022794956530552 6 2.226147426417468 2.2302580548312627 0.9981568821576992 7 1.8724335807709964 1.8694505760993239 1.001595658483733 8 1.5598919908074356 1.5545330104647772 1.0034473248921592 10 5.4807183747103965 5.476743768137176 1.0007257243978336 11 4.761537546570885 4.7675649177997785 0.9987357547652072 12 4.121772914340678 4.128961894445791 0.9982588892101951 13 3.563349307700689 3.5578474161711995 1.001546410198617 14 3.0563384329392687 3.0500572640304964 1.0020593609775288 15 2.6002600073919195 2.59817936182737 1.000800809056957 16 2.2022932254441625 2.196708500319393 1.0025423150700048 17 1.8470104694503924 1.843965542639693 1.0016512926843202 18 1.5354822924410063 1.5359493521691818 0.999695914629271 20 5.447983684666091 5.446791938029558 1.0002187978997714 21 4.691874796295272 4.709464347899645 0.9962650632205724 22 4.051455215605245 4.057001534092247 0.9986329020483734 23 3.478413068549578 3.4800297107604585 0.9995354516066683 24 2.9740099813862546 2.974560234322525 0.9998150136850749 25 2.5356667532748176 2.5299461517537427 1.0022611554467709 26 2.142232551219016 2.140204532546692 1.000947581710759 27 1.8042748015143624 1.7993287532556201 1.0027488296676152 28 1.5039300349664328 1.5035741138862109 1.0002367166852202 30 5.322909801132186 5.326657041733677 0.9992965117573495 31 4.580948366140474 4.581103865131712 0.9999660564362179 32 3.9282382336343367 3.927965494799325 1.0000694351402457 33 3.3592332623819354 3.35697856438951 1.0006716450371005 34 2.8644918288282435 2.8619804509160254 1.0008774965291654 35 2.432702701818389 2.431681746097864 1.0004198558146697 36 2.0548357854823047 2.058143741949571 0.9983927476007419 37 1.7372993484005628 1.7337265019297263 1.002060790134349 38 1.4557416827878473 1.4537402670533837 1.0013767354319218 42 5.095984950919991 5.091725079969108 1.000836626267911 43 4.363184050264116 4.366174317563032 0.9993151287416795 44 3.7362843550111835 3.7318307004182816 1.0011934235367115 45 3.1854249728904893 3.1831090725354816 1.00072755922032 46 2.7079501149330105 2.708095527293655 0.9999463045674796 47 2.3058999729773717 2.3010835821505435 1.002093096862795 48 1.9505716579493524 1.948086408931246 1.001275738594917 49 1.6468135561755077 1.6444489335125705 1.0014379422885977 50 1.3840385058261189 1.3828495825134448 1.0008597632943659 54 4.743375789074989 4.7438115798517355 0.9999081348891264 55 4.0581536194238925 4.0607616838479785 0.9993577400923429 56 3.469998286723662 3.466442331031852 1.0010258228328155 57 2.9511805448912454 2.9511803728720944 1.0000000582882538 58 2.514414837387943 2.5108844968183104 1.0014060147227426 59 2.1375020376846967 2.1322678803775172 1.002454737209779 60 1.8079916652853671 1.8056246845872135 1.0013108929659402 61 1.5243994551489402 1.525125452362275 0.9995239754131632 62 1.2897369168719202 1.2852822698516397 1.003465890042033 66 4.2662591596166575 4.289243259113009 0.9946414558214857 67 3.666308071157606 3.6663871267221486 0.9999784377476217 68 3.1242293836247734 3.1274087819047565 0.9989833761744294 69 2.6595774320150767 2.661352293800126 0.9993330977679339 70 2.2641690572028743 2.2624068938049775 1.0007788888032132 71 1.9171138186879628 1.9215126359877979 0.9977107528634211 72 1.62974314542598 1.6269943247504046 1.001689508459716 73 1.3780147461714223 1.373941541928328 1.002964612480803 74 1.1545423612552905 1.1566476863823476 0.9981798043156582 78 3.7495348118038443 3.7484480694672455 1.0002899179384264 79 3.187317224678795 3.204636134761223 0.9945956703494143 80 2.725118681936777 2.7274406444226553 0.9991486661714798 81 2.3235934817640036 2.31588277563865 1.0033294889562048 82 1.9734298541067328 1.9695033324073008 1.001993660855924 83 1.6759557901830093 1.6703999324182106 1.0033260644094708 84 1.416728052896304 1.4135607811410367 1.0022406335811826 85 1.1888394984860804 1.190687091255299 0.9984482969683742 86 1.006125171908575 1.0021913323133538 1.0039252380941481 90 3.1377141692936923 3.1455643975405962 0.9975043498543403 91 2.6844011354975317 2.684896083009364 0.9998156548720956 92 2.275787859638898 2.2781925206100273 0.9989444873734878 93 1.9389391038296546 1.9373245599164417 1.000833388450556 94 1.6445878759130688 1.6430461835866612 1.000938313445969 95 1.3959299495506834 1.395157703603685 1.0005535187491734 96 1.1765734687898948 1.1767819097800933 0.9998228720305214 97 0.9949453002121362 0.9913858933810864 1.0035903343539725 98 0.8344829874314138 0.8340508790387636 1.0005180839724648 103 2.146338456836429 2.1545099135470847 0.9962072782031426 104 1.832557257066124 1.8278512135979887 1.0025746315855062 105 1.5442005192160548 1.5464156899493398 0.9985675450994955 106 1.3158603513137357 1.31329900179275 1.0019503171155155 107 1.1133303322106982 1.112509936562661 1.0007374277038568 108 0.9395518277945049 0.9358498179682271 1.0039557734106472 109 0.7887035484669978 0.7877046803179686 1.0012680744116258 110 0.6656840475575747 0.6653069202722249 1.0005668470804354 ```

DYE886R

``` data vp pine ratio 0 1.0575337673933138 1.057509727136195 1.000022732894556 1 1.081730248434689 1.081983093264858 0.9997663135110493 2 1.1027612358774566 1.1027859334832713 0.9999776043517923 3 1.1225336839613018 1.122840718860291 0.9997265552505963 4 1.14334021950913 1.1436327862475362 0.9997441777273927 5 1.1516665382797153 1.1515064730292144 1.0001390050809527 6 1.1555183000864855 1.1555936858036338 0.9999347645127571 7 1.1557930841553183 1.1557961464618363 0.9999973504786919 8 1.1563856915478905 1.1566913355419486 0.9997357601075875 9 1.154800618188727 1.1550146337343867 0.9998147075029102 10 1.151044899626755 1.1512556200206583 0.999816964720746 11 1.144923041493198 1.1448929404112247 1.0000262916128755 12 1.1345644688922316 1.1341946172599031 1.000326091859986 13 1.0872104710993686 1.0872169414124813 0.9999940487377761 14 1.0565369932338926 1.0572378712692545 0.9993370668471037 ```

E905

Not sure how to begin as the grid in the NNPDF repos have a very different number of points and not sure how they match with the runcards in vrap.

[cschwan]

@scarlehoff can you paste the actual numbers here so we have them for reference?

[scarlehoff]

You were too fast! I'm trying to put them in beautiful foldable markdown because there are ~100 of them (plus cuts being applied)

[scarlehoff]

Updated comment.

[cschwan]

@scarlehoff is E905 supposed to be E906 or are those two different experiments?

[scarlehoff]

I'm just bad at numbers.

[cschwan]

OK, does that mean with https://github.com/scarlehoff/Hawaiian_vrap/pull/7 closed this Issue can be closed as well? Is there anything else to be done?

[scarlehoff]

Yes. Positivity, but that I'll do directly in the runcard runner. Now the question is how to make this public.