Notes on BaseQA with HFT hits
Updated on Tue, 2015-09-15 11:24. Originally created by jwebb on 2015-09-15 11:18.
Goal: compare TPC only tracks which would have been reconstructed with HFT hits to corresponding HFT tracks.
Method: Reconstruct events with HFT in tracking. Create a TTree associating ID truth with number of reco hits in each layer. Process TPC only tracking sample and apply cuts based on ID truth / TTree lookup.
Figure 1 -- Shows corrspondence between the number of hits in the TPC (0-49), SST( 50-54), IST (55-59) and PXL (60-65).
z-axis is cut... < 100 entries are eliminated. Note that the z-axis is scaled by 1000. Looks like a good (enough) correlation to proceed.
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAhgAAAIoCAYAAAAr/pK9AAAABHNCSVQICAgIfAhkiAAAIABJREFUeF7tvVusJdtVphl5fDCIi31cth9ASD5uEBiXhI/KhnZLSJlZolVgU/hW3U8NHGxskGxKllKi+y0zH3wpiRRG1UiAO+1jqlotUSowIFRVaqTc+UZLtjp5QCAEpeMSEgifYxubO4Jd8cfe/86xxpozbisu8/KHcuVeMWPOMcf4xlwxxxoxI9a1Jrydh4tVKgIiIAIiIAIiIAJHBK75El/QBRY3b968qnft2mGVU/ch+FQZvv0aMtXH46HiWQzt5+qPXPX2/tjKjvPzw+8hS+/DjqVlenm59pGr3qXwL8WOJcfRgwcPIA7bVdBgo4fzZ555hhUuarng4uAgJQ3UCZ38rJyh45226kP+cIPv1DGxxLhaQsaQHSn3ETrJWjcNHUfdoTqnHl+ijyVkDNmhPh6PHLEazyJFVo8ePYIBXWzBAOMouHhsot6JgAiIgAiIgAiIwDgCDDKeGFddtURABERABERABERgPAFkMJS9GM9LNUVABERABERABAYIIIvxZF+dj33sY32HB4998IMfHKyjCiIgAiIgAiIgAuUR6A0wvLlf+tKXmrOz55tHj750cOiZZ55qbtx4unnqqad8E+2LgAhsTODv//7vm7/5m79p/uEf/uGo52/4hm9oXvKSlxyVq0AEREAEliYwOsD49KcfNZ/+9Jea559vumACL2537py1xx41b3/7U+3r8E6UpRWWPBEQgX4CCCy+53u+p/mxH/ux5ju/8zubl73sZV2D7/3e723+7M/+rL+xjoqACIjAQgRGBRgIHj72sYusxdnZs22m4ksNshncnn/+2ebZZ8/aOs+3x54/CD4W0lNiREAEDAFkKP72b//26DbPJ598svm6r/u65s1vfnP3+uM//uOrAEMARUAEyiDwla98ZZQhyFhO2XBO+bu/+7uuiW/LPl/60pc2X/3VXz1K7KgAoy+4QC/Pt2mNj33smTaLwczG86M6VyUREIHpBP7yL/+yCyze//73N+985zubb/7mb27++q//ur10+aj50R/90e65Mb/8y7/cfPzjH29e85rXNL/yK78yvRO1EAERSJrA/fv3e/V7z3ve03s8dBDBBeWiPYMMBBe2fGyAMXibKtZccHv66eYgc4FjuDyCtRcX6zPOupOcNhEQgXUI4PIH1lj8zM/8THPr1q3mta99bfOnf/qnzec///nussjLX/7yruMvfvGLV99E1tFEUkVABEomgIACgYUNLqbaO5jBwJoLbHfu3OgyFdzwHos9z85ut+W/dhlg3GiebqOQp5466/5qEwERWJYAAoxv//Zvb37gB36gE4w7vZ577rkua6HP3LKsJU0EUicQy1IMZTdidiFjAZls7+XYrEZMhi0fDDBiQp59FkcuFnRicWfT4PV0e7Jr2tfTbeARa6lyERCBUwi84Q1vuGqOEwDvCrFfAE6Rr7YiIAL5EIitlZhrgQ8yKGdqcIF2XYARup3tsXJnTXP2XJvBaEuebl9n15pn3v6GNoDANycEE8huXO+q37nzuXaB52u6V9M8HJDbNdEmAiIwgUDot0uw4BOLO7k98cTFlc9//Md/PJIcKjuqpAIREIGkCXDOxl8/f/cdG2uUl4l2ob6G5HVnpdhJ59FTj5rnPvlc03yqrYQfSrvbvu6cN4+ee9R88JkPtvmLD7YF19s7SD53dfnkIsB42PUbk9sd1CYCIjCLwO/+7u9etcO3ip/92Z/t7ijBJZLPfe5zDVZ5Y7FW6IeQsEAUQQoWael5GLPwq5EIjCbwF3/xF6Pqfv3Xf/2oeqyEzzZffp7tOzamE+iMy67+/PHJT36yneufbabo2pvBQIDRBRXPX6r1bPv3ufbV/n101gYfH3ymy1ogNXuxyBOXSS6CC7QIRUEo1yYCIjCfwB/90R81v/mbv9m89a1vbfC0XNxJgizGt37rtzave93rus/du971ruYDH/jAwe1kv/RLv9QtEP3Qhz7UZiDP5iugliIgAqMIYPL/1KfwDT2+/ciP/MjkuRJyEQDg75e//OUj4Tw2dQ7GFxDoy+ACumGjDQgyUIZb4cdsvRmM5s6FiGfbP59s+/nRmxexRndjSfs6e/tZg594R3CBZ2VcXC650V4iuWjnI6uLUv0vAiJwKoGf+qmfav7wD/+wedvb3tbdpooA47d/+7e7rAQ+d694xSu6W1Tt9k3f9E3dLu40QRZDn89TvaD2ItBPgIFAXy3UmfpZpFxkGkIbA4y5ciHzh3/4h68CiR/6oR9q8AUF2xR9r7X1z7/lW74lpGO3khS3nT7V3Xr6XPvv6TZ58XQXVOCFQOIiuMDdJBdP+Pz0p99+datqbIVrsDMVioAIiIAIiEBBBP7qr/5qlDVf+7VfO6oeK60lFxlOvLB5ndjnV33VVzV4DW3ItPZnMFoJCCTa/5rrZ88172gfA37jS890wcSjR893GQtueHS4DS5QPjV6eixN70RABERABEQgbwJf8zVfM8qAqXPlWnKRAeXaLK+T7dMfixnZuwaDjZCGef7p9kbU9oFab3/70+1isufb140rmfixM7z8Q7amXv+JKalyERABERABERCBvAgMZjAuzPlUc+P59qE+7aUSrCJ99tmnnZX4EbQvHVk+Nso5aqgCERABERABERCBrAl0azC+8Ru/MWsjpLwIiIAIiIAIiEA6BP7kT/5keA0G1MU99t3WrjxHRNK+efy+LcN+98f8vap7caBbtX7Rrqvl2mP38njXja/7uE3Xy4zjF7pRe/RPfS/KLvS15eOOWztJ59D2x2w6s49sDx+3dpJt2PZLm6Kc43Zc6HJ5/EI50L3k24G+xIS/ZIP62kRABESgJXCOf/iv+7/7e/EHf7uDF8fte9MmWPdC0NWtklfPY+j6uJDZSb3s6+J9J+lx2VHdi+NWz0NZx8cvJTrb0AfrXtrmbB46/pgJRDtOzvYLfQ3HqcePmABThGNXfnH8kGmIzeO6fT98NvISSdtBN9lcTj7XnugmIb66Y6bM7l+8f3y8d5+T76W8x3Uv2vft9+nAgIQ69uqAifrKtvaJiMZOtD/cR91LFgE2gzoZbr06HXAxfV7qdqVDYH+UDs2FDb06XHLBUNMmAiIgAh0BTErd6+K5DHb/8XtMRu1TZa/qXtTv228u5V3IxgR7IT+0zwkzpEM3gRsd+/pEe3v8qk9O0kc6PLb7cT/HHLwOfv9iMn/MsU/Hi2OPWfTtT+FypFOvry780clvp6O+pRD4+tsFLdpEQAREQAREQAREYCkCXQYDm57stxRSyREBERABERCBegncuHzaZverSBeXEOqFIctFQAREQAREQASWIcCYQgHGMjwlRQREQAREQAREoCWgAEPDQAREQAREQAREYHECCjAWRyqBIiACIiACIiACCjA0BkRABERABERABBYnoABjcaQSKAIiIAIiIAIiwACju001dBcJbzOxqIZuZWWboXpz8I/VB/XW6N/r7Pvx+77+FvuW0RYMtraJ/c21bc3xuQUL9SECIiACcwj4+cnvh2Secr4clcHAidy+rBJ2MmP53BN/yLhQWZ8+ofooC+kZKou1n1u+RR9WNw4YMhrTf1+dvmNzmcxpR3vQ1r6fK2tOO7URAREQgdoInDKfjwow+oCe0nmf3KWP7aXnXv0uxS93/ZfiIDkiIAIiIALTCBxcIpnW9HFWYOwkZL8Nj20zRafYt22W2z5ZFjo2pc++ujHZIQ4o8/otycizoew+DnP0D7UJ2dvHbeqxmPxYeUz+1PoxOSoXAREQgdwITDn/2bqwc2iuiq7BQONYxxDqO4pBRb2lJtCQPiH51CWkJ8uGwMTs8bJD8kL9hvQ8VYchHa2ufE89Qnp726z8Pv3tMb7vqz9W7756Mfmx8pisqfVjclQuAiIgAikTwLnOb3POf2Pm88FFnlBkqQkwZJg3dMz+UvqM6WuNOn0cvKPX6H8tmSG/2CAL70N11tJnitxc9Jxik+qKgAiIgCdgz8F2Luqbl7yMsfujAoyxwobqpTq5DOm99PHaONDe1IOnXPRcejxKngiIgAisMS+dvMhzrlvWiJbm6oJ2p+jDb79zHHRKv1PtRV98xXRdWh8rj5ym6r1F/Vz03IKF+hABEaibwJh5YMx8wgDjWovz/GUve1nz67/+6wdkIWTMZGTr+JM1BcbKp7gyZDj79sdC5d4WtvHlU3RC3RCnmL195ZyET9EnJD+kn7UxxCEkh7ayrWXsx8AY/ntx9vrH7Jqqn+qLgAiIQKoE/Dxg908538fs/cEf/MHmy1/+ctMFGC9/+cubX/u1X4vVVXnGBPzAytgUqS4CIiACIrAjgbHzydve9rbmz//8z5voXSQ2qrH2nPIN+xQu0qefXoxPf6v0j8bs2mscpk9MGoqACIjAIYGtz6MHl0he8YpXNL/6q78qn4iACIiACIiACIjASQTe8Y53NF/84hebJyCF0cZJEtVYBERABERABESgegKHd5FUj0MAREAEREAEREAEliCAxZ3YugxGm8K43NUfERABERABERABETiBwGVMcbXI8zd+4zdOkKamIiACIiACIiACtRP46Z/+6atlFwd3kdy7d692NrJfBERABERABERgBgHe3Rd9VPj5+fkMsWoiAiIgAiIgAiJQM4GHDx925u/2qPCa4ct2ERABERABESidQDSDETJct7GGqNRVpsxWXf6WtSIgAiIwl8CkAAOdaIKZizr/dgow8/ehLBABERCBrQjoEslWpNWPCIiACIiACFREQAFGRc6WqSIgAiIgAiKwFQEXYGzVrfoRAREQAREQAREomQCf3XnxWyTdr7ZrEwEREAEREAEREIHTCFy7jCkWeVQ4fgoWrzt37hxpxWNb/bw2+onpcqTcTgXg1KcjbUCdvm1ITl9bHRMBERABERCBVQj4R4Wf0gkfroG/DCgoj8c4aZ7Sz5i26Ad94hUKeMbIWLMO+JBJrB/aEDvO8lA9yn/w4EHnC20iIAIiIAIisCWBybepjlXu5s2buqW1BxaDi9u3b58cAIUCqK0Duh5TdUgEREAERKBCAqveRdJ3OSR2CYOZD7a16X++Rx1/3JZ5P7Iv1AlNxqiPY4ARqsN+IYeyfB92n3X65LE+ZfbJwzEv0/cHHVEHm7XRy6ctMd2G9NBxERABERABERhDgAEG6p6/9rWvPb916xZ+hKR9ntbxFitnTRzH6/r1691fW5/77Tf2rjr++jooZ1mb2u/qWVk81veX7ax8Xx8yuaG+P459W4c6WF3Yz5WgyzcxeZDJzfdn+7Ly+myw8qx+aB+TH2MZ69/bZvv0x7QvAiIgAiIgAiDQfrHtQCCmwLzxRPvf1Q+T4P0pG7IA7aTViYhlDObIx3oCvOzmy/gt3taBLq2tDS5HYMPlA/9tn3VQj3W87miHepADG0MbLg1hozyrL9uwD9TD8ZDOXjb6tDb443bfykc7yudlE/SJOtRtaC1IX186JgIiIAIiIAIhAqtdIuGkdvfu3VC/k8swYWOCtsFLqCw0Wdtgwgc+nFxtMME6XhbKUeYDDxpj6/M99F1iImeftm+v31iokIEXdEOgYQOSsTJUTwREQAREQAT6CKwWYKBTZgz4zb1Pka2OURebxUDfyDwABl4MOvw3+yE7YhP+ULutbGfgBLsQ+NHemN5b6aV+REAEREAEyiPAAONJmManbi1lJr4lYyLzE/VS8k+Rw8mWMhgMWZlLBQapTODUg36hrQiukGVZyt5T/KK2IiACIiACZRBY/Umefs3EHtjsBM9LNphM7YSK9/bSwRw9rbzQ5ZI5MpdqA31oH/76SyOpBEFL2Ss5IiACIiAC+xLgkzwvMxjLPyocky6yBT6LgXI72QODr7MUGnxDR4bCrgfBJIuNurEOyrxeY/WwtlIeJm7aFcqSjJU9tx76Z/BAu7BvgyHI9vtz+1M7ERABERABEQCBgzUYi18juWTMydwi52SMMkzAePnLFku6yAYXNquCyZb9og7r9d0p0qeXzQRAlg0uQhz6ZJ1yjMEM+ke/4G3LlrD1FP3UVgREQAREoHACl9dIkLo4f93rXte89a1vbe7duxe8swDRyNJ3HPAbNibANb9Fc+KP9WEDg1idKUNhaXlT+mbdGNtY+VAfa/h/qE8dFwEREAERyIsAEwbf8R3f0fz+7/9+95Nn59h5y1vesmmAkRe2urVVgFG3/2W9CIiACIwhwADj9a9/ffN7v/d7zRNoxOslYwSojgiIgAiIgAiIgAjECByswVCAEcOkchEQAREQAREQgSkEFGBMoaW6IiACIiACIiACowgowBiFSZVEQAREQAREQASmEHABxpSmqisCIiACIiACIiACYQJ8kufoB21pnUYYpEpFQAREQAREQAQeE2C8cBFgdHerxreln4ER70lHREAEREAEREAELAE8w4gPbMTzmrZ8eOMcT/BR4bpNdQ49tREBERABERCBDQgguMBPUPDp0Hga8xIPhVxT9YMMxlqPCl/TAMkWAREQAREQgdIJ8OnQvJKAyZvBRrK2Xy7CGL0GI1lDpJgIiIAIiIAIFEqAvyllzVvz97uWwDjpNlVU1isPBksMDskQAREQARFIhwDXXXDiTn4NxmUGY/QaDKRn9EqbQTofB2kiAiIgAiIAAp/61KeaD3zgA1EYv/iLv9i86U1vunr9+I//eLAuggr+IjjWZKS8TcpgpGyIdBMBERABERCBVAk8++yzzc/93M91gYbfEEzg9dnPfvbqhYCDEzTq+ztIbt++3Ynh2gwvM4V9BRgpeEE6iIAIiIAIFEkAWQsbKISMRDDB7X3ve1/zxje+8WqfmQwEErhzBBkMvkellO8kcQFGyHSViYAIiIAIiIAITCHw3d/93V1ggaxF32YvhfzCL/xCg9dnPvOZqyYMPrjeAkEGL40wi9Enf89jfJLn6DUYIWUZUYWOoQzHCYd1mO5ZI72ztmxEjH0PObH2UhfPxjMZYujba18EREAERGA7Apj4kV2IbThuMw8IML7ru76re/VtuCzCzcq37xlkYP0j1l/ghfd+Xu3rZ49jNnNz/j+9+c3nt27dOm8VaXU/3vrK29tljhu0JS2ITh5eeM+tjby6MvxdeltLNu3wf61d4MDjsIv7tg7KbR0yIkPu+zZjOcX8NLa96omACIiACBwTwLn1ve9979GBn//5n+/O6aFjqPz+97//6pz/3HPPHbS384k90AYYV23a4OWoz5QL2i/WnXqIKWDfSRmMVsCoLfUVr31G8DoX7jtuuXUvpqdCduH42A2yUd9nc/z+WHmqJwIiIAIisDwBnKc//vGPH2QykLn4iZ/4iaYNLhq7luLU3m02xGY5TpW7ZfvFFnniiWIQxlcsdRMrZzsaj3ooY30ex2Rs++jrExO0PW4nbMrncR5jG14GwXFsfGKal8EHnUAe2rAe241xpu0T7xmwcEEPZFi7LZcx8lVHBERABERgGQI2yFgruFhG0/2lcB7sMhjtbHySRphs+a0ek6PdcAwvXz6nQ9tHX5+Y7FHXZxoQDEAPtG0vQ3QqYFK3wQPa4oU6DHIYTFideQxt+R7H2SfrQr4NdmJ2I5BgP+wbcqkL9CVHq29MnspFQAREQASWJcAgY43MRUhTm80IHU+2bOqDtvoM8ZOsr8sJEZPo3M334fetXDtBc9JGferBrAODAT9h+8sWQ3rb4MAGG9QJOvDVZz/bxvqDngg0Ysf7ZOuYCIiACIjAaQSQubAbzsV8+Xkk1BPqhuYI1sWxMZdF0Jftu09mSI+1yw4yGFPS+l6x0Ld7Xwf7mMzxbXyME0Ltp5TZCdi+52UMZhWYVbE6WXvYNqRzqCykIwIC1OUrVCdWhkEDfaA3dU5tIMV0V7kIiIAIlETAXhaBXViT8W3f9m3d+ZnnaDvfhGxHPcw7mHN9dgJlOGYDDNzK6s/5mEswH7BfKzPU5x5liwUYY5UnJE7yY9vNqWcn/1AgwMWa/Buqg345WEI6Mzjxzp+jb18b6AY9mW2BLmv32aePjomACIhAbQRscIHAghveY5EnN5yfY/MJ6tgvsPZ2VPtMDJ/B8MsL7M0FmBdCX4r39s/mAQYM5rqHkPGYNOEYDzNUd6iMToZMBgcIFjhJc4JGGUD0DQi2QT20w4vw4NihiHVI177j1I/9+vUdfW11TAREQAREYBkCoTUXOB9zTYYNMji/8K/VwM41NqgI3YViAxDKsu3RP8pRxiAj9GV4GQLTpLgAY1rjubVDkzGDDgQWiMzGXnIZ0gGyGKywDzgD8pmigjOGggS04cSOdpSJdn2ByZB+fcfRB/rloKK+7JvlfTJ0TAREQARE4HQCyF7wVlR7zud8hiDDZx3Qa2x+sF8U/ZoOaovgwh6jLCvTzgN2bo31ezqJ8RJ43whuHzm/efNG88/+2Rube/fudRGZ3xCNhMp9vVP2ASUUgKwhkw6Y2t/cdlNt8CzsgOnTeQs/TbVF9UVABESgJAKxecBmtlGH+7Ddzp8IDPyX3z/4gz+4ClJsYOFl4PzPLIWVif546YSZjT2Y80v7P//nN9srFmfNk1DCGrGHUuizb+Kcq1NMZqx8qJ+57Ybk+uO+H7/v62tfBERABERgGwKh87HNJgxlmNGeAQa/TPJyiV9GgEx5Kpc9ptBlTHERYDSnPQdjSsequy6BFILFdS2UdBEQARE4JDCYYb/M2c+d6frkI2BgEDB0yf1Uv+USbLSP3uxMTSaDcSp4tb8g0PdBECMREAERKI3AlC9VxwsAhmn0BSW27zmXJkLZkD6NcslokMsTMGaKg/qM1zEREAEREAERKJ2AX2OBSxtDl0bIxK6p85z8sVwyFt6Og0skYx4VriDEI9T+XgQwFkvP1MjGvUbXsv3Kj8vyTEGaXVAJfWLnoli2wQYRPoPhAwzay7sr7eUYrt9AnT6ZuzC7vCQ1+hJJDOIuyqvTIIFagsAaxqJsDA7x7Arlx+xcNqiwf9CVz1wgCMAL5ayL96zXl5Wwx6xcBiL4ax9XwMCCZVDeBy2DBq1Q4SCDUcvEtAJHiRQBERABEaiEgM8w2ImdCFDHZhdQjnq+rn0eBi6xMBgJzccMNhA8MDOCYMTXTeVhjNRLazAq+WCUZKb/UJVkG22RjWV4VX6c7sf3tKsqQ6/pkpZv4QOMoR78baes7xeEInCIBQdeRkwHBB426zGk25rHOe67B219//d/X/P61//TXR+0taaxtciGU2tIydbiT9kpAjUSQHAR2u5HbgEZdd7joyVDggfKukky8ADKgWYHh5nVQOFQEGCzFX2XO6xM1OurO0XXU+ryQVtvecv3N//pP/1n3aZ6Cky1FQEREAEREIEhAlMCgKEAhH1NkTmk39LHmcEYvcgzpEAsYgKg2LGQnFgZAY4FHpMzVA75vGbGvmx0aNtTJxxn3ZC9oYjVtoEca9cceUN2lXp81DeWzI2XjZk78FJ9+bEMP8qKaQQOA4yZT/K0Ty+b1v3FSle0PzX1NLVfX58gcP2KC3GgEyb/0GpfBgZ2wQ4W5/hralyww9uLEEDYMrRnEGKPQT/Uw3U3G6R5ed6Omvb3HjNbsJaNW1Bevw/5cX3G6iE9AnyS5yKLPO1EifecGDFJ85s5yvDebzYLYOujHtrYSZay7Td/ykMZ6oaO+T79PoIDyOZiGrzHCycHvtgGfTC4sMdswBHSgbajDd8zgGF9K8/KCMnzNmhfBERABERABFIgwC/u0OX8B//lvzy/desWltC0c9zx1leOY+3EfPCCBJThmH+xri1H/VA9lLWTf6dQmwno6vCv1cm3RZ2xG2WiPvqycinDl2MfdnBj/3af9ltdfBseC7X39oXkeRtDuvs6JezXYKdsLGGkXpzXyrAkbsXSY/Xd7TQUesU0GNU/5raZr1HyY8pVVt5+ee4sRkwBbietwWgFdBvT99xv5fNtd88uMwKoh/f4Ro7y0CWSdiLtMhH8lk9BqMvLEDgGWfabPY+hzLe9UibwBu2QfbCXSnw1HEc9bqHsAo9Dd9hmMy9sR7vIy8pBm9DWJy9Uv4YyO75KtVc2luHZ2vwYuwNkijdjd4tMkaG6+xLgfNpdIhnzqPA+dREU2Jety4mWf4cm/9DEzDYMBDhBM1hBfzxmy/p0xjHURTtM7jgRIEhAIGMnfupjyygXxxh8MLBhe+qMfWsz2oAVNh+YeX3Rbkieb6N9ERABERABEdiVwNRHhfcpGwoK+urPPYYAINQXggMGFpiQMXGP+ebAiZ/BA/4iYLABAeSFsguM0JhxgU1s5wMHyGUAgr+wgZkTtkE/sS0kz+oYa1dquVbml+FZ+VF+LIOArPAEDjIYdkGGrzhmnxOonUjHtBtbh0EFJlW+t5dbqD+OMxgYMwFTlp38oZMPYmL7vPxj7Udgwxdk2UtECF7YF96zL6szjmNjMOXl4dgY2zohhW5jgsfcTZeNuXvwQn/5sQw/yoppBGxMcf7Od77jpEWebdfdgg6+2snxnIsS8Z4bjreTabeLctbHvn2Pfd/e1rdybFvKYB9dRwNbn1zqgL928/Z63Vl3SE/PxsoNqe3l+To4rk0EREAE9iIQWpw5tWyq7qPOezMXeLYTUzc3aRtHoP3y21VETAFu3VNQ3/WudzZPP/3azR8Vzm/iPkPQ6hTdbBbDVpojy7efokdUwYEDfXr2HRsQ2x2uIeVci501+FI2jvlUp1/H+nGPRZ6jxtHOjwpP34vLaMglBf/qX72r+Y//8Vf2fVT4nAk91iZWPhbbqe2X6GcrHcbqmmq9NkBOVbXF9JKNi6HcVZD8uCt+db4TAV4iWeQ21Z1sULciIAIiIAI7EJiSqfgfu0T58fb/Xdvwi8IbjvsfXfI7o2uq4iUBBhjdbap8rKfoiEAOBDh4c9B1ro6ycS65tNrV4MdPzPypibQ8JW2WJLDoo8KXVEyyRGCIgNLOQ4TyOC4/5uGnIS3bRZxDVXS8MgKHGYwTFsBUxk10ZoSnAAAgAElEQVTmioAIiIAIiIAI9BA4CDBOfZJnTz86JAKLE6gh7SwbFx82uwiswY+6RLLL0Eq708ukxcUaDGUw0naWtDsgoNR6GQNCfizDj7pEUoYfl7SCgbXuIlmSqmSJgAiIQMUEPnH+D8fWb3m3yHHvKtmBgFuDsYMG6lIEZhKoIe0sG2cOjsSa1eBHXSJJbNAloA4viugSSQLOkArTCCi1Po1XqrXlx1Q9M00vXSKZxquG2idfIsFjrfnDXHgCJd/zx7wAEe9rejqlZcBBRB6WUWiAxbhRpm0TKgvJVJkIiIAIiIAIbE3g5AADCod+YhyBR6h8awO37g+TPn8hlQEAIePXUnEMr9C3NrZFPf9z85Rp7ak9wADXEMetfb5mf7JxTbrbya7Bj7hEoizGdmMqh5449y3yJE8bUMSCC0yK+BbPH/QiJGZC7Dd4HMM+y/q+/bM95drMylaOQJ8+EKA++Nl1vMdfbCy3uqEtf9YdEyfrsk77a65XPwFf+sQ6xmc1MJCNY0ZC+nVq8KOCi/TH4dYaLvYkT0yM2DBxcvJkGY1CNMMJ+ObNm1eXTRBAYJ/f7vEeZdhsGYIW7Icut6AMx9AW/fMv+97iL3RGUODtRt/UmX99gMF9BlGQ5euyjGy2sEl9iIAIiIAIiMAcAocZjBOeg2EnQztZUilOivgWjuO8DMDj2EeUz0jfT8D2WCw7wm/8CC6weRlzAI1tg76gl5/8yYXl/riXz4CKwZI9TrtjQZaXVfo+B2/JdsrGMrxbgx/tXSTvvvaSxr/K8GTcCswBfed3HMd8wFdfXRyzdWNz2RSZcc3XO8Jx/yS6OPVDwPQ+ZPlv8QTEyZ8m0SmAeWr/cAozJP7ywnoILyTDLl4GsbbBLpRDrzH2kSG4MJMDuxCYQRY2/I0FWWvbmZL8GtLOsjGlETdflxr8WPslEs5tDAzsaLFzE8uZkfdjIzRPcH6xQQnnCNtPTOb8kXtaS9rSrcE49VHhnPhgJCdDrx4mSvuiM9CGawx8m7H7NspjoDG27RL1mHXg5M8Bh0GBQUS70Zfnw33/N6QX68Si2lAblYmACIiACCxPAOf3UFDAnnCetvMRvnDaL8C2rQ0g8GXT1vNzmv2y7i/N+/lleatHSlzyUeHWKG+gnRTxHiAJiBMyy0eqflSN8ugU66yjygsX8BIO/jJ7w8gUA4j2UyfsY+B5mzGIWA4VLSu0tYOVMmOm2IFb4nvYVKJd8Cftko3XroZ3zr4u2Y+8NIK/9jLJUu+n+j12Ply6HOdf6OYnft+P/SKIOQHncbxCWXYrC+1QD19MuXH+sDIhB+Uo49yTSoabvusyGCcswejstxOen/wAAMbzUgEAEDD/WmdNBUTwhA19hhxPp639FzrBHgRA1m70i30OFg4k1gMvBhi0h0GUHXQx/W3qrcT3DOpov2y8IJEbB/nx8c+c5+Y7jDheGsFfe5lkqfdTmcTOh0uX49yMc7RfDjC2Hz9H2nY2+LD1OFfYAINzH9qH6o7VZ416jCnwNeH8ve/9seZlL3t5c+/eveDzBRAAWGfPUYhgPNxY+Zw+Um0DG73dIV1j9cYyWsJPIb1UJgIiIAKWwHseJ5gWB3P/cdw1Svao894z8xW+9jvtJNlmIEIbJnl+obXr5VAX522bXWe2gWV+3R3a2C/K2LdZHOiAeYRfwq1Oob5C+q5dBt1g1/ve997m4x//v5rLDMZ8+GMVBpjQJBsrHys3h3ohu0N6x+rVwCjEI1ZmP3SxOrmXy8bcPXihfw1+1G+RhMcqztvMODOD74ML39LPAXOzJF7u1vsc90+i4xo+BFsDVn/rEYh9m1ivx+0ly8btma/RY6l+tFmG++0lkuS377vQ8Nq/2VZTBhS+V3zTR0bDXubwdUL7U5cQhGRsUXYYYLSLdLSJgAiIgAiIQMkEzv/36dbNDUpsNsJePkE5AgVkNXzGYkg7ZDRyCDIWe5LnEBAdF4GlCdSQcZONS4+afeTJj/twT6FXBgJctE+dbNaC6+t4zO/nEEyEWHPcb7YGI6SEykRgDoFS086WhWycMzLSayM/pueT1DSyWQwfYFBXrsWI1bXtpmZF1uChAGMNqpIpAiIgAiIgAoYAshA2AIjdXoomNmMRqmeDB3vcPppBAYaGnwicQEBp5xPgJdRUfkzIGSeoUoMf5+CxzyzCYk9wwssGEQwGbF3Ws0EDgwnUZzYDcliX+oUe4jVH91PbHGQwTn7S1qnaqL0ITCCgtPMEWAlXlR8Tds4E1Wrw4wQcV1URDMQmfAQJlltfXf9wxb7LKFPvSplj16g2l0/a6h609ZM/+YHmpS/96lUftDVKKVU6iQCiRn3YT0KoxiIgApkRGHXe+z/m3ymJu0hOPa8iKGBggGCi7zKGzVb01ZsicyuXIquC4Olf/+ufbP7tv/0/m9nPwYBxFgTf4y9B4n0foK2MVj9lERh1QsncZNmYuQMv1Zcfy/DjqVYMBRVW/tgsxBSZp+o/tf3BJRLuTBWCaMUvYEFwwfKp8nKvj4HhBweDLF/ubWU9DBoGaKiD9xxIQzK8zFL3T/02kQMX2ZiDl4Z1lB+HGalGeQQWCTCIxS5ase8tNk6gdvLEcezjGF72mC3DBBubXFnP9+X7WduF0AOLcuzCHEDmPh8VG9KDbXEM/OzT3/gDaDg258Esof5UJgIiIAIiIAJrEVgswOCKVkzonNT989PtRIsJk5dNMLFinxMz3jOQsGWYdPsmVxxj33ayXguel4u+bWCB49QHi3zwnot9WG5loC2Y4Ri+8bAuWWCRD46hTiyA8zqVvD8345YTE9mYk7fiusqPcTY6Ui6BwwDjhEeFM1jABMjJk2XANzRJcjUtU4l+AkY5j4UmV8q3gQn6tTqs7Ub0jaDAB1ZWD8vJ6mOZQQ5erIu/kLulLWuzWkK+0s5LUNxfhvy4vw+W0KAGPy7BqSYZfFT4kzD61Cib374hy0+ynED9j76gnJPpqf2jXxt8xG4NWsPBsAN9468NBPgeNtLWvv5tBgTv8aGFDMphP55vn0wdEwERqJfAlJ9Un/oT6fVSleVjCCx2iQSdYRLEJIsXJ0SvBNL89sXJE21QfkoUzICCfTOb4XVYYx+BEy+DUD6DKpTDPsAOZV+sPszk8J5nawPeox8byK1hSy4ylwhIU7dVNqbuoXH6yY/jOKlWWQQWDzCIxwcY9hs43nOyRH1Ouiyfi5iT8dAkPlf+UDtkHOxiTGZroBcCJwZWkBPjw/LQccj3QcyQTiUfPyUYzYWLbMzFU/16yo/9fHS0TAKbBRiYZPHNG5Mkv8kz48C/KOclgrlBAi8dbHl5BEMDJxC+qANPKrCLAQODIOzzcorNUsB+e5mF9SwP1LdtyhyaskoEREAERCBnAgwwujUYc9Z4YgK00bl97ydCTJy8bGC/oaMe9205VPKRv9/38NEek3FKEzCCHWY3oK8NfmzggAwHsh7MfCBQgT2WiV2jkZKN3g9b7GPwDo2HLfRYsw/ZuCbd7WTLj9uxVk8JEbh8cOoiizzHmOUDCLaJlY+RaWVgwt57ASSDKOrFQMtmJnAMNtsJkvu+npc3hUnJdUsPLuA72VjGCJYfy/CjrJhG4CCDwZ1pItKqjexAqt/sxwZRY+ulRV7aiIAIiEDiBJ5JXL/C1CsqwNA3/cJG54A5SjsPAMrksPy4jKNit6PeD1z7fk9zvkynRkoNflwcWuECGWA8ATtLyGAU7i+ZZwgo7VzGcJAf5ccyCMgKT+AwwAhEur6B9kVABERABERABERgiACf5KkMxhApHU+OQA0ZN9mY3LCbpZD8OAubGmVOQJdIMndgzeortV6G9+VH+bEMArLCE1CA4YloXwREQAREQARE4GQCCjBORigBexFQ2nkv8sv2Kz8uy3MvaTX4cS+2ufbLMXH5oK1czZDeNRJQar0Mr8uP6/pxjVtSQxrX4MeQ3SqLE2gfttxt3SLP9j7VeM2eI3gwVOjFJnwUuH+AlG1j6yLq4SvVh2b14Ag+6CvGwMthvZjdKMdLz/zw5LQvAiIgAiKQFIHLmOKkR4Xz9zT4iG7766iYCPFiGSdIQPA/aIZj/K0NyMJx/vhXLhMqgibozeAJdjJNRJuwH4r2bT3YjZetBwb2t0h8wJbUwNpAmRjHDbrerAvZuBnqVTuSH1fFK+GJEuCctshtqgwmOCn6AAIMGCiEAgaWoT3eD8mx3+RRP/atfwv26BswQzajf/7MOn7QDFtMV1/PcsKPoO39OytbsBzbRyhIG9s2l3qyMRdP9espP/bz0dEyCSwaYAwh4jd41AsFGGyPb+Y22Ah9OFGHv1CKuph8+2QO6XbqcQQM0DP0M/EoG5ttoA38y3YMSGKByan6q70IiIAIiIAILElg0QADwviCkvy2ToXtJIsJ1H8b56SKLAACBsqKBQ6czPnz5rF6SwKbIwtBAW2nrqFAAfbAdtiN4Mnywb7nOUeXktpw8JZkk7dFNnoiee7Lj3n6bYrWmH9C53Urg3MB5gO8YnPW2HpoT1n4O9T/FHuWqHsYYJz4qHBMkHhxYvTwOMminOsUvBHIAuCFydQHEL6uhRnKHPj6e+7DZsIOBQo4zqCCmRAw4kCj7p7pnjbt3Xcos7W3Tkv3LxuXJrqPPPlxH+5b9oovj1wz6Pvl+R/HcV7nK5R55xdMX88HD5CJ9rYe5KcUzPJR4eBx/uEPfej81q1b+Jm99vNwvPWV+2PYbwONTgj+Yr+dWLu/3Mcx7LNtGySc42U37LPtwYF2h/KsDF9n6/2QvtbumD6hdmRFXrSTfyE3tJFn6JjKREAERKBEAqPOe/9PO7fNfMXk89zdd162527Oc7aM/rCycN73sq3fvEw7T3Du3cvPbfDTdY2YAnqedBdJK2D01hoeXQjJzAaEMdvBuya4bzvi5YbWCV3kiAjPR3mjFVuxIvW0KSzawwwF9q0NzFT4tBfqox5sDjFZ0YzkRCNSb8dwcnotqZBsXJLmfrLkx/3Yr9Uzzr/IHgxtNuuM7DXP2/7cDjmc7/DezgH20jrmACsTcwHnPeo0Rq8hvZc4fniJZOZzMGKKhIwkCAQafiM0QAZQQg1dUqAcC9c6x8vea98OBOjHF8vBCO8xMGAL9plqAyPaSf1rDyqsH0sPLmCrbNzrk7tsv/LjsjxTkIZzMc7RobnM6mfnAHv+xnscs8fZDnMBN9uGdW0bO0eE6u7Jyl6uOf83H/3IrEska6RhkP6PXQJYo7+UZJ5qdzugUjJHuoiACIjA6gRGnfdmXh7BZZU++fZShj9/20sXOMZ9/LV18R594AV5dmM5dbAybb0+Gas7wHTQBkDdHmIK6PxE+19Si0MQidX6bb1WuzEGp2wpLWaaoveUurJxCq1068qP6fpmbc1sJp+LMtEns9U+S41jfg4YypKsbcNc+Rz33RqMNsKYK0ftRGBzAm2AvHmfW3coG7cmvk5/8uM6XCdL/RcXLa79k8ktF2vQZhk6WVwCwLV3UzoILT+Y0n6zuks8KnwzZdWRCIiACIiACJxI4PwL0wUsEZTYQBOBBoOM0DqMPg2R0cghyFh1kWcfIB0TgVMJKO18KsE02suPafjhVC1q8OMcRrHLG30LMn3AkUMwEWKjACNERWVZEFDaOQs3DSopPw4iyqJCDX6c4wi/noIybBDh1xz6AINtGKzEghMvc46+S7Y5DDBOfJLnkopJlgiIgAiIgAjkTsAGA3ZBZ+z2UthrMxahejGZ9lENscBmS558kmdyd5FsCUF95UmghpSsbMxzbHqt5UdPpJ59O9EjAMBYwItBhL2EwgWgoMN6NmhgsAGZbAc5rEuq9jkae5I+zGDoLpI9faG+JxKoISUrGycOikSry4+JOmYjtWL+R5DgL2vEggMbfEDtvssooVtfNzL1oBsGGJs9KnwPI9WnCIiACIiACKxFABP60KSOIANBAQMDv+6CulGWzVbELneMlbmW3UNyXYAxVD18PGY8QQIU39uoy7azdW1KyD4KPNy7SmslgMEb+2ZQChPZWIYn5ccy/HiqFbGgIiR3KGBhmykyQ/2sWcaLIietwcA1ILsohfsMIBA8sMxCYxnb4hiDC15fwn4sgFkTzBKyQwMEZWMGBOuFZEA3MLXB2hL65iaj9OAC/pCNuY3KsL7yY5iLSssmsOgaDE56/DDZoIMYOSmGJkeWMe3TJwcTr598sR+Su4cLEUQgOLL6ADYDKC7MCelm63FRkK0HO/GAllRsDdmgMhEQAREQgboJHAQYaz8q3D59rG9yxOTM4wgyYtG/ncAx6dpLK3u5FXoAaii4gk645APbuGDHB0nU29dDG7xs8LGXjan0y8Gbij5r6CEb16C6vUz5cXvm6jEBApfXSE66REIz8CHiC2V+1au91IHJ0j/hjEEFJmd8Q6cslltcnJj5l8HF3pdToA8CotBKYJSN1Y828y8vq0C255rAMNpFhVjguYsyK3UqG1cCu7FY+XFj4OouCQIMrBcJMDCB4sXAwQcGnFxRjiAiNNkyY4FJlJM0n9ceImYzBaFJPdRmrzIEH7SZNoUyGLCDl1AQOPlALMRtL5vUrwiIgAiIgAiECCwaYGCyxMt/+/Yd++wDj7M99jGJYp9Bgw9WUIfHOOGGJmvf9977sIPQQ5kIHGdQwUwIgo0cbNuardLOWxNfpz/5cR2uW0utwY9bM829P46JJ2EIH+u5plF2HYbvh5kNlDNo6Lv0gUkXx2PrHbz8vfdhHzIX/uEqVi8GUjYIs2tN9rYhpf6Vdk7JG/N1KdGP77l2yOPdzXmDsvuhn2NoL3uWsJXoxxL8sqcNqz4qPDTxc+L0aX9A4OSKCRUTMS8jhL7pExrlpH55BPrSHmZnwAI248Uym41BOXmxfM/Bor5FQAREQAREYCyBwwzGzEeFhyJXW8bAgUrZY74t99lmaGLF8RwuIVgG/m4X2gA+CCgQLDHIQhmCKAYaYx1bQz0MXj9+SrNbNpbh0U+0mQtkMUreahirJftvDdtcgLFGF/NkDgUWkMqJOZQNmdfrcq0QENigALr2TYb2GNsysxHSqk9WqH6JZTUwkI1ljNzSgwt4qYaxWsZo3M6KRZ7kuZ26xz3h277PkBzXyrNkTJCVp2XSWgREQAREoHQCi1wi2QtSqYHFXjxz67eGlKxszG1UhvXVJZIwl61L/9srtu6x7v6yDjDqdp2sryElKxvzHOe/49R+Y7v+AmXvCazDuJ+niUda1zBWj4xWQS8BBhjdg7bWflR4ryY6KAIiIAIiIAIiUA6BqY8KR0SiV9oMODoZPWK/xPcchyXbKxsfP1AixzH82cvnXuAv32O84pJJaeM2lbF6BVZvdifAz+yoR4UjBbbHi5Rwt4hPw1l97PMy8N4eg4w9dN+rT9pLdpZbKe/JVjYeruDPzb8l+xGXRrDhL99j395Vkpu/YueWVPzI84H+7k/gMMAIPWVufx2vNAg9uMuqp0WfCTlLqoiACIiACFRNYNUneU4hi+DAPzvCt+fzLmLPvbABRuwWT5tm9fJL2a/BRviqBjtlYxmfSnt5pAyLjq2oYaweW62SPgKTLpH0CVriGJ5eiZd9QJV/WBX6iR0fynCgrU1HLqFzijJqsFG+THHkzdOphvFqL4/Mo5R+qxr8mL4X0tKQAQZWHJ3/+3/3S83//+h3mnv37u0yEQ9FwHYA99XFw7dsEALkqK8PQFqDT9qIgAiIwFIExpzj/9sJywBe086SmkPGeQtf9nGl4f/+9/+u+d9+6IebUYs8x4meX6vvR828Y2N1Q8GF1agvMJmveVota7ARxGuwUzam9dmaq438OJec2uVMgOP+4ufaZ/7Y2VIA+HsdWEvB9RQoC62nCNX1WYuQXj5QCdXJvawGG+GjGuyUjbl/Gi/0lx/L8KOsmEYgqQCDqseCipBpU+qG2qtMBERABERABERgeQIMMJK4RLK8eccSlao8ZpJriXyZq+cO9ZYf5ccyCMgKT+Agg1HDo8KVqvRDIN99+TJf31nN5Uf5sQQCvLQ/5lJ97E5IcsBx+9gF7IeWCvDxDmyHOmP634z35bKLJNZgbGa0OhIBERABERCBBQncvHmzkzZ02R7H7SMVfEAQyuhBtr+BAcEF+6QZkItHPaQStOsSyYIDLBVRoQGaim5L6lGDnbJxyRGznyz5cT/2a/eMAGGsfxEU9D2vyQYbuM0TQQU3BA52s8EF6tkHUIayHWtzCMk/DDBOuEc4JDy1MhqLv/aVmp6n6pNK9HqqHUPta7BTNg6NgjyOy495+GmKlpjEMY/4ib9PBoOC2NOorSxe/rCPZGAAYi+fMLOBMsrtC2L69Fv6WDKPCl/asJA8fMhDr1BdlYmACIiACIhAjAACDEzosWDBt2NwYLMSvg73bR2bjWBg4ddnsF2obqyPLcoPMxiPf0F4i7536WNsKmsX5RbqtAYbgaoGO2XjQh+KncXIjzs7YIXuETBgosdrTNDA7IRfc0HVbNAQU5eZiVjd9AKMC0uquU1VqcrY0M2vXL7Mz2chjeXHEJX8ymrw41yvcOKPPYHay/VrKMZmSbycvfcZWOsukr09of5FQAREQARWJfCfL6X/+IbZemQb+NscPnCYa2wqayyG9K8uwIDBpUfaNdiIgV2DnbJx6BSWx3H5MS0//cL5dH3mBiVc2MlLKr5nBCBTAw9kNHIIMqoLMEoPLjB4a7CxFjtr8KVs9FNOnvs1+PEUz/hnVlAWn3Hh10/Y/RyCiRAbBhjdGowanuQZgqAyERABERABEdiTgA8wQrpwLUasrl38OTUrEurv5LLLJ3lWs8iTEdXJ4BIWUIONwF+DnbIx4Q/aBNXkxwmwCquKzI5/2UWbOObvLLEZC3uMQYMNHuxx+xyNFAIMXSIpbDDDnFpSlTXYKRvL+IDKj2X4cW0rcJcJL6WEglIGEwgeuA4DwYivO+a22bVtgXzqdZHBKPxJnlsAVR8iIAIiIAIiMIcAAodYcOBvcY09CwOBh8+IzNFliTZ8kmc1t6kioir920QNNmLw12CnbFziNLe/DPlxfx+sqQEm9CmTeiw4gI6UZbMVscsdmMsgi/JQL1Z3TftjspnBuAwwYtXKKS89uICnarCxFjtr8KVsLOP8WoMft/TU2IAltaDCMrpc49lUs8hzywGivkRABERABESgVgKHazAYbhRKA8aGXqWZS6eWZpe3pwY7ZaP3ep778mOefpPWpxFwl0g2fH7qaXrPal1LCk92zhoeSTaqwZeyMcmhN1mpGvw4GUrlDarKYFTua5kvAiIgAiIgApsRqC7AUKpys7G1ekfy5eqIN+lAftwE8+qd1ODH1SEW1sFBgFHDo8JrSOPVYCM+hzXYKRvLOOPKj2X4UVZMJHC5rrOa52BMxKPqIiACmRN4T2Rp2f0Zv6iZOQqpLwKbEji8RFLBkzxrSOPVYCM+JTXYKRs3PR+u1pn8uBpaCU6YAJ/kWc1zMJSqTHg0TlRNvpwILNHq8mOijpmoVg1+nIik+uqHGYxIKrF6SgIgAiIgAiIgAiIwiUB1T/JUqnLS+Ei6snyZtHtGKyc/jkaVdMUa/Ji0AxJUjmOimkWeNaTxarARn6Ua7JSNCZ41Z6gkP86AtkKTRyvIlMg4AXeJRNdI4qh0RAREQAREQAREYCyB6jIYMLj0bxM12IgBXoOdsnHsqaxpUr4dVX4c70fVLIdAVQEGjeVfurG0gKM0e2IftxrslI0x7+dVLj/m5S9puwyBqgKMGj7kywwLSREBERABERCB0wgwwOieg1HDo8J99uI0fGm2rsHGi+Fa/pqhGnwpG9M8j0zVqgY/TmVSff3L+1T1oK2CRkItmZoa7JSNZXww5ccy/CgrphE4yGDwsZ7TRKi2CIiACIiACIiACBwScI8KLx9PDWm8GmzESK3BTtlYxjlJfizDj7JiGgE+yVMP2prGLenaNaRj4YAa7JSNxx+1lG9HPdb2okR+jJFReckEDi+RMNwo2WLZJgIiIAIiIAIisDqB6gIMpSpXH1ObdSBfboZ61Y7kx1Xxbia8Bj9uBrOQjqoLMJSqLGTktmbIl2X4Un6UH8sgICs8geoCDA9A+yIgAiIgAiJwKoGzs7Pmzp07UTE4fuPGjatXX10cs3XRNrRNkRlqv3YZA4xqFnnC4NK/MdVgIz4YNdgpG9c+BW4jX37chvOevdy8ebPrnoGB1QVlDx8+PFAP+3fv3m0ePHjQteHGSdlWhuzbt28fBDAILtgn61JmKnNcdRmMVMAfjLSFd2qwEchqsLNmG3G3SOh1/7xpQq+FP0aLiqvZj4uCTFAYsg2hoICqIhCwwQUChevXr19ZYoMEm9VAHdTlhmDEbradl2kDlj2RkUuXwSj9UeGxQVDDh3/PQaa+RUAERKA0AqGsRMhGe3nDzjUIJhg08FKHDSLYDv0wmEAbvKxMm9mgTj5bEtJrk7KDR4U3Zf+2A52Lv/a1CegNO4kFUhuqsElXNdgpGzcZSqt3Ij+ujnjzDjCZI8tgsxEhJWwwYI/3ZRls5sLWoywr02Y9QnVDOm1V5p7kWXaAAag1ZCtqsFG+3OoUsX4/NYxX2bj+ONq6B2YSMNnbgMDrgeNYZ4GX3Xzg4fe9HOwzMxGrm1yAcfhjZyGTVCYCIiACIiACIjCXACZ+n7Gwl0P8Mb8/lCWZq9fa7ap7VDhSlaV/m6jBRnwwarBTNq59CtxGvvy4DeehXv7LZYX/umOyHtkHu0DTZzaGbMDxZNZYDCjLS4PV3KZaenABf9dgYy121uBL2Thwls7kcE5+/B/aO5GmbksEJXZhJ/r3t6iO1QkZjRyCjKQDDER6vNYUSjHhmF/gYvfHOkv1REAE9iUQ+gEz3IqqTQRKIYA5jEEBAoTYOgrYi2P2MkkOwUTIT8kGGADsHyJigftIEMaNeciIUpWhYZBnmXyZp9+81vDju5uyowmNVe/1uvZtcKebCSwAACAASURBVOEfmEUSdn6LBR9ci2Hl2WDEtrPy9qLNAOMJKMCdvZSx/frgwutkF8j4h4z0ZTFySuN5m8fu12AjWNRgp2wcO+rTric/pu2fNbXDpO8zEJij7Mv3b+v7LD3q+i/bbN+3cNT3scV+khkMAo1dZ7LA7TUsGgPIfUHGFmDVhwiIgAiIgAj4ucg/kROEEDDghfmMX65DX/gpC3U5PyIY8XX7bpvd0iNJZjDoAO8YgomlgSzUWIrJO2JL2Fv1VYONYFmDnbJxq0/Nuv3Ij+vyTVm6z1706YrAIRYc+LtNYnMcAo/Y3NnX9xrHOO6TeVQ4Uz+AzPfe8CkO822VqvRE8t2XL/P1ndUcfgwt8izDugsrNFZL8uaxLZjQY5P6VN9Tls1WxOZCyEagwWAD9WJ1j7XeoOTyQRgXt6nu/KhwXquaG4EBrH+2+wYI1YUIiIAIiIAILEogFrD4TpILKoyCfFT45XMwvOrb7tsfdFmj5zlpSkafaJvLe7CDrjnpPIct/TmnLRjlwKdEG3nHyCfaLzR4j78cs7n4ZaqeJfrRf+72trEbRNqSIsAnee5+FwmjNd6GY9M+IOb3hyiG0kT4QEx9sR9+mLCf+nvql7qep7KlL0+Vk7JPS7PR3o7K9/irsXoxinPmkMpYHZobdHw7Agw6k3mSJ9ZXhG5RZbn9AHpMvA7ly7UvAiIgAiIgAiKwLQEGGLtnMKaYzSyHb2MDjFAGA/VpsG9b0n4NNsqX5YxYXiIpx6JjS2r4TNZg47FnVdJHIJkAA5dImGKzf6k87iph9sIufrHvx9xd0pcB6QOV07EabIQ/arCzBhtLf4qnxmpOZ0/puiSBZAKMKUbZ7ATuGoERNnqO3Uc8pQ/VFQEREAEREAERmE+A83IyazDGmmKfeGbbxJ7zzjowuPRvhTXYCH/WYGeuNsaeaxH6ATPYeF+/RTL21JdsvVzHarJAC1As+QAjFgwgi4Fj9u6SMfcNx+QV4MsrE2qwEcbWYKdsLOOTKT+W4UdZMY1A8gHGkDkpP2RkSHcdFwEREAEREIFSCTDA6O4iaZPOpdp5ZRcNLtnQGmzsRiuf4lKwM2VjGc6VH8vwo6yYSuAipkjiSZ5TVZ9TX6nKOdTSbCNfpumXqVrJj1OJpVk/Bz/+1+YbToD3lRPa1tk0mSd51olfVouACIiACIhAmQQOLpEojVeGk2vwIzxVg52yUZ/JXAjUMFZz8UUqenJMZHeb6lyAOaTx5trGdjXYCFtrsDN1G6fcjhob16nbGNN7SrlsnEJLdUshUFWAEYuwa/jwlzJgZYcIiIAIiEAeBDjnZvVbJHPRMpDAX/uaKy/VdrFAKlV95+pVg52yce7oSKud/JiWP6TNNgSqCjCAtIZsRQ02ypfbnCC26KWG8SobtxhJ6iM1AtUFGKk5QPqIgAiIgAiIQIkEqgswlKosZxjLl2X4Un6UH8sgICs8geoCDKUq/RDId1++zNd3VnP5UX4sg4Cs8AQYYHS3qWoTAREQgRCBJW5HDclVmQiIQPkELn+LpHxDlY4tx8fyZRm+lB/lxzIIyIoYgWoCDKVjY0Mgv3L5Mj+fhTSWH0NU8iurwY/5eSUNjasJMNLALS1EQAREQARKInB2dtbcuXMnahKO37hx4+rVVxfHbF20DW1TZIbab1WG31Q9f/GFzzcf/shHm3v37hX3vAikYRFh8+9WYPfopwYbwbUGO1Oxcc01GKnYuOZnVTauSfdC9hjG15qXzVfk2ld650X0j+3BgwddcGA3BAx3794N9u0zP5TjK9++ffsggEFwcfPmTV+t2/cyg5VWLHz48GFz/fr15gsvvtC88lWvbqrJYOwNfkWfXomuwcYUPkTy5TIEahivsnGZsZKiFAQPsaCA+trgAoECJl9uNpNh36MO6oZkoMwGF16mD3D25lZNgLE3aPUvAiIgAiKQPwFM4ggsYpkJWmiDBmQ3sG8vedj29j0vuaCNl2XbM7OBMgYuyCCktFVxm2osyizt28WYVGFKg2+uLjXYuYeNocsh98/nemm43R42Dmu1bA3ZuCzPFKTZLEHfhG6DAdsGgQEDCq6loF02c2HbUJaVaQMY1KUuXuaezKrIYCCQCL32BL9G36UFTDFGNdgpG2Pez6tcfszLX2O0ZSYCE7kNCHzbvuDD1rVBg5fBfRs8hOqEgpFQva3Lqggwtoaq/kRABERABEQgRKAvGPBrKOyajZCs1MuquEQCJyhVmfpQHK+ffDmeVco15ceUvTNetzz8+H0XBl37D+MNS7Dm2KxIKqpXE2AoVZnKkDtdD/nydIYpSJAfU/DC6Tpk5cfz/2W6wSsGJT5jMaQcMho5BRm6RDLkUR0XAREQAREQgYUI9K258MdyCiZCeKrJYOSRxgu5aHxZDTaCRg12rmlj6G4RcF3zjpHQKF7TxlB/e5TJxj2op9FnLNtggwifwfABBi3hWozY3SJ9MvekUU0GI6s03swRUYONQFODnbJx5ocgsWbyY2IO2VCd2IO0+rIS9pi/DRWq24DEHrfP0fBBy4YmH3VVTYBxZLkKREAEREAERGAlAnaiRwCAbBZe3OwtrvahWqxngwYGE5BpH6rVJ3MlsyaJrSbAsI6dRCijyjXYCHfUYKdszOiD16Oq/NgDp4JDNnCw5vIpnCxD4BB7poaX0XcZxWY1UsBbzY+dpQBbOohACgRSWYORAgvpkD8BBHFDl6KuNf/rfEPbu0iG5A8JR1DAwGAoCLDZir7LHVYm6vXVHdJvqeO4xIMMC3/srJpFnksBlBwREAEREAERmEJgSgAwFICw3ykyp+i6ZF1dIlmS5s6yakjHAnENdsrGnT9MC3UvPy4EUmKyJFBFBoMfcv9hPzXtlZrHS7MnxrcGO5ewMfVLIUvYGBsjqZTLxlQ8IT32IFBFgFHDh3yPwaM+RUAEREAERCBGQJdIYmQyLPcZmgxNGKVyDXbKxlFDIflK8mPyLpKCKxKoJsCoIYtRg434LNRgp2xc8ay3oWj5cUPY6io5AtUEGMmRl0IiIAIiIAIiUDCBagIMpSrLGcXyZRm+lB/lxzIIyIoYgWoCDKUqY0Mgv3L5Mj+fhTSWH0NU8iurwY/5eSUNjau4iyQN1NJCBNYh8KbHP29w0MFnztfpT1JFQAREYAyBajIYSseOGQ551JEv8/DTkJby4xChPI7X4Mc8PJGeltUEGDWk8WqwER+hGuyUjemdLOdoJD/OoaY2pRDQJZJSPCk7REAEREAEIgSeiZSPKf4PYyqpToBANRmMGtJ4NdiIMVyDnbIxcLbKsEh+zNBpUnkxAtUEGEpVLjZmdhckX+7ugkUUkB8Xwbi7kBr8uDvkTBWo4hJJ7FuEPhiZjlqpLQIiIAIikDyBKgIMBBIIMkoPKGqwEZ+oGuwM2Rj7ddRcb0cN2Zj8GXOigrJxIjBVL4qALpEU5M7SAyi6qgY7ZWMZH0z5sQw/yop5BKoJMObhUSsREAEREAEREIE5BKoJMGLrMOZAS7VNDTaCfQ12ysZUP2XT9JIfp/FS7bIIVBNgKFVZzsCVL8vwpfwoP5ZBQFbECFQTYMQAqFwEREAEREAERGB5AkndRXLnzp3m7Ozsykrs37hx48hq1MExbqhj948atAVazR2ikmdZLb58d3P4a2X3C/vxslr8WHqmpgY/5nmm3F/rJAIMBAw3b948ooGy27dvHwQPCCTu3r17UPfhw4ddWd8Hue/YUceZFtRgI1xTg50+uMh0SPaqXYMfZWPvENDBwgkkcYnEBhcIKPDi5oMJu496169fv6o7lMUo3JcyTwREQAREQASSIbB7BsNeEvHZCgYTqOMvgzx48ODq8glXaqN+LMioIY1Xg4345NRg5yeaa03pWYwa/Cgbk5nrFlcEc83Sl/THylzcmJUE7h5gWLvsegu899kL60xbF4GJD0Y8L6UqPZF892vwZenBBUZfDX6UjfmeZ2Kar3VJn1+Ubb+hZQIxvVIs3/0SCQIFfAjxskGDzUSwHGsttImACIiACIjAXgSGLunbL8JjL+nb+Q6X/fuWCexl95x+dw8wvNIIJhDJMZgY8w3ABibWuVZ2KDr0fee+X4ON8FENduISSelbDX6UjeWOYl7SR3Bg5ynOQTZowCV9f/nDBh/2PdqjLtpws7JyIppcgOGzFDZ4mAuWH3L8HftiX/YEkfp7DvLU9TyVLewsxUb8gBlswV/7HjYiyMCtqXjZgKMU20vyI8Z0yC+y8XGgHOIT4zanfO78MLfd0NzUd0mfffovxDZzMeaL81zdt2qXXICBDyRevDsEAceU6C3kdMqc8pcOQBu9P7xeLiYXI2IpDna9hX2/lHzJWdZf4pkmz60mTc5NuFQSylZwvvJflmP6+SAjVG+srFDbPcuSCzAIw0IfcsDQcci00fOewNfsuwYb5cs1R9C2smsYr7Jx2zG1RW92vkGQAR/zEoe9rBHTpS8z4b8g28cwxOSlXL77XSSI9qxzPGDAY/QG2KFIzjo81B4y7LeOlB1yim412ChfnjJC0mpbw3iVjamMuX9xoci1N52sUF9GHQFHDT4fC3H3DIYNCKzj7HtGcbHjoaBjLADVEwEREAERqITA+Wfab5sTXwYNvszyCzHmJQQTeNnMRexLbojwlLqh9qmXJRVgIFBAusmmnACQgYV1BpzMuoRsF8h48EpVeiL57suX+frOai4/yo+5EYhdusfcZNcN9tnVd0nfH8v9y/PuAQYc0ZdSsk/sRN3YNS7/FFDv4L4+fN1c92uwcWi85Oo7r3cNvpSN3ut57tfgxzGe8cFAbP2Ev6Rvvzj7AIP9xmSN0WvPOruvwaDxGKSAS8CAHkofoczX7bsmtidc9S0CJIBbUENbab+QGrJRZSJQCgHMP7xEgvecr2xgYC/p86FcmKM4T/lAxLKxx+y8FpoLc2CK0975iy98vvnwRz7a3Lt3rzebkINBXkekYRGQ8K8/XtJ+DTbCXznaOTXAyNHGqZ8l2TiVWJr19/bjmP6vNZ+dD69dGGqzNOivbwv9Tlaovs26I0CxTwj19XPJEiFAQoD1hRdfaF75qlc3SVwi8TDX2M/FQafYXoON4FODnbLxlE9COm3lx3R8sZQmfT6de0kfGYrYGsLYsoCl7FlTTjKXSNY0UrJFQAREQAREYCkCCDKWvqTPyyj2poZcL42QczUBxpg02lKDby85NdgItjXYKRv3+hQt26/8uCzPlKRh8h8bAEypW9KaQl0iSWnEnqhLX+ruRNFJNa/BTtmY1JCbrYz8OBudGhZAoIoMRmxRTg0f/gLGaFYmTF3MmZVxUlYEREAEJhCoIoPBQAJ/7WsCpyyqxgKpLJSfoGQNdsrGCQMi4aryY8LOkWqrE6giwADFGrIVNdgoX65+TtisgxrGq2zcbDipowQJVBNgJMheKomACIiACIhAsQSqCTCUqixnDMuXZfhSfpQfyyAgK2IEqgkwlKqMDYH8yuXL/HwW0lh+DFHJr6wGP+bnlTQ0ribASAO3tBABERABERCBOghUcZsqXIl0bOmRdg02puLLtW9HrcGXsrGMSaYGP5bhqe2tqCaDUXpwgaFTg4212FmDL2Xj9if8NXqswY9rcKtBZjUBRg3OlI0iIAIiIAIikAoBXSJJxRML6FFLqrIGO2XjAh+IBETIjwk4oVPhmVQUqUqPajIYNaTxarARn84a7JSNZZyH5ccy/Cgr5hGoJsCYh0etREAEREAEREAE5hCoIsBAmjL0mgMs5TawsYatBjtlYxkjWX4sw4+yYh6BKtZg1JCmhPtl57wPwVCr0C2p98+HWp12vAZfysbTxkgqrWvwYyqsc9OjigxGbk6RviIgAiIgAiKQO4FqAgylKnMfqo/1ly/L8KX8KD+WQUBWxAhUE2DUkMarwUYM5BrslI2xU1Ze5fJjXv6StssSqCbAWBabpImACIiACIiACPQRqCbAUDq2bxjkdUy+zMtfMW3lxxiZvMpr8GNeHklH22oCDKUq0xl0p2oiX55KMI328mMafjhVixr8eCqjWttXcZtqrc6V3dMIhG5HhYS1b0mdpqVqi4AIiEAeBKrJYNSQxqvBRnysarBTNuZxAh3SUn4cIqTjJROoJsCoIY1Xg434MNZgp2ws47QrP5bhx5AVZ2dnzY0bN65ed+7cCVVrxtZDY8iwMtE2502XSHL2nnQXAREQARHYnACCgIcPHx70i30EBDYoQMBw9+7do3oo88FnKNt18+bN5vbt213gkeNWTQYj5LwcHdancw02wv4a7JSNfSM9n2PyYz6+GqspJnsGF9evX+8CAPzFxiCDsmxwYevhuA0a7HvKDMkYq2Mq9arJYPhoMRUHLKlHDTaCVw12ysYlPxn7yZIf92O/Vs82aLDZCmQ17GaDhgcPHnSXPrAx6IQc1gnJRH1kMLChnpVn+0n5fRUBRuxbRA0f/pQH31666W6RvcirXxEohwAyEnbz6yViwQfaMaBAHRuYWJm23MvOhWIVl0gYSOCvfeXipLF6xgKpse1zqVeDnbIxl9HYr6f82M8nt6N+okcQAB/j5TMMfo1GzFYvM1RvrKxQ2z3LqggwALiGbEUNNsqXe54ulu27hvEqG5cdMylJQxbCTvzYHxNQ9mUm7DHYyrUdKdk9RZdqAowpUFRXBERABESgIAK/9ZKmwetauypg6stg8NkGXNLA+gobCPhMRkEUJ5tSxRoMUEFkWfq3iRpslC8nf8aTbVDDeJWNiQ2///d8ukL/87VgG3v7KAIPZi98EBJsfFnoMxZ9dXM8Vk0Go/TgAoOvBhtrsbMGX8rGHKeMY51r8COttgGBDw7srarHlB6X9AUg/liuay9obTUZjD6H65gIiIAIiIAIDBHwQYWt74MBBBy+DPVtEOHl+QCD8nNdi1FNgKFU5dBHJ5/jY32Z8+2oY23Mx2vHmsrGYyY5ltTgx5Bf8IwKZm9sYMBbTbEWI/Qci1DQQfn2mF3L4QORkD4plukSSYpemalTLanKGuyUjTM/BIk1kx8Tc8gC6mBRJzcEV3gxkEA5AwMbFPAOE67TQD37zIuQTPvwrVwXjlYTYCwwriRCBERABESgcgIIHPxDtojEBgoo8/usZxeIomyKzJzwY3ns+YsvfL758Ec+2ty7d6+4hYJM39WQxqvBRny4xtqpSyRpn4rG+jFtK/q1k439fJY4Oobxtd86oaf2LpJYJspmK/ouY+ASCi+jDGUjxso8waLVmuISD9aLfOHFF5pXvurVTTVrMGIDZDXSOwiuwUZgrcFO2bjDB2iFLuXHFaAmJHIoWKCqCD76AhBr0liZCWGIqqJLJFE0OiACIiACIiACIjCXQBUBBlJooddcaKm2g401bDXYKRvLGMnyYxl+lBXzCFRxiaSGNCXcX6udOa+1iH1sa/ClbIx5P6/yGvyYl0fS0baKDEY6uKWJCIiACIiACNRBoJoAQ6nKcga0fFmGL+VH+bEMArIiRqCaAKOGNF4NNmIg12CnbIydsvIqlx/z8pe0XZZANQHGstgkTQREQAREQAREoI9AUgEG7v/l/cL4G/vhF5TbemPuG1Y6tm8Y5HVMvszLXzFt5ccYmbzKa/BjXh5JR9sk7iJBwGCf5U48KMNTwWyggWDCPqMddfH0MJT1pSP7jqXjjtM0qcFGEKrBTtl42mchldbyYyqekB57EEgig2GDCwQU9jnvCB5sgGGDC9SzP2M7JpOxB2T1uRwB3JLqX/fPmyb0Wq5XSRIBERABEZhKIIkMBpW2PwCDYIGpN7xHkGEDCPyIDB+9ynoIPmJBBuqU/m2iBhsxVj7RXGve3bRRRcFbDb6UjWUM4Cz8+KgM1rlZsXsGw2Ynhp7VHqtrMx62jnVG6cEFbK3BRthZenBRiy9rGK+yMbcpUfouSWD3DAaCCv6k7VCAgcsl2kRABERABERABNInsHsGA4h4R4jFZS91xC57sL4NTGIZjBpWOtdgI3yOSySlbzX4UjaWMYpr8GMZntreiiQCDG82AgYu5sQizqHMhm/v9/kBwN+xL8qwH57U3zMdm7qep7LFJRIbZJRoL3xZol3wPe2SjY8D5Zx9nYof/Xlf+/sT2P0SiUdgP2h20aevF9sPBSOnXAe1bfX+gvoWHPgDZggm7HvcLYLtvlnkad9voRvHnvrabjyI+eEaK429x2PPzhkcJ/qbBoFkMhi4tGEHCtZlDF0aIcLYZRGLuIZBWION8GkNdsrGNE6Qp2ohP55KUO1zJpBEBgMBgn0WRizjgMsloYWesbtLrGNiMnN2nte9Bhthcw12ykY/uvPclx/z9Ju0XoZAEgGGDS5Cl0W4CBQZDdbFe2Y4QkHHMngkRQREQAREQAREYA6B3QMMf3nDPwYcRqEOXnZ9Ber5uvZ5GB4GUpWlf5uowUb4tQY7ZaP/BOe5Lz/m6TdpvQyB3ddg+ABjyCw+M8PXC2U+bJ3SgwvYWoONtdhZgy9loz+L5blfgx/z9Mz+Wu+ewbCXOsbgQBYDA5pZDbQZuxh0jHzVEQEREAEREAEROJ3A7gHGXBO4LmNse6Uqx5Lath5vQfW98nZUX459+TJEJb8y+TE/n4U0rsGPIbtVNkwg2wBj2LTDGjWk8WqwEV6twU7ZOPUTnmZ9+TFNvyyplc2mh7744rjNsqNOLOuOcrtsAPt27eGSem8hq5oAYwuY6kMEREAERKAeAggG7F2QsNwGBAgQ/M0IuOsRZT74RCbIb5A9tL7Qt0lpf/dFnlvBCDlvq7636qcGG8GyBjtl41afmnX7kR/X5bu3dB9ceH1scIFAAc9y4mazGPY96tg7In2A4vtIeb+aAMNHiyk7Za5uNdgINjXYKRvnfgrSaic/puWPJbVhUGCDBivfBg18MrW9/GEDB/uel1TsHZOxSypL2rOGrCoCDHyLCL3WACqZIiACIiAC5RNgUBCb/G0wYS+b2OyErQNi9pht4+vlQreKAIPfIvDXvnJx0lg9a0jHgkUNdsrGsaM+7XryY9r+masdJ38EBLFFmGOfMD0meBgra649a7WrZpGnUpVrDaFxcufcjhqTLF/GyORVLj/m5a+YtjX40dqOgAATPi6NxLIXMVYoR0DC7Adk2QDFByux39/qk5/SsWoCjJSgSxcREAEREIENCfyXy75+6/hOjala2N/Dmtq2tvrVBBhIVZYeaddgIz6gNdgpG8s4FcuPifnxe8+nK2SCEr+w01/ewL7PSvR16DMWfXVzPFZNgFF6cIHBV4ONtdhZgy9lY45TxrHONfjRW41LJKFbVFnex8QHJVa2D05yXXtBm6oJMPwA0b4IiIAIiIAIrEUgtn7CBhg+gxELPmK3wq6l+1Jyq7iLBLC0mnupIbO/HPlyfx8soYH8uATF/WXU4EdSxiUSZCf8i8dxVwmzF3YBqH3fl5Wwx2wbH4js7/VxGlSTwehLWY1DlX6tFGxc8m6RGPEU7IzptlS5bFyK5L5y5Md9+e/Zuw0KcNeIfyKnfeYFHqrFSy6hgG3O3Sp72s6+q8lgpABbOoiACIiACNRDwD6N01rtf18EwYgNOGzdmIwcKFaTwdBq7hyG4zgd5ctxnFKvJT+m7qFx+tXgxyESsUwVAgcc490lkBPLRqCcL9RD21wvjZBXNQFGbAAMDZycjtdgI/xRg52yMadPXlxX+THOpqYjU4KFWACSIy9dIsnRa9JZBERABERABBInUE2AEVo4k7hvJqtXg42AUoOdsnHy8E+ygfyYpFuk1EYEqrhEwg+5/7CXlr4szZ7YZ6AGO2VjzPt5lcuPeflL2i5LoIoAo4YP+bLDYpy00C2p92c8iXdcb6olAiIgAiKQEwFdIsnJWwO6+gzNQPVsD9dgp2zMdngeKC4/luFHWTGPQDUBRg1ZjBpsxDCvwU7ZOO+Ellor+TE1j0ifLQlUE2BsCVV9iYAIiIAIiEDtBKoJMJSqLGeoy5dl+FJ+lB/LICArYgSqCTCUqowNgfzK5cv8fBbSWH4MUcmvrAY/5ueVNDSuJsBIA7e0EAEREAEREIE6CFRxmypciXRs6ZH2WjaGbkcF071uSV3LzpQ+8rIxJW/M10V+nM9OLfMnUE2AUXpwgaFYg4212FmDL2Vj/hNINp/HR2Wwzs0KXSLJzWPSVwREQAREQAQyIFBNgKEV6xmMxpEqypcjQSVeTX5M3EEj1avBjyNRqJojUE2AoXRsOWNfvizDl/Kj/FgGAVkRI1DFGoxYhF3DCS7meJWLgAiIgAiIwJoEqshgMJDAX/taE+wesmOB1B66rNlnDXbKxjVH0Hay5cftWKun9AhUkcEA9hqyFafamNrtqLGPy6l2xuSmVC4bU/LGfF3kx/ns1DJ/AlVkMPJ3kywQAREQAREQgbwIVBNgKFWZ18Ds01a+7KOTzzH5MR9f9Wlagx/77NexOIFqAgylKuODILcj8mVuHgvrKz+GueRWWoMfc/NJKvpWE2CkAlx6iIAIiIAIiEANBKoJMGpI49VgIz6UNdgpG8s4/cqPZfhRVswjoLtI5nFLstXYVGUud4vEII+1M9Y+h3LZmIOXhnWUH4cZqUa5BKrJYJTrQlkmAiIgAiKwNYE7d+40N27cuHqdnZ0FVUC5rYd2sW2szFj71MqryWAgVVn6t4kabMQHqAY7ZWNqp8p5+siP87il3AoBw82bN49URNn169cbG2ggYLh79+5B3YcPH3Zlfj4KXU6DzNu3bzd9QcmRIgkVVJPB8M5MyAeLqVKDjYBVg52ycbGPxa6C5Mdd8a/SuQ0uEFAgAOCG4MEGGDa4QD3U52aDBvvey/QByipGrSS0mgBjJX4SKwIiIAIiUCEBBAwIJhAc2ECSwYINGh48eNDViwUfNoigTLThZmXlhLqKAAOpp9ArJ0eN0TWUYhvTLrc6NdgpG3MblWF95ccwl1xLbYCAdRV9W6yuzXjYOpBlj1n5vl5fvykdq2INRg1pSgwq2ZnSx82DZQAAETVJREFUR+s0XWrwpWw8bYyk0roGP5I1Jn1mFoYCDFwuGbONCR7GyhrT35Z1qggwtgSaUl+5346aEkvpIgIiIAIgEAos7CWMocsZaM9LIggurDwvG+sxcg0uwKqaAEOrucs5OciXZfhSfpQftyLwxhcuevpsc23xLhEUMAhAQOCDhMU7zEhgNQFGDWm8GmzEZ6sGO2VjRmfRHlXlxx44Oxx6Y3M+ude+oMSusZlzO2npwUgVizwnjyg1EAEREAEREIEIAVzasMEF7xKJVD8o7ltz4Y/lfHkERicZYABy33UsXrdC9IdXX116Vqu5xwz9POrIl3n4aUhL+XGIUB7Ha/Cj9QTmH/ssDGSpQpkI+8wL3577nMO47wMMlsdkpT5CkrxEQud5+ICJYMI/eCT2ZDQLX6nK1IfieP3ky/GsUq4pP6bsnfG61eBHS8MGF6HLIvaLL+ti3uIX4b6shD1mvziHApjxHtqvZlIBRih48Gj8k9EQ8dEp1om+nfZFQAREQARE4BQCPsPgv+xCNuowy86+UM/Xtc+8wCUWBiOhjJANNk7Rf+u2SQQYiM76ojpCsZDhEEZ1dAgcGHME6pQeaXsb709fz7T1+JvVn7dzlpDEG8nGxB00Uj35cSSoTKr5AGNIbRs42Lo+84G5DGU+CEEb+0TPof5SO55MgEEwfYGGda5NGVnH+MiRcksPLmBnDTbWYmcNvpSNqU0H8/SpwY8kMzVLjnkKfJjVgJzYl2DK5nFeapnnlTRaJRFgWOB4H4rigKsv+EgDp7QQAREQAREQgUMCU4KFWACSI9Mk7yKZCtJmM2IprNB1ran9pF6/BhvhgxrslI2pf9rG6Sc/juOkWmUSKCLAGHINP+T4O/ZFmfYEkfp7pipT1/NUtrBTNl5QzJmD/Pj4qZLy4+njeWge0PHtCSRxiWRJs202g3JPuUZo2+r9BVFxEAeNAY2BVMaADc6WnEsk63QCRWQwYpdFLJ4aBmENNsKnNdgpG08/uaUgQX5MwQvSYS8CWQUYsaeZxe4usVBPyWLs5Zyp/dZgI5jUYKdsnDr606wvP6bpF2m1DYGsAgx/twkR6e6SbQaLehEBERABERCBsQSyCjDs+grcyor0o01B2iejeQBKVXoi+e7Ll/n6zmouP8qPZRCQFTECWQUYMCL2VDP/ZDRvsFKVnki++/Jlvr6zmsuP8mMZBGRFjEByd5HwaWYxhac8GS0mQ+UiIAIiIAIiIALrEkguwBhr7pQno0Em0rGlf2OqwUb5cuwnJP16NYxX2Zj+OJSG6xHI7hLJXBSlBxfgUoONtdhZgy9l49yzWVrtavBjWsTz0SbbDEY+iKWpCIiACIjAngTecELnnz2hbe1Nq8lgaMV6OUNdvizDl/Kj/FgGAVkRI1BNgFFDGq8GGzGQa7BTNsZOWXmVy495+UvaLkugmgBjWWySJgIiIAIiIAIi0EegigADqdjQqw9MjsdqSDnDLzXYKRtz/AQe6yw/HjNRST0EqljkWUOaEkNWdpbzwa3Bl7KxjPFagx/L8NT2VlSRwdgeq3oUAREQAREQgboJVBNgKFVZzkCXL8vwpfwoP5ZBQFbECFQTYNSQxqvBRgzkGuyUjbFTVl7l8mNe/pK2yxKoJsDQt6VlB86e0uTLPekv17f8uBzLPSXV4Mc9+ebcdzUBRs5Oku4iIAIiIAIikBsBBRi5eUz6ioAIiIAIiEAGBBRgZOAkqSgCIiACIiACuRFQgJGbx6SvCIiACIiACGRAQAFGxElTFi5NqYvuptSfUjdiSrR4iuwpdWXjMfKp/I4lxEumyp5Sf626cWvCR6bosfb4m6pL2KLj0qlyp9SfUncqv2NL+kum6DKlbn+vOroHAQUYe1BXnyIgAiIgAiJQOAEFGIU7WOaJgAiIgAiIwB4EqvgtEqbZfLqthofg7DGo1KcIiIAIiIAIVBFgIJBAcKGAQgNeBERABERABLYhUEWAEctgDCH2GY+++lPqQs6U+mvVXVOPXGVPYZ2rjWvqLX7HZ4kpTKbUXdOPOcs+9kDT3D8PlY4r+8S1cfVU65iA1mAcM1GJCIiACIiACHQElPmePxCKz2BocMwfHGopAiIgAiIgAnMJKIMxl5zaiYAIiIAIiIAIRAkowIii0QEREAEREAEREIG5BIq/RHLnzp3m7Ozsig/2b9y4MZdXEu1gD16wJbT5Y7A3VjfUfu+ysT7L2c4abOQ44njFPsai//zl6se+z5Q/lquN8NlY3cfW2/v8ov63JXD+4gufP7916xbW2bZLFsrZYE/o9eDBg6yNpE0hO27fvh20OQffwp6Qv1B2/fr1A5/lamefjbDJbrnaaG3w9pZkY2ysotx+NnP2Iz53ITtL+TxmPREkqHwbZHZaIaa4HDdNkQGG/WDgvf+QJ+ibQZW8DaEAw54MUN9y8Cf3wQ43rmB1D/nM2purnV5v71OL3NfNyZe0w09OfgzmaqMPnLydJYxVOzb5ebRjsAQbNz7FFd9dNQEGP/D4QHCzHxh/okvZ8/ZDbU9kPsCw9sU+/CnbSdu8b7wvc7XTTkrWxpA9oTL4zvo/ZV9CN9oQC4xyttHa1ueHnG2MjTX4Ey+eY3K2sc93OjadgA8wilzkiWuB3Oz1Xntd1Na5qpzoG9jQfqC7V98Ws7s9AVw1S9XumO4he2N1c7AzNC79mgTUKcHGu3fvdub69QhkUIKNId/ZMVuCjfZzxbEJu2h7CTZan+n9cgSqCjCWw7atJC4IxAfZf9itJg8fPtxWsQV7w8mq/UbUvYZO2rnaCbuQhMDL2mgnYJbnaiOHBO3AeI35M2cbOaniL56+yZe3NVcbbdAAn8Iu2ugDxlxtXPD0JVERAkUGGBFbD4pr+lDYk54/cYxhtVUd6OlP0PZk5k9sXq9c7KTePGlzLCLwGNpysBFjDDYh4zbks5C9OdhIn/nzCPbHPO47BxvpG2SirJ3YL83G0DhU2ekEqg0wTkcnCWsTwEmYaXZMVj74WLv/teX7yakU+27evNmhmxNcrM18afkYl8xK2UuYudvuv4ggE4XsYkk2Lj0WJO+YgAKMYyZFl+QyieEbEidgnNz8CW/ISTnY6Scm2DtlYkrRRurPiQh+s77z+7n6kb7zttGeKeM1RT9av+DzB79Cz5i9ufpxSG8dP41AkQHGmA/s0ILJ07Cm1XrKyW5vzaGrTb/iW9PYSTcnOy3nKSftXGxEsIRMBl+0l+V94ywXG/ts8NkpXzd1G+051J9Pee7M3UbvE+0vT6DIACN28rbl/kOzPNrtJcaCJnsyS9lu6Mn0Ouj5xZAkmqudCJS4UC42wfCknauNU0Z9rjbCd/BlX+BL23K1se884QOLXG2cMlZVdz6BIh+0hfmJL97N234QrsrscyJ4PIe/bboyakPsOQshFinaavWEnf5Fn+Vqp9UbY5Gb9SnLc7UxNq7oW9jKLVcbY3pbP9LOWN0cPpMhHUP2hMrg41D72PhQeRkE2uC7M6T4J3naD7sd6D7oyM2t1i58sP0WsjV0cvft9t63J6mYDXZSjtVBuZ3E9rbL99+nN45Zn/bVTdlGbzP2Y2MwVxv79MYxu/XVTdmPQ5/JEmwMjVWVzSfgA4xiL5EgfRlL3V2c79qPfYFbe1IIWtWeyHpTusFGGxbGLhnEVMjVzr6xB5tsajpXG2M+C5XnamNMb9jofRyrm/pnEmMROoY2b5PfZ5vUbQzZprLlCFzD56FNZzQf/shHm3v37h19OJbraj9JvFaKD0zftcX9NFy+Z0zYnLT7rhUv3/O2EnO10+o9NC5ztXHKSMjVxlr8OPYcmqsfp4xV1Y0TwPocfLH/wosvNK981aubKgKMOA4dEQEREAEREAERWIKADzCKvUSyBCzJEAEREAEREAERmEdAAcY8bmolAiIgAiIgAiLQQ0ABRg8cHRIBERABERABEZhHQAHGPG5qJQIiIAIiIAIi0ENAAUYPHB0SAREQAREQARGYR0ABxjxuaiUCIiACIiACItBDQAFGDxwdEgEREAEREAERmEdAAcY8bmolAiIgAiIgAiLQQ0ABRg8cHRIBERABERABEZhHQAHGPG5qJQIiIAIiIAIi0ENAAUYPHB0SAREQAREQARGYR0ABxjxuaiUCIiACIiACItBDQAFGDxwdEgEREAEREAERmEdAAcY8bmolAiIgAiIgAiLQQ0ABRg8cHRIBERABERABEZhHQAHGPG5qJQIiIAIiIAIi0EPgyZ5jV4du3LgRrXZ2dhY9NvXAnTt3mrt37141u337doOyNTbqbW3zfeFYn+2o79tgn7LHsFlCD8jAy+pi92GD15NM++rZY2jvWYT6DfUTKlvDp5IpAiIgAiKQFoHzF1/4/PmtW7fOW7XOQ9v169fP8cJxvlgWqj+nrA0mDmTbfubI62vDvvCX24MHDw7sQ//2uJdHHmjHzerM976d3V9CD6t3yBbqib+hzevMeiEe1tZQv5Bvx8gYBiGdVCYCIiACIpAfgfZLZ6c0Ygqc/0dlMPBNFRv+3rx5s2knoe49N35Djf3lN13Uj2UFKK/V7UrutWvXmocPH17tx+TYb9d8H/rWTP2vBAbe2P4DhxvItVkWX4dZF/aF+njRPth/ih5oy0wC+IQ29kVbPEe2IaM2WLjyC3nzmJVBO2L9Ui7lhXRTmQiIgAiIQD0EBjMYjKP4rdV/G25RHWQ/7P7QN2HKtt+07Tdl3zdl4y/r2TK+78s+UCdbh9kEZmb62kMn1re6ooz7vg/PCDJ8HSs3pAcZ2T7RJiab7HDc+4z9hxiEZIb6iJVR91Cf1El/RUAEREAEyiLQfrntDGIGo50jmlUDDD8J+X2Ll8fsX06mvp3d7zsWcx/a2MnVBjiU1xdkhAKMkC0sQ/2QvCl6gEVMBuT4jQEMjvmgJFaXQQHa2ADBM0Z7X2b7I8+QXr5v7YuACIiACORPwAcYm91FgrT6UGodcxZe7UTVtBNpOzc13SUZu42Rw7a8VHAgILKD9D/6RRvogK3vUkhETNeedkIeN14qibWz9WJ64NIIL10MycFx1KcOnqNtD5mhS19j+rB12B/4gUMbZEwVofoiIAIiIAKFEBi1BmMJW+1kG5LHiRN/MVHhhQ2TvA0UhuSEZI8tY59j6/t60HOpidrLnrIPPfAiSwRc5OhtxD7WXaCOD17s+pcp/bMuZUMX3+8ceWojAiIgAiKQD4FFAwwbCBABvsViouIEg2/3KPN1sc8JjXWZQcD+GDmcUG27sa7gN3xmD9Bu6jdwyoC+nKzxni/I9HZ7/fr0gEy0p52+LfdRxzKI8bDM0ZY64y95ow51ZmYo1i/aoS/Ug82xfmPtVS4CIiACIlAWgZPXYGBdQIvk6IUrSu2EfVSOstAWksG6fXL62oX6QRnaQG9uIfmxtiinzX362T6oo5c5RQ+ua/D8QrI9E2srdbDrJGx9HvcyQrqjjt18G6+rl6F9ERABERCBMgi0X0Y7Q7jIE/c5Yqf58Ec+2ty7d+9q/UFbPnnDN118cw1t/BYcO27b9NUNHbNZkT4dQnr5spB8X2eL/SX02EvGEv1uwVh9iIAIiIAILEcAVyGQ/f7Ciy80r3zVq5tFA4zl1JwmKXbZZZoU1RYBERABERABEZhLwAcYT8wVlFI7RExjMiMp6SxdREAEREAERKBkAosu8twLFFPye/WvfkVABERABERABA4JFJHBkFNFQAREQAREQATSInCQwUAm4NRnH6RlnrQRAREQAREQARHYg8DVIs9/8spX7dG/+hQBERABERABESiIAO8iucpgoECbCIiACIiACIiACCxBQGswlqAoGSIgAiIgAiIgAgcEFGBoQIiACIiACIiACCxOoFuDsbhUCRQBERABERABEaiawH8HREuB4wWYcjAAAAAASUVORK5CYII=)
NOTE -- in debugging these plots, I found that the base QA code has an issue with respect to clone tracks. TPC only tracks are not considered when counting clones -- i.e. a track may be split into two. One is HFT, one is TPC only, and the base QA doesn't register this as a clone. The reason for this is that TPC only tracks are eliminated from study before the clone classification is made. Will need to revisit this.
Method: Reconstruct events with HFT in tracking. Create a TTree associating ID truth with number of reco hits in each layer. Process TPC only tracking sample and apply cuts based on ID truth / TTree lookup.
Figure 1 -- Shows corrspondence between the number of hits in the TPC (0-49), SST( 50-54), IST (55-59) and PXL (60-65).
z-axis is cut... < 100 entries are eliminated. Note that the z-axis is scaled by 1000. Looks like a good (enough) correlation to proceed.
NOTE -- in debugging these plots, I found that the base QA code has an issue with respect to clone tracks. TPC only tracks are not considered when counting clones -- i.e. a track may be split into two. One is HFT, one is TPC only, and the base QA doesn't register this as a clone. The reason for this is that TPC only tracks are eliminated from study before the clone classification is made. Will need to revisit this.
»
- jwebb's blog
- Login or register to post comments