import numpy as np from numpy import nan add_legend_handle = [ 'alpha0.yoda', '_alpha0.yoda' ] xpoints = { 'alpha0.yoda' : [-8.774222, -6.3226665, -3.871111, -1.4195554000000001, 1.0320001, 3.4835555, 5.935111, 8.3866665, 10.838222, 13.28978, 15.741335000000001, 18.19289, 20.644444999999997, 23.096, 25.547555000000003, 27.99911, 30.450665, 32.90222, 35.35378, 37.805335, 40.25689, 42.708445, 45.16, 47.611554999999996, 50.06311, 52.514665, 54.96622, 57.41778, 59.869335, 62.32089, 64.772445, 67.224, 69.675555, 72.12711, 74.578665, 77.03022, 79.48178, 81.933335, 84.38489, 86.836445, 89.288, 91.739555, 94.19111, 96.642665, 99.09422], '_alpha0.yoda' : [-8.774222, -6.3226665, -3.871111, -1.4195554000000001, 1.0320001, 3.4835555, 5.935111, 8.3866665, 10.838222, 13.28978, 15.741335000000001, 18.19289, 20.644444999999997, 23.096, 25.547555000000003, 27.99911, 30.450665, 32.90222, 35.35378, 37.805335, 40.25689, 42.708445, 45.16, 47.611554999999996, 50.06311, 52.514665, 54.96622, 57.41778, 59.869335, 62.32089, 64.772445, 67.224, 69.675555, 72.12711, 74.578665, 77.03022, 79.48178, 81.933335, 84.38489, 86.836445, 89.288, 91.739555, 94.19111, 96.642665, 99.09422], } xedges = { 'alpha0.yoda' : [-10.0, -7.548444, -5.096889, -2.645333, -0.1937777999999999, 2.257778, 4.709333, 7.160889, 9.612444, 12.064, 14.51556, 16.96711, 19.41867, 21.87022, 24.32178, 26.77333, 29.22489, 31.67644, 34.128, 36.57956, 39.03111, 41.48267, 43.93422, 46.38578, 48.83733, 51.28889, 53.74044, 56.192, 58.64356, 61.09511, 63.54667, 65.99822, 68.44978, 70.90133, 73.35289, 75.80444, 78.256, 80.70756, 83.15911, 85.61067, 88.06222, 90.51378, 92.96533, 95.41689, 97.86844, 100.32], '_alpha0.yoda' : [-10.0, -7.548444, -5.096889, -2.645333, -0.1937777999999999, 2.257778, 4.709333, 7.160889, 9.612444, 12.064, 14.51556, 16.96711, 19.41867, 21.87022, 24.32178, 26.77333, 29.22489, 31.67644, 34.128, 36.57956, 39.03111, 41.48267, 43.93422, 46.38578, 48.83733, 51.28889, 53.74044, 56.192, 58.64356, 61.09511, 63.54667, 65.99822, 68.44978, 70.90133, 73.35289, 75.80444, 78.256, 80.70756, 83.15911, 85.61067, 88.06222, 90.51378, 92.96533, 95.41689, 97.86844, 100.32], } ref_xerrs = [ [abs(xpoints['alpha0.yoda'][i] - xedges['alpha0.yoda'][i]) for i in range(len(xpoints['alpha0.yoda']))], [abs(xedges['alpha0.yoda'][i+1] - xpoints['alpha0.yoda'][i]) for i in range(len(xpoints['alpha0.yoda']))] ] yvals = { 'alpha0.yoda' : [0.0, 0.0, 0.0, 0.0, 2.5326545698042033, 0.8602181064671198, 0.5414708046644661, 0.44509586772477067, 0.4147818773056785, 0.4143822708805821, 0.43384022353205104, 0.4617802542054861, 0.498902327099182, 0.5437794710306905, 0.5895633374803693, 0.643996067809885, 0.701284901388917, 0.7682035112336633, 0.8397648844001367, 0.9142350757683925, 0.9992996296235863, 1.0881593277722237, 1.1846640506453066, 1.2925361505985984, 1.4172335166179897, 1.5497607635985413, 1.7015092430941925, 1.8689960678098838, 2.0597515857314783, 2.2905831388993128, 2.556853011360155, 2.868539215846236, 3.255160612673618, 3.742010393382172, 4.368546429809716, 5.215630047806293, 6.41062017654065, 8.292859619424453, 11.730657214181988, 20.241014868144664, 75.4442477442934, 416.62540025698075, 0.0, 0.0, 0.0], '_alpha0.yoda' : [0.0, 0.0, 0.0, 0.0, 2.528346, 0.8586234, 0.5404824, 0.4436655, 0.4136297, 0.4150145, 0.4330568, 0.4616596, 0.4994486, 0.5416651, 0.5904333, 0.6437876, 0.7037288, 0.7684984, 0.8381973, 0.9144926, 1.001988, 1.08963, 1.190752, 1.296746, 1.412251, 1.547844, 1.696953, 1.867383, 2.063406, 2.28981, 2.552768, 2.868955, 3.257372, 3.741224, 4.365709, 5.208933, 6.410612, 8.300188, 11.7348, 20.23518, 75.45669, 416.6254, 0.0, 0.0, 0.0], } xerrs = { 'alpha0.yoda' : [ [1.225778, 1.2257775000000004, 1.225778, 1.2257775999999998, 1.2257779000000002, 1.2257775, 1.225778, 1.2257774999999995, 1.225778, 1.2257800000000003, 1.2257750000000005, 1.2257800000000003, 1.2257750000000023, 1.2257800000000003, 1.2257749999999987, 1.2257799999999968, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257750000000058, 1.2257800000000003, 1.2257750000000058, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.225775000000013, 1.225779999999986], [1.225778, 1.2257774999999995, 1.225778, 1.2257776, 1.2257778999999998, 1.2257775, 1.225778, 1.2257775000000004, 1.225778, 1.2257800000000003, 1.2257750000000005, 1.2257799999999968, 1.2257749999999987, 1.2257800000000003, 1.2257750000000023, 1.2257800000000003, 1.2257750000000023, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257799999999932, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257750000000058, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003], ], '_alpha0.yoda' : [ [1.225778, 1.2257775000000004, 1.225778, 1.2257775999999998, 1.2257779000000002, 1.2257775, 1.225778, 1.2257774999999995, 1.225778, 1.2257800000000003, 1.2257750000000005, 1.2257800000000003, 1.2257750000000023, 1.2257800000000003, 1.2257749999999987, 1.2257799999999968, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257750000000058, 1.2257800000000003, 1.2257750000000058, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.225775000000013, 1.225779999999986], [1.225778, 1.2257774999999995, 1.225778, 1.2257776, 1.2257778999999998, 1.2257775, 1.225778, 1.2257775000000004, 1.225778, 1.2257800000000003, 1.2257750000000005, 1.2257799999999968, 1.2257749999999987, 1.2257800000000003, 1.2257750000000023, 1.2257800000000003, 1.2257750000000023, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257799999999932, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257750000000058, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003, 1.2257749999999987, 1.2257800000000003], ], } yerrs = { 'alpha0.yoda' : [ [0.0, 0.0, 0.0, 0.0, 0.21407399665430032, 0.07272354967452557, 0.04578829623047529, 0.0376382803462929, 0.03508005328587472, 0.03504977707849241, 0.03669240506776881, 0.03906038156967019, 0.04219171561882855, 0.04597939642755021, 0.04985220741571528, 0.054454144697553906, 0.05929082780163272, 0.06494959299453266, 0.07100014274932817, 0.07729255754314958, 0.08448143165218135, 0.09198801357932859, 0.1001412908027922, 0.109260795045065, 0.11980042219786104, 0.13100121834835385, 0.14382409236715624, 0.157987628258901, 0.1741074922801707, 0.19361359650501606, 0.2161235204803589, 0.24246271526221516, 0.2751383033026405, 0.3162837064064915, 0.369235645402662, 0.4408257913958774, 0.5418202837662557, 0.7008977070306374, 0.9914450252680916, 1.710699176539073, 6.376219945510411, 35.21126927155089, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.21407399665430032, 0.07272354967452557, 0.04578829623047529, 0.0376382803462929, 0.03508005328587472, 0.03504977707849241, 0.03669240506776881, 0.03906038156967019, 0.04219171561882855, 0.04597939642755021, 0.04985220741571528, 0.054454144697553906, 0.05929082780163272, 0.06494959299453266, 0.07100014274932817, 0.07729255754314958, 0.08448143165218135, 0.09198801357932859, 0.1001412908027922, 0.109260795045065, 0.11980042219786104, 0.13100121834835385, 0.14382409236715624, 0.157987628258901, 0.1741074922801707, 0.19361359650501606, 0.2161235204803589, 0.24246271526221516, 0.2751383033026405, 0.3162837064064915, 0.369235645402662, 0.4408257913958774, 0.5418202837662557, 0.7008977070306374, 0.9914450252680916, 1.710699176539073, 6.376219945510411, 35.21126927155089, 0.0, 0.0, 0.0], ], '_alpha0.yoda' : [ [0.0, 0.0, 0.0, 0.0, 0.001724063, 0.001004696, 0.0007971206, 0.000722206, 0.0006973311, 0.0006984975, 0.0007135192, 0.0007367059, 0.0007662644, 0.0007979922, 0.0008331413, 0.0008699705, 0.0009095695, 0.0009505057, 0.0009926733, 0.001036868, 0.001085337, 0.001131808, 0.001183161, 0.001234698, 0.001288514, 0.001348953, 0.001412434, 0.001481665, 0.001557491, 0.001640714, 0.001732363, 0.001836518, 0.001956892, 0.002097201, 0.002265482, 0.002474615, 0.002745257, 0.003123757, 0.00371425, 0.004877381, 0.009418508, 0.02213126, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.001724063, 0.001004696, 0.0007971206, 0.000722206, 0.0006973311, 0.0006984975, 0.0007135192, 0.0007367059, 0.0007662644, 0.0007979922, 0.0008331413, 0.0008699705, 0.0009095695, 0.0009505057, 0.0009926733, 0.001036868, 0.001085337, 0.001131808, 0.001183161, 0.001234698, 0.001288514, 0.001348953, 0.001412434, 0.001481665, 0.001557491, 0.001640714, 0.001732363, 0.001836518, 0.001956892, 0.002097201, 0.002265482, 0.002474615, 0.002745257, 0.003123757, 0.00371425, 0.004877381, 0.009418508, 0.02213126, 0.0, 0.0, 0.0], ], } variation_yvals = { } # lists for ratio plot ratio0_yvals = { 'alpha0.yoda' : [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0], '_alpha0.yoda' : [1.0, 1.0, 1.0, 1.0, 0.9982987929520383, 0.9981461603108203, 0.9981745928756426, 0.9967863828256093, 0.9972222091448094, 1.001525714693523, 0.9981942118559849, 0.9997387194355166, 1.0010949495946317, 0.9961117122963785, 1.0014756048490885, 0.9996762902441407, 1.0034848869642607, 1.0003838680271888, 0.9981333056082042, 1.0002816827295662, 1.0026902545510064, 1.001351522879271, 1.0051389669091226, 1.0032570457696304, 0.9964843361665058, 0.9987631874263673, 0.9973222342971778, 0.9991369335454119, 1.0017742014589706, 0.9996624707105439, 0.9984023284318633, 1.0001449463027967, 1.0006793481457634, 0.9997898473548971, 0.9993504865164407, 0.9987159657136514, 0.9999987245320383, 1.0008836976521804, 1.0003531588846533, 0.9997117304550847, 1.0001649198723375, 0.9999999993831852, 1.0, 1.0, 1.0], } ratio0_yerrs = { 'alpha0.yoda' : [ [1.0, 1.0, 1.0, 1.0, 0.08452554059547494, 0.08454082648085401, 0.08456281638092931, 0.08456218777921096, 0.0845747010784226, 0.08458319658321763, 0.08457584861321202, 0.084586513203938, 0.08456908963353225, 0.08455522666274978, 0.08455784857445484, 0.08455664159991329, 0.08454599219832819, 0.08454737845474003, 0.08454764430885342, 0.08454341732424459, 0.08454064141303004, 0.08453542714894023, 0.08453138317841548, 0.08453209992963386, 0.08453118049574947, 0.08452996193049131, 0.0845273647209887, 0.08453074406091884, 0.08452839336853332, 0.08452589788906445, 0.08452715878469247, 0.08452480409638995, 0.08452372587436056, 0.0845224019061642, 0.08452139660988908, 0.08452014183431009, 0.08451916801263946, 0.08451821677879573, 0.08451743215797546, 0.0845164725031339, 0.08451565409097353, 0.08451541660645763, 1.0, 1.0, 1.0], [1.0, 1.0, 1.0, 1.0, 0.08452554059547494, 0.08454082648085401, 0.08456281638092931, 0.08456218777921096, 0.0845747010784226, 0.08458319658321763, 0.08457584861321202, 0.084586513203938, 0.08456908963353225, 0.08455522666274978, 0.08455784857445484, 0.08455664159991329, 0.08454599219832819, 0.08454737845474003, 0.08454764430885342, 0.08454341732424459, 0.08454064141303004, 0.08453542714894023, 0.08453138317841548, 0.08453209992963386, 0.08453118049574947, 0.08452996193049131, 0.0845273647209887, 0.08453074406091884, 0.08452839336853332, 0.08452589788906445, 0.08452715878469247, 0.08452480409638995, 0.08452372587436056, 0.0845224019061642, 0.08452139660988908, 0.08452014183431009, 0.08451916801263946, 0.08451821677879573, 0.08451743215797546, 0.0845164725031339, 0.08451565409097353, 0.08451541660645763, 1.0, 1.0, 1.0], ], '_alpha0.yoda' : [ [1.0, 1.0, 1.0, 1.0, 0.0006807335751804817, 0.0011679549551987982, 0.0014721395745315442, 0.0016225852729025627, 0.0016811995368015884, 0.0016856355811643666, 0.0016446589350129425, 0.0015953603327356126, 0.0015359006330064008, 0.001467492324576854, 0.0014131497788865495, 0.0013508941179697159, 0.0012970042534761104, 0.0012373097572459366, 0.001182084793542051, 0.0011341371901844146, 0.0010860976706344042, 0.001040112390817931, 0.0009987312431364085, 0.0009552521988868068, 0.0009091755062883644, 0.0008704266049862651, 0.0008301065690548295, 0.0007927598273313845, 0.000756154776522184, 0.00071628659625444, 0.0006775371882165586, 0.0006402276077854548, 0.0006011660353658284, 0.000560447668373382, 0.0005185894293216152, 0.0004744613742381573, 0.0004282357906721932, 0.00037668031817193564, 0.00031662761362676175, 0.00024096523972599954, 0.00012484063771067843, 5.312028500026429e-05, 1.0, 1.0, 1.0], [1.0, 1.0, 1.0, 1.0, 0.0006807335751804817, 0.0011679549551987982, 0.0014721395745315442, 0.0016225852729025627, 0.0016811995368015884, 0.0016856355811643666, 0.0016446589350129425, 0.0015953603327356126, 0.0015359006330064008, 0.001467492324576854, 0.0014131497788865495, 0.0013508941179697159, 0.0012970042534761104, 0.0012373097572459366, 0.001182084793542051, 0.0011341371901844146, 0.0010860976706344042, 0.001040112390817931, 0.0009987312431364085, 0.0009552521988868068, 0.0009091755062883644, 0.0008704266049862651, 0.0008301065690548295, 0.0007927598273313845, 0.000756154776522184, 0.00071628659625444, 0.0006775371882165586, 0.0006402276077854548, 0.0006011660353658284, 0.000560447668373382, 0.0005185894293216152, 0.0004744613742381573, 0.0004282357906721932, 0.00037668031817193564, 0.00031662761362676175, 0.00024096523972599954, 0.00012484063771067843, 5.312028500026429e-05, 1.0, 1.0, 1.0], ], } ratio0_variation_vals = { } ratio_band_edges = { }