import numpy as np from numpy import nan add_legend_handle = [ 'alpha0.yoda', '_alpha0.yoda' ] xpoints = { 'alpha0.yoda' : [-8.7644445, -6.2933335, -3.8222224999999996, -1.3511113, 1.1200002, 3.5911115, 6.0622225, 8.5333335, 11.0044445, 13.475555, 15.946665, 18.417775, 20.888885000000002, 23.36, 25.831114999999997, 28.302225, 30.773335000000003, 33.244445, 35.715554999999995, 38.186665000000005, 40.657775, 43.128885, 45.6, 48.071115, 50.542225, 53.013335, 55.484445, 57.955555000000004, 60.426665, 62.897774999999996, 65.368885, 67.84, 70.311115, 72.78222500000001, 75.25333499999999, 77.724445, 80.195555, 82.666665, 85.137775, 87.608885, 90.08000000000001, 92.55111500000001, 95.02222499999999, 97.493335, 99.96444500000001], '_alpha0.yoda' : [-8.7644445, -6.2933335, -3.8222224999999996, -1.3511113, 1.1200002, 3.5911115, 6.0622225, 8.5333335, 11.0044445, 13.475555, 15.946665, 18.417775, 20.888885000000002, 23.36, 25.831114999999997, 28.302225, 30.773335000000003, 33.244445, 35.715554999999995, 38.186665000000005, 40.657775, 43.128885, 45.6, 48.071115, 50.542225, 53.013335, 55.484445, 57.955555000000004, 60.426665, 62.897774999999996, 65.368885, 67.84, 70.311115, 72.78222500000001, 75.25333499999999, 77.724445, 80.195555, 82.666665, 85.137775, 87.608885, 90.08000000000001, 92.55111500000001, 95.02222499999999, 97.493335, 99.96444500000001], } xedges = { 'alpha0.yoda' : [-10.0, -7.528889, -5.057778, -2.586667, -0.11555559999999998, 2.355556, 4.826667, 7.297778, 9.768889, 12.24, 14.71111, 17.18222, 19.65333, 22.12444, 24.59556, 27.06667, 29.53778, 32.00889, 34.48, 36.95111, 39.42222, 41.89333, 44.36444, 46.83556, 49.30667, 51.77778, 54.24889, 56.72, 59.19111, 61.66222, 64.13333, 66.60444, 69.07556, 71.54667, 74.01778, 76.48889, 78.96, 81.43111, 83.90222, 86.37333, 88.84444, 91.31556, 93.78667, 96.25778, 98.72889, 101.2], '_alpha0.yoda' : [-10.0, -7.528889, -5.057778, -2.586667, -0.11555559999999998, 2.355556, 4.826667, 7.297778, 9.768889, 12.24, 14.71111, 17.18222, 19.65333, 22.12444, 24.59556, 27.06667, 29.53778, 32.00889, 34.48, 36.95111, 39.42222, 41.89333, 44.36444, 46.83556, 49.30667, 51.77778, 54.24889, 56.72, 59.19111, 61.66222, 64.13333, 66.60444, 69.07556, 71.54667, 74.01778, 76.48889, 78.96, 81.43111, 83.90222, 86.37333, 88.84444, 91.31556, 93.78667, 96.25778, 98.72889, 101.2], } 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, 694.8152402343949, 29.510353844890016, 17.542813738435864, 12.602048228509364, 9.832917258674335, 8.086325578383805, 6.868868646074028, 5.968750885229716, 5.268344994759442, 4.706800964744732, 4.253987074634477, 3.8680528183690686, 3.5368805921225692, 3.2433372047379567, 2.985939921735573, 2.7365787034976155, 2.497902562006552, 2.219123389893608, 1.0632223445239408, 0.18628527261028477, 0.13072384474993004, 0.09770819591195846, 0.07541465171522121, 0.05822666736810576, 0.04689645543905378, 0.03690407144967237, 0.02958012391192627, 0.022903319142736904, 0.018130868314239275, 0.013492244376009184, 0.010056654701733246, 0.00711657918910936, 0.004907984670856397, 0.0031052450113511795, 0.001628130273439874, 0.0007361979029666799, 5.1911198161157715e-05, 0.0, 0.0, 0.0, 0.0], '_alpha0.yoda' : [0.0, 0.0, 0.0, 0.0, 694.8691, 29.51814, 17.5317, 12.5871, 9.841658, 8.086543, 6.863058, 5.9545, 5.262598, 4.704198, 4.250012, 3.864797, 3.528831, 3.244839, 2.979187, 2.739133, 2.495132, 2.225684, 1.06552, 0.1865824, 0.129779, 0.09751308, 0.07454638, 0.05908119, 0.04698766, 0.0367806, 0.02904225, 0.02298808, 0.01781876, 0.01351292, 0.01012976, 0.007284411, 0.005024582, 0.003014749, 0.001679247, 0.000664462, 6.578831e-05, 0.0, 0.0, 0.0, 0.0], } xerrs = { 'alpha0.yoda' : [ [1.2355554999999994, 1.2355555000000003, 1.2355554999999998, 1.2355557, 1.2355558, 1.2355554999999998, 1.2355555000000003, 1.2355555000000003, 1.2355555000000003, 1.2355549999999997, 1.2355550000000015, 1.2355550000000015, 1.235554999999998, 1.2355599999999995, 1.2355550000000015, 1.2355550000000015, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.2355599999999995, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.2355549999999909, 1.2355599999999924, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.2355599999999924, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355549999999909], [1.2355555000000003, 1.2355555000000003, 1.2355555000000003, 1.2355557, 1.2355558, 1.2355554999999998, 1.2355555000000003, 1.2355554999999994, 1.2355555000000003, 1.2355549999999997, 1.2355549999999997, 1.235554999999998, 1.2355550000000015, 1.2355599999999995, 1.235554999999998, 1.2355550000000015, 1.2355550000000015, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.2355599999999995, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.2355600000000067, 1.235555000000005, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355600000000067, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005], ], '_alpha0.yoda' : [ [1.2355554999999994, 1.2355555000000003, 1.2355554999999998, 1.2355557, 1.2355558, 1.2355554999999998, 1.2355555000000003, 1.2355555000000003, 1.2355555000000003, 1.2355549999999997, 1.2355550000000015, 1.2355550000000015, 1.235554999999998, 1.2355599999999995, 1.2355550000000015, 1.2355550000000015, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.2355599999999995, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.2355549999999909, 1.2355599999999924, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.2355599999999924, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355549999999909], [1.2355555000000003, 1.2355555000000003, 1.2355555000000003, 1.2355557, 1.2355558, 1.2355554999999998, 1.2355555000000003, 1.2355554999999994, 1.2355555000000003, 1.2355549999999997, 1.2355549999999997, 1.235554999999998, 1.2355550000000015, 1.2355599999999995, 1.235554999999998, 1.2355550000000015, 1.2355550000000015, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.2355599999999995, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.235554999999998, 1.235554999999998, 1.235555000000005, 1.2355600000000067, 1.235555000000005, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005, 1.2355600000000067, 1.235555000000005, 1.2355549999999909, 1.235555000000005, 1.235555000000005], ], } yerrs = { 'alpha0.yoda' : [ [0.0, 0.0, 0.0, 0.0, 58.72261317154907, 2.4941049602724314, 1.4826700401205029, 1.065097736977935, 0.8310635033999701, 0.6834491503558526, 0.5805565568567868, 0.5044794613361472, 0.4452888253149391, 0.3978348503403543, 0.3595630679555053, 0.32694009773089483, 0.29895577181435007, 0.27414381294053425, 0.25239433121826815, 0.23131239523363478, 0.21114723825155182, 0.1875842274522165, 0.08989460691236732, 0.015772852859989413, 0.011074034949852905, 0.00828677756741927, 0.00639507863426902, 0.004946298407903836, 0.003993985085787852, 0.003147356287995667, 0.002534173246865579, 0.0019669543414642713, 0.00155841575197379, 0.0011673078352676016, 0.0008792321878576671, 0.0006293745228084631, 0.0004404322588720385, 0.0002945638673988023, 0.0001627274868927708, 8.624685101482389e-05, 1.565180967302425e-05, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 58.72261317154907, 2.4941049602724314, 1.4826700401205029, 1.065097736977935, 0.8310635033999701, 0.6834491503558526, 0.5805565568567868, 0.5044794613361472, 0.4452888253149391, 0.3978348503403543, 0.3595630679555053, 0.32694009773089483, 0.29895577181435007, 0.27414381294053425, 0.25239433121826815, 0.23131239523363478, 0.21114723825155182, 0.1875842274522165, 0.08989460691236732, 0.015772852859989413, 0.011074034949852905, 0.00828677756741927, 0.00639507863426902, 0.004946298407903836, 0.003993985085787852, 0.003147356287995667, 0.002534173246865579, 0.0019669543414642713, 0.00155841575197379, 0.0011673078352676016, 0.0008792321878576671, 0.0006293745228084631, 0.0004404322588720385, 0.0002945638673988023, 0.0001627274868927708, 8.624685101482389e-05, 1.565180967302425e-05, 0.0, 0.0, 0.0, 0.0], ], '_alpha0.yoda' : [ [0.0, 0.0, 0.0, 0.0, 0.03380616, 0.006967691, 0.005369779, 0.004549955, 0.004023264, 0.003646917, 0.003359721, 0.003129443, 0.002942012, 0.002781552, 0.002643866, 0.002521202, 0.002409128, 0.002310154, 0.00221357, 0.002122516, 0.002025775, 0.00191327, 0.001324033, 0.000553962, 0.0004620048, 0.0004004754, 0.0003501528, 0.0003117231, 0.0002779945, 0.0002459539, 0.0002185544, 0.0001944445, 0.0001711919, 0.0001490799, 0.0001290755, 0.0001094565, 9.090638e-05, 7.041578e-05, 5.25535e-05, 3.305822e-05, 1.040205e-05, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.03380616, 0.006967691, 0.005369779, 0.004549955, 0.004023264, 0.003646917, 0.003359721, 0.003129443, 0.002942012, 0.002781552, 0.002643866, 0.002521202, 0.002409128, 0.002310154, 0.00221357, 0.002122516, 0.002025775, 0.00191327, 0.001324033, 0.000553962, 0.0004620048, 0.0004004754, 0.0003501528, 0.0003117231, 0.0002779945, 0.0002459539, 0.0002185544, 0.0001944445, 0.0001711919, 0.0001490799, 0.0001290755, 0.0001094565, 9.090638e-05, 7.041578e-05, 5.25535e-05, 3.305822e-05, 1.040205e-05, 0.0, 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, 1.000077516672759, 1.0002638448576695, 0.9993664791406003, 0.9988138254799289, 1.0008889265612353, 1.000026887566434, 0.9991540606796507, 0.9976124174883926, 0.998909146085695, 0.9994469779444193, 0.9990655649477218, 0.9991582797541939, 0.9977240984215023, 1.0004630401241812, 0.9977384267893616, 1.0009333904773576, 0.9988908446435452, 1.0029563971684814, 1.0021610300873502, 1.001595012775577, 0.9927722080716224, 0.9980030752780015, 0.9884866972733632, 1.0146757949667975, 1.0019448071307808, 0.9966542594130634, 0.9818163739432677, 1.003700811080475, 0.9827858043624884, 1.0015324080571486, 1.0072693455661905, 1.0235832141301087, 1.023756661229192, 0.9708570463778631, 1.031395968365681, 0.9025589414509285, 1.2673240520428934, 1.0, 1.0, 1.0, 1.0], } ratio0_yerrs = { 'alpha0.yoda' : [ [1.0, 1.0, 1.0, 1.0, 0.08451543629316347, 0.08451626752365453, 0.08451723094294789, 0.0845178273931999, 0.08451850875352665, 0.08451912351672242, 0.08451996780992599, 0.08452010664149733, 0.0845215766541256, 0.0845234063050997, 0.0845237800790452, 0.08452317304931382, 0.08452526570452845, 0.08452522683736288, 0.08452759862347275, 0.08452612561005277, 0.08452981371776903, 0.08453077837245, 0.08454920776953549, 0.08467042315785621, 0.084713198047664, 0.08481148884261668, 0.08479888839662013, 0.08494902132443219, 0.08516603330455112, 0.08528479824476409, 0.08567148854450328, 0.08588075506462266, 0.08595372957123466, 0.08651695023722027, 0.0874278986337408, 0.08843778816817036, 0.08973790433522018, 0.09486010486194431, 0.09994746092950174, 0.11715172057305819, 0.3015112389514376, 1.0, 1.0, 1.0, 1.0], [1.0, 1.0, 1.0, 1.0, 0.08451543629316347, 0.08451626752365453, 0.08451723094294789, 0.0845178273931999, 0.08451850875352665, 0.08451912351672242, 0.08451996780992599, 0.08452010664149733, 0.0845215766541256, 0.0845234063050997, 0.0845237800790452, 0.08452317304931382, 0.08452526570452845, 0.08452522683736288, 0.08452759862347275, 0.08452612561005277, 0.08452981371776903, 0.08453077837245, 0.08454920776953549, 0.08467042315785621, 0.084713198047664, 0.08481148884261668, 0.08479888839662013, 0.08494902132443219, 0.08516603330455112, 0.08528479824476409, 0.08567148854450328, 0.08588075506462266, 0.08595372957123466, 0.08651695023722027, 0.0874278986337408, 0.08843778816817036, 0.08973790433522018, 0.09486010486194431, 0.09994746092950174, 0.11715172057305819, 0.3015112389514376, 1.0, 1.0, 1.0, 1.0], ], '_alpha0.yoda' : [ [1.0, 1.0, 1.0, 1.0, 4.865489131843963e-05, 0.00023611004587146007, 0.00030609565147664703, 0.0003610488483694838, 0.0004091628042990787, 0.0004509980416505689, 0.0004891229070045302, 0.000524304508627446, 0.0005584319179792696, 0.0005909644407814585, 0.0006215030637410138, 0.0006518013373620486, 0.0006811448498899486, 0.0007122768476325135, 0.0007413310575630624, 0.0007756093392407157, 0.0008109904008316105, 0.0008621737793912074, 0.0012453020826916853, 0.0029737294432228556, 0.00353420449715045, 0.004098687896774338, 0.004643034105921985, 0.005353613972603032, 0.005927836067723267, 0.006664681980561892, 0.007388555932041992, 0.00848979568368204, 0.009442013313038768, 0.011049303277153932, 0.012834834627239817, 0.015380493505574057, 0.01852214016474051, 0.022676400651992386, 0.03227843671806817, 0.0449039855544063, 0.20038162031450238, 1.0, 1.0, 1.0, 1.0], [1.0, 1.0, 1.0, 1.0, 4.865489131843963e-05, 0.00023611004587146007, 0.00030609565147664703, 0.0003610488483694838, 0.0004091628042990787, 0.0004509980416505689, 0.0004891229070045302, 0.000524304508627446, 0.0005584319179792696, 0.0005909644407814585, 0.0006215030637410138, 0.0006518013373620486, 0.0006811448498899486, 0.0007122768476325135, 0.0007413310575630624, 0.0007756093392407157, 0.0008109904008316105, 0.0008621737793912074, 0.0012453020826916853, 0.0029737294432228556, 0.00353420449715045, 0.004098687896774338, 0.004643034105921985, 0.005353613972603032, 0.005927836067723267, 0.006664681980561892, 0.007388555932041992, 0.00848979568368204, 0.009442013313038768, 0.011049303277153932, 0.012834834627239817, 0.015380493505574057, 0.01852214016474051, 0.022676400651992386, 0.03227843671806817, 0.0449039855544063, 0.20038162031450238, 1.0, 1.0, 1.0, 1.0], ], } ratio0_variation_vals = { } ratio_band_edges = { }