Skip to content

Instantly share code, notes, and snippets.

@anaderi
Created June 12, 2013 15:13
Show Gist options
  • Save anaderi/5766173 to your computer and use it in GitHub Desktop.
Save anaderi/5766173 to your computer and use it in GitHub Desktop.
{
"metadata": {
"name": "Untitled0"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": "import pandas",
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 12
},
{
"cell_type": "code",
"collapsed": false,
"input": "import util",
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 13
},
{
"cell_type": "code",
"collapsed": false,
"input": "util.mogo_connection()",
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 14,
"text": "Connection('localhost', 27017)"
}
],
"prompt_number": 14
},
{
"cell_type": "code",
"collapsed": false,
"input": "util.set_of_tracks()",
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 15,
"text": "{None,\n 1,\n u'0test',\n u'2mesh_features-1',\n u'3_opts',\n u'3_reweight',\n u'None',\n u'bdt13',\n u'bdt13_pca',\n u'binary_facrors',\n u'pca1',\n u'pipi_options'}"
}
],
"prompt_number": 15
},
{
"cell_type": "code",
"collapsed": false,
"input": "l = list(util.Formula.find({\"analysis\":\"bsmumu\"}))",
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 18
},
{
"cell_type": "code",
"collapsed": false,
"input": "len(l)",
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 19,
"text": "467"
}
],
"prompt_number": 19
},
{
"cell_type": "code",
"collapsed": false,
"input": "l[0]",
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 20,
"text": "FORMULA (517db25abadd84374b5218c4)\nFOM : 36.143005092\n_id : 517db25abadd84374b5218c4\nanalysis : bsmumu\ncorr : {u'corr_bg_abs_mass': -0.0011597901546226638, u'corr_bg': -0.1912023049340575, u'corr_signal_abs_mass': 0.014105375934108001, u'corr_signal_mass': -0.012702035149387399, u'corr_bg_abs': -0.03730395491344927}\ncount_lines : None\ncreate_at : 2013-06-12 18:44:10.431060\ndir : /mnt/w76/notebook/f16_Weight_binB__B_n3_p1_medium.dir\nfeatures : B_s0_TAU_1000,sqrt_min_IPCHI2,isolation_sum_D1_D2,B_s0_PT,B_s0_CDFiso,min_muPT,B_s0_cosnk,B_s0_doca,B_s0_IP_OWNPV,B_s0_otherB_ang,abs_delta_phi,abs_delta_eta,B_s0_phi,B_s0_OWNPV_NDOF,B_s0_DIRA_OWNPV,B_s0_FD_OWNPV\nhost : w78.cern.yandex.net\noptions : -i 5000 -w 0.02 -x 64 -E -D 0\npredict : is_signal\ntest : /mnt/w76/data/bsmumu5_a.t\ntrack : 3_reweight\ntrain : /mnt/w76/notebook/bsmumu5_b.t.weight\nweight : Weight_binB__B_n3_p1"
}
],
"prompt_number": 20
},
{
"cell_type": "code",
"collapsed": false,
"input": "util.formula2table(l, [\"FOM\", \"dir\"], count=10)",
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": "Showing 10 of 467\n"
},
{
"html": "<table border=\"1\" cellpadding=\"3\" cellspacing=\"0\" style=\"border:1px solid black;border-collapse:collapse;\"><tr><td style=\"background-color:LightGray;\"><b>FOM</b></td><td style=\"background-color:LightGray;\"><b>dir</b></td></tr><tr><td style=\"background-color:Ivory;\">36.1430</td><td style=\"background-color:Ivory;\">/mnt/w76/notebook/f16_Weight_binB__B_n3_p1_medium.dir</td></tr><tr><td style=\"background-color:AliceBlue;\">40.4409</td><td style=\"background-color:AliceBlue;\">/mnt/w76/notebook/f16_Weight_binS__S_n3_p1_medium.dir</td></tr><tr><td style=\"background-color:Ivory;\">41.7615</td><td style=\"background-color:Ivory;\">/mnt/w76/notebook/f16_Weight_binT__T_n3_p1_medium.dir</td></tr><tr><td style=\"background-color:AliceBlue;\">36.4802</td><td style=\"background-color:AliceBlue;\">/mnt/w76/notebook/f16_Weight_binB__S_S_n3_p1_medium.dir</td></tr><tr><td style=\"background-color:Ivory;\">27.2398</td><td style=\"background-color:Ivory;\">/mnt/w76/notebook/f16_Weight_binB__T_T_n3_p1_medium.dir</td></tr><tr><td style=\"background-color:AliceBlue;\">41.1464</td><td style=\"background-color:AliceBlue;\">/mnt/w76/notebook/f16_Weight_binB__B_n10_p1_medium.dir</td></tr><tr><td style=\"background-color:Ivory;\">41.6093</td><td style=\"background-color:Ivory;\">/mnt/w76/notebook/f16_Weight_binS__S_n10_p1_medium.dir</td></tr><tr><td style=\"background-color:AliceBlue;\">41.7568</td><td style=\"background-color:AliceBlue;\">/mnt/w76/notebook/f16_Weight_binT__T_n10_p1_medium.dir</td></tr><tr><td style=\"background-color:Ivory;\">41.8000</td><td style=\"background-color:Ivory;\">/mnt/w76/notebook/f16_Weight_binB__S_S_n10_p1_medium.dir</td></tr><tr><td style=\"background-color:AliceBlue;\">41.0312</td><td style=\"background-color:AliceBlue;\">/mnt/w76/notebook/f16_Weight_binB__T_T_n10_p1_medium.dir</td></tr></table>",
"metadata": {},
"output_type": "pyout",
"prompt_number": 21,
"text": "<ipy_table.IpyTable at 0x1144ca1d0>"
}
],
"prompt_number": 21
},
{
"cell_type": "code",
"collapsed": false,
"input": "array = [f.FOM for f in l]",
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 22
},
{
"cell_type": "code",
"collapsed": false,
"input": "len(array)",
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 23,
"text": "467"
}
],
"prompt_number": 23
},
{
"cell_type": "code",
"collapsed": false,
"input": "hist(array, bins=100)",
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 25,
"text": "(array([ 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1,\n 1, 2, 1, 0, 0, 0, 0, 0, 2, 0, 0, 2, 0, 0, 2, 0, 1,\n 3, 2, 2, 4, 1, 3, 1, 3, 1, 0, 2, 3, 3, 5, 3, 2, 3,\n 1, 0, 4, 4, 1, 5, 3, 1, 2, 4, 5, 5, 4, 6, 2, 2, 3,\n 6, 2, 4, 4, 7, 4, 11, 4, 11, 12, 5, 15, 6, 14, 26, 12, 65,\n 18, 10, 13, 23, 9, 5, 0, 0, 0, 11, 60, 0, 0, 1, 1]),\n array([ 18.74621728, 19.01819801, 19.29017875, 19.56215948,\n 19.83414022, 20.10612095, 20.37810168, 20.65008242,\n 20.92206315, 21.19404388, 21.46602462, 21.73800535,\n 22.00998609, 22.28196682, 22.55394755, 22.82592829,\n 23.09790902, 23.36988976, 23.64187049, 23.91385122,\n 24.18583196, 24.45781269, 24.72979343, 25.00177416,\n 25.27375489, 25.54573563, 25.81771636, 26.0896971 ,\n 26.36167783, 26.63365856, 26.9056393 , 27.17762003,\n 27.44960077, 27.7215815 , 27.99356223, 28.26554297,\n 28.5375237 , 28.80950444, 29.08148517, 29.3534659 ,\n 29.62544664, 29.89742737, 30.1694081 , 30.44138884,\n 30.71336957, 30.98535031, 31.25733104, 31.52931177,\n 31.80129251, 32.07327324, 32.34525398, 32.61723471,\n 32.88921544, 33.16119618, 33.43317691, 33.70515765,\n 33.97713838, 34.24911911, 34.52109985, 34.79308058,\n 35.06506132, 35.33704205, 35.60902278, 35.88100352,\n 36.15298425, 36.42496499, 36.69694572, 36.96892645,\n 37.24090719, 37.51288792, 37.78486866, 38.05684939,\n 38.32883012, 38.60081086, 38.87279159, 39.14477232,\n 39.41675306, 39.68873379, 39.96071453, 40.23269526,\n 40.50467599, 40.77665673, 41.04863746, 41.3206182 ,\n 41.59259893, 41.86457966, 42.1365604 , 42.40854113,\n 42.68052187, 42.9525026 , 43.22448333, 43.49646407,\n 43.7684448 , 44.04042554, 44.31240627, 44.584387 ,\n 44.85636774, 45.12834847, 45.40032921, 45.67230994, 45.94429067]),\n <a list of 100 Patch objects>)"
},
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD9CAYAAAChtfywAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFipJREFUeJzt3X9MVff9x/HXbfH7h1OEO7kHJiY0a4lCEe5s65a13e3o\npdkaKFbDala86XRZmi2Z3TKmfyyj+6O9rG6p2v61LOZuS9qRLGW3jZLq7HVuJpI627q0K8uqAw33\nZu5eKmgNCuf7h+Mi5cL9zeUDz0dCcj33nnPefC68enj38znXYdu2LQCAsW4rdAEAgOwQ5ABgOIIc\nAAxHkAOA4QhyADAcQQ4Ahksa5B9++KHcbnf8a9WqVdq/f7+i0ai8Xq+qq6vV1NSk4eHh+agXAPAp\njnTmkU9MTGjNmjXq6+vTgQMHtHr1anV0dKirq0uxWEx+vz+ftQIAEkirtXL06FHdeeedWrt2rYLB\noHw+nyTJ5/Opp6cnLwUCAOaWVpC/+uqr2rZtmyQpEonIsixJkmVZikQiua8OAJBUyq2VsbExrVmz\nRu+//77KyspUWlqqWCwWf97pdCoajU4/uMOR22oBYIlI5+4pKV+RHz58WBs3blRZWZmkm1fh4XBY\nkjQ0NCSXyzVrMaZ+/fSnPy14DUuxduov/Bf1F/YrXSkH+SuvvBJvq0hSS0uLAoGAJCkQCKi1tTXt\nkwMAspdSkF+5ckVHjx7V448/Ht+2e/duHTlyRNXV1Tp27Jh2796dtyIBALMrSuVFn/nMZ3Tp0qVp\n25xOp44ePZqXohYKj8dT6BIyZnLtEvUXGvWbJa155Gkf3OHIqN8DAEtZutnJEn0AMBxBDgCGI8gB\nwHAEOQAYjiAHAMMR5ABgOIIcAAxHkAOA4QhyADAcQQ4AhiPIAeRUcbFTDodDDodDxcXOQpezJHCv\nFQA5dfMDZSZ/78mATHCvFQBYYghyADAcQQ4AhiPIAcBwBDkAGI4gBwDDEeQAYDiCHAAMR5ADgOEI\ncgAwXEpBPjw8rK1bt2r9+vWqqanRqVOnFI1G5fV6VV1draamJg0PD+e7VgBAAikF+fe//319/etf\n1wcffKD33ntP69atk9/vl9frVX9/vxobG+X3+/NdKwAggaQ3zfr444/ldrv10UcfTdu+bt06HT9+\nXJZlKRwOy+Px6B//+Mf0g3PTLGDJ4aZZ2Us3O4uSveDcuXMqKyvTU089pXfffVcbN27Uiy++qEgk\nIsuyJEmWZSkSiSTcv7OzM/7Y4/HI4/GkXBwApKK42KmRkZhWrizV5cvRQpeTtlAopFAolPH+Sa/I\n3377bX3pS1/SyZMnde+992rXrl1auXKlXnrpJcVisfjrnE6notHpA8gVObD0FOKKfOqciyNzcn4b\n28rKSlVWVuree++VJG3dulV/+9vfVF5ernA4LEkaGhqSy+XKsGQAQDaSBnl5ebnWrl2r/v5+SdLR\no0dVW1ur5uZmBQIBSVIgEFBra2t+KwUAJJTSJwS9++672rlzp8bGxvT5z39eBw8e1Pj4uNra2jQw\nMKCqqip1d3erpKRk+sFprQBLDq2V7KWbnXzUG4CcIsizx0e9AcASQ5ADgOEIcgAwHEEOAIYjyAHA\ncAQ5ABiOIAcAwxHkAGA4ghwADEeQA4DhCHIAMBxBDgCGI8gBwHAEOQAYjiAHAMMR5ABgOIIcAAxH\nkAOA4QhyADAcQQ4AhiPIAcBwBDkAGI4gBwDDFaXyoqqqKhUXF+v222/XsmXL1NfXp2g0qm984xv6\n97//raqqKnV3d6ukpCTf9QIAPiWlK3KHw6FQKKQzZ86or69PkuT3++X1etXf36/Gxkb5/f68FgoA\nSCzl1opt29P+HQwG5fP5JEk+n089PT25rQwAkJKUWisOh0MPP/ywbr/9dn3nO9/Rt7/9bUUiEVmW\nJUmyLEuRSCThvp2dnfHHHo9HHo8n66IBYDEJhUIKhUIZ7++wP32pncDQ0JAqKir0n//8R16vVwcO\nHFBLS4tisVj8NU6nU9FodPrBHY4ZV/IAFjeHwyFp8vd+fjJg6pyLI3PSzc6UWisVFRWSpLKyMm3e\nvFl9fX2yLEvhcFjSzaB3uVwZlAsAyFbSIL969apGRkYkSVeuXNGbb76puro6tbS0KBAISJICgYBa\nW1vzWykAIKGkrZVz585p8+bNkqQbN27om9/8pvbs2aNoNKq2tjYNDAzMOv2Q1gqw9NBayV662ZlS\nj3y+igFgPoI8e3npkQMAFi6CHAAMR5ADgOEIcgAwHEEOAIYjyAHAcAQ5ABiOIAcAwxHkAGA4ghwA\nDEeQA4DhCHIAMBxBDgCGI8gBwHAEOQAYjiAHAMMR5ABgOIIcAAxHkAOA4QhyADAcQQ4AhiPIAcBw\nBDkAGC6lIB8fH5fb7VZzc7MkKRqNyuv1qrq6Wk1NTRoeHs5rkQCA2aUU5Pv27VNNTY0cDockye/3\ny+v1qr+/X42NjfL7/XktEgAwu6RBfuHCBR06dEg7d+6UbduSpGAwKJ/PJ0ny+Xzq6enJb5UAgFkV\nJXvBM888oxdeeEGXL1+Ob4tEIrIsS5JkWZYikcis+3d2dsYfezweeTyezKsFgEUoFAopFAplvL/D\nnrzMTuCNN97Q4cOH9fLLLysUCukXv/iFXn/9dZWWlioWi8Vf53Q6FY1GZx7c4dAchwewCN1swU7+\n3s9PBkydc3FkTrrZOecV+cmTJxUMBnXo0CFdu3ZNly9fVnt7uyzLUjgcVnl5uYaGhuRyubIuHACQ\nmTmvyG91/Phx7d27V6+//ro6Ojr02c9+Vj/+8Y/l9/s1PDyc8H94ckUOLD1ckWcv3exMax755KyV\n3bt368iRI6qurtaxY8e0e/fu9KoEAORMylfkGR2cK3JgyeGKPHt5vSIHACw8BDkAGI4gBwDDEeQA\nYDiCHAAMR5ADgOEIcgAwHEEOAIYjyAHAcAQ5ABiOIAcAwxHkAGA4ghwADEeQA4DhCHIAMBxBDgCG\nI8gBwHAEOQAYjiAHAMMR5ABgOIIcAAxHkAOA4QhyADDcnEF+7do1bdq0SQ0NDaqpqdGePXskSdFo\nVF6vV9XV1WpqatLw8PC8FAsAmMlh27Y91wuuXr2q5cuX68aNG7r//vu1d+9eBYNBrV69Wh0dHerq\n6lIsFpPf7595cIdDSQ4PYJFxOBySJn/v5ycDps65ODIn3exM2lpZvny5JGlsbEzj4+MqLS1VMBiU\nz+eTJPl8PvX09GRYLgAgW0XJXjAxMaEvfOEL+te//qWnn35atbW1ikQisixLkmRZliKRyKz7d3Z2\nxh97PB55PJ6siwaAxSQUCikUCmW8f9LWyqSPP/5YjzzyiJ5//nk9/vjjisVi8eecTqei0ejMg9Na\nAZYcWivZy3lrZdKqVav06KOP6vTp07IsS+FwWJI0NDQkl8uVfqUAgJyYM8gvXboUn5HyySef6MiR\nI3K73WppaVEgEJAkBQIBtba25r9SAEBCc7ZWzp49K5/Pp4mJCU1MTKi9vV0/+tGPFI1G1dbWpoGB\nAVVVVam7u1slJSUzD05rBVhyaK1kL93sTLlHPh/FADAfQZ69vPXIAQALE0EOAIYjyAHAcAQ5ABiO\nIAcAwxHkAGA4ghwADEeQA8hYcbFTDodDDodDxcXOQpezZLEgCEDGEi3+YUFQ9lgQBABLDEEOAIYj\nyAHAcAQ5ABiOIAcAwxHkAGA4ghwADEeQAyiYyQVFLCbKDguCAGQs2wVBuVrIw4IgAIDRCHIAMBxB\nDgCGI8gBwHAEOQAYLmmQDw4O6qGHHlJtba3uvvtu7d+/X5IUjUbl9XpVXV2tpqYmDQ8P571YAMBM\nSacfhsNhhcNhNTQ0aHR0VBs3blRPT48OHjyo1atXq6OjQ11dXYrFYvL7/dMPzvRDYFFj+mF+5Hz6\nYXl5uRoaGiRJK1as0Pr163Xx4kUFg0H5fD5Jks/nU09PT4YlAwCyUZTOi8+fP68zZ85o06ZNikQi\nsixLkmRZliKRSMJ9Ojs74489Ho88Hk/GxQLAYhQKhRQKhTLeP+WVnaOjo/rKV76in/zkJ2ptbVVp\naalisVj8eafTqWg0Ov3gtFaARY3WSn7kZWXn9evXtWXLFrW3t6u1tVXSzavwcDgsSRoaGpLL5cqg\nXABAtpIGuW3b2rFjh2pqarRr16749paWFgUCAUlSIBCIBzwAYH4lba385S9/0YMPPqgNGzb8788X\n6fnnn9d9992ntrY2DQwMqKqqSt3d3SopKZl+cForwKJGayU/0s1O7n4IIGMEeX5w90MAC0gR9xuf\nB2lNPwSA9NyQZGtkxFHoQhY1rsgBwHAEOQAYjiAHAMMR5ABgOIIcAAxHkAOA4QhyADAcQQ4AhiPI\nAcyquNgph8PB6swFjnutAJhVsvumpHKvlbnugcK9VhLjXisAsMQQ5ADScmu7BQsDQQ4gLSMjMd1s\nY5jfwlgsCHIAMBxBDmBe0ZrJPYIcwLyiNZN7BDkAGI4gBwDDEeQAJOViFWcRfe8CYWUnAEmprdKc\nfeXmrY/nXtmZzsrP9GtfHJnDyk4AWGKSBvm3vvUtWZalurq6+LZoNCqv16vq6mo1NTVpeHg4r0UC\nyF6i1onJUwFNrj3Xkgb5U089pd7e3mnb/H6/vF6v+vv71djYKL/fn7cCAeTGrdP+bj42eyqgybXn\nWtIgf+CBB1RaWjptWzAYlM/nkyT5fD719PTkpzoAQFJFmewUiURkWZYkybIsRSKRWV/b2dkZf+zx\neOTxeDI5JQCjMaNlLqFQSKFQKOP9U5q1cv78eTU3N+vs2bOSpNLSUsVisfjzTqdT0Wh05sGZtQIs\nGOnOQMn1rJVk+6SbFfmY/bJQzMusFcuyFA6HJUlDQ0NyuVyZHAYAkAMZBXlLS4sCgYAkKRAIqLW1\nNadFAVhqivhIuSwkba1s27ZNx48f16VLl2RZln72s5/pscceU1tbmwYGBlRVVaXu7m6VlJTMPDit\nFWDBWOitlbk+Ui6V72cpt1ZY2QksEQS5OVjZCQBLDEEOAIYjyAHAcAQ5sAhlf0tamCSjlZ0AFrap\n+5BIIyOsqFzsuCIHAMMR5ABgOIIcKLB0+tmTry1M35sbXy1ULAgCCizRopvkr031dVPHzMWCoEwW\n97AgKH0sCAKAJYYgx6KTz/ZDJtP6bt3H4fg/Q6YFFrKNUmTA+CwstFaw6OTzE9XTaYPMtk8qLY/k\nx8p/ayVXbZJs9knne6S1AgAwFkGOeVPYGRepyWUbZPJYs8tN+2Kq5ql6k5+HGSiLCa0VzJt8tjxy\ndZ58tUFyNWMj0feW6XlM2IfWSmq4IgcAwxHkyJlsZ3RM7pPJApnEbYXEHx+WuxbP/MyuuHU8gERo\nrSBnkrUd5m4LTO2Tu/ZF4uPM1XrJdIFMohpz1VpJdk5aK6nvYwpaKwCwxBDkBsjlvaXnaiskO096\ndRTN0fLITmathqJ5aU9k3wZJddzm5/uBGWitGCCTRSjJj5X8T/hsn0/lT+ZMWivpnCfZn/O5bq0U\n6t4ki3UfWiup4YocAAyXVZD39vZq3bp1uuuuu9TV1ZWrmvLiy1/2yulcK6dzrb773R+mtE8oFMpv\nUXlkcu03hQpdQJZChS5gSTP/5z89GQf5+Pi4vve976m3t1fvv/++XnnlFX3wwQe5rC2n/v73vysW\ne02x2LM6e7Y/pX2y/WHIZGpdopWF0yWeUpfr2lOV2371rasNQ9O25WeKX2pjObO22Y8zJZSbEpGR\nT//8J5qmupBXGKcr4yDv6+vTnXfeqaqqKi1btkxPPPGE/vjHP+aytjxYI6ls3s429bmJ9v8eJ942\n2z7SdU31ACfdmHP/+Ta93lRNfQ+Jt8/clp/vNZ2xTFTbzONgYZr6OZ38nVoYvz+5kvGHL1+8eFFr\n166N/7uyslKnTp3KSVH5UFR0m1au3K7x8WEtW/a5QpcDADmT8ayVP/zhD+rt7dWvfvUrSdLvfvc7\nnTp1SgcOHJg6OFOjACAj6URzxlfka9as0eDgYPzfg4ODqqyszLgQAEBmMu6R33PPPfrnP/+p8+fP\na2xsTL///e/V0tKSy9oAACnI+Iq8qKhIL730kh555BGNj49rx44dWr9+fS5rAwCkIKt55F/72tf0\n4Ycf6sEHH9SLL76ourq6+HOdnZ2qrKyU2+2W2+1Wb29v1sXmy+DgoB566CHV1tbq7rvv1v79+yVJ\n0WhUXq9X1dXVampq0vDwcIErTWy2+k15D65du6ZNmzapoaFBNTU12rNnjyRzxn+2+k0Zf+nmdGK3\n263m5mZJ5oz9pE/Xb9LYV1VVacOGDXK73brvvvskpT/+OVmif+LECa1YsULbt2/X2bNnJUnPPvus\nVq5cqR/84AfZHj7vwuGwwuGwGhoaNDo6qo0bN6qnp0cHDx7U6tWr1dHRoa6uLsViMfn9/kKXO8Ns\n9Xd3dxvzHly9elXLly/XjRs3dP/992vv3r0KBoNGjL+UuP4//elPxoz/L3/5S50+fVojIyMKBoPq\n6OgwZuylmfWblD933HGHTp8+Ladzal57uuOfkyX6DzzwgEpLS2dsN+V/dpaXl6uhoUGStGLFCq1f\nv14XL15UMBiUz+eTJPl8PvX09BSyzFnNVr9kznuwfPlySdLY2JjGx8dVWlpqzPhLieuXzBj/Cxcu\n6NChQ9q5c2e8XpPGPlH9tm0bMfaTPl1ruuOf13utHDhwQPX19dqxY8eC/9Ns0vnz53XmzBlt2rRJ\nkUhElmVJkizLUiQSKXB1yU3W/8UvflGSOe/BxMSEGhoaZFlWvE1k0vgnql8yY/yfeeYZvfDCC7rt\ntqk4MGnsE9XvcDiMGHvpZq0PP/yw7rnnnvh07nTHP29B/vTTT+vcuXN65513VFFRoR/+MLX7mxTS\n6OiotmzZon379mnlypXTnjPhlqGjo6PaunWr9u3bpxUrVhj1Htx222165513dOHCBf35z3/WW2+9\nNe35hT7+n64/FAoZMf5vvPGGXC6X3G73rFewC3nsZ6vfhLGf9Ne//lVnzpzR4cOH9fLLL+vEiRPT\nnk9l/PMW5C6XK17Azp071dfXl69T5cT169e1ZcsWtbe3q7W1VdLN/xKGw2FJ0tDQkFwuVyFLnNNk\n/U8++WS8ftPeA0latWqVHn30UZ0+fdqo8Z80Wf/bb79txPifPHlSwWBQd9xxh7Zt26Zjx46pvb3d\nmLFPVP/27duNGPtJFRUVkqSysjJt3rxZfX19aY9/3oJ8aGgo/vi1116bNqNlobFtWzt27FBNTY12\n7doV397S0qJAICBJCgQC8YBcaGar35T34NKlS/E/fT/55BMdOXJEbrfbmPGfrf7JX0Rp4Y7/c889\np8HBQZ07d06vvvqqvvrVr+q3v/2tMWOfqP7f/OY3xvzsX716VSMjI5KkK1eu6M0331RdXV3642/n\nwBNPPGFXVFTYy5YtsysrK+1f//rXdnt7u11XV2dv2LDBfuyxx+xwOJyLU+XFiRMnbIfDYdfX19sN\nDQ12Q0ODffjwYfu///2v3djYaN9111221+u1Y7FYoUtNKFH9hw4dMuY9eO+992y3223X19fbdXV1\n9s9//nPbtm1jxn+2+k0Z/0mhUMhubm62bducsb/VW2+9Fa//ySefNGLsP/roI7u+vt6ur6+3a2tr\n7eeee8627fTHP6+fEAQAyD8+IQgADEeQA4DhCHIAMBxBDgCGI8gBwHAEOQAY7v8BvXl1Fvr1G0YA\nAAAASUVORK5CYII=\n",
"text": "<matplotlib.figure.Figure at 0x1142fdf50>"
}
],
"prompt_number": 25
},
{
"cell_type": "code",
"collapsed": false,
"input": "f = util.Formula.find_one({\"analysis\":\"bsmumu\", \"track\":\"bdt13_pca\"})",
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 26
},
{
"cell_type": "code",
"collapsed": false,
"input": "f",
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 27,
"text": "FORMULA (518badbdbadd842c84f191e2)\nFOM : 39.5212346138\n_id : 518badbdbadd842c84f191e2\nanalysis : bsmumu\ncreate_at : 2013-06-12 18:47:26.833806\ndir : /mnt/w76/notebook/bsmumu/bdt13_pca/f26_i1500.000000_w0.300000_x32.000000\nfeatures : [u'B_s0_TAU_1000', u'sqrt_min_IPCHI2', u'isolation_sum_D1_D2', u'B_s0_PT', u'B_s0_CDFiso', u'min_muPT', u'B_s0_cosnk', u'B_s0_doca', u'B_s0_IP_OWNPV', u'abs_delta_eta', u'abs_delta_phi', u'B_s0_otherB_ang', u'B_s0_otherB_boo_ang', u'B_s0_TAU_1000_pca', u'sqrt_min_IPCHI2_pca', u'isolation_sum_D1_D2_pca', u'B_s0_PT_pca', u'B_s0_CDFiso_pca', u'min_muPT_pca', u'B_s0_cosnk_pca', u'B_s0_doca_pca', u'B_s0_IP_OWNPV_pca', u'abs_delta_eta_pca', u'abs_delta_phi_pca', u'B_s0_otherB_ang_pca', u'B_s0_otherB_boo_ang_pca']\nfeatures_csv : B_s0_TAU_1000,sqrt_min_IPCHI2,isolation_sum_D1_D2,B_s0_PT,B_s0_CDFiso,min_muPT,B_s0_cosnk,B_s0_doca,B_s0_IP_OWNPV,abs_delta_eta,abs_delta_phi,B_s0_otherB_ang,B_s0_otherB_boo_ang,B_s0_TAU_1000_pca,sqrt_min_IPCHI2_pca,isolation_sum_D1_D2_pca,B_s0_PT_pca,B_s0_CDFiso_pca,min_muPT_pca,B_s0_cosnk_pca,B_s0_doca_pca,B_s0_IP_OWNPV_pca,abs_delta_eta_pca,abs_delta_phi_pca,B_s0_otherB_ang_pca,B_s0_otherB_boo_ang_pca\nmatrixnet_bin : /mnt/w76/notebook/bsmumu/bdt13_pca/f26_i1500.000000_w0.300000_x32.000000/data/matrixnet.bin\nmatrixnet_inc : /mnt/w76/notebook/bsmumu/bdt13_pca/f26_i1500.000000_w0.300000_x32.000000/data/matrixnet.inc\noptions : -i 1500.000000 -w 0.300000 -x 32.000000\ntest : /mnt/w76/data/bsmumu5_b.t.pca_bdt13\ntest_MN : /mnt/w76/notebook/bsmumu/bdt13_pca/f26_i1500.000000_w0.300000_x32.000000/data/test.MN\ntest_mass : None\ntrack : bdt13_pca\ntrain : /mnt/w76/data/bsmumu5_a.t.pca_bdt13\ntrain_MN : /mnt/w76/notebook/bsmumu/bdt13_pca/f26_i1500.000000_w0.300000_x32.000000/data/train.MN\nweight : None"
}
],
"prompt_number": 27
},
{
"cell_type": "code",
"collapsed": false,
"input": "f.plot_corr()",
"language": "python",
"metadata": {},
"outputs": [
{
"ename": "TypeError",
"evalue": "'NoneType' object is not iterable",
"output_type": "pyerr",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mTypeError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m<ipython-input-28-ba71dcd0665f>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_corr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
"\u001b[0;32m/Users/anaderi/notebook/util.py\u001b[0m in \u001b[0;36mplot_corr\u001b[0;34m(self, T_flat, bins)\u001b[0m\n\u001b[1;32m 476\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mplot_corr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mT_flat\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mNone\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbins\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m30\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 477\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mT_flat\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 478\u001b[0;31m \u001b[0mT_flat\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpandas_load_csvs\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mtest_MN_mass\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 479\u001b[0m \u001b[0;32massert\u001b[0m \u001b[0mMASS_COL\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mT_flat\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdtypes\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 480\u001b[0m \u001b[0mmn_flat\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0m_\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mflatten_classifier2\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mT_flat\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mMN\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mT_flat\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mis_signal\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/anaderi/notebook/util.py\u001b[0m in \u001b[0;36mpandas_load_csvs\u001b[0;34m(fnames)\u001b[0m\n\u001b[1;32m 1062\u001b[0m \u001b[0mfh_map\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m{\u001b[0m\u001b[0;34m}\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1063\u001b[0m \u001b[0mresult\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpandas\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mDataFrame\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1064\u001b[0;31m \u001b[0;32mfor\u001b[0m \u001b[0mf\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mfnames\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 1065\u001b[0m \u001b[0mfh\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1066\u001b[0m \u001b[0mnames\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfh\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mreadline\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mreplace\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"# \"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstrip\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msplit\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'\\t'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mTypeError\u001b[0m: 'NoneType' object is not iterable"
]
}
],
"prompt_number": 28
},
{
"cell_type": "code",
"collapsed": false,
"input": "",
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment