{
"cells": [
{
"cell_type": "markdown",
"id": "3d755697",
"metadata": {},
"source": [
"# ratioplot5\n",
"Example that shows how you can set the colors of the confidence interval bands by using\n",
"the method `TRatioPlot::SetConfidenceIntervalColors`.\n",
"Inspired by the tutorial of Paul Gessinger.\n",
"\n",
"\n",
"\n",
"\n",
"**Author:** Alberto Ferro \n",
"This notebook tutorial was automatically generated with ROOTBOOK-izer from the macro found in the ROOT repository on Tuesday, March 19, 2024 at 07:11 PM."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "28051247",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:41.138113Z",
"iopub.status.busy": "2024-03-19T19:11:41.137738Z",
"iopub.status.idle": "2024-03-19T19:11:44.682603Z",
"shell.execute_reply": "2024-03-19T19:11:44.681458Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3dzY7lRv3/8aq/uBBugkg9PWqXCckKgoTIBcAiI7FMoswIZWIboiiDCDsizWzy2xNGIo1YZADbrXR6EbgcVgj5v6hpj9tPx8/+Vvn90ijp9vE5XfbxsT+nnqyLolAAAABt/t/eBQAAAHIRFAAAQCeCAgAA6ERQAAAAnQgKAACgE0EBAAB0IigAh5Zl2d5FACAaQQE4tDAMex7NsswYs1VZAEhEUACO6P/+7/9ee+21/nWyLCtjxBdffLF+oQBIRFAAjujzzz//7rvvyl+zLIvjuNoMYYyppoRf/vKXG5cQgBCaKZyBA3rttde+++67oii01kEQ2IV5nkdRFMdxuZptd8jzXCnFuQI4JmoUgKMzxmRZlmVZFEVJkpTLv/jii//85z82JQA4LIICcGhffPHF3/72t9dee+211177/ve/X33o888//9e//rVTuQBI8b29CwBgN7YBovz1F7/4RbXT4nfffWebHv7zn/+QGIDDokYBOK7vvvvuiy++KIqiKIpf/epXSqnPP/9870IBkIWgABzXr371q1/84hf25z/+8Y/7FgaATAQF4LgIBwBOoo8CcGj9gx7tzAo/+MEPNioNAHmoUQAAAJ0ICgAAoBNBAQAAdCIoAACATgQFAADQiZtCAQCATtQoAACATgQFAADQiaAAAAA6ERQAAEAnggIAAOhEUAAAAJ0ICgAAoJOsoBDHcetCY0zzoSzL4jhufQoAAFiEoAmXsiwLwzBNU2NMuVBrrZQKgiDPc1W5JW4cx0mSBEGglMrzvPYsAACwjEKANE3tJV8plaZpudwuLH9VSkVRVP5crhkEQRAE1RfccX8CAI5s3evlHr43YS9orYMgsC0CS+1ZY4wxJkmS2vIyQNifsyxTty0U5V+P4zgMw9oTC+ICAGBbthbcM1P6KKRpqpQKw1BrbYyxF+85bBeEZm+DOI7zPI/j2HZHsD+3Pl0pNb8YAACgZkpQsOGgKIooitRtYlijU6ExJgiCJEnCMEySJIqiMhNUaxpa6QEWLzAAAJ6ZNerBfte3cSFJElvBsEy5bl/fdlQsiiJN0yRJhseRIe0uCxYVAAAvTQ8Kto+C1tqOPkjT1DZJLPhNPUmScjiDMSaKItuJwRhjB0FYttGBUQ8AACxuYtODzQfqduhBlmW2N6K9Zq/dXaCWCeidAADASiYGhWo+qD264JQGQRCUwxmyLCsnTijrGOxDtvvCIn8RAABUCZpwSSmlta7ljGorQzk8Ut3OztRcXr6OqO0CAByBl1cfNzapteqih5dvFQBAOC+vPh5ukuroUOnllgIA5PAyKEyZmdEJ/r1VAABsT9bdIwEAgCje1igAaKX1s5PrFMU7G5QEgBMICsDh9OeAIUkCwHHQ9AAAADoRFAAAQCeCAgAA6ORtH4XmVAoMmAQAYCxvgwKxAACA+Wh6AAAAnQgKAACgE0EBAAB0IigAAIBOBAUAANCJoAAAADp5OzySeRQAAJjP26BALAAAYD5ZTQ9xHDcXZlkWx3Hzoa7lAABgKYKCQpZlSZJkWVZdaIwJw9A+VG1NiOPYLs+yTGtdexYAAFiEllBFb+sG8jxXSqVpaowpl4dhWJZQax0Egc0EWutyTfvfalbQWsR2AQJp/awo3ulf4eSL9L8CcFheXn2k9FEwxhhjkiSpLozjOAiC8tdy79vmhjJP2NqFbcoJeO9kCBiSJAB4Q0RQsClBKVULCnmeR1Fk2xdURw8GValRKKMDAABYhKA+Cq2SJKn2UbCJIcuyak1DKz3AFhsAAIDLRNQo9CtbHGzHxoHNP/61EgEr6srNfI6Aw5NeoxBFUflz2bJgjLE9Hy1bzUC7AzCF1p0p4eSjAA5AdFDoal+oZQLGRgJT1EJAUdT/da0J4EhEB4U4jsuZFWw3BVvBUB0VqZQqlwMYqhkRmppxAcDxiO6jYIyJoqgc+hgEQTnwIU3TMAxth8TqcgCnlZf8IV0Q7Dr2KS+f+HSdYgGQyI2pIcYOffRyygtgAbWKhM2eCxyDl1cfDzdJefpWAXMtcqUf8CLM7YjD8vLqI7rpYQ5uMw3ccfuJ0OrprIt0UbzKClp3ZYX5s0QDEEJ0Z8Y5ioa9SwTsZ1SnhJMvpp6+eh16OAK+87ZGAcBLi6aE25d8ptTTQj1QVA8AviMoAF5pXrZtOtDqqVroiv6qWUE/UEoV6gF9GwGPERQA39zpH3BbnbBK58Gyv0J3ZwUArvO2jwKANRod6uisAPiOoAB4aoOUsM3rA9gVQQHw2pZXcSoVAB9520eBeRRwaBtfs+msAPjL26BALMBxbdboUFWdiAmAR2h6ADy1V1YmLgB+ISgAXrGTIO30txkBAXiIoAD4aK/qBJr8AO8QFACPyPkqL6ckAOYhKAC+2KUPYxMNEIBfvB31wPBIHJNWT/c/0BkBAXjE26BALMCxyLwwM60C4D6aHgCPyLkqyykJgHlkBYU4jnseyrKsuiTLsjiOe54CHIXM6gRLctkADCAoKGRZliRJLQ10PRTHcRiGWZZlWaa1bn0WcCzSvsRLKw+ASUQEhSzLjDFhGHat0HwoSZI0TW1QCIKAegUcl/yv7PJLCKCbiKCglDLGRFHU9VAQBEEQlEtsLDDGlL/meb56EQHJZH59l1kqAGOICArGmK7eBjYE9Lcs2MRA6wOOyJUv666UE0CD9OGRtomhttA2N/Q/sTmPQhNDKOEJyUcycyoAjhMdFLTWURSVTQyjEALgP7cuwMypALhJRNNDK9uUYPs5GmPyPE+SxIYG+2ttzWl5AnCe/Kuv/BIC6Ca3RqHWvTHP8yAIyqCQJEn5EL0TcERuVSdYVCoADpIbFNTd+Zds1UJ1vIMxxkaEJEm6RkwAnnPluktPBcBZcpse+qVpmue51lprzTwKOBwuugC2IqtGoacHYq19wRhDd0UcnYsfAVofANe4WqMAHJej1QnkA8BN3gYF3bB3iYBFuXvd5cMIOEVW08OCaJiAn5y+ytKlEXCQtzUKgM8cz8GFerB3EQAMRVAA3OHB13HHIw5wQAQFwDV+XGs9CD3AMRAUAEd4c2X1I+gAh0FQAJzi01XWm+gDeI2gALjAr2uqVk/3LgKAobwdHtmcOIEBk3Cef8cwEzUC4nkbFIgF8Idf1QkvMacC4AiaHgBH+Jp9iQuAbAQFADvxNfoAfiEoALId4Qv3EbYRcBZBAXCBr1++fd0uwCMEBUCw43zVPs6WAq7xdtQDwyPhIq2fVX+1h6xWT1VleVG8s22hVsbwB0A2b4MCsQCOepUDbi+f1WRQSxIAsDZvgwLgiUbk9TYrMPkSIJKsoBDHcRzHzYVZlhljag9lWZZlmV1hm+IBu/Ot3cGi9QEQTFBnxizLkiSx1/6S1jpJEmNMkiRa6/LROI7DMLRZoboc8MQxL5zH3GpANhFBwVYYhGFYW26MUUoVRRHHcVEUQRCU6yRJkqapDQpBEFCpAD8dpyr+OFsKuEZEUFBKGWOiKKotzPO8utDmBnXb1lD9Nc/z9csIbOXIX6yPvO2ASCL6KBhj7FU/SZLq8trIha72BftcWy2xSvmAXRztSzY9FQCRRASFIYwxeZ6naaqUss0N/es351FoYgglJDrGxbJ17EZRecjPbpuAgxwICnEc25qGNE2H1xkQAuA2rw/gzhCgHyilCvVAq6ebFghANyl9FLrY8Q5RFBVFUaYEW7tQrmObJGh3AJzndTwCHCW6RsH2UmzWDdj0UP7K2Eh4o1AP9i4CANyhRVXRa62r7Qta6+bQR/uofchGBK11FEXV1bSWtV3AULaDwpGP3rKLxpF3Apzl5dVH1iY1g0JzHVvgLMvKORXKxFB9HVHbBQxFUFDsBDjMy6uPh5ukPH2r4D++TFsEBTjLy6uP9M6MwOF4d5aZ6BjDRAH5RHdmnKPZbOFfyoNXuC6WmHkJkMTboEAsgIu0esqB+wo3ngYEoOkBgDhMuATIQVAABKCmHYBUBAVADKrZm4hQwN4ICgBEIjYBMhAUgL3xpRmAYN6OemB4JBzD8dmFsQ/ArrwNCsQCwHlMqAAIQNMDsCsuhEOwl4D9EBQAAagA68KeAfbmbdMDIJDWz6q/FupB63IAkIOgAGyqKN559Yt+YBfxrfk0ujQCO6HpAYBs5ANgVwQFYCd00BuLPQbswdumB+ZRgBs4LIdgnCSwH2+DArEAonHZA+AImh6A/RBnxyJgAZuTFRTiOG5daIxpPpRlWRzHrU8B4BtCFbATQUEhy7IkSbIsqy7UWidJopRKkqTa7SCO4zAMsyzLskxrXXsWIBpfi+dg7wHbEhEUsiwzxoRhWFtuawuKosiyzPY5KOsPkiRJ09QGhSAIqFeAe/iKPBZ7DNiDiKCglDLGRFFUW5gkSRAE5a9BENiaAxsLjDF2eRzHeZ5vUkxgNr4QA3CKiKBguyC01gqUacD+3BoI7Dq0PsAlfDmeg7AFbMjJ4ZG2uaF/neY8Ck0MoQQcw4QKwOacDApDEAIgUHkXKJzUdaOs4vbRO3fNALAaJ4OCMcYOhbBso0O1kQIQjRR7Sl8I0A/Uy8hFUAC2IKKPQpey96JV9m2sZQJ6JwAHQswCtiU6KNjhDDYHNMc7lHEhSZLmiAlAHBrXAThIdNODHTNZzq8QRVEZDtI0DcPQ9lhkHgW4hC/ES9GanQlsQHvZ6U9rP7cLbrM1ChyZi2BnQiQvrz6iaxTm4DbTkIV2BwBu8jYoEAsgkFZPOS6XROsDsD7RnRkBoJVWT/cuAnAUBAVgfbQ7AHAWQQHYCpXkayCEASsjKABwE8EL2ARBAVgZX3kBuIygAGyCr7/rIYoBa/J2eCTzKAD+467TwPq8DQrEAojAZWwbTKgArIamB2B9XMPWw74FVkZQAFZDdQIA9xEUgJXxlXcbxDJgHQQFAI4jigFrIigA6+AL7vbY58AKCArAmviyuw32M7Aab4dHMo8C9sRXWwC+8DYoEAuwPw7C7TGhArA0B4JCdiuOY2NM9aE4jrMsM8bEcbxP4YAKrZ/ZH4q7v2ILzNIIrEML/+Ydx3GSJEEQKKXyPA+CIMsy+5BtXAiCIM9zdbcKQWvp2wUvaf2sKN55dbniIFzNy13dWPryB/Y8duLl1Ud6Z8YkSaIosjUKaZrmeW6Dgq1CKIoiyzL7rlCpAFm8O1k4gH0OrEB6UFBKlc0N1XaHsprBqtY0AHui9huAX6QHhSAIwjC0NQo2KLTmBmOMbYAAROCr7b6Ia8BypHdmzLJMax2Gof01TdOBT2wOj2zyryUJODq6NAJLkx4UtNZls0Icx2EYpmlaG/vQihCA7RXqwd5FOJCuQSXlkJOW3o4AxhMdFGw+KDsf2BEQZRsEIBQhdX19IUA/UC9DG0EBWID0Pgo1Ze1CrfdirW8jgOMiqAGLEh0Uar0X4zjO89wOg7Q/V6scGB6JndE0DsBHopselFJpmoZhWPZMjKKoTA9RFJWdHMvlwM74OisH0zkDS/BwDinl6dxYEI05AaWx7whvB7bl5dVHdNMD4Batnu5dBNxFexAwm7dBQTfsXSIA2yG0AUuR3kdhMv8qfyAXMRSAv7wNCsCyeu4Z/XKGH77CykSXRmAeggIwVOckP3aGH+YBlIbpnIEleNtHAdgIlyIAXiMoAEugclsywhwwA0EBgL8IcMBsBAVgBr6quoJ3CpjK286MzYkTGDCJtXBoSUaXRmAeb4MCsQCr4/LjFsZJApPQ9ADMw7VHPt4jYAaCAjAJ1QkAjoGgAMzAV1W3EO+A8QgKAA6ASAdMRVAAxuOLqbt474CRvB31wPBIrI4jyi2MkwQm8TYoEAuwFi42rmOcJDAGTQ/AJFxpXMS7BoznRlCI49gYE8fxwOXAWqhOAHAwWn4Vve1tEARBnudBEGRZ1lyu7rY1aO3AdsEtWj8rindeBgWOLvFevl9tD7z8gTcRK/Dy6iO9RsEYEwRBURRZlqVpmue5DQq2FsEut+8K9QoATvPuJA6sTXpnxjzP0zS1PxtjyqSWJEkQBOVq1ZoGYC20O/iELo3AMKKDgr32V3shVKsNjDHVn5Mk2bBoODCuLq5jnCQwhuigYGmtbeVBnudJkgxs/mnOo9DkX0sS1lOoB3sXAQB2IL2PglIqiqIsy8b2RSgGWLfc8BKHjU+oVwAGcCAoVJMBfRGwD64oniHwAYOJDgq2F0I1GeR5bhfWEkOtbyOwCq4uAI5Heh+FIAjCMKw2OtigEMdxGIZZlhljqgMmAcDS+lnPo8XtSuQ/oJ8DU0MYY+yUSkqpNE3LwQ5xHJcjHaIoqgYFL6e8wG6YosdXTJ+FpXl59fFwk5SnbxV2o7VSSqun7TP9wV1EQCzNy6uP6D4KwP7oxugvrZ7uXQTAAd4GBd2wd4ngMu++IgDAQNI7M07mX+UPdkC+PAi6NALdvA0KwCitPeTtpUOrp6q3/zwcxnTOwCkEBeClel/F2+sHfRgPgUoFoIO3fRSAZXDx8B5vMdCLoAC0oTr6gHjTgTYEBaAb3zUPgjca6EZQABr4ZnlYvPVAg7edGZsTJzBgEuNwwBwKwx+ADt4GBWIBJuJqcSTVYbHlbaKqMzYy5gXwNigAsxA0D6ARAt6xMbFc3n//SeAg6KMAVFCdAMVhANxBUAAaqE44LN56oIGgAAANVCoAtwgKwC2uDVBUKgB13gYFbjONibhOwOKkASilPB71wPBIjFKoB3sXAWIwpwJQob28oGrt53ZhRfbCwGEDqwwKHBIYw8urj0tND8aY2pI4jo0xcRzvUBr4hK+PqPHuXA9M5kz2McbkeZ6maRkXbLeDIAjyPFd32xq8zHRYEdUJaKJSAeN5efVxo0YhyzKbBkq2FqEoiizL7LtCvQImojoBrbw73QPTuBEUwjCMoqi6JEmSIAjKX4MgyLJs62LBI9Xp/YE7iJI4NgeCgjEmiqJmhUG1y4JtmNiyVPAE1wB0Iz4CSv7wyDiO8zyfUFswZOIE/1qSMFFRKG7/gx53bylZxe0l4T3RQSHLsiRJpl3OCQE4jeoE9CqKd5q3lKzi9pI4AtFBwTY3VJsYwjCkOwIWRqbEEFpzqOCYpAeFaibI8zwIApseanGh1rcROI3qBAxUTtRIVsAhuTTiU2tdzqOQZVkYhvbX6s/lmg5tF/Zxd+4ErZ/R2Iw+bZNtcNigxsurj+gahR52KEQYhvbXKIqa8zYCnahOwDRUKuB4PMw+ytNMh8W0zbjHV0Oc1jhyOGxQ4+XVx4F5FIBVePdhxuo4ZnBIrjY9nNScR8G/lIcpaHTAfDRA4Ei8DQrEAvTh8MA0jIDA8dD0gCOhOgHzkQ9wMN7WKABVWj8r1IOXP6unTNiMBWituBkEDoCggIOhUQrzlQ0QwAF4OJBDeTpABbO0DYkEpuOIQhsvrz70UcABcE7H4jiWcBjeNj0wPBJ1HABYA8Mf4DtvgwKxAC/RloyV0FMBx+BtUACqtHpKcsRatNbdwx+Y4xmuIyjAa3zhw6puKxUK9aC1AUIzEBfuozMj/EUfRmyAowu+IyjAU6QEbIzqK3iKoACvkRKwgfIwIyvARwQF+IjzNTZGVoC/vO3MyDwKh1LtMsY9HbAPRkvCU94GBWLB0bwchFbpmsARgH0wBRP84m1QgE9GjzHjNI1dlJUKZAV4xIGgkN0yxsRxXH0ojuPW5fDPoFlrqPjF7sgK8I70zoxxHIdhmGWZUipJkmrPA611kiTN5TgoxkNCCDo2wi/Sg0KSJFEU2RoF2+hsKw/sf4uiqC3HQZESIMrtcVh2rQXcJf3O2bV7extjlFJZlmmtgyCwNQ3V5a3Pguu0ftbX9EBKgEwcmcfj5dVHeo1CbY/neW4zgboNB+XPeZ5vWC6IwbkYYtEGAS9IDwolW4ugBjcx6AHWLTE2QEqAbK/uKskJB85yIygYY8IwDIJgeJVOMcCqZcbqSAlwAvUKcJwDwyPtV/80TattDTg6UgIcofUzpZ6+7NWo9as6hluDhv4C+5He7aLWabHU7L1YXc3L7iRHVu/MSEqAi+xxe/egPdFRF67x8uojukbBXviNMbWgYGdYsvMrlI8yPPIoSAlwGhMxwTUOBIUkSezESpatOTDGRFEUhqFdGEURDRPHwqkWzmHSRrjJw0oS5Wnlz5G9qp6lOgGuu3sM0/TgGS+vPqJrFObgNtMeIiXAA9QrwDUeZh/laaY7tGrs452FBzikPeXl1cfbGgX4g1Mq/FPWK6j2MZN316VtAnsiKEAwIgI8Zg9prZW9d1THEa71sy0LBTS5MTMjjqh6S/He71uAw5i3EeIRFCBStd8idQnwG1kBshEUIIzWjG7A4VSzAnEBwhAUIAmdEnBY1QOerABJvO3MyDwKjiEiAJXuja03hgB24W1QIBY4oxbpeONwcHdHThZKKcXwSOzJ26AABxARgFbVqgVF7QJ2Rh8F7KTW1sBJEKipfS7ouICdEBR20+xFIdPy5ayNa1giIhx3Z66AQi5ofjm1err2mIjj7ExMQ9MDNkRbAzCe1s+UelqoB+Xv6u4sZMzxjFURFLAJIgIwSSUE3LnT+svcUBTM8Yy1eXifK+XI/bucKKQaUM6e89Sr70DlylMnY+7/zuTNzpSAQi5orXI2KuHn3Fbq6DtzUU4UcixvaxSYR2FL9dNQa1NiUSileA+ABdSGRVRDeeNER5UDZvI2KBAL9tGMCLwRwEpabxLBWEoszeFKkizLsixTShljjDHVhwbUlvuwgohiaK1b6wkqT9m/kMd5O1hh8ApCirHwCm2VefrUFycJWyGkGNtsqXNcrVHIsiwMwyAIlFJJkkRRFMfx3oU6hpMjlLz7kADOaDRJvHR3CfdtxyiuZh+tdRAEtkYhjuMkSaobIiE2epLBuyoMqqu48H3Fk7eDFZZbQUgxNlqh52EZH+Ft/oqEFVzk6iZprdM0LVscmr/ufjS48dEamAPqS2pfRx7sviskrCCkGKwwcAUhxdhhhZYREwt8GVAnailOnCVayunpCi5ysumh7JpQW1hbsjYhfYnLYjTHIpZrzP0bd4/7em7QHX8XgEDNy9iQU8SpdVrOP3dqeZc5Swg56x6Nk0GhlU0PJX3qsJ6/guq6MG9ahlfF6Fr15EucXmH9DfFjBSHFYIWBKwgphogV+h+etkKtY8QSb4eEs+6w07JX/AkKd6sTXlaCdX7J9lTb8Tuk19Kx9hKAkUafRnyrfB/A4y6i/gSFqvqkp4cx9cN5rL0EYAVHP414nI2cvHukrTyotTVs3EEBAIAjcDIoKKWCIAjD0P5sIwJBAQCAxTk8kMMYk+e5/dndrQAAQDKHgwIAAFibq00PAABgAwQFAADQiaAAAAA6ERQAAEAnggIAAOhEUAAAAJ0ICgAAoBNBAQAAdCIoAACATgQFAADQiaAAAAA6ERQAAEAnggIAAOhEUAAAAJ0ICgAAoBNBAQAAdCIoAACATt/buwCr0FrvXQQAwBEVRbF3ERbmZ1BQMt6qMq9IKIw3tNbV/Vn7FaNU994ae5I361B4f5WnX1MlBoXLy8vr62ul1Pn5+U9+8pOu1R49evTNN9/cv3//008/3bB0I1RPwbUlgDRFUXCix2QcPB4T99ZeXl6+9dZb5+fnSqnr6+uHDx+25gB76T0/P7eR4vLy8sc//nHt0RoJW0pimK95PuIMNc3ae5J36lB4cy0v94O4TdJan5+ff/PNN0qpR48ePXnypFnC2vLqU8ol0rarhsQwGZefpWwfFBb/E5CDd9bycj9IHPXw8OFD+4OtS7i8vOxf31Y/1OiGxcs5R3FLZvHm8GlbAAzh5dURJVlBwWaCWr8E27hQZQPEo0eP/vrXvz569Oj6+vr+/fu1dYqGNQs+nceJAc6xx+HepQAgi8TOjDXVNoXSw4cPnzx58uTJE6XU+fl5sx/DkPOdqPTQ7PmohJUQNU5/i3K68AC2JKtGoVWztsD2Ubi8vCyKwg6RGFKj4Fwdg4vVDPbys16B5VzeHHpTdiHnnTogDk4sy4Gg0PTkyZOHDx/aYQ4//vGPy+GU/nE3McyitRqzpVSYL4idua+ldv6WbyKh0HuygoLtnVDrvdjaV/Ek4Z0ZxzpQYig3TfA22jMj11Qsa5HDics2FicrKCilzs/P33rrLfuzbVCw6SHLsjiO7XLbQeGvf/2r/bV1UiZXGhrG8jwx1DZH5NZ5cCLecRNIV/2W2j/sZyxI4inv/v37ZVNCWbwsy8IwTNPUGFNbpzaJghI84dIa5EzJUF5+Jl6Hus5rA+YA3uzKV/tDjoaG/mIvslE9L+LoTlvb3I/P0q8z9s/B8nKHiKtRsB0Ozs/PbYvDo0eP7HJjTFEUNiUopcIwtCtfXFxcX1+/ePGi9jq+1ig01eoY9i6OUtO+zfSsL2OjZLHdOIb8G2/+l1FZ58o195U01T1PpQKWIunzrJQaNjOjXe3rr79+4403lFJBEFxdXXnwPW8RO2579U+PK8aQ01nvl6RBf67nrwwrautf2W6HL3fS10r1l/jkCnVjZnjcYo/tdYGcul1L1VQ1n7j23j7yybaVlztE3CZprb/66quy20HtV+vx48cff/xx/5mouVDalq7EvaAw8Jx+qja1/tDkS8XICvN1d/gKF7whIWB0UGg+fWxQOMh332GHkCtBwcuL4kxe7hNZEy51zcxYW3J1dXVxcfHixYurqyul1MXFha1aqPLvrRqo2OkegAv3S2hdc9kqipPPnf8Nz9mLXzE7K5zYdmf3zFxtR9dyr91ylO51QoBPxPVRaGrOzHh1dXV1dfXmm2+WPzx+/Li2TnN4pGcDJn2wYD+G6jpLvbO11uuehm1/27yxltuDpOfqvtSforMCZpJVo9CqOeuiVX66bEvEb3/729ZHD0jCd4jTZVjjzLXO2bD/u/UCX779cuS9MXbbX65vj1vXTlm7n2SwGQeCQtPFxUXPr1YzQXNMSzHzct5sg1jz25JnlxE20ncAABTRSURBVD0Jm+NltLKH4KxNu3sYF/ZK3LP+mBOahC8PcJespoeBMzO2JoOa4wyPbCW3snGRUlVnbyQlrKO4vfjhJHuc2EPF7rchu26BXiCVg58cgPWIO7bsTEq2VNWfsyyrTs6otf7www9tc4PW+uLiIs/z8kX4zKhtd8LQcYOzL+r959bFr+sDX9ChPDGqqNO2a4NRFXKUFQmtD80cgzp0Lw2eXmnBcwLn2C5e7hmJmzRkZsYXL168+eab9qFaSlDHHh5Ztdd8hS3Ll7rVzYZBYYNr6vbW3ij/olWXnohQXadrhcV3VP+Q1FerLXRO8PJyuAgv94yspgel1PPnz6+vr8/Ozs7OzpRS77//vl1em5nxjTfeKBsU/ve//zVf5+BND7JskhJwEjtwKdW2hh5dzTd7VX0B04jLPlrrs7Ozb7/9Vin1/vvvf/bZZycnVjr4vR767Tvf+4neWGP/yuLzCS73ak6cqderIRi7vhO7q9WQioTmU4pTS4Y/9/RqmzRAePm9eRFe7hlxNQqqUovw+9//Xin1/Pnz1tUePnyoum9CTY3CZjo/GDJ7Uw7j7pUM6xlYkVBT69446tCa0qV0/c+dl9dC9JAVFGwm+NnPflZdaGsXar766qvf/e53HKwnLTb8YeytdFxOCVjb0YZUFGtuckvy6P2Qyh0SBakcmEehNSj89Kc//ctf/tLzrCGfBHLGCaPmV26z8Yj5HQfoy58bYFrxRm2X8D0w0/ytKyr/3VrzEzr4TmyAA0Hh3r17tSX3798/Pz9/6623ep5FCCiNm2tl768autrYeqyvnUCLNTLo0Ndsu8EK7Q4H5EBQaLKDJ+3UzuXPH3zwQTU6MDPjFONTwoKnMN2sQFWFIi70au60GvZeyct9tXotjpvTS2NZsoKC7Z3w/PnzajeFZo3CBx98UP5sg0KzgoFYULXPnReGOXn6nhYX5lwVJFehn9yu7vX12OcOf6KES+zUrRv+Cm7sh1WMuncrvCOuEunevXs3NzdRFCmlvv7665ubm9aZGZVSjx8/vrq6spMoMDzypJXGJswZsjj+zH56sKX9+jPyZVtfpphw0q9t7PxL13Jmhp81hj1KzmNNy5R2cpLo+vNzijVxuOyxT6Qnedk0I26Tsiz79a9/fXNzY3+NosiGg9rMjDYKXFxc2KAQhmH17pFevlXz1XfLErUIE+brnXH5fPlKtbOtrr38QkFhkydudrHcICiI3fZF7F7aVwUoj/9t0l/9KZxae/l59WnON7AvpVQQBPZnW6/QXOfDDz+sLq/9al9kzTK66tVuuR2vNf+fGrPOJn9t+GprvMKoJ6rb/84v8Db7pP9FJvyJDTZczg5cowBbb9Krp6CDl1cfWfMoWGX7QlmXUFvh6uqqegPJrttM16xTWJcUdvz0crtizMC5QlI9fL+NvzvaE6/849OVcq5h9+oE1bXz53yyxr6dd/YCZ9QjkdWZ0WaC8oYO5cLaktotoK6urpovVfhX+TOf1sXm5zx9rKvL8B1cW63cS2u8OQu+56uWExOVWWHT3pQMiDgMWUGhVbNGoerx48cff/yxbX2oYsKlO/aI/+5UIQhhd5f8y7Ar5fRPLYa2Tcn4sj1i27hwnBPpUTkQFGrVCaXyTtNff/31G2+8UXv0QCGgX3M+iSVHIgwccbd9RYZwPXujOLWCELVyyi/wZO5t2qhPse7dwEGZg6oF30nsozDE48eP33zzzQ8//LAoimZKgFKq70YMzXVVMb6x070TqDtstJLfZONKOTFReWY4/Wmn14K/xA3k0FqXYyCbv1q2LqG1IqF8VnOhtC1d0YBP7BLzDWw8aH6bP7dBOSXvt41JLltJWiH10lV043rVnK5jOM6Zto2XwyPFNT0EQRCGod3RNh/Y/1YnXPr444/tyi9evCifWAsN/r1Vg5yKCNV5+tYuy1HRzgK3jDtiT3eDoNeCdyRmH2NMOa6hLF51wqWTFQaHqFEYU9HX0aaw5Xwtx6lR6H+R7WclEp5aKN5Yi9coqAGv1rlCZ2Lw7Hw7mJc1Cq5u0pdffvn222+fnZ0ppW5ubt59993PPvusfNS3t2pS49+ADgcbT+zmRC4hKGyM4k2wRqn6X/P0Xxwx1MKnk3ODb1cfpZS7QUFrfXZ29u233yql3nvvvT/84Q8+1CjM7g004e4JG15QCQrTXnz7iXq3RPGEmBsUStMHZ8o/RQ/gZVBwddSDUurdd9+1P9i6hC+//LL6aHMSyh2K2MMOSaj961r3VcfjE/+23AIfLXVVKOgC4oXjpIR+4/bD9HNR61mRwRQCOBkUbCb4+c9/Xl1Y3kdKnLZDf/i1f+XLv6+XNIHbtctVh0sdBlr4IzPqFHfiREeA2Ju4UQ+T1YLC8Js7lJUN5VOmLbk7Y1p1SdFY0lznIEv0jL1RLlyvhGrAOs4t0cLK42IJj7Ok9b1QjSWr/HU99lnzztjrLfGPP0HBdmwsTWhraD5l2pL6d7g7R4++/d+pZ223pLh7dhjyrHL9LcvcXLj432pNJ0u98vZ7rPnOyjnqyiViS9hVGSOnhIsvaX0vqgt3L+Hd3xc6Yy+7xMu44E9QqGm+W7t1U2j7u8UCgxuxrGVr6bviF2rYUXIs+V4seb8JaT3MjsfJPgq2d0Kt92KzRkF0Z8aiaPnXta7SA/9NLo28Fn2fzDz58u7AJdPPRa1nRWmn7kOSOJDj8vLy+vpaKXV+fv6Tn/ykdZ179+7d3NwEQRAEwYsXL25ubmrDIwVu13RrzaNwd/UR60++8s25ZE547jYbdfIF57/yLjt8MzILKbNUa6tt9YidwDwKlm9XH6WUwKBweXn51ltvnZ+fK6Wur68fPnz46aefNlez8yiUHRg/+uijJEmqjzafIm1L51q+8WLgSWHmCXTLa56EoLDIyxIUNiawSJup9Wc8sR+YmbHGy6DQUkW/L6XU+fm5/fnhw4etJfzoo4+qy2u/2hdZtZByKdX/78SDJ549arXFnz7ticOfNXO7ul5zkZdde9v3/SewnAKLtMu2d+6HE69yYF5efST2UbD5QCll6xIuLy9rK+R5HgRB+Wv155JuWK28kpxq0uvt0FCcagvX/n7TWm+7vNxdOK7TXRD8+z59eLKCgs0EtX4Jtr9Clb2TZPlreQepqmYmWr64Yg3rAXSqF6Ru/Cv2u+xNvpCfTD+rOtJR5w9f0/BA5Ufmzn4Y1GOavoeecmB45DfffNPzaBRFv/nNb2zrQ9WQKgTP04PdumFVKeXnX/fNJQCsoTj8tVkuuihC7RIUavUBpTiOW9e/f/9+6/J//vOfr7/+ulLqH//4xw9/+MPao56HgOGKYlS3x8pe65kbqvaU9tfveCJXhQnYaRhh5hwGWil71hj6OpxsfSer6WG4KIpef/31jz76qCiKZkrAHfPqAydP3jBvaocF7dv6ALe4HcgW+dCNi6WkhAMQN5BDa/3VV1+V3RRqv1q2LqG1IqF8VnOhtC3dx06dOu/WLqw90HHs0524Nmy/07Ykp6hySjLCPomcM2obL4dHituk+/fvX19f21JVf7YNFrZ5whiT5/k//vGP6hOrocHLt2ox+w0AuY0Lu1zzel7BiWsDQWEbckpy2p41dpxgO3h59ZG4STYf2J/L4mVZFoZhmqbGmJMVBtQonLZnXFAzr3m18+Ps6aScuDaMLaQTG1WSU9qdSyKjta4XJ9JeXgYFcX0U7PzN5+fndnLGR48e2eXGmKIojDFK1edWuri4aL4xhx4eOcSuA5kmt6S2PnHY3S5G91QYfouNyf/G7wBXObIfdkgJUndFB06khyQu+2itz8/P7ZDIR48ePXnypKeEtubg4uKiNpWCl5luLXtULdydJHZQfcD8t7PrIiDny2yTru4nwTUKS3Szn2uJm6yeKIgDF/L1cEYdxsurj7gaBTVgZkbLzp1wcXHR+uhBZ2acYO85UoaMnljvr0tOCerO102lRtZVbFAjstT34EXGpSyyFY7UfGzOuysfRpEVFAbOzKiU+vvf//7b3/62J7jR9DDO3nFBtTUrLHt2ZpQkMJqAMwN25+rMjG+88caLFy96nsXMjFOMnJ1plSJwNccehNct7YMzJJRSjs7MGATBxcXFj370o56/QgiYqLnfTkaH1l29d+BoxQSHOCjOh5jBgRqFpqurK3V708jy58ePH1ejQ7NGgegw0bT9JqBy4iRfQ4OL20WGWxLnOixqh6BgjLGjHJts74TLy8tqNwU7TrLq8ePH5c9lUKhVMBAL9jfmrlQAWWEZnPqwNHEDOYbMzFhl6xUYHilaW1YYfklY/OJRvqBzl6WBBXZuu0p7ldzdPXYHJz0BvLz6iGt6+Oabb+7fv182HFT3eJIk1dqIL7/88ubmput1aHoQxIVmCEhApcJ0nN+wGlnDI5VSz58/v76+Pjs7Ozs7U0q9//77dnl1Zkal1Hvvvff222/f3Nz897//vbq6+vLLL2uvw/BIWSQNsmKcpGS8O6NJ+nDBS+IqSbTWZ2dn3377rVLq/fff/+yzz1pLqLX+05/+9POf/1wpde/evZubG+714Ibbt2Zu08O8DhD2ac4dEN43PZQ23gSH9xinNWG8bHoQt0la6z//+c8/+9nPWn+13nvvvT/84Q/9UztL2y68orWaGRRqb+6kuODiteE4QUFtuBWu7i5OcSJ5efWR1Ufh+fPnSqlaLPj2229rS25ubs7Ozso+CmdnZ7ZqAW6YORqi+SGc9IK+fZRvuXrZg3dXF3hDXB+FJtsMUXVzc3Nzc2P7KNgf3nvvvdo6zXs9cPcHUYqimHT7yO6TKS21fjlKZwV73HLoQjAHZma8d+9e6/Kyese2RHz22Wetj8IfQ97TBSdv2OAQGl8LcqgKA5+3lxMU3CGr6WEgOyCi61eL4ZFuGD5yctTbN2dA5pbHycm/5Xq91+ydWZSNvuvsiq2DCGchOEjWzIy2L8Lz58+rnRKaNQpnZ2c9MyhYxAJnrHQNGHXTCrFHS7VgroSGlXamhFAl9jgB1iSuf6Yd6xhFkVLq66+/Lsc91mZm1Fq/++67//73v40xSZKUIyrLR5uvLG1Lcad7cPdZXvPeVZz8el3/iuzFrpvYk3xAdLizu7zYV9gXox628Mknn/z6179OksT+ahODVZ2Z8U9/+tPbb7+tlMrz/Ozs7JNPPqm9jvy3ypXjaaNyzqtXONzO7HkRrWde8ATuzKIoarO1DirkkK2YvbtOvby4ndnkRCGVO+X0j7j9rrUOgsD2dozjOEmSZglry6tPKZdI264mJwqp1ixnyys3s0JRDCnAcXbmyVfY4E/sa9nuR2tvrPCdaTlRSOVIOZ0o5FgSh0eW7Qv2h9YhElX2vlA1DIaUr/o1sVzU9ytO8fIkVVPOyF7Ozj75M36E3QXMJyso2ExQ6+rYDAo2QMRxbHst5Hne7B3JvR5cVb5TvGVtWtLV4fEBB1Ylro9CU2uNQhRFSZLYrgxBEDTnYBhyMuXkIhTvC8YrhrVSARjLgQmXmrUFto9CmqbGmCzLwjC0P1TX4XzhB079ALAvWU0PAyVJEkWRDRDGmDRN8zwf9Qon6xs2WOGkIa/gRDkPUsiB66xdDD925rRCVttltnlDPd6Z25fBiXIes+FP1oRLdnmWZdUVulbux8yMAADMJ65GIQiCMAztz2WdgVKqOtuS7aBQtjWU61fRmREAgPnEBYUsy4IgsIOd8jyvXuDLcBDHsc0TdrUgCJo5gOGRTqAP/wRdO43+HIojCliB22eWWiNFaf6kNE6sIKQYM1ewp/WuFexzdy/kkBW2LEbrag7tq1VXGLIf2FfSVhBSjG221DkODI/s0dN9gT4KwDFRqQAsy8Pso47aMRUAsDv/rqp+BgUAALAIcZ0ZAQCAHAQFAADQye3OjP1qc0J3TREtRHMWalGyW3EcT5sCawNlIY0xwt9uK45jUeUsp1fvmRVNDml7r8aJo9GJz3WV5POkW1eccZoTE/khiiKHttTeKTtN070L0s7uzCAIbDntxBXS1Aop/B0viiJNU1Fvui1PuQOjKNq7RH2k7b0aJ45GJz7XVZLPk25dccbyamOq7NG/dykGsac8sR+AoiiqZbOfh12L0652bZN8qUvTtLx+yHnTq5cKse9yIXXv1ThxNDrxuS4JP086dMWZwOc+Ck7UpCmlwjBsplE5yrpo+6utT5NZ+1et6wuCQGYhLWOMwDe93IGS32Ulde/VCD8aHfpcW8LPk8qdK84UeyeVtSilyq8dQRDIDKFFUQRBYL9qKKlJuar8Mrd3QU5TIr/D1ch50+3XteoS+TtQzt47SfjOlP+5ln+edOWKM43PnRmVUvb0F8dxGIaFvBkj4jjO81xyiq8yxtjbeZd1gDJlWWbvE+ZVZ6I9uHJkSib/aJT/uXboPCn8ijOZw0Gh7KFdYz+Q1TcpyzKt9S59pHsKmWVZkiRCDqb+nVmuo24/A2mabl/PNqSQ9qy3e03vkKLK53NV6iaEHI39dv9c9xN1nuwh5IqzEp/7KFQJ/KzaY8jcUkqFYSjtU2pVd53ktkx7x9E0TWUWD4ci/2h04nPt0HmySuAVZw6HaxR6hnrbYcHV9ynP8106wvQUslnCIAj2SqD95RSS6PsH99sbjgv5cDoxD0GVLW3tdqxubYIooo7GVnI+1/1EnSe7yLnirGW/7hHrUo6M9SopqZ10irtdsezOlFZU2zQYRVF6197lOkHUnqx2ZxPetc0StfeqXDka5X+um8QW0rkrzihebUxVrWOOzGOrSnIh5e/M1vAuf1iztJ1Zdtt24jQnbe+VXDka5X+um8SW08WdORx3jwQAAJ2O0pkRAABMQFAAAACdCAoAAKATQQEAAHQiKAAAgE4EBQAA0ImgAAAAOhEUAAAYx975qfzVuVnbR2HCJQAARivvkGJvJu7xxZQaBQAARsuyLM9ze4dur24B1UCNAgAAU9i6BOF3Cp2PGgUAANCJoAAAwBS20cE2QOxdlhXR9AAAwGjH6cxIUAAAYJw4jpMkKS+gZWjYsUjrISgAAIBO9FEAAACdCAoAAKATQQEAAHQiKAAAgE4EBQAA0ImgAAAAOhEUAABAJ4ICAADoRFAAAACdCAoAAKATQQEAAHT6/8ogUYaLu5lgAAAAAElFTkSuQmCC\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"import ROOT\n",
"\n",
"ROOT.gStyle.SetOptStat(0)\n",
"\n",
"c1 = ROOT.TCanvas(\"c1\", \"fit residual simple\")\n",
"h1 = ROOT.TH1D(\"h1\", \"h1\", 50, -5, 5)\n",
"h1.FillRandom(\"gaus\", 2000)\n",
"h1.Fit(\"gaus\")\n",
"h1.GetXaxis().SetTitle(\"x\")\n",
"h1.GetYaxis().SetTitle(\"y\")\n",
"\n",
"rp1 = ROOT.TRatioPlot(h1)\n",
"rp1.SetConfidenceIntervalColors(ROOT.kBlue, ROOT.kRed)\n",
"rp1.Draw()\n",
"c1.Update()"
]
},