import numpy as np from numpy import nan, inf add_legend_handle = [ 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz', 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' ] xpoints = { 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz' : [2.361292e-01, 3.083874e-01, 3.806457e-01, 4.529039e-01, 5.251622e-01, 5.974204e-01, 6.696788e-01, 7.419370e-01, 8.141953e-01, 8.864535e-01, 9.587118e-01, 1.030970e+00, 1.103228e+00, 1.175486e+00, 1.247745e+00, 1.320003e+00, 1.392262e+00, 1.464520e+00, 1.536778e+00, 1.609036e+00, 1.681294e+00, 1.753553e+00, 1.825811e+00, 1.898069e+00, 1.970328e+00, 2.042586e+00, 2.114844e+00, 2.187102e+00, 2.259361e+00, 2.331619e+00, 2.403877e+00, 2.476135e+00, 2.548394e+00, 2.620652e+00, 2.692910e+00, 2.765169e+00, 2.837427e+00, 2.909685e+00, 2.981943e+00, 3.054202e+00, 3.126460e+00, 3.198718e+00, 3.270976e+00, 3.343234e+00, 3.415493e+00, 3.487751e+00, 3.560010e+00, 3.632268e+00, 3.704526e+00, 3.767752e+00], 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' : [2.361292e-01, 3.083874e-01, 3.806457e-01, 4.529039e-01, 5.251622e-01, 5.974204e-01, 6.696788e-01, 7.419370e-01, 8.141953e-01, 8.864535e-01, 9.587118e-01, 1.030970e+00, 1.103228e+00, 1.175486e+00, 1.247745e+00, 1.320003e+00, 1.392262e+00, 1.464520e+00, 1.536778e+00, 1.609036e+00, 1.681294e+00, 1.753553e+00, 1.825811e+00, 1.898069e+00, 1.970328e+00, 2.042586e+00, 2.114844e+00, 2.187102e+00, 2.259361e+00, 2.331619e+00, 2.403877e+00, 2.476135e+00, 2.548394e+00, 2.620652e+00, 2.692910e+00, 2.765169e+00, 2.837427e+00, 2.909685e+00, 2.981943e+00, 3.054202e+00, 3.126460e+00, 3.198718e+00, 3.270976e+00, 3.343234e+00, 3.415493e+00, 3.487751e+00, 3.560010e+00, 3.632268e+00, 3.704526e+00, 3.767752e+00], } xedges = { 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz' : [2.000000e-01, 2.722583e-01, 3.445165e-01, 4.167748e-01, 4.890331e-01, 5.612913e-01, 6.335496e-01, 7.058079e-01, 7.780661e-01, 8.503244e-01, 9.225827e-01, 9.948409e-01, 1.067099e+00, 1.139357e+00, 1.211616e+00, 1.283874e+00, 1.356132e+00, 1.428391e+00, 1.500649e+00, 1.572907e+00, 1.645165e+00, 1.717424e+00, 1.789682e+00, 1.861940e+00, 1.934198e+00, 2.006457e+00, 2.078715e+00, 2.150973e+00, 2.223231e+00, 2.295490e+00, 2.367748e+00, 2.440006e+00, 2.512265e+00, 2.584523e+00, 2.656781e+00, 2.729039e+00, 2.801298e+00, 2.873556e+00, 2.945814e+00, 3.018072e+00, 3.090331e+00, 3.162589e+00, 3.234847e+00, 3.307105e+00, 3.379364e+00, 3.451622e+00, 3.523880e+00, 3.596139e+00, 3.668397e+00, 3.740655e+00, 3.794849e+00], 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' : [2.000000e-01, 2.722583e-01, 3.445165e-01, 4.167748e-01, 4.890331e-01, 5.612913e-01, 6.335496e-01, 7.058079e-01, 7.780661e-01, 8.503244e-01, 9.225827e-01, 9.948409e-01, 1.067099e+00, 1.139357e+00, 1.211616e+00, 1.283874e+00, 1.356132e+00, 1.428391e+00, 1.500649e+00, 1.572907e+00, 1.645165e+00, 1.717424e+00, 1.789682e+00, 1.861940e+00, 1.934198e+00, 2.006457e+00, 2.078715e+00, 2.150973e+00, 2.223231e+00, 2.295490e+00, 2.367748e+00, 2.440006e+00, 2.512265e+00, 2.584523e+00, 2.656781e+00, 2.729039e+00, 2.801298e+00, 2.873556e+00, 2.945814e+00, 3.018072e+00, 3.090331e+00, 3.162589e+00, 3.234847e+00, 3.307105e+00, 3.379364e+00, 3.451622e+00, 3.523880e+00, 3.596139e+00, 3.668397e+00, 3.740655e+00, 3.794849e+00], } ref_xerrs = [ [abs(xpoints['Analysis_CSS_FragPythia8_UEAmisic.yoda.gz'][i] - xedges['Analysis_CSS_FragPythia8_UEAmisic.yoda.gz'][i]) for i in range(len(xpoints['Analysis_CSS_FragPythia8_UEAmisic.yoda.gz']))], [abs(xedges['Analysis_CSS_FragPythia8_UEAmisic.yoda.gz'][i+1] - xpoints['Analysis_CSS_FragPythia8_UEAmisic.yoda.gz'][i]) for i in range(len(xpoints['Analysis_CSS_FragPythia8_UEAmisic.yoda.gz']))] ] yvals = { 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz' : [0.000000e+00, 7.084530e+06, 0.000000e+00, 0.000000e+00, 2.719099e+07, 1.979641e+08, 4.177450e+08, 1.186238e+09, 2.449463e+09, 4.097705e+09, 6.965177e+09, 9.380214e+09, 1.030414e+10, 8.773104e+09, 6.999046e+09, 4.933843e+09, 3.179420e+09, 2.012352e+09, 1.214388e+09, 6.759855e+08, 3.778831e+08, 1.961168e+08, 1.016287e+08, 5.326736e+07, 2.613100e+07, 1.271004e+07, 6.352079e+06, 3.110735e+06, 1.441137e+06, 6.981758e+05, 3.218412e+05, 1.442525e+05, 6.756390e+04, 2.772353e+04, 1.168412e+04, 4.712228e+03, 1.852861e+03, 6.865412e+02, 2.698332e+02, 9.631203e+01, 2.653128e+01, 9.329103e+00, 2.404662e+00, 5.118160e-01, 8.239576e-02, 1.049287e-02, 3.612273e-03, 1.576703e-04, 1.276666e-04, 0.000000e+00], 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' : [0.000000e+00, 0.000000e+00, 4.938158e+06, 5.664372e+06, 4.350548e+07, 1.464427e+08, 3.537703e+08, 8.939869e+08, 1.836276e+09, 4.067714e+09, 6.672275e+09, 9.396881e+09, 1.042217e+10, 9.251763e+09, 7.691634e+09, 5.420845e+09, 3.733349e+09, 2.342582e+09, 1.331723e+09, 7.549961e+08, 4.176863e+08, 2.199200e+08, 1.128819e+08, 5.438267e+07, 2.825770e+07, 1.377997e+07, 6.653831e+06, 3.159396e+06, 1.534629e+06, 7.190900e+05, 3.164055e+05, 1.462373e+05, 6.441120e+04, 2.636309e+04, 1.149050e+04, 4.650048e+03, 1.762575e+03, 6.968032e+02, 2.371174e+02, 8.477580e+01, 2.599675e+01, 7.768919e+00, 1.812148e+00, 4.261418e-01, 1.115167e-01, 1.930765e-02, 9.038090e-04, 0.000000e+00, 0.000000e+00, 0.000000e+00], } xerrs = { 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz' : [ [3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612905e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 2.709675e-02], [3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612905e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 2.709675e-02], ], 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' : [ [3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612905e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 2.709675e-02], [3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612915e-02, 3.612915e-02, 3.612910e-02, 3.612905e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 3.612950e-02, 3.612900e-02, 3.612900e-02, 2.709675e-02], ], } yerrs = { 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz' : [ [0.000000e+00, 7.084530e+06, 0.000000e+00, 0.000000e+00, 1.159999e+07, 4.369441e+07, 4.508829e+07, 8.492373e+07, 1.151544e+08, 1.413181e+08, 1.966814e+08, 1.995938e+08, 1.975207e+08, 1.613992e+08, 1.230754e+08, 8.596318e+07, 5.657020e+07, 3.224283e+07, 1.979142e+07, 1.025605e+07, 5.762062e+06, 2.725463e+06, 1.382292e+06, 7.701343e+05, 3.666004e+05, 1.910873e+05, 8.727902e+04, 4.320291e+04, 2.246166e+04, 1.057962e+04, 5.350414e+03, 2.326314e+03, 1.343141e+03, 5.355260e+02, 2.245219e+02, 9.722047e+01, 4.931246e+01, 2.051969e+01, 8.629665e+00, 3.879162e+00, 1.198978e+00, 5.535158e-01, 2.090691e-01, 5.333588e-02, 1.278607e-02, 3.355911e-03, 1.512857e-03, 1.576703e-04, 1.276666e-04, 0.000000e+00], [0.000000e+00, 7.084530e+06, 0.000000e+00, 0.000000e+00, 1.159999e+07, 4.369441e+07, 4.508829e+07, 8.492373e+07, 1.151544e+08, 1.413181e+08, 1.966814e+08, 1.995938e+08, 1.975207e+08, 1.613992e+08, 1.230754e+08, 8.596318e+07, 5.657020e+07, 3.224283e+07, 1.979142e+07, 1.025605e+07, 5.762062e+06, 2.725463e+06, 1.382292e+06, 7.701343e+05, 3.666004e+05, 1.910873e+05, 8.727902e+04, 4.320291e+04, 2.246166e+04, 1.057962e+04, 5.350414e+03, 2.326314e+03, 1.343141e+03, 5.355260e+02, 2.245219e+02, 9.722047e+01, 4.931246e+01, 2.051969e+01, 8.629665e+00, 3.879162e+00, 1.198978e+00, 5.535158e-01, 2.090691e-01, 5.333588e-02, 1.278607e-02, 3.355911e-03, 1.512857e-03, 1.576703e-04, 1.276666e-04, 0.000000e+00], ], 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' : [ [0.000000e+00, 0.000000e+00, 4.938158e+06, 4.435657e+06, 1.905215e+07, 2.781074e+07, 4.412101e+07, 7.049922e+07, 1.006189e+08, 1.556924e+08, 1.788569e+08, 2.036093e+08, 1.978146e+08, 1.691288e+08, 1.411196e+08, 9.746518e+07, 7.327826e+07, 4.455295e+07, 2.263520e+07, 1.484917e+07, 6.784496e+06, 3.680978e+06, 2.136235e+06, 7.914542e+05, 4.367073e+05, 2.348154e+05, 9.818717e+04, 4.501224e+04, 2.292104e+04, 1.126376e+04, 5.145770e+03, 2.397670e+03, 1.172337e+03, 4.498262e+02, 2.830895e+02, 1.064169e+02, 4.127258e+01, 2.416971e+01, 7.279222e+00, 3.199095e+00, 1.162852e+00, 4.850638e-01, 1.475298e-01, 4.395583e-02, 1.634030e-02, 4.705511e-03, 6.422896e-04, 0.000000e+00, 0.000000e+00, 0.000000e+00], [0.000000e+00, 0.000000e+00, 4.938158e+06, 4.435657e+06, 1.905215e+07, 2.781074e+07, 4.412101e+07, 7.049922e+07, 1.006189e+08, 1.556924e+08, 1.788569e+08, 2.036093e+08, 1.978146e+08, 1.691288e+08, 1.411196e+08, 9.746518e+07, 7.327826e+07, 4.455295e+07, 2.263520e+07, 1.484917e+07, 6.784496e+06, 3.680978e+06, 2.136235e+06, 7.914542e+05, 4.367073e+05, 2.348154e+05, 9.818717e+04, 4.501224e+04, 2.292104e+04, 1.126376e+04, 5.145770e+03, 2.397670e+03, 1.172337e+03, 4.498262e+02, 2.830895e+02, 1.064169e+02, 4.127258e+01, 2.416971e+01, 7.279222e+00, 3.199095e+00, 1.162852e+00, 4.850638e-01, 1.475298e-01, 4.395583e-02, 1.634030e-02, 4.705511e-03, 6.422896e-04, 0.000000e+00, 0.000000e+00, 0.000000e+00], ], } variation_yvals = { } # lists for ratio plot ratio0_yvals = { 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz' : [1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00], 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' : [1.000000e+00, 0.000000e+00, nan, nan, 1.599996e+00, 7.397435e-01, 8.468570e-01, 7.536323e-01, 7.496647e-01, 9.926809e-01, 9.579477e-01, 1.001777e+00, 1.011455e+00, 1.054560e+00, 1.098955e+00, 1.098706e+00, 1.174223e+00, 1.164102e+00, 1.096621e+00, 1.116882e+00, 1.105332e+00, 1.121373e+00, 1.110728e+00, 1.020938e+00, 1.081386e+00, 1.084179e+00, 1.047505e+00, 1.015643e+00, 1.064874e+00, 1.029955e+00, 9.831108e-01, 1.013760e+00, 9.533376e-01, 9.509285e-01, 9.834284e-01, 9.868045e-01, 9.512721e-01, 1.014947e+00, 8.787553e-01, 8.802203e-01, 9.798528e-01, 8.327616e-01, 7.535980e-01, 8.326075e-01, 1.353427e+00, 1.840073e+00, 2.502051e-01, 0.000000e+00, 0.000000e+00, 1.000000e+00], } ratio0_yerrs = { 'Analysis_CSS_FragPythia8_UEAmisic.yoda.gz' : [ [1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 4.266115e-01, 2.207188e-01, 1.079325e-01, 7.159084e-02, 4.701208e-02, 3.448713e-02, 2.823782e-02, 2.127817e-02, 1.916907e-02, 1.839705e-02, 1.758460e-02, 1.742317e-02, 1.779262e-02, 1.602246e-02, 1.629744e-02, 1.517199e-02, 1.524826e-02, 1.389714e-02, 1.360140e-02, 1.445790e-02, 1.402933e-02, 1.503435e-02, 1.374023e-02, 1.388833e-02, 1.558607e-02, 1.515324e-02, 1.662439e-02, 1.612669e-02, 1.987957e-02, 1.931666e-02, 1.921598e-02, 2.063153e-02, 2.661423e-02, 2.988851e-02, 3.198147e-02, 4.027703e-02, 4.519109e-02, 5.933215e-02, 8.694325e-02, 1.042091e-01, 1.551787e-01, 3.198277e-01, 4.188101e-01, 1.000000e+00, 1.000000e+00, 1.000000e+00], [1.000000e+00, 1.000000e+00, 1.000000e+00, 1.000000e+00, 4.266115e-01, 2.207188e-01, 1.079325e-01, 7.159084e-02, 4.701208e-02, 3.448713e-02, 2.823782e-02, 2.127817e-02, 1.916907e-02, 1.839705e-02, 1.758460e-02, 1.742317e-02, 1.779262e-02, 1.602246e-02, 1.629744e-02, 1.517199e-02, 1.524826e-02, 1.389714e-02, 1.360140e-02, 1.445790e-02, 1.402933e-02, 1.503435e-02, 1.374023e-02, 1.388833e-02, 1.558607e-02, 1.515324e-02, 1.662439e-02, 1.612669e-02, 1.987957e-02, 1.931666e-02, 1.921598e-02, 2.063153e-02, 2.661423e-02, 2.988851e-02, 3.198147e-02, 4.027703e-02, 4.519109e-02, 5.933215e-02, 8.694325e-02, 1.042091e-01, 1.551787e-01, 3.198277e-01, 4.188101e-01, 1.000000e+00, 1.000000e+00, 1.000000e+00], ], 'Analysis_Alaric_FragPythia8_UEAmisic.yoda.gz' : [ [1.000000e+00, 0.000000e+00, nan, nan, 7.006787e-01, 1.404838e-01, 1.056171e-01, 5.943095e-02, 4.107796e-02, 3.799503e-02, 2.567873e-02, 2.170625e-02, 1.919759e-02, 1.927810e-02, 2.016269e-02, 1.975441e-02, 2.304768e-02, 2.213974e-02, 1.863917e-02, 2.196669e-02, 1.795395e-02, 1.876931e-02, 2.102000e-02, 1.485814e-02, 1.671223e-02, 1.847479e-02, 1.545749e-02, 1.446997e-02, 1.590484e-02, 1.613313e-02, 1.598854e-02, 1.662134e-02, 1.735154e-02, 1.622543e-02, 2.422856e-02, 2.258313e-02, 2.227505e-02, 3.520504e-02, 2.697674e-02, 3.321594e-02, 4.382945e-02, 5.199469e-02, 6.135160e-02, 8.588210e-02, 1.983148e-01, 4.484484e-01, 1.778076e-01, 0.000000e+00, 0.000000e+00, 1.000000e+00], [1.000000e+00, 0.000000e+00, nan, nan, 7.006787e-01, 1.404838e-01, 1.056171e-01, 5.943095e-02, 4.107796e-02, 3.799503e-02, 2.567873e-02, 2.170625e-02, 1.919759e-02, 1.927810e-02, 2.016269e-02, 1.975441e-02, 2.304768e-02, 2.213974e-02, 1.863917e-02, 2.196669e-02, 1.795395e-02, 1.876931e-02, 2.102000e-02, 1.485814e-02, 1.671223e-02, 1.847479e-02, 1.545749e-02, 1.446997e-02, 1.590484e-02, 1.613313e-02, 1.598854e-02, 1.662134e-02, 1.735154e-02, 1.622543e-02, 2.422856e-02, 2.258313e-02, 2.227505e-02, 3.520504e-02, 2.697674e-02, 3.321594e-02, 4.382945e-02, 5.199469e-02, 6.135160e-02, 8.588210e-02, 1.983148e-01, 4.484484e-01, 1.778076e-01, 0.000000e+00, 0.000000e+00, 1.000000e+00], ], } ratio0_variation_vals = { } ratio_band_edges = { }