qn_upmu_signal_calc.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "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>Event</th>\n",
       "      <th>NmHitDAQ1</th>\n",
       "      <th>MinFracDay</th>\n",
       "      <th>JulDay</th>\n",
       "      <th>SSDB</th>\n",
       "      <th>eventDateTime</th>\n",
       "      <th>Hit1.ch1</th>\n",
       "      <th>PW1</th>\n",
       "      <th>Hit1.ch2</th>\n",
       "      <th>PW2</th>\n",
       "      <th>...</th>\n",
       "      <th>Hit1.ch4</th>\n",
       "      <th>PW4</th>\n",
       "      <th>Hit2.ch1</th>\n",
       "      <th>PW1.1</th>\n",
       "      <th>Hit2.ch2</th>\n",
       "      <th>PW2.1</th>\n",
       "      <th>Hit2.ch3</th>\n",
       "      <th>PW3.1</th>\n",
       "      <th>Hit2.ch4</th>\n",
       "      <th>PW4.1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>0.500012</td>\n",
       "      <td>2459747</td>\n",
       "      <td>43201.03</td>\n",
       "      <td>6/17/2022 0:0:1</td>\n",
       "      <td>0.00</td>\n",
       "      <td>17.50</td>\n",
       "      <td>2.50</td>\n",
       "      <td>25.00</td>\n",
       "      <td>...</td>\n",
       "      <td>3.0</td>\n",
       "      <td>25.01</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>4</td>\n",
       "      <td>0.500012</td>\n",
       "      <td>2459747</td>\n",
       "      <td>43201.07</td>\n",
       "      <td>6/17/2022 0:0:1</td>\n",
       "      <td>0.75</td>\n",
       "      <td>22.49</td>\n",
       "      <td>0.75</td>\n",
       "      <td>23.75</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>23.75</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2</td>\n",
       "      <td>0.500013</td>\n",
       "      <td>2459747</td>\n",
       "      <td>43201.16</td>\n",
       "      <td>6/17/2022 0:0:1</td>\n",
       "      <td>-1.00</td>\n",
       "      <td>-1.00</td>\n",
       "      <td>-1.00</td>\n",
       "      <td>-1.00</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>32.50</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3 rows × 22 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "   Event  NmHitDAQ1  MinFracDay   JulDay      SSDB    eventDateTime  Hit1.ch1  \\\n",
       "0      1          4    0.500012  2459747  43201.03  6/17/2022 0:0:1      0.00   \n",
       "1      2          4    0.500012  2459747  43201.07  6/17/2022 0:0:1      0.75   \n",
       "2      3          2    0.500013  2459747  43201.16  6/17/2022 0:0:1     -1.00   \n",
       "\n",
       "     PW1  Hit1.ch2    PW2  ...  Hit1.ch4    PW4  Hit2.ch1  PW1.1  Hit2.ch2  \\\n",
       "0  17.50      2.50  25.00  ...       3.0  25.01      -1.0   -1.0      -1.0   \n",
       "1  22.49      0.75  23.75  ...       0.0  23.75      -1.0   -1.0      -1.0   \n",
       "2  -1.00     -1.00  -1.00  ...       0.0  32.50      -1.0   -1.0      -1.0   \n",
       "\n",
       "   PW2.1  Hit2.ch3  PW3.1  Hit2.ch4  PW4.1  \n",
       "0   -1.0      -1.0   -1.0      -1.0   -1.0  \n",
       "1   -1.0      -1.0   -1.0      -1.0   -1.0  \n",
       "2   -1.0      -1.0   -1.0      -1.0   -1.0  \n",
       "\n",
       "[3 rows x 22 columns]"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "pd.options.mode.chained_assignment = None  # default='warn'\n",
    "import math\n",
    "import numpy as np\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "from IPython.core.interactiveshell import InteractiveShell\n",
    "InteractiveShell.ast_node_interactivity = \"all\"\n",
    "import math\n",
    "\n",
    "\n",
    "#read the file\n",
    "    \n",
    "data = pd.read_table('6119-hybridOut-June17-27.txt', skiprows=1)\n",
    "\n",
    "#drop columns with NaN\n",
    "\n",
    "data = data.dropna(axis='columns', how ='all')\n",
    "\n",
    "#fix the headers\n",
    "\n",
    "data = data.rename({'Event': 'Event', 'MinFracDay': 'NmHitDAQ1', \\\n",
    "\n",
    "                   'JulDay':'MinFracDay','Unnamed: 6':'JulDay', \\\n",
    "\n",
    "                   'Unnamed: 8':'SSDB', 'Unnamed: 10':'Hit1.ch1', \\\n",
    "\n",
    "                   'Unnamed: 13':'Hit1.ch2','Unnamed: 16':'Hit1.ch3', \\\n",
    "\n",
    "                   'Unnamed: 19':'Hit1.ch4','Unnamed: 22':'Hit2.ch1', \\\n",
    "\n",
    "                   'Unnamed: 25':'Hit2.ch2','Unnamed: 28':'Hit2.ch3', \\\n",
    "\n",
    "                   'Unnamed: 31':'Hit2.ch4'\n",
    "\n",
    "                   },\n",
    "\n",
    "                  axis=1)\n",
    "cols = data.columns\n",
    "\n",
    "cols = cols.map(lambda x: x.replace('#', '_') if isinstance(x, (str or unicode)) else x)\n",
    "\n",
    "data.columns = cols\n",
    "\n",
    "data.head(3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "data['Filter_1'] = data.iloc[:, 6] != -1  \n",
    "data['Filter_2'] = data.iloc[:, 8] != -1\n",
    "data['Filter_3'] = data.iloc[:, 10] != -1 \n",
    "data['Filter_4'] = data.iloc[:, 12] != -1 \n",
    "\n",
    "#data.head(5)\n",
    "\n",
    "filter_to_data = data.query('Filter_1 == True and Filter_3 == True and Filter_2 == True and Filter_4 == True')\n",
    "\n",
    "#Change the Hit Columns \n",
    "filter_to_data['Delta_T'] = (filter_to_data['Hit1.ch3'] - filter_to_data['Hit1.ch2'])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Limit 1 -3.3157375574237733 (nanoseconds)\n",
      "Limit 2 -13.050569036356471 (nanoseconds)\n"
     ]
    }
   ],
   "source": [
    "channel_1 = 2.179\n",
    "channel_2 = 1.697\n",
    "channel_3 = 0.032\n",
    "channel_4 = 0.758\n",
    "speed_of_light_constant = 0.298\n",
    "\n",
    "#Change channel numbers\n",
    "channel_distance = abs(channel_3 - channel_2)\n",
    "\n",
    "\n",
    "#Filter out outliers in calculations\n",
    "no_outlier_data = list(filter(lambda x: x < 30 and x > -30, filter_to_data['Delta_T'])) \n",
    "mean = np.mean(no_outlier_data)\n",
    "sigma = np.std(no_outlier_data)\n",
    "\n",
    "#Limit Calculations\n",
    "calculation_1 =(channel_distance / speed_of_light_constant)\n",
    "calculation_2 = round(calculation_1 + calculation_1)\n",
    "calculation_3 = mean - calculation_2\n",
    "\n",
    "\n",
    "#You Can Change Sigma Cut\n",
    "limit_1 = calculation_3 + (2 * sigma)\n",
    "limit_2 = calculation_3 - (2 * sigma)\n",
    "\n",
    "print(\"Limit 1 \" + str(limit_1) + \" (nanoseconds)\")\n",
    "print(\"Limit 2 \" + str(limit_2) + \" (nanoseconds)\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  },
  "vscode": {
   "interpreter": {
    "hash": "d254c79ca0780151aaa7e4cd16056ebdec783c71a8789cb97d8cd0ee929d8311"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}