{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Interactive Usage of 3rd Party Libraries: [Fastjet](http://fastjet.fr) Example\n",
"
\n",
"This notebook illustrates how to interactively use 3rd party libraries. The example is built around [Fastjet](http://fastjet.fr), a software package for jet finding which is widely adopted in High Energy Physics. Nevertheless, the concepts discussed are valid for any C++ library.\n",
"## Some Context\n",
"Jets are the experimental signatures of quarks and gluons produced in high-energy processes such as head-on proton-proton collisions. As quarks and gluons have a net colour charge and cannot exist freely due to colour-confinement, they are not directly observed in Nature. Instead, they come together to form colour-neutral hadrons, a process called hadronisation that leads to a collimated spray of hadrons called a jet (from http://cms.web.cern.ch/news/jets-cms-and-determination-their-energy-scale).\n",
"The kinematic properties of originating quarks and gluons can be inferred by the ones of the jets. In this example we will build jets starting from an event which is an emulation of a collision of protons at the [LHC](http://home.cern/topics/large-hadron-collider).\n",
"## Preparation of the environment\n",
"ROOT needs to be told where to find\n",
"1. The interfaces of the C++ entites present in the 3rd party libraries\n",
"2. The implementations of such entities\n",
"It is therefore necessary to update the include path and load the fastjet library."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"gSystem->Load(\"libfastjet\");"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Preparation of the Input Particles\n",
"Here we convert an ASCII file in ROOT format to then fill the particles necessary to Fastjet.\n",
"Notice the inclusion of the header file. Given that we instructed ROOT about its location and the Fastjet library has been loaded, the function can be smoothly JIT-ted."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"%%cpp -d\n",
"#include \"fastjet/ClusterSequence.hh\"\n",
"void fillInputParticles(vector& input_particles, const char* inputFileName) \n",
"{\n",
" auto ntupleFormat = \"px:py:pz:E\";\n",
" TNtuple input_particles_ntuple(\"InputParticles\",\"Input Particles\",ntupleFormat);\n",
" auto n_particles = input_particles_ntuple.ReadFile(inputFileName,ntupleFormat);\n",
" input_particles.reserve(n_particles);\n",
"\n",
" for (auto i : ROOT::TSeqI(n_particles)) {\n",
" input_particles_ntuple.GetEntry(i);\n",
" auto v = input_particles_ntuple.GetArgs();\n",
" fastjet::PseudoJet particle(v[0],v[1],v[2],v[3]);\n",
" particle.set_user_index(i);\n",
" input_particles.emplace_back(particle); \n",
" }\n",
"}"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Preparation of the Histogram Displaying the Jets' Components\n",
"We now need to:\n",
"1. Cluster the jets with fastjet\n",
"2. Extract the clustered particles from the jet and display them as 2D histograms\n",
"This function is flexible: the jet algorithm, its radius and the jets' minimum transverse momentum can be set by the user. This is useful to visualise the behaviour of the different algorithms."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"%%cpp -d\n",
"THStack* getJetsComponents(const vector& input_particles,\n",
" fastjet::JetAlgorithm jetAlgorithm = fastjet::antikt_algorithm, \n",
" double R = .6,\n",
" double ptMin = 14.)\n",
"{\n",
" // create a jet definition: \n",
" // a jet algorithm with a given radius parameter\n",
" fastjet::JetDefinition jet_def(jetAlgorithm, R);\n",
" fastjet::ClusterSequence clust_seq(input_particles, jet_def);\n",
"\n",
" // get the resulting jets ordered in pt\n",
" auto inclusive_jets = sorted_by_pt(clust_seq.inclusive_jets(ptMin));\n",
" \n",
" auto hs = new THStack(\"hs\",\"Jets;rapidity;#phi\");\n",
" int i=0;\n",
" for (auto&& jet : inclusive_jets) {\n",
" // get the constituents of the jet\n",
" auto constituents = jet.constituents();\n",
" auto h2Name = TString::Format(\"JetHist_%d\",i);\n",
" if (auto oldH2 = (TH2F*) gDirectory->GetObjectChecked(h2Name,\"TH2F\")) {\n",
" delete oldH2;\n",
" }\n",
" auto h2 = new TH2F(h2Name,\"JetHist\",48, -5, 5, 48, 0, TMath::TwoPi());\n",
" h2->SetFillColor(i+1);\n",
" h2->SetLineWidth(1);\n",
" h2->SetLineColor(kBlack);\n",
" for (auto&& constituent : constituents){\n",
" h2->Fill(constituent.rap(), constituent.phi_02pi(), constituent.perp());\n",
" }\n",
" hs->Add(h2);\n",
" i++;\n",
" }\n",
" hs->Draw();// This creates the axes\n",
" hs->GetXaxis()->SetTitleOffset(1.7);\n",
" hs->GetYaxis()->SetTitleOffset(1.7);\n",
" return hs;\n",
"}"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Clustering and Visualisation\n",
"We now prepare all the elements needed to build our plot: a canvas to hold the graphics primitives, the containers and the stack of histograms containing the jets' components. First of all we uncompress the input file if needed."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"%%python\n",
"inputFileName = 'Pythia-Z2jets-lhc-pileup-1ev.dat'\n",
"import os\n",
"if not os.path.exists(inputFileName):\n",
" import urllib2\n",
" response = urllib2.urlopen('https://raw.githubusercontent.com/dpiparo/swanExamples/master/notebooks/Pythia-Z2jets-lhc-pileup-1ev.dat')\n",
" filecontent = response.read() \n",
" with open(inputFileName,\"w\") as f_out:\n",
" f_out.write(filecontent)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"TCanvas c;\n",
"c.SetTheta(90);\n",
"c.SetPhi(37);\n",
"unique_ptr jetComponents(nullptr);\n",
"\n",
"vector input_particles;\n",
"fillInputParticles(input_particles, \"Pythia-Z2jets-lhc-pileup-1ev.dat\");"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We are ready to go now: we can play with the parameters of the clustering and visualise the result. Each jet is visualised with a different colour in the Rapidity-Phi plane."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"#--------------------------------------------------------------------------\n",
"# FastJet release 3.1.1\n",
"# M. Cacciari, G.P. Salam and G. Soyez \n",
"# A software package for jet finding and analysis at colliders \n",
"# http://fastjet.fr \n",
"#\t \n",
"# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package\n",
"# for scientific work and optionally PLB641(2006)57 [hep-ph/0512210]. \n",
"# \n",
"# FastJet is provided without warranty under the terms of the GNU GPLv2.\n",
"# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code\n",
"# and 3rd party plugin jet algorithms. See COPYING file for details.\n",
"#--------------------------------------------------------------------------\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElE\nQVR4nO29Qa9syVXvubblKf4Ar74AE0/oO/AruXzy0jQGXLLEc7VFSW2JRm0XzQReYQnXK+N7roXB\nTwILJjRVlp6eZImirbKQiirACHHzuEolBiUmnvAF/AU8ZOB8gzgZJzIi1n+tiNg7c2fm/zeoumef\n2BE782Tm+ueKtf572u12QgghhBBS42OnvgBCCCGErBcKBUIIIYSoUCgQQgghRIVCgRBCCCEqFAqE\nEEIIUaFQIIQQQogKhQIhV8R2u91ut6e+CkLIOTHRR4GQ62GaJhHhu54Q4ocZBUIIYpqmIC8IIdcJ\nhQIhhBBCVD5+6gsghJyeWLiw2WyWO4UQco4wo0DIVXN7eztN0+M90zRFBRB+Ff49TVMUBJvNRjuF\nEHJ5UCgQcr3c3t4+ffpURJ48efLs2bMnT56IyOPHj0Pg32w24UgYcHt7G065u7uLp9zc3IRTTvMA\nCCHLw64HQq6IrOuhbILYbrePHz++ubmJSQLzlM1mc3d39+zZM+5BEHKRsEaBkCslSIGYMwhsNpub\nm5uQM2iahxByqVAoEHKlxABfjfTb7baaIQgyYpqmm5ub29tbZhEIuXi49UDIFZFuHIQtA21k3Eoo\n9xpiZUM5mBByebCYkZDLJPQshApEwLNnz3Y1QOC/vb3d7Xah+DEWM3IDgpBLhUKBkCslSIEywOP7\nQcTfbjab29vb7Xb77NkzETEVCSHkTKFQIOSSSUN+Fv5DaM82EUQkuCNoE+LfEkIuDwoFQi6ToAPu\n7u7C9/7Q9yiHbQ5h4yA6Jm2321CREJIEgM1mE08J0zKjQMilwmJGQi6WGMUjqUFCoCxpjN5K2YBw\nbjlneQoh5JKgUCDkwkmrCrQSxRjmq/E+HExPT+sYKBEIuWwoFAghhBCiwhoFQgghhKhQKBBCCCFE\nhUKBEEIIISoUCoQQQghRoVAghBBCiAqFAiGEEEJUKBQIIYQQokKhQAghhBAVCgVCCCGEqFAoEEII\nIUSFQoEQQgghKhQKhBBCCFGhUCCEEEKICoUCIYQQQlQoFAghhBCiQqFACCGEEBUKBUIIIYSoUCgQ\nQgghRIVCgRBCCCEqFAqEEEIIUaFQIIQQQogKhQIhhBBCVCgUCCGnZ5qmU18CIaQOhQIh5JRM0zRN\n0263o1YgZJ1Mu93u1NdACLlGgjJIP4KCYjjdFRFCKvBtSQg5NqVESH/FDyVCVgXfk4SQ4wEkQtMY\nQsjRoFAghByD1vDP1AIhK4FvRULIsnRnCKgVCFkDfB8SQpZifBOBWoGQk8M3ISFkfmasM6BWIOS0\n8B1ICJmTJUoRqRUIOSF8+xFCZiDaJS30kUKtQMipoDMjIWSIYK0oshMxmh5HvBdp3UjIqaBIJ4R0\nso/c2WdIxW8xGTb0mcO8AiHHh+86QkgzikQ4GJLkAAwl0bE0P7gIORoUCoSQBpLYP+G9BmsAUwuE\nnAesUSCEeDmsRQhaoT5wrxJAVcFQmGfJAiFHg6qcEGKj7zVkaYNyWDWvMM/2AfMKhBwBvs0IIQhP\nOUKSPKgOy8XEjB871AqELA3fY4SQOmnZ4H7ToT5QRHz1CovUIVIrELIofIMRQnKqnQW11MJB36OV\ndViwVYFagZDl4LuLEPKA2Xy4Ty2ofY/zNkP6oVYgZCE+fuoLIISsghZ/Ai15kKUWjmp4EPogqBUI\nmR2+rwi5dpwSITNYPOFGA4B2TITMDjMKhFwvXRJB9v9W+x4tB4UFiaWX1AqEzAXfToRcIwMS4eD3\n0DLhlB8v1AqEzAXfS4RcF+0SwdxlACYK1AqEnD18IxFyLQxkEUZMFKgVCDlvWKNAyOUzvNFQzSv4\naxtPBlshCBmHbyFCLpmsVUF7vzt8miXre3SYKKylAYFagZAR+P4h5DJRYn8evFubHmtzZmOWsmoe\ngVqBkG745iHk0vDcxim5TXObuyK86YOsUCJEqBUI6YPvHEIuin34F/c+Av6t86YP9ZFrg3ZMhHTA\nYkZCLoQsCupf/T19j2j7oEgtnEH0DU/O3/7bh8LUAiGN8A1DyNmjfVEutIKz79Eb+M/iC3oqESK/\n/gvPr/yyCVkPFAqEnDFjN3u8H7IeO8V5qUqECLUCIU4+duoLIIT0ME1TSKHjaLfbxS0GY5dBZLoY\nlRCeHE0iBP723z5M6jkIISoX8rlAyPXQkfA/31aFVsoswq//wvMgqSCX8sAJWQ4KBULOhpGagPNt\nVXACNhpKrZBKBNY2EoLhO4SQM2CussHzalXwY240RK1QzSJQKxAC4NuDkFUze2fBWbQqdIC1QtAH\noj9wagVCNPjeIGSlXGpEXw6tE1IabojFJ5yQHAoFQlYHJcIIMbXQJBGaTiHkqqBQIGRFzCsRrjad\n3nRb7b4MBCHXAy2cCVkFs0sEEfnKV75ynVrBaStZrWn423/78Nd/4fnrfN4IqcI3AyEnZolyxa98\n5SvxxzfffJNv85SYSKDFAiEeKBQIORlNGXL/zRdSlRCgVgj47ZiCwTPzCoQItx4IOQlNEiEEfhy0\nNIkQuNo9iIhmx1TmFdJcArUCIUKhQMiR6ZAIgRDsy3PTkW+++aamFeS6axuBxULUCtXtBmoFQvgG\nIORINO0ggHgf9xGqWYQ333wzO5gewXsQF9yW6bFjws/MRT4thHjgq5+Qxem7jRPWCqJvNMg+tVAd\nFg7itMRFfizQjomQPigUCFmKwehS1Qox9uNdBrG2Icy0xKV+MtCOiZBWKBQImZ97ibATmYbiStno\nCH5Mj8t+c707LXHZWkFox0SIGwoFQmZmmqaDmznrWsETsUKw14J6KR2yCUFawqMkLuzzIXVQMJ92\nWiwQEqBQIGQ2HhIJ+S/qNQHxhs8eraD9Ni1HqM6D/ZfA5BcmFLLYr2kF2jERkkGhQMgMqBLhYYSk\nNQH5yHAMVt2DcI7PlcO0RDmynNwz5xmhmShkWqHJjin842KeIkIAFAqEDJHlBlStMO3/4RATYKGF\n+h6jVrgwiSCOxsiwBSM1JSGFVggSgbWN5KqgUCCkE5AbyGoUHo4AJSGuykdQsjDS93jZ7X9mwQFQ\nErLXCtWR1ArkGqBQIKQZz0ZDaHmoDNMPOt+M43ZMZb794u2YqtsKIg/7QaYdk5mWmPFqCVkVtHAm\npIG07xGx05MHWWqhRSLcTwC7FaLTc3VAetuIPpfoM/UoTAVBmQbQdh+yjQagFc5aRRGCoVAgxMtB\n3yOQAvvYn/dJpux6JMLD2YpWSO2YwOlBEJhKQvRbTpxpRAQ3bshSC2UWobr7kCmJ831mCAHwZU2I\nDeh7LMsR8sR+9Sxx2TGZd4ykHdPsgMLGAEsWyLVBoUAIoqkcQe3L3xXjkx9BN78pJpazYzJdoi9b\nK3jKG7VfCbUCuSwoFAip09T3aCYG1NpGyaVARZpY3o7z2jF50hJyiY2UKZpWiDqA1o3kemCNAiE5\nlVBdDfPJEXtzGkiN3UF1YWXYrjJ/1q3Q/dU/7kFU0xLVvEI60twckfMMmaUUyMI/1gpytlWfhJTw\npUzIPdFLYKTvEVs1G3ZM+L146O0Ignd5xNP3aKYlqqvImH/DytG6JLIB8UfaMZGLhEKBEHc4F0dE\nd1g1a/WPqEsinuveaOgI53iXwVy6yb/hXDATA6mYYG0juUi49UCuGu8ug7gNFnew/jE9vchAqLsP\nmZiAaPHe0/dopiXw0k7/BjzJ2jAvGOxBhIPn+KgJSeErmFwp3nYG0XUD7Ht07iOo16bbMZlbD2bJ\nQnl6NgCnJbqXlnPehsCYtY3C1AI5WygUyNVx4K5ohfP7kdpvYd8jsGPy9D1Ky22cyuDd3feonV7+\n1r90tvpFfuxgl2hhaoGcLXzhkuvCcDVIj5vuivF0ve9RvTvU/kj1ps/3wxx2TAI9kbrtmDwmCk1L\nS7GjcamfPLj+kVqBnCN81ZJrYWZ3xWRwv5KQmomCz0Ehju/bRBBfOB+xYwJLy6VrBZwKutQHTi4V\nvmTJ5TO/u6K4axvdfY9Ax3RrhZFyhHSA88empeUStULcfcD3k6RWIOcFX6/kkkl3+s3EgMtdUZTY\nP9736Gi5BNfWV2DoDOf+tIR/abksoVAtUND2lX7wzx9+4Zd4Z2pyNlAokMsE3C5ZrRtwVBfmp+eD\n1MyEbasA0xLm+9SsbRTY9zjeJeFfOh68jA8fcBOpTCtM0/SDf34YRq1AzgUKBXJRxHBuNx/O7q6o\nnF5ZWtDquGACzSwip+h7dNoxpWLiYj52PLePivtKqUoIUCuQs4CGS+RC2IfzXfgBjES9DPdnH7gS\ntdoxgR1oVIsgh1qkUSI4PZGwHRPIDQA7Jr8TlPMRnRHmHR+kSCSk/OCfP2S9Alk/fI2Ss+dAIiRH\nzds3D9YN4NyAK6+gT94qESJmqt9slChPzwZot6P0OEFd5AdOdfchdEgGifCFX3pe0wpf+CXaMZG1\nQ6FAzpi6REh+bd6+edBdMZ/w8Lc9WqHFjqk71V/VCk12TCNOUOX1XAYxtZBKhEgQBFmNgiRKQi70\naSEXAIUCOUseJMI0qUJB7rUCyvYPuysCMbGDQqNi/ST56lrZvDOcix57Mq3QZMfkqawU3x2kLgyt\nFiESUgulaBBqBbJiKBTI+TFl4gBoBdNBQSxPJIe7opaWcAVUsLTkWgFH9/S4OMySw2xNdkx+iXDN\nN30AQkHgNoSwvJGsEgoFck6oew2lVshqG8fLEWA4b7JTTK7RpWNiUqTPXdGzTVAFlCN4Zrvmmz5U\npUBMJFArkPOCQoGcB0Y5giRaQaltFOm3MaicW45pb1Xwt1wuZ5ZsGjt2l1XKsI45X0rLBCkKFKrd\nkkI7JrI+KBTI2sn7HvEug5hiIvwjjMx+Wxx0KgnplQhg6cM5R3yapV0r+FsVnGWVS9z0wbzB5mkJ\nWqFajhDItEL542ofGrk26KNA1kslixCqF7ODbgcFEZgeKGwMZnRQMIwdMy1SyI7q7oPfXRETJgcO\nCsDm0untCH4b5m8KimHp3/7mXwp8zk9LeFbBFkPcg6iKCVoskPXAFyJZI96NBrTL0FuyII7kgZWW\nUNsyPeYNlrdjh7viyPZBOlXqreR0ifabN/gvMkiEyF9943ecp2sTLvcxiLUCyDfEAfyIJieHGQWy\nOkyPPxHYGJltUmhpCb3vEVk3ilXYmKQW+owdBT769Nt/RtVd0R9mHi5YIaYWzNXLpbV6zFRemF+g\np2nKJEKgL7XgT0uMiAnNvzkzUQBaYeU7LOQaYEaBrIjsMzFvg0zGibhMFO6Hab8tyhHyTIDe9+hx\nV1zI2FHm+4LunznOiZeWhe2YNK0Q+Ktv/M68aYl0pHNysGjqrUQ7JnJGUCiQVaB9bapYJkhRoFDt\nlhSfkoCtCqUnkq0kJIn9Iy7RXVqhKYfv92KSQkKd0I6pqhVCmP/tb/6lGc6dUmP2PY44J+2YyNnB\nrQdyYnBmdf/dXSlHkJoaSH+0tQL68H3IHNTERCWvkMV+TSukSkJd28iHl1/une0DMU5r88eZtcBf\nXbo6Ulsa37yqqSozSoTw429/8y+1xxVj/19943c0rRBOl0IixN9KV/lklAjxzg5VgMVCOMjyRnIS\n+LIjJyPrezT6BfALFYsJy45pZH/a467Y4xItdlJBDmsbO77Ka9rCsy/etHQcny7dnZYIWYFMIqRk\neYUy9le1QpqWAFkHEfn//vD/Bb8tr9Z0UBDaMZEVQ6FATkA99hfByeugIKJKhHSAVttYW7pc3dIx\n6soiezWgiYlSKzS6Moxk+zOtkDU9zlLi12TH5De1xLE8aAWQHsjUAP4xPd50kVLbaygtE4R2TGTF\nUCiQ45F4CaCInvQL6GmA6sFhO6Z0f9qpYz768Psi8uj5L85402p/9P2Tr3xRRF578/tYx3jKD0FX\nAnZQGKy7bEpLZDPjagNRJEI6BmQmSunQdJFmYyTtmMi5QKFAjsGRcwNl/eODAjCXtlaPMTVIhMij\n578ognIDqFHC4aCQXMKDRIi89ub3paZjopIwuyQGHRSanB/952JAyWFsVRgRE6mS6LhIoBVMBwWx\nxAS1AjkaFApkWebPDWinZwP2JQtoQ6FjaRGZpkwipDykFmq5gco+xSSSpCXMb+dSSISUkFoASgJb\nNeNwXp6ejTH3KYCt5GDnYZQFAnVD9SAQE4NtDlLTCmn/Au2YyFlAoUCWwuuuKLpuUPoeH1oh4PJW\nBWLj0vvBH334/UfPf9HQCu76x2paAp8LhIKIvPbm97GSMG0M+lyinUWI2FZyxF1RHLkBcYuJTEnM\n6KBQ1Y59dkyxh4If42RRKBTI/ITPvpd//NJbn3zb7lYQ7Il0sIlQ2Bio4TxkFMzvuP6l5TCoV7VC\nCPNBSZghE0uNDq0QcgZ/8pUvmlpBem/jZB4cL1nodlc0WxWaNhpmd4nGp9OOiawZCgUyM9M0vfzj\nl+KPqlZo8ESqbyJIqRVqBQquvIK1dDWoZ1qh/BF0HnrSEn6tECVCeqTUClEimHeQct7kSdMN4IGb\nq7/xxhvgwrLZynCOWxU8vRLjzo/Vi/Q4QQntmMgqoVAgsxETCdnxXCtU6w1ht4Id7EF6oFsr6BIh\nEoJ9TCSUv9W6JNLTy7PEkZYIWqGUCJFMK2T1jKaNgcclGtQ0aC2XYPUmSwZQVYBbFbRznSYKTgcF\nTceYWsFT3gh+y89zsgQUCmQGNIkQudcKMJz3eSI57ZjatMLh0mY5gpkYqHZJxAEdaQlxFDbKXito\new2mjQFwiTYTA6DlUlu99cs61gpgsyA7t1qgoNkxjeuYPq0QcwZm/SM/0snsUCiQIVKJ8NYn39a0\nwluffFvE564oqCei247JzvpaS1erDmO3gqYVTCUhjWkJKfoe5ypH0GIqyASYjRKeDY7uMkacGxj3\ndsykw4w6pkkrlNsKVa2QKgl+qpN54UuKdFLNIgRBkNcoLGNj4G+59JsT4KUfvJWgbiiPeLokPPUK\noO8RHAQmCoOtClLTCumcnv2LWRojA622iX47pnl1jMdJM5YjaAIO+C9RK5B54euJtBFLuPFGw33L\nw+w2Bo0tl82lZw4dgyO6lhswD5ppCbEcFKIsKEeWqYUQpE1jx4DHXVGUwA88mpyrY9LcwBLejifU\nMZ7EgyYmqBXIjPDFRLw4dxkCb33y7SEbA3zQEdGHOtn0rYfx3EBHWiJ+Ozd3GcQSEzE3YBo7Zpju\nis7T+1YH+G85AU6f93bS2fwjaQk8M56TWoHMBV9JxEbbZahqhWoiQZtzITsmvLRJpmNmzA1op2cD\nYloCVCdEsl2GQTumDndF/1fk8uL9q2sTms2HHiWReTvOHtGXkAjRRMH8k/FDngxCoUAMMl+ElEwr\nOCVCOvNydkyDgD6FgJkb0PoeQ42FWd7oMVHQdhkG7Zj8JQsdqX5Tx/gLSsw0gN/GwOOJtBBNsRzY\nODK1QJbj46e+ALJeUoPFqlaIv2qVCIEdcGKeJhH5h+9991e/9GVwvkA7phEeSiYVQJ+C1PyU0tgP\ntEKcE3y475S7OQRKKZDpCU0rxGE4rozU1YdzsY7Bk1clQiAcnJJy13TYb3/zL6shORMT+Prnjbjp\n0p6KBFHsmMJB8wVDrUC6oVAgFbK9BqAVxCpHwFS0wl4ihJ/+4Xvf/dWx1sduPOE8lhmWRK1Q3UQo\nJ89kh6YVokSIRQklUQpodkxVMZEqCS2uxNXNcI6DFtYxHpGkEQSBqSSkpjnCueCBgwFxWF9GxFw6\n1i1qYKcmagUyAl865ABgnVTdaIi1jYPF4TExECVCyq9+6ct9fY/je7Q4nIujWwHHFZyZKL0dTZ/m\n+CtpsWOqJhjKNH62OtYx5u4JtpXEp49YNeObPmSbFGVQr+5ixGEem4SmpbM2SGzHBNYV7kGQXvi6\nIQ+AcoRAutFQdVAYKT6vSoTIg1ZokQimsaPz2sxwjmsbTa0AVk+9HavhvLu20RQT4KbV6YD0sZhK\nIo7ES+PTxdF5WI3HYdhHH3306NEj854OYJujjOged6ZUTGCtADYaMq3glwiyT0vwM5+0wq0HIpJ8\n83ZqhfJ4ODi4e421Qo9E2M8+klowyw/TioRST4RNBG11T20jqAFMNxHKwK99a/eLCVx+GLcJqoG/\nuomQjjR7NMw9CNB5WOYV/uobv/PRRx+Ff5taAeck4k5BVUyU+wjZMKwVJmjSHII9tmMCE/7gnz9k\nXoG0wlfMtVPGsOruQ7bRgK2aB5vaM7kQ6hnDwV/90pftmq9lChpwOI/FCn03hKw2WWRpCbzLIMN2\nTFrdwCxpCVGyCB4dM973GBMJ2W9LrZDuSph3khR9/0J8aQktKYK1QpNEkFpmgnkF0gSFwpVi9oPF\n1ALYaFiuZCHKAlF0g1EYv5hWkDEHBTOqmXZMI+FcHGLCqST8q8+lY0otGwL/o0eP7Ge1kAiRqBWq\ngd88iDcRynPLMVpSpK8TEp+eTcUPf+KEQuHqSJsePR1ZTqvmUkmMljcqhY0BnFqoa4V9mPcEbE9V\nQatVs//zHauNJcK5dno2oCMtEb8ij+uY+AU9C/yPHj0SSztqWiGcK1ZuoElMtKYl8JVHYyU8Mj1F\nEjsmTSvEHgqGAGJCoXBdZCUIWjhPxURfyQKY3H+pZnljg1YomhfAAzfFhCecp2KiSSLgtITAbQhp\nD+fSHtGdaQlQndCxdBiMcwP2nyw5PUiEmJYww7k4xER3WsK88u4sQhAEZf1jq/4g1wyFwrWgpQey\ncF4Ow1bNy1ksiK4V4paESysoQb3sPJSi73GwW0FaHPc8LZepf8OMuYEyogM7Jo9LtNkY2bR0tnq3\nVpB9aiGVCOm5ooRzj5gQR/LATEsMborhjYZY/9i9nUGuFgqFy8ezgxCTuq0OCtUB2cgZtUJZtQC0\ngud7/0M222FjUM5vhnNnj0ZVx5j+DVo438luEvvbucAKFWDY4HGJ7v6KXC4txY5GtSxR9r0MM5Ys\nlGIC+DfgXgZx1z+Ciwfg+keB2xDCkgUCoVC4ZNLYb7YqdDsopAOqc9p3krS2aaOXs5ZgqGa5Y0Uk\nNkQavO2CGc5HdIzp35D7FkiSI1G0Qiomuj2RsEu0jN082mPHlG0iZD8OagWptUjE32rlCPGI5t9g\nKglZRivERAK1AumDQuEyMVscyyO4HGFQTFRTC5mOwQHVWa8AGiy1iO652aOpFfzuinIoEUxjR0+9\nQoz9+bVJ5dt5Jia6PZGwsaPMoRXAgKAVtKDepxViUaQn6wDSA2Xfo6kk5FBMzKUVqgUK1W5Joc0z\n0eHL4qKITY/drQr44KCYiGoA6JiybiAmBnATRPgHHqNFdM+do7sTA+npZoNldWlbxxQS4eHyZAJK\nQmpawemJ5HSJXkgrhNXH6xW08WajhL/vsZpgMNMSg1qhWo4QKL0dswfOoEAy+Jq4EJy7DAEz3jc5\nKIhbTIRh5enZmFg3AJyXsoMxDJtiYuTO0VLTMf5w7lm9ueVSdhKCvaIVJpkEKglJtILTE6m15XKw\nRq96PbHAEIdzf21jObLUCq2NEn1dEnHA0vUKWhkjtQLJ4Avi7NG+nXe3KogvN4DFRHdaIgzA8T61\nYwLVCekp4i5ZkN76R09awhQTDS2Xh7G/FATpEaAk4khZ0o5psKA1bZQAwbs84skreAwYWhsl5rJj\nWkgrmJ0O1Aokha+GMwa3M/S1KjiVRHXC6upOJZHOGUoXzdyAWf8oMAnRl+o3B3jSEiBX4TGKxiE/\nVCmWw8yDuFFiFjumvk+buKGGNxpw3+OoHed83o6ZVfNCTRDxyoGtAq5e9Ps3kIuHQuEscRoidbQq\nzOWgIG4lkQ7O+h5LrWBKhIjH2xGk+kHf4866jVNzbqBIilS7JePS9j6CpSQE6oZF7ZiaPnCyagkQ\nrQVuQ4gjtYAvw+ySEPdGQzUHNq+tQjZ/6q3UaqLA1AIRCoVzpHRXXK7vcSQtsbPuRZkurfU9Oi0T\nqni8HcE2RJlyyIvdlHDudFcU/QO67LrMy82UskRxiAmzaiGmFlrtmJz+DX6dV0qWasBOTRSW0woC\ndz08dkw4HqdpiXk/lj12TNQKBMBXwDkB3BXrHQThLKvaoDpnNnlTWiKLaljHmLE8mih0vFY93o5Y\nK4CIju0UPe6Kzq7LerlZWY5giYkoEXDWoTphtnq3f4PTjgm0WmjlCPHIoB3TSMnCSN/joql+s7zx\nv/zv/xmfzkhxzfDPfx543BXTmoDsjxr2eEdaFfxpiWpUw0uPlCOYeLwd8UYDnhyYKHjcFXHXpfE9\nL5YjKDomt0zQf0yPS6Mdk9O/wW/HZG5wAAcFqfkveeyYYs6g9Z4R2dLL9T0Ogu2YPBaiDBZXC//2\na+c+YIR/wy3/8A/w55yS2CwL2DHhj1dzadyq0B1QZa8VNG/HUig4HRTi6oPuitiOafCBa7WNUpMO\n2TBTKzj9GzrsmIBWMB0UBNY2Ss0yQYq0RFXvim+DY7m+x0Gqdkz+i6FWuFr4h18vqUR4OCgitXAe\nlQQWClKcnjFox4Sjmrm02aoA7JjMcG4WNorbQUGKcD6LSzTouiwf+EPfo5nN7uqSeBjQpRWqdZEZ\nrVohK6vstmGWvRQA6QHsv2T6NyzX9zhIasfUt5fHkHGFfPzUF0AqVCVC4F4QpKn+w99OybCHCdNz\nrdW15EGWWih1w8s/fqn8mps5QeGlwQ0dwsHwOVUG/n/43ner37D93o4gooeD8SOyDOfVnEFTyyXI\nTKSrl4F/JzvwwJ0mCuC3YX4tNlQvKZUI5q0mzcm1skotzGcFhkArYKnx0Ucfpc95NrKaV2iqbTwV\nuzGf5sHTyZnCP/m6iBIB5wbEkTxIZcGMJQviqFoAt6M07wThaVXAuYE0nJvejk4vJlnSXXEuq+b0\ngc9ux+Tf45jRjikmogY9kYAd02Baos+O6dw/dakVrg3+vddCmUXQpECM/Z6NhpGShe6+R+308rfa\nzaJAq0JY3ax/NL0dnV5Msg/SHjEhy9kx+ayaxRQTA3ZM3a0K0m7HlDkomP+mEDcAACAASURBVDJF\nxuyY+tIS0UHBb8d0MZ+3tGO6KigUTs/D16byV4V0EIeYcCqJONLvieTpe2yyYwJfJcvv/dnqmh2T\naewojrREGs5LMWH2PS5q1VwP5/t0wtJ2THiPY9yOSXNQKJfOLmPp3MBIWgJf/PnC1MKVwD/zKQES\n4WGMXnmQDnD+mB6X/czO2kaz77HJjimEc/zRjx0USjFhGjuKW0zgtESmFaoFCjPYMenJg4NwXkTu\npe2Y8B6HZscUmiDM8sZuY8fBkgWnmOhOS1zqJy21wjXAv/FpyCQC+Opv7iCIJSY8aQncdWn2PTbZ\nMTnNFmfpezQdFLrTEsABOoD9l1x2TL11A0eyY9J1TGaZAH5Mj4tDTHztjf8fXJg4ShawHdOgSzRe\n+lI/bKkVLh7+gY+NlkWoxm85SCnrc9YmzAaYaQlQYAgwnaDEsmoutYK/71Gs4kdQjpAOcP6YnSiO\nvsdROyYczk9rx9RrolBKBynsmFSXaLf9MygswHZMgy7ReOlL/bylVrhs+Nc9Hs0bDfXN4Hy8zJeW\nAAWGlfGHfY9mOcJgq4K5SaHNjJcWS0yUWiGVCLOUIxh2TDicn9aOqddEAdy0Oh2Q/tgkETxf/YEd\n06BL9EL3m1g51AoXDP+0xyBzV7RzA8Y3xWRk9XTliGdpvHq8Bpm17zGEXvGJCXBJoO9xvEtClCzC\nSN+j144Jh3PsoLC0HROsrMTuip56BWcioequiLsVupVEObkzLRGXvtRPXWqFS4V/12Ux0wD5QffX\npr6NBvPg5PgaOtL3KAPh3BMtcN8jdok2dcyMfY/HdlDosmMyCxvT07WFBm9a3SERIqYnUreSkMa0\nRLn0pX7wUitcJPyjLoUznEuXRPD0PTYsXfwWX4+pFczTO8K5//nx5wacdkwPnsGnsGMy+x69DgpL\n2jGJ/gcC7orxekbsmJzlh+u0Y+r7+F25hwG1wuXBv+gieMoRpEUiSFHP6Ol7HK1/bNQKqYPCeAWi\n9EqEWTYaqktXjZazcD6jHRNIA9wf6QjnXRsNI3ZMTnfFDjumJhuD1doxNX0Chwf+4osvisi77767\n2k/vlUsZ0gqFwsxk7xCtW+EgnPuSq+UgT99j9SxxKAlp1ArYnSk9S9w2Bt1ZhJG+R3HUQ+BwPm7H\n5G1n6Ajnvn0E/FunHZO/WiIOdtoxzRjOtdOzAYumJfyv8yARIu+++67z9JPA1MLFwD/kbGgiuswE\nSBnOce06WNTR9wikgyst0WuiUDo5itvGYHyjoa/v0bk6/sKUaYUmOybnAx8J54N9j+Wc+fy1O0WJ\nTyuY5Y0zhnNZQEwsYcdUlQgpK08trPbaiB/+FWfAzLMFrWCEc70ZcrzvEazuSksM1Cvsajd7jJTx\nu0kimGmJVjsmp0Tw2Eq22jG1pqC7w3lf32OTHRNauksrhN2HQRsDcYgJre8xvIwXtWMS5UMgbjQA\nrfB3f/d34MJOC7XCBcA/4RDOrTgzMSBJPHY2SrT2PY6mJbq0gumgILAmoLqWONIS8Vd49aalpSit\ncO5QaKS1jR1vQ+Su6PBpbup7bLNjgqu3aoUoEeKRcXdFUcRE1UEhy4otZ8ckSWpB22gotcLKdx8C\n1ArnDv9+nXRIBE84l96I7up7hEsP1itIcZcHSaTPeDmCPy0hhxLBVBIj2xwjWqFbIqTzA/8ls1HC\nkxhos2NyrO68aXV0UJjRxsATztMdBK3IZjk7pjhSSx5kWmHNmw6B8Mb5o7vvff3mSyu/VAKgUGhm\nJIswEs6107MBO33kDGkJh8WC00FBinBudlSCpUf6Hk1MoWNKKNAogdc1h03QqrmM1q2NEt60hBJT\nUY+G9bhmtDFwhnNRtgCq1za7HVO6x4ErEl588cVzSST80d334o/UCucLhUID/p4frdnh/reN4Vzc\nEf0YaQkzbjnifYdfgictUc6ZjWnaaHCu7pkQN0pUx4vIN//Hn4rIN37rq6YQMbskRAn8I40SztyA\n0aOB8zHLVBemA6oTeqoLZ9cxWT4GJBXkHCSCiKQqIUCtcKZQKLjoaAs2tYKcKDdgDnAt3fuy8Xgi\n4R6Qvl0GacwieFZvyoiEfzRJhMg3fuur5bnO5sPE0sPIDYhbTPhzA9X2h3zpxpCchlUgJgbdFUF1\n4exdElo+5ryaISNZIiGDWuEc+fipL2DtdDuH7Eu6igmT36LToZhIBUE5bKdLgcka4F26tzrJfOD/\n8L3vapNnFQ9lWgIbPg5+PKWrV5euvlSy3ABeYpqmTCIEwsH4tJTRdye7cvVMSRhL6xsNaWqhHFZd\nulwdaQX9tZSpgTKsVp0V/E2PYEA4mD3n6UhPM+TIDaLSPYhzkQiBcOVAK3R/qJJTwYyCyiyv5tJE\noSj4qp0lIlY4f5jQ96tqEkLVMc6lB76dAxsDT/OhtFg1z/sid5YsaLkBM1df1QrxdHF8QddsDPDB\nGeyYFB1TXb21ZEEc+xSnclfssGPyl9A6B6+Ncvfh6zdfkuSxsBXijOCfqsK8gvfARKFebpb8GIYd\nHukWE6AcoTo5/rG+9ExaoVpvOFjeKIt9wnoKJrR436cV4u4D7nsUTziHuwyL2zH1liyYmDUN3eHc\nU1lZzlmu3ufDCFob1v/VPKQWMomQ/nbNF08i3Ho4YIk33v2Hu55cnaAdUykFspEeraBeGxQTrqXH\n9iCAIRLYg0hPL3+VpiU6LsyDtnqUCGCX4Zv/40/xkxYmT7VCqi2AVvCEc2yfUJ6bDdYmT5fWHhe+\ntnv5MhA5dtATCRc2pqkFbYPDbFUA1xambZUIgRdffFHbVwrDzCfthPHY9elHrbB6KBTuWVSb42m1\naoaHAVZuAEf0kf0L19IDb3VcftiqFTLNgU8fpLp6jO5BK4BNBKdWqJYxllphxnAez9X0RFVMmEuL\nX8csoBVSGwNwOgjn5cyZnjDFxEcffWT2BInioJAKAk1JaBWReMARcH36USusGwqF+/fSD944qBc7\nMqZWkIHcgBbm/WLCXrrxeUuLDMxbNHm0Ql9awnORZjiPq2eyoKoVQuDPihPB5FYiav5wXg6urA7F\nxAxpiVm1QlqCACyPzH2BOLM2UhMT6erVhxYrJbGOCfFeUxLZzJmeOK1W0Ih1DCu8NpJy1X+eVCJE\nvvBK5TvcMa9nJJyXp2cDdvrIanmEf2nnM1atQzRdklp7JjOct5Cozukpq8Tlh2nLg6fvsek6Z3RX\nlJbaxpFyhOrqlR+H6xW6bQzwzGae3/RvyCJ6djFVKRAbHzStIPtqBpCZWI+TY7XUcSXXRkquVChU\nJULKF14xKs4GV8expymcizuix3SzKSaal3Z/yPpbFcwWR39aYvym1eXpqUQAuwyxWAGLiYW0gjOc\nayPNqkOw+jwu0WOmFzNWF7aubhY/xohe1TGZT7MUr3wgJoCro6xAK5QSIUKtsFqubushSoSQOdAA\nGmJ89Z/85CdmVlmUcJ6pAedGQ6YkADtdatSXbnljgwrEALiNUziYbtP6HRTSc7WlQWYiO73MIsRi\ngqwCUZKAqq07SPJq6QnnAiN6llooGxOqq8/Vctm3B5FG3xEbgxFwWWXArLsEPs1lXiEdCVIOcelT\nxeMJ+itwD2K1XNFfpcwifOGVr1YFwRK7D1EixCPPPfectlsJovXDSBEZyA1UT4yDXXscY9/2tDsv\nYCUhjuRBq0t0a2bC3GgARorzbj1k82vhXB5eVz3hPFY+Wmkw2PfYtcfRlFSofkEftzGYN6vR5N9g\n7nFoYqK6+5DmG06YV8BaITRSyok2f4nGVQiF+H2ulAWZVjiORIhkWiGrZxwM52CG8lfVJER57sHS\nYx1ZaRlg+TEH+h5H0hJxQPpFvxQTuO8RCAVx3JohLVlY4pW2YDg3/+Lz7XE02SrMvtGQpSXGqyWk\n0b/BE849f+7U2HE9htDV3QfaMa2ZC/9jAIkQCVphoRrGaZqqEiEStIK219AXzkVcYgKkHLTV8x/1\njqxBMTHV+h4zUaWF85/+9Kef+MQnPHZM4D5PWt+j6BUJztg/+D11kMFwDv/is9oxte84eMoCPPOI\n774P/gsL/3DqGE84TxMGph1TOWc2/8lTC7RjWj8XW6OQSQRtlyGwXOliEAFYK2SJhIPTrXBeSoFs\npEcrqBcPV98V7+RM7pQDspHbv39bG5B1HmrbtFk4/+lPfxr+DbRCmpbQHnhY19/32CQRsB3TcnjC\nuW3HVP+L79LTR+2Yut6JoCzAGeZBWiK774OTLC0BRsaihPJXmaVSuadg2jEFqaFxwrbJhy9IuAeY\nWmEFXOCfQcsihJxBudEQlcSibQ5S7D4899xz8eBzzz0HYrn4thiacgNxpKscAS8Nqy+n2kfY9u/f\njkc2v/YS+Erk37+IKiGSaQVnWiId3FeOoF2qZra4NNUCBRzOzXs6VCsoq5PjH/PVx56TCXoigbME\nWiykE46kJbAdk6YVAmZHA+iNfPfdd3GXxMlbIQDUCmvg0v4GZcViRrrRMK+DgmfLMMoCUXTDQuFc\noJgwD+KchGvp/UdYKhEiQCtgwpylRIgErQDaI0utkFVW4r5Hf9io3r5huRfbwbD5wrk4dhDi6U7/\nhlkkwsNs0BNJGx9/9IRzs24AlCgCOyZRTBTScG62PgKLhbJLIjuytlgQSxnYNnlyLk0oiMg0TVgo\nyNwOCmnCoNrLkA4TpbAxEFILs4dz7fRswE4fWS2PaFq6KhEirVoha1XQtMInPvEJcZRKgFtNlmG+\nqQ4RSA3p0grp9YDTgUR4GLMSO6YFCoO6v/eLL5xr/UpmWiItR9AUSRa/wY/pcdnLcawVQGZiPVqh\nrHakVjgtF1ijEL+8gq0HULLwgzeMe/aklHsKwSNBavn2NJ2A5lSCbhaSy8jtURjdYiI9vXXpIBE2\nv/YS1gpNn++ZiUIQBKlcSI/gmfFeaZpa6GhV2BV3eCofTtOLLfNv0F5s9y7RvzVmY6DXNnq6Jzxb\nGODauhkpR5Ak2GvhvHRiTtMSTq2g/TYzUcgieplXyKodPVpBW3oNNs/VhgihxcKpucynvqxklJpu\nqHZLishPfvITkBhIlzByA/sIBKoTsoMrzw14ri1dOlMGm197CR/E9QrhH7jvMZQxSrEf4UktAAZb\nFbS0hJkYkOSBm5mJ6h7HvbjpCucPDgpWbkBqYqJ6cKEsgpPwFN09+xsRuXn8G9jGYNBdEbRO4B6N\navtDSpob0DbUtI0GrBVOm1SYoMWCMK9wOi5NKKQf6GZvZGmiYBoiSSIRnnvuObP1ESuJNMdQhvNy\n7fuI7ruDVDlnOb+mJKonxsGDGw0xtQB0g1b/uPm1l8y+R1yy0JfqD3Pi0z2ll1UPR1G2M7JdBiwU\nyjmzAT3h3L3LUM6ZjVlio6FJvaUSIXLz+DekZmNgKglxiImoFaoFE6AP04zoHguEKbFjAs5LTXMe\nAawVaMd0Ki5KKFQ/qc2ShVjbWAb10hBJCndFzUlJHGKiKhEeliu+92cXg8L52O0oy6WlvBhtaZHJ\nUZEglpjQ6h81rZD2Pc6iFaqVktXT40hzcmfJAsgNAIVhTi7+cG4pCWkXEwtIhPBvo4I4/COTCClB\nEKQSIf2VjLkrCnzg06EdkxnRm9wVwR6EJEJkJRIhUt19oB3TabmKpxtoBU0iRMAOQhyQSQdxiIlU\nSXi/91cNSQ4HS5OYcGgFkJmo6JjDH7UGBxHZ/v3bZskCGJBphWrfY7VbUtzhXPRmiuz0bC1t8hj7\nx3MDwNuxqhXinDGvoOG89/QS3o5NFSrZwCAGql/ZY24ACAVzQEwtNLkr+h+OuDcaWtMAuF4h/GOd\nIWCiHdOauJbnutQKaTPkpHsiVesJyjEh5FdHepQEiNZifRAAJREH5BMeHsFiAoDFxFT4JUhxpFtM\nBK2A+x6z+O0J59rpGeF0TU+USkKKAgUtnJvFj2K1XZSGDbl/Q+2v6jdRkLGSBS0fs9sZuQG5jw3g\ntwd7jmVuoCoFQs7g7tnfmFpBFnCJ9tQNmAM8Ts8jaYnTFpSMi0syzhU90VErVP0SqtsKkrQqmDsI\nZjkCGFZqhRB9nS2X9lu9MTcgiYQyyxHwgBDvRU8wjIgJT7Cv1jbKmFbQ5swmr+4gBCrx+1BbAJfo\n+9wA/ADViiGk0ArVAoVuO6bqhNlv86RXMlDTCodiQl96EtF3GTIpECWCNiAbhksWpNGOqcMTaaQI\nEZcsgOecdkxErkooiE+iahHdPAjEhDMtES4rlQjZDOX31HiK+ZWiITdQk1BVHbNLztWmFUe9QreY\nMN0XBusVytPT7glTSZiJAdByWSYGpEhLYPlo1jaCiD7orjiZd5R4aEktL772Os/FROUs2SsJnBgI\nIV8UPVGKiexHYMcUZm76cp/i8UTCdkzdJQvZudUkxGojBbXCEbiup9j5/dvfqmCKida0RDlnNman\n+CWbMd7TqmDucYjUswjVHQ1JAiquW5QBMWGmFsxw7tQK1cHl7kMqETzlCN2JAalpBWejhFnYKHs1\n0O2uCLTCfrC+9CRASUihFcofQVJBYGGjWGIi0wplCWRVK6TfzvvKEdIB4EczN4C1ApAy1ArXzLU8\nv1nb5KBWkN6NhpG0RByJ4326TWuGc3FIhPth+6orY3UoWbCDwoiYwKkFLZw7mx6l9oiyMXGPw9+q\n4HRXdN7V2tkokS89XI4AOhpKrZBKBLCJENMDgHC6NrLUCmng95Qj9HVJxAGZZvK4K4rjhg5iiQlP\nbgDvceClVxsvqBUW5fKf3GoWoU8rOJsepdFBoUNMDNYNHOQGfBLBk5YQxyaI6aDQJyYGSxbMbMF4\n32PaqmCaKDgTA/Hccs5ydbR0PZzHIhWcGDCeloPaxl32W8kOpkc85Qh4QJQFovQ9Vg/GBz6YlgBf\n0Mto3eSJVJ0zHWDmBs7RjsmEWmE5LvmZxRsNznqF+GPZvDDuoDAiJsrrjoEfVyBKl0RIz62e1jTn\noB3T0iUL1XYGc58CaAVnONfqHz1piTntmOobWMVGg6/pMfzD2mioB354MGYH1WkDZsmCHEqE7OLL\n+kdTScSRfhsDvyeS046pOy2xZjsmjCf5R/q4TKHgqUWII/tqG2UOBwUzLVFdN/42qwkodohrj0hE\neiXCwVrtc6aTj9sxgZIFsE0bMgpYK4ijnUH7bbUTUnzh3BzgSUsMigmrSAXaMVXlxUy5gTgPXqic\n0KxtDP/ADzyWLJhKQoo9DtPGYEYxkYZzULqoTQiWXnmYmBKPJs13gYxwHkLh9vZ2u91uNpvb21tz\ncGsCyqMVtN96+h670xLh+6X3VpO1+qmOcB4lgln8CJY20bTCXHZMZRofeCrEg+JQEiLycz/3c+C3\nsk8taJEb9D3uLBOFVjsmz+pN98v23OSpWoQ4khsQRSLk1wbTEjic2xsolpIQuJ3RneoP/+i774PW\nhykOMbFyO6ZIKhEi1AqzcwZCIbwUbm5u7u7uZJk/f7dWWMiOKYtzeHWzVcEfzkttAZIKzjnBWova\nMe10QySp+S+ZSkJa7iyFtyHwTaur57baMflXb1bVtlUzON3IDVTLEsuR6rUNlCyARJSz7xGnLjpS\n/X47JnxtIDdQrXYsh62QqkRI4R2kZmTtQmGz2YjIdrsNP07T9OzZs3BwRjz1CnIUOybt08qobeyt\nw8hHlsdh6wQuXfQsuqgdE04MBDUA7JjKLgn/MxlmxtEafJWfChMFU0nIoR2Ts7axLxuk2jFZhkhi\nBf4s2JtejfEsz+qmHVOZiHL2PZrGjgK3IcSxy2DaMZkfX6aYOBeJEJl4t8ljsXahsJAySOeXxrbJ\n1laFOHl3YiA93VwofVymkhBdIjwMyLokit+aF4+XHrFjKoVCVts4WI4QlUSTREhnwBsNeLZBOybD\nRGFM4eX+S3l0r54l8hDOnZsIvp0ya3W/HVNMRDn7Hk1jRzn0duyuG/DYMWk6xmnHtPJYUEXTCnH3\n4Uwf19pY9ZO43W4fP3682+2aahScdLRNlnmFjPzO0bUsH95owEubq++U+1eB63nIDWhL7/9hiInG\nDLbM56BgdkmMuCv2SYR0ErzRgOc0myywHZPnZtYdS8teK2h7DaUDklTC+VAarFoJoa2e/ThX32PV\njsl0icbuioN2TJmOmdFdcbURN9MKZYHCaq/8jFjwGdxutyHAj8zw+PFjEbm5uRGR5WoUUqqvquwL\nuu2JZFU8tEqEufY4DsrNsoVEsoPhyEPnIbjIcK774czuoFCeno3pcFf02zGZexxNqf60dWLcjgnc\narJ76TjYLEcAHQ2lVigaJYysA7w2tPrSdkymSzR2Vxzse9xB/4a+zSZZ900fglYANYzUCoPM/PRt\nNpu7u7ubm5vNZrPdbkdUguyFwpMnT2IiYZqm9MeFSF9V2hd0re/RXXXVJhGyVZb1djyUCJF7U6Ni\nvPhsDMyNhj4HhSh9zLSEWAJuxI7J3ONolQiRWeyYmiQCWLo8S/ul2fQoiRqAjRKuaobW1c3aRnwQ\niAlnWgK7K/a1KohDTLQG+3nTEsvh2The55WfBbM9d/Hbf3w1j0f0uPUQj2S1jcthpvrLboWm2O+p\nGxixY1J1DNxluF/dNEveDxPfnRdm32gwV6lOeAQ7JnCw1YAho8+OqVudlEvj07WNBk91oVgR3dqk\nqMyZVB0ZS/vtmEwx0ZqW6LYxqCYM/HZMzs/8edMSa8Bf3E0yZhMK6d8gpBNmUXBZMePRhEJcfaRb\nIZtKigIFsMcxbseUr5WuomiFuIPg+YLuufOCFryrZYmepcWRt8BpCS0zKTPZMaV7HKVuMGsOtN9K\nox1Txx4HmPwP/+/fN6fSUv3mQY+Y8GQ1tJZLgTrGY8fUt9HgSUuM2BhMh32PTXZMGE0ipFOdb7hl\naqGDeZ6y8NU/RPS46ZAlA/pIlUG6yugVW0TdM9itIDWJkM6QpZTndYl+0DHlVdXKEbKL0cK50xBJ\ndDFR9j2aS2erd+9xjNgxhdeDWZEgupjo0wr+vkcZ2+MYtWOyShbEoRuqB82SBXyRpo5ZumRBG6nZ\nMTUVOPfZMYE5Z0xLrBNqhVbmFAq73S40JgS5cHd3N8vk8f0gIscpUJAiZI7YGJiJAa1VQXpdotPB\nznIE7Xs2MC0YMUQSuIlQXVqKAgUtH+NZetdox5S9HoAdk5mWwMaOU1ffY6okzO0VWdSOaaxuAIsJ\n7S++s1ouPde20B2ktNOzAUuE84AzqJvbGemRC4iy1ApNzLn18OTJkyARbm9vnz59OmNQDxmF4yQS\ntJA5YmOw0E2r0wHa9Uz6DaI8rQoT9ETqNkQSqVcz+JcWKw9hLu23Y6rGJ61LIj29nLM6Wzl5U9+j\nM7pHibC4HVNxUmpjYJY3ii4mMjWQraVphUMxYSw9cgeparek7BNRnq5LMCCuIo12TM6U54xpiXOB\nWsHPbM9UiOK3t7e3t7d3d3dH+Op/ZDwFhlpuYKRVwe8SLcpuSNkMOdjNmIINkaQWsJ0FhubScpgb\nWMKOqa9LIs4A0hIYZ72CJiZKJ0c5sh1THr/VH9Pj4gvn8S9eVSTlAzdXj8UNnnAusBEDODlqk6d2\nTOb2SpMdk18izJWWWIJFwzm1gpPZth6Ca0IQBxcmESJmbaOMtSpIsa0gLWkJMbdpeyWC+e3ceQMn\n0PdYfr6f1o7JbFUAOxfpVE3ejtn8o3eOVrok4oDuRgkTUBMgNekgjnCuDa5NriqJ6urzhnPtVpPl\n5NW+CU/1YpWOcgSPf4PTJXoJ9pvOuxm/0FZXoVYwmeE5GjdWOiPMhIGZGyhPz8Z0uER7rJrF/ZYw\n+wzLI55wXs6ZjdFyA0e2YzJNFFKJMJeDQhXQDBn3lXBiYF47Jj8gTu8HIEMkqWkFZ6OEqSSK1ecM\n556yStOOqU8rjDc19LlEL8Q05fXW1AonhE9QM1n9gcyaGzAHeNISg5Gpb6NhJJxrp5dLexbqW1qs\nvket/tHv0VTiaeyOaqAazqvbCuJQEnGk5zr9hOt54403ROSVV14ZD+fa4GoeQhxKIhk5ZzgPD9y8\nb3U6TKNVK2QOCt37F9LuEr0ESSIh/w21wqngs9NG2oIhA7kB0PeIezLFt8cxohWc+wj4tycUE61L\n+3MDfgcFp0TweDGJQ09M+h2kQFoiCeed34/LwWHOyCuvvCL1cL6vtLXCufRmJkZaLpvCeXVwufuQ\nKQmzS8KjFYCDQvVl4Ol7LCcsxywasHGr1sJLX1S15oxQKHjxtE1m4G0IfHeo6uROBwUR+U//6T+B\n32I0reD0RBrPDQz6N/iX7nBXHN9oAGmJ4RJCQ0nIoUSI7MN5HtVMJVEOrhJSC9Xv8dVNivbcABrg\nSEvU6x+nOeyYdsodpDx2TIO5gVTHNPU9msUQi0YNqBUWD1hMLVThk+LCbJvM8ER009hxgv5LZloC\nPyJM2WooxefpMe2Y/P4NTh0z4q6o9T2GukUzpnpcoptI+x7N0GKF8x1MDKDJwcxiqY0pd1fMfgsO\nejMTffWP1Z7JI9sxgbSHedMH6e17BPLiCCGjtvvg/bo/nhigVijhMzJKGc7F54nkiehTcqtJpx3T\nXH/QKbEx0HTMouG81b+hSccMuisCywQt2MdEgpmWaOqilKJAoVqdJ3slYYZzU0mAi6mea4qM9CJb\ncgN5YCvPdbZcppNrvz2tHRPeIum76YM4An+2x7FAuauZqdq1SoS9vBgKbdQKGXw6ZgCHc4EGi7Pc\nHWpeiZCubr6TZw/n6eBuY8dy8vJiRtwV077HaoJB65JITwdLO8Oq1veYrZ4pCRzOcW6gSSs4JUJ6\nuhnOtQmnwrxBfC2Xfv8GWasdk6dkoTs3MPvOfVM4b2rUKjYshq6cWiGFz8U8OOsV/K0KWW/F4P0m\nFqKpS8ITzqWx73FQx3S4K2qDywFal0R1cr87E5AIEWyikGmFaoFCqRXCsDfeeMOjFVolQnq6I5yr\nz4xW2Bgn6RATUZ14wrms0o5JrNtRHgctnIMMnGmioLdIGPN7oFaIQ+9slQAAIABJREFU8ImYjbla\nFapZBLPl8sh/x1QizBLORckimH2Pzsm1y2tyVxR3baOpJKQxLSH75/y//83/FJE/+I3fHLRjwuG8\nzA1kP1oFB4NZ3/RHkUo4V0OLZ4uhSUz4w7knD3QqOyZZo0R4+H2ZgTNNFFIlAYVC/XQn1AoBPgsz\nELNzIMZHEbCoHdPR6oxmDOf7N/xQ3yOef/ymD+LoiSiPjLdcptoxSITIH/zGb8qYHZOpsaKYqCYY\nRj588eq+3EAeWmIFw0hho+RiwhXOnU+F54EvYce0hts0wF4GCZEI6olSTBwoCe2s/YTUCkPwKRgl\nexllWwYC9xQCTiURBy/koGDStNHgCefOtIT4ShZMLQJ0jBnOxdEo4VQS8WD8i+O0hBQSISWmFprs\nmDx9j84uCe23YM4koqPtFUc4R60K4KBPTLjCeZNEcPY9zmjHtAaJEIFawdxBEEtMlNIhG0at0M/H\nT30B5035AgovZfD9PhUE5bCd7rbkabn8yU9+stxrOnmX1gGNkeFgvLYycnuMkoCY2P792+CBx0Wr\nq2z//m2PgwJ44KInD7LUQumXUC1lyJQEXhqkB9LUQumg8MYbb1S/46YSAfdB+F9sZUTf7SrfsP25\ngf3qlePZ6aVu0NIVh2ICvZxaH3ga+++e/Y05OdYK5unrvNnj/qWe1yjIwUvdyHLpA+IehKYkGv5q\nypWv6/k8JtculJbDfFGCAVNXy2WqJBbNK6zZjmm8ZEHcuwzSvtGAvwtqd5ACWw/+cOXJDVSzCOXu\nQ3rEk5YQmB7ozQ3cX7A2c9RX1tKuhVqJD1wL+Z7ixBE7ptV+sCfBXhWpVTGRnDWalmi53sOTV/zE\nLsqVPuzjMPKqmtpbLrMfj6YV1mbHNKgV5KR2TPjTP1YySmPSO1QbYK0g1kaDVrJQLYrM0hKesgDz\nt5qDghTVjvGIZ5dBHDoGXX0+/mCXAQgFcZSJdNsxrRlPrUYqJqBuKI8YxRB4aZPr1ArX+JiPSfer\nqloXmQFqG+UoWmEhOyZTTJh2TB1aIUoEnBiQ5e2YzC+arRIhHqlqhbTvcRY7Ji0tUb3qtPlwsG7A\ns6FQPejY43A952YaoDyCEwPanNn8F/kZbjVKCNxosNMSVn1Mv6C/SK7uAR+Zjp2trBzyVHeHMvHs\nrbTaMXnERHVwOaBJK5SeCqe1Yxr8q4E6xLLR0VQS4hYTOC2RaYUy9le1gsdBAZyezoN3GaCOMVb3\nhHOtnUGL9HFO7BINZjh3zJIFR/LASEuU70RTScTBF/mca1zXoz0mqURwvqqqWQTzTtbYjumEf9+R\nLol0QHZEHLWNrXZM1XB+WjumkUSUJzEA+h61coTqgGyYI+uA0gNZsMeuBulZ8hDOtcXrEgGslQ2u\nrt4Uzp31CqXs0M512jGdlsEdWND3qFRH5sN62zLt/ZHVPuezc0UP9WholeTmNxLTQUEcuiE9cqo/\nbioRRhID4qhtFNj3aGqF7sSA1LTCvHZMHYmoWMSwUDlCOsD5Y7m0JzGgBfVJd1BIBlTOEt8GB1i6\nuro/nON2BnCfyerkrXZMJyEJwMann69e4f6AVDIB5ekPSw+mJagVhEJhXvBGQ/VV5dxoaLVjWoNE\niGjhPKZbzOrCcs5szIgd07i7ojZ43I7Jn4iq9kSA/gUznM+YliiXdtRhqEuDLolkQD4+nR8kFfDS\n2erV2zdgOyZNK5i3hhK4fxEHrOTzvPYd3XBXdGgF9dO1JibAj+lxOy1hvR6uIoZexYM8As5ahPRV\n5dxocO4yxMGrkggp4L4Pprw4gh3TSDiXBeyYYm7AvHJsxxRjs7j7HqOAG09L4KVB/wLYRHCHc3WP\no9QKrY0SznIELaqB2saR+0zKCrQCLEI03BUHCwyhiUJ/WqLpU/1SufxH2MF2u729vc2OgPFNL5TY\nYN1t1dx90+qlwbkBcQTspl0GaclM3PzqF8yL7wjnc1k1p2mJJr+E2DNZJZwr7o0G0OKYjhd3WsJc\nOm2UAMG7POIrR8BJ4/oq1sH7Cx7ve1zoptWnLksy8vx6TYCU8chfYNjUJeFMS7R+ql9wMKUzY4Xt\ndnt3d3dzc+Mc3/T6SD8ZNYD9YvmrNUiEwE53W/JUF8ZPEODtmJUspBNqq3uWjqf73RXTkR6FAcRE\nPL0a8sNBTYzGJ60qL1LfBYC2F5ClFsrcQHX3IVMSeGmwg5AKgjJya8UEh2ICh5bK6eXSUpMd1d2H\nJmMDLXmQpRbKObWNjzXYKjiqC7GYeDBPrAX+uj1iJia6Cw6Si2+O+k1F6+fIxT6wETabjVhZhEFA\n9WIM/BO8g9TJyxE0ys5DablIT1oCTFiu3loSOGNuwG/H9N//5n/+wW/8JsgNiGMbIpUFHSUL5te1\nETsmsPWACwzFkTxw5AayJHNDWkLgq3dqvO9DGuNNOybP0rIOiZACmggc4Tx+g7Ljvb4QWtpRS9v/\nTF6qVrjMRzXINE1PnjzJdh9mnFySj0icG9C0wjolQmSy7Jg8p5fHw4Q/+9nPPvaxj5lbGCOra+Ec\nZB2ykX47pjSia1ohxv6lSxa0crNZ7JjSpcsvhVpZouyTB46NBjRAq390pCU8bRptEiEyix3TOj8H\nksyBqRvKI56mx3LOdICx9HJP2kVqhQt8SONM03Rzc3N3dxd+fPbsWcgxjE8rtY9InBuoiomPPvy+\niDx6/our/fONC/MyLfGzn/0sHgFaYfzTM1MDpUXSXHZMZVDPtEKZHqhqhSgRPGkJMTyRDh7p7HZM\n4K+mlSPEI+WpTiWhnV7+VtMN+q52XN0QcPPaMfk5YdxqqhtwKIns4DxpiSW4PK1waY9nFoJQCBmF\n29vbu7u7RZ8l81WViokgESKPnv+irPUrxSBpYiCVCJFSK4QPpjB4UElM+i2apKYkpNGOCX/vDyFf\nG5ZphayeEaclfJ5IO22vYdyOySwwtGwMvD+mx2U4LVH2WEpFx9jbBBp+OyY/4SL/23fe+ONXm28I\nPkja92iWNzqUhCyQlpAlwl/8bLkwrXBRD2Yhgm5YumTB/EaSSYSUNacWukmjvkZUA9XBH/vYx6S2\nP20qCXGICdAlEdDsmOIFa1oBSIR0jFbGKEpaoskTaVE7JkvAaSuHMWoZo9Skg8yXlgD7F+kA/XTk\noBDC+VxaIUqEeORoWqGWRdACtlNMDKYljKXnKlkoPzEuSStcziNZjiPUNoryqkolwqPnv6hphQvL\nK2ThHGiFIAUE6okgCJxKolzd1ArggWS1jdWFtI0GvIlgiok0LeEJ53JcO6bqU+G51aSvHGH+tIRn\naWnUClEixCOaVvALhVIipJMvnRkVET0wa9/7qwOyYSNpCXtpsK/kFBOibxxrvzovKBRygolCKguO\nkFGIC8U/RzWLEARBejA9gvMKZ/GSBeEcxPhpmkwxYSoJvDr4EuyxYzK7FUCrAu57xGkJZ/lhkx2T\nqSS0CbMB6XOu9EQcnJLZGAyHc29aIs750UcficijR4/MrANc+qFLQgvnfS2XQCKkky/0CWB2RcLv\n/ekA54/pcbHEhL10mXf0iwlzo+ECUgtn/wCWIO162G63jx8/nque0Vw3/htvNHz04fdL0SBKaiFL\nS6z2L27G+/gtX2AaID1FHJkJf1pCu2xZ2F3R7HsEaQkznMuSfY/DdkwihxIhUq0SiEfm6JJ4GBkk\nQuTRo0fl6dWLVCafzHDe0XIZ5vzjV18Bk//xq/a9NrqxtIK5gyBWRO9LS3iX3iEnKKmKCf9Gw7lr\nhfO++oW4vb19+vRp/HG5VslI9rKbpgkIBYHbEJKoAS0tsdo/uic3gN+NqZgwMxNOJREnx79tdVdM\nJcKMGw0z3vRB3CULcigRImXCoHXpcs5sfi034Ot71Ka9/0cmEVJiasEvEdLLniWcV7MI4XRwEOQV\nBvOOSqB9mHNYTHSkJeLIGXRMUqRS2SZe7efqOJf82AYJew1LJxLAy64qBWIiwdQKoqclzk4rmBIh\nPV2s3EBwYiiHaVrBubonN5AKuDKiV8sSZZ/bNMsbZ7zpg7hzA6JIhJRwet/SGLP+0ex7zA5mSgIL\nhfCPwZaE9GBTOC9Pz6aK82gLlV+RzaU9JGqg8uEGxIQjog+mJUaWlurDuRIoFNZLphWqBQqlFIi7\nDzgtcUZaAVcUpmfJ3o7J7JUA2xBlysEvEZw2BqCwAFgmaFohtWMy0xIjuQGt79H51+leWptQ3HUD\n0lv/CLYekqWNogQsoUbCuZmWKOfMxsSvyObSTXiahlIx4YjoM+YGWpeWLCky+ydn692Fjg/v9bBe\nYrDH6YFIVqBgaoXV5srClcfv/c44HQM80AoxkaANiL9qymHE8AyidRrOtani6dWtijKvkMkObXXP\n0pLcuKH6qMPqzkbHlLj7ADYa8NLVCdPg/dFHH01TRSvEMB+//VfZ6Y2R+9TCoygLakujSjcwQPbP\nqijhPBwMp5fDQkdleW4a4/F+2X/7zhvaJki6NJhBwzzLKghIBUE5zJMw6BYT2em54gG7D9203l3o\n+Kw0VJCImRiISkL7tAKtEyCvcPIuCc+HFNhoqOYGsiMjGw1x9T53Rdyt0Ncloa2e/Wi6Kzo3EZx/\nnbnsmKQmEVKyvEK2caBlHWKFAeikEKmnFg4nz7/3F6vbeQUNM3mQpiXKVJDWXiH7Bw4m/9Z//Qq4\nsD7SzkPfPoL529nTEnLMjYbjdOCPQKFwBgCt4HFQiKc7GyXOokvCrEUQWI6QDkh/lBY7pkF3xabb\nSfuXNlfPtEK1DXJQPlabHuPkfXZMcWazwFAL6qWSkErherni/WyPHj2ylt4BKdOnFWLOYKSdIZu8\nuqNRTr5Ei0S18xC0Kuzs21GO1D8aS8ux0q7TkncXmgUKhTMg5g9HHBTEarl8+Jhbt0t0KhE8dkzd\nYsJjxzTorpj2PTZJhPEuCdNdEfQvmEpiXjumcgnwzV4crQqgT2EqbjmRJQa6l5ZGrVBuK1S1gjOc\nh8nB4HL15SXCwy/L7/3FV5fyxFnSEsbSchStMC1zd6EZoVBYNekXuBjFux0UurskZMnUgn+Po5pF\nMDshB8WEx45pxF1RerMIR7BjSh+42ZiQKYlZ7Jg0Sq2QFRjicO4oP0RpCcVWYVe9sKbV03CuJRj6\nwrnnXWaKiW6srICk3/v1gg+Qh0C5AUtMGEsfgem4dxfqgEJhpVTf2NVv/Cmmg4J28OFj7ohaIft2\nbn6K+e2YZhQTTjumEXfF8T0Oz0LVpcftmHbKHaTMZkjTJdqpFapf5cuDHokQZy4nzCZPl66/SXsb\nJZz1Cn1eC1hYLFqW5NlBcFSoaJsIUs0NtO5xrCQaTsfyAvazlqeG+DHDuTiMHUVPQmgtlyGZMcsL\nRvt2PqgVyjmzAcvZMYFtCPGF82q8ceYGyjnL1VvtmKJ2nNGOqdUl2rWhZkX0PolgpiXwnB4d4ywp\nyJjFjqn6YpvFRMG8HtB52HV6NhVQElIVE2VW7OQxcYW1jRQKZwkI52ZXpCgSQZs827MY1Aq4CLFP\nK8RiAqwkxC0mOuyYwlmDdQNpNsiZG4inzGjHBKoTIrN7O/qdl5zhXIazCFpawpkYCP8w6x/B6X12\nTNrp2VTpi21eEwVMklrokwh4E0GgkhBZwUZDxgnvLuSHQuFcKcN5Jo1NY0fj+5Du3zCuFfruqlA9\nPRuvTW46KKQjm+yYsiOD7orizg2U4XPcjslsjBQ9NwDsmHBaQjs9/e3BhpoVzmWBjQb/Kn4dI+5t\nCKcdU/XcKn47JnDx3bSWJTn7HpMXib3HsZ6NBpGDrodj3l3Iz4qeLNJKGs7Lv2OpJKTom8BfaJar\nVwBawWNjMEE7plJJSIuDAhYTnpbL5dwVywlTJmjHJNAlOizt2CRWcwO40VHTCqmJgpmW6K4bqBIj\nuiecexIDmo7RuiXj6rPkBpznZnZMuB5iCRMFDzEf49xoqGo1c49jPVrh+HcXamUtzxRpIkpy/FrH\niQGpxXtno8R4z2SpFVptDPCAbgeF6gAwsjpG+94vveFcHF/64+Qjdkzm/M5NhOoDN/0bsB2TGc6r\n62oPxNP3mIZzc48D1DRo+xfxSPXpEoeDgji2IUbsmMC6S+DeZbgfbn1KGHsc69EKcqy7C/WxoqeJ\nOMle3OC1biYGpL1Rwu/fYDI5bvZY1TGmmMCVEHF+ICZKoZBej6esEn+EtYZzPGE6szZnNjmwY+rW\nCiN3h4ozYDsmEM6bGgXNvsdqOK+eJbEwqGv/Ig7I3tSmg4Ic2jF1C+tyLTB4OfQ0gNqqEDIKnkoR\nPIZB0ITP0ZlRfVmXB52JAdkHfrE2Gpr8G5zEt/GgjUFZtzyXg0KqY4Bkqc5p6hh/OG+VCON2TK1a\nwdmtgDcv4lTYjqmv6VF0iZDOZoZzkAkAWsFvx6TtNZQOClLYMeG0hJl1WJlEePi92feoPfB43My8\nMg5i+ARdCPG17kwMiLQ5KJSnZwOaXkit4VwcYsKM3NWDHgcFc48DL1TqGH84784imH2YoLLSucdh\n9j2WD3yuLommpkc5lAjj7opg9VIrtNoxmfd0wHZMOC2x2mDpqEDc1STCw4BUE4jyp1ntw18/fHYu\nB08FInZQKE/3i4n/7T//n/6L9Idzf6p/OQcFT27TqWNOa8dU/YuXexxNNgbrtGPKrrO173HQXbHP\njsm8eTSQCOkYnJbo/sD3tyoMzA+SCuLoexRZTdPjhUGhcGl4cgNm/aM2sjq5f/fBjPeDqf5xMdFd\nLbFyOyZnwYTfxmBGBwVt5DHtmPwtl6aJgnNfPBUTwC+hPOK5O9S8TY+p7FiuZ1LqWqHBQSF+qVg0\nqG02m1U5HBwHCoULBIRzj7uiWfo04r80XgaIcwODdkwj1RKgS8JpYzBix9R3L8qA2cuQLu3PDcSz\n1mzH5G+5NE0UsqU9OqZvo8G0VQB9j02RvpqZOKYdk9NBQWp/iCWucLPZ3N3drc3k4AhQKFwUadvk\nou6KwL/BczoI557TBXo7LmfHhNMS1ZZLp43BbtiOqfVelBmzhHPRgzfuezyJHdNcYmKn3F3aY8c0\nmBvAfY9SS0s0bdmc0I7JXbJwvL2GYIUkJ7q742uvvfbpT3/6c5/73JHXDVAoXAiliJ4Wc1ccF+wT\ndFfE63rC+ex2TH4lEZd29j1mf7JF3RUxIJzvfO6KYPWpMFEwHRT8q7faMXk8kZwOCtrp6W9NHYOt\nmkfu+zDSzgDSEnH15fMK9d+H/+2Oa7A4TdOTJ0+ePn16EqEwTdPXvva1P/mTPznyuverUyicOyBs\nm/UKfbuVpomCR0lM0F2xOqHHEEnc+xdNdkxNS4/0PU4LuytiynCe6RjTXRE/7bjvsSxHMFcPw8LH\nt8eOye+JdDQ7JklSC9qeQrUsURx9jzK2Z69pBf/q3Sh5haNuNKRsNpvNZnN7eztN03JC4b333vvg\ngw+qauCFF1749Kc/TaFAmnHG41ncFc2eiHKYpxiiQyJEsCeSOe0SdkyxfmLwrgoLuSs6XzCD7op4\n9W5jx+rqT548iT9qWmEWT6QROyZnWkL0VH9pouDvexyk7LEUh3/DjKsfliycpqkhWCzHzNDsQuFT\nn/rUv/7rv4rIe++99+KLL+52u6AYvv3tb7/77rthu+GFF14Qkffff3/Gdf1QKKwOZ1Wt/4tCGeP9\n7ooP1eNWZsKpJJrwlx/6JYJz8llaLvsKDM2mR2l3Vwxz/tN73xORX/7cl8y/eHdiQGpawdko4el7\nTFdPVUIg0wrVesNuE4UOOyZ/WiK8g7pNFI6gFfz+DbOvHv5x5I2GSHaXphmFwg9/+MPPfvazYc6/\n/uu/fvnll2X/eD/96U+///77L7zwwgcffBC0wmuvvfbBBx9QKBCRxapqYxT3uytmuYFj2jFll+Gx\nY+qbf6Tl0ukE1VRg6Gx6lPYuiSARIr/8uS+V19MkJsx6Bdwo0WfsGCcsJUIkaAXgdlBqBWdtY6sd\nkzMtkemYbhOFpXcBTlKvUC29OnLMCp/JNzc3InJ3dxcOznIb6GmaXn/99T/6oz+KR1577bVvf/vb\nkjzkqBVEJCQbBhftvFQKhfWwaFVt9Rt/SutNH6TITIzbMYEr776DlGf+ETumWXSM1vSID3ocFDKJ\nkBLVgF9MpBdsVheK41aTrX2PQSI8ffpU0wrhLnz2fpzV9zhix+RMS2g6pq+2cYmuxVQimGWVMy6d\nZhHmmrOb7XabaoKnT5/e3NyEeoXBmb/+9a9/61vfCo/xU5/61Kc+9ak///M/F5FpmkJGIQwLWuFr\nX/vat7/9bQoFsmBVrdY2meHJDaR2TKVuGLFjwnjC+WAfx4gd0+CdMAftmLQ7SE3ThIVC+IdHTGhX\nXo3xMZzjXYby9GxMqmMyZRAEQVajkC7t2lBbzI4JpyVMHYNrGwXaMY3sK2UjW/0bZmFVKqFk3hqF\nmFT4vd/7vb/4i78IDzlWKsRhQSvI6Z4TCoW1sFBVrdY2mQ1L7ZjMu02KLibmNWzI8ITzvj4Os+9R\nHMmDvjthzrXRoHVSmNkCLCbEvQ1RrRsIwRu4K3rSEnijIWhrbWmclljajgmkJTw6ZsSOqfrAY12g\nuXq3f8MI2YfVSSoSTDo+n99++20Reemll9KDX//613/0ox+JyPvvvx/fib/7u78bkgpl6WJarzDD\nw2hkjX+JK2ShqlrtnTZZdkxVHeAUE+mtJud6dTWFc/+iU+GpsJwdU3pt2kYD7ns0yxvxp38sPpBa\nFKlqhTD4n977nqf+EYdz7K6I0xJgZnNA+rYqLzKIiUXtmHBaol4YNJMdU/rAq22GIFdkLg1O70BL\ndaxTK7Ty+7//+9/5znfK7rDXX3/9Rz/60fvvv68lFTJZME0ThcL1slxVLWCy7Ji0LonqgOxcmftD\npDWcd3ddVt0VnXZMO59LtDjcFT0/psfFHc5BbEi1QpQI6RE8uRnO/caOZTivTh4TCaZWKOfMBpyd\nHVM8y47o4K2gpxbA1sMSZRCzbJecC2lpQvZjllTIthum03kuUSicnuWqajFmvQJQElLTCv6uS/8V\nSm84x5+AfeUI2urZdwVTxzg3EaoJhrIcwQzncffBDPZhjChVC31aIW4BeOofQYmiVo5QHZAN86Ql\nxN30KMWewqB/g7Z0HICtmjU7pnubY/wutLYhFtplkEMFYGYOLiC18Pzzz3/44Yci8pnPfEZEwr5D\noFqpICIvvPDC1772tZBCqOYYjsbZP/sXwHJVtSZAK5iFjeKobZTFTBQ84bya13X2PfrtmKqB2dQx\nuFvBqSQ84bxMFWAJNVivUIbz7GKwu2JfOUI6wPljdpE4N9DhoCAO/4ZWOybnTsEfv/pKftek6nnR\n+dBRsnCELIIntXDW0WqapldfffXP/uzPSqGgJRVSXnvtNRGhMyO55zhbD3Gt8I9yoyG+as0uCYF6\nYiETBWc4T/O6s/c94g81U8eM9D16wrkW+Mtgnzko9DVBxKliONfqY7C7Iu579NQraCOraYkmq2aw\ntKfvsduOyRPOY21jRRZkWqE0R65N/lDcAJdu4vL2EQCxXDH8N1Yq/PCHP/yVX/mV9EkIQuGFF174\n0Y9+9NZbbwXnpVVBobA6jlajIIkakMSOqV5gpWcLPGJixsbFNEibWYfwj+P3PZobHFLYMc3e9zjo\noNBkqxDxRILJclcEnZDj5QhRSbRuNHgKG80H3mrHZIbzOHM+OB+R3HBRSTDk1Xa+pZ00bTRcAJ/5\nzGeCtWLoXPjHf/zHz372szGpUPZ3vP766yKSmi+tisv/g5GS8o3qyfuluwx+MfHyS/8uIm+9/fMj\nr7RJv4OUmS0Yb1Vo6nt0bnBIomNm73s0cwNOBwXQKFGSdR4OdknEecqR5kEzLdHtruhxUPBoBc/S\nnnAumUQAFQlAIhxOrmqOXq2gbTTMG3q2223Yqz3Opi0gJBJi06PsRcNut0vbH8LTEpREyCWc8JpN\nKBTIA/jd2yomgkSIvPX2z+PTB6+tyV1R5hYT8RPWv3T8vBi8g1Qazs3cgN9BIYqJJokQ8dgY4HBe\nzpmNaRITmZJY2kHBeJtAOyZn+aHa91geEYeSEIeYaNQKR9toyOrBj1AJDkj3EeLBWK4YkwqylxSv\nv/56uOPDmqFQIAcMKv3w0ZBJhJTB1IK59EhuQOt7nN2OqRo+++4gFU83cwPm/R3S8eK76YPozQ4B\nj42BVv9ozhz+MWLHNOKuWJ6ejTE1t7rR4AnnApMHu2Qk0A3Vg1Zawv+a3B3rNk7pdm3WbX4Ebm9v\nt9vt3d3dv/zLvwQP/mmaMqEQCxhLT4Wz4OOnvgCyLrrf21EihMyBxssv/fvsnx1RIsRkfhVcXViq\ngfRHoBU8YgK0SEiRWuiQCPHbfxXQGJmlFsoMxD+99z389wovGC1aP3nyBGw0pKmFMgMRkzR4l0G7\nMIGdFOEgyA1Udx8yJQGWFqi5H4p7qhsNWjHBdHAuXLt2uiQHgZiorp4caco7xj/icqExJA+iLDia\nPojL3d3dhZfTL/7iL4aH+frrr3/rW99Kh33mM58JR0Iu4exgRoFUaIrlZRbhrbd/vppUCBoiiIkZ\n66jndVcEemLQjsnzWJwjpdg40CoSmioMQE1DR14hdVAwuxU8XZdabmDcjmm5mz6Yz5tdgSj1r/KV\nrQf/LoNAMZGubkmW7HpOHlBC5D7OZaSOuiIyTVPc9ciSCmWO4bygUCAVnO95sNGQaYUoEdIjWheZ\nZ+k4cq5wrg2uTt5txzQLoO9RK0fwzDZLY2T8MYv9mlbw3OwxXuRydkyD7oq4ZGFUK5jGhSCcV8+L\ngz1iwlGycOSNBo0gEWSZu++WC+12u7DpEJRBEA2y/2OFSgXZlyvKCiTUCBQK5ABnF5NZiyB7rVBK\nhHRA2S/g7JIwGyObwrm4axvH7ZgGMSsQtXIEz2xmo4RHK2ixv1QSnouMEmFpOybsoDBy04dOrZCE\n836pMekph+qA2upNOubkrY8hYC+kFUI/RTl/LI/IkgoisvJ85jAgAAAfpElEQVS+RycUCuQeZxdT\nFs5x3aI4xETcxWztkgBawR/OtcGaHVOUUCN2TIMArdAkEcSRlpCW2kZPXwy2Y6qOjz8ubcfUYdW8\nrzDYidX32BbsiwAP0m9GOAc7CMpa+cjq5Hqq4+RaIbxslmiSDJWSWVdFKg6CTHny5MnNzc1//Md/\nZMZK5wuFAjGIb3sQzsHBQTGBUwvl7sNydkxVwTRixzTeXaLdxskTzv0OCh7nJaeJQuu+0vHtmJqt\nmtPHYvU9urSCFtFLQyRpzA2gtVvExL6MciWxIysUkCWFghzKAhF59uxZrGGUfWohdkCce2lChEKB\n2HiaHqMsKEeOiAlP2ePUZcc0S27A499Qfr5/9zvfFJEvv/qNwXffpN/GyWyGHLdjMvseywd+YXZM\nUj6KQiscjHT0C5g1AWgYyA04dhBcxY9hoLU1eUzCt/yoDJZujyznD9IhJqvSZEPp1nymrOIvTc6C\naZoGcwNNYqLJoMls2TJtDPpyA9rp2VTx8z1IhMiXX/1GeT1+PLmBNJw77ZicN30I/zBNFJxKQgox\ncXw7phjOzXKEikRIHsZB/2Q2Ur8f0sNZYBNBBnID5kG8fxEGWluTJyEWEgaWSycEqkkFbd0f/vCH\n6/dTMlnFn5mcC1WtkDY9mhUJ4hATrRJBHC2Xnu1znBuoVueJyHe/880vv/oNsx4ikwgpHakFZ25A\nfPd9wGkJXLIwmOo30xLAXXFGO6ZqOAeWi9M0qVphegiq2tJyuJ2nJCHS8fJwZDA3MCXzlCOrBQrh\n2GL7aHORGSosutCRPZ1ODoUCacC0TPA4KHjKG7svBodze5NYzw1kWmGapjT2a1oh5jCy8Rn/z3/9\nQ3Bh5UWWwbtalhiXNu82OViycJzbOM24y/Aw2J0bUGwMDmoUHo4AJSEPYsJIXSyZG/A6KPj6HleS\nWjgasc3hxNdxLK7rr0u6ydomQd+jx0EBiAl/pZuz5TIL56B63MwNgE0EKbRCtUChepa0PGpQWJBt\nImRLAzsmj4kCWFqS1ILfE8lpiKSdXv62w44pdCuYEd2yMdhJKRpE7INwj0PwNoT05gY8aYn94LLK\nBD8V1xNNjpa9WAlX9KclfWhtk84dBL8d06BESGfD4VyrLowDcNcl3kTADgpO4aIxYqJQWiZIUaCg\npSVMJSFKaiFi2hiMuys22THlmYBaWeL9MMeXaUNqxByDVHQD8G9w2hj0uytWTywkAiEUCqQToBWc\nDgpN5Qit/g1mTYCmJKSmFVKJ8OVXv6FN7kkPBK3QXcY4YqIAuiTiDEA6YJfoMHO3J9K4u2JnOUIg\nC/bFj6NaQVtXcq1QL1AoT42x35MbANsQRcrBs9FArg2+IEg/1SoBaQnn3VkE0DoR90dwOAdKQhKt\nUM0ilDO07iAMmiiMeCJ56hU0MVF1ic6OdHsijbsrGo0S+DkHmwjSqxWStIRZDIHuMV1GdFNJiE9M\nJCmHpo2GVrbbbegI2Gw2i7YkkCWgUCBDpE0H0uWgYH7++lsuq+VmZTg3EwNxpFiZiabcgNNEwfMZ\nPen3exQrtWDuX5QTZgPMtMSIjcGguyLowxz63i/tWgH/WFkazA03EaQmHaRBTCwdBUL74s3NjYiE\n2zEw7pwXFApklFnsmPCXYLzHYZ4ew7mZGJDD3ABuVeiTCOB0vx3TtC+b77uNE7BViA8cb3CYaYkO\nG4MZuyTkUCJEKq0K4u5W8Ak4nJawlgZTi0hvOUI64HDC42w0ZFaJizonkiWgUCBDxK+/43ZMZlwc\nsWMyEwOiBH4Q4/0Rva9gAlyP30HBsw1h2jimE85ox7RE32N0ULBbFcTRmFAecbkrwryFvXR2ijwc\nNPse3WmJTKEuFwiC60A6/7X1Fl4AFAqkkzJDPm7H5NmGaJII6bkjuYF4upmEaF1aHJsgTe6K0liy\nAK5cYP2jaaswb27AmZbwJ3hEHBFdD+coVzFsx9TR95hsr6gTaxsNR65eZEbh7KBQIM2A6JJphQ47\npuV8XYCNgd+Oqc9gUdMKUWF4dMyIJ9K8tZNg6Xk9kZrqHxskghnOpaFkwd7OKI/47JicGw2KZCnG\nr6AyIGQX1nAlpAkKBdKAs85u0I5puddkaWPQ5GrgyQ04tUK1QAHrGI+NASgs6GvFFLcd07yeSDH2\nT4d3mq5O3iwRkqPGwWEx0WXHFP5t7yNY2yv7eZJ8wwk/8Debzd3dXXaPZnIWfPzUF0BORmvDkvNT\nBtcrxD0ITUm8/NK/d3+cmTomfnev7iB89zvfxEuDr/5xQm2GbOlyTDl5qmO++51vYhuDWOhQ5Z/e\n+17Ts5ptc2haIXVQALPhpsfSBSGN/UArxDnxQ0N1A5kaKEfulM7GtO8RoJ2uLp0OMFoVZIfekvH1\nkA6oHjwCIZFwc3NzVfdHuCSYUbhSoroXkSVk/rgdU6tlobS0XI7coqmaG3BmJsxoDXRMmLmsLszS\nEoNJBWAUbdoxDborgk2ETCtUCxSqZ8WlbRMF0XMDfX2PbhMFEcvnIP1R8iPaA1/PZ/t0eLtFcnZQ\nKFwp0zRFdb/QzdBObsfkKU4s8dgq43DunESbWZszzgx8mqWmFcz2h3TpQTumDnfFdDD+k5l2TJlm\nyi9GD+f3X7X7NhGkJh2yYbYdk7qyp+8xPvBUIqzEYDF+vGTHmVo4I1bxSiJHpmxYSnXDLKRp8IXs\nmDz+DU2Ni4O+CBkd93GQ+VyiQd+jtvRcN61uclcUd7dCTA+AMUEraCNzrVBr2zG6FbqVhNTERP6j\nNrOIOPoew6iisPHkn/DBbSk7yATDeXH6lxE5OWEbYq5XQpb2XMKOKdtowHsc/lR/q0TweDsO2jEB\nW0mPHVNfFsHTcontmJx3g+yQCONdEjGca/mYztpGU0mIJSZKreDse4QdDWvbhiDnCIXCVRMkgojM\nlU7QvsHMZcekbTQMukTPGM7FHdFjOHe6RHfYMTkTD1XMlkup2THFP5MZzsuLr16hLGDHZPfBarWN\n1YMzd0mEf0h5nZW+xySRYFZ08qOedMNXD7nPDS5dkDxuxyS6mDDTEurXR4eJgqduoMmOyZ+WMAd4\n7JhMP6Xqr2IjQ9+NrUGM9+uY2e2Y4imefMwCfY++tIT1J3soWWgxWKRWIN3wpUNE9h/KS3uljdsx\ndaclMq1QLVDQqseDkjDDuQxYNeOWS88eR99NH6SmMPxWzZ7J03DeoWPmsmPyFEbMa8c0kpaYZV+J\nkLmgj8I1ElIIx/+giRYL1aBeSoFsmKYVUiWhLR3tGbT0QGmBkHZGgGidhnPw2MP81RnKyTPZoa0e\nhr3zzjuf//znwdLYRCEmJKrtDMBtKQ72dPNXZwgHzdM1O6YoAjRAT2Y4OAF3RZgwSOoGlDcRLmxM\ntUht6WphQTzI3AA5MnzBXSOh6yGmEBZqj9Qw6xWAHZNYeQizttHTqqCJCRzOywHlSNzOYLZcZmu9\n88478cfPf/7zOJyPlCyUlglSJCG0fIxp7Cgi/8ev/V/4wkCPpSYF/P4NIg25gUrwNvoeHeWN+m9T\nHePfZRjh9vaWd2EgGRQKV0rWs3Sce7SkX4kG7Zj8LtFS1Daaafy+6sJ0AJjQFBPa0tnqqUoIZFrB\nGc5l/3cJmQmsFUB7ZDZ5JjtMW8nqoulsoOuylA7i8G94kAiecI7bCoy+x9qJLXZM2l9tiU/v2Tul\nyQVAoXDVhFbmoyUSyioB7c7R410SZsslKDnEdkzlueUY7K6IuyRMHVNKhEiI9E3hPJvQzExoS8fJ\ntb2GalrC+fljmiiAcoR0QPqjP5y7ujOMvkelHAGsftxahGDoPm8PFLkYKBTIMah++4lawTRLiDi9\nHWOxAhYTaatCdXXNjgkoCXGnJbSRZuvENE2aVgjFCv5wXp2nmplISxnAzHh10CWh4TRRMAsb5VBM\nOMN500XOZsdUOJEc4VN6u92Grw1HaIAiZweFAjklphtMTC04vR1b0xLm6podk2mrYLorlqdnY0wd\nk2UC4hGwgyD7cA7SEnKYmQBWDdmcppgYzyKYfY8zlCN0S4RxO6Y9JzRY5NYDKaFQIKfH7P8Wx0ZD\nd1pixI4pphacuYHWXYbyEWWrx5aHLPBXUwtZOMdbGObSwCXabLm0LY/G+h6xHRMO50NZhG47ptU0\nPVIokBIKBbIKzHg8crdJT1qi+43gtGPq22XwrG6WLAgM51WFETMTpq2kKSaAsePSDgrayEpqod3G\nYDY7puSZXEPfI4UCKTn965IQadmDyHCGeU9aovW90GTH1L3L4LmMwdxAlAVS0w2LukSbIbka41/+\n8Usi8tYn357rvg/ALCG7njntmGpiy3wjLA2FAimhUCAnJv1kbMortGYCBtMS5WxzOShIr0RIL6ba\nLSn7pkdPqh8USGKXaNByaZooYAcFSVILqUSIVLWC00FB8DaE5DYGdm4gO9jb93haKBRICYUCORma\n/ZypFbo3C3BawrzpQ3rZc9kxDUqE9Kqy2sbsR1MrmGkJUKJYfic2TRT8hY332aBDiRDJtEKTHdN9\n98d4bkAbqfc9rmSjoYRCgZSs8ZVKrhzwATqemK2mJcybPqSnz2vH1FQ65wn21R0E6dUKcefCs4kA\n3JpLE4WmP+I0TZpQkL1WaLJjypWNZmOw89oxoSpFpe9RjtvO4IRCgZSs7mVKiNQ+QNMig456gmyq\nmJbQwnmfVmi1Y/J/nzarC+NgnBhoOj0bb/Y9mv4NrSYK2bWBpIK4uyTUjZKswBD8mB4Xr7tioOom\nwg9hsnJ4UyiyRtLEbFmHGO7wJAOpBWzkXN4gqrw26e17FHcioSwdeOedd/ADD9eGPRIA8fRqCWQ1\nr5C6Mvzy55B/g7QnEsrVyxqFeJ1mowRY+mEPopoeKKVANgxoBagSZMV7EIREKBTISonxuFpVEA62\nfsKmE375VaQV8OSpP2OZRcDNkM4aCFGqC8NBj44p7Zg8zxWOW6lWKI2b3nnnnV/+nOrfYN4r0hR/\nqVZIJUL6W9AlgV8tRr1ClALaRoMmJhwei9QKZOXw1UlWzVxtjcB5Cbsrmv17M/Y9ZrsMI02Pcliy\nMOPbHOuYeHmpmDBNFNJyxbc++bbHxgCoDdwlAc4VrV5h/7uwsHbu/Ril71Echbr8NCbrhC9NslKy\ntsnqXSJlf6sIM2T6b/rgSQOkEmGWvkfT+6h60FNzMLtE8OsYse778JA0coRzp5iISgIXP9paRNto\n8JUj8HOVXBIUCmR1aG2TWbdC9iPo5g9iwrzvA473acFEKSZG+h7NCsQmQ6Q4pzhMFJxKolvHmCUL\nnnDuFBOZkugWCg+zaRsNlkfTQrmBcINHEbm9vWVXAjkmFArkbEi7FaoJhkrY0LVFPCh7JWFWIIou\nJrr7HuPVml/QTUOkOJUUBQrVYntx2zGN6Bhw0+qXf/ySGdHFkRsASqJ16yHFtliIOQYl/zHjp+t2\nu338+PHNzY2I3N3dPXnyJIgGQo4AhQI5J8xNBFwCWeYhsh+Bu2KY2Wmi0Nf+h90Vzb5HYLBYigmP\nHVMqJkZ0THl3qDR+V7VCDPNYSWinl79tkggRjx3TcT5Cp2m6ubkJd4K+vb3FTRyEzAuFAjkzgFbQ\nkg3ZGJyWKJsesxCLTRQGDR6wuyIIxuK4czQQE560BLZjMmsbS4kQyYI9/jE9Hlc30xLdfxdgx7To\nRkN5GakPUqjZZFKBHAcKBXJ+VLcVJLFjMssRzNpGEPVLrdBUsehvVagGZlwT0P29X3xpCWzHhHWM\nJzGg7TWU0kEKAYfTEoMCDrgrygIbDRlh3yF7vBQK5GhQKJCzJK1SFKXvERwEYsJzI4moFZokgqfl\nUoZrAvDB8bQEaLksrzxberAcISoJTcDhtES3f0McEzj+Z2ZVKMSdCEKWhkKBnBnxY91jsRBlQTmy\nPNh0rynnl0itSwKfa347Fys30CQmWtMSnisHS4OqQ1Bw4NlBCFoBpCUGWy5P5YxEoUBOC50ZydmQ\nxWZgsRABGw1paqHjdpTOrsJqQUM4OOKumH7BrQLSA6kgKIdpGx/pwdYrz5f+5OfTAsM0olfDfHrE\njNNAarz845emzBf8cOlS/JVH1uOiyA5JcjQoFMgZAL6+V1ML2S4DmNnj19QNbpQw7yghuruiFo9b\nxYSpJKSmOfAtJ+K1AcBGQ5paqCoJzWZDHGkJqUmEbOnwR8EPED+62QmaYLvdpuKAQoEcjVVIY0I0\nmjL8oGRB65aUhbWCKHmF6MrgKVnwp/qr4XPp+kftwrSlPeHcHOCxY9LSEh7zhvVkDgJRKwjbI8nR\nWdE7gVw80Vpus9l4CrabPqlxyUJpmSBFEuJoWsHpzuR3V5SxuoFu/waxbKSzpctwDvoeZ7Fj0koW\nnHZMq9IKaZ4mbZUkZGlW9DYgl034GhSt5WSBFK5HK4CqhSNoBc3eMdMKS7grip4bwIYNmlZIxYSp\nY0CJova9vzogG3YEO6ZVaYWQUaBEIEdmRe8Bctlknd9LNII/3BAIbjSAGZbTCtX2h5SgFU7oruj3\nb6j2TWAd0/e9Px3g/DE9Lm47JjFvKcn7PJErhkKBHIOyvyvdcx0n/SgftGNaKM9h3moyJBuO764Y\nf4VXx0qiunq2NO57PJUdUxxmZg5WlVog5JjwpU9Ow1wZBVQD327HNO/bAZgoaEbRptuSzOquKO7a\nxllcovscFMQhJprsmLRhnmoPfmCSK4Sve3JsQnZBFs7lttoxHUEipIA7SHn6HsfdFbXB43ZMeIDm\noOA5OG7HFP/Nzz1C/FAokKOy2Wzu7u6OZipn3kFKFosZwEFBWswZnX2PreFcejcawEGnjhFHbkDc\nYiJTEh53RWFugJAW+G4hRyIkEm5ubm5vb49Ztq2VLBzhld/n35ydboqJ2cO5dno2oEPH+HcZZKzv\nEbsrxuPH/PS7vb3lPZzImUKhQI7ECd3ps5KFI7/mW+8g1Tr5jOFc2sWEU8e02hh4WhVMO6ZVdStM\nh/eJJuSMoFAgxyCkE549e5YdX/pzM7uD1Kle7Uus7u97xEubexzYjsksC+i2MTAtF/3WjSf8lAsm\nY8E4hEKBnCm81wM5BiGREGoYI0snGNII4bkBwXIsJBHCj5pWSMM5vjbzxg1aj2W46YPh3/BJtTES\n9ylkqQXtDlJm3+PJZeJms9lsNk+fPj3J6oSMw4wCuUyqAezcS9j8fY/ljyMPHBRDlJNXCxSwu6Kn\nW8Hj7ejpfTgV3Hog58t5f24S0sqIVjj5/oWn79EZzluXrs6ZTY51DHZXBGkJ500fIuv8TKNQIOcL\ntx7IddF3V8AQsX7wxp/K6dIS5S5ACahXANsEgDT242IIfG1xD6JahxhuHl2mJeIwoBVKD+bj/4G2\n2211E41tDuQyoFAgV0fTpnUqEQI/eONPR1ILI+dq39r9Jgqt15lOVU1XNC0NEgNRK2gbDaVW0G7T\n0KcFCSEafDuR6wWEk1j5mEqEjC+88tWmt0+qOVrPLafqcFf0C6Nyzmz+Dv+GeMFz3fRB1rrLUIVb\nD+R8YUaBXCNp22S1Ky+GczAJ0BDV5WZMS+CCgPT7fatESPUHQEseZKmFUjeEPQiB7oqmVpCkJ+KM\ntAIhZwrfZuS6KMNzDDZlOBeRL7zy1aogCBrCfPtU58zmGakx9NylyTOPdOUGxu/70HrTh/PtZGFG\ngZwv5/EeI2QWcNO/Fs4zrdAqETSpEfkvX/l968LRKlqXgf8i/d6OpW7A/g2mjum2YyofBT/KCFkI\nCgVCRHxawSkRwmyZthhJS5hr9ZUjVE8vwbkB078B6JggJsy+R+cuw7mkFgg5O/jWIuQeoBWaJEJ1\nkr60hJPUEKljTpyW8LhEY/8GrUsiHjHdFeN1lgezAfxAI2R2+L4i5IEyzIeI7mxVyBIJGa1piSYG\nY2Q1LeHJDVQHlwOAkpBDrbBmd0VCrhMKBUIOSGsLpFbbKNaXWrB/gc89LdinWWpawXmrSVNJxDEB\nGiEQsir4biTkgfTGUd2tClpaYv3vNb/zo/O+1U1KQopulHKJJQg3eBSRzWZDL0VCSigUCBFR2iZB\nbsDTULDcRsNyYK1QjejVGZqUBJ5wUW5vb58+fXpzcyMi4WbQZ/THIuQ4UCgQgtomW7cePNOuHI9L\n9EizYqok5NQbDdM0PXnyJCYSsh8JIUJnRkJEj2ohhqUujX6JIHv7xdVW6WurZzsFVXfFqhpIFQZ+\naFm+4VRaIdzJKZUFNzc31ds7EXLNUCgQgmi6g5QUBQqaVnAqCVlGTMSIDib3uERnppammFibM9Jm\ns8ku5u7u7smTJ6e6HkLWCYUCIQYxZGI0B4VSDaQ7GkAr+MWEnyyia1ohDsM3fYhqwBQT8QkEuqTp\ngczOdrt9/Pix8N7QhBRQKBBig3Pj5g0dYrCvjqwqCXGkJZrQInqpFdJ6RuC2lNoxmUvH+U+1xVDd\nUIiaYLPZ3N3dcd+BkCoUCoQY+G81ie/pANoosJKQYa2AexmiVqiKiVIrZFULmpiodjS0buUsTUgk\n3Nzc8I5NhGicZVU2Iceh41aT+CAQE9VzyzELaYVZ3BXj8ZM3PfqZpomJBEIwFAqEVADfej13jk4b\nJUwxkfZTmCbQ2lV58PQ9zuiuKKdu6zAJ6YRnz55lx5laICRl1W9jQo6PJzFu3moy/KPvvg9mrmIk\nryCw77F6cMRdUdatFYLbUnaQCQZCMtb7HiZkXm5vb82Cdn9U84TzkdxA5t8w4vukXbxzo2HcXXFV\nRQmEkFYoFMi1ME3T7AVrnnDeHeM9exwLlSyIw34xc1cUS2atObVACAGw64FcOOGWP8HGfwlAEWKW\nWuiQCHEXowqufMSE7gPQ94ivs9q8gJtIeU9IQs4Uvm/JhRN76J8+fTpjRiENk7gCUdrtn7NSx+rk\ns9xuKtMKHokg3Eog5MqgUCDXwlxbD9ptDkA4D7kBZwAu58m0wrx3pAQ3e9SukLkBQq4Kbj0Q0oZ2\nI6UyEyBud0VcjhB0xkI3rXZG/UweHVkrhP0jEdlsNrRYJuTIUCgQMg9RK2hVC1Wt4Pd2HOyKBHh6\nQcsxR9MK0V9ZRJ4+far5MRNCFoJCgZDZAB7MgVQrlCOrIiM9YrYVyKz5hpXcyenu7i7uGQWLpO12\nS08kQo4GhQIhcxK0As4NADGRphZK0RB0hkCPplm+5WernLAiISQPoiygPiDk+FAoEDIzmlbI7Jjw\nJJrUyNRANS0xklqonnvC6sXNZpOuG4QC5QIhx4RCgZA5qbZNVnMDph2TebfJ8vT4W2lPA2B5cfIe\nh1CpICLlrRkIIYvCHidC5sHcFKiC7ZiqWiFtuZzlbpOZuDnJZ4JWopj1OIS7M/CW0IQcE2YUCJkH\n0DYJzsJNj9Wuy7TlEtsx4Z5MqYmbquXieghCgcWMhBwTCgVClqUM9pE0N2CeXt2qKLWC078BqIFT\n1SJUY39QButULYRcCdx6IOQYYDsmgdsEzv0LcNeJdPI1bDT4Cf2QT548CXsQ4UduPRByTNb+MUHI\nxYDtmKSmFZx2TGDObIzUyiBW/iEQkgrxxygaCCHHYe2fEYRcDLFYwWO/WM0iYDsmU0mAjMX6Pwcy\nQwVCyNE4gw8IQs6dLNsvulaI3/v7Nho0JSGWF8JZaAVCyEngpwMhC+LvmXTmBsSx0ZAqiWx1agVC\nSCv8aCDkNAA7JnzQtGPSyhHEcbcIfiAQQjL4uUDIyZjdjglIBEII6YNCgZBTAhIDAbMOMbt5xGn7\nHjebDe8BTciFQaFAyInBdkzmOzR1fszKETynz0i4HQNNDgi5MOjMSMgpCeG86q7YLRG0I4uy3W7D\nTZsIIRfGx059AYRcKdM0hQ2CQBAHX3jlq8FKAYf57NxjXTIi+Cee+ioIIfPDrQdCTkDr/RfMMae9\nk1O4U8Pt7e00Tdx6IOTC4NYDIScA3I1Jq0OMGw3ruZNT4Pb29u7ujjWMhFwqFAqErIvSQjHLFpyk\no+G999774IMPyuOf/exnnz59+uzZsyNfDyHkaLBGgZDVEbSCHNYilL9dA+H+TLe3t5vNJlzV48eP\nufVAyCXBGgVCVspZ+CRut9t00+Hp06c3NzehXuFk10QImZUz+CQi5NqIew1noRVSWMxIyOXBGgVC\nVkRWjnCOWoEQcmHwM4iQVQD6G6kVCCEnhB9AhJyYEfsEQghZGnY9EDIn2+02uA85fQXKpoYqYcx6\nmh0IIdcDhQIhs7Hdbh8/fhz+/fjxY0/lf1OSgFqBEHJ8uPVAyGxM03RzcxNyCbe3t0+fPuX7ixBy\n7lAoEDIbWXPgNE1Pnjw5uaPAdrvNroF2y4QQP9x6IGQeQvRdoYUAbwBNCBmBPgqELMgavrtvt9u4\nIUIIIa0wo0DIhXN3d7fCPAch5FygUCBkQVYSobfb7bSHqQVCSBMUCoTMQ9AEWRheiVAQkWfPnj17\n9uzm5iY2cBJCiAd2PRAyG6lWOHJ7ZHYXx0jZc5H2cBJCiAmLGQmZjZjhDz8+e/bstNdT5ebm5tSX\nQAg5J5hRIGRmVtUnGUwU0vwBMwqEkCYoFAi5cFLfp2AynbpCEUIIhkKBkAsnVEvEH9dgFkkIOSMo\nFAi5Cla1IUIIOSMoFAghhBCiQh8FQgghhKhQKBBCCCFEhUKBEEIIISoUCoQQQghRoVAghBBCiAqF\nAiGEEEJUKBQIIYQQokKhQAghhBAVCgVCCCGEqFAoEEIIIUSFQoEQQgghKhQKhBBCCFGhUCCEEEKI\nCoUCIYQQQlQoFAghhBCiQqFACCGEEBUKBUIIIYSoUCgQQgghRIVCgRBCCCEqFAqEEEIIUaFQIIQQ\nQogKhQIhhBBCVCgUCCGEEKJCoUAIIYQQFQoFQgghhKhQKBBCCCFEhUKBEEIIISr/C7sfambVs/Sl\nAAAAAElFTkSuQmCC\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"jetComponents.reset(getJetsComponents(input_particles, fastjet::antikt_algorithm, .8, 14.));\n",
"jetComponents->Draw(\"LEGO1 0\");\n",
"c.Draw();"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "ROOT C++",
"language": "c++",
"name": "root"
},
"language_info": {
"codemirror_mode": "text/x-c++src",
"file_extension": ".C",
"mimetype": " text/x-c++src",
"name": "c++"
}
},
"nbformat": 4,
"nbformat_minor": 0
}