
{
"cell_type": "markdown",
"id": "f47ad160",
"metadata": {},
"source": [
"Draw all canvases "
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "865d301e",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:44.687504Z",
"iopub.status.busy": "2024-03-19T19:11:44.687176Z",
"iopub.status.idle": "2024-03-19T19:11:44.845506Z",
"shell.execute_reply": "2024-03-19T19:11:44.830300Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3dzY7lRv3/8aq/uBBugkg9PWqXCckKgoTIBcAiI7FMoswIZWIboiiDCDsizWzy2xNGIo1YZADbrXR6EbgcVgj5v6hpj9tPx8/+Vvn90ijp9vE5XfbxsT+nnqyLolAAAABt/t/eBQAAAHIRFAAAQCeCAgAA6ERQAAAAnQgKAACgE0EBAAB0IigAh5Zl2d5FACAaQQE4tDAMex7NsswYs1VZAEhEUACO6P/+7/9ee+21/nWyLCtjxBdffLF+oQBIRFAAjujzzz//7rvvyl+zLIvjuNoMYYyppoRf/vKXG5cQgBCaKZyBA3rttde+++67oii01kEQ2IV5nkdRFMdxuZptd8jzXCnFuQI4JmoUgKMzxmRZlmVZFEVJkpTLv/jii//85z82JQA4LIICcGhffPHF3/72t9dee+211177/ve/X33o888//9e//rVTuQBI8b29CwBgN7YBovz1F7/4RbXT4nfffWebHv7zn/+QGIDDokYBOK7vvvvuiy++KIqiKIpf/epXSqnPP/9870IBkIWgABzXr371q1/84hf25z/+8Y/7FgaATAQF4LgIBwBOoo8CcGj9gx7tzAo/+MEPNioNAHmoUQAAAJ0ICgAAoBNBAQAAdCIoAACATgQFAADQiZtCAQCATtQoAACATgQFAADQiaAAAAA6ERQAAEAnggIAAOhEUAAAAJ0ICgAAoJOsoBDHcetCY0zzoSzL4jhufQoAAFiEoAmXsiwLwzBNU2NMuVBrrZQKgiDPc1W5JW4cx0mSBEGglMrzvPYsAACwjEKANE3tJV8plaZpudwuLH9VSkVRVP5crhkEQRAE1RfccX8CAI5s3evlHr43YS9orYMgsC0CS+1ZY4wxJkmS2vIyQNifsyxTty0U5V+P4zgMw9oTC+ICAGBbthbcM1P6KKRpqpQKw1BrbYyxF+85bBeEZm+DOI7zPI/j2HZHsD+3Pl0pNb8YAACgZkpQsOGgKIooitRtYlijU6ExJgiCJEnCMEySJIqiMhNUaxpa6QEWLzAAAJ6ZNerBfte3cSFJElvBsEy5bl/fdlQsiiJN0yRJhseRIe0uCxYVAAAvTQ8Kto+C1tqOPkjT1DZJLPhNPUmScjiDMSaKItuJwRhjB0FYttGBUQ8AACxuYtODzQfqduhBlmW2N6K9Zq/dXaCWCeidAADASiYGhWo+qD264JQGQRCUwxmyLCsnTijrGOxDtvvCIn8RAABUCZpwSSmlta7ljGorQzk8Ut3OztRcXr6OqO0CAByBl1cfNzapteqih5dvFQBAOC+vPh5ukuroUOnllgIA5PAyKEyZmdEJ/r1VAABsT9bdIwEAgCje1igAaKX1s5PrFMU7G5QEgBMICsDh9OeAIUkCwHHQ9AAAADoRFAAAQCeCAgAA6ORtH4XmVAoMmAQAYCxvgwKxAACA+Wh6AAAAnQgKAACgE0EBAAB0IigAAIBOBAUAANCJoAAAADp5OzySeRQAAJjP26BALAAAYD5ZTQ9xHDcXZlkWx3Hzoa7lAABgKYKCQpZlSZJkWVZdaIwJw9A+VG1NiOPYLs+yTGtdexYAAFiEllBFb+sG8jxXSqVpaowpl4dhWJZQax0Egc0EWutyTfvfalbQWsR2AQJp/awo3ulf4eSL9L8CcFheXn2k9FEwxhhjkiSpLozjOAiC8tdy79vmhjJP2NqFbcoJeO9kCBiSJAB4Q0RQsClBKVULCnmeR1Fk2xdURw8GValRKKMDAABYhKA+Cq2SJKn2UbCJIcuyak1DKz3AFhsAAIDLRNQo9CtbHGzHxoHNP/61EgEr6srNfI6Aw5NeoxBFUflz2bJgjLE9Hy1bzUC7AzCF1p0p4eSjAA5AdFDoal+oZQLGRgJT1EJAUdT/da0J4EhEB4U4jsuZFWw3BVvBUB0VqZQqlwMYqhkRmppxAcDxiO6jYIyJoqgc+hgEQTnwIU3TMAxth8TqcgCnlZf8IV0Q7Dr2KS+f+HSdYgGQyI2pIcYOffRyygtgAbWKhM2eCxyDl1cfDzdJefpWAXMtcqUf8CLM7YjD8vLqI7rpYQ5uMw3ccfuJ0OrprIt0UbzKClp3ZYX5s0QDEEJ0Z8Y5ioa9SwTsZ1SnhJMvpp6+eh16OAK+87ZGAcBLi6aE25d8ptTTQj1QVA8AviMoAF5pXrZtOtDqqVroiv6qWUE/UEoV6gF9GwGPERQA39zpH3BbnbBK58Gyv0J3ZwUArvO2jwKANRod6uisAPiOoAB4aoOUsM3rA9gVQQHw2pZXcSoVAB9520eBeRRwaBtfs+msAPjL26BALMBxbdboUFWdiAmAR2h6ADy1V1YmLgB+ISgAXrGTIO30txkBAXiIoAD4aK/qBJr8AO8QFACPyPkqL6ckAOYhKAC+2KUPYxMNEIBfvB31wPBIHJNWT/c/0BkBAXjE26BALMCxyLwwM60C4D6aHgCPyLkqyykJgHlkBYU4jnseyrKsuiTLsjiOe54CHIXM6gRLctkADCAoKGRZliRJLQ10PRTHcRiGWZZlWaa1bn0WcCzSvsRLKw+ASUQEhSzLjDFhGHat0HwoSZI0TW1QCIKAegUcl/yv7PJLCKCbiKCglDLGRFHU9VAQBEEQlEtsLDDGlL/meb56EQHJZH59l1kqAGOICArGmK7eBjYE9Lcs2MRA6wOOyJUv666UE0CD9OGRtomhttA2N/Q/sTmPQhNDKOEJyUcycyoAjhMdFLTWURSVTQyjEALgP7cuwMypALhJRNNDK9uUYPs5GmPyPE+SxIYG+2ttzWl5AnCe/Kuv/BIC6Ca3RqHWvTHP8yAIyqCQJEn5EL0TcERuVSdYVCoADpIbFNTd+Zds1UJ1vIMxxkaEJEm6RkwAnnPluktPBcBZcpse+qVpmue51lprzTwKOBwuugC2IqtGoacHYq19wRhDd0UcnYsfAVofANe4WqMAHJej1QnkA8BN3gYF3bB3iYBFuXvd5cMIOEVW08OCaJiAn5y+ytKlEXCQtzUKgM8cz8GFerB3EQAMRVAA3OHB13HHIw5wQAQFwDV+XGs9CD3AMRAUAEd4c2X1I+gAh0FQAJzi01XWm+gDeI2gALjAr2uqVk/3LgKAobwdHtmcOIEBk3Cef8cwEzUC4nkbFIgF8Idf1QkvMacC4AiaHgBH+Jp9iQuAbAQFADvxNfoAfiEoALId4Qv3EbYRcBZBAXCBr1++fd0uwCMEBUCw43zVPs6WAq7xdtQDwyPhIq2fVX+1h6xWT1VleVG8s22hVsbwB0A2b4MCsQCOepUDbi+f1WRQSxIAsDZvgwLgiUbk9TYrMPkSIJKsoBDHcRzHzYVZlhljag9lWZZlmV1hm+IBu/Ot3cGi9QEQTFBnxizLkiSx1/6S1jpJEmNMkiRa6/LROI7DMLRZoboc8MQxL5zH3GpANhFBwVYYhGFYW26MUUoVRRHHcVEUQRCU6yRJkqapDQpBEFCpAD8dpyr+OFsKuEZEUFBKGWOiKKotzPO8utDmBnXb1lD9Nc/z9csIbOXIX6yPvO2ASCL6KBhj7FU/SZLq8trIha72BftcWy2xSvmAXRztSzY9FQCRRASFIYwxeZ6naaqUss0N/es351FoYgglJDrGxbJ17EZRecjPbpuAgxwICnEc25qGNE2H1xkQAuA2rw/gzhCgHyilCvVAq6ebFghANyl9FLrY8Q5RFBVFUaYEW7tQrmObJGh3AJzndTwCHCW6RsH2UmzWDdj0UP7K2Eh4o1AP9i4CANyhRVXRa62r7Qta6+bQR/uofchGBK11FEXV1bSWtV3AULaDwpGP3rKLxpF3Apzl5dVH1iY1g0JzHVvgLMvKORXKxFB9HVHbBQxFUFDsBDjMy6uPh5ukPH2r4D++TFsEBTjLy6uP9M6MwOF4d5aZ6BjDRAH5RHdmnKPZbOFfyoNXuC6WmHkJkMTboEAsgIu0esqB+wo3ngYEoOkBgDhMuATIQVAABKCmHYBUBAVADKrZm4hQwN4ICgBEIjYBMhAUgL3xpRmAYN6OemB4JBzD8dmFsQ/ArrwNCsQCwHlMqAAIQNMDsCsuhEOwl4D9EBQAAagA68KeAfbmbdMDIJDWz6q/FupB63IAkIOgAGyqKN559Yt+YBfxrfk0ujQCO6HpAYBs5ANgVwQFYCd00BuLPQbswdumB+ZRgBs4LIdgnCSwH2+DArEAonHZA+AImh6A/RBnxyJgAZuTFRTiOG5daIxpPpRlWRzHrU8B4BtCFbATQUEhy7IkSbIsqy7UWidJopRKkqTa7SCO4zAMsyzLskxrXXsWIBpfi+dg7wHbEhEUsiwzxoRhWFtuawuKosiyzPY5KOsPkiRJ09QGhSAIqFeAe/iKPBZ7DNiDiKCglDLGRFFUW5gkSRAE5a9BENiaAxsLjDF2eRzHeZ5vUkxgNr4QA3CKiKBguyC01gqUacD+3BoI7Dq0PsAlfDmeg7AFbMjJ4ZG2uaF/neY8Ck0MoQQcw4QKwOacDApDEAIgUHkXKJzUdaOs4vbRO3fNALAaJ4OCMcYOhbBso0O1kQIQjRR7Sl8I0A/Uy8hFUAC2IKKPQpey96JV9m2sZQJ6JwAHQswCtiU6KNjhDDYHNMc7lHEhSZLmiAlAHBrXAThIdNODHTNZzq8QRVEZDtI0DcPQ9lhkHgW4hC/ES9GanQlsQHvZ6U9rP7cLbrM1ChyZi2BnQiQvrz6iaxTm4DbTkIV2BwBu8jYoEAsgkFZPOS6XROsDsD7RnRkBoJVWT/cuAnAUBAVgfbQ7AHAWQQHYCpXkayCEASsjKABwE8EL2ARBAVgZX3kBuIygAGyCr7/rIYoBa/J2eCTzKAD+467TwPq8DQrEAojAZWwbTKgArIamB2B9XMPWw74FVkZQAFZDdQIA9xEUgJXxlXcbxDJgHQQFAI4jigFrIigA6+AL7vbY58AKCArAmviyuw32M7Aab4dHMo8C9sRXWwC+8DYoEAuwPw7C7TGhArA0B4JCdiuOY2NM9aE4jrMsM8bEcbxP4YAKrZ/ZH4q7v2ILzNIIrEML/+Ydx3GSJEEQKKXyPA+CIMsy+5BtXAiCIM9zdbcKQWvp2wUvaf2sKN55dbniIFzNy13dWPryB/Y8duLl1Ud6Z8YkSaIosjUKaZrmeW6Dgq1CKIoiyzL7rlCpAFm8O1k4gH0OrEB6UFBKlc0N1XaHsprBqtY0AHui9huAX6QHhSAIwjC0NQo2KLTmBmOMbYAAROCr7b6Ia8BypHdmzLJMax2Gof01TdOBT2wOj2zyryUJODq6NAJLkx4UtNZls0Icx2EYpmlaG/vQihCA7RXqwd5FOJCuQSXlkJOW3o4AxhMdFGw+KDsf2BEQZRsEIBQhdX19IUA/UC9DG0EBWID0Pgo1Ze1CrfdirW8jgOMiqAGLEh0Uar0X4zjO89wOg7Q/V6scGB6JndE0DsBHopselFJpmoZhWPZMjKKoTA9RFJWdHMvlwM74OisH0zkDS/BwDinl6dxYEI05AaWx7whvB7bl5dVHdNMD4Batnu5dBNxFexAwm7dBQTfsXSIA2yG0AUuR3kdhMv8qfyAXMRSAv7wNCsCyeu4Z/XKGH77CykSXRmAeggIwVOckP3aGH+YBlIbpnIEleNtHAdgIlyIAXiMoAEugclsywhwwA0EBgL8IcMBsBAVgBr6quoJ3CpjK286MzYkTGDCJtXBoSUaXRmAeb4MCsQCr4/LjFsZJApPQ9ADMw7VHPt4jYAaCAjAJ1QkAjoGgAMzAV1W3EO+A8QgKAA6ASAdMRVAAxuOLqbt474CRvB31wPBIrI4jyi2MkwQm8TYoEAuwFi42rmOcJDAGTQ/AJFxpXMS7BoznRlCI49gYE8fxwOXAWqhOAHAwWn4Vve1tEARBnudBEGRZ1lyu7rY1aO3AdsEtWj8rindeBgWOLvFevl9tD7z8gTcRK/Dy6iO9RsEYEwRBURRZlqVpmue5DQq2FsEut+8K9QoATvPuJA6sTXpnxjzP0zS1PxtjyqSWJEkQBOVq1ZoGYC20O/iELo3AMKKDgr32V3shVKsNjDHVn5Mk2bBoODCuLq5jnCQwhuigYGmtbeVBnudJkgxs/mnOo9DkX0sS1lOoB3sXAQB2IL2PglIqiqIsy8b2RSgGWLfc8BKHjU+oVwAGcCAoVJMBfRGwD64oniHwAYOJDgq2F0I1GeR5bhfWEkOtbyOwCq4uAI5Heh+FIAjCMKw2OtigEMdxGIZZlhljqgMmAcDS+lnPo8XtSuQ/oJ8DU0MYY+yUSkqpNE3LwQ5xHJcjHaIoqgYFL6e8wG6YosdXTJ+FpXl59fFwk5SnbxV2o7VSSqun7TP9wV1EQCzNy6uP6D4KwP7oxugvrZ7uXQTAAd4GBd2wd4ngMu++IgDAQNI7M07mX+UPdkC+PAi6NALdvA0KwCitPeTtpUOrp6q3/zwcxnTOwCkEBeClel/F2+sHfRgPgUoFoIO3fRSAZXDx8B5vMdCLoAC0oTr6gHjTgTYEBaAb3zUPgjca6EZQABr4ZnlYvPVAg7edGZsTJzBgEuNwwBwKwx+ADt4GBWIBJuJqcSTVYbHlbaKqMzYy5gXwNigAsxA0D6ARAt6xMbFc3n//SeAg6KMAVFCdAMVhANxBUAAaqE44LN56oIGgAAANVCoAtwgKwC2uDVBUKgB13gYFbjONibhOwOKkASilPB71wPBIjFKoB3sXAWIwpwJQob28oGrt53ZhRfbCwGEDqwwKHBIYw8urj0tND8aY2pI4jo0xcRzvUBr4hK+PqPHuXA9M5kz2McbkeZ6maRkXbLeDIAjyPFd32xq8zHRYEdUJaKJSAeN5efVxo0YhyzKbBkq2FqEoiizL7LtCvQImojoBrbw73QPTuBEUwjCMoqi6JEmSIAjKX4MgyLJs62LBI9Xp/YE7iJI4NgeCgjEmiqJmhUG1y4JtmNiyVPAE1wB0Iz4CSv7wyDiO8zyfUFswZOIE/1qSMFFRKG7/gx53bylZxe0l4T3RQSHLsiRJpl3OCQE4jeoE9CqKd5q3lKzi9pI4AtFBwTY3VJsYwjCkOwIWRqbEEFpzqOCYpAeFaibI8zwIApseanGh1rcROI3qBAxUTtRIVsAhuTTiU2tdzqOQZVkYhvbX6s/lmg5tF/Zxd+4ErZ/R2Iw+bZNtcNigxsurj+gahR52KEQYhvbXKIqa8zYCnahOwDRUKuB4PMw+ytNMh8W0zbjHV0Oc1jhyOGxQ4+XVx4F5FIBVePdhxuo4ZnBIrjY9nNScR8G/lIcpaHTAfDRA4Ei8DQrEAvTh8MA0jIDA8dD0gCOhOgHzkQ9wMN7WKABVWj8r1IOXP6unTNiMBWituBkEDoCggIOhUQrzlQ0QwAF4OJBDeTpABbO0DYkEpuOIQhsvrz70UcABcE7H4jiWcBjeNj0wPBJ1HABYA8Mf4DtvgwKxAC/RloyV0FMBx+BtUACqtHpKcsRatNbdwx+Y4xmuIyjAa3zhw6puKxUK9aC1AUIzEBfuozMj/EUfRmyAowu+IyjAU6QEbIzqK3iKoACvkRKwgfIwIyvARwQF+IjzNTZGVoC/vO3MyDwKh1LtMsY9HbAPRkvCU94GBWLB0bwchFbpmsARgH0wBRP84m1QgE9GjzHjNI1dlJUKZAV4xIGgkN0yxsRxXH0ojuPW5fDPoFlrqPjF7sgK8I70zoxxHIdhmGWZUipJkmrPA611kiTN5TgoxkNCCDo2wi/Sg0KSJFEU2RoF2+hsKw/sf4uiqC3HQZESIMrtcVh2rQXcJf3O2bV7extjlFJZlmmtgyCwNQ3V5a3Pguu0ftbX9EBKgEwcmcfj5dVHeo1CbY/neW4zgboNB+XPeZ5vWC6IwbkYYtEGAS9IDwolW4ugBjcx6AHWLTE2QEqAbK/uKskJB85yIygYY8IwDIJgeJVOMcCqZcbqSAlwAvUKcJwDwyPtV/80TattDTg6UgIcofUzpZ6+7NWo9as6hluDhv4C+5He7aLWabHU7L1YXc3L7iRHVu/MSEqAi+xxe/egPdFRF67x8uojukbBXviNMbWgYGdYsvMrlI8yPPIoSAlwGhMxwTUOBIUkSezESpatOTDGRFEUhqFdGEURDRPHwqkWzmHSRrjJw0oS5Wnlz5G9qp6lOgGuu3sM0/TgGS+vPqJrFObgNtMeIiXAA9QrwDUeZh/laaY7tGrs452FBzikPeXl1cfbGgX4g1Mq/FPWK6j2MZN316VtAnsiKEAwIgI8Zg9prZW9d1THEa71sy0LBTS5MTMjjqh6S/He71uAw5i3EeIRFCBStd8idQnwG1kBshEUIIzWjG7A4VSzAnEBwhAUIAmdEnBY1QOerABJvO3MyDwKjiEiAJXuja03hgB24W1QIBY4oxbpeONwcHdHThZKKcXwSOzJ26AABxARgFbVqgVF7QJ2Rh8F7KTW1sBJEKipfS7ouICdEBR20+xFIdPy5ayNa1giIhx3Z66AQi5ofjm1err2mIjj7ExMQ9MDNkRbAzCe1s+UelqoB+Xv6u4sZMzxjFURFLAJIgIwSSUE3LnT+svcUBTM8Yy1eXifK+XI/bucKKQaUM6e89Sr70DlylMnY+7/zuTNzpSAQi5orXI2KuHn3Fbq6DtzUU4UcixvaxSYR2FL9dNQa1NiUSileA+ABdSGRVRDeeNER5UDZvI2KBAL9tGMCLwRwEpabxLBWEoszeFKkizLsixTShljjDHVhwbUlvuwgohiaK1b6wkqT9m/kMd5O1hh8ApCirHwCm2VefrUFycJWyGkGNtsqXNcrVHIsiwMwyAIlFJJkkRRFMfx3oU6hpMjlLz7kADOaDRJvHR3CfdtxyiuZh+tdRAEtkYhjuMkSaobIiE2epLBuyoMqqu48H3Fk7eDFZZbQUgxNlqh52EZH+Ft/oqEFVzk6iZprdM0LVscmr/ufjS48dEamAPqS2pfRx7sviskrCCkGKwwcAUhxdhhhZYREwt8GVAnailOnCVayunpCi5ysumh7JpQW1hbsjYhfYnLYjTHIpZrzP0bd4/7em7QHX8XgEDNy9iQU8SpdVrOP3dqeZc5Swg56x6Nk0GhlU0PJX3qsJ6/guq6MG9ahlfF6Fr15EucXmH9DfFjBSHFYIWBKwgphogV+h+etkKtY8QSb4eEs+6w07JX/AkKd6sTXlaCdX7J9lTb8Tuk19Kx9hKAkUafRnyrfB/A4y6i/gSFqvqkp4cx9cN5rL0EYAVHP414nI2cvHukrTyotTVs3EEBAIAjcDIoKKWCIAjD0P5sIwJBAQCAxTk8kMMYk+e5/dndrQAAQDKHgwIAAFibq00PAABgAwQFAADQiaAAAAA6ERQAAEAnggIAAOhEUAAAAJ0ICgAAoBNBAQAAdCIoAACATgQFAADQiaAAAAA6ERQAAEAnggIAAOhEUAAAAJ0ICgAAoBNBAQAAdCIoAACATt/buwCr0FrvXQQAwBEVRbF3ERbmZ1BQMt6qMq9IKIw3tNbV/Vn7FaNU994ae5I361B4f5WnX1MlBoXLy8vr62ul1Pn5+U9+8pOu1R49evTNN9/cv3//008/3bB0I1RPwbUlgDRFUXCix2QcPB4T99ZeXl6+9dZb5+fnSqnr6+uHDx+25gB76T0/P7eR4vLy8sc//nHt0RoJW0pimK95PuIMNc3ae5J36lB4cy0v94O4TdJan5+ff/PNN0qpR48ePXnypFnC2vLqU8ol0rarhsQwGZefpWwfFBb/E5CDd9bycj9IHPXw8OFD+4OtS7i8vOxf31Y/1OiGxcs5R3FLZvHm8GlbAAzh5dURJVlBwWaCWr8E27hQZQPEo0eP/vrXvz569Oj6+vr+/fu1dYqGNQs+nceJAc6xx+HepQAgi8TOjDXVNoXSw4cPnzx58uTJE6XU+fl5sx/DkPOdqPTQ7PmohJUQNU5/i3K68AC2JKtGoVWztsD2Ubi8vCyKwg6RGFKj4Fwdg4vVDPbys16B5VzeHHpTdiHnnTogDk4sy4Gg0PTkyZOHDx/aYQ4//vGPy+GU/nE3McyitRqzpVSYL4idua+ldv6WbyKh0HuygoLtnVDrvdjaV/Ek4Z0ZxzpQYig3TfA22jMj11Qsa5HDics2FicrKCilzs/P33rrLfuzbVCw6SHLsjiO7XLbQeGvf/2r/bV1UiZXGhrG8jwx1DZH5NZ5cCLecRNIV/2W2j/sZyxI4inv/v37ZVNCWbwsy8IwTNPUGFNbpzaJghI84dIa5EzJUF5+Jl6Hus5rA+YA3uzKV/tDjoaG/mIvslE9L+LoTlvb3I/P0q8z9s/B8nKHiKtRsB0Ozs/PbYvDo0eP7HJjTFEUNiUopcIwtCtfXFxcX1+/ePGi9jq+1ig01eoY9i6OUtO+zfSsL2OjZLHdOIb8G2/+l1FZ58o195U01T1PpQKWIunzrJQaNjOjXe3rr79+4403lFJBEFxdXXnwPW8RO2579U+PK8aQ01nvl6RBf67nrwwrautf2W6HL3fS10r1l/jkCnVjZnjcYo/tdYGcul1L1VQ1n7j23j7yybaVlztE3CZprb/66quy20HtV+vx48cff/xx/5mouVDalq7EvaAw8Jx+qja1/tDkS8XICvN1d/gKF7whIWB0UGg+fWxQOMh332GHkCtBwcuL4kxe7hNZEy51zcxYW3J1dXVxcfHixYurqyul1MXFha1aqPLvrRqo2OkegAv3S2hdc9kqipPPnf8Nz9mLXzE7K5zYdmf3zFxtR9dyr91ylO51QoBPxPVRaGrOzHh1dXV1dfXmm2+WPzx+/Li2TnN4pGcDJn2wYD+G6jpLvbO11uuehm1/27yxltuDpOfqvtSforMCZpJVo9CqOeuiVX66bEvEb3/729ZHD0jCd4jTZVjjzLXO2bD/u/UCX779cuS9MXbbX65vj1vXTlm7n2SwGQeCQtPFxUXPr1YzQXNMSzHzct5sg1jz25JnlxE20ncAABTRSURBVD0Jm+NltLKH4KxNu3sYF/ZK3LP+mBOahC8PcJespoeBMzO2JoOa4wyPbCW3snGRUlVnbyQlrKO4vfjhJHuc2EPF7rchu26BXiCVg58cgPWIO7bsTEq2VNWfsyyrTs6otf7www9tc4PW+uLiIs/z8kX4zKhtd8LQcYOzL+r959bFr+sDX9ChPDGqqNO2a4NRFXKUFQmtD80cgzp0Lw2eXmnBcwLn2C5e7hmJmzRkZsYXL168+eab9qFaSlDHHh5Ztdd8hS3Ll7rVzYZBYYNr6vbW3ij/olWXnohQXadrhcV3VP+Q1FerLXRO8PJyuAgv94yspgel1PPnz6+vr8/Ozs7OzpRS77//vl1em5nxjTfeKBsU/ve//zVf5+BND7JskhJwEjtwKdW2hh5dzTd7VX0B04jLPlrrs7Ozb7/9Vin1/vvvf/bZZycnVjr4vR767Tvf+4neWGP/yuLzCS73ak6cqderIRi7vhO7q9WQioTmU4pTS4Y/9/RqmzRAePm9eRFe7hlxNQqqUovw+9//Xin1/Pnz1tUePnyoum9CTY3CZjo/GDJ7Uw7j7pUM6xlYkVBT69446tCa0qV0/c+dl9dC9JAVFGwm+NnPflZdaGsXar766qvf/e53HKwnLTb8YeytdFxOCVjb0YZUFGtuckvy6P2Qyh0SBakcmEehNSj89Kc//ctf/tLzrCGfBHLGCaPmV26z8Yj5HQfoy58bYFrxRm2X8D0w0/ytKyr/3VrzEzr4TmyAA0Hh3r17tSX3798/Pz9/6623ep5FCCiNm2tl768autrYeqyvnUCLNTLo0Ndsu8EK7Q4H5EBQaLKDJ+3UzuXPH3zwQTU6MDPjFONTwoKnMN2sQFWFIi70au60GvZeyct9tXotjpvTS2NZsoKC7Z3w/PnzajeFZo3CBx98UP5sg0KzgoFYULXPnReGOXn6nhYX5lwVJFehn9yu7vX12OcOf6KES+zUrRv+Cm7sh1WMuncrvCOuEunevXs3NzdRFCmlvv7665ubm9aZGZVSjx8/vrq6spMoMDzypJXGJswZsjj+zH56sKX9+jPyZVtfpphw0q9t7PxL13Jmhp81hj1KzmNNy5R2cpLo+vNzijVxuOyxT6Qnedk0I26Tsiz79a9/fXNzY3+NosiGg9rMjDYKXFxc2KAQhmH17pFevlXz1XfLErUIE+brnXH5fPlKtbOtrr38QkFhkydudrHcICiI3fZF7F7aVwUoj/9t0l/9KZxae/l59WnON7AvpVQQBPZnW6/QXOfDDz+sLq/9al9kzTK66tVuuR2vNf+fGrPOJn9t+GprvMKoJ6rb/84v8Db7pP9FJvyJDTZczg5cowBbb9Krp6CDl1cfWfMoWGX7QlmXUFvh6uqqegPJrttM16xTWJcUdvz0crtizMC5QlI9fL+NvzvaE6/849OVcq5h9+oE1bXz53yyxr6dd/YCZ9QjkdWZ0WaC8oYO5cLaktotoK6urpovVfhX+TOf1sXm5zx9rKvL8B1cW63cS2u8OQu+56uWExOVWWHT3pQMiDgMWUGhVbNGoerx48cff/yxbX2oYsKlO/aI/+5UIQhhd5f8y7Ar5fRPLYa2Tcn4sj1i27hwnBPpUTkQFGrVCaXyTtNff/31G2+8UXv0QCGgX3M+iSVHIgwccbd9RYZwPXujOLWCELVyyi/wZO5t2qhPse7dwEGZg6oF30nsozDE48eP33zzzQ8//LAoimZKgFKq70YMzXVVMb6x070TqDtstJLfZONKOTFReWY4/Wmn14K/xA3k0FqXYyCbv1q2LqG1IqF8VnOhtC1d0YBP7BLzDWw8aH6bP7dBOSXvt41JLltJWiH10lV043rVnK5jOM6Zto2XwyPFNT0EQRCGod3RNh/Y/1YnXPr444/tyi9evCifWAsN/r1Vg5yKCNV5+tYuy1HRzgK3jDtiT3eDoNeCdyRmH2NMOa6hLF51wqWTFQaHqFEYU9HX0aaw5Xwtx6lR6H+R7WclEp5aKN5Yi9coqAGv1rlCZ2Lw7Hw7mJc1Cq5u0pdffvn222+fnZ0ppW5ubt59993PPvusfNS3t2pS49+ADgcbT+zmRC4hKGyM4k2wRqn6X/P0Xxwx1MKnk3ODb1cfpZS7QUFrfXZ29u233yql3nvvvT/84Q8+1CjM7g004e4JG15QCQrTXnz7iXq3RPGEmBsUStMHZ8o/RQ/gZVBwddSDUurdd9+1P9i6hC+//LL6aHMSyh2K2MMOSaj961r3VcfjE/+23AIfLXVVKOgC4oXjpIR+4/bD9HNR61mRwRQCOBkUbCb4+c9/Xl1Y3kdKnLZDf/i1f+XLv6+XNIHbtctVh0sdBlr4IzPqFHfiREeA2Ju4UQ+T1YLC8Js7lJUN5VOmLbk7Y1p1SdFY0lznIEv0jL1RLlyvhGrAOs4t0cLK42IJj7Ok9b1QjSWr/HU99lnzztjrLfGPP0HBdmwsTWhraD5l2pL6d7g7R4++/d+pZ223pLh7dhjyrHL9LcvcXLj432pNJ0u98vZ7rPnOyjnqyiViS9hVGSOnhIsvaX0vqgt3L+Hd3xc6Yy+7xMu44E9QqGm+W7t1U2j7u8UCgxuxrGVr6bviF2rYUXIs+V4seb8JaT3MjsfJPgq2d0Kt92KzRkF0Z8aiaPnXta7SA/9NLo28Fn2fzDz58u7AJdPPRa1nRWmn7kOSOJDj8vLy+vpaKXV+fv6Tn/ykdZ179+7d3NwEQRAEwYsXL25ubmrDIwVu13RrzaNwd/UR60++8s25ZE547jYbdfIF57/yLjt8MzILKbNUa6tt9YidwDwKlm9XH6WUwKBweXn51ltvnZ+fK6Wur68fPnz46aefNlez8yiUHRg/+uijJEmqjzafIm1L51q+8WLgSWHmCXTLa56EoLDIyxIUNiawSJup9Wc8sR+YmbHGy6DQUkW/L6XU+fm5/fnhw4etJfzoo4+qy2u/2hdZtZByKdX/78SDJ549arXFnz7ticOfNXO7ul5zkZdde9v3/SewnAKLtMu2d+6HE69yYF5efST2UbD5QCll6xIuLy9rK+R5HgRB+Wv155JuWK28kpxq0uvt0FCcagvX/n7TWm+7vNxdOK7TXRD8+z59eLKCgs0EtX4Jtr9Clb2TZPlreQepqmYmWr64Yg3rAXSqF6Ru/Cv2u+xNvpCfTD+rOtJR5w9f0/BA5Ufmzn4Y1GOavoeecmB45DfffNPzaBRFv/nNb2zrQ9WQKgTP04PdumFVKeXnX/fNJQCsoTj8tVkuuihC7RIUavUBpTiOW9e/f/9+6/J//vOfr7/+ulLqH//4xw9/+MPao56HgOGKYlS3x8pe65kbqvaU9tfveCJXhQnYaRhh5hwGWil71hj6OpxsfSer6WG4KIpef/31jz76qCiKZkrAHfPqAydP3jBvaocF7dv6ALe4HcgW+dCNi6WkhAMQN5BDa/3VV1+V3RRqv1q2LqG1IqF8VnOhtC3dx06dOu/WLqw90HHs0524Nmy/07Ykp6hySjLCPomcM2obL4dHituk+/fvX19f21JVf7YNFrZ5whiT5/k//vGP6hOrocHLt2ox+w0AuY0Lu1zzel7BiWsDQWEbckpy2p41dpxgO3h59ZG4STYf2J/L4mVZFoZhmqbGmJMVBtQonLZnXFAzr3m18+Ps6aScuDaMLaQTG1WSU9qdSyKjta4XJ9JeXgYFcX0U7PzN5+fndnLGR48e2eXGmKIojDFK1edWuri4aL4xhx4eOcSuA5kmt6S2PnHY3S5G91QYfouNyf/G7wBXObIfdkgJUndFB06khyQu+2itz8/P7ZDIR48ePXnypKeEtubg4uKiNpWCl5luLXtULdydJHZQfcD8t7PrIiDny2yTru4nwTUKS3Szn2uJm6yeKIgDF/L1cEYdxsurj7gaBTVgZkbLzp1wcXHR+uhBZ2acYO85UoaMnljvr0tOCerO102lRtZVbFAjstT34EXGpSyyFY7UfGzOuysfRpEVFAbOzKiU+vvf//7b3/62J7jR9DDO3nFBtTUrLHt2ZpQkMJqAMwN25+rMjG+88caLFy96nsXMjFOMnJ1plSJwNccehNct7YMzJJRSjs7MGATBxcXFj370o56/QgiYqLnfTkaH1l29d+BoxQSHOCjOh5jBgRqFpqurK3V708jy58ePH1ejQ7NGgegw0bT9JqBy4iRfQ4OL20WGWxLnOixqh6BgjLGjHJts74TLy8tqNwU7TrLq8ePH5c9lUKhVMBAL9jfmrlQAWWEZnPqwNHEDOYbMzFhl6xUYHilaW1YYfklY/OJRvqBzl6WBBXZuu0p7ldzdPXYHJz0BvLz6iGt6+Oabb+7fv182HFT3eJIk1dqIL7/88ubmput1aHoQxIVmCEhApcJ0nN+wGlnDI5VSz58/v76+Pjs7Ozs7U0q9//77dnl1Zkal1Hvvvff222/f3Nz897//vbq6+vLLL2uvw/BIWSQNsmKcpGS8O6NJ+nDBS+IqSbTWZ2dn3377rVLq/fff/+yzz1pLqLX+05/+9POf/1wpde/evZubG+714Ibbt2Zu08O8DhD2ac4dEN43PZQ23gSH9xinNWG8bHoQt0la6z//+c8/+9nPWn+13nvvvT/84Q/9UztL2y68orWaGRRqb+6kuODiteE4QUFtuBWu7i5OcSJ5efWR1Ufh+fPnSqlaLPj2229rS25ubs7Ozso+CmdnZ7ZqAW6YORqi+SGc9IK+fZRvuXrZg3dXF3hDXB+FJtsMUXVzc3Nzc2P7KNgf3nvvvdo6zXs9cPcHUYqimHT7yO6TKS21fjlKZwV73HLoQjAHZma8d+9e6/Kyese2RHz22Wetj8IfQ97TBSdv2OAQGl8LcqgKA5+3lxMU3CGr6WEgOyCi61eL4ZFuGD5yctTbN2dA5pbHycm/5Xq91+ydWZSNvuvsiq2DCGchOEjWzIy2L8Lz58+rnRKaNQpnZ2c9MyhYxAJnrHQNGHXTCrFHS7VgroSGlXamhFAl9jgB1iSuf6Yd6xhFkVLq66+/Lsc91mZm1Fq/++67//73v40xSZKUIyrLR5uvLG1Lcad7cPdZXvPeVZz8el3/iuzFrpvYk3xAdLizu7zYV9gXox628Mknn/z6179OksT+ahODVZ2Z8U9/+tPbb7+tlMrz/Ozs7JNPPqm9jvy3ypXjaaNyzqtXONzO7HkRrWde8ATuzKIoarO1DirkkK2YvbtOvby4ndnkRCGVO+X0j7j9rrUOgsD2dozjOEmSZglry6tPKZdI264mJwqp1ixnyys3s0JRDCnAcXbmyVfY4E/sa9nuR2tvrPCdaTlRSOVIOZ0o5FgSh0eW7Qv2h9YhElX2vlA1DIaUr/o1sVzU9ytO8fIkVVPOyF7Ozj75M36E3QXMJyso2ExQ6+rYDAo2QMRxbHst5Hne7B3JvR5cVb5TvGVtWtLV4fEBB1Ylro9CU2uNQhRFSZLYrgxBEDTnYBhyMuXkIhTvC8YrhrVSARjLgQmXmrUFto9CmqbGmCzLwjC0P1TX4XzhB079ALAvWU0PAyVJEkWRDRDGmDRN8zwf9Qon6xs2WOGkIa/gRDkPUsiB66xdDD925rRCVttltnlDPd6Z25fBiXIes+FP1oRLdnmWZdUVulbux8yMAADMJ65GIQiCMAztz2WdgVKqOtuS7aBQtjWU61fRmREAgPnEBYUsy4IgsIOd8jyvXuDLcBDHsc0TdrUgCJo5gOGRTqAP/wRdO43+HIojCliB22eWWiNFaf6kNE6sIKQYM1ewp/WuFexzdy/kkBW2LEbrag7tq1VXGLIf2FfSVhBSjG221DkODI/s0dN9gT4KwDFRqQAsy8Pso47aMRUAsDv/rqp+BgUAALAIcZ0ZAQCAHAQFAADQye3OjP1qc0J3TREtRHMWalGyW3EcT5sCawNlIY0xwt9uK45jUeUsp1fvmRVNDml7r8aJo9GJz3WV5POkW1eccZoTE/khiiKHttTeKTtN070L0s7uzCAIbDntxBXS1Aop/B0viiJNU1Fvui1PuQOjKNq7RH2k7b0aJ45GJz7XVZLPk25dccbyamOq7NG/dykGsac8sR+AoiiqZbOfh12L0652bZN8qUvTtLx+yHnTq5cKse9yIXXv1ThxNDrxuS4JP086dMWZwOc+Ck7UpCmlwjBsplE5yrpo+6utT5NZ+1et6wuCQGYhLWOMwDe93IGS32Ulde/VCD8aHfpcW8LPk8qdK84UeyeVtSilyq8dQRDIDKFFUQRBYL9qKKlJuar8Mrd3QU5TIr/D1ch50+3XteoS+TtQzt47SfjOlP+5ln+edOWKM43PnRmVUvb0F8dxGIaFvBkj4jjO81xyiq8yxtjbeZd1gDJlWWbvE+ZVZ6I9uHJkSib/aJT/uXboPCn8ijOZw0Gh7KFdYz+Q1TcpyzKt9S59pHsKmWVZkiRCDqb+nVmuo24/A2mabl/PNqSQ9qy3e03vkKLK53NV6iaEHI39dv9c9xN1nuwh5IqzEp/7KFQJ/KzaY8jcUkqFYSjtU2pVd53ktkx7x9E0TWUWD4ci/2h04nPt0HmySuAVZw6HaxR6hnrbYcHV9ynP8106wvQUslnCIAj2SqD95RSS6PsH99sbjgv5cDoxD0GVLW3tdqxubYIooo7GVnI+1/1EnSe7yLnirGW/7hHrUo6M9SopqZ10irtdsezOlFZU2zQYRVF6197lOkHUnqx2ZxPetc0StfeqXDka5X+um8QW0rkrzihebUxVrWOOzGOrSnIh5e/M1vAuf1iztJ1Zdtt24jQnbe+VXDka5X+um8SW08WdORx3jwQAAJ2O0pkRAABMQFAAAACdCAoAAKATQQEAAHQiKAAAgE4EBQAA0ImgAAAAOhEUAAAYx975qfzVuVnbR2HCJQAARivvkGJvJu7xxZQaBQAARsuyLM9ze4dur24B1UCNAgAAU9i6BOF3Cp2PGgUAANCJoAAAwBS20cE2QOxdlhXR9AAAwGjH6cxIUAAAYJw4jpMkKS+gZWjYsUjrISgAAIBO9FEAAACdCAoAAKATQQEAAHQiKAAAgE4EBQAA0ImgAAAAOhEUAABAJ4ICAADoRFAAAACdCAoAAKATQQEAAHT6/8ogUYaLu5lgAAAAAElFTkSuQmCC\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from ROOT import gROOT \n",
"gROOT.GetListOfCanvases().Draw()"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.10"
}
},
"nbformat": 4,
"nbformat_minor": 5
}