dbscan_grid_search_multi_weight.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To speed up training here instead of doing separate DBSCAN iterations for each set of provided weights we will alternate which weights we use during the existing iterations."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "_cell_guid": "e081740e-8169-4481-b1df-f5dd5488314f",
    "_uuid": "0bee86255243664f24e4bcf48af2228a3100a8b7"
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import os\n",
    "import math\n",
    "# from scipy.spatial import KDTree\n",
    "from sklearn.neighbors import KDTree\n",
    "from sklearn.neighbors import NearestNeighbors\n",
    "\n",
    "from trackml.dataset import load_event, load_dataset\n",
    "from trackml.score import score_event\n",
    "\n",
    "from multiprocessing import Pool\n",
    "\n",
    "from itertools import product\n",
    "\n",
    "def create_one_event_submission(event_id, hits, labels):\n",
    "    sub_data = np.column_stack(([event_id]*len(hits), hits.hit_id.values, labels))\n",
    "    submission = pd.DataFrame(data=sub_data, columns=[\"event_id\", \"hit_id\", \"track_id\"]).astype(int)\n",
    "    return submission\n",
    "\n",
    "# disable chain assignment warnings\n",
    "pd.options.mode.chained_assignment = None"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "def extend(submission, hits, limit=0.04, num_neighbours=18, threshold=0.01):\n",
    "\n",
    "    df = submission.merge(hits,  on=['hit_id'], how='left')\n",
    "    df = df.assign(d = np.sqrt( df.x**2 + df.y**2 + df.z**2 ))\n",
    "    df = df.assign(r = np.sqrt( df.x**2 + df.y**2))\n",
    "    df = df.assign(arctan2 = np.arctan2(df.z, df.r))\n",
    "\n",
    "    for angle in range(-90,90,1):\n",
    "\n",
    "        df1 = df.loc[(df.arctan2>(angle-1.5)/180*np.pi) & (df.arctan2<(angle+1.5)/180*np.pi)]\n",
    "\n",
    "        min_num_neighbours = len(df1)\n",
    "        if min_num_neighbours<3: continue\n",
    "\n",
    "        hit_ids = df1.hit_id.values\n",
    "        x,y,z = df1[['x', 'y', 'z']].values.T\n",
    "        r  = (x**2 + y**2)**0.5\n",
    "        r  = r/1000\n",
    "        a  = np.arctan2(y,x)\n",
    "        c = np.cos(a)\n",
    "        s = np.sin(a)\n",
    "\n",
    "        tree = KDTree(np.column_stack([c, s, r]), metric='euclidean')\n",
    "\n",
    "        track_ids = list(df1.track_id.unique())\n",
    "        num_track_ids = len(track_ids)\n",
    "        min_length=3\n",
    "\n",
    "        for i in range(num_track_ids):\n",
    "            p = track_ids[i]\n",
    "            if p==0: continue\n",
    "\n",
    "            idx = np.where(df1.track_id==p)[0]\n",
    "            if len(idx)<min_length: continue\n",
    "\n",
    "            if angle>0:\n",
    "                idx = idx[np.argsort( z[idx])]\n",
    "            else:\n",
    "                idx = idx[np.argsort(-z[idx])]\n",
    "\n",
    "\n",
    "            ## start and end points  ##\n",
    "            idx0,idx1 = idx[0],idx[-1]\n",
    "            a0 = a[idx0]\n",
    "            a1 = a[idx1]\n",
    "            r0 = r[idx0]\n",
    "            r1 = r[idx1]\n",
    "            c0 = c[idx0]\n",
    "            c1 = c[idx1]\n",
    "            s0 = s[idx0]\n",
    "            s1 = s[idx1]\n",
    "\n",
    "            da0 = a[idx[1]] - a[idx[0]]  #direction\n",
    "            dr0 = r[idx[1]] - r[idx[0]]\n",
    "            direction0 = np.arctan2(dr0,da0)\n",
    "\n",
    "            da1 = a[idx[-1]] - a[idx[-2]]\n",
    "            dr1 = r[idx[-1]] - r[idx[-2]]\n",
    "            direction1 = np.arctan2(dr1,da1)\n",
    "\n",
    "            ## extend start point\n",
    "            ns = tree.query([[c0, s0, r0]], k=min(num_neighbours, min_num_neighbours), return_distance=False)\n",
    "            ns = np.concatenate(ns)\n",
    "\n",
    "            direction = np.arctan2(r0 - r[ns], a0 - a[ns])\n",
    "            diff = 1 - np.cos(direction - direction0)\n",
    "            ns = ns[(r0 - r[ns] > threshold) & (diff < (1 - np.cos(limit)))]\n",
    "            for n in ns: df.loc[df.hit_id == hit_ids[n], 'track_id'] = p\n",
    "\n",
    "            ## extend end point\n",
    "            ns = tree.query([[c1, s1, r1]], k=min(num_neighbours, min_num_neighbours), return_distance=False)\n",
    "            ns = np.concatenate(ns)\n",
    "\n",
    "            direction = np.arctan2(r[ns] - r1, a[ns] - a1)\n",
    "            diff = 1 - np.cos(direction - direction1)\n",
    "            ns = ns[(r[ns] - r1 > threshold) & (diff < (1 - np.cos(limit)))]\n",
    "            for n in ns:  df.loc[df.hit_id == hit_ids[n], 'track_id'] = p\n",
    "\n",
    "    df = df[['event_id', 'hit_id', 'track_id']]\n",
    "    return df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "_cell_guid": "572fcbb6-8c7b-4a09-8916-8ec76689130f",
    "_uuid": "63414de98667e95f60407c9155899a25a321cffc"
   },
   "outputs": [],
   "source": [
    "# Change this according to your directory preferred setting\n",
    "path_to_train = \"./data/train_1\"\n",
    "\n",
    "# This event is in Train_1\n",
    "event_prefix = \"event000001000\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "weights_0 = np.array([1.1, 1.1, .375, .25, 0.05, 0.05, 0])\n",
    "weights_2 = np.array([1.05, 1.05, .37, .25, .01, .045, 0])\n",
    "weights_3 = np.array([1.08, 1.08, .38,.25, 0.009, 0.0001, 0])\n",
    "weights_4 = np.array([1.1, 1.1, .5, .25, 0.008, 0.001, 0])\n",
    "weights_5 = np.array([1.01, 1.01, .4048, .2, 2e-16, 2e-16, 0])\n",
    "weights_9 = np.array([1.02, 1.02, .38, .25, 0.009, 0.02, 0])\n",
    "\n",
    "weights_1 = np.array([1.08, 1.08, .38,.25, 0.0001, 0.0001, 0])\n",
    "weights_6 = np.array([1.08, 1.08, .38, .25, 2e-4, 0, 0])\n",
    "weights_8 = np.array([1.07, 1.07, .37, .25, 0.0002, 0.0001, 0])\n",
    "weights_12 = np.array([1.01, 1.01, .42, .25, 0.0023, 0.001, 0])\n",
    "weights_13 = np.array([1.0, 1.0, .40, .20, 0.0023, 0.0001, 0])\n",
    "weights_14 = np.array([1.02, 1.02, .40, .24, 0.0002, 0, 0])\n",
    "weights_7 = np.array([1.08, 1.08, .40, .20, 0.0023, 0, 0])\n",
    "\n",
    "weights_10 = np.array([1.02, 1.02, .42, .24, 0.0002, 0, 0])\n",
    "weights_11 = np.array([1.08, 1.08, .40, .20, 0.0023, 0.0001, 0])\n",
    "weights_15 = np.array([1.01, 1.01, .42, .25, 0.0002, 0, 0])\n",
    "weights_16 = np.array([1.005, 1.005, .40, .25, 0.0002, 0, 0])\n",
    "weights_17 = np.array([1.02, 1.02, .42, .24, 0.02, 0, 0])\n",
    "weights_18 = np.array([1.033, 1.033, .376, .2, 0.0002, 0, 0])\n",
    "weights_19 = np.array([0.9734, 0.9734, .4072, .2217, 0.0, 0.0, 0])\n",
    "weights_20 = np.array([1.0, 1.0, .42, .24, 0.0, 0, 0])\n",
    "weights_21 = np.array([1.02, 1.02, .42, .24, 0.0, 0, 0])\n",
    "weights_22 = np.array([1.05, 1.05, .4, .25, 0.0, 0, 0])\n",
    "weights_23 = np.array([1.05, 1.05, .5, .25, 0.0, 0, 0])\n",
    "weights_24 = np.array([1.05, 1.05, .4, .20, 0.0, 0, 0])\n",
    "weights_26 = np.array([1.2, 1.2, .45, .15, 0.0, 0, 0])\n",
    "weights_27 = np.array([1.2, 1.2, .45, .20, 0.0, 0, 0])\n",
    "weights_28 = np.array([1.2, 1.2, .45, .25, 0.0, 0, 0])\n",
    "weights_29 = np.array([1.09, 1.09, .5, .20, 0.0, 0, 0])\n",
    "weights_25 = np.array([1.09, 1.09, .45, .20, 0.0, 0, 0])\n",
    "weights_30 = np.array([0.9, .9, .35, .15, 0.0, 0.007, 0.007])\n",
    "weights_31 = np.array([1.09, 1.09, .45, .20, 0.0, 0.005, 0.005])\n",
    "weights_34 = np.array([1.09, 1.09, .45, .20, 0.0, 0.01, 0.01])\n",
    "weights_36 = np.array([1.0, 1.0, .35, .2, 0.0, 0.02, 0.02])\n",
    "weights_32 = np.array([0.9, .9, .35, .2, 0.0, 0.01, 0.01])\n",
    "weights_33 = np.array([1.0, 1.0, .35, .2, 0.0, 0.01, 0.01])\n",
    "weights_37 = np.array([0.9, .9, .4, .2, 0.0, 0.015, 0.015]) # raise w4\n",
    "weights_38 = np.array([0.9, .9, .4, .25, 0.0, 0.01, 0.01])  # raise w3\n",
    "weights_39 = np.array([0.9, .9, .45, .2, 0.0, 0.01, 0.01])\n",
    "weights_40 = np.array([0.95, .95, .4, .2, 0.0, 0.01, 0.01])\n",
    "weights_41 = np.array([0.95, .95, .45, .2, 0.0, 0.01, 0.01])  # 38 + 40\n",
    "weights_42 = np.array([0.95, .95, .4, .2, 0.0, 0.008, 0.008]) # w4 lowered\n",
    "weights_43 = np.array([1.0, 1.0, .4, .2, 0.0, 0.01, 0.01])    # w1 raised\n",
    "weights_44 = np.array([0.95, .95, .4, .2, 0.0, 0.012, 0.012]) # w4 raised\n",
    "weights_45 = np.array([1.0, 1.0, .4, .2, 0.0, 0.01, 0.01])  # raise w1 from 35\n",
    "weights_46 = np.array([0.9, .9, .4, .2, 0.0, 0.009, 0.009]) # lower w4\n",
    "weights_47 = np.array([0.85, .85, .4, .2, 0.0, 0.01, 0.01]) # lower w1\n",
    "weights_49 = np.array([0.9, .9, .4, .2, 0.0, 0.007, 0.007]) # still lower w4\n",
    "weights_48 = np.array([0.9, .9, .4, .2, 0.0, 0.008, 0.008]) # lower w4\n",
    "weights_51 = np.array([0.9, .9, .4, .2, 0.0, 0.009, 0.009]) # lower w4\n",
    "weights_50 = np.array([0.9, .9, .4, .2, 0.0, 0.011, 0.011]) # higher w4\n",
    "weights_52 = np.array([0.9, .9, .4, .2, 0.0, 0.012, 0.012]) # higher w4\n",
    "weights_53 = np.array([0.9, .9, .4, .2, 0.0, 0.013, 0.013]) # higher w4\n",
    "weights_54 = np.array([0.9, .9, .4, .2, 0.0, 0.014, 0.014]) # higher w4\n",
    "weights_55 = np.array([0.9, .9, .4, .2, 0.0, 0.015, 0.015]) # higher w4\n",
    "weights_56 = np.array([0.9, .9, .5, .2, 0.0, 0.01, 0.01])     # raise w2\n",
    "weights_60 = np.array([0.7, .7, .5, .13, 0.0, 0.01, 0.01])     # raise w2, lower w1, w3\n",
    "weights_59 = np.array([0.9, .9, .55, .1, 0.0, 0.01, 0.01])     # raise w2, lower w3\n",
    "weights_35 = np.array([0.9, .9, .4, .2, 0.0, 0.01, 0.01])     \n",
    "weights_57 = np.array([0.9, .9, .5, .15, 0.0, 0.01, 0.01]) # second best\n",
    "weights_61 = np.array([0.9, .9, .4, .1, 0.0, 0.01, 0.01])  # lower w3 some more\n",
    "weights_62 = np.array([0.9, .9, .45, .1, 0.0, 0.01, 0.01])  # lower w3 and raise w2\n",
    "weights_63 = np.array([0.9, .9, .5, .1, 0.0, 0.01, 0.01])  # raise w2 some more\n",
    "weights_64 = np.array([0.9, .9, .5, .07, 0.0, 0.01, 0.01])  # raise w2 and lower w3 some more\n",
    "weights_65 = np.array([0.9, .9, .35, .15, 0.0, 0.01, 0.01])   # lower w2 a bit\n",
    "weights_66 = np.array([0.9, .9, .35, .2, 0.0, 0.01, 0.01])   # lower w2 a bit and raise w3\n",
    "weights_67 = np.array([0.95, .95, .3, .2, 0.0, 0.012, 0.012])# weights from r optimization\n",
    "weights_68 = np.array([0.96, .96, .35, .2, 0.0, 0.0085, 0.0085])# weights from r optimization\n",
    "weights_69 = np.array([0.9, .9, .4, .15, 0.0, 0.0085, 0.0085])   # w4 down\n",
    "weights_70 = np.array([1.0, 1.0, .4, .1214, 0.0, 0.008, 0.008])   # w4 down\n",
    "weights_71 = np.array([0.9, 0.9, .5, .15, 0.0, 0.0085, 0.0085])   # w4 down\n",
    "weights_73 = np.array([0.9, .9, .45, .15, 0.0, 0.01, 0.01])   # w2 up a bit\n",
    "weights_74 = np.array([0.9, .9, .35, .2, 0.0, 0.011, 0.011])   # w2 down and w3 up\n",
    "weights_72 = np.array([0.9, .9, .4, .15, 0.0, 0.012, 0.012])   # w4 up a bit\n",
    "weights_75 = np.array([0.9, .9, .4, .15, 0.0, 0.014, 0.014]) \n",
    "weights_76 = np.array([0.9, .9, .4, .20, 0.0, 0.014, 0.014]) \n",
    "weights_77 = np.array([0.85, 0.85, .4, .20, 0.0, 0.012, 0.012]) \n",
    "weights_78 = np.array([0.9, .9, .4, .15, 0.005, 0.01, 0.01]) \n",
    "weights_79 = np.array([0.9, .9, .4, .15, 0.0005, 0.01, 0.01]) \n",
    "weights_80 = np.array([0.9, .9, .4, .15, 0.01, 0.01, 0.01]) \n",
    "weights_82 = np.array([0.96, .96, .45, .126, 0.0, 0.007, 0.007])  \n",
    "weights_84 = np.array([0.96, .96, .45, .15, 0.0, 0.0108, 0.0108])  \n",
    "weights_86 = np.array([0.9, .9, .4, .126, 0.0, 0.0108, 0.0108])  \n",
    "weights_85 = np.array([0.96, .96, .4, .15, 0.0, 0.0108, 0.0108])  \n",
    "weights_87 = np.array([1, 1, .4, .13, 0.0, 0.0108, 0.0108])  \n",
    "weights_88 = np.array([0.96, .96, .4, .13, 0.0, 0.0108, 0.0108])  \n",
    "weights_58 = np.array([0.9, .9, .4, .15, 0.0, 0.01, 0.01])   \n",
    "weights_81 = np.array([0.96, .96, .45, .126, 0.0, 0.0108, 0.0108])  \n",
    "weights_90 = np.array([0.96, .96, .4, .126, 0.15, 0.0108, 0.0108])  \n",
    "weights_91 = np.array([0.96, .96, .4, .126, 0.20, 0.0108, 0.0108]) \n",
    "weights_92 = np.array([0.96, .96, .2, .126, 0.30, 0.0108, 0.0108]) \n",
    "weights_93 = np.array([0.96, .96, .2, .126, 0.40, 0.0108, 0.0108])  \n",
    "weights_94 = np.array([0.96, .96, .3, .126, 0.25, 0.0108, 0.0108])  \n",
    "weights_95 = np.array([0.96, .96, .3, .126, 0.20, 0.0108, 0.0108]) \n",
    "weights_96 = np.array([0.96, .96, .3, .05, 0.21, 0.005, 0.005]) \n",
    "weights_97 = np.array([0.96, .96, .3, .05, 0.21, 0.008, 0.008]) \n",
    "weights_83 = np.array([0.96, .96, .4, .126, 0.0, 0.0108, 0.0108])  \n",
    "weights_100 = np.array([1.0, 1.0, .4, .05, 0.20, 0.0108, 0.0108]) \n",
    "weights_102 = np.array([0.96, .96, .4, .0005, 0.20, 0.0108, 0.0108]) \n",
    "weights_99 = np.array([0.96, .96, .35, .05, 0.16, 0.0108, 0.0108]) \n",
    "\n",
    "weights_101 = np.array([1.0, 1.0, .4, 0.005, 0.20, 0.0108, 0.0108]) \n",
    "weights_89 = np.array([0.96, .96, .4, .126, 0.1, 0.0108, 0.0108])  \n",
    "weights_98 = np.array([0.96, .96, .4, .05, 0.20, 0.0108, 0.0108]) \n",
    "weights_103 = np.array([1.2, 1.2, .1, .2, 0.30, 0.0108, 0.0108])  \n",
    "weights_104 = np.array([0.96, .96, .35, .2, 0.25, 0.008, 0.008]) \n",
    "weights_105 = np.array([0.96, .96, .12, .4, 0.25, 0.008, 0.008]) \n",
    "weights_106 = np.array([1.0, 1.0, .05, .3, 0.15, 0.012, 0.012]) \n",
    "weights_107 = np.array([0.95, 0.95, 0.05, .1, 0.35, 0.0083, 0.0083]) \n",
    "\n",
    "weights_arr = np.vstack([weights_0, weights_1, weights_2, weights_3, weights_4, weights_5, weights_6, weights_7, weights_8, weights_9, weights_10, weights_11, weights_12, weights_13, weights_14,  weights_15,  weights_16,  weights_17,  weights_18,  weights_19,  weights_20,  weights_21,  weights_22,  weights_23,  weights_24,  weights_25,  weights_26,  weights_27,  weights_28,  weights_29,  weights_30,  weights_31,  weights_32,  weights_33,  weights_34,  weights_35,  weights_36,  weights_37,  weights_38,  weights_39,  weights_40,  weights_41,  weights_42,  weights_43,  weights_44,  weights_45,  weights_46,  weights_47,  weights_48,  weights_49,  weights_50,  weights_51,  weights_52,  weights_53,  weights_54,  weights_55,  weights_56,  weights_57,  weights_58,  weights_59,  weights_60,  weights_61,  weights_62,  weights_63,  weights_64,  weights_65,  weights_66,  weights_67,  weights_68,  weights_69,  weights_70,  weights_71,  weights_72,  weights_73,  weights_74,  weights_75,  weights_76,  weights_77,  weights_78,  weights_79,  weights_80,  weights_81,  weights_82,  weights_83,  weights_84,  weights_85,  weights_86,  weights_87,  weights_88,  weights_89,  weights_90,  weights_91,  weights_92,  weights_93,  weights_94,  weights_95,  weights_96,  weights_97,  weights_98,  weights_99,  weights_100,  weights_101,  weights_102,  weights_103,  weights_104,  weights_105,  weights_106,  weights_107])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "_cell_guid": "e06d1ed7-5091-4d67-abb4-5984b137e2e6",
    "_uuid": "c2f70ae63abffcc09a534bb17fb89df8ffddb722",
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "from sklearn.preprocessing import StandardScaler\n",
    "import hdbscan\n",
    "from scipy import stats\n",
    "from tqdm import tqdm_notebook as tqdm\n",
    "from sklearn.cluster import DBSCAN\n",
    "\n",
    "class Clusterer(object):\n",
    "    def __init__(self,rz_scales=[0.65, 0.965, 1.5], eps=0.00285, dz0 = 0, stepdz = 1.5e-7, stepeps = 0.001205, num_loops=350, final_cluster=1, weight=[98, 83, 100, 101], weight_arr=weights_arr, max_size=16, step_z=0.0, size_incr=0.95, min_points=1, max_cluster_size=22, final_cluster_size=2, extension_loops=16, ext_threshold=0.01, threshold_delta=1.05, z_rounds=4, loop_mod=3, start_loop_mod=0):\n",
    "        self.max_cluster_size = max_cluster_size\n",
    "        self.rz_scales=rz_scales\n",
    "        self.epsilon = eps\n",
    "        self.dz0 = dz0\n",
    "        self.stepdz = stepdz\n",
    "        self.stepeps = stepeps / num_loops\n",
    "        self.num_loops = num_loops\n",
    "        self.final_cluster = final_cluster\n",
    "        self.weight_arr = weight_arr\n",
    "        self.weights = weight\n",
    "        self.max_size = max_size\n",
    "        self.step_z = step_z\n",
    "        if size_incr != 0:\n",
    "            self.size_incr = (self.max_cluster_size - max_size) / (num_loops * size_incr)\n",
    "        else:\n",
    "            self.size_incr = 0\n",
    "            \n",
    "        self.min_points = min_points\n",
    "        self.final_cluster_size = final_cluster_size\n",
    "        self.extension_loops = extension_loops\n",
    "        self.ext_threshold=ext_threshold\n",
    "        self.threshold_delta=threshold_delta\n",
    "        self.z_rounds = z_rounds\n",
    "        self.loop_mod = loop_mod\n",
    "        self.start_loop_mod = start_loop_mod\n",
    "        \n",
    "    # remove outliers\n",
    "    def _eliminate_outliers(self,labels,M):\n",
    "        norms=np.zeros((len(labels)),np.float32)\n",
    "        indices=np.zeros((len(labels)),np.float32)\n",
    "        for i, cluster in tqdm(enumerate(labels),total=len(labels)):\n",
    "            if cluster == 0:\n",
    "                continue\n",
    "            index = np.argwhere(self.clusters==cluster)\n",
    "            index = np.reshape(index,(index.shape[0]))\n",
    "            indices[i] = len(index)\n",
    "            x = M[index]\n",
    "            norms[i] = self._test_quadric(x)\n",
    "        threshold1 = np.percentile(norms,90)*5\n",
    "        threshold2 = 25\n",
    "        threshold3 = 5\n",
    "        for i, cluster in enumerate(labels):\n",
    "            if norms[i] > threshold1 or indices[i] > threshold2 or indices[i] < threshold3:\n",
    "                self.clusters[self.clusters==cluster]=0   \n",
    "    \n",
    "    # not sure what this function does?\n",
    "    def _test_quadric(self,x):\n",
    "        if x.size == 0 or len(x.shape)<2:\n",
    "            return 0\n",
    "        Z = np.zeros((x.shape[0],10), np.float32)\n",
    "        Z[:,0] = x[:,0]**2\n",
    "        Z[:,1] = 2*x[:,0]*x[:,1]\n",
    "        Z[:,2] = 2*x[:,0]*x[:,2]\n",
    "        Z[:,3] = 2*x[:,0]\n",
    "        Z[:,4] = x[:,1]**2\n",
    "        Z[:,5] = 2*x[:,1]*x[:,2]\n",
    "        Z[:,6] = 2*x[:,1]\n",
    "        Z[:,7] = x[:,2]**2\n",
    "        Z[:,8] = 2*x[:,2]\n",
    "        Z[:,9] = 1\n",
    "        v, s, t = np.linalg.svd(Z,full_matrices=False)        \n",
    "        smallest_index = np.argmin(np.array(s))\n",
    "        T = np.array(t)\n",
    "        T = T[smallest_index,:]        \n",
    "        norm = np.linalg.norm(np.dot(Z,T), ord=2)**2\n",
    "        return norm\n",
    "\n",
    "    # standard scale our data\n",
    "    def _preprocess(self, hits):\n",
    "        # old preprocess, may work better?\n",
    "        x = hits.x.values\n",
    "        y = hits.y.values\n",
    "        z = hits.z.values\n",
    "\n",
    "        r = np.sqrt(x**2 + y**2 + z**2)\n",
    "        hits['x2'] = x/r\n",
    "        hits['y2'] = y/r\n",
    "\n",
    "        r = np.sqrt(x**2 + y**2)\n",
    "        hits['z2'] = z/r\n",
    "\n",
    "        ss = StandardScaler()\n",
    "        X = ss.fit_transform(hits[['x2', 'y2', 'z2']].values)\n",
    "        for i, rz_scale in enumerate(self.rz_scales):\n",
    "            X[:,i] = X[:,i] * rz_scale\n",
    "       \n",
    "        return X\n",
    "    \n",
    "    def preprocess_outliers(self, hits):\n",
    "        hits2 = hits.copy()\n",
    "    \n",
    "        x,y,z = hits[['x', 'y', 'z']].values.T\n",
    "        r  = (x**2 + y**2)**0.5\n",
    "        r  = r/1000\n",
    "        a  = np.arctan2(y,x)\n",
    "        c = np.cos(a)\n",
    "        s = np.sin(a)\n",
    "        x2 = x / (r * 1000)\n",
    "        y2 = y / (r * 1000)\n",
    "        z2 = z / (r * 1000)\n",
    "        \n",
    "        hits2['c'] = c\n",
    "        hits2['s'] = s\n",
    "        hits2['r'] = r\n",
    "        hits2['x2'] = x2\n",
    "        hits2['y2'] = y2\n",
    "        hits2['z2'] = z2\n",
    "\n",
    "        return hits2\n",
    "    \n",
    "    def assign_outliers(self, hits):\n",
    "        # add the clusters to our hits\n",
    "        hits['track_id'] = self.clusters\n",
    "        \n",
    "        # get the size of each cluster\n",
    "        hits['N1'] = hits.groupby('track_id')['track_id'].transform('count')\n",
    "        \n",
    "        # if the cluster has size 1 or is 0\n",
    "        outliers_mask = ( (hits['N1'] == 1) | (hits['track_id'] == 0) )\n",
    "        outliers = hits[outliers_mask]\n",
    "        \n",
    "        # mask for even tracks\n",
    "        even_tracks_mask = ((hits['N1'] % 2) == 0)\n",
    "        even_tracks = hits[even_tracks_mask]\n",
    "        \n",
    "        # make sure we have outliers and even tracks to use\n",
    "        if (len(outliers) > 0) and (len(even_tracks) > 0):\n",
    "            print(\"Event:\", self.event_id, \"assigning outliers\")\n",
    "            \n",
    "            # preprocess the even tracks and fit nearest neighbors to it\n",
    "            X = self.preprocess_outliers(even_tracks)\n",
    "            outliers_proc = self.preprocess_outliers(outliers)\n",
    "\n",
    "            nbrs = NearestNeighbors(n_neighbors=1, algorithm='kd_tree').fit(X[['c', 's', 'r']])\n",
    "            indices = nbrs.kneighbors(outliers_proc[['c', 's', 'r']], n_neighbors=1, return_distance=False)\n",
    "            outliers['nbr'] = indices\n",
    "\n",
    "            nbrs = NearestNeighbors(n_neighbors=1, algorithm='kd_tree').fit(X)\n",
    "\n",
    "            # assign the outliers to the track of the nearest neighbor\n",
    "            for index, row in outliers.iterrows():\n",
    "                outliers.loc[index, 'track_id'] = X.iloc[int(row['nbr'])]['track_id']\n",
    "\n",
    "            # merge the outliers back into hits\n",
    "            hits.loc[outliers_mask, 'track_id'] = outliers['track_id']\n",
    "        \n",
    "        self.clusters = hits['track_id']\n",
    "        \n",
    "        return self.clusters\n",
    "\n",
    "    def _init(self,dfh):\n",
    "        dfh['s1'] = dfh.hit_id\n",
    "        dfh['N1'] =1\n",
    "        dfh['e1'] = self.epsilon\n",
    "        dfh['group_locked'] = 0\n",
    "        \n",
    "        mm = 1\n",
    "        dz0 = self.dz0\n",
    "        dfh['stepped_z'] = dfh.z\n",
    "        \n",
    "        # features which don't depend on z\n",
    "        dfh['rt'] = np.sqrt(dfh['x'].values**2+dfh['y'].values**2)\n",
    "        dfh['a0'] = np.arctan2(dfh['y'].values,dfh['x'].values)\n",
    "            \n",
    "        # try the origins in a different order?\n",
    "        if dz0 == 1:\n",
    "            z_shifts = [2, 3, 0, -2, -3]\n",
    "        elif dz0 == 2:\n",
    "            z_shifts = [3, 2, 0, -2, -3]\n",
    "        else:\n",
    "            z_shifts = [0, 2, 3, -2, -3]\n",
    "        \n",
    "        if type(self.weights) == int:\n",
    "            weights = [self.weights]\n",
    "        else:\n",
    "            weights = self.weights\n",
    "            \n",
    "        cx = self.weight_arr[self.weights]\n",
    "        start_loop = self.num_loops // 5\n",
    "        num_loops = start_loop + (self.num_loops // 2)\n",
    "        \n",
    "        # init the clusters\n",
    "        dfh['s1'] = 0\n",
    "        dfh['N1'] = 0\n",
    "        \n",
    "        for j, z_shift in enumerate(z_shifts[:self.z_rounds]):\n",
    "            dfh['stepped_z'] = dfh.z - z_shift\n",
    "\n",
    "            if j == 0:\n",
    "                num_loops = self.num_loops\n",
    "                start_loop = 0\n",
    "            else:\n",
    "                if self.start_loop_mod == 0:\n",
    "                    start_loop = 0\n",
    "                else:\n",
    "                    start_loop = self.num_loops // np.max([self.start_loop_mod - j, 2])\n",
    "\n",
    "                num_loops = self.num_loops - start_loop\n",
    "            \n",
    "            # features which depend on z\n",
    "            dfh['r'] = np.sqrt(dfh['x'].values**2+dfh['y'].values**2+dfh['stepped_z'].values**2)\n",
    "            dfh['z1'] = dfh['stepped_z'].values/dfh['rt'].values\n",
    "            dfh['z2'] = dfh['stepped_z']/dfh['r']\n",
    "            dfh['x2'] = dfh['x'].values/dfh['r'].values\n",
    "            dfh['y2'] = dfh['y'].values/dfh['r'].values\n",
    "            dfh['z3'] = np.arctan2(dfh['stepped_z'].values,dfh['rt'].values)\n",
    "#             dfh['t1'] = np.arcsin(dfh['rt']/(2*dfh['r']))\n",
    "#             dfh['t0'] = dfh['t1'] - dfh['a0']\n",
    "            \n",
    "            for ii in range(start_loop, num_loops):\n",
    "                # features which depend on ii\n",
    "                mm = mm*(-1)  \n",
    "                max_size = np.max([self.max_size + (ii * self.size_incr), self.max_cluster_size])\n",
    "                z_step = mm * self.step_z * ii\n",
    "                dz = mm*((ii*self.stepdz))\n",
    "\n",
    "                dfh['stepped_z'] = dfh['z'] + z_step    \n",
    "\n",
    "                dfh['a1'] = dfh['a0']+mm*(dfh['rt']+(dz)*dfh['rt']**2)/1000*(ii/self.loop_mod)/180*math.pi\n",
    "                dfh['sina1'] = np.sin(dfh['a1'].values)\n",
    "                dfh['cosa1'] = np.cos(dfh['a1'].values)\n",
    "\n",
    "                ss = StandardScaler()\n",
    "                dfs = ss.fit_transform(dfh[['sina1','cosa1','z1','z2','z3','x2','y2']].values)\n",
    "                weight = self.weights[ii % len(self.weights)]\n",
    "                cx = self.weight_arr[weight]\n",
    "                dfs = np.multiply(dfs, cx)\n",
    "\n",
    "                clusters=DBSCAN(eps=self.epsilon+(ii*self.stepeps),min_samples=self.min_points,metric='euclidean',n_jobs=1).fit(dfs).labels_ \n",
    "\n",
    "                # put our new clusters to another feature\n",
    "                dfh['s2'] = clusters\n",
    "                dfh['e2'] = self.epsilon+(ii*self.stepeps)\n",
    "\n",
    "                # get the count of those clusters\n",
    "                dfh['N2'] = dfh.groupby('s2')['s2'].transform('count')\n",
    "                maxs1 = dfh['s1'].max()\n",
    "\n",
    "                # if our new clusters are bigger, but less than our max size, use the new ones instead\n",
    "                cond = np.where((dfh['N2'].values > dfh['N1'].values) & (dfh['N2'].values < max_size))\n",
    "\n",
    "                s1 = dfh['s1'].values\n",
    "                s1[cond] = dfh['s2'].values[cond]+maxs1\n",
    "\n",
    "                # write the new clusters back to our dataframe\n",
    "                dfh['s1'] = s1\n",
    "                dfh['s1'] = dfh['s1'].astype('int64')\n",
    "                dfh['N1'] = dfh.groupby('s1')['s1'].transform('count')\n",
    "            \n",
    "        # for debugging\n",
    "        # return dfh\n",
    "        \n",
    "        # return our clusters\n",
    "        return dfh['s1'].values    \n",
    "    \n",
    "    def predict(self, hits, event_id):    \n",
    "        original_hits = hits.copy()\n",
    "        self.event_id = event_id\n",
    "        \n",
    "        # init our clusters\n",
    "        self.clusters = self._init(hits) \n",
    "        \n",
    "        # if we have unassigned points assign them with HDBSCAN\n",
    "        mask = self.clusters == 0\n",
    "        if (np.sum(mask) > 0) & (self.final_cluster >= 1):\n",
    "            # preprocess our data\n",
    "            X = self._preprocess(hits) \n",
    "\n",
    "            # create our last clusterer\n",
    "            cl = hdbscan.HDBSCAN(min_samples=1, min_cluster_size=self.final_cluster_size, metric='braycurtis', cluster_selection_method='leaf', algorithm='best', leaf_size=20)\n",
    "\n",
    "            # labels = unique clusters\n",
    "            labels = np.unique(self.clusters)\n",
    "\n",
    "            # remove outliers\n",
    "            self._eliminate_outliers(labels,X)\n",
    "            mask = self.clusters == 0\n",
    "            \n",
    "            n_labels = 0\n",
    "            # assign unassigned points with HDBSCANS\n",
    "            max_len = np.max(self.clusters)\n",
    "            \n",
    "            if np.sum(mask) > 3:\n",
    "                try:\n",
    "                    self.clusters[mask] = cl.fit_predict(X[mask])+max_len\n",
    "                except:\n",
    "                    print(\"Error with HDBSCAN\")\n",
    "                \n",
    "        # if we still have unassigned points assign them to their nearest neighbor\n",
    "        if self.final_cluster >= 2:\n",
    "            self.clusters = self.assign_outliers(hits)\n",
    "        \n",
    "        # create a submission\n",
    "        one_submission = create_one_event_submission(event_id, original_hits, self.clusters)\n",
    "        \n",
    "        threshold = self.ext_threshold\n",
    "        \n",
    "        # extend the submission\n",
    "        for i in range(self.extension_loops): \n",
    "            one_submission = extend(one_submission, original_hits, threshold=threshold)\n",
    "            \n",
    "            # decay the threshold\n",
    "            threshold = threshold * self.threshold_delta\n",
    "        \n",
    "        self.clusters = one_submission['track_id']\n",
    "        \n",
    "        return one_submission"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "results_df = pd.read_pickle(\"gs_results.pkl\")\n",
    "# results_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Length: 1\n"
     ]
    }
   ],
   "source": [
    "# defaults: \n",
    "# eps=0.0035, dz0 = -0.00070, stepdz = 0.00001, stepeps = 0.000005\n",
    "# 0.0035, 0.0005, 5e-05, 1e-09, 150, True, 10, True\n",
    "\n",
    "# create our params to iterate over\n",
    "eps_vals = [  0.00285,  ] \n",
    "stepdz_vals = [ 1.5e-7 ] \n",
    "stepeps_vals = [ 0.001155,  ] # \n",
    "loop_vals = [ 355, ] # 225\n",
    "weight_vals = [ [98, 89, 100, 101], ] # [98, 89, 101 ]\n",
    "step_z_vals = [ 0 ] # , 1.5e-5\n",
    "size_vals = [ 16 ]\n",
    "max_size_vals = [ 22 ]\n",
    "size_inc_vals = [ 0.95 ]\n",
    "final_cluster_size_vals = [ 1 ]\n",
    "final_cluster = [ 0  ] # unused\n",
    "ext_delta_vals = [ 1.05 ]\n",
    "z_round_vals = [ 5 ]\n",
    "loop_mod_vals = [ 3 ]\n",
    "start_loop_mod_vals = [ 0 ]\n",
    "\n",
    "dz0_vals = [ 1 ] # unused\n",
    "min_size_vals = [ 1 ] # unused\n",
    "min_size_change_vals = [ 0 ] # unused\n",
    "dz_incr = [ True ] # unused\n",
    "\n",
    "track_ext_round_vals = [ 8 ]\n",
    "\n",
    "step = \"eps\"\n",
    "\n",
    "foo = product(eps_vals, dz0_vals, stepdz_vals, stepeps_vals, loop_vals, final_cluster, weight_vals, dz_incr, size_vals, step_z_vals, size_inc_vals, max_size_vals, final_cluster_size_vals, track_ext_round_vals, ext_delta_vals, z_round_vals, loop_mod_vals, start_loop_mod_vals)\n",
    "\n",
    "offset = np.max(results_df.index.values) + 1\n",
    "\n",
    "iter_list = []\n",
    "for i, item in enumerate(foo):\n",
    "    iter_list.append([i + offset, item])\n",
    "    \n",
    "print(\"Length:\", len(iter_list))    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "def grid_search_loop(params):\n",
    "    counter = params[0]\n",
    "    eps, dz0, stepdz, stepeps, loops, final_cluster, weights, dz_incr, max_size, step_z, size_incr, max_cluster_size, final_cluster_size, track_ext_rounds, ext_delta, z_rounds, loop_mod, start_loop_mod = params[1]\n",
    "    step = \"eps\"\n",
    "    dataset_scores = []\n",
    "    print(\" - \", counter, \"Eps:\", eps, \"stepeps:\", stepeps, \"loops:\", loops, \"weight:\", weights, \"max size:\", max_size, \"step z:\", step_z, \"max:\", max_cluster_size, \"ext rounds:\", track_ext_rounds)\n",
    "    \n",
    "    for event_id, hits, cells, particles, truth in load_dataset(\"./data/train_1\", skip=50, nevents=2):\n",
    "        # Track pattern recognition\n",
    "        model = Clusterer(eps=eps, dz0=dz0, stepdz=stepdz, stepeps=stepeps, num_loops=loops, final_cluster=final_cluster, weight=weights, weight_arr=weights_arr, max_size=max_size, step_z=step_z, size_incr=size_incr, max_cluster_size=max_cluster_size, final_cluster_size=final_cluster_size, extension_loops=track_ext_rounds, threshold_delta=ext_delta, z_rounds= z_rounds, loop_mod=loop_mod, start_loop_mod=start_loop_mod)\n",
    "        one_submission = model.predict(hits.copy(), event_id)\n",
    "\n",
    "        # Score for the event\n",
    "        score = score_event(truth, one_submission)\n",
    "        dataset_scores.append(score)\n",
    "        print(event_id, params[1], \"score:\", score)\n",
    "        \n",
    "    mean_score = np.mean(dataset_scores)\n",
    "    \n",
    "    # create and return our results\n",
    "    result = [counter, [eps, dz0, stepdz, stepeps, loops, step, weights, final_cluster, max_size, step_z, size_incr, max_cluster_size, final_cluster_size, track_ext_rounds, ext_delta, z_rounds, loop_mod, start_loop_mod, mean_score]]\n",
    "    print(result)\n",
    "    \n",
    "#     return one_submission\n",
    "    return result\n",
    "\n",
    "def grid_search(iter_list, results_df=None, start=0, end=5):\n",
    "    pool = Pool(processes=8)\n",
    "    results = pool.map(grid_search_loop, iter_list[start:end])\n",
    "    pool.close()\n",
    "    \n",
    "    return results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " -  1801 Eps: 0.00285 stepeps: 0.001155 loops: 355 weight: [98, 89, 100, 101] max size: 16 step z: 0 max: 22 ext rounds: 8\n",
      "1050 (0.00285, 1, 1.5e-07, 0.001155, 355, 0, [98, 89, 100, 101], True, 16, 0, 0.95, 22, 1, 8, 1.05, 5, 3, 0) score: 0.5974128512048524\n"
     ]
    }
   ],
   "source": [
    "results = grid_search(iter_list, results_df, start=0, end=12)\n",
    "\n",
    "for item in results:\n",
    "    results_df.loc[item[0]] = item[1]\n",
    "    \n",
    "results_df.sort_values(\"acc\", ascending=False).to_pickle(\"gs_results.pkl\")\n",
    "results_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "results_df.sort_values(\"acc\", ascending=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# get one submission to play with\n",
    "submission = grid_search_loop(iter_list[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "results_df.sort_values(\"acc\", ascending=False).to_pickle(\"gs_results.pkl\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Scratch"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 95,
   "metadata": {},
   "outputs": [],
   "source": [
    "def preprocess_outliers(hits):\n",
    "    hits2 = hits.copy()\n",
    "    \n",
    "    x,y,z = hits[['x', 'y', 'z']].values.T\n",
    "    r  = (x**2 + y**2)**0.5\n",
    "    r  = r/1000\n",
    "    a  = np.arctan2(y,x)\n",
    "    c = np.cos(a)\n",
    "    s = np.sin(a)\n",
    "\n",
    "    hits2['c'] = c\n",
    "    hits2['s'] = s\n",
    "    hits2['r'] = r\n",
    "    \n",
    "    return hits2\n",
    "    \n",
    "def assign_outliers(hits):\n",
    "    # add the clusters to our hits\n",
    "#     hits['track_id'] = self.clusters\n",
    "\n",
    "    # get the size of each cluster\n",
    "    hits['N1'] = hits.groupby('track_id')['track_id'].transform('count')\n",
    "\n",
    "    # if the cluster has size 1 or is 0\n",
    "    mask = (hits['N1'] == 1 | hits['track_id'] == 0)\n",
    "    outliers = hits[mask]\n",
    "\n",
    "    # mask for even tracks\n",
    "    even_tracks_mask = ((hits['N1'] % 2) == 0)\n",
    "    even_tracks = hits[even_tracks_mask]\n",
    "\n",
    "    # preprocess the even tracks and fit nearest neighbors to it\n",
    "    X = preprocess_outliers(even_tracks)\n",
    "    nbrs = NearestNeighbors(n_neighbors=1, algorithm='kd_tree').fit(X)\n",
    "\n",
    "    # preprocess the outliers and get the nearest neighbor for each\n",
    "    outliers = preprocess_outliers(outliers)\n",
    "    indices = nbrs.kneighbors(outliers, n_neighbors=1, return_distance=False)\n",
    "    outliers['nbr'] = indices\n",
    "    \n",
    "    # assign the outliers to the track of the nearest neighbor\n",
    "    for row in outliers.iterrows():\n",
    "        neighbor = X.iloc[row['nbr']]\n",
    "        \n",
    "\n",
    "    return hits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 96,
   "metadata": {},
   "outputs": [],
   "source": [
    "for event_id, hits, cells, particles, truth in load_dataset(\"./data/train_1\", skip=50, nevents=1):\n",
    "    hits['track_id'] = submission['track_id']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 97,
   "metadata": {},
   "outputs": [],
   "source": [
    "hits['N1'] = hits.groupby('track_id')['track_id'].transform('count')\n",
    "\n",
    "mask = ( (hits['N1'] == 1) | (hits['track_id'] == 0) )\n",
    "outliers = hits[mask]\n",
    "\n",
    "# mask for even tracks\n",
    "even_tracks_mask = ((hits['N1'] % 2) == 0)\n",
    "even_tracks = hits[even_tracks_mask]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 98,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/skooch/.local/lib/python3.5/site-packages/ipykernel_launcher.py:6: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy\n",
      "  \n"
     ]
    }
   ],
   "source": [
    "X = preprocess_outliers(even_tracks)\n",
    "outliers_proc = preprocess_outliers(outliers)\n",
    "\n",
    "nbrs = NearestNeighbors(n_neighbors=1, algorithm='kd_tree').fit(X[['c', 's', 'r']])\n",
    "indices = nbrs.kneighbors(outliers_proc[['c', 's', 'r']], n_neighbors=1, return_distance=False)\n",
    "outliers['nbr'] = indices"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 99,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>hit_id</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>volume_id</th>\n",
       "      <th>layer_id</th>\n",
       "      <th>module_id</th>\n",
       "      <th>track_id</th>\n",
       "      <th>N1</th>\n",
       "      <th>nbr</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>32648</th>\n",
       "      <td>32649</td>\n",
       "      <td>62.141201</td>\n",
       "      <td>98.091499</td>\n",
       "      <td>-128.850006</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>295</td>\n",
       "      <td>1949622</td>\n",
       "      <td>1</td>\n",
       "      <td>17112</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33495</th>\n",
       "      <td>33496</td>\n",
       "      <td>-27.344801</td>\n",
       "      <td>114.413002</td>\n",
       "      <td>-51.299801</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>354</td>\n",
       "      <td>18409903</td>\n",
       "      <td>1</td>\n",
       "      <td>17360</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33551</th>\n",
       "      <td>33552</td>\n",
       "      <td>-67.969498</td>\n",
       "      <td>95.379700</td>\n",
       "      <td>-17.507700</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>357</td>\n",
       "      <td>11685994</td>\n",
       "      <td>1</td>\n",
       "      <td>16835</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34974</th>\n",
       "      <td>34975</td>\n",
       "      <td>-110.622002</td>\n",
       "      <td>36.811901</td>\n",
       "      <td>115.222000</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>466</td>\n",
       "      <td>9636902</td>\n",
       "      <td>1</td>\n",
       "      <td>2567</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41705</th>\n",
       "      <td>41706</td>\n",
       "      <td>-60.984100</td>\n",
       "      <td>-160.223999</td>\n",
       "      <td>369.039001</td>\n",
       "      <td>8</td>\n",
       "      <td>8</td>\n",
       "      <td>951</td>\n",
       "      <td>18521710</td>\n",
       "      <td>1</td>\n",
       "      <td>20502</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>42091</th>\n",
       "      <td>42092</td>\n",
       "      <td>171.774002</td>\n",
       "      <td>6.021760</td>\n",
       "      <td>468.938995</td>\n",
       "      <td>8</td>\n",
       "      <td>8</td>\n",
       "      <td>1054</td>\n",
       "      <td>18521332</td>\n",
       "      <td>1</td>\n",
       "      <td>18734</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58239</th>\n",
       "      <td>58240</td>\n",
       "      <td>-560.581970</td>\n",
       "      <td>-30.300100</td>\n",
       "      <td>-2948.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>2</td>\n",
       "      <td>3</td>\n",
       "      <td>18523761</td>\n",
       "      <td>1</td>\n",
       "      <td>30068</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59251</th>\n",
       "      <td>59252</td>\n",
       "      <td>-526.934998</td>\n",
       "      <td>441.579010</td>\n",
       "      <td>-2945.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>2</td>\n",
       "      <td>152</td>\n",
       "      <td>18521978</td>\n",
       "      <td>1</td>\n",
       "      <td>30562</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61271</th>\n",
       "      <td>61272</td>\n",
       "      <td>292.415009</td>\n",
       "      <td>534.953003</td>\n",
       "      <td>-2148.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>6</td>\n",
       "      <td>116</td>\n",
       "      <td>3937267</td>\n",
       "      <td>1</td>\n",
       "      <td>31021</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61513</th>\n",
       "      <td>61514</td>\n",
       "      <td>-324.071991</td>\n",
       "      <td>266.165009</td>\n",
       "      <td>-2151.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>6</td>\n",
       "      <td>151</td>\n",
       "      <td>18523119</td>\n",
       "      <td>1</td>\n",
       "      <td>46446</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62766</th>\n",
       "      <td>62767</td>\n",
       "      <td>-596.731018</td>\n",
       "      <td>345.621002</td>\n",
       "      <td>-1795.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>8</td>\n",
       "      <td>158</td>\n",
       "      <td>6087058</td>\n",
       "      <td>1</td>\n",
       "      <td>30574</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72418</th>\n",
       "      <td>72419</td>\n",
       "      <td>-361.653992</td>\n",
       "      <td>2.315940</td>\n",
       "      <td>-886.200012</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>112</td>\n",
       "      <td>18521053</td>\n",
       "      <td>1</td>\n",
       "      <td>35623</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72881</th>\n",
       "      <td>72882</td>\n",
       "      <td>-356.710999</td>\n",
       "      <td>-43.158298</td>\n",
       "      <td>-606.119995</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>226</td>\n",
       "      <td>18520989</td>\n",
       "      <td>1</td>\n",
       "      <td>35966</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72907</th>\n",
       "      <td>72908</td>\n",
       "      <td>-271.240997</td>\n",
       "      <td>-238.298004</td>\n",
       "      <td>-578.099976</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>231</td>\n",
       "      <td>7302158</td>\n",
       "      <td>1</td>\n",
       "      <td>30708</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74722</th>\n",
       "      <td>74723</td>\n",
       "      <td>-29.566700</td>\n",
       "      <td>-355.957001</td>\n",
       "      <td>-60.400002</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>518</td>\n",
       "      <td>7220777</td>\n",
       "      <td>1</td>\n",
       "      <td>35823</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76207</th>\n",
       "      <td>76208</td>\n",
       "      <td>54.579601</td>\n",
       "      <td>358.360992</td>\n",
       "      <td>216.199997</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>713</td>\n",
       "      <td>5601208</td>\n",
       "      <td>1</td>\n",
       "      <td>36490</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76360</th>\n",
       "      <td>76361</td>\n",
       "      <td>-298.911987</td>\n",
       "      <td>-194.498001</td>\n",
       "      <td>291.600006</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>734</td>\n",
       "      <td>6568774</td>\n",
       "      <td>1</td>\n",
       "      <td>37581</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79306</th>\n",
       "      <td>79307</td>\n",
       "      <td>-495.473999</td>\n",
       "      <td>57.894001</td>\n",
       "      <td>-746.000000</td>\n",
       "      <td>13</td>\n",
       "      <td>6</td>\n",
       "      <td>311</td>\n",
       "      <td>5200397</td>\n",
       "      <td>1</td>\n",
       "      <td>38948</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79962</th>\n",
       "      <td>79963</td>\n",
       "      <td>430.200989</td>\n",
       "      <td>-255.901001</td>\n",
       "      <td>-368.200012</td>\n",
       "      <td>13</td>\n",
       "      <td>6</td>\n",
       "      <td>501</td>\n",
       "      <td>18525024</td>\n",
       "      <td>1</td>\n",
       "      <td>39158</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82254</th>\n",
       "      <td>82255</td>\n",
       "      <td>-363.391998</td>\n",
       "      <td>-339.842987</td>\n",
       "      <td>327.600006</td>\n",
       "      <td>13</td>\n",
       "      <td>6</td>\n",
       "      <td>1024</td>\n",
       "      <td>8096766</td>\n",
       "      <td>1</td>\n",
       "      <td>39734</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>84339</th>\n",
       "      <td>84340</td>\n",
       "      <td>155.970993</td>\n",
       "      <td>-642.041016</td>\n",
       "      <td>-1078.199951</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>30</td>\n",
       "      <td>2100380</td>\n",
       "      <td>1</td>\n",
       "      <td>42108</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>85129</th>\n",
       "      <td>85130</td>\n",
       "      <td>-404.596008</td>\n",
       "      <td>-526.367004</td>\n",
       "      <td>-608.143005</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>423</td>\n",
       "      <td>18525192</td>\n",
       "      <td>1</td>\n",
       "      <td>41325</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86527</th>\n",
       "      <td>86528</td>\n",
       "      <td>281.730988</td>\n",
       "      <td>-596.346008</td>\n",
       "      <td>-131.199997</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>951</td>\n",
       "      <td>18526337</td>\n",
       "      <td>1</td>\n",
       "      <td>41215</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89115</th>\n",
       "      <td>89116</td>\n",
       "      <td>-121.209999</td>\n",
       "      <td>650.249023</td>\n",
       "      <td>867.062012</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>1916</td>\n",
       "      <td>18524296</td>\n",
       "      <td>1</td>\n",
       "      <td>42495</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89755</th>\n",
       "      <td>89756</td>\n",
       "      <td>-230.985992</td>\n",
       "      <td>-497.566010</td>\n",
       "      <td>1218.500000</td>\n",
       "      <td>14</td>\n",
       "      <td>2</td>\n",
       "      <td>31</td>\n",
       "      <td>18522458</td>\n",
       "      <td>1</td>\n",
       "      <td>31318</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100122</th>\n",
       "      <td>100123</td>\n",
       "      <td>-592.258972</td>\n",
       "      <td>-565.510986</td>\n",
       "      <td>-737.200012</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>375</td>\n",
       "      <td>18525052</td>\n",
       "      <td>1</td>\n",
       "      <td>49452</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100497</th>\n",
       "      <td>100498</td>\n",
       "      <td>-33.920502</td>\n",
       "      <td>822.763977</td>\n",
       "      <td>-612.599976</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>571</td>\n",
       "      <td>18525130</td>\n",
       "      <td>1</td>\n",
       "      <td>49413</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>101337</th>\n",
       "      <td>101338</td>\n",
       "      <td>-732.999023</td>\n",
       "      <td>360.893005</td>\n",
       "      <td>-260.399994</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>952</td>\n",
       "      <td>18526139</td>\n",
       "      <td>1</td>\n",
       "      <td>49306</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>103330</th>\n",
       "      <td>103331</td>\n",
       "      <td>-600.182983</td>\n",
       "      <td>558.385010</td>\n",
       "      <td>542.000000</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>1906</td>\n",
       "      <td>18525814</td>\n",
       "      <td>1</td>\n",
       "      <td>49968</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104496</th>\n",
       "      <td>104497</td>\n",
       "      <td>126.666000</td>\n",
       "      <td>-1010.719971</td>\n",
       "      <td>-1046.199951</td>\n",
       "      <td>17</td>\n",
       "      <td>4</td>\n",
       "      <td>42</td>\n",
       "      <td>18523429</td>\n",
       "      <td>1</td>\n",
       "      <td>51786</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104709</th>\n",
       "      <td>104710</td>\n",
       "      <td>711.646973</td>\n",
       "      <td>-730.331970</td>\n",
       "      <td>-954.000000</td>\n",
       "      <td>17</td>\n",
       "      <td>4</td>\n",
       "      <td>209</td>\n",
       "      <td>18524734</td>\n",
       "      <td>1</td>\n",
       "      <td>50721</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>106893</th>\n",
       "      <td>106894</td>\n",
       "      <td>1003.309998</td>\n",
       "      <td>204.781998</td>\n",
       "      <td>-48.599998</td>\n",
       "      <td>17</td>\n",
       "      <td>4</td>\n",
       "      <td>1601</td>\n",
       "      <td>18525994</td>\n",
       "      <td>1</td>\n",
       "      <td>52238</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>109465</th>\n",
       "      <td>109466</td>\n",
       "      <td>-572.379028</td>\n",
       "      <td>-761.374023</td>\n",
       "      <td>1222.500000</td>\n",
       "      <td>18</td>\n",
       "      <td>2</td>\n",
       "      <td>15</td>\n",
       "      <td>3714269</td>\n",
       "      <td>1</td>\n",
       "      <td>52537</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111716</th>\n",
       "      <td>111717</td>\n",
       "      <td>-522.752014</td>\n",
       "      <td>-656.690002</td>\n",
       "      <td>2947.500000</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>14</td>\n",
       "      <td>18522615</td>\n",
       "      <td>1</td>\n",
       "      <td>47105</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111919</th>\n",
       "      <td>111920</td>\n",
       "      <td>616.473022</td>\n",
       "      <td>455.427002</td>\n",
       "      <td>2947.500000</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>58</td>\n",
       "      <td>18523047</td>\n",
       "      <td>1</td>\n",
       "      <td>53218</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        hit_id            x            y            z  volume_id  layer_id  \\\n",
       "32648    32649    62.141201    98.091499  -128.850006          8         6   \n",
       "33495    33496   -27.344801   114.413002   -51.299801          8         6   \n",
       "33551    33552   -67.969498    95.379700   -17.507700          8         6   \n",
       "34974    34975  -110.622002    36.811901   115.222000          8         6   \n",
       "41705    41706   -60.984100  -160.223999   369.039001          8         8   \n",
       "42091    42092   171.774002     6.021760   468.938995          8         8   \n",
       "58239    58240  -560.581970   -30.300100 -2948.500000         12         2   \n",
       "59251    59252  -526.934998   441.579010 -2945.500000         12         2   \n",
       "61271    61272   292.415009   534.953003 -2148.500000         12         6   \n",
       "61513    61514  -324.071991   266.165009 -2151.500000         12         6   \n",
       "62766    62767  -596.731018   345.621002 -1795.500000         12         8   \n",
       "72418    72419  -361.653992     2.315940  -886.200012         13         4   \n",
       "72881    72882  -356.710999   -43.158298  -606.119995         13         4   \n",
       "72907    72908  -271.240997  -238.298004  -578.099976         13         4   \n",
       "74722    74723   -29.566700  -355.957001   -60.400002         13         4   \n",
       "76207    76208    54.579601   358.360992   216.199997         13         4   \n",
       "76360    76361  -298.911987  -194.498001   291.600006         13         4   \n",
       "79306    79307  -495.473999    57.894001  -746.000000         13         6   \n",
       "79962    79963   430.200989  -255.901001  -368.200012         13         6   \n",
       "82254    82255  -363.391998  -339.842987   327.600006         13         6   \n",
       "84339    84340   155.970993  -642.041016 -1078.199951         13         8   \n",
       "85129    85130  -404.596008  -526.367004  -608.143005         13         8   \n",
       "86527    86528   281.730988  -596.346008  -131.199997         13         8   \n",
       "89115    89116  -121.209999   650.249023   867.062012         13         8   \n",
       "89755    89756  -230.985992  -497.566010  1218.500000         14         2   \n",
       "100122  100123  -592.258972  -565.510986  -737.200012         17         2   \n",
       "100497  100498   -33.920502   822.763977  -612.599976         17         2   \n",
       "101337  101338  -732.999023   360.893005  -260.399994         17         2   \n",
       "103330  103331  -600.182983   558.385010   542.000000         17         2   \n",
       "104496  104497   126.666000 -1010.719971 -1046.199951         17         4   \n",
       "104709  104710   711.646973  -730.331970  -954.000000         17         4   \n",
       "106893  106894  1003.309998   204.781998   -48.599998         17         4   \n",
       "109465  109466  -572.379028  -761.374023  1222.500000         18         2   \n",
       "111716  111717  -522.752014  -656.690002  2947.500000         18        12   \n",
       "111919  111920   616.473022   455.427002  2947.500000         18        12   \n",
       "\n",
       "        module_id  track_id  N1    nbr  \n",
       "32648         295   1949622   1  17112  \n",
       "33495         354  18409903   1  17360  \n",
       "33551         357  11685994   1  16835  \n",
       "34974         466   9636902   1   2567  \n",
       "41705         951  18521710   1  20502  \n",
       "42091        1054  18521332   1  18734  \n",
       "58239           3  18523761   1  30068  \n",
       "59251         152  18521978   1  30562  \n",
       "61271         116   3937267   1  31021  \n",
       "61513         151  18523119   1  46446  \n",
       "62766         158   6087058   1  30574  \n",
       "72418         112  18521053   1  35623  \n",
       "72881         226  18520989   1  35966  \n",
       "72907         231   7302158   1  30708  \n",
       "74722         518   7220777   1  35823  \n",
       "76207         713   5601208   1  36490  \n",
       "76360         734   6568774   1  37581  \n",
       "79306         311   5200397   1  38948  \n",
       "79962         501  18525024   1  39158  \n",
       "82254        1024   8096766   1  39734  \n",
       "84339          30   2100380   1  42108  \n",
       "85129         423  18525192   1  41325  \n",
       "86527         951  18526337   1  41215  \n",
       "89115        1916  18524296   1  42495  \n",
       "89755          31  18522458   1  31318  \n",
       "100122        375  18525052   1  49452  \n",
       "100497        571  18525130   1  49413  \n",
       "101337        952  18526139   1  49306  \n",
       "103330       1906  18525814   1  49968  \n",
       "104496         42  18523429   1  51786  \n",
       "104709        209  18524734   1  50721  \n",
       "106893       1601  18525994   1  52238  \n",
       "109465         15   3714269   1  52537  \n",
       "111716         14  18522615   1  47105  \n",
       "111919         58  18523047   1  53218  "
      ]
     },
     "execution_count": 99,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "outliers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 106,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "18524844.0\n",
      "8719220.0\n",
      "1575617.0\n",
      "11582980.0\n",
      "2789076.0\n",
      "18520929.0\n",
      "6648030.0\n",
      "18521324.0\n",
      "18522229.0\n",
      "1347411.0\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/skooch/.local/lib/python3.5/site-packages/pandas/core/indexing.py:543: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy\n",
      "  self.obj[item] = s\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "165095.0\n",
      "5919807.0\n",
      "18523733.0\n",
      "1864044.0\n",
      "18525541.0\n",
      "18526747.0\n",
      "2552722.0\n",
      "18521714.0\n",
      "18526105.0\n",
      "6737515.0\n",
      "18525952.0\n",
      "18525680.0\n",
      "3946074.0\n",
      "18526114.0\n",
      "18524018.0\n",
      "4025041.0\n",
      "3906428.0\n",
      "7227919.0\n",
      "18524913.0\n",
      "18525458.0\n",
      "18525977.0\n",
      "5489690.0\n",
      "18522829.0\n",
      "18524172.0\n",
      "18523786.0\n"
     ]
    }
   ],
   "source": [
    "for index, row in outliers.iterrows():\n",
    "    print(X.iloc[int(row['nbr'])]['track_id'])\n",
    "    outliers.loc[index, 'track_id'] = X.iloc[int(row['nbr'])]['track_id']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 108,
   "metadata": {},
   "outputs": [],
   "source": [
    "hits.loc[mask, 'track_id'] = outliers['track_id']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 109,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>hit_id</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>volume_id</th>\n",
       "      <th>layer_id</th>\n",
       "      <th>module_id</th>\n",
       "      <th>track_id</th>\n",
       "      <th>N1</th>\n",
       "      <th>nbr</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>32648</th>\n",
       "      <td>32649</td>\n",
       "      <td>62.141201</td>\n",
       "      <td>98.091499</td>\n",
       "      <td>-128.850006</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>295</td>\n",
       "      <td>18524844.0</td>\n",
       "      <td>1</td>\n",
       "      <td>17112</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33495</th>\n",
       "      <td>33496</td>\n",
       "      <td>-27.344801</td>\n",
       "      <td>114.413002</td>\n",
       "      <td>-51.299801</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>354</td>\n",
       "      <td>8719220.0</td>\n",
       "      <td>1</td>\n",
       "      <td>17360</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33551</th>\n",
       "      <td>33552</td>\n",
       "      <td>-67.969498</td>\n",
       "      <td>95.379700</td>\n",
       "      <td>-17.507700</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>357</td>\n",
       "      <td>1575617.0</td>\n",
       "      <td>1</td>\n",
       "      <td>16835</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34974</th>\n",
       "      <td>34975</td>\n",
       "      <td>-110.622002</td>\n",
       "      <td>36.811901</td>\n",
       "      <td>115.222000</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>466</td>\n",
       "      <td>11582980.0</td>\n",
       "      <td>1</td>\n",
       "      <td>2567</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41705</th>\n",
       "      <td>41706</td>\n",
       "      <td>-60.984100</td>\n",
       "      <td>-160.223999</td>\n",
       "      <td>369.039001</td>\n",
       "      <td>8</td>\n",
       "      <td>8</td>\n",
       "      <td>951</td>\n",
       "      <td>2789076.0</td>\n",
       "      <td>1</td>\n",
       "      <td>20502</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>42091</th>\n",
       "      <td>42092</td>\n",
       "      <td>171.774002</td>\n",
       "      <td>6.021760</td>\n",
       "      <td>468.938995</td>\n",
       "      <td>8</td>\n",
       "      <td>8</td>\n",
       "      <td>1054</td>\n",
       "      <td>18520929.0</td>\n",
       "      <td>1</td>\n",
       "      <td>18734</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58239</th>\n",
       "      <td>58240</td>\n",
       "      <td>-560.581970</td>\n",
       "      <td>-30.300100</td>\n",
       "      <td>-2948.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>2</td>\n",
       "      <td>3</td>\n",
       "      <td>6648030.0</td>\n",
       "      <td>1</td>\n",
       "      <td>30068</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59251</th>\n",
       "      <td>59252</td>\n",
       "      <td>-526.934998</td>\n",
       "      <td>441.579010</td>\n",
       "      <td>-2945.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>2</td>\n",
       "      <td>152</td>\n",
       "      <td>18521324.0</td>\n",
       "      <td>1</td>\n",
       "      <td>30562</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61271</th>\n",
       "      <td>61272</td>\n",
       "      <td>292.415009</td>\n",
       "      <td>534.953003</td>\n",
       "      <td>-2148.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>6</td>\n",
       "      <td>116</td>\n",
       "      <td>18522229.0</td>\n",
       "      <td>1</td>\n",
       "      <td>31021</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61513</th>\n",
       "      <td>61514</td>\n",
       "      <td>-324.071991</td>\n",
       "      <td>266.165009</td>\n",
       "      <td>-2151.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>6</td>\n",
       "      <td>151</td>\n",
       "      <td>1347411.0</td>\n",
       "      <td>1</td>\n",
       "      <td>46446</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62766</th>\n",
       "      <td>62767</td>\n",
       "      <td>-596.731018</td>\n",
       "      <td>345.621002</td>\n",
       "      <td>-1795.500000</td>\n",
       "      <td>12</td>\n",
       "      <td>8</td>\n",
       "      <td>158</td>\n",
       "      <td>165095.0</td>\n",
       "      <td>1</td>\n",
       "      <td>30574</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72418</th>\n",
       "      <td>72419</td>\n",
       "      <td>-361.653992</td>\n",
       "      <td>2.315940</td>\n",
       "      <td>-886.200012</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>112</td>\n",
       "      <td>5919807.0</td>\n",
       "      <td>1</td>\n",
       "      <td>35623</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72881</th>\n",
       "      <td>72882</td>\n",
       "      <td>-356.710999</td>\n",
       "      <td>-43.158298</td>\n",
       "      <td>-606.119995</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>226</td>\n",
       "      <td>18523733.0</td>\n",
       "      <td>1</td>\n",
       "      <td>35966</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72907</th>\n",
       "      <td>72908</td>\n",
       "      <td>-271.240997</td>\n",
       "      <td>-238.298004</td>\n",
       "      <td>-578.099976</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>231</td>\n",
       "      <td>1864044.0</td>\n",
       "      <td>1</td>\n",
       "      <td>30708</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74722</th>\n",
       "      <td>74723</td>\n",
       "      <td>-29.566700</td>\n",
       "      <td>-355.957001</td>\n",
       "      <td>-60.400002</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>518</td>\n",
       "      <td>18525541.0</td>\n",
       "      <td>1</td>\n",
       "      <td>35823</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76207</th>\n",
       "      <td>76208</td>\n",
       "      <td>54.579601</td>\n",
       "      <td>358.360992</td>\n",
       "      <td>216.199997</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>713</td>\n",
       "      <td>18526747.0</td>\n",
       "      <td>1</td>\n",
       "      <td>36490</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76360</th>\n",
       "      <td>76361</td>\n",
       "      <td>-298.911987</td>\n",
       "      <td>-194.498001</td>\n",
       "      <td>291.600006</td>\n",
       "      <td>13</td>\n",
       "      <td>4</td>\n",
       "      <td>734</td>\n",
       "      <td>2552722.0</td>\n",
       "      <td>1</td>\n",
       "      <td>37581</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79306</th>\n",
       "      <td>79307</td>\n",
       "      <td>-495.473999</td>\n",
       "      <td>57.894001</td>\n",
       "      <td>-746.000000</td>\n",
       "      <td>13</td>\n",
       "      <td>6</td>\n",
       "      <td>311</td>\n",
       "      <td>18521714.0</td>\n",
       "      <td>1</td>\n",
       "      <td>38948</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79962</th>\n",
       "      <td>79963</td>\n",
       "      <td>430.200989</td>\n",
       "      <td>-255.901001</td>\n",
       "      <td>-368.200012</td>\n",
       "      <td>13</td>\n",
       "      <td>6</td>\n",
       "      <td>501</td>\n",
       "      <td>18526105.0</td>\n",
       "      <td>1</td>\n",
       "      <td>39158</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82254</th>\n",
       "      <td>82255</td>\n",
       "      <td>-363.391998</td>\n",
       "      <td>-339.842987</td>\n",
       "      <td>327.600006</td>\n",
       "      <td>13</td>\n",
       "      <td>6</td>\n",
       "      <td>1024</td>\n",
       "      <td>6737515.0</td>\n",
       "      <td>1</td>\n",
       "      <td>39734</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>84339</th>\n",
       "      <td>84340</td>\n",
       "      <td>155.970993</td>\n",
       "      <td>-642.041016</td>\n",
       "      <td>-1078.199951</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>30</td>\n",
       "      <td>18525952.0</td>\n",
       "      <td>1</td>\n",
       "      <td>42108</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>85129</th>\n",
       "      <td>85130</td>\n",
       "      <td>-404.596008</td>\n",
       "      <td>-526.367004</td>\n",
       "      <td>-608.143005</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>423</td>\n",
       "      <td>18525680.0</td>\n",
       "      <td>1</td>\n",
       "      <td>41325</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86527</th>\n",
       "      <td>86528</td>\n",
       "      <td>281.730988</td>\n",
       "      <td>-596.346008</td>\n",
       "      <td>-131.199997</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>951</td>\n",
       "      <td>3946074.0</td>\n",
       "      <td>1</td>\n",
       "      <td>41215</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89115</th>\n",
       "      <td>89116</td>\n",
       "      <td>-121.209999</td>\n",
       "      <td>650.249023</td>\n",
       "      <td>867.062012</td>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>1916</td>\n",
       "      <td>18526114.0</td>\n",
       "      <td>1</td>\n",
       "      <td>42495</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89755</th>\n",
       "      <td>89756</td>\n",
       "      <td>-230.985992</td>\n",
       "      <td>-497.566010</td>\n",
       "      <td>1218.500000</td>\n",
       "      <td>14</td>\n",
       "      <td>2</td>\n",
       "      <td>31</td>\n",
       "      <td>18524018.0</td>\n",
       "      <td>1</td>\n",
       "      <td>31318</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100122</th>\n",
       "      <td>100123</td>\n",
       "      <td>-592.258972</td>\n",
       "      <td>-565.510986</td>\n",
       "      <td>-737.200012</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>375</td>\n",
       "      <td>4025041.0</td>\n",
       "      <td>1</td>\n",
       "      <td>49452</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100497</th>\n",
       "      <td>100498</td>\n",
       "      <td>-33.920502</td>\n",
       "      <td>822.763977</td>\n",
       "      <td>-612.599976</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>571</td>\n",
       "      <td>3906428.0</td>\n",
       "      <td>1</td>\n",
       "      <td>49413</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>101337</th>\n",
       "      <td>101338</td>\n",
       "      <td>-732.999023</td>\n",
       "      <td>360.893005</td>\n",
       "      <td>-260.399994</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>952</td>\n",
       "      <td>7227919.0</td>\n",
       "      <td>1</td>\n",
       "      <td>49306</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>103330</th>\n",
       "      <td>103331</td>\n",
       "      <td>-600.182983</td>\n",
       "      <td>558.385010</td>\n",
       "      <td>542.000000</td>\n",
       "      <td>17</td>\n",
       "      <td>2</td>\n",
       "      <td>1906</td>\n",
       "      <td>18524913.0</td>\n",
       "      <td>1</td>\n",
       "      <td>49968</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104496</th>\n",
       "      <td>104497</td>\n",
       "      <td>126.666000</td>\n",
       "      <td>-1010.719971</td>\n",
       "      <td>-1046.199951</td>\n",
       "      <td>17</td>\n",
       "      <td>4</td>\n",
       "      <td>42</td>\n",
       "      <td>18525458.0</td>\n",
       "      <td>1</td>\n",
       "      <td>51786</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104709</th>\n",
       "      <td>104710</td>\n",
       "      <td>711.646973</td>\n",
       "      <td>-730.331970</td>\n",
       "      <td>-954.000000</td>\n",
       "      <td>17</td>\n",
       "      <td>4</td>\n",
       "      <td>209</td>\n",
       "      <td>18525977.0</td>\n",
       "      <td>1</td>\n",
       "      <td>50721</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>106893</th>\n",
       "      <td>106894</td>\n",
       "      <td>1003.309998</td>\n",
       "      <td>204.781998</td>\n",
       "      <td>-48.599998</td>\n",
       "      <td>17</td>\n",
       "      <td>4</td>\n",
       "      <td>1601</td>\n",
       "      <td>5489690.0</td>\n",
       "      <td>1</td>\n",
       "      <td>52238</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>109465</th>\n",
       "      <td>109466</td>\n",
       "      <td>-572.379028</td>\n",
       "      <td>-761.374023</td>\n",
       "      <td>1222.500000</td>\n",
       "      <td>18</td>\n",
       "      <td>2</td>\n",
       "      <td>15</td>\n",
       "      <td>18522829.0</td>\n",
       "      <td>1</td>\n",
       "      <td>52537</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111716</th>\n",
       "      <td>111717</td>\n",
       "      <td>-522.752014</td>\n",
       "      <td>-656.690002</td>\n",
       "      <td>2947.500000</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>14</td>\n",
       "      <td>18524172.0</td>\n",
       "      <td>1</td>\n",
       "      <td>47105</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>111919</th>\n",
       "      <td>111920</td>\n",
       "      <td>616.473022</td>\n",
       "      <td>455.427002</td>\n",
       "      <td>2947.500000</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>58</td>\n",
       "      <td>18523786.0</td>\n",
       "      <td>1</td>\n",
       "      <td>53218</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        hit_id            x            y            z  volume_id  layer_id  \\\n",
       "32648    32649    62.141201    98.091499  -128.850006          8         6   \n",
       "33495    33496   -27.344801   114.413002   -51.299801          8         6   \n",
       "33551    33552   -67.969498    95.379700   -17.507700          8         6   \n",
       "34974    34975  -110.622002    36.811901   115.222000          8         6   \n",
       "41705    41706   -60.984100  -160.223999   369.039001          8         8   \n",
       "42091    42092   171.774002     6.021760   468.938995          8         8   \n",
       "58239    58240  -560.581970   -30.300100 -2948.500000         12         2   \n",
       "59251    59252  -526.934998   441.579010 -2945.500000         12         2   \n",
       "61271    61272   292.415009   534.953003 -2148.500000         12         6   \n",
       "61513    61514  -324.071991   266.165009 -2151.500000         12         6   \n",
       "62766    62767  -596.731018   345.621002 -1795.500000         12         8   \n",
       "72418    72419  -361.653992     2.315940  -886.200012         13         4   \n",
       "72881    72882  -356.710999   -43.158298  -606.119995         13         4   \n",
       "72907    72908  -271.240997  -238.298004  -578.099976         13         4   \n",
       "74722    74723   -29.566700  -355.957001   -60.400002         13         4   \n",
       "76207    76208    54.579601   358.360992   216.199997         13         4   \n",
       "76360    76361  -298.911987  -194.498001   291.600006         13         4   \n",
       "79306    79307  -495.473999    57.894001  -746.000000         13         6   \n",
       "79962    79963   430.200989  -255.901001  -368.200012         13         6   \n",
       "82254    82255  -363.391998  -339.842987   327.600006         13         6   \n",
       "84339    84340   155.970993  -642.041016 -1078.199951         13         8   \n",
       "85129    85130  -404.596008  -526.367004  -608.143005         13         8   \n",
       "86527    86528   281.730988  -596.346008  -131.199997         13         8   \n",
       "89115    89116  -121.209999   650.249023   867.062012         13         8   \n",
       "89755    89756  -230.985992  -497.566010  1218.500000         14         2   \n",
       "100122  100123  -592.258972  -565.510986  -737.200012         17         2   \n",
       "100497  100498   -33.920502   822.763977  -612.599976         17         2   \n",
       "101337  101338  -732.999023   360.893005  -260.399994         17         2   \n",
       "103330  103331  -600.182983   558.385010   542.000000         17         2   \n",
       "104496  104497   126.666000 -1010.719971 -1046.199951         17         4   \n",
       "104709  104710   711.646973  -730.331970  -954.000000         17         4   \n",
       "106893  106894  1003.309998   204.781998   -48.599998         17         4   \n",
       "109465  109466  -572.379028  -761.374023  1222.500000         18         2   \n",
       "111716  111717  -522.752014  -656.690002  2947.500000         18        12   \n",
       "111919  111920   616.473022   455.427002  2947.500000         18        12   \n",
       "\n",
       "        module_id    track_id  N1    nbr  \n",
       "32648         295  18524844.0   1  17112  \n",
       "33495         354   8719220.0   1  17360  \n",
       "33551         357   1575617.0   1  16835  \n",
       "34974         466  11582980.0   1   2567  \n",
       "41705         951   2789076.0   1  20502  \n",
       "42091        1054  18520929.0   1  18734  \n",
       "58239           3   6648030.0   1  30068  \n",
       "59251         152  18521324.0   1  30562  \n",
       "61271         116  18522229.0   1  31021  \n",
       "61513         151   1347411.0   1  46446  \n",
       "62766         158    165095.0   1  30574  \n",
       "72418         112   5919807.0   1  35623  \n",
       "72881         226  18523733.0   1  35966  \n",
       "72907         231   1864044.0   1  30708  \n",
       "74722         518  18525541.0   1  35823  \n",
       "76207         713  18526747.0   1  36490  \n",
       "76360         734   2552722.0   1  37581  \n",
       "79306         311  18521714.0   1  38948  \n",
       "79962         501  18526105.0   1  39158  \n",
       "82254        1024   6737515.0   1  39734  \n",
       "84339          30  18525952.0   1  42108  \n",
       "85129         423  18525680.0   1  41325  \n",
       "86527         951   3946074.0   1  41215  \n",
       "89115        1916  18526114.0   1  42495  \n",
       "89755          31  18524018.0   1  31318  \n",
       "100122        375   4025041.0   1  49452  \n",
       "100497        571   3906428.0   1  49413  \n",
       "101337        952   7227919.0   1  49306  \n",
       "103330       1906  18524913.0   1  49968  \n",
       "104496         42  18525458.0   1  51786  \n",
       "104709        209  18525977.0   1  50721  \n",
       "106893       1601   5489690.0   1  52238  \n",
       "109465         15  18522829.0   1  52537  \n",
       "111716         14  18524172.0   1  47105  \n",
       "111919         58  18523786.0   1  53218  "
      ]
     },
     "execution_count": 109,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "outliers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 110,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>hit_id</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>volume_id</th>\n",
       "      <th>layer_id</th>\n",
       "      <th>module_id</th>\n",
       "      <th>track_id</th>\n",
       "      <th>N1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "Empty DataFrame\n",
       "Columns: [hit_id, x, y, z, volume_id, layer_id, module_id, track_id, N1]\n",
       "Index: []"
      ]
     },
     "execution_count": 110,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "hits['N1'] = hits.groupby('track_id')['track_id'].transform('count')\n",
    "\n",
    "mask = ( (hits['N1'] == 1) | (hits['track_id'] == 0) )\n",
    "\n",
    "hits[mask]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "112098"
      ]
     },
     "execution_count": 91,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### EDA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [],
   "source": [
    "for event_id, hits, cells, particles, truth in load_dataset(\"./data/train_1\", skip=40, nevents=1):\n",
    "    # Track pattern recognition\n",
    "    model = Clusterer(eps=0.0035, dz0=0, stepdz=1.5e-7, stepeps=7e-8, num_loops=100, final_cluster=True, weight=58, weight_arr=weights_arr, max_size=20, step_z=0.000045, size_incr=0)\n",
    "    labels = model._init(hits)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>hit_id</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>volume_id</th>\n",
       "      <th>layer_id</th>\n",
       "      <th>module_id</th>\n",
       "      <th>s1</th>\n",
       "      <th>N1</th>\n",
       "      <th>stepped_z</th>\n",
       "      <th>...</th>\n",
       "      <th>z1</th>\n",
       "      <th>z2</th>\n",
       "      <th>z3</th>\n",
       "      <th>a1</th>\n",
       "      <th>x2</th>\n",
       "      <th>y2</th>\n",
       "      <th>sina1</th>\n",
       "      <th>cosa1</th>\n",
       "      <th>s2</th>\n",
       "      <th>N2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>-48.674801</td>\n",
       "      <td>-5.761740</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-30.654200</td>\n",
       "      <td>-0.999468</td>\n",
       "      <td>-0.032622</td>\n",
       "      <td>-2.981392</td>\n",
       "      <td>-0.032379</td>\n",
       "      <td>-0.003833</td>\n",
       "      <td>-0.159516</td>\n",
       "      <td>-0.987195</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>-75.548500</td>\n",
       "      <td>-9.682310</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>1822912</td>\n",
       "      <td>15</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-19.726601</td>\n",
       "      <td>-0.998718</td>\n",
       "      <td>-0.050693</td>\n",
       "      <td>-2.948250</td>\n",
       "      <td>-0.050217</td>\n",
       "      <td>-0.006436</td>\n",
       "      <td>-0.192141</td>\n",
       "      <td>-0.981367</td>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>-69.968399</td>\n",
       "      <td>-5.581860</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-21.406033</td>\n",
       "      <td>-0.998911</td>\n",
       "      <td>-0.046716</td>\n",
       "      <td>-3.001281</td>\n",
       "      <td>-0.046517</td>\n",
       "      <td>-0.003711</td>\n",
       "      <td>-0.139852</td>\n",
       "      <td>-0.990172</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>-35.324799</td>\n",
       "      <td>-2.704890</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-42.409836</td>\n",
       "      <td>-0.999722</td>\n",
       "      <td>-0.023579</td>\n",
       "      <td>-3.034546</td>\n",
       "      <td>-0.023504</td>\n",
       "      <td>-0.001800</td>\n",
       "      <td>-0.106842</td>\n",
       "      <td>-0.994276</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>-99.261803</td>\n",
       "      <td>-11.197500</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>5588395</td>\n",
       "      <td>4</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-15.041381</td>\n",
       "      <td>-0.997797</td>\n",
       "      <td>-0.066483</td>\n",
       "      <td>-2.942832</td>\n",
       "      <td>-0.065919</td>\n",
       "      <td>-0.007436</td>\n",
       "      <td>-0.197455</td>\n",
       "      <td>-0.980312</td>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>6</td>\n",
       "      <td>-67.558197</td>\n",
       "      <td>-6.094340</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>5</td>\n",
       "      <td>14</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-22.150208</td>\n",
       "      <td>-0.998982</td>\n",
       "      <td>-0.045146</td>\n",
       "      <td>-2.992965</td>\n",
       "      <td>-0.044918</td>\n",
       "      <td>-0.004052</td>\n",
       "      <td>-0.148081</td>\n",
       "      <td>-0.988975</td>\n",
       "      <td>5</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>7</td>\n",
       "      <td>-63.736301</td>\n",
       "      <td>3.384740</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>6</td>\n",
       "      <td>11</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-23.540590</td>\n",
       "      <td>-0.999099</td>\n",
       "      <td>-0.042480</td>\n",
       "      <td>3.143731</td>\n",
       "      <td>-0.042382</td>\n",
       "      <td>0.002251</td>\n",
       "      <td>-0.002139</td>\n",
       "      <td>-0.999998</td>\n",
       "      <td>6</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>8</td>\n",
       "      <td>-59.426899</td>\n",
       "      <td>-12.626500</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>4726944</td>\n",
       "      <td>12</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-24.731169</td>\n",
       "      <td>-0.999184</td>\n",
       "      <td>-0.040435</td>\n",
       "      <td>-2.879700</td>\n",
       "      <td>-0.039520</td>\n",
       "      <td>-0.008397</td>\n",
       "      <td>-0.258909</td>\n",
       "      <td>-0.965902</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>9</td>\n",
       "      <td>-78.164299</td>\n",
       "      <td>-1.813290</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>8</td>\n",
       "      <td>2</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-19.217215</td>\n",
       "      <td>-0.998649</td>\n",
       "      <td>-0.052037</td>\n",
       "      <td>-3.050773</td>\n",
       "      <td>-0.051952</td>\n",
       "      <td>-0.001205</td>\n",
       "      <td>-0.090695</td>\n",
       "      <td>-0.995879</td>\n",
       "      <td>8</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>10</td>\n",
       "      <td>-83.351303</td>\n",
       "      <td>0.286217</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>3127933</td>\n",
       "      <td>16</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-18.026058</td>\n",
       "      <td>-0.998465</td>\n",
       "      <td>-0.055475</td>\n",
       "      <td>3.210259</td>\n",
       "      <td>-0.055390</td>\n",
       "      <td>0.000190</td>\n",
       "      <td>-0.068612</td>\n",
       "      <td>-0.997643</td>\n",
       "      <td>9</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>11</td>\n",
       "      <td>-91.341797</td>\n",
       "      <td>-14.536400</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>5502020</td>\n",
       "      <td>2</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-16.244827</td>\n",
       "      <td>-0.998111</td>\n",
       "      <td>-0.061558</td>\n",
       "      <td>-2.903757</td>\n",
       "      <td>-0.060678</td>\n",
       "      <td>-0.009657</td>\n",
       "      <td>-0.235600</td>\n",
       "      <td>-0.971850</td>\n",
       "      <td>10</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>12</td>\n",
       "      <td>-66.740898</td>\n",
       "      <td>-10.744500</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>1422236</td>\n",
       "      <td>14</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-22.226318</td>\n",
       "      <td>-0.998989</td>\n",
       "      <td>-0.044992</td>\n",
       "      <td>-2.923513</td>\n",
       "      <td>-0.044375</td>\n",
       "      <td>-0.007144</td>\n",
       "      <td>-0.216355</td>\n",
       "      <td>-0.976315</td>\n",
       "      <td>11</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>13</td>\n",
       "      <td>-60.177299</td>\n",
       "      <td>-3.808200</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>12</td>\n",
       "      <td>12</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-24.918114</td>\n",
       "      <td>-0.999196</td>\n",
       "      <td>-0.040131</td>\n",
       "      <td>-3.026254</td>\n",
       "      <td>-0.040019</td>\n",
       "      <td>-0.002533</td>\n",
       "      <td>-0.115083</td>\n",
       "      <td>-0.993356</td>\n",
       "      <td>12</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>14</td>\n",
       "      <td>-96.132599</td>\n",
       "      <td>-7.941960</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>237016</td>\n",
       "      <td>12</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-15.576433</td>\n",
       "      <td>-0.997946</td>\n",
       "      <td>-0.064200</td>\n",
       "      <td>-2.975710</td>\n",
       "      <td>-0.063850</td>\n",
       "      <td>-0.005275</td>\n",
       "      <td>-0.165123</td>\n",
       "      <td>-0.986273</td>\n",
       "      <td>13</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>15</td>\n",
       "      <td>-93.830498</td>\n",
       "      <td>-15.685200</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>3210493</td>\n",
       "      <td>2</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-15.793812</td>\n",
       "      <td>-0.998002</td>\n",
       "      <td>-0.063316</td>\n",
       "      <td>-2.893654</td>\n",
       "      <td>-0.062325</td>\n",
       "      <td>-0.010419</td>\n",
       "      <td>-0.245406</td>\n",
       "      <td>-0.969420</td>\n",
       "      <td>14</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>16</td>\n",
       "      <td>-87.631104</td>\n",
       "      <td>-1.254320</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>15</td>\n",
       "      <td>15</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-17.144032</td>\n",
       "      <td>-0.998303</td>\n",
       "      <td>-0.058329</td>\n",
       "      <td>-3.051466</td>\n",
       "      <td>-0.058224</td>\n",
       "      <td>-0.000833</td>\n",
       "      <td>-0.090005</td>\n",
       "      <td>-0.995941</td>\n",
       "      <td>15</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>17</td>\n",
       "      <td>-82.778000</td>\n",
       "      <td>-13.260600</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>7803580</td>\n",
       "      <td>6</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-17.922501</td>\n",
       "      <td>-0.998447</td>\n",
       "      <td>-0.055796</td>\n",
       "      <td>-2.910231</td>\n",
       "      <td>-0.055008</td>\n",
       "      <td>-0.008812</td>\n",
       "      <td>-0.229303</td>\n",
       "      <td>-0.973355</td>\n",
       "      <td>16</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>18</td>\n",
       "      <td>-84.366203</td>\n",
       "      <td>-10.476600</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>2145544</td>\n",
       "      <td>14</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-17.673569</td>\n",
       "      <td>-0.998403</td>\n",
       "      <td>-0.056582</td>\n",
       "      <td>-2.944505</td>\n",
       "      <td>-0.056061</td>\n",
       "      <td>-0.006962</td>\n",
       "      <td>-0.195814</td>\n",
       "      <td>-0.980641</td>\n",
       "      <td>17</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>19</td>\n",
       "      <td>-36.901402</td>\n",
       "      <td>-9.198610</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>18</td>\n",
       "      <td>1</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-39.507759</td>\n",
       "      <td>-0.999680</td>\n",
       "      <td>-0.025311</td>\n",
       "      <td>-2.864422</td>\n",
       "      <td>-0.024552</td>\n",
       "      <td>-0.006120</td>\n",
       "      <td>-0.273636</td>\n",
       "      <td>-0.961833</td>\n",
       "      <td>18</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>20</td>\n",
       "      <td>-62.858200</td>\n",
       "      <td>-12.043900</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>2307828</td>\n",
       "      <td>12</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-23.476032</td>\n",
       "      <td>-0.999094</td>\n",
       "      <td>-0.042597</td>\n",
       "      <td>-2.896937</td>\n",
       "      <td>-0.041798</td>\n",
       "      <td>-0.008009</td>\n",
       "      <td>-0.242222</td>\n",
       "      <td>-0.970221</td>\n",
       "      <td>19</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>21</td>\n",
       "      <td>-51.928200</td>\n",
       "      <td>-0.367017</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>20</td>\n",
       "      <td>1</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-28.933544</td>\n",
       "      <td>-0.999403</td>\n",
       "      <td>-0.034562</td>\n",
       "      <td>-3.089627</td>\n",
       "      <td>-0.034540</td>\n",
       "      <td>-0.000244</td>\n",
       "      <td>-0.051943</td>\n",
       "      <td>-0.998650</td>\n",
       "      <td>20</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>22</td>\n",
       "      <td>-74.732399</td>\n",
       "      <td>-4.602630</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>21</td>\n",
       "      <td>1</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-20.067104</td>\n",
       "      <td>-0.998761</td>\n",
       "      <td>-0.049833</td>\n",
       "      <td>-3.015324</td>\n",
       "      <td>-0.049677</td>\n",
       "      <td>-0.003060</td>\n",
       "      <td>-0.125934</td>\n",
       "      <td>-0.992039</td>\n",
       "      <td>21</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>23</td>\n",
       "      <td>-103.884003</td>\n",
       "      <td>-12.836400</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2307831</td>\n",
       "      <td>13</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-14.349346</td>\n",
       "      <td>-0.997580</td>\n",
       "      <td>-0.069690</td>\n",
       "      <td>-2.928079</td>\n",
       "      <td>-0.068996</td>\n",
       "      <td>-0.008526</td>\n",
       "      <td>-0.211895</td>\n",
       "      <td>-0.977292</td>\n",
       "      <td>22</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>24</td>\n",
       "      <td>-155.029999</td>\n",
       "      <td>-3.939040</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2064563</td>\n",
       "      <td>15</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.685350</td>\n",
       "      <td>-0.994712</td>\n",
       "      <td>-0.103249</td>\n",
       "      <td>-2.981902</td>\n",
       "      <td>-0.102670</td>\n",
       "      <td>-0.002609</td>\n",
       "      <td>-0.159013</td>\n",
       "      <td>-0.987276</td>\n",
       "      <td>23</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>25</td>\n",
       "      <td>-122.973999</td>\n",
       "      <td>-8.463410</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2798314</td>\n",
       "      <td>10</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-12.185176</td>\n",
       "      <td>-0.996649</td>\n",
       "      <td>-0.082067</td>\n",
       "      <td>-2.966190</td>\n",
       "      <td>-0.081599</td>\n",
       "      <td>-0.005616</td>\n",
       "      <td>-0.174505</td>\n",
       "      <td>-0.984656</td>\n",
       "      <td>24</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>26</td>\n",
       "      <td>-164.014999</td>\n",
       "      <td>-10.461000</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>4812712</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.139155</td>\n",
       "      <td>-0.994067</td>\n",
       "      <td>-0.109419</td>\n",
       "      <td>-2.935565</td>\n",
       "      <td>-0.108550</td>\n",
       "      <td>-0.006923</td>\n",
       "      <td>-0.204573</td>\n",
       "      <td>-0.978851</td>\n",
       "      <td>25</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>27</td>\n",
       "      <td>-141.460007</td>\n",
       "      <td>0.942888</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>1983944</td>\n",
       "      <td>4</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.617637</td>\n",
       "      <td>-0.995594</td>\n",
       "      <td>-0.094183</td>\n",
       "      <td>3.257399</td>\n",
       "      <td>-0.093766</td>\n",
       "      <td>0.000625</td>\n",
       "      <td>-0.115548</td>\n",
       "      <td>-0.993302</td>\n",
       "      <td>26</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>28</td>\n",
       "      <td>-156.419006</td>\n",
       "      <td>-4.854170</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2064563</td>\n",
       "      <td>15</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.597821</td>\n",
       "      <td>-0.994616</td>\n",
       "      <td>-0.104190</td>\n",
       "      <td>-2.975054</td>\n",
       "      <td>-0.103579</td>\n",
       "      <td>-0.003214</td>\n",
       "      <td>-0.165770</td>\n",
       "      <td>-0.986164</td>\n",
       "      <td>27</td>\n",
       "      <td>9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>29</td>\n",
       "      <td>-150.740997</td>\n",
       "      <td>-9.847230</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>4898332</td>\n",
       "      <td>11</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.942947</td>\n",
       "      <td>-0.994980</td>\n",
       "      <td>-0.100574</td>\n",
       "      <td>-2.945559</td>\n",
       "      <td>-0.099856</td>\n",
       "      <td>-0.006523</td>\n",
       "      <td>-0.194781</td>\n",
       "      <td>-0.980847</td>\n",
       "      <td>28</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>30</td>\n",
       "      <td>-120.870003</td>\n",
       "      <td>-13.952400</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>29</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-12.344637</td>\n",
       "      <td>-0.996735</td>\n",
       "      <td>-0.081007</td>\n",
       "      <td>-2.921360</td>\n",
       "      <td>-0.080210</td>\n",
       "      <td>-0.009259</td>\n",
       "      <td>-0.218456</td>\n",
       "      <td>-0.975847</td>\n",
       "      <td>29</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124355</th>\n",
       "      <td>124356</td>\n",
       "      <td>-976.395996</td>\n",
       "      <td>245.266998</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>94</td>\n",
       "      <td>2334520</td>\n",
       "      <td>6</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>2.932759</td>\n",
       "      <td>0.946491</td>\n",
       "      <td>0.340976</td>\n",
       "      <td>3.778244</td>\n",
       "      <td>-0.313006</td>\n",
       "      <td>0.078626</td>\n",
       "      <td>-0.594506</td>\n",
       "      <td>-0.804091</td>\n",
       "      <td>90959</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124356</th>\n",
       "      <td>124357</td>\n",
       "      <td>-904.656982</td>\n",
       "      <td>230.473999</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>94</td>\n",
       "      <td>2648289</td>\n",
       "      <td>8</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.162641</td>\n",
       "      <td>0.953473</td>\n",
       "      <td>0.316191</td>\n",
       "      <td>3.709849</td>\n",
       "      <td>-0.292148</td>\n",
       "      <td>0.074429</td>\n",
       "      <td>-0.538163</td>\n",
       "      <td>-0.842841</td>\n",
       "      <td>90960</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124357</th>\n",
       "      <td>124358</td>\n",
       "      <td>-809.007019</td>\n",
       "      <td>322.921997</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>94</td>\n",
       "      <td>6640453</td>\n",
       "      <td>4</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.389486</td>\n",
       "      <td>0.959128</td>\n",
       "      <td>0.295030</td>\n",
       "      <td>3.524101</td>\n",
       "      <td>-0.262809</td>\n",
       "      <td>0.104902</td>\n",
       "      <td>-0.373249</td>\n",
       "      <td>-0.927731</td>\n",
       "      <td>90942</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124358</th>\n",
       "      <td>124359</td>\n",
       "      <td>-779.510986</td>\n",
       "      <td>165.393005</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>95</td>\n",
       "      <td>7766413</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.698869</td>\n",
       "      <td>0.965343</td>\n",
       "      <td>0.270353</td>\n",
       "      <td>3.629106</td>\n",
       "      <td>-0.255300</td>\n",
       "      <td>0.054168</td>\n",
       "      <td>-0.468430</td>\n",
       "      <td>-0.883501</td>\n",
       "      <td>90961</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124359</th>\n",
       "      <td>124360</td>\n",
       "      <td>-791.286987</td>\n",
       "      <td>106.193001</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>95</td>\n",
       "      <td>7605044</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.691841</td>\n",
       "      <td>0.965218</td>\n",
       "      <td>0.270868</td>\n",
       "      <td>3.706116</td>\n",
       "      <td>-0.259123</td>\n",
       "      <td>0.034775</td>\n",
       "      <td>-0.535013</td>\n",
       "      <td>-0.844844</td>\n",
       "      <td>90962</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124360</th>\n",
       "      <td>124361</td>\n",
       "      <td>-780.950989</td>\n",
       "      <td>158.154999</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>95</td>\n",
       "      <td>7766293</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.699145</td>\n",
       "      <td>0.965348</td>\n",
       "      <td>0.270333</td>\n",
       "      <td>3.638314</td>\n",
       "      <td>-0.255773</td>\n",
       "      <td>0.051798</td>\n",
       "      <td>-0.476546</td>\n",
       "      <td>-0.879150</td>\n",
       "      <td>90963</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124361</th>\n",
       "      <td>124362</td>\n",
       "      <td>-830.892029</td>\n",
       "      <td>226.938995</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>95</td>\n",
       "      <td>6751146</td>\n",
       "      <td>5</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.422043</td>\n",
       "      <td>0.959856</td>\n",
       "      <td>0.292223</td>\n",
       "      <td>3.628619</td>\n",
       "      <td>-0.270581</td>\n",
       "      <td>0.073903</td>\n",
       "      <td>-0.468000</td>\n",
       "      <td>-0.883728</td>\n",
       "      <td>90946</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124362</th>\n",
       "      <td>124363</td>\n",
       "      <td>-746.534973</td>\n",
       "      <td>171.251999</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>95</td>\n",
       "      <td>8402232</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.848280</td>\n",
       "      <td>0.967856</td>\n",
       "      <td>0.259856</td>\n",
       "      <td>3.585336</td>\n",
       "      <td>-0.245136</td>\n",
       "      <td>0.056233</td>\n",
       "      <td>-0.429324</td>\n",
       "      <td>-0.903151</td>\n",
       "      <td>53912</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124363</th>\n",
       "      <td>124364</td>\n",
       "      <td>-813.125000</td>\n",
       "      <td>209.641006</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>95</td>\n",
       "      <td>7785159</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.510113</td>\n",
       "      <td>0.961733</td>\n",
       "      <td>0.284891</td>\n",
       "      <td>3.623775</td>\n",
       "      <td>-0.265313</td>\n",
       "      <td>0.068403</td>\n",
       "      <td>-0.463713</td>\n",
       "      <td>-0.885985</td>\n",
       "      <td>90964</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124364</th>\n",
       "      <td>124365</td>\n",
       "      <td>-991.192017</td>\n",
       "      <td>219.682007</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>8148851</td>\n",
       "      <td>1</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>2.911117</td>\n",
       "      <td>0.945756</td>\n",
       "      <td>0.343511</td>\n",
       "      <td>3.813816</td>\n",
       "      <td>-0.317181</td>\n",
       "      <td>0.070298</td>\n",
       "      <td>-0.622727</td>\n",
       "      <td>-0.782439</td>\n",
       "      <td>90965</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124365</th>\n",
       "      <td>124366</td>\n",
       "      <td>-950.107971</td>\n",
       "      <td>157.542999</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>7605020</td>\n",
       "      <td>1</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.068792</td>\n",
       "      <td>0.950793</td>\n",
       "      <td>0.325861</td>\n",
       "      <td>3.821214</td>\n",
       "      <td>-0.305653</td>\n",
       "      <td>0.050682</td>\n",
       "      <td>-0.628498</td>\n",
       "      <td>-0.777811</td>\n",
       "      <td>90966</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124366</th>\n",
       "      <td>124367</td>\n",
       "      <td>-914.064026</td>\n",
       "      <td>179.987000</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>7427479</td>\n",
       "      <td>2</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.172440</td>\n",
       "      <td>0.953740</td>\n",
       "      <td>0.315215</td>\n",
       "      <td>3.763165</td>\n",
       "      <td>-0.294969</td>\n",
       "      <td>0.058082</td>\n",
       "      <td>-0.582314</td>\n",
       "      <td>-0.812964</td>\n",
       "      <td>90967</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124367</th>\n",
       "      <td>124368</td>\n",
       "      <td>-975.681030</td>\n",
       "      <td>245.492996</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>2334520</td>\n",
       "      <td>6</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>2.937601</td>\n",
       "      <td>0.946653</td>\n",
       "      <td>0.340414</td>\n",
       "      <td>3.777285</td>\n",
       "      <td>-0.312513</td>\n",
       "      <td>0.078632</td>\n",
       "      <td>-0.593734</td>\n",
       "      <td>-0.804661</td>\n",
       "      <td>90959</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124368</th>\n",
       "      <td>124369</td>\n",
       "      <td>-863.390015</td>\n",
       "      <td>168.121002</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>8330865</td>\n",
       "      <td>1</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.360022</td>\n",
       "      <td>0.958453</td>\n",
       "      <td>0.297617</td>\n",
       "      <td>3.719129</td>\n",
       "      <td>-0.279993</td>\n",
       "      <td>0.054521</td>\n",
       "      <td>-0.545961</td>\n",
       "      <td>-0.837810</td>\n",
       "      <td>90968</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124369</th>\n",
       "      <td>124370</td>\n",
       "      <td>-901.502991</td>\n",
       "      <td>134.832993</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>7785186</td>\n",
       "      <td>2</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.242345</td>\n",
       "      <td>0.955584</td>\n",
       "      <td>0.308419</td>\n",
       "      <td>3.791294</td>\n",
       "      <td>-0.291478</td>\n",
       "      <td>0.043595</td>\n",
       "      <td>-0.604949</td>\n",
       "      <td>-0.796264</td>\n",
       "      <td>90969</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124370</th>\n",
       "      <td>124371</td>\n",
       "      <td>-904.025024</td>\n",
       "      <td>232.617996</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>2648289</td>\n",
       "      <td>8</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.166128</td>\n",
       "      <td>0.953568</td>\n",
       "      <td>0.315843</td>\n",
       "      <td>3.707385</td>\n",
       "      <td>-0.291677</td>\n",
       "      <td>0.075052</td>\n",
       "      <td>-0.536085</td>\n",
       "      <td>-0.844164</td>\n",
       "      <td>90960</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124371</th>\n",
       "      <td>124372</td>\n",
       "      <td>-910.455017</td>\n",
       "      <td>198.906006</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>7622978</td>\n",
       "      <td>2</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.171374</td>\n",
       "      <td>0.953711</td>\n",
       "      <td>0.315321</td>\n",
       "      <td>3.742775</td>\n",
       "      <td>-0.293795</td>\n",
       "      <td>0.064185</td>\n",
       "      <td>-0.565618</td>\n",
       "      <td>-0.824668</td>\n",
       "      <td>90970</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124372</th>\n",
       "      <td>124373</td>\n",
       "      <td>-933.367981</td>\n",
       "      <td>245.300995</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>96</td>\n",
       "      <td>1757465</td>\n",
       "      <td>10</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.062487</td>\n",
       "      <td>0.950605</td>\n",
       "      <td>0.326532</td>\n",
       "      <td>3.730296</td>\n",
       "      <td>-0.300208</td>\n",
       "      <td>0.078899</td>\n",
       "      <td>-0.555283</td>\n",
       "      <td>-0.831661</td>\n",
       "      <td>90954</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124373</th>\n",
       "      <td>124374</td>\n",
       "      <td>-795.552979</td>\n",
       "      <td>45.108200</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>97</td>\n",
       "      <td>7566778</td>\n",
       "      <td>1</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.695258</td>\n",
       "      <td>0.965279</td>\n",
       "      <td>0.270617</td>\n",
       "      <td>3.781511</td>\n",
       "      <td>-0.260802</td>\n",
       "      <td>0.014788</td>\n",
       "      <td>-0.597130</td>\n",
       "      <td>-0.802144</td>\n",
       "      <td>90971</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124374</th>\n",
       "      <td>124375</td>\n",
       "      <td>-817.885986</td>\n",
       "      <td>22.400299</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>97</td>\n",
       "      <td>5783904</td>\n",
       "      <td>7</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.598780</td>\n",
       "      <td>0.963495</td>\n",
       "      <td>0.277872</td>\n",
       "      <td>3.829668</td>\n",
       "      <td>-0.267628</td>\n",
       "      <td>0.007330</td>\n",
       "      <td>-0.635051</td>\n",
       "      <td>-0.772470</td>\n",
       "      <td>90972</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124375</th>\n",
       "      <td>124376</td>\n",
       "      <td>-865.489990</td>\n",
       "      <td>91.181000</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>97</td>\n",
       "      <td>7966217</td>\n",
       "      <td>6</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.383390</td>\n",
       "      <td>0.958990</td>\n",
       "      <td>0.295562</td>\n",
       "      <td>3.798213</td>\n",
       "      <td>-0.281881</td>\n",
       "      <td>0.029697</td>\n",
       "      <td>-0.610443</td>\n",
       "      <td>-0.792060</td>\n",
       "      <td>90973</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124376</th>\n",
       "      <td>124377</td>\n",
       "      <td>-897.953003</td>\n",
       "      <td>72.925102</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>97</td>\n",
       "      <td>8128579</td>\n",
       "      <td>3</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.268359</td>\n",
       "      <td>0.956242</td>\n",
       "      <td>0.305964</td>\n",
       "      <td>3.849301</td>\n",
       "      <td>-0.291616</td>\n",
       "      <td>0.023683</td>\n",
       "      <td>-0.650094</td>\n",
       "      <td>-0.759854</td>\n",
       "      <td>90974</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124377</th>\n",
       "      <td>124378</td>\n",
       "      <td>-880.408020</td>\n",
       "      <td>22.589800</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>97</td>\n",
       "      <td>1918612</td>\n",
       "      <td>10</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.343367</td>\n",
       "      <td>0.958063</td>\n",
       "      <td>0.299100</td>\n",
       "      <td>3.886759</td>\n",
       "      <td>-0.286462</td>\n",
       "      <td>0.007350</td>\n",
       "      <td>-0.678094</td>\n",
       "      <td>-0.734975</td>\n",
       "      <td>90975</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124378</th>\n",
       "      <td>124379</td>\n",
       "      <td>-781.200012</td>\n",
       "      <td>105.078003</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>97</td>\n",
       "      <td>6358094</td>\n",
       "      <td>1</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.735554</td>\n",
       "      <td>0.965986</td>\n",
       "      <td>0.267698</td>\n",
       "      <td>3.696844</td>\n",
       "      <td>-0.256285</td>\n",
       "      <td>0.034472</td>\n",
       "      <td>-0.527157</td>\n",
       "      <td>-0.849768</td>\n",
       "      <td>90976</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124379</th>\n",
       "      <td>124380</td>\n",
       "      <td>-786.262024</td>\n",
       "      <td>27.843599</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>97</td>\n",
       "      <td>1924892</td>\n",
       "      <td>4</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.742583</td>\n",
       "      <td>0.966108</td>\n",
       "      <td>0.267195</td>\n",
       "      <td>3.793843</td>\n",
       "      <td>-0.257978</td>\n",
       "      <td>0.009136</td>\n",
       "      <td>-0.606977</td>\n",
       "      <td>-0.794720</td>\n",
       "      <td>90977</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124380</th>\n",
       "      <td>124381</td>\n",
       "      <td>-890.322998</td>\n",
       "      <td>22.828899</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>98</td>\n",
       "      <td>1918612</td>\n",
       "      <td>10</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.315118</td>\n",
       "      <td>0.957391</td>\n",
       "      <td>0.301648</td>\n",
       "      <td>3.895570</td>\n",
       "      <td>-0.288701</td>\n",
       "      <td>0.007403</td>\n",
       "      <td>-0.684544</td>\n",
       "      <td>-0.728972</td>\n",
       "      <td>73690</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124381</th>\n",
       "      <td>124382</td>\n",
       "      <td>-897.323975</td>\n",
       "      <td>77.195702</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>98</td>\n",
       "      <td>7908973</td>\n",
       "      <td>9</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.278225</td>\n",
       "      <td>0.956488</td>\n",
       "      <td>0.305043</td>\n",
       "      <td>3.844278</td>\n",
       "      <td>-0.290696</td>\n",
       "      <td>0.025008</td>\n",
       "      <td>-0.646269</td>\n",
       "      <td>-0.763110</td>\n",
       "      <td>90978</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124382</th>\n",
       "      <td>124383</td>\n",
       "      <td>-865.239014</td>\n",
       "      <td>90.266998</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>98</td>\n",
       "      <td>7966217</td>\n",
       "      <td>6</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.393927</td>\n",
       "      <td>0.959229</td>\n",
       "      <td>0.294644</td>\n",
       "      <td>3.798922</td>\n",
       "      <td>-0.281105</td>\n",
       "      <td>0.029327</td>\n",
       "      <td>-0.611005</td>\n",
       "      <td>-0.791627</td>\n",
       "      <td>90973</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124383</th>\n",
       "      <td>124384</td>\n",
       "      <td>-971.604004</td>\n",
       "      <td>55.957901</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>98</td>\n",
       "      <td>4401397</td>\n",
       "      <td>2</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.033758</td>\n",
       "      <td>0.949735</td>\n",
       "      <td>0.329624</td>\n",
       "      <td>3.937011</td>\n",
       "      <td>-0.312538</td>\n",
       "      <td>0.018000</td>\n",
       "      <td>-0.714156</td>\n",
       "      <td>-0.699986</td>\n",
       "      <td>90979</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124384</th>\n",
       "      <td>124385</td>\n",
       "      <td>-906.614014</td>\n",
       "      <td>95.154503</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>98</td>\n",
       "      <td>6907385</td>\n",
       "      <td>1</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.238828</td>\n",
       "      <td>0.955494</td>\n",
       "      <td>0.308754</td>\n",
       "      <td>3.835241</td>\n",
       "      <td>-0.293401</td>\n",
       "      <td>0.030794</td>\n",
       "      <td>-0.639347</td>\n",
       "      <td>-0.768918</td>\n",
       "      <td>90980</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>124385 rows × 23 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        hit_id           x           y       z  volume_id  layer_id  \\\n",
       "0            1  -48.674801   -5.761740 -1502.5          7         2   \n",
       "1            2  -75.548500   -9.682310 -1502.5          7         2   \n",
       "2            3  -69.968399   -5.581860 -1502.5          7         2   \n",
       "3            4  -35.324799   -2.704890 -1502.5          7         2   \n",
       "4            5  -99.261803  -11.197500 -1502.5          7         2   \n",
       "5            6  -67.558197   -6.094340 -1502.5          7         2   \n",
       "6            7  -63.736301    3.384740 -1502.5          7         2   \n",
       "7            8  -59.426899  -12.626500 -1502.5          7         2   \n",
       "8            9  -78.164299   -1.813290 -1502.5          7         2   \n",
       "9           10  -83.351303    0.286217 -1502.5          7         2   \n",
       "10          11  -91.341797  -14.536400 -1502.5          7         2   \n",
       "11          12  -66.740898  -10.744500 -1502.5          7         2   \n",
       "12          13  -60.177299   -3.808200 -1502.5          7         2   \n",
       "13          14  -96.132599   -7.941960 -1502.5          7         2   \n",
       "14          15  -93.830498  -15.685200 -1502.5          7         2   \n",
       "15          16  -87.631104   -1.254320 -1502.5          7         2   \n",
       "16          17  -82.778000  -13.260600 -1502.5          7         2   \n",
       "17          18  -84.366203  -10.476600 -1502.5          7         2   \n",
       "18          19  -36.901402   -9.198610 -1502.5          7         2   \n",
       "19          20  -62.858200  -12.043900 -1502.5          7         2   \n",
       "20          21  -51.928200   -0.367017 -1502.5          7         2   \n",
       "21          22  -74.732399   -4.602630 -1502.5          7         2   \n",
       "22          23 -103.884003  -12.836400 -1502.0          7         2   \n",
       "23          24 -155.029999   -3.939040 -1502.0          7         2   \n",
       "24          25 -122.973999   -8.463410 -1502.0          7         2   \n",
       "25          26 -164.014999  -10.461000 -1502.0          7         2   \n",
       "26          27 -141.460007    0.942888 -1502.0          7         2   \n",
       "27          28 -156.419006   -4.854170 -1502.0          7         2   \n",
       "28          29 -150.740997   -9.847230 -1502.0          7         2   \n",
       "29          30 -120.870003  -13.952400 -1502.0          7         2   \n",
       "...        ...         ...         ...     ...        ...       ...   \n",
       "124355  124356 -976.395996  245.266998  2952.5         18        12   \n",
       "124356  124357 -904.656982  230.473999  2952.5         18        12   \n",
       "124357  124358 -809.007019  322.921997  2952.5         18        12   \n",
       "124358  124359 -779.510986  165.393005  2947.5         18        12   \n",
       "124359  124360 -791.286987  106.193001  2947.5         18        12   \n",
       "124360  124361 -780.950989  158.154999  2947.5         18        12   \n",
       "124361  124362 -830.892029  226.938995  2947.5         18        12   \n",
       "124362  124363 -746.534973  171.251999  2947.5         18        12   \n",
       "124363  124364 -813.125000  209.641006  2947.5         18        12   \n",
       "124364  124365 -991.192017  219.682007  2955.5         18        12   \n",
       "124365  124366 -950.107971  157.542999  2955.5         18        12   \n",
       "124366  124367 -914.064026  179.987000  2955.5         18        12   \n",
       "124367  124368 -975.681030  245.492996  2955.5         18        12   \n",
       "124368  124369 -863.390015  168.121002  2955.5         18        12   \n",
       "124369  124370 -901.502991  134.832993  2955.5         18        12   \n",
       "124370  124371 -904.025024  232.617996  2955.5         18        12   \n",
       "124371  124372 -910.455017  198.906006  2955.5         18        12   \n",
       "124372  124373 -933.367981  245.300995  2955.5         18        12   \n",
       "124373  124374 -795.552979   45.108200  2944.5         18        12   \n",
       "124374  124375 -817.885986   22.400299  2944.5         18        12   \n",
       "124375  124376 -865.489990   91.181000  2944.5         18        12   \n",
       "124376  124377 -897.953003   72.925102  2944.5         18        12   \n",
       "124377  124378 -880.408020   22.589800  2944.5         18        12   \n",
       "124378  124379 -781.200012  105.078003  2944.5         18        12   \n",
       "124379  124380 -786.262024   27.843599  2944.5         18        12   \n",
       "124380  124381 -890.322998   22.828899  2952.5         18        12   \n",
       "124381  124382 -897.323975   77.195702  2952.5         18        12   \n",
       "124382  124383 -865.239014   90.266998  2952.5         18        12   \n",
       "124383  124384 -971.604004   55.957901  2952.5         18        12   \n",
       "124384  124385 -906.614014   95.154503  2952.5         18        12   \n",
       "\n",
       "        module_id       s1  N1    stepped_z ...         z1        z2  \\\n",
       "0               1        0   2 -1502.495605 ... -30.654200 -0.999468   \n",
       "1               1  1822912  15 -1502.495605 ... -19.726601 -0.998718   \n",
       "2               1        2   1 -1502.495605 ... -21.406033 -0.998911   \n",
       "3               1        3   1 -1502.495605 ... -42.409836 -0.999722   \n",
       "4               1  5588395   4 -1502.495605 ... -15.041381 -0.997797   \n",
       "5               1        5  14 -1502.495605 ... -22.150208 -0.998982   \n",
       "6               1        6  11 -1502.495605 ... -23.540590 -0.999099   \n",
       "7               1  4726944  12 -1502.495605 ... -24.731169 -0.999184   \n",
       "8               1        8   2 -1502.495605 ... -19.217215 -0.998649   \n",
       "9               1  3127933  16 -1502.495605 ... -18.026058 -0.998465   \n",
       "10              1  5502020   2 -1502.495605 ... -16.244827 -0.998111   \n",
       "11              1  1422236  14 -1502.495605 ... -22.226318 -0.998989   \n",
       "12              1       12  12 -1502.495605 ... -24.918114 -0.999196   \n",
       "13              1   237016  12 -1502.495605 ... -15.576433 -0.997946   \n",
       "14              1  3210493   2 -1502.495605 ... -15.793812 -0.998002   \n",
       "15              1       15  15 -1502.495605 ... -17.144032 -0.998303   \n",
       "16              1  7803580   6 -1502.495605 ... -17.922501 -0.998447   \n",
       "17              1  2145544  14 -1502.495605 ... -17.673569 -0.998403   \n",
       "18              1       18   1 -1502.495605 ... -39.507759 -0.999680   \n",
       "19              1  2307828  12 -1502.495605 ... -23.476032 -0.999094   \n",
       "20              1       20   1 -1502.495605 ... -28.933544 -0.999403   \n",
       "21              1       21   1 -1502.495605 ... -20.067104 -0.998761   \n",
       "22              2  2307831  13 -1501.995605 ... -14.349346 -0.997580   \n",
       "23              2  2064563  15 -1501.995605 ...  -9.685350 -0.994712   \n",
       "24              2  2798314  10 -1501.995605 ... -12.185176 -0.996649   \n",
       "25              2  4812712   1 -1501.995605 ...  -9.139155 -0.994067   \n",
       "26              2  1983944   4 -1501.995605 ... -10.617637 -0.995594   \n",
       "27              2  2064563  15 -1501.995605 ...  -9.597821 -0.994616   \n",
       "28              2  4898332  11 -1501.995605 ...  -9.942947 -0.994980   \n",
       "29              2       29   1 -1501.995605 ... -12.344637 -0.996735   \n",
       "...           ...      ...  ..          ... ...        ...       ...   \n",
       "124355         94  2334520   6  2952.504395 ...   2.932759  0.946491   \n",
       "124356         94  2648289   8  2952.504395 ...   3.162641  0.953473   \n",
       "124357         94  6640453   4  2952.504395 ...   3.389486  0.959128   \n",
       "124358         95  7766413   1  2947.504395 ...   3.698869  0.965343   \n",
       "124359         95  7605044   1  2947.504395 ...   3.691841  0.965218   \n",
       "124360         95  7766293   1  2947.504395 ...   3.699145  0.965348   \n",
       "124361         95  6751146   5  2947.504395 ...   3.422043  0.959856   \n",
       "124362         95  8402232   1  2947.504395 ...   3.848280  0.967856   \n",
       "124363         95  7785159   1  2947.504395 ...   3.510113  0.961733   \n",
       "124364         96  8148851   1  2955.504395 ...   2.911117  0.945756   \n",
       "124365         96  7605020   1  2955.504395 ...   3.068792  0.950793   \n",
       "124366         96  7427479   2  2955.504395 ...   3.172440  0.953740   \n",
       "124367         96  2334520   6  2955.504395 ...   2.937601  0.946653   \n",
       "124368         96  8330865   1  2955.504395 ...   3.360022  0.958453   \n",
       "124369         96  7785186   2  2955.504395 ...   3.242345  0.955584   \n",
       "124370         96  2648289   8  2955.504395 ...   3.166128  0.953568   \n",
       "124371         96  7622978   2  2955.504395 ...   3.171374  0.953711   \n",
       "124372         96  1757465  10  2955.504395 ...   3.062487  0.950605   \n",
       "124373         97  7566778   1  2944.504395 ...   3.695258  0.965279   \n",
       "124374         97  5783904   7  2944.504395 ...   3.598780  0.963495   \n",
       "124375         97  7966217   6  2944.504395 ...   3.383390  0.958990   \n",
       "124376         97  8128579   3  2944.504395 ...   3.268359  0.956242   \n",
       "124377         97  1918612  10  2944.504395 ...   3.343367  0.958063   \n",
       "124378         97  6358094   1  2944.504395 ...   3.735554  0.965986   \n",
       "124379         97  1924892   4  2944.504395 ...   3.742583  0.966108   \n",
       "124380         98  1918612  10  2952.504395 ...   3.315118  0.957391   \n",
       "124381         98  7908973   9  2952.504395 ...   3.278225  0.956488   \n",
       "124382         98  7966217   6  2952.504395 ...   3.393927  0.959229   \n",
       "124383         98  4401397   2  2952.504395 ...   3.033758  0.949735   \n",
       "124384         98  6907385   1  2952.504395 ...   3.238828  0.955494   \n",
       "\n",
       "              z3        a1        x2        y2     sina1     cosa1     s2  N2  \n",
       "0      -0.032622 -2.981392 -0.032379 -0.003833 -0.159516 -0.987195      0   2  \n",
       "1      -0.050693 -2.948250 -0.050217 -0.006436 -0.192141 -0.981367      1   2  \n",
       "2      -0.046716 -3.001281 -0.046517 -0.003711 -0.139852 -0.990172      2   1  \n",
       "3      -0.023579 -3.034546 -0.023504 -0.001800 -0.106842 -0.994276      3   1  \n",
       "4      -0.066483 -2.942832 -0.065919 -0.007436 -0.197455 -0.980312      4   1  \n",
       "5      -0.045146 -2.992965 -0.044918 -0.004052 -0.148081 -0.988975      5   1  \n",
       "6      -0.042480  3.143731 -0.042382  0.002251 -0.002139 -0.999998      6   2  \n",
       "7      -0.040435 -2.879700 -0.039520 -0.008397 -0.258909 -0.965902      7   2  \n",
       "8      -0.052037 -3.050773 -0.051952 -0.001205 -0.090695 -0.995879      8   2  \n",
       "9      -0.055475  3.210259 -0.055390  0.000190 -0.068612 -0.997643      9   2  \n",
       "10     -0.061558 -2.903757 -0.060678 -0.009657 -0.235600 -0.971850     10   2  \n",
       "11     -0.044992 -2.923513 -0.044375 -0.007144 -0.216355 -0.976315     11   2  \n",
       "12     -0.040131 -3.026254 -0.040019 -0.002533 -0.115083 -0.993356     12   1  \n",
       "13     -0.064200 -2.975710 -0.063850 -0.005275 -0.165123 -0.986273     13   1  \n",
       "14     -0.063316 -2.893654 -0.062325 -0.010419 -0.245406 -0.969420     14   1  \n",
       "15     -0.058329 -3.051466 -0.058224 -0.000833 -0.090005 -0.995941     15   2  \n",
       "16     -0.055796 -2.910231 -0.055008 -0.008812 -0.229303 -0.973355     16   6  \n",
       "17     -0.056582 -2.944505 -0.056061 -0.006962 -0.195814 -0.980641     17   1  \n",
       "18     -0.025311 -2.864422 -0.024552 -0.006120 -0.273636 -0.961833     18   1  \n",
       "19     -0.042597 -2.896937 -0.041798 -0.008009 -0.242222 -0.970221     19   2  \n",
       "20     -0.034562 -3.089627 -0.034540 -0.000244 -0.051943 -0.998650     20   1  \n",
       "21     -0.049833 -3.015324 -0.049677 -0.003060 -0.125934 -0.992039     21   1  \n",
       "22     -0.069690 -2.928079 -0.068996 -0.008526 -0.211895 -0.977292     22   2  \n",
       "23     -0.103249 -2.981902 -0.102670 -0.002609 -0.159013 -0.987276     23   2  \n",
       "24     -0.082067 -2.966190 -0.081599 -0.005616 -0.174505 -0.984656     24   1  \n",
       "25     -0.109419 -2.935565 -0.108550 -0.006923 -0.204573 -0.978851     25   1  \n",
       "26     -0.094183  3.257399 -0.093766  0.000625 -0.115548 -0.993302     26   2  \n",
       "27     -0.104190 -2.975054 -0.103579 -0.003214 -0.165770 -0.986164     27   9  \n",
       "28     -0.100574 -2.945559 -0.099856 -0.006523 -0.194781 -0.980847     28   1  \n",
       "29     -0.081007 -2.921360 -0.080210 -0.009259 -0.218456 -0.975847     29   1  \n",
       "...          ...       ...       ...       ...       ...       ...    ...  ..  \n",
       "124355  0.340976  3.778244 -0.313006  0.078626 -0.594506 -0.804091  90959   2  \n",
       "124356  0.316191  3.709849 -0.292148  0.074429 -0.538163 -0.842841  90960   2  \n",
       "124357  0.295030  3.524101 -0.262809  0.104902 -0.373249 -0.927731  90942   3  \n",
       "124358  0.270353  3.629106 -0.255300  0.054168 -0.468430 -0.883501  90961   1  \n",
       "124359  0.270868  3.706116 -0.259123  0.034775 -0.535013 -0.844844  90962   1  \n",
       "124360  0.270333  3.638314 -0.255773  0.051798 -0.476546 -0.879150  90963   1  \n",
       "124361  0.292223  3.628619 -0.270581  0.073903 -0.468000 -0.883728  90946   2  \n",
       "124362  0.259856  3.585336 -0.245136  0.056233 -0.429324 -0.903151  53912   2  \n",
       "124363  0.284891  3.623775 -0.265313  0.068403 -0.463713 -0.885985  90964   1  \n",
       "124364  0.343511  3.813816 -0.317181  0.070298 -0.622727 -0.782439  90965   1  \n",
       "124365  0.325861  3.821214 -0.305653  0.050682 -0.628498 -0.777811  90966   1  \n",
       "124366  0.315215  3.763165 -0.294969  0.058082 -0.582314 -0.812964  90967   1  \n",
       "124367  0.340414  3.777285 -0.312513  0.078632 -0.593734 -0.804661  90959   2  \n",
       "124368  0.297617  3.719129 -0.279993  0.054521 -0.545961 -0.837810  90968   1  \n",
       "124369  0.308419  3.791294 -0.291478  0.043595 -0.604949 -0.796264  90969   1  \n",
       "124370  0.315843  3.707385 -0.291677  0.075052 -0.536085 -0.844164  90960   2  \n",
       "124371  0.315321  3.742775 -0.293795  0.064185 -0.565618 -0.824668  90970   1  \n",
       "124372  0.326532  3.730296 -0.300208  0.078899 -0.555283 -0.831661  90954   2  \n",
       "124373  0.270617  3.781511 -0.260802  0.014788 -0.597130 -0.802144  90971   1  \n",
       "124374  0.277872  3.829668 -0.267628  0.007330 -0.635051 -0.772470  90972   1  \n",
       "124375  0.295562  3.798213 -0.281881  0.029697 -0.610443 -0.792060  90973   2  \n",
       "124376  0.305964  3.849301 -0.291616  0.023683 -0.650094 -0.759854  90974   1  \n",
       "124377  0.299100  3.886759 -0.286462  0.007350 -0.678094 -0.734975  90975   1  \n",
       "124378  0.267698  3.696844 -0.256285  0.034472 -0.527157 -0.849768  90976   1  \n",
       "124379  0.267195  3.793843 -0.257978  0.009136 -0.606977 -0.794720  90977   1  \n",
       "124380  0.301648  3.895570 -0.288701  0.007403 -0.684544 -0.728972  73690   2  \n",
       "124381  0.305043  3.844278 -0.290696  0.025008 -0.646269 -0.763110  90978   1  \n",
       "124382  0.294644  3.798922 -0.281105  0.029327 -0.611005 -0.791627  90973   2  \n",
       "124383  0.329624  3.937011 -0.312538  0.018000 -0.714156 -0.699986  90979   1  \n",
       "124384  0.308754  3.835241 -0.293401  0.030794 -0.639347 -0.768918  90980   1  \n",
       "\n",
       "[124385 rows x 23 columns]"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "labels"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [],
   "source": [
    "cond = np.where((labels['N2'].values > labels['N1'].values) & (labels['N2'].values < 20))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "cond = np.where((labels['N2'].values > labels['N1'].values) & ( ((labels['N2'].values < 20) & ( (labels['x'] > 275) | (labels['y'] > 625))) | (labels['N2'] < 15)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [
    {
     "ename": "TypeError",
     "evalue": "ufunc 'bitwise_or' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-46-e3bad33eb4b7>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0;34m(\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlabels\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'N2'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m \u001b[0;34m<\u001b[0m \u001b[0;36m20\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m&\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mlabels\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'x'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m \u001b[0;34m>\u001b[0m \u001b[0;36m275\u001b[0m \u001b[0;34m|\u001b[0m \u001b[0mlabels\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'y'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m \u001b[0;34m>\u001b[0m \u001b[0;36m625\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m|\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mlabels\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'N2'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m \u001b[0;34m<\u001b[0m \u001b[0;36m15\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m: ufunc 'bitwise_or' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''"
     ]
    }
   ],
   "source": [
    "( ((labels['N2'].values < 20) & (labels['x'].values > 275 | labels['y'].values > 625)) | (labels['N2'].values < 15))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([False, False, False, ..., False, False, False])"
      ]
     },
     "execution_count": 56,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(labels['x'].values > 275) | (labels['y'].values > 625)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([False, False, False, ..., False, False, False])"
      ]
     },
     "execution_count": 55,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "labels['y'].values > 625"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>hit_id</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>z</th>\n",
       "      <th>volume_id</th>\n",
       "      <th>layer_id</th>\n",
       "      <th>module_id</th>\n",
       "      <th>s1</th>\n",
       "      <th>N1</th>\n",
       "      <th>stepped_z</th>\n",
       "      <th>...</th>\n",
       "      <th>z1</th>\n",
       "      <th>z2</th>\n",
       "      <th>z3</th>\n",
       "      <th>a1</th>\n",
       "      <th>x2</th>\n",
       "      <th>y2</th>\n",
       "      <th>sina1</th>\n",
       "      <th>cosa1</th>\n",
       "      <th>s2</th>\n",
       "      <th>N2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>319</th>\n",
       "      <td>320</td>\n",
       "      <td>-60.398602</td>\n",
       "      <td>-127.692001</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>21</td>\n",
       "      <td>8348561</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.633213</td>\n",
       "      <td>-0.995607</td>\n",
       "      <td>-0.094045</td>\n",
       "      <td>-1.890321</td>\n",
       "      <td>-0.040035</td>\n",
       "      <td>-0.084641</td>\n",
       "      <td>-0.949385</td>\n",
       "      <td>-0.314115</td>\n",
       "      <td>241</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>343</th>\n",
       "      <td>344</td>\n",
       "      <td>-23.916300</td>\n",
       "      <td>-97.369598</td>\n",
       "      <td>-1502.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>23</td>\n",
       "      <td>8348576</td>\n",
       "      <td>2</td>\n",
       "      <td>-1502.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-14.985513</td>\n",
       "      <td>-0.997781</td>\n",
       "      <td>-0.066731</td>\n",
       "      <td>-1.724901</td>\n",
       "      <td>-0.015882</td>\n",
       "      <td>-0.064661</td>\n",
       "      <td>-0.988149</td>\n",
       "      <td>-0.153496</td>\n",
       "      <td>256</td>\n",
       "      <td>13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>363</th>\n",
       "      <td>364</td>\n",
       "      <td>-32.573200</td>\n",
       "      <td>-138.384003</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>24</td>\n",
       "      <td>8348592</td>\n",
       "      <td>2</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.565152</td>\n",
       "      <td>-0.995550</td>\n",
       "      <td>-0.094651</td>\n",
       "      <td>-1.678889</td>\n",
       "      <td>-0.021590</td>\n",
       "      <td>-0.091723</td>\n",
       "      <td>-0.994164</td>\n",
       "      <td>-0.107882</td>\n",
       "      <td>272</td>\n",
       "      <td>13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>370</th>\n",
       "      <td>371</td>\n",
       "      <td>-35.410400</td>\n",
       "      <td>-138.871994</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>24</td>\n",
       "      <td>8348597</td>\n",
       "      <td>2</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.480405</td>\n",
       "      <td>-0.995479</td>\n",
       "      <td>-0.095416</td>\n",
       "      <td>-1.696383</td>\n",
       "      <td>-0.023469</td>\n",
       "      <td>-0.092040</td>\n",
       "      <td>-0.992124</td>\n",
       "      <td>-0.125257</td>\n",
       "      <td>277</td>\n",
       "      <td>7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>391</th>\n",
       "      <td>392</td>\n",
       "      <td>-11.331600</td>\n",
       "      <td>-113.625000</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>26</td>\n",
       "      <td>8348613</td>\n",
       "      <td>6</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-13.118681</td>\n",
       "      <td>-0.997107</td>\n",
       "      <td>-0.076227</td>\n",
       "      <td>-1.571377</td>\n",
       "      <td>-0.007543</td>\n",
       "      <td>-0.075631</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-0.000580</td>\n",
       "      <td>293</td>\n",
       "      <td>16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>393</th>\n",
       "      <td>394</td>\n",
       "      <td>-11.398400</td>\n",
       "      <td>-113.919998</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>27</td>\n",
       "      <td>8348613</td>\n",
       "      <td>6</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-13.119224</td>\n",
       "      <td>-0.997108</td>\n",
       "      <td>-0.076224</td>\n",
       "      <td>-1.571441</td>\n",
       "      <td>-0.007567</td>\n",
       "      <td>-0.075626</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-0.000645</td>\n",
       "      <td>293</td>\n",
       "      <td>16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>492</th>\n",
       "      <td>493</td>\n",
       "      <td>57.278702</td>\n",
       "      <td>-166.098007</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>34</td>\n",
       "      <td>8348688</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-8.548839</td>\n",
       "      <td>-0.993228</td>\n",
       "      <td>-0.116975</td>\n",
       "      <td>-1.086531</td>\n",
       "      <td>0.037877</td>\n",
       "      <td>-0.109835</td>\n",
       "      <td>-0.885017</td>\n",
       "      <td>0.465559</td>\n",
       "      <td>368</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>649</th>\n",
       "      <td>650</td>\n",
       "      <td>134.445007</td>\n",
       "      <td>-66.787003</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>46</td>\n",
       "      <td>8348804</td>\n",
       "      <td>7</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.005369</td>\n",
       "      <td>-0.995042</td>\n",
       "      <td>-0.099946</td>\n",
       "      <td>-0.331070</td>\n",
       "      <td>0.089067</td>\n",
       "      <td>-0.044245</td>\n",
       "      <td>-0.325055</td>\n",
       "      <td>0.945695</td>\n",
       "      <td>484</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>667</th>\n",
       "      <td>668</td>\n",
       "      <td>58.875198</td>\n",
       "      <td>-23.493900</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>47</td>\n",
       "      <td>8348816</td>\n",
       "      <td>1</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-23.631674</td>\n",
       "      <td>-0.999106</td>\n",
       "      <td>-0.042316</td>\n",
       "      <td>-0.324867</td>\n",
       "      <td>0.039267</td>\n",
       "      <td>-0.015669</td>\n",
       "      <td>-0.319183</td>\n",
       "      <td>0.947693</td>\n",
       "      <td>496</td>\n",
       "      <td>11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>688</th>\n",
       "      <td>689</td>\n",
       "      <td>134.128006</td>\n",
       "      <td>-66.574799</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>48</td>\n",
       "      <td>8348804</td>\n",
       "      <td>7</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.003934</td>\n",
       "      <td>-0.995041</td>\n",
       "      <td>-0.099961</td>\n",
       "      <td>-0.331070</td>\n",
       "      <td>0.089094</td>\n",
       "      <td>-0.044222</td>\n",
       "      <td>-0.325056</td>\n",
       "      <td>0.945695</td>\n",
       "      <td>484</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>701</th>\n",
       "      <td>702</td>\n",
       "      <td>141.481995</td>\n",
       "      <td>-52.666901</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>49</td>\n",
       "      <td>8257314</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.949241</td>\n",
       "      <td>-0.994987</td>\n",
       "      <td>-0.100510</td>\n",
       "      <td>-0.225641</td>\n",
       "      <td>0.093723</td>\n",
       "      <td>-0.034889</td>\n",
       "      <td>-0.223731</td>\n",
       "      <td>0.974651</td>\n",
       "      <td>520</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>702</th>\n",
       "      <td>703</td>\n",
       "      <td>148.559006</td>\n",
       "      <td>-42.171799</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>49</td>\n",
       "      <td>8348841</td>\n",
       "      <td>2</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.726197</td>\n",
       "      <td>-0.994756</td>\n",
       "      <td>-0.102815</td>\n",
       "      <td>-0.142873</td>\n",
       "      <td>0.098389</td>\n",
       "      <td>-0.027930</td>\n",
       "      <td>-0.142388</td>\n",
       "      <td>0.989811</td>\n",
       "      <td>521</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>733</th>\n",
       "      <td>734</td>\n",
       "      <td>148.192993</td>\n",
       "      <td>-41.968700</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>51</td>\n",
       "      <td>8348841</td>\n",
       "      <td>2</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.725962</td>\n",
       "      <td>-0.994756</td>\n",
       "      <td>-0.102818</td>\n",
       "      <td>-0.142607</td>\n",
       "      <td>0.098408</td>\n",
       "      <td>-0.027869</td>\n",
       "      <td>-0.142124</td>\n",
       "      <td>0.989849</td>\n",
       "      <td>521</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>771</th>\n",
       "      <td>772</td>\n",
       "      <td>140.257004</td>\n",
       "      <td>-3.155920</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>54</td>\n",
       "      <td>8348896</td>\n",
       "      <td>6</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.677723</td>\n",
       "      <td>-0.995643</td>\n",
       "      <td>-0.093653</td>\n",
       "      <td>0.098959</td>\n",
       "      <td>0.093221</td>\n",
       "      <td>-0.002098</td>\n",
       "      <td>0.098798</td>\n",
       "      <td>0.995108</td>\n",
       "      <td>576</td>\n",
       "      <td>10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1202</th>\n",
       "      <td>1203</td>\n",
       "      <td>-10.810400</td>\n",
       "      <td>93.706299</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>82</td>\n",
       "      <td>8349210</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-15.923242</td>\n",
       "      <td>-0.998034</td>\n",
       "      <td>-0.062801</td>\n",
       "      <td>1.767261</td>\n",
       "      <td>-0.007183</td>\n",
       "      <td>0.062265</td>\n",
       "      <td>0.980763</td>\n",
       "      <td>-0.195203</td>\n",
       "      <td>890</td>\n",
       "      <td>16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1225</th>\n",
       "      <td>1226</td>\n",
       "      <td>-4.662660</td>\n",
       "      <td>172.847000</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>83</td>\n",
       "      <td>8349237</td>\n",
       "      <td>1</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-8.663498</td>\n",
       "      <td>-0.993404</td>\n",
       "      <td>-0.115427</td>\n",
       "      <td>1.747532</td>\n",
       "      <td>-0.003092</td>\n",
       "      <td>0.114624</td>\n",
       "      <td>0.984423</td>\n",
       "      <td>-0.175818</td>\n",
       "      <td>917</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1234</th>\n",
       "      <td>1235</td>\n",
       "      <td>-15.301700</td>\n",
       "      <td>133.399994</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>84</td>\n",
       "      <td>8349243</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-11.186055</td>\n",
       "      <td>-0.996028</td>\n",
       "      <td>-0.089397</td>\n",
       "      <td>1.801239</td>\n",
       "      <td>-0.010147</td>\n",
       "      <td>0.088462</td>\n",
       "      <td>0.973565</td>\n",
       "      <td>-0.228408</td>\n",
       "      <td>923</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1243</th>\n",
       "      <td>1244</td>\n",
       "      <td>-19.023701</td>\n",
       "      <td>63.438000</td>\n",
       "      <td>-1497.5</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>85</td>\n",
       "      <td>8349246</td>\n",
       "      <td>1</td>\n",
       "      <td>-1497.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-22.611006</td>\n",
       "      <td>-0.999023</td>\n",
       "      <td>-0.044226</td>\n",
       "      <td>1.919416</td>\n",
       "      <td>-0.012691</td>\n",
       "      <td>0.042321</td>\n",
       "      <td>0.939845</td>\n",
       "      <td>-0.341601</td>\n",
       "      <td>926</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1299</th>\n",
       "      <td>1300</td>\n",
       "      <td>-44.287102</td>\n",
       "      <td>129.695999</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>88</td>\n",
       "      <td>8349287</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-10.959623</td>\n",
       "      <td>-0.995863</td>\n",
       "      <td>-0.091244</td>\n",
       "      <td>2.018493</td>\n",
       "      <td>-0.029363</td>\n",
       "      <td>0.085991</td>\n",
       "      <td>0.901447</td>\n",
       "      <td>-0.432891</td>\n",
       "      <td>967</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1355</th>\n",
       "      <td>1356</td>\n",
       "      <td>-57.589699</td>\n",
       "      <td>74.781303</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>93</td>\n",
       "      <td>8349329</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-15.913335</td>\n",
       "      <td>-0.998031</td>\n",
       "      <td>-0.062840</td>\n",
       "      <td>2.308702</td>\n",
       "      <td>-0.038266</td>\n",
       "      <td>0.049690</td>\n",
       "      <td>0.739879</td>\n",
       "      <td>-0.672740</td>\n",
       "      <td>1009</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1375</th>\n",
       "      <td>1376</td>\n",
       "      <td>-110.156998</td>\n",
       "      <td>127.799004</td>\n",
       "      <td>-1502.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>94</td>\n",
       "      <td>8349345</td>\n",
       "      <td>1</td>\n",
       "      <td>-1501.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-8.902238</td>\n",
       "      <td>-0.993750</td>\n",
       "      <td>-0.112331</td>\n",
       "      <td>2.428321</td>\n",
       "      <td>-0.072882</td>\n",
       "      <td>0.084554</td>\n",
       "      <td>0.654312</td>\n",
       "      <td>-0.756225</td>\n",
       "      <td>1025</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1517</th>\n",
       "      <td>1518</td>\n",
       "      <td>-129.451004</td>\n",
       "      <td>31.396900</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>105</td>\n",
       "      <td>8349453</td>\n",
       "      <td>1</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-11.245934</td>\n",
       "      <td>-0.996070</td>\n",
       "      <td>-0.088921</td>\n",
       "      <td>3.018956</td>\n",
       "      <td>-0.086076</td>\n",
       "      <td>0.020877</td>\n",
       "      <td>0.122329</td>\n",
       "      <td>-0.992490</td>\n",
       "      <td>1133</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1564</th>\n",
       "      <td>1565</td>\n",
       "      <td>-155.727997</td>\n",
       "      <td>9.427130</td>\n",
       "      <td>-1498.0</td>\n",
       "      <td>7</td>\n",
       "      <td>2</td>\n",
       "      <td>108</td>\n",
       "      <td>8349479</td>\n",
       "      <td>1</td>\n",
       "      <td>-1497.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.601788</td>\n",
       "      <td>-0.994620</td>\n",
       "      <td>-0.104147</td>\n",
       "      <td>3.216228</td>\n",
       "      <td>-0.103398</td>\n",
       "      <td>0.006259</td>\n",
       "      <td>-0.074567</td>\n",
       "      <td>-0.997216</td>\n",
       "      <td>1159</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1571</th>\n",
       "      <td>1572</td>\n",
       "      <td>-79.428703</td>\n",
       "      <td>-10.840300</td>\n",
       "      <td>-1302.5</td>\n",
       "      <td>7</td>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "      <td>8349484</td>\n",
       "      <td>2</td>\n",
       "      <td>-1302.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-16.247789</td>\n",
       "      <td>-0.998111</td>\n",
       "      <td>-0.061547</td>\n",
       "      <td>-2.936612</td>\n",
       "      <td>-0.060866</td>\n",
       "      <td>-0.008307</td>\n",
       "      <td>-0.203548</td>\n",
       "      <td>-0.979065</td>\n",
       "      <td>1164</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1603</th>\n",
       "      <td>1604</td>\n",
       "      <td>-134.376007</td>\n",
       "      <td>-6.137480</td>\n",
       "      <td>-1302.0</td>\n",
       "      <td>7</td>\n",
       "      <td>4</td>\n",
       "      <td>2</td>\n",
       "      <td>8348343</td>\n",
       "      <td>1</td>\n",
       "      <td>-1301.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-9.679172</td>\n",
       "      <td>-0.994705</td>\n",
       "      <td>-0.103315</td>\n",
       "      <td>-2.979505</td>\n",
       "      <td>-0.102661</td>\n",
       "      <td>-0.004689</td>\n",
       "      <td>-0.161379</td>\n",
       "      <td>-0.986892</td>\n",
       "      <td>23</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1617</th>\n",
       "      <td>1618</td>\n",
       "      <td>-79.156403</td>\n",
       "      <td>-10.750700</td>\n",
       "      <td>-1298.0</td>\n",
       "      <td>7</td>\n",
       "      <td>4</td>\n",
       "      <td>4</td>\n",
       "      <td>8349484</td>\n",
       "      <td>2</td>\n",
       "      <td>-1297.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-16.248793</td>\n",
       "      <td>-0.998112</td>\n",
       "      <td>-0.061543</td>\n",
       "      <td>-2.937507</td>\n",
       "      <td>-0.060868</td>\n",
       "      <td>-0.008267</td>\n",
       "      <td>-0.202672</td>\n",
       "      <td>-0.979247</td>\n",
       "      <td>1164</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1802</th>\n",
       "      <td>1803</td>\n",
       "      <td>-121.322998</td>\n",
       "      <td>-100.648003</td>\n",
       "      <td>-1298.0</td>\n",
       "      <td>7</td>\n",
       "      <td>4</td>\n",
       "      <td>13</td>\n",
       "      <td>8349635</td>\n",
       "      <td>4</td>\n",
       "      <td>-1297.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-8.234158</td>\n",
       "      <td>-0.992706</td>\n",
       "      <td>-0.121445</td>\n",
       "      <td>-2.312562</td>\n",
       "      <td>-0.092787</td>\n",
       "      <td>-0.076975</td>\n",
       "      <td>-0.737277</td>\n",
       "      <td>-0.675591</td>\n",
       "      <td>1315</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1858</th>\n",
       "      <td>1859</td>\n",
       "      <td>-102.922997</td>\n",
       "      <td>-136.904999</td>\n",
       "      <td>-1298.0</td>\n",
       "      <td>7</td>\n",
       "      <td>4</td>\n",
       "      <td>16</td>\n",
       "      <td>8349670</td>\n",
       "      <td>2</td>\n",
       "      <td>-1297.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-7.578352</td>\n",
       "      <td>-0.991406</td>\n",
       "      <td>-0.131955</td>\n",
       "      <td>-2.067089</td>\n",
       "      <td>-0.078612</td>\n",
       "      <td>-0.104567</td>\n",
       "      <td>-0.879354</td>\n",
       "      <td>-0.476168</td>\n",
       "      <td>1350</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1906</th>\n",
       "      <td>1907</td>\n",
       "      <td>-103.083000</td>\n",
       "      <td>-137.432999</td>\n",
       "      <td>-1302.0</td>\n",
       "      <td>7</td>\n",
       "      <td>4</td>\n",
       "      <td>18</td>\n",
       "      <td>8349670</td>\n",
       "      <td>2</td>\n",
       "      <td>-1301.995605</td>\n",
       "      <td>...</td>\n",
       "      <td>-7.578773</td>\n",
       "      <td>-0.991407</td>\n",
       "      <td>-0.131947</td>\n",
       "      <td>-2.065536</td>\n",
       "      <td>-0.078492</td>\n",
       "      <td>-0.104648</td>\n",
       "      <td>-0.880092</td>\n",
       "      <td>-0.474803</td>\n",
       "      <td>1350</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000</th>\n",
       "      <td>2001</td>\n",
       "      <td>-20.154900</td>\n",
       "      <td>-84.700798</td>\n",
       "      <td>-1302.5</td>\n",
       "      <td>7</td>\n",
       "      <td>4</td>\n",
       "      <td>23</td>\n",
       "      <td>8348576</td>\n",
       "      <td>2</td>\n",
       "      <td>-1302.495605</td>\n",
       "      <td>...</td>\n",
       "      <td>-14.960009</td>\n",
       "      <td>-0.997773</td>\n",
       "      <td>-0.066845</td>\n",
       "      <td>-1.729089</td>\n",
       "      <td>-0.015440</td>\n",
       "      <td>-0.064884</td>\n",
       "      <td>-0.987498</td>\n",
       "      <td>-0.157633</td>\n",
       "      <td>256</td>\n",
       "      <td>13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123820</th>\n",
       "      <td>123821</td>\n",
       "      <td>-821.000977</td>\n",
       "      <td>286.739014</td>\n",
       "      <td>2544.5</td>\n",
       "      <td>18</td>\n",
       "      <td>10</td>\n",
       "      <td>93</td>\n",
       "      <td>8438885</td>\n",
       "      <td>2</td>\n",
       "      <td>2544.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>2.925941</td>\n",
       "      <td>0.946261</td>\n",
       "      <td>0.341770</td>\n",
       "      <td>3.566593</td>\n",
       "      <td>-0.305318</td>\n",
       "      <td>0.106634</td>\n",
       "      <td>-0.412321</td>\n",
       "      <td>-0.911039</td>\n",
       "      <td>90565</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123827</th>\n",
       "      <td>123828</td>\n",
       "      <td>-828.724976</td>\n",
       "      <td>295.891998</td>\n",
       "      <td>2552.5</td>\n",
       "      <td>18</td>\n",
       "      <td>10</td>\n",
       "      <td>94</td>\n",
       "      <td>8438885</td>\n",
       "      <td>2</td>\n",
       "      <td>2552.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>2.900681</td>\n",
       "      <td>0.945397</td>\n",
       "      <td>0.344747</td>\n",
       "      <td>3.568825</td>\n",
       "      <td>-0.306944</td>\n",
       "      <td>0.109593</td>\n",
       "      <td>-0.414353</td>\n",
       "      <td>-0.910116</td>\n",
       "      <td>90565</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123867</th>\n",
       "      <td>123868</td>\n",
       "      <td>-880.921021</td>\n",
       "      <td>-187.223999</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>3</td>\n",
       "      <td>8421026</td>\n",
       "      <td>2</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.278377</td>\n",
       "      <td>0.956492</td>\n",
       "      <td>0.305029</td>\n",
       "      <td>-2.143711</td>\n",
       "      <td>-0.285384</td>\n",
       "      <td>-0.060653</td>\n",
       "      <td>-0.840325</td>\n",
       "      <td>-0.542084</td>\n",
       "      <td>72706</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123873</th>\n",
       "      <td>123874</td>\n",
       "      <td>-881.361023</td>\n",
       "      <td>-186.447998</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>4</td>\n",
       "      <td>8421026</td>\n",
       "      <td>2</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.268516</td>\n",
       "      <td>0.956246</td>\n",
       "      <td>0.305949</td>\n",
       "      <td>-2.144416</td>\n",
       "      <td>-0.286228</td>\n",
       "      <td>-0.060550</td>\n",
       "      <td>-0.839942</td>\n",
       "      <td>-0.542676</td>\n",
       "      <td>72706</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123890</th>\n",
       "      <td>123891</td>\n",
       "      <td>-729.818970</td>\n",
       "      <td>-435.549011</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>10</td>\n",
       "      <td>8422100</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.468029</td>\n",
       "      <td>0.960853</td>\n",
       "      <td>0.288348</td>\n",
       "      <td>-1.860004</td>\n",
       "      <td>-0.237913</td>\n",
       "      <td>-0.141984</td>\n",
       "      <td>-0.958470</td>\n",
       "      <td>-0.285193</td>\n",
       "      <td>73780</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123912</th>\n",
       "      <td>123913</td>\n",
       "      <td>-471.006989</td>\n",
       "      <td>-591.406982</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>14</td>\n",
       "      <td>8423147</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.898552</td>\n",
       "      <td>0.968642</td>\n",
       "      <td>0.256505</td>\n",
       "      <td>-1.582836</td>\n",
       "      <td>-0.154788</td>\n",
       "      <td>-0.194355</td>\n",
       "      <td>-0.999928</td>\n",
       "      <td>-0.012040</td>\n",
       "      <td>74827</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123929</th>\n",
       "      <td>123930</td>\n",
       "      <td>-421.453003</td>\n",
       "      <td>-639.172974</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>16</td>\n",
       "      <td>8420050</td>\n",
       "      <td>1</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.845927</td>\n",
       "      <td>0.967819</td>\n",
       "      <td>0.260015</td>\n",
       "      <td>-1.484769</td>\n",
       "      <td>-0.138526</td>\n",
       "      <td>-0.210088</td>\n",
       "      <td>-0.996302</td>\n",
       "      <td>0.085921</td>\n",
       "      <td>71730</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123936</th>\n",
       "      <td>123937</td>\n",
       "      <td>-412.959015</td>\n",
       "      <td>-847.940002</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>17</td>\n",
       "      <td>8438567</td>\n",
       "      <td>1</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.133634</td>\n",
       "      <td>0.952668</td>\n",
       "      <td>0.319118</td>\n",
       "      <td>-1.197764</td>\n",
       "      <td>-0.133112</td>\n",
       "      <td>-0.273323</td>\n",
       "      <td>-0.931226</td>\n",
       "      <td>0.364441</td>\n",
       "      <td>90247</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123979</th>\n",
       "      <td>123980</td>\n",
       "      <td>94.224503</td>\n",
       "      <td>-802.755981</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>26</td>\n",
       "      <td>8419050</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.646686</td>\n",
       "      <td>0.964397</td>\n",
       "      <td>0.274222</td>\n",
       "      <td>-0.747281</td>\n",
       "      <td>0.030830</td>\n",
       "      <td>-0.262655</td>\n",
       "      <td>-0.679647</td>\n",
       "      <td>0.733540</td>\n",
       "      <td>70730</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124014</th>\n",
       "      <td>124015</td>\n",
       "      <td>425.089996</td>\n",
       "      <td>-711.791016</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>34</td>\n",
       "      <td>8420272</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.555207</td>\n",
       "      <td>0.962644</td>\n",
       "      <td>0.281278</td>\n",
       "      <td>-0.307351</td>\n",
       "      <td>0.138833</td>\n",
       "      <td>-0.232469</td>\n",
       "      <td>-0.302535</td>\n",
       "      <td>0.953138</td>\n",
       "      <td>71952</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124037</th>\n",
       "      <td>124038</td>\n",
       "      <td>656.642029</td>\n",
       "      <td>-488.907013</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>40</td>\n",
       "      <td>8402142</td>\n",
       "      <td>1</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.596713</td>\n",
       "      <td>0.963455</td>\n",
       "      <td>0.278032</td>\n",
       "      <td>0.075864</td>\n",
       "      <td>0.214857</td>\n",
       "      <td>-0.159973</td>\n",
       "      <td>0.075792</td>\n",
       "      <td>0.997124</td>\n",
       "      <td>53822</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124045</th>\n",
       "      <td>124046</td>\n",
       "      <td>669.901001</td>\n",
       "      <td>-372.571991</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>42</td>\n",
       "      <td>8422440</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.845215</td>\n",
       "      <td>0.967808</td>\n",
       "      <td>0.260063</td>\n",
       "      <td>0.162218</td>\n",
       "      <td>0.219961</td>\n",
       "      <td>-0.122334</td>\n",
       "      <td>0.161507</td>\n",
       "      <td>0.986872</td>\n",
       "      <td>74120</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124086</th>\n",
       "      <td>124087</td>\n",
       "      <td>765.083008</td>\n",
       "      <td>-32.949402</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>48</td>\n",
       "      <td>8419264</td>\n",
       "      <td>1</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.845032</td>\n",
       "      <td>0.967805</td>\n",
       "      <td>0.260076</td>\n",
       "      <td>0.626081</td>\n",
       "      <td>0.251470</td>\n",
       "      <td>-0.010830</td>\n",
       "      <td>0.585973</td>\n",
       "      <td>0.810330</td>\n",
       "      <td>70944</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124092</th>\n",
       "      <td>124093</td>\n",
       "      <td>760.989990</td>\n",
       "      <td>95.395401</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>50</td>\n",
       "      <td>8376287</td>\n",
       "      <td>2</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.843160</td>\n",
       "      <td>0.967775</td>\n",
       "      <td>0.260203</td>\n",
       "      <td>0.794847</td>\n",
       "      <td>0.249862</td>\n",
       "      <td>0.031322</td>\n",
       "      <td>0.713756</td>\n",
       "      <td>0.700394</td>\n",
       "      <td>27967</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124101</th>\n",
       "      <td>124102</td>\n",
       "      <td>761.505981</td>\n",
       "      <td>95.986801</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>52</td>\n",
       "      <td>8376287</td>\n",
       "      <td>2</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.836318</td>\n",
       "      <td>0.967665</td>\n",
       "      <td>0.260667</td>\n",
       "      <td>0.796045</td>\n",
       "      <td>0.250258</td>\n",
       "      <td>0.031545</td>\n",
       "      <td>0.714595</td>\n",
       "      <td>0.699538</td>\n",
       "      <td>27967</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124140</th>\n",
       "      <td>124141</td>\n",
       "      <td>700.924011</td>\n",
       "      <td>597.921021</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>60</td>\n",
       "      <td>8420560</td>\n",
       "      <td>2</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.207943</td>\n",
       "      <td>0.954690</td>\n",
       "      <td>0.311726</td>\n",
       "      <td>1.513101</td>\n",
       "      <td>0.226414</td>\n",
       "      <td>0.193142</td>\n",
       "      <td>0.998336</td>\n",
       "      <td>0.057663</td>\n",
       "      <td>72240</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124211</th>\n",
       "      <td>124212</td>\n",
       "      <td>150.177002</td>\n",
       "      <td>814.348999</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>71</td>\n",
       "      <td>8439182</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.559431</td>\n",
       "      <td>0.962728</td>\n",
       "      <td>0.280944</td>\n",
       "      <td>2.112639</td>\n",
       "      <td>0.049052</td>\n",
       "      <td>0.265987</td>\n",
       "      <td>0.856760</td>\n",
       "      <td>-0.515716</td>\n",
       "      <td>90862</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124217</th>\n",
       "      <td>124218</td>\n",
       "      <td>82.925400</td>\n",
       "      <td>782.651978</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>73</td>\n",
       "      <td>8423709</td>\n",
       "      <td>1</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.741261</td>\n",
       "      <td>0.966085</td>\n",
       "      <td>0.267290</td>\n",
       "      <td>2.153130</td>\n",
       "      <td>0.027208</td>\n",
       "      <td>0.256787</td>\n",
       "      <td>0.835181</td>\n",
       "      <td>-0.549975</td>\n",
       "      <td>75389</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124230</th>\n",
       "      <td>124231</td>\n",
       "      <td>-76.138397</td>\n",
       "      <td>772.674988</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>75</td>\n",
       "      <td>8382488</td>\n",
       "      <td>3</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.796278</td>\n",
       "      <td>0.967013</td>\n",
       "      <td>0.263416</td>\n",
       "      <td>2.347528</td>\n",
       "      <td>-0.024979</td>\n",
       "      <td>0.253499</td>\n",
       "      <td>0.713208</td>\n",
       "      <td>-0.700952</td>\n",
       "      <td>34168</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124232</th>\n",
       "      <td>124233</td>\n",
       "      <td>-71.489998</td>\n",
       "      <td>783.401978</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>75</td>\n",
       "      <td>8382488</td>\n",
       "      <td>3</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.746862</td>\n",
       "      <td>0.966181</td>\n",
       "      <td>0.266890</td>\n",
       "      <td>2.349362</td>\n",
       "      <td>-0.023434</td>\n",
       "      <td>0.256797</td>\n",
       "      <td>0.711921</td>\n",
       "      <td>-0.702259</td>\n",
       "      <td>34168</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124290</th>\n",
       "      <td>124291</td>\n",
       "      <td>-415.937988</td>\n",
       "      <td>680.382019</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>83</td>\n",
       "      <td>8402206</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.696158</td>\n",
       "      <td>0.965295</td>\n",
       "      <td>0.270551</td>\n",
       "      <td>2.816610</td>\n",
       "      <td>-0.136218</td>\n",
       "      <td>0.222823</td>\n",
       "      <td>0.319293</td>\n",
       "      <td>-0.947656</td>\n",
       "      <td>53886</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124294</th>\n",
       "      <td>124295</td>\n",
       "      <td>-615.476990</td>\n",
       "      <td>755.770996</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>84</td>\n",
       "      <td>8421896</td>\n",
       "      <td>2</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.032272</td>\n",
       "      <td>0.949689</td>\n",
       "      <td>0.329786</td>\n",
       "      <td>3.108490</td>\n",
       "      <td>-0.197771</td>\n",
       "      <td>0.242852</td>\n",
       "      <td>0.033097</td>\n",
       "      <td>-0.999452</td>\n",
       "      <td>73576</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124301</th>\n",
       "      <td>124302</td>\n",
       "      <td>-506.084991</td>\n",
       "      <td>758.119019</td>\n",
       "      <td>2955.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>84</td>\n",
       "      <td>8439238</td>\n",
       "      <td>1</td>\n",
       "      <td>2955.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.242386</td>\n",
       "      <td>0.955585</td>\n",
       "      <td>0.308415</td>\n",
       "      <td>2.957568</td>\n",
       "      <td>-0.163630</td>\n",
       "      <td>0.245119</td>\n",
       "      <td>0.182988</td>\n",
       "      <td>-0.983115</td>\n",
       "      <td>90918</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124305</th>\n",
       "      <td>124306</td>\n",
       "      <td>-472.179993</td>\n",
       "      <td>604.210999</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>85</td>\n",
       "      <td>8420846</td>\n",
       "      <td>4</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.839839</td>\n",
       "      <td>0.967722</td>\n",
       "      <td>0.260428</td>\n",
       "      <td>2.904179</td>\n",
       "      <td>-0.155184</td>\n",
       "      <td>0.198577</td>\n",
       "      <td>0.235190</td>\n",
       "      <td>-0.971950</td>\n",
       "      <td>72526</td>\n",
       "      <td>7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124306</th>\n",
       "      <td>124307</td>\n",
       "      <td>-517.020996</td>\n",
       "      <td>606.385010</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>85</td>\n",
       "      <td>8363872</td>\n",
       "      <td>2</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.695043</td>\n",
       "      <td>0.965275</td>\n",
       "      <td>0.270633</td>\n",
       "      <td>2.973414</td>\n",
       "      <td>-0.169492</td>\n",
       "      <td>0.198787</td>\n",
       "      <td>0.167387</td>\n",
       "      <td>-0.985891</td>\n",
       "      <td>15552</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124307</th>\n",
       "      <td>124308</td>\n",
       "      <td>-615.684021</td>\n",
       "      <td>756.870972</td>\n",
       "      <td>2952.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>86</td>\n",
       "      <td>8421896</td>\n",
       "      <td>2</td>\n",
       "      <td>2952.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.026139</td>\n",
       "      <td>0.949500</td>\n",
       "      <td>0.330454</td>\n",
       "      <td>3.108817</td>\n",
       "      <td>-0.197999</td>\n",
       "      <td>0.243404</td>\n",
       "      <td>0.032770</td>\n",
       "      <td>-0.999463</td>\n",
       "      <td>73576</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124316</th>\n",
       "      <td>124317</td>\n",
       "      <td>-639.760010</td>\n",
       "      <td>494.737000</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>87</td>\n",
       "      <td>8439251</td>\n",
       "      <td>2</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.644561</td>\n",
       "      <td>0.964358</td>\n",
       "      <td>0.274381</td>\n",
       "      <td>3.190426</td>\n",
       "      <td>-0.209316</td>\n",
       "      <td>0.161867</td>\n",
       "      <td>-0.048814</td>\n",
       "      <td>-0.998808</td>\n",
       "      <td>90931</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124325</th>\n",
       "      <td>124326</td>\n",
       "      <td>-640.460999</td>\n",
       "      <td>494.403015</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>89</td>\n",
       "      <td>8439251</td>\n",
       "      <td>2</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.639274</td>\n",
       "      <td>0.964260</td>\n",
       "      <td>0.274780</td>\n",
       "      <td>3.191593</td>\n",
       "      <td>-0.209737</td>\n",
       "      <td>0.161906</td>\n",
       "      <td>-0.049979</td>\n",
       "      <td>-0.998750</td>\n",
       "      <td>90931</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124345</th>\n",
       "      <td>124346</td>\n",
       "      <td>-776.440002</td>\n",
       "      <td>288.593994</td>\n",
       "      <td>2944.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>93</td>\n",
       "      <td>8422986</td>\n",
       "      <td>1</td>\n",
       "      <td>2944.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.554698</td>\n",
       "      <td>0.962634</td>\n",
       "      <td>0.281318</td>\n",
       "      <td>3.510165</td>\n",
       "      <td>-0.253839</td>\n",
       "      <td>0.094349</td>\n",
       "      <td>-0.360284</td>\n",
       "      <td>-0.932843</td>\n",
       "      <td>74666</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124362</th>\n",
       "      <td>124363</td>\n",
       "      <td>-746.534973</td>\n",
       "      <td>171.251999</td>\n",
       "      <td>2947.5</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>95</td>\n",
       "      <td>8402232</td>\n",
       "      <td>1</td>\n",
       "      <td>2947.504395</td>\n",
       "      <td>...</td>\n",
       "      <td>3.848280</td>\n",
       "      <td>0.967856</td>\n",
       "      <td>0.259856</td>\n",
       "      <td>3.585336</td>\n",
       "      <td>-0.245136</td>\n",
       "      <td>0.056233</td>\n",
       "      <td>-0.429324</td>\n",
       "      <td>-0.903151</td>\n",
       "      <td>53912</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>4280 rows × 23 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        hit_id           x           y       z  volume_id  layer_id  \\\n",
       "319        320  -60.398602 -127.692001 -1502.0          7         2   \n",
       "343        344  -23.916300  -97.369598 -1502.5          7         2   \n",
       "363        364  -32.573200 -138.384003 -1502.0          7         2   \n",
       "370        371  -35.410400 -138.871994 -1502.0          7         2   \n",
       "391        392  -11.331600 -113.625000 -1498.0          7         2   \n",
       "393        394  -11.398400 -113.919998 -1502.0          7         2   \n",
       "492        493   57.278702 -166.098007 -1502.0          7         2   \n",
       "649        650  134.445007  -66.787003 -1502.0          7         2   \n",
       "667        668   58.875198  -23.493900 -1498.0          7         2   \n",
       "688        689  134.128006  -66.574799 -1498.0          7         2   \n",
       "701        702  141.481995  -52.666901 -1502.0          7         2   \n",
       "702        703  148.559006  -42.171799 -1502.0          7         2   \n",
       "733        734  148.192993  -41.968700 -1498.0          7         2   \n",
       "771        772  140.257004   -3.155920 -1498.0          7         2   \n",
       "1202      1203  -10.810400   93.706299 -1502.0          7         2   \n",
       "1225      1226   -4.662660  172.847000 -1498.0          7         2   \n",
       "1234      1235  -15.301700  133.399994 -1502.0          7         2   \n",
       "1243      1244  -19.023701   63.438000 -1497.5          7         2   \n",
       "1299      1300  -44.287102  129.695999 -1502.0          7         2   \n",
       "1355      1356  -57.589699   74.781303 -1502.0          7         2   \n",
       "1375      1376 -110.156998  127.799004 -1502.0          7         2   \n",
       "1517      1518 -129.451004   31.396900 -1498.0          7         2   \n",
       "1564      1565 -155.727997    9.427130 -1498.0          7         2   \n",
       "1571      1572  -79.428703  -10.840300 -1302.5          7         4   \n",
       "1603      1604 -134.376007   -6.137480 -1302.0          7         4   \n",
       "1617      1618  -79.156403  -10.750700 -1298.0          7         4   \n",
       "1802      1803 -121.322998 -100.648003 -1298.0          7         4   \n",
       "1858      1859 -102.922997 -136.904999 -1298.0          7         4   \n",
       "1906      1907 -103.083000 -137.432999 -1302.0          7         4   \n",
       "2000      2001  -20.154900  -84.700798 -1302.5          7         4   \n",
       "...        ...         ...         ...     ...        ...       ...   \n",
       "123820  123821 -821.000977  286.739014  2544.5         18        10   \n",
       "123827  123828 -828.724976  295.891998  2552.5         18        10   \n",
       "123867  123868 -880.921021 -187.223999  2952.5         18        12   \n",
       "123873  123874 -881.361023 -186.447998  2944.5         18        12   \n",
       "123890  123891 -729.818970 -435.549011  2947.5         18        12   \n",
       "123912  123913 -471.006989 -591.406982  2947.5         18        12   \n",
       "123929  123930 -421.453003 -639.172974  2944.5         18        12   \n",
       "123936  123937 -412.959015 -847.940002  2955.5         18        12   \n",
       "123979  123980   94.224503 -802.755981  2947.5         18        12   \n",
       "124014  124015  425.089996 -711.791016  2947.5         18        12   \n",
       "124037  124038  656.642029 -488.907013  2944.5         18        12   \n",
       "124045  124046  669.901001 -372.571991  2947.5         18        12   \n",
       "124086  124087  765.083008  -32.949402  2944.5         18        12   \n",
       "124092  124093  760.989990   95.395401  2947.5         18        12   \n",
       "124101  124102  761.505981   95.986801  2944.5         18        12   \n",
       "124140  124141  700.924011  597.921021  2955.5         18        12   \n",
       "124211  124212  150.177002  814.348999  2947.5         18        12   \n",
       "124217  124218   82.925400  782.651978  2944.5         18        12   \n",
       "124230  124231  -76.138397  772.674988  2947.5         18        12   \n",
       "124232  124233  -71.489998  783.401978  2947.5         18        12   \n",
       "124290  124291 -415.937988  680.382019  2947.5         18        12   \n",
       "124294  124295 -615.476990  755.770996  2955.5         18        12   \n",
       "124301  124302 -506.084991  758.119019  2955.5         18        12   \n",
       "124305  124306 -472.179993  604.210999  2944.5         18        12   \n",
       "124306  124307 -517.020996  606.385010  2944.5         18        12   \n",
       "124307  124308 -615.684021  756.870972  2952.5         18        12   \n",
       "124316  124317 -639.760010  494.737000  2947.5         18        12   \n",
       "124325  124326 -640.460999  494.403015  2944.5         18        12   \n",
       "124345  124346 -776.440002  288.593994  2944.5         18        12   \n",
       "124362  124363 -746.534973  171.251999  2947.5         18        12   \n",
       "\n",
       "        module_id       s1  N1    stepped_z ...         z1        z2  \\\n",
       "319            21  8348561   1 -1501.995605 ... -10.633213 -0.995607   \n",
       "343            23  8348576   2 -1502.495605 ... -14.985513 -0.997781   \n",
       "363            24  8348592   2 -1501.995605 ... -10.565152 -0.995550   \n",
       "370            24  8348597   2 -1501.995605 ... -10.480405 -0.995479   \n",
       "391            26  8348613   6 -1497.995605 ... -13.118681 -0.997107   \n",
       "393            27  8348613   6 -1501.995605 ... -13.119224 -0.997108   \n",
       "492            34  8348688   1 -1501.995605 ...  -8.548839 -0.993228   \n",
       "649            46  8348804   7 -1501.995605 ... -10.005369 -0.995042   \n",
       "667            47  8348816   1 -1497.995605 ... -23.631674 -0.999106   \n",
       "688            48  8348804   7 -1497.995605 ... -10.003934 -0.995041   \n",
       "701            49  8257314   1 -1501.995605 ...  -9.949241 -0.994987   \n",
       "702            49  8348841   2 -1501.995605 ...  -9.726197 -0.994756   \n",
       "733            51  8348841   2 -1497.995605 ...  -9.725962 -0.994756   \n",
       "771            54  8348896   6 -1497.995605 ... -10.677723 -0.995643   \n",
       "1202           82  8349210   1 -1501.995605 ... -15.923242 -0.998034   \n",
       "1225           83  8349237   1 -1497.995605 ...  -8.663498 -0.993404   \n",
       "1234           84  8349243   1 -1501.995605 ... -11.186055 -0.996028   \n",
       "1243           85  8349246   1 -1497.495605 ... -22.611006 -0.999023   \n",
       "1299           88  8349287   1 -1501.995605 ... -10.959623 -0.995863   \n",
       "1355           93  8349329   1 -1501.995605 ... -15.913335 -0.998031   \n",
       "1375           94  8349345   1 -1501.995605 ...  -8.902238 -0.993750   \n",
       "1517          105  8349453   1 -1497.995605 ... -11.245934 -0.996070   \n",
       "1564          108  8349479   1 -1497.995605 ...  -9.601788 -0.994620   \n",
       "1571            1  8349484   2 -1302.495605 ... -16.247789 -0.998111   \n",
       "1603            2  8348343   1 -1301.995605 ...  -9.679172 -0.994705   \n",
       "1617            4  8349484   2 -1297.995605 ... -16.248793 -0.998112   \n",
       "1802           13  8349635   4 -1297.995605 ...  -8.234158 -0.992706   \n",
       "1858           16  8349670   2 -1297.995605 ...  -7.578352 -0.991406   \n",
       "1906           18  8349670   2 -1301.995605 ...  -7.578773 -0.991407   \n",
       "2000           23  8348576   2 -1302.495605 ... -14.960009 -0.997773   \n",
       "...           ...      ...  ..          ... ...        ...       ...   \n",
       "123820         93  8438885   2  2544.504395 ...   2.925941  0.946261   \n",
       "123827         94  8438885   2  2552.504395 ...   2.900681  0.945397   \n",
       "123867          3  8421026   2  2952.504395 ...   3.278377  0.956492   \n",
       "123873          4  8421026   2  2944.504395 ...   3.268516  0.956246   \n",
       "123890         10  8422100   1  2947.504395 ...   3.468029  0.960853   \n",
       "123912         14  8423147   1  2947.504395 ...   3.898552  0.968642   \n",
       "123929         16  8420050   1  2944.504395 ...   3.845927  0.967819   \n",
       "123936         17  8438567   1  2955.504395 ...   3.133634  0.952668   \n",
       "123979         26  8419050   1  2947.504395 ...   3.646686  0.964397   \n",
       "124014         34  8420272   1  2947.504395 ...   3.555207  0.962644   \n",
       "124037         40  8402142   1  2944.504395 ...   3.596713  0.963455   \n",
       "124045         42  8422440   1  2947.504395 ...   3.845215  0.967808   \n",
       "124086         48  8419264   1  2944.504395 ...   3.845032  0.967805   \n",
       "124092         50  8376287   2  2947.504395 ...   3.843160  0.967775   \n",
       "124101         52  8376287   2  2944.504395 ...   3.836318  0.967665   \n",
       "124140         60  8420560   2  2955.504395 ...   3.207943  0.954690   \n",
       "124211         71  8439182   1  2947.504395 ...   3.559431  0.962728   \n",
       "124217         73  8423709   1  2944.504395 ...   3.741261  0.966085   \n",
       "124230         75  8382488   3  2947.504395 ...   3.796278  0.967013   \n",
       "124232         75  8382488   3  2947.504395 ...   3.746862  0.966181   \n",
       "124290         83  8402206   1  2947.504395 ...   3.696158  0.965295   \n",
       "124294         84  8421896   2  2955.504395 ...   3.032272  0.949689   \n",
       "124301         84  8439238   1  2955.504395 ...   3.242386  0.955585   \n",
       "124305         85  8420846   4  2944.504395 ...   3.839839  0.967722   \n",
       "124306         85  8363872   2  2944.504395 ...   3.695043  0.965275   \n",
       "124307         86  8421896   2  2952.504395 ...   3.026139  0.949500   \n",
       "124316         87  8439251   2  2947.504395 ...   3.644561  0.964358   \n",
       "124325         89  8439251   2  2944.504395 ...   3.639274  0.964260   \n",
       "124345         93  8422986   1  2944.504395 ...   3.554698  0.962634   \n",
       "124362         95  8402232   1  2947.504395 ...   3.848280  0.967856   \n",
       "\n",
       "              z3        a1        x2        y2     sina1     cosa1     s2  N2  \n",
       "319    -0.094045 -1.890321 -0.040035 -0.084641 -0.949385 -0.314115    241   2  \n",
       "343    -0.066731 -1.724901 -0.015882 -0.064661 -0.988149 -0.153496    256  13  \n",
       "363    -0.094651 -1.678889 -0.021590 -0.091723 -0.994164 -0.107882    272  13  \n",
       "370    -0.095416 -1.696383 -0.023469 -0.092040 -0.992124 -0.125257    277   7  \n",
       "391    -0.076227 -1.571377 -0.007543 -0.075631 -1.000000 -0.000580    293  16  \n",
       "393    -0.076224 -1.571441 -0.007567 -0.075626 -1.000000 -0.000645    293  16  \n",
       "492    -0.116975 -1.086531  0.037877 -0.109835 -0.885017  0.465559    368   3  \n",
       "649    -0.099946 -0.331070  0.089067 -0.044245 -0.325055  0.945695    484   8  \n",
       "667    -0.042316 -0.324867  0.039267 -0.015669 -0.319183  0.947693    496  11  \n",
       "688    -0.099961 -0.331070  0.089094 -0.044222 -0.325056  0.945695    484   8  \n",
       "701    -0.100510 -0.225641  0.093723 -0.034889 -0.223731  0.974651    520   2  \n",
       "702    -0.102815 -0.142873  0.098389 -0.027930 -0.142388  0.989811    521   3  \n",
       "733    -0.102818 -0.142607  0.098408 -0.027869 -0.142124  0.989849    521   3  \n",
       "771    -0.093653  0.098959  0.093221 -0.002098  0.098798  0.995108    576  10  \n",
       "1202   -0.062801  1.767261 -0.007183  0.062265  0.980763 -0.195203    890  16  \n",
       "1225   -0.115427  1.747532 -0.003092  0.114624  0.984423 -0.175818    917   3  \n",
       "1234   -0.089397  1.801239 -0.010147  0.088462  0.973565 -0.228408    923   3  \n",
       "1243   -0.044226  1.919416 -0.012691  0.042321  0.939845 -0.341601    926   3  \n",
       "1299   -0.091244  2.018493 -0.029363  0.085991  0.901447 -0.432891    967   2  \n",
       "1355   -0.062840  2.308702 -0.038266  0.049690  0.739879 -0.672740   1009   3  \n",
       "1375   -0.112331  2.428321 -0.072882  0.084554  0.654312 -0.756225   1025   3  \n",
       "1517   -0.088921  3.018956 -0.086076  0.020877  0.122329 -0.992490   1133   2  \n",
       "1564   -0.104147  3.216228 -0.103398  0.006259 -0.074567 -0.997216   1159   3  \n",
       "1571   -0.061547 -2.936612 -0.060866 -0.008307 -0.203548 -0.979065   1164   4  \n",
       "1603   -0.103315 -2.979505 -0.102661 -0.004689 -0.161379 -0.986892     23   2  \n",
       "1617   -0.061543 -2.937507 -0.060868 -0.008267 -0.202672 -0.979247   1164   4  \n",
       "1802   -0.121445 -2.312562 -0.092787 -0.076975 -0.737277 -0.675591   1315   6  \n",
       "1858   -0.131955 -2.067089 -0.078612 -0.104567 -0.879354 -0.476168   1350   3  \n",
       "1906   -0.131947 -2.065536 -0.078492 -0.104648 -0.880092 -0.474803   1350   3  \n",
       "2000   -0.066845 -1.729089 -0.015440 -0.064884 -0.987498 -0.157633    256  13  \n",
       "...          ...       ...       ...       ...       ...       ...    ...  ..  \n",
       "123820  0.341770  3.566593 -0.305318  0.106634 -0.412321 -0.911039  90565   3  \n",
       "123827  0.344747  3.568825 -0.306944  0.109593 -0.414353 -0.910116  90565   3  \n",
       "123867  0.305029 -2.143711 -0.285384 -0.060653 -0.840325 -0.542084  72706   5  \n",
       "123873  0.305949 -2.144416 -0.286228 -0.060550 -0.839942 -0.542676  72706   5  \n",
       "123890  0.288348 -1.860004 -0.237913 -0.141984 -0.958470 -0.285193  73780   2  \n",
       "123912  0.256505 -1.582836 -0.154788 -0.194355 -0.999928 -0.012040  74827   2  \n",
       "123929  0.260015 -1.484769 -0.138526 -0.210088 -0.996302  0.085921  71730   3  \n",
       "123936  0.319118 -1.197764 -0.133112 -0.273323 -0.931226  0.364441  90247   2  \n",
       "123979  0.274222 -0.747281  0.030830 -0.262655 -0.679647  0.733540  70730   2  \n",
       "124014  0.281278 -0.307351  0.138833 -0.232469 -0.302535  0.953138  71952   3  \n",
       "124037  0.278032  0.075864  0.214857 -0.159973  0.075792  0.997124  53822   3  \n",
       "124045  0.260063  0.162218  0.219961 -0.122334  0.161507  0.986872  74120   3  \n",
       "124086  0.260076  0.626081  0.251470 -0.010830  0.585973  0.810330  70944   3  \n",
       "124092  0.260203  0.794847  0.249862  0.031322  0.713756  0.700394  27967   5  \n",
       "124101  0.260667  0.796045  0.250258  0.031545  0.714595  0.699538  27967   5  \n",
       "124140  0.311726  1.513101  0.226414  0.193142  0.998336  0.057663  72240   3  \n",
       "124211  0.280944  2.112639  0.049052  0.265987  0.856760 -0.515716  90862   2  \n",
       "124217  0.267290  2.153130  0.027208  0.256787  0.835181 -0.549975  75389   2  \n",
       "124230  0.263416  2.347528 -0.024979  0.253499  0.713208 -0.700952  34168   5  \n",
       "124232  0.266890  2.349362 -0.023434  0.256797  0.711921 -0.702259  34168   5  \n",
       "124290  0.270551  2.816610 -0.136218  0.222823  0.319293 -0.947656  53886   3  \n",
       "124294  0.329786  3.108490 -0.197771  0.242852  0.033097 -0.999452  73576   3  \n",
       "124301  0.308415  2.957568 -0.163630  0.245119  0.182988 -0.983115  90918   3  \n",
       "124305  0.260428  2.904179 -0.155184  0.198577  0.235190 -0.971950  72526   7  \n",
       "124306  0.270633  2.973414 -0.169492  0.198787  0.167387 -0.985891  15552   4  \n",
       "124307  0.330454  3.108817 -0.197999  0.243404  0.032770 -0.999463  73576   3  \n",
       "124316  0.274381  3.190426 -0.209316  0.161867 -0.048814 -0.998808  90931   4  \n",
       "124325  0.274780  3.191593 -0.209737  0.161906 -0.049979 -0.998750  90931   4  \n",
       "124345  0.281318  3.510165 -0.253839  0.094349 -0.360284 -0.932843  74666   2  \n",
       "124362  0.259856  3.585336 -0.245136  0.056233 -0.429324 -0.903151  53912   2  \n",
       "\n",
       "[4280 rows x 23 columns]"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "labels.loc[cond]"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}