import numpy as np from numpy import nan add_legend_handle = [ 'alpha2.yoda', '_alpha2.yoda' ] xpoints = { 'alpha2.yoda' : [-9.271111000000001, -7.813333, -6.3555555, -4.897778000000001, -3.44, -1.982222, -0.5244443000000001, 0.9333332000000001, 2.391111, 3.8488890000000002, 5.306667, 6.7644445, 8.222221999999999, 9.680000499999998, 11.13778, 12.595555, 14.053329999999999, 15.511109999999999, 16.968890000000002, 18.42667, 19.884445, 21.34222, 22.8, 24.25778, 25.715555000000002, 27.17333, 28.63111, 30.08889, 31.54667, 33.004445000000004, 34.46222, 35.92, 37.37778, 38.835555, 40.29333, 41.75111, 43.20889, 44.666669999999996, 46.124445, 47.58222, 49.04, 50.49778, 51.955555000000004, 53.41333, 54.87111], '_alpha2.yoda' : [-9.271111000000001, -7.813333, -6.3555555, -4.897778000000001, -3.44, -1.982222, -0.5244443000000001, 0.9333332000000001, 2.391111, 3.8488890000000002, 5.306667, 6.7644445, 8.222221999999999, 9.680000499999998, 11.13778, 12.595555, 14.053329999999999, 15.511109999999999, 16.968890000000002, 18.42667, 19.884445, 21.34222, 22.8, 24.25778, 25.715555000000002, 27.17333, 28.63111, 30.08889, 31.54667, 33.004445000000004, 34.46222, 35.92, 37.37778, 38.835555, 40.29333, 41.75111, 43.20889, 44.666669999999996, 46.124445, 47.58222, 49.04, 50.49778, 51.955555000000004, 53.41333, 54.87111], } xedges = { 'alpha2.yoda' : [-10.0, -8.542222, -7.084444, -5.626667, -4.168889, -2.711111, -1.253333, 0.20444439999999997, 1.662222, 3.12, 4.577778, 6.035556, 7.493333, 8.951111, 10.40889, 11.86667, 13.32444, 14.78222, 16.24, 17.69778, 19.15556, 20.61333, 22.07111, 23.52889, 24.98667, 26.44444, 27.90222, 29.36, 30.81778, 32.27556, 33.73333, 35.19111, 36.64889, 38.10667, 39.56444, 41.02222, 42.48, 43.93778, 45.39556, 46.85333, 48.31111, 49.76889, 51.22667, 52.68444, 54.14222, 55.6], '_alpha2.yoda' : [-10.0, -8.542222, -7.084444, -5.626667, -4.168889, -2.711111, -1.253333, 0.20444439999999997, 1.662222, 3.12, 4.577778, 6.035556, 7.493333, 8.951111, 10.40889, 11.86667, 13.32444, 14.78222, 16.24, 17.69778, 19.15556, 20.61333, 22.07111, 23.52889, 24.98667, 26.44444, 27.90222, 29.36, 30.81778, 32.27556, 33.73333, 35.19111, 36.64889, 38.10667, 39.56444, 41.02222, 42.48, 43.93778, 45.39556, 46.85333, 48.31111, 49.76889, 51.22667, 52.68444, 54.14222, 55.6], } ref_xerrs = [ [abs(xpoints['alpha2.yoda'][i] - xedges['alpha2.yoda'][i]) for i in range(len(xpoints['alpha2.yoda']))], [abs(xedges['alpha2.yoda'][i+1] - xpoints['alpha2.yoda'][i]) for i in range(len(xpoints['alpha2.yoda']))] ] yvals = { 'alpha2.yoda' : [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 774.8899111757391, 121.33394010170002, 32.06335944156106, 18.53795982653051, 12.765873816177777, 9.549780247596168, 7.501272484562125, 6.081467081087049, 5.031587070751417, 4.246489501087277, 3.628286161149141, 3.1278951556476335, 2.7114983056428206, 2.3732181810698467, 2.079562619617635, 1.8388762364691522, 1.6233382266185574, 1.4357084059323084, 1.2816301611365304, 1.130883946823252, 1.0043998408539014, 0.8845696881559635, 0.7812090987666179, 0.6877484102430409, 0.5995205723771765, 0.514792149707089, 0.4441524098286437, 0.37191340197699313, 0.30494066320020863, 0.2368784041487742, 0.17036692779431745, 0.094346540630273, 0.0071135364289291346, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], '_alpha2.yoda' : [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 774.9692, 121.323, 32.08979, 18.51516, 12.7485, 9.53548, 7.499842, 6.074973, 5.028292, 4.24147, 3.621975, 3.120045, 2.714622, 2.372609, 2.084361, 1.836855, 1.625827, 1.438208, 1.275633, 1.132048, 1.00417, 0.8853146, 0.782019, 0.6901463, 0.5997965, 0.5200745, 0.4448663, 0.3723675, 0.3028513, 0.2373025, 0.1694153, 0.09482274, 0.007221134, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], } xerrs = { 'alpha2.yoda' : [ [0.7288890000000006, 0.7288889999999997, 0.7288885, 0.7288890000000006, 0.7288890000000001, 0.7288889999999999, 0.7288887, 0.7288888000000001, 0.7288890000000001, 0.7288890000000001, 0.7288889999999997, 0.7288885, 0.7288890000000006, 0.7288895000000011, 0.7288899999999998, 0.728885, 0.7288900000000016, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000018, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288900000000069, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998], [0.7288889999999988, 0.7288890000000006, 0.7288885, 0.7288889999999997, 0.7288890000000001, 0.7288889999999999, 0.7288887000000001, 0.7288888, 0.7288889999999999, 0.7288890000000001, 0.7288889999999997, 0.7288885, 0.7288889999999988, 0.7288894999999993, 0.7288899999999998, 0.728885, 0.7288899999999998, 0.728889999999998, 0.7288900000000034, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000018, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000053, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000053, 0.7288899999999998, 0.7288899999999998], ], '_alpha2.yoda' : [ [0.7288890000000006, 0.7288889999999997, 0.7288885, 0.7288890000000006, 0.7288890000000001, 0.7288889999999999, 0.7288887, 0.7288888000000001, 0.7288890000000001, 0.7288890000000001, 0.7288889999999997, 0.7288885, 0.7288890000000006, 0.7288895000000011, 0.7288899999999998, 0.728885, 0.7288900000000016, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000018, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288900000000069, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998], [0.7288889999999988, 0.7288890000000006, 0.7288885, 0.7288889999999997, 0.7288890000000001, 0.7288889999999999, 0.7288887000000001, 0.7288888, 0.7288889999999999, 0.7288890000000001, 0.7288889999999997, 0.7288885, 0.7288889999999988, 0.7288894999999993, 0.7288899999999998, 0.728885, 0.7288899999999998, 0.728889999999998, 0.7288900000000034, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000018, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000053, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288849999999982, 0.7288899999999998, 0.7288899999999998, 0.7288899999999998, 0.7288850000000053, 0.7288899999999998, 0.7288899999999998], ], } yerrs = { 'alpha2.yoda' : [ [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 65.4901626835493, 10.254623085191456, 2.7098941182220226, 1.5667950841716012, 1.07894963060142, 0.8071471540843923, 0.6340103511104622, 0.5140145097433302, 0.42529045468427407, 0.3589302869757196, 0.3066828794517672, 0.26440871552629797, 0.22920534123708444, 0.2006239774941092, 0.1757988747708189, 0.15544790573193534, 0.1372431454954325, 0.12138803999063782, 0.10836134235494697, 0.09562706861333001, 0.08492844385213841, 0.07480507801453508, 0.06606568283662693, 0.05816465920154102, 0.0507118220122281, 0.043558785261313036, 0.03757911163627051, 0.03148081367936391, 0.025810981436328467, 0.0200642283191673, 0.014445446779228653, 0.00802538936797656, 0.0006430815830673818, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 65.4901626835493, 10.254623085191456, 2.7098941182220226, 1.5667950841716012, 1.07894963060142, 0.8071471540843923, 0.6340103511104622, 0.5140145097433302, 0.42529045468427407, 0.3589302869757196, 0.3066828794517672, 0.26440871552629797, 0.22920534123708444, 0.2006239774941092, 0.1757988747708189, 0.15544790573193534, 0.1372431454954325, 0.12138803999063782, 0.10836134235494697, 0.09562706861333001, 0.08492844385213841, 0.07480507801453508, 0.06606568283662693, 0.05816465920154102, 0.0507118220122281, 0.043558785261313036, 0.03757911163627051, 0.03148081367936391, 0.025810981436328467, 0.0200642283191673, 0.014445446779228653, 0.00802538936797656, 0.0006430815830673818, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], ], '_alpha2.yoda' : [ [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03974849, 0.01572715, 0.008088386, 0.006143872, 0.005098096, 0.004409097, 0.00391025, 0.003519254, 0.003201759, 0.002940604, 0.002717385, 0.00252208, 0.002352519, 0.002199336, 0.002061413, 0.001935155, 0.001820605, 0.001712337, 0.001612654, 0.001519186, 0.00143081, 0.001343468, 0.001262662, 0.001186176, 0.00110581, 0.001029701, 0.0009523425, 0.0008712931, 0.0007857658, 0.000695552, 0.0005876987, 0.0004396779, 0.0001213336, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03974849, 0.01572715, 0.008088386, 0.006143872, 0.005098096, 0.004409097, 0.00391025, 0.003519254, 0.003201759, 0.002940604, 0.002717385, 0.00252208, 0.002352519, 0.002199336, 0.002061413, 0.001935155, 0.001820605, 0.001712337, 0.001612654, 0.001519186, 0.00143081, 0.001343468, 0.001262662, 0.001186176, 0.00110581, 0.001029701, 0.0009523425, 0.0008712931, 0.0007857658, 0.000695552, 0.0005876987, 0.0004396779, 0.0001213336, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], ], } variation_yvals = { } # lists for ratio plot ratio0_yvals = { 'alpha2.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], '_alpha2.yoda' : [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0001023226952852, 0.9999098347775498, 1.0008243228064455, 0.9987701005534666, 0.9986390421503494, 0.9985025574175105, 0.9998093010799076, 0.9989321522256948, 0.9993451229790753, 0.9988179645596693, 0.9982605668712905, 0.9974902753266971, 1.001152017816378, 0.9997433101285395, 1.0023073988429583, 0.9989008306110729, 1.0015331206649563, 1.0017410179235304, 0.9953206772761869, 1.0010293303570341, 0.9997711659793713, 1.0008421177596412, 1.0010367278551424, 1.0034865798615393, 1.0004602471300181, 1.0102611321791068, 1.0016073090127593, 1.0012209778421348, 0.9931482958740833, 1.0017903525344567, 0.9944142457304487, 1.0050473431939932, 1.0151257496388562, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0], } ratio0_yerrs = { 'alpha2.yoda' : [ [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.08451544114722205, 0.08451570167915265, 0.08451684930773076, 0.08451820474490888, 0.08451827474858024, 0.08451997147133979, 0.08452037336535598, 0.08452146544406702, 0.08452411708355097, 0.08452400197476495, 0.08452554893151962, 0.08453247387428876, 0.08453088123274569, 0.08453667643978183, 0.08453646603973997, 0.08453418596044979, 0.08454377728867535, 0.08454922983599296, 0.08454961941505322, 0.0845595773836515, 0.08455640910888194, 0.084566630550589, 0.08456850149458348, 0.08457258255382433, 0.0845872924946498, 0.08461431528452308, 0.0846085956187173, 0.08464554789373074, 0.0846426356047579, 0.0847026489868014, 0.08479020527193235, 0.08506288958093977, 0.09040251490835349, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0], [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.08451544114722205, 0.08451570167915265, 0.08451684930773076, 0.08451820474490888, 0.08451827474858024, 0.08451997147133979, 0.08452037336535598, 0.08452146544406702, 0.08452411708355097, 0.08452400197476495, 0.08452554893151962, 0.08453247387428876, 0.08453088123274569, 0.08453667643978183, 0.08453646603973997, 0.08453418596044979, 0.08454377728867535, 0.08454922983599296, 0.08454961941505322, 0.0845595773836515, 0.08455640910888194, 0.084566630550589, 0.08456850149458348, 0.08457258255382433, 0.0845872924946498, 0.08461431528452308, 0.0846085956187173, 0.08464554789373074, 0.0846426356047579, 0.0847026489868014, 0.08479020527193235, 0.08506288958093977, 0.09040251490835349, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0], ], '_alpha2.yoda' : [ [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 5.129566074707268e-05, 0.00012961871992962375, 0.00025226258698006856, 0.0003314211519224045, 0.00039935346952429913, 0.00046169617370094354, 0.0005212782242009509, 0.0005786850365341352, 0.0006363318283035992, 0.0006924788108500171, 0.0007489445096963789, 0.000806318586301145, 0.0008676085082200644, 0.0009267314811352656, 0.000991272386103491, 0.0010523573917708099, 0.0011215192066242122, 0.0011926774217693994, 0.0012582834337871095, 0.0013433615396765699, 0.0014245422408505973, 0.0015187814120114023, 0.0016162919786693543, 0.0017247237250331434, 0.0018444904994924871, 0.0020002266945715634, 0.002144179517943894, 0.0023427311179657325, 0.002576782616505644, 0.0029363250841691357, 0.0034496055520207756, 0.004660244001134266, 0.017056720129605837, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0], [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 5.129566074707268e-05, 0.00012961871992962375, 0.00025226258698006856, 0.0003314211519224045, 0.00039935346952429913, 0.00046169617370094354, 0.0005212782242009509, 0.0005786850365341352, 0.0006363318283035992, 0.0006924788108500171, 0.0007489445096963789, 0.000806318586301145, 0.0008676085082200644, 0.0009267314811352656, 0.000991272386103491, 0.0010523573917708099, 0.0011215192066242122, 0.0011926774217693994, 0.0012582834337871095, 0.0013433615396765699, 0.0014245422408505973, 0.0015187814120114023, 0.0016162919786693543, 0.0017247237250331434, 0.0018444904994924871, 0.0020002266945715634, 0.002144179517943894, 0.0023427311179657325, 0.002576782616505644, 0.0029363250841691357, 0.0034496055520207756, 0.004660244001134266, 0.017056720129605837, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0], ], } ratio0_variation_vals = { } ratio_band_edges = { }