Warning, /SimTracker/TrackerMaterialAnalysis/test/dEdxWeightsHFnose.ipynb is written in an unsupported language. File is not indexed.
0001 {
0002 "cells": [
0003 {
0004 "cell_type": "markdown",
0005 "metadata": {},
0006 "source": [
0007 "# dEdx weights calculation"
0008 ]
0009 },
0010 {
0011 "cell_type": "markdown",
0012 "metadata": {},
0013 "source": [
0014 "## Input"
0015 ]
0016 },
0017 {
0018 "cell_type": "markdown",
0019 "metadata": {},
0020 "source": [
0021 "This notebook expects to find a txt file called \"VolumesZPosition_HFnose.txt\". This file is coming from the _SimTracker/TrackerMaterialAnalysis_ package. The HFnose material is adapted for absorber/material as in V9. Just for the record we correct for the Materials with spaces in their name: \n",
0022 "\n",
0023 "sed -i -e 's/M_NEMA\\ FR4\\ plate/M_NEMA_FR4_plate/g' VolumesZPosition_HFnose.txt\n",
0024 "\n",
0025 "but ok now we use HGC_G10-FR4."
0026 ]
0027 },
0028 {
0029 "cell_type": "code",
0030 "execution_count": 1,
0031 "metadata": {},
0032 "outputs": [],
0033 "source": [
0034 "inputfile = \"VolumesZPosition_HFnose.txt\""
0035 ]
0036 },
0037 {
0038 "cell_type": "markdown",
0039 "metadata": {},
0040 "source": [
0041 "## Some imports"
0042 ]
0043 },
0044 {
0045 "cell_type": "code",
0046 "execution_count": 2,
0047 "metadata": {},
0048 "outputs": [],
0049 "source": [
0050 "from collections import OrderedDict\n",
0051 "import numpy as np\n",
0052 "import pandas as pd\n",
0053 "from IPython.display import display\n",
0054 "import matplotlib.pylab as plt\n",
0055 "import seaborn as sns"
0056 ]
0057 },
0058 {
0059 "cell_type": "markdown",
0060 "metadata": {},
0061 "source": [
0062 "## Materials Properties"
0063 ]
0064 },
0065 {
0066 "cell_type": "markdown",
0067 "metadata": {},
0068 "source": [
0069 "### Materials radiation lengths from Chris, Geant (in mm)"
0070 ]
0071 },
0072 {
0073 "cell_type": "code",
0074 "execution_count": 3,
0075 "metadata": {},
0076 "outputs": [
0077 {
0078 "name": "stdout",
0079 "output_type": "stream",
0080 "text": [
0081 "\u001b[1m Radiation Lengths in mm\n"
0082 ]
0083 },
0084 {
0085 "data": {
0086 "text/html": [
0087 "<div>\n",
0088 "<style scoped>\n",
0089 " .dataframe tbody tr th:only-of-type {\n",
0090 " vertical-align: middle;\n",
0091 " }\n",
0092 "\n",
0093 " .dataframe tbody tr th {\n",
0094 " vertical-align: top;\n",
0095 " }\n",
0096 "\n",
0097 " .dataframe thead th {\n",
0098 " text-align: right;\n",
0099 " }\n",
0100 "</style>\n",
0101 "<table border=\"1\" class=\"dataframe\">\n",
0102 " <thead>\n",
0103 " <tr style=\"text-align: right;\">\n",
0104 " <th></th>\n",
0105 " <th>FromChris</th>\n",
0106 " <th>FromGeant</th>\n",
0107 " </tr>\n",
0108 " </thead>\n",
0109 " <tbody>\n",
0110 " <tr>\n",
0111 " <th>Polyethylene</th>\n",
0112 " <td>503.100</td>\n",
0113 " <td>NaN</td>\n",
0114 " </tr>\n",
0115 " <tr>\n",
0116 " <th>Copper</th>\n",
0117 " <td>14.360</td>\n",
0118 " <td>14.3559</td>\n",
0119 " </tr>\n",
0120 " <tr>\n",
0121 " <th>H_Scintillator</th>\n",
0122 " <td>425.400</td>\n",
0123 " <td>425.3930</td>\n",
0124 " </tr>\n",
0125 " <tr>\n",
0126 " <th>Cables</th>\n",
0127 " <td>0.000</td>\n",
0128 " <td>66.7220</td>\n",
0129 " </tr>\n",
0130 " <tr>\n",
0131 " <th>HGC_G10-FR4</th>\n",
0132 " <td>175.000</td>\n",
0133 " <td>175.0560</td>\n",
0134 " </tr>\n",
0135 " <tr>\n",
0136 " <th>Silicon</th>\n",
0137 " <td>93.660</td>\n",
0138 " <td>93.6762</td>\n",
0139 " </tr>\n",
0140 " <tr>\n",
0141 " <th>Other</th>\n",
0142 " <td>0.000</td>\n",
0143 " <td>0.0000</td>\n",
0144 " </tr>\n",
0145 " <tr>\n",
0146 " <th>Air</th>\n",
0147 " <td>300000.000</td>\n",
0148 " <td>301522.0000</td>\n",
0149 " </tr>\n",
0150 " <tr>\n",
0151 " <th>StainlessSteel</th>\n",
0152 " <td>17.350</td>\n",
0153 " <td>17.3555</td>\n",
0154 " </tr>\n",
0155 " <tr>\n",
0156 " <th>WCu</th>\n",
0157 " <td>5.122</td>\n",
0158 " <td>5.1225</td>\n",
0159 " </tr>\n",
0160 " <tr>\n",
0161 " <th>Lead</th>\n",
0162 " <td>5.612</td>\n",
0163 " <td>5.6118</td>\n",
0164 " </tr>\n",
0165 " </tbody>\n",
0166 "</table>\n",
0167 "</div>"
0168 ],
0169 "text/plain": [
0170 " FromChris FromGeant\n",
0171 "Polyethylene 503.100 NaN\n",
0172 "Copper 14.360 14.3559\n",
0173 "H_Scintillator 425.400 425.3930\n",
0174 "Cables 0.000 66.7220\n",
0175 "HGC_G10-FR4 175.000 175.0560\n",
0176 "Silicon 93.660 93.6762\n",
0177 "Other 0.000 0.0000\n",
0178 "Air 300000.000 301522.0000\n",
0179 "StainlessSteel 17.350 17.3555\n",
0180 "WCu 5.122 5.1225\n",
0181 "Lead 5.612 5.6118"
0182 ]
0183 },
0184 "execution_count": 3,
0185 "metadata": {},
0186 "output_type": "execute_result"
0187 }
0188 ],
0189 "source": [
0190 "#In mm\n",
0191 "#-------- \n",
0192 "FromChrisMatXo = OrderedDict()\n",
0193 "FromChrisMatXo['Polyethylene'] = 503.100\n",
0194 "FromChrisMatXo['Copper'] = 14.36\n",
0195 "FromChrisMatXo['H_Scintillator'] = 425.4\n",
0196 "FromChrisMatXo['Cables'] = 0.\n",
0197 "#FromChrisMatXo['M_NEMA_FR4_plate'] = 175.\n",
0198 "#Just duplicate the above in this case\n",
0199 "FromChrisMatXo['HGC_G10-FR4'] = 175.\n",
0200 "FromChrisMatXo['Silicon'] = 93.66\n",
0201 "FromChrisMatXo['Other'] = 0.\n",
0202 "FromChrisMatXo['Air'] = 300000.\n",
0203 "FromChrisMatXo['StainlessSteel'] = 17.35\n",
0204 "FromChrisMatXo['WCu'] = 5.122\n",
0205 "FromChrisMatXo['Lead'] = 5.612\n",
0206 "#-------- \n",
0207 "FromGeantMatXo = OrderedDict()\n",
0208 "FromChrisMatXo['Polyethylene'] = 503.100\n",
0209 "FromGeantMatXo['Copper'] = 14.3559\n",
0210 "FromGeantMatXo['H_Scintillator'] = 425.393\n",
0211 "FromGeantMatXo['Cables'] = 66.722\n",
0212 "#FromGeantMatXo['M_NEMA_FR4_plate'] = 175.056\n",
0213 "FromGeantMatXo['HGC_G10-FR4'] = 175.056\n",
0214 "FromGeantMatXo['Silicon'] = 93.6762\n",
0215 "FromGeantMatXo['Other'] = 0.\n",
0216 "FromGeantMatXo['Air'] = 301522.\n",
0217 "FromGeantMatXo['StainlessSteel'] = 17.3555\n",
0218 "FromGeantMatXo['WCu'] = 5.1225\n",
0219 "FromGeantMatXo['Lead'] = 5.6118\n",
0220 "#-------- \n",
0221 "df = pd.DataFrame.from_dict([FromChrisMatXo, FromGeantMatXo])\n",
0222 "df = df.transpose()\n",
0223 "df.columns = ['FromChris', 'FromGeant']\n",
0224 "print( '\\033[1m' + ' Radiation Lengths in mm')\n",
0225 "df"
0226 ]
0227 },
0228 {
0229 "cell_type": "markdown",
0230 "metadata": {},
0231 "source": [
0232 "### Materials nuclear interaction lengths from Chris, Geant (in mm)"
0233 ]
0234 },
0235 {
0236 "cell_type": "code",
0237 "execution_count": 4,
0238 "metadata": {},
0239 "outputs": [
0240 {
0241 "name": "stdout",
0242 "output_type": "stream",
0243 "text": [
0244 "\u001b[1m Nuclear interaction Lengths in mm\n"
0245 ]
0246 },
0247 {
0248 "data": {
0249 "text/html": [
0250 "<div>\n",
0251 "<style scoped>\n",
0252 " .dataframe tbody tr th:only-of-type {\n",
0253 " vertical-align: middle;\n",
0254 " }\n",
0255 "\n",
0256 " .dataframe tbody tr th {\n",
0257 " vertical-align: top;\n",
0258 " }\n",
0259 "\n",
0260 " .dataframe thead th {\n",
0261 " text-align: right;\n",
0262 " }\n",
0263 "</style>\n",
0264 "<table border=\"1\" class=\"dataframe\">\n",
0265 " <thead>\n",
0266 " <tr style=\"text-align: right;\">\n",
0267 " <th></th>\n",
0268 " <th>FromChris</th>\n",
0269 " <th>FromGeant</th>\n",
0270 " </tr>\n",
0271 " </thead>\n",
0272 " <tbody>\n",
0273 " <tr>\n",
0274 " <th>Polyethylene</th>\n",
0275 " <td>881.8</td>\n",
0276 " <td>NaN</td>\n",
0277 " </tr>\n",
0278 " <tr>\n",
0279 " <th>Copper</th>\n",
0280 " <td>155.1</td>\n",
0281 " <td>155.880</td>\n",
0282 " </tr>\n",
0283 " <tr>\n",
0284 " <th>H_Scintillator</th>\n",
0285 " <td>701.3</td>\n",
0286 " <td>700.034</td>\n",
0287 " </tr>\n",
0288 " <tr>\n",
0289 " <th>Cables</th>\n",
0290 " <td>0.0</td>\n",
0291 " <td>393.710</td>\n",
0292 " </tr>\n",
0293 " <tr>\n",
0294 " <th>HGC_G10-FR4</th>\n",
0295 " <td>484.2</td>\n",
0296 " <td>483.429</td>\n",
0297 " </tr>\n",
0298 " <tr>\n",
0299 " <th>Silicon</th>\n",
0300 " <td>457.5</td>\n",
0301 " <td>456.628</td>\n",
0302 " </tr>\n",
0303 " <tr>\n",
0304 " <th>Other</th>\n",
0305 " <td>0.0</td>\n",
0306 " <td>0.000</td>\n",
0307 " </tr>\n",
0308 " <tr>\n",
0309 " <th>Air</th>\n",
0310 " <td>700000.0</td>\n",
0311 " <td>704083.000</td>\n",
0312 " </tr>\n",
0313 " <tr>\n",
0314 " <th>StainlessSteel</th>\n",
0315 " <td>166.0</td>\n",
0316 " <td>166.272</td>\n",
0317 " </tr>\n",
0318 " <tr>\n",
0319 " <th>WCu</th>\n",
0320 " <td>119.9</td>\n",
0321 " <td>120.105</td>\n",
0322 " </tr>\n",
0323 " <tr>\n",
0324 " <th>Lead</th>\n",
0325 " <td>182.6</td>\n",
0326 " <td>182.472</td>\n",
0327 " </tr>\n",
0328 " </tbody>\n",
0329 "</table>\n",
0330 "</div>"
0331 ],
0332 "text/plain": [
0333 " FromChris FromGeant\n",
0334 "Polyethylene 881.8 NaN\n",
0335 "Copper 155.1 155.880\n",
0336 "H_Scintillator 701.3 700.034\n",
0337 "Cables 0.0 393.710\n",
0338 "HGC_G10-FR4 484.2 483.429\n",
0339 "Silicon 457.5 456.628\n",
0340 "Other 0.0 0.000\n",
0341 "Air 700000.0 704083.000\n",
0342 "StainlessSteel 166.0 166.272\n",
0343 "WCu 119.9 120.105\n",
0344 "Lead 182.6 182.472"
0345 ]
0346 },
0347 "execution_count": 4,
0348 "metadata": {},
0349 "output_type": "execute_result"
0350 }
0351 ],
0352 "source": [
0353 "#In mm\n",
0354 "#-------- \n",
0355 "MatNucIntLength = OrderedDict()\n",
0356 "FromChrisMatNucIntLength = OrderedDict()\n",
0357 "FromChrisMatNucIntLength['Polyethylene'] = 881.8\n",
0358 "FromChrisMatNucIntLength['Copper'] = 155.1\n",
0359 "FromChrisMatNucIntLength['H_Scintillator'] = 701.3\n",
0360 "FromChrisMatNucIntLength['Cables'] = 0.\n",
0361 "#FromChrisMatNucIntLength['M_NEMA_FR4_plate'] = 484.2\n",
0362 "#Just duplicate the above in this case\n",
0363 "FromChrisMatNucIntLength['HGC_G10-FR4'] = 484.2\n",
0364 "FromChrisMatNucIntLength['Silicon'] = 457.5\n",
0365 "FromChrisMatNucIntLength['Other'] = 0.\n",
0366 "FromChrisMatNucIntLength['Air'] = 700000\n",
0367 "FromChrisMatNucIntLength['StainlessSteel'] = 166\n",
0368 "FromChrisMatNucIntLength['WCu'] = 119.9\n",
0369 "FromChrisMatNucIntLength['Lead'] = 182.6\n",
0370 "#-------- \n",
0371 "FromGeantMatNucIntLength = OrderedDict()\n",
0372 "FromChrisMatNucIntLength['Polyethylene'] = 881.8\n",
0373 "FromGeantMatNucIntLength['Copper'] = 155.88\n",
0374 "FromGeantMatNucIntLength['H_Scintillator'] = 700.034\n",
0375 "FromGeantMatNucIntLength['Cables'] = 393.71\n",
0376 "#FromGeantMatNucIntLength['M_NEMA_FR4_plate'] = 483.429\n",
0377 "FromGeantMatNucIntLength['HGC_G10-FR4'] = 483.429\n",
0378 "FromGeantMatNucIntLength['Silicon'] = 456.628\n",
0379 "FromGeantMatNucIntLength['Other'] = 0.\n",
0380 "FromGeantMatNucIntLength['Air'] = 704083\n",
0381 "FromGeantMatNucIntLength['StainlessSteel'] = 166.272\n",
0382 "FromGeantMatNucIntLength['WCu'] = 120.105\n",
0383 "FromGeantMatNucIntLength['Lead'] = 182.472\n",
0384 "#-------- \n",
0385 "df = pd.DataFrame.from_dict([FromChrisMatNucIntLength, FromGeantMatNucIntLength])\n",
0386 "df = df.transpose()\n",
0387 "df.columns = ['FromChris', 'FromGeant']\n",
0388 "print( '\\033[1m' + ' Nuclear interaction Lengths in mm')\n",
0389 "df"
0390 ]
0391 },
0392 {
0393 "cell_type": "markdown",
0394 "metadata": {},
0395 "source": [
0396 "### Materials dEdx from PDG, Chris, Geant (in MeV/mm)"
0397 ]
0398 },
0399 {
0400 "cell_type": "code",
0401 "execution_count": 5,
0402 "metadata": {},
0403 "outputs": [
0404 {
0405 "name": "stdout",
0406 "output_type": "stream",
0407 "text": [
0408 "\u001b[1m dEdx in MeV/mm\n"
0409 ]
0410 },
0411 {
0412 "data": {
0413 "text/html": [
0414 "<div>\n",
0415 "<style scoped>\n",
0416 " .dataframe tbody tr th:only-of-type {\n",
0417 " vertical-align: middle;\n",
0418 " }\n",
0419 "\n",
0420 " .dataframe tbody tr th {\n",
0421 " vertical-align: top;\n",
0422 " }\n",
0423 "\n",
0424 " .dataframe thead th {\n",
0425 " text-align: right;\n",
0426 " }\n",
0427 "</style>\n",
0428 "<table border=\"1\" class=\"dataframe\">\n",
0429 " <thead>\n",
0430 " <tr style=\"text-align: right;\">\n",
0431 " <th></th>\n",
0432 " <th>FromPDG</th>\n",
0433 " <th>FromChris</th>\n",
0434 " <th>FromGeantWithGetDEDX</th>\n",
0435 " <th>FromGeantWithComputeTotalDEDX</th>\n",
0436 " </tr>\n",
0437 " </thead>\n",
0438 " <tbody>\n",
0439 " <tr>\n",
0440 " <th>Polyethylene</th>\n",
0441 " <td>0.197027</td>\n",
0442 " <td>NaN</td>\n",
0443 " <td>NaN</td>\n",
0444 " <td>NaN</td>\n",
0445 " </tr>\n",
0446 " <tr>\n",
0447 " <th>Copper</th>\n",
0448 " <td>1.25709</td>\n",
0449 " <td>1.26</td>\n",
0450 " <td>1.14494 - 1.19191</td>\n",
0451 " <td>1.25937 - 1.51496</td>\n",
0452 " </tr>\n",
0453 " <tr>\n",
0454 " <th>H_Scintillator</th>\n",
0455 " <td>0.200135</td>\n",
0456 " <td>0.395</td>\n",
0457 " <td>0.155602 - 0.156968</td>\n",
0458 " <td>0.212843 - 0.230636</td>\n",
0459 " </tr>\n",
0460 " <tr>\n",
0461 " <th>Silicon</th>\n",
0462 " <td>0.387546</td>\n",
0463 " <td>0.388</td>\n",
0464 " <td>0.298729 - 0.349232</td>\n",
0465 " <td>0.38849 - 0.442262</td>\n",
0466 " </tr>\n",
0467 " <tr>\n",
0468 " <th>HGC_G10-FR4</th>\n",
0469 " <td>0.320491</td>\n",
0470 " <td>0.322</td>\n",
0471 " <td>0.200079 - 0.203959</td>\n",
0472 " <td>0.319344 - 0.364818</td>\n",
0473 " </tr>\n",
0474 " <tr>\n",
0475 " <th>Other</th>\n",
0476 " <td>0</td>\n",
0477 " <td>0</td>\n",
0478 " <td>0</td>\n",
0479 " <td>0</td>\n",
0480 " </tr>\n",
0481 " <tr>\n",
0482 " <th>Air</th>\n",
0483 " <td>0</td>\n",
0484 " <td>0</td>\n",
0485 " <td>0.000162487 - 0.000200878</td>\n",
0486 " <td>0.000221838 - 0.000303024</td>\n",
0487 " </tr>\n",
0488 " <tr>\n",
0489 " <th>StainlessSteel</th>\n",
0490 " <td>1.16574</td>\n",
0491 " <td>1.14</td>\n",
0492 " <td>1.17868 - 1.24005</td>\n",
0493 " <td>1.18145 - 1.40497</td>\n",
0494 " </tr>\n",
0495 " <tr>\n",
0496 " <th>WCu</th>\n",
0497 " <td>1.81171</td>\n",
0498 " <td>1.81</td>\n",
0499 " <td>1.4676 - 1.54044</td>\n",
0500 " <td>1.85741 - 2.18222</td>\n",
0501 " </tr>\n",
0502 " <tr>\n",
0503 " <th>Lead</th>\n",
0504 " <td>1.27347</td>\n",
0505 " <td>1.27</td>\n",
0506 " <td>1.28276 - 1.39512</td>\n",
0507 " <td>1.28388 - 1.52941</td>\n",
0508 " </tr>\n",
0509 " <tr>\n",
0510 " <th>Cables</th>\n",
0511 " <td>0.426449</td>\n",
0512 " <td>0</td>\n",
0513 " <td>0.334163 - 0.343928</td>\n",
0514 " <td>0.458855 - 0.507352</td>\n",
0515 " </tr>\n",
0516 " </tbody>\n",
0517 "</table>\n",
0518 "</div>"
0519 ],
0520 "text/plain": [
0521 " FromPDG FromChris FromGeantWithGetDEDX \\\n",
0522 "Polyethylene 0.197027 NaN NaN \n",
0523 "Copper 1.25709 1.26 1.14494 - 1.19191 \n",
0524 "H_Scintillator 0.200135 0.395 0.155602 - 0.156968 \n",
0525 "Silicon 0.387546 0.388 0.298729 - 0.349232 \n",
0526 "HGC_G10-FR4 0.320491 0.322 0.200079 - 0.203959 \n",
0527 "Other 0 0 0 \n",
0528 "Air 0 0 0.000162487 - 0.000200878 \n",
0529 "StainlessSteel 1.16574 1.14 1.17868 - 1.24005 \n",
0530 "WCu 1.81171 1.81 1.4676 - 1.54044 \n",
0531 "Lead 1.27347 1.27 1.28276 - 1.39512 \n",
0532 "Cables 0.426449 0 0.334163 - 0.343928 \n",
0533 "\n",
0534 " FromGeantWithComputeTotalDEDX \n",
0535 "Polyethylene NaN \n",
0536 "Copper 1.25937 - 1.51496 \n",
0537 "H_Scintillator 0.212843 - 0.230636 \n",
0538 "Silicon 0.38849 - 0.442262 \n",
0539 "HGC_G10-FR4 0.319344 - 0.364818 \n",
0540 "Other 0 \n",
0541 "Air 0.000221838 - 0.000303024 \n",
0542 "StainlessSteel 1.18145 - 1.40497 \n",
0543 "WCu 1.85741 - 2.18222 \n",
0544 "Lead 1.28388 - 1.52941 \n",
0545 "Cables 0.458855 - 0.507352 "
0546 ]
0547 },
0548 "execution_count": 5,
0549 "metadata": {},
0550 "output_type": "execute_result"
0551 }
0552 ],
0553 "source": [
0554 "#First in MeV gr^-1 cm^2 and at the end will convert to MeV/mm\n",
0555 "dEdx = OrderedDict()\n",
0556 "#In gr/cm^3 \n",
0557 "rho = OrderedDict()\n",
0558 "#--------\n",
0559 "#Some elements necessary to build our materials\n",
0560 "#Note: Whenever in PDG values (2018) are given in more than one state (e.g gas, liquid)\n",
0561 "#we go to the materials.xml file and choose the corresponding value. \n",
0562 "dEdx['Fe'] = 1.451\n",
0563 "dEdx['Mn'] = 1.428\n",
0564 "dEdx['Cr'] = 1.456\n",
0565 "dEdx['Ni'] = 1.468\n",
0566 "dEdx['C'] = 1.745\n",
0567 "dEdx['0'] = 1.801 \n",
0568 "dEdx['H'] = 4.034\n",
0569 "dEdx['Br'] = 1.380\n",
0570 "dEdx['W'] = 1.145\n",
0571 "dEdx['N'] = 1.813\n",
0572 "dEdx['Ar'] = 1.519\n",
0573 "\n",
0574 "#-------- \n",
0575 "rho['Fe'] = 7.874 \n",
0576 "rho['Mn'] = 7.440 \n",
0577 "rho['Cr'] = 7.180\n",
0578 "rho['Ni'] = 8.902 \n",
0579 "rho['C'] = 2.265 \n",
0580 "rho['0'] = 1.43**-3 #Here we choose the materials.xml value.\n",
0581 "rho['H'] = 0.07080\n",
0582 "rho['Br'] = 3.103 \n",
0583 "rho['W'] = 19.30 \n",
0584 "rho['N'] = 0.8070 \n",
0585 "rho['Ar'] = 1.639**-3 #Here we choose the materials.xml value. \n",
0586 "\n",
0587 "#--------\n",
0588 "#Detector materials (From PDG and Geometry/CMSCommonData/data/materials.xml)\n",
0589 "#You need a measurement of the density for a compound which we take from the same file\n",
0590 "dEdx['Polyethylene'] = 0.85628451*dEdx['C'] + 0.14371549*dEdx['H'] \n",
0591 "rho['Polyethylene'] = 0.95\n",
0592 "\n",
0593 "dEdx['Copper'] = 1.403\n",
0594 "rho['Copper'] = 8.960\n",
0595 "\n",
0596 "dEdx['H_Scintillator'] = 0.91512109*dEdx['C'] + 0.084878906*dEdx['H']\n",
0597 "rho['H_Scintillator'] = 1.032\n",
0598 "\n",
0599 "dEdx['Silicon'] = 1.664\n",
0600 "rho['Silicon'] = 2.329\n",
0601 "\n",
0602 "#dEdx['M_NEMA_FR4_plate'] = 0.18077359*dEdx['Silicon'] + 0.4056325*dEdx['0'] + 0.27804208*dEdx['C'] + 0.068442752*dEdx['H'] + 0.067109079*dEdx['Br']\n",
0603 "#rho['M_NEMA_FR4_plate'] = 1.025\n",
0604 "dEdx['HGC_G10-FR4'] = 0.18077359*dEdx['Silicon'] + 0.4056325*dEdx['0'] + 0.27804208*dEdx['C'] + 0.068442752*dEdx['H'] + 0.067109079*dEdx['Br']\n",
0605 "rho['HGC_G10-FR4'] = 1.70\n",
0606 "\n",
0607 "\n",
0608 "dEdx['Other'] = 0.\n",
0609 "rho['Other'] = 0.\n",
0610 "\n",
0611 "#dEdx['Air'] = 0.7494*dEdx['N'] + 0.2369*dEdx['0'] + 0.0129*dEdx['Ar'] + 0.0008*dEdx['H']\n",
0612 "dEdx['Air'] = 0.\n",
0613 "rho['Air'] = 1.214**-3\n",
0614 "\n",
0615 "dEdx['StainlessSteel'] = 0.6996*dEdx['Fe']+0.01*dEdx['Mn']+0.19*dEdx['Cr']+0.1*dEdx['Ni']+0.0004*dEdx['C'];\n",
0616 "rho['StainlessSteel'] = 8.02\n",
0617 "\n",
0618 "dEdx['WCu'] = 0.75*dEdx['W']+0.25*dEdx['Copper']\n",
0619 "rho['WCu'] = 14.979\n",
0620 "\n",
0621 "dEdx['Lead'] = 1.122 #Pb\n",
0622 "rho['Lead'] = 11.35 #Pb\n",
0623 "\n",
0624 "dEdx['Cables'] = 0.586*dEdx['Copper'] + 0.259*dEdx['C'] + 0.138*dEdx['0'] + 0.017*dEdx['H']\n",
0625 "rho['Cables'] = 2.68\n",
0626 "#--------\n",
0627 "# Now to the calculation. First, we will loop through the composite materials, not \n",
0628 "# to mess with the elements. \n",
0629 "composite_materials = ['Cables','StainlessSteel','M_NEMA_FR4_plate','H_Scintillator', 'Air', 'WCu']\n",
0630 "for element in dEdx: \n",
0631 " if element not in composite_materials: continue\n",
0632 " dEdx[element] = (rho[element] * dEdx[element]) / 10.\n",
0633 "#And for the rest\n",
0634 "for element in dEdx: \n",
0635 " if element in composite_materials: continue\n",
0636 " dEdx[element] = (rho[element] * dEdx[element]) / 10.\n",
0637 "\n",
0638 "\n",
0639 "#--------\n",
0640 "FromChrisdEdx = OrderedDict()\n",
0641 "FromChrisdEdx['Copper'] = 1.26\n",
0642 "FromChrisdEdx['H_Scintillator'] = 0.395\n",
0643 "FromChrisdEdx['Cables'] = 0.\n",
0644 "#FromChrisdEdx['M_NEMA_FR4_plate'] = 0.322\n",
0645 "#Just duplicate the above\n",
0646 "FromChrisdEdx['HGC_G10-FR4'] = 0.322\n",
0647 "FromChrisdEdx['Silicon'] = 0.388\n",
0648 "FromChrisdEdx['Other'] = 0.\n",
0649 "FromChrisdEdx['Air'] = 0.\n",
0650 "FromChrisdEdx['StainlessSteel'] = 1.14\n",
0651 "FromChrisdEdx['WCu'] = 1.81\n",
0652 "FromChrisdEdx['Lead'] = 1.27\n",
0653 "#-------- \n",
0654 "#The range below is from commands like: \n",
0655 "# Always recheck columns in txt, so that array is pointing to desired column. \n",
0656 "#array=($(cat VolumesZPosition_HFnose.txt | grep Copper | awk '{print $7}'))\n",
0657 "#IFS=$'\\n'\n",
0658 "#Highest value\n",
0659 "#echo \"${array[*]}\" | sort -nr | head -n1\n",
0660 "#Lowest value\n",
0661 "#echo \"${array[*]}\" | sort -nr | tail -n1\n",
0662 "\n",
0663 "FromGeantdEdxWithGetDEDX = OrderedDict()\n",
0664 "FromGeantdEdxWithGetDEDX['Copper'] = \"1.14494 - 1.19191\"\n",
0665 "FromGeantdEdxWithGetDEDX['H_Scintillator'] = \"0.155602 - 0.156968\"\n",
0666 "FromGeantdEdxWithGetDEDX['Cables'] = \"0.334163 - 0.343928\"\n",
0667 "#FromGeantdEdxWithGetDEDX['M_NEMA_FR4_plate'] = \"0.200079 - 0.203959\"\n",
0668 "FromGeantdEdxWithGetDEDX['HGC_G10-FR4'] = \"0.200079 - 0.203959\"\n",
0669 "FromGeantdEdxWithGetDEDX['Silicon'] = \"0.298729 - 0.349232\"\n",
0670 "FromGeantdEdxWithGetDEDX['Other'] = 0.\n",
0671 "FromGeantdEdxWithGetDEDX['Air'] = \"0.000162487 - 0.000200878\"\n",
0672 "FromGeantdEdxWithGetDEDX['StainlessSteel'] = \"1.17868 - 1.24005\"\n",
0673 "FromGeantdEdxWithGetDEDX['WCu'] = \"1.4676 - 1.54044\"\n",
0674 "FromGeantdEdxWithGetDEDX['Lead'] = \"1.28276 - 1.39512\"\n",
0675 "#-------- \n",
0676 "#The range below is from commands like: \n",
0677 "# Always recheck columns in txt, so that array is pointing to desired column. \n",
0678 "#array=($(cat VolumesZPosition_HFnose.txt | grep Copper | awk '{print $9}'))\n",
0679 "#IFS=$'\\n'\n",
0680 "#Highest value\n",
0681 "#echo \"${array[*]}\" | sort -nr | head -n1\n",
0682 "#Lowest value\n",
0683 "#echo \"${array[*]}\" | sort -nr | tail -n1\n",
0684 "FromGeantdEdxWithComputeTotalDEDX = OrderedDict()\n",
0685 "FromGeantdEdxWithComputeTotalDEDX['Copper'] = \"1.25937 - 1.51496\"\n",
0686 "FromGeantdEdxWithComputeTotalDEDX['H_Scintillator'] = \"0.212843 - 0.230636\"\n",
0687 "FromGeantdEdxWithComputeTotalDEDX['Cables'] = \"0.458855 - 0.507352\"\n",
0688 "#FromGeantdEdxWithComputeTotalDEDX['M_NEMA_FR4_plate'] = \"0.319344 - 0.364818\"\n",
0689 "FromGeantdEdxWithComputeTotalDEDX['HGC_G10-FR4'] = \"0.319344 - 0.364818\"\n",
0690 "FromGeantdEdxWithComputeTotalDEDX['Silicon'] = \"0.38849 - 0.442262\"\n",
0691 "FromGeantdEdxWithComputeTotalDEDX['Other'] = 0.\n",
0692 "FromGeantdEdxWithComputeTotalDEDX['Air'] = \"0.000221838 - 0.000303024\"\n",
0693 "FromGeantdEdxWithComputeTotalDEDX['StainlessSteel'] = \"1.18145 - 1.40497\"\n",
0694 "FromGeantdEdxWithComputeTotalDEDX['WCu'] = \"1.85741 - 2.18222\"\n",
0695 "FromGeantdEdxWithComputeTotalDEDX['Lead'] = \"1.28388 - 1.52941\"\n",
0696 "#-------- \n",
0697 "df = pd.DataFrame.from_dict([dEdx,FromChrisdEdx, FromGeantdEdxWithGetDEDX, FromGeantdEdxWithComputeTotalDEDX])\n",
0698 "df = df.transpose()[11:]\n",
0699 "df.columns = ['FromPDG','FromChris', 'FromGeantWithGetDEDX', 'FromGeantWithComputeTotalDEDX']\n",
0700 "print( '\\033[1m' + ' dEdx in MeV/mm')\n",
0701 "df"
0702 ]
0703 },
0704 {
0705 "cell_type": "markdown",
0706 "metadata": {},
0707 "source": [
0708 "### Choice for radiation lengths values -> Geant values"
0709 ]
0710 },
0711 {
0712 "cell_type": "code",
0713 "execution_count": 6,
0714 "metadata": {},
0715 "outputs": [],
0716 "source": [
0717 "# Current choice -> Geant values\n",
0718 "#In mm\n",
0719 "MatXo = OrderedDict()\n",
0720 "MatXo['Polyethylene'] =503.100\n",
0721 "MatXo['Copper'] = 14.3559\n",
0722 "MatXo['H_Scintillator'] = 425.393\n",
0723 "MatXo['Cables'] = 66.722\n",
0724 "#MatXo['M_NEMA_FR4_plate'] = 175.056\n",
0725 "MatXo['HGC_G10-FR4'] = 175.056\n",
0726 "MatXo['Silicon'] = 93.6762\n",
0727 "MatXo['Other'] = 0.\n",
0728 "MatXo['Air'] = 301522.\n",
0729 "MatXo['StainlessSteel'] = 17.3555\n",
0730 "MatXo['WCu'] = 5.1225\n",
0731 "MatXo['Lead'] = 5.6118"
0732 ]
0733 },
0734 {
0735 "cell_type": "code",
0736 "execution_count": 7,
0737 "metadata": {},
0738 "outputs": [],
0739 "source": [
0740 "# Let's read the input file\n",
0741 "# Mat is prepoint material, Z is post point material start, so \n",
0742 "# upper edge of prepoint material. Etable,Efull is the energy loss in the \n",
0743 "# prepoint volume with GetDEDX and ComputeTotalDEDX. \n",
0744 "matZ = pd.read_csv(inputfile, sep=\" \", header=None, names=[\"Mat\", \"Z\", \"Eta\", \"R\", \"Etable\", \"Efull\"], index_col=False)"
0745 ]
0746 },
0747 {
0748 "cell_type": "code",
0749 "execution_count": 8,
0750 "metadata": {},
0751 "outputs": [],
0752 "source": [
0753 "# We will make a new column with the physical thickness of the volumes in mm\n",
0754 "matZ[\"PhysThickInmm\"] = abs(matZ[\"Z\"].shift(-1) - matZ[\"Z\"])\n",
0755 "matZ[\"PhysThickInmm\"] = matZ[\"PhysThickInmm\"].shift(1)"
0756 ]
0757 },
0758 {
0759 "cell_type": "code",
0760 "execution_count": 9,
0761 "metadata": {},
0762 "outputs": [],
0763 "source": [
0764 "# We will add a column that indicates the track the relevant volume belongs to. \n",
0765 "# The logic is that right before the next track \"PhysThickInmm\" column will be \n",
0766 "# very large. \n",
0767 "matZ[\"trackflag\"] = matZ.apply(lambda row: True if row[\"PhysThickInmm\"] < 400. else False ,axis=1)\n",
0768 "matZ[\"tracknum\"] = (( matZ[\"trackflag\"] == False) & (matZ[\"trackflag\"].shift() == True)).cumsum()\n",
0769 "matZ[\"tracknum\"] = matZ[\"tracknum\"].shift(1)\n",
0770 "matZ = matZ.drop('trackflag', 1)"
0771 ]
0772 },
0773 {
0774 "cell_type": "code",
0775 "execution_count": 10,
0776 "metadata": {},
0777 "outputs": [],
0778 "source": [
0779 "# Now that we have the tracknum we will not let the last Copper volume to \n",
0780 "# its huge thickness due to track change. According to TDR BH back is at 5137.7 mm \n",
0781 "# so we put 6.0 mm for that Copper. In any case this do *not* effect the \n",
0782 "# dEdx weights calculation since it is after the sensitive material of the last layer. \n",
0783 "matZ.loc[ matZ[\"PhysThickInmm\"] > 400. , \"PhysThickInmm\" ] = 6.0"
0784 ]
0785 },
0786 {
0787 "cell_type": "code",
0788 "execution_count": 11,
0789 "metadata": {},
0790 "outputs": [],
0791 "source": [
0792 "# The line below is when comparing with Chris which has no Cables. \n",
0793 "# matZ = matZ.query('Mat != \"Cables\"')\n",
0794 "\n",
0795 "# Again, adding a new column with the physical thickness of the volumes in radiation lengths\n",
0796 "matZ[\"PhysThickInXo\"] = matZ.apply(lambda row: row[\"PhysThickInmm\"] / MatXo[row[\"Mat\"]],axis=1)"
0797 ]
0798 },
0799 {
0800 "cell_type": "code",
0801 "execution_count": 12,
0802 "metadata": {},
0803 "outputs": [],
0804 "source": [
0805 "# Adding a new column with the dEdx of the material that the volumes is build\n",
0806 "matZ[\"dEdx\"] = matZ.apply(lambda row: dEdx[row[\"Mat\"]],axis=1)"
0807 ]
0808 },
0809 {
0810 "cell_type": "code",
0811 "execution_count": 13,
0812 "metadata": {},
0813 "outputs": [],
0814 "source": [
0815 "# Another column with the dEdx times thickness to help us with the calculation of the \n",
0816 "# final dEdx weights \n",
0817 "matZ[\"dEdxtimesdx\"] = matZ[\"dEdx\"] * matZ[\"PhysThickInmm\"]"
0818 ]
0819 },
0820 {
0821 "cell_type": "code",
0822 "execution_count": 14,
0823 "metadata": {},
0824 "outputs": [],
0825 "source": [
0826 "# And here a column with the cumulative sum\n",
0827 "matZ[\"dEdxtimesdxCum\"] = matZ.groupby('tracknum')[\"dEdxtimesdx\"].cumsum()"
0828 ]
0829 },
0830 {
0831 "cell_type": "code",
0832 "execution_count": 15,
0833 "metadata": {},
0834 "outputs": [],
0835 "source": [
0836 "# And here a column with the cumulative sum for Etable\n",
0837 "matZ[\"EtableCum\"] = matZ.groupby('tracknum')[\"Etable\"].cumsum()"
0838 ]
0839 },
0840 {
0841 "cell_type": "code",
0842 "execution_count": 16,
0843 "metadata": {},
0844 "outputs": [],
0845 "source": [
0846 "# And here a column with the cumulative sum for Efull\n",
0847 "matZ[\"EfullCum\"] = matZ.groupby('tracknum')[\"Efull\"].cumsum()"
0848 ]
0849 },
0850 {
0851 "cell_type": "code",
0852 "execution_count": 17,
0853 "metadata": {},
0854 "outputs": [],
0855 "source": [
0856 "# And here a column for the cumulative sum in radiation length \n",
0857 "matZ[\"PhysThickInXoCum\"] = matZ.groupby('tracknum')[\"PhysThickInXo\"].cumsum()"
0858 ]
0859 },
0860 {
0861 "cell_type": "code",
0862 "execution_count": 18,
0863 "metadata": {},
0864 "outputs": [],
0865 "source": [
0866 "# We will add a column that indicates the layer the relevant volume belongs to. \n",
0867 "# The logic is that if the previous material is Silicon or Scintillator\n",
0868 "# we change layer.\n",
0869 "matZ[\"layerflag\"] = matZ.apply(lambda row: False if row[\"Mat\"] == \"Silicon\" or row[\"Mat\"] == \"H_Scintillator\" else True ,axis=1)\n",
0870 "matZ[\"layer\"] = ( matZ[\"layerflag\"] == True) & (matZ[\"layerflag\"].shift(1) == False) \n",
0871 "matZ[\"layer\"] = matZ.groupby('tracknum')[\"layer\"].cumsum()\n",
0872 "#The convention is that layers starts from 1\n",
0873 "matZ[\"layer\"] = matZ.apply(lambda row: row[\"layer\"] + 1,axis=1)"
0874 ]
0875 },
0876 {
0877 "cell_type": "code",
0878 "execution_count": 19,
0879 "metadata": {},
0880 "outputs": [],
0881 "source": [
0882 "# Drop auxillary columns\n",
0883 "#We need layerflag for not counting the silicon/scintillator in the dedx. \n",
0884 "#matZ = matZ.drop('layerflag', 1)"
0885 ]
0886 },
0887 {
0888 "cell_type": "code",
0889 "execution_count": 20,
0890 "metadata": {},
0891 "outputs": [],
0892 "source": [
0893 "# This was my misunderstanding. There should be 53 in the layers column since \n",
0894 "# after the last scintillator or silicon the index will increase, there is \n",
0895 "# material there. In case we want to filter them out uncomment the following. \n",
0896 "# Be careful! Filter chooses and not disgards!\n",
0897 "# matZ = matZ.groupby('tracknum').filter(lambda g: ~(g['layer'] == 53.0).any() ) "
0898 ]
0899 },
0900 {
0901 "cell_type": "code",
0902 "execution_count": 21,
0903 "metadata": {},
0904 "outputs": [
0905 {
0906 "data": {
0907 "text/html": [
0908 "<div>\n",
0909 "<style scoped>\n",
0910 " .dataframe tbody tr th:only-of-type {\n",
0911 " vertical-align: middle;\n",
0912 " }\n",
0913 "\n",
0914 " .dataframe tbody tr th {\n",
0915 " vertical-align: top;\n",
0916 " }\n",
0917 "\n",
0918 " .dataframe thead th {\n",
0919 " text-align: right;\n",
0920 " }\n",
0921 "</style>\n",
0922 "<table border=\"1\" class=\"dataframe\">\n",
0923 " <thead>\n",
0924 " <tr style=\"text-align: right;\">\n",
0925 " <th></th>\n",
0926 " <th>Mat</th>\n",
0927 " <th>Z</th>\n",
0928 " <th>Eta</th>\n",
0929 " <th>R</th>\n",
0930 " <th>Etable</th>\n",
0931 " <th>Efull</th>\n",
0932 " <th>PhysThickInmm</th>\n",
0933 " <th>tracknum</th>\n",
0934 " <th>PhysThickInXo</th>\n",
0935 " <th>dEdx</th>\n",
0936 " <th>dEdxtimesdx</th>\n",
0937 " <th>dEdxtimesdxCum</th>\n",
0938 " <th>EtableCum</th>\n",
0939 " <th>EfullCum</th>\n",
0940 " <th>PhysThickInXoCum</th>\n",
0941 " <th>layerflag</th>\n",
0942 " <th>layer</th>\n",
0943 " </tr>\n",
0944 " </thead>\n",
0945 " <tbody>\n",
0946 " <tr>\n",
0947 " <th>71</th>\n",
0948 " <td>Polyethylene</td>\n",
0949 " <td>10484.0</td>\n",
0950 " <td>3.43512</td>\n",
0951 " <td>676.324</td>\n",
0952 " <td>10.293100</td>\n",
0953 " <td>0</td>\n",
0954 " <td>0.0</td>\n",
0955 " <td>1.0</td>\n",
0956 " <td>0.000000</td>\n",
0957 " <td>0.197027</td>\n",
0958 " <td>0.000000</td>\n",
0959 " <td>0.000000</td>\n",
0960 " <td>10.293100</td>\n",
0961 " <td>0</td>\n",
0962 " <td>0.000000</td>\n",
0963 " <td>True</td>\n",
0964 " <td>1.0</td>\n",
0965 " </tr>\n",
0966 " <tr>\n",
0967 " <th>72</th>\n",
0968 " <td>StainlessSteel</td>\n",
0969 " <td>10485.0</td>\n",
0970 " <td>3.43512</td>\n",
0971 " <td>676.388</td>\n",
0972 " <td>1.414020</td>\n",
0973 " <td>0</td>\n",
0974 " <td>1.0</td>\n",
0975 " <td>1.0</td>\n",
0976 " <td>0.057619</td>\n",
0977 " <td>1.165737</td>\n",
0978 " <td>1.165737</td>\n",
0979 " <td>1.165737</td>\n",
0980 " <td>11.707120</td>\n",
0981 " <td>0</td>\n",
0982 " <td>0.057619</td>\n",
0983 " <td>True</td>\n",
0984 " <td>1.0</td>\n",
0985 " </tr>\n",
0986 " <tr>\n",
0987 " <th>73</th>\n",
0988 " <td>Lead</td>\n",
0989 " <td>10506.0</td>\n",
0990 " <td>3.43512</td>\n",
0991 " <td>677.743</td>\n",
0992 " <td>30.729100</td>\n",
0993 " <td>0</td>\n",
0994 " <td>21.0</td>\n",
0995 " <td>1.0</td>\n",
0996 " <td>3.742115</td>\n",
0997 " <td>1.273470</td>\n",
0998 " <td>26.742870</td>\n",
0999 " <td>27.908607</td>\n",
1000 " <td>42.436220</td>\n",
1001 " <td>0</td>\n",
1002 " <td>3.799733</td>\n",
1003 " <td>True</td>\n",
1004 " <td>1.0</td>\n",
1005 " </tr>\n",
1006 " <tr>\n",
1007 " <th>74</th>\n",
1008 " <td>StainlessSteel</td>\n",
1009 " <td>10507.0</td>\n",
1010 " <td>3.43512</td>\n",
1011 " <td>677.807</td>\n",
1012 " <td>1.413140</td>\n",
1013 " <td>0</td>\n",
1014 " <td>1.0</td>\n",
1015 " <td>1.0</td>\n",
1016 " <td>0.057619</td>\n",
1017 " <td>1.165737</td>\n",
1018 " <td>1.165737</td>\n",
1019 " <td>29.074344</td>\n",
1020 " <td>43.849360</td>\n",
1021 " <td>0</td>\n",
1022 " <td>3.857352</td>\n",
1023 " <td>True</td>\n",
1024 " <td>1.0</td>\n",
1025 " </tr>\n",
1026 " <tr>\n",
1027 " <th>75</th>\n",
1028 " <td>Air</td>\n",
1029 " <td>10510.3</td>\n",
1030 " <td>3.43512</td>\n",
1031 " <td>678.024</td>\n",
1032 " <td>0.000685</td>\n",
1033 " <td>0</td>\n",
1034 " <td>3.3</td>\n",
1035 " <td>1.0</td>\n",
1036 " <td>0.000011</td>\n",
1037 " <td>0.000000</td>\n",
1038 " <td>0.000000</td>\n",
1039 " <td>29.074344</td>\n",
1040 " <td>43.850045</td>\n",
1041 " <td>0</td>\n",
1042 " <td>3.857363</td>\n",
1043 " <td>True</td>\n",
1044 " <td>1.0</td>\n",
1045 " </tr>\n",
1046 " <tr>\n",
1047 " <th>76</th>\n",
1048 " <td>HGC_G10-FR4</td>\n",
1049 " <td>10511.9</td>\n",
1050 " <td>3.43512</td>\n",
1051 " <td>678.127</td>\n",
1052 " <td>0.523805</td>\n",
1053 " <td>0</td>\n",
1054 " <td>1.6</td>\n",
1055 " <td>1.0</td>\n",
1056 " <td>0.009140</td>\n",
1057 " <td>0.320491</td>\n",
1058 " <td>0.512786</td>\n",
1059 " <td>29.587131</td>\n",
1060 " <td>44.373850</td>\n",
1061 " <td>0</td>\n",
1062 " <td>3.866503</td>\n",
1063 " <td>True</td>\n",
1064 " <td>1.0</td>\n",
1065 " </tr>\n",
1066 " <tr>\n",
1067 " <th>77</th>\n",
1068 " <td>Silicon</td>\n",
1069 " <td>10512.1</td>\n",
1070 " <td>3.43512</td>\n",
1071 " <td>678.135</td>\n",
1072 " <td>0.142698</td>\n",
1073 " <td>0</td>\n",
1074 " <td>0.2</td>\n",
1075 " <td>1.0</td>\n",
1076 " <td>0.002135</td>\n",
1077 " <td>0.387546</td>\n",
1078 " <td>0.077509</td>\n",
1079 " <td>29.664640</td>\n",
1080 " <td>44.516548</td>\n",
1081 " <td>0</td>\n",
1082 " <td>3.868638</td>\n",
1083 " <td>False</td>\n",
1084 " <td>1.0</td>\n",
1085 " </tr>\n",
1086 " <tr>\n",
1087 " <th>78</th>\n",
1088 " <td>Silicon</td>\n",
1089 " <td>10512.2</td>\n",
1090 " <td>3.43512</td>\n",
1091 " <td>678.146</td>\n",
1092 " <td>0.214046</td>\n",
1093 " <td>0</td>\n",
1094 " <td>0.1</td>\n",
1095 " <td>1.0</td>\n",
1096 " <td>0.001068</td>\n",
1097 " <td>0.387546</td>\n",
1098 " <td>0.038755</td>\n",
1099 " <td>29.703394</td>\n",
1100 " <td>44.730594</td>\n",
1101 " <td>0</td>\n",
1102 " <td>3.869705</td>\n",
1103 " <td>False</td>\n",
1104 " <td>1.0</td>\n",
1105 " </tr>\n",
1106 " <tr>\n",
1107 " <th>79</th>\n",
1108 " <td>WCu</td>\n",
1109 " <td>10513.6</td>\n",
1110 " <td>3.43512</td>\n",
1111 " <td>678.236</td>\n",
1112 " <td>2.590920</td>\n",
1113 " <td>0</td>\n",
1114 " <td>1.4</td>\n",
1115 " <td>1.0</td>\n",
1116 " <td>0.273304</td>\n",
1117 " <td>1.811710</td>\n",
1118 " <td>2.536394</td>\n",
1119 " <td>32.239788</td>\n",
1120 " <td>47.321514</td>\n",
1121 " <td>0</td>\n",
1122 " <td>4.143010</td>\n",
1123 " <td>True</td>\n",
1124 " <td>2.0</td>\n",
1125 " </tr>\n",
1126 " <tr>\n",
1127 " <th>80</th>\n",
1128 " <td>Copper</td>\n",
1129 " <td>10519.6</td>\n",
1130 " <td>3.43512</td>\n",
1131 " <td>678.624</td>\n",
1132 " <td>8.913600</td>\n",
1133 " <td>0</td>\n",
1134 " <td>6.0</td>\n",
1135 " <td>1.0</td>\n",
1136 " <td>0.417947</td>\n",
1137 " <td>1.257088</td>\n",
1138 " <td>7.542528</td>\n",
1139 " <td>39.782316</td>\n",
1140 " <td>56.235114</td>\n",
1141 " <td>0</td>\n",
1142 " <td>4.560956</td>\n",
1143 " <td>True</td>\n",
1144 " <td>2.0</td>\n",
1145 " </tr>\n",
1146 " <tr>\n",
1147 " <th>81</th>\n",
1148 " <td>Air</td>\n",
1149 " <td>10523.0</td>\n",
1150 " <td>3.43512</td>\n",
1151 " <td>678.840</td>\n",
1152 " <td>0.000685</td>\n",
1153 " <td>0</td>\n",
1154 " <td>3.4</td>\n",
1155 " <td>1.0</td>\n",
1156 " <td>0.000011</td>\n",
1157 " <td>0.000000</td>\n",
1158 " <td>0.000000</td>\n",
1159 " <td>39.782316</td>\n",
1160 " <td>56.235799</td>\n",
1161 " <td>0</td>\n",
1162 " <td>4.560967</td>\n",
1163 " <td>True</td>\n",
1164 " <td>2.0</td>\n",
1165 " </tr>\n",
1166 " <tr>\n",
1167 " <th>82</th>\n",
1168 " <td>StainlessSteel</td>\n",
1169 " <td>10524.0</td>\n",
1170 " <td>3.43512</td>\n",
1171 " <td>678.904</td>\n",
1172 " <td>1.412470</td>\n",
1173 " <td>0</td>\n",
1174 " <td>1.0</td>\n",
1175 " <td>1.0</td>\n",
1176 " <td>0.057619</td>\n",
1177 " <td>1.165737</td>\n",
1178 " <td>1.165737</td>\n",
1179 " <td>40.948053</td>\n",
1180 " <td>57.648269</td>\n",
1181 " <td>0</td>\n",
1182 " <td>4.618586</td>\n",
1183 " <td>True</td>\n",
1184 " <td>2.0</td>\n",
1185 " </tr>\n",
1186 " <tr>\n",
1187 " <th>83</th>\n",
1188 " <td>Lead</td>\n",
1189 " <td>10545.0</td>\n",
1190 " <td>3.43512</td>\n",
1191 " <td>680.259</td>\n",
1192 " <td>30.698500</td>\n",
1193 " <td>0</td>\n",
1194 " <td>21.0</td>\n",
1195 " <td>1.0</td>\n",
1196 " <td>3.742115</td>\n",
1197 " <td>1.273470</td>\n",
1198 " <td>26.742870</td>\n",
1199 " <td>67.690923</td>\n",
1200 " <td>88.346769</td>\n",
1201 " <td>0</td>\n",
1202 " <td>8.360701</td>\n",
1203 " <td>True</td>\n",
1204 " <td>2.0</td>\n",
1205 " </tr>\n",
1206 " <tr>\n",
1207 " <th>84</th>\n",
1208 " <td>StainlessSteel</td>\n",
1209 " <td>10546.0</td>\n",
1210 " <td>3.43512</td>\n",
1211 " <td>680.323</td>\n",
1212 " <td>1.411610</td>\n",
1213 " <td>0</td>\n",
1214 " <td>1.0</td>\n",
1215 " <td>1.0</td>\n",
1216 " <td>0.057619</td>\n",
1217 " <td>1.165737</td>\n",
1218 " <td>1.165737</td>\n",
1219 " <td>68.856661</td>\n",
1220 " <td>89.758379</td>\n",
1221 " <td>0</td>\n",
1222 " <td>8.418320</td>\n",
1223 " <td>True</td>\n",
1224 " <td>2.0</td>\n",
1225 " </tr>\n",
1226 " <tr>\n",
1227 " <th>85</th>\n",
1228 " <td>Air</td>\n",
1229 " <td>10549.3</td>\n",
1230 " <td>3.43512</td>\n",
1231 " <td>680.539</td>\n",
1232 " <td>0.000684</td>\n",
1233 " <td>0</td>\n",
1234 " <td>3.3</td>\n",
1235 " <td>1.0</td>\n",
1236 " <td>0.000011</td>\n",
1237 " <td>0.000000</td>\n",
1238 " <td>0.000000</td>\n",
1239 " <td>68.856661</td>\n",
1240 " <td>89.759063</td>\n",
1241 " <td>0</td>\n",
1242 " <td>8.418330</td>\n",
1243 " <td>True</td>\n",
1244 " <td>2.0</td>\n",
1245 " </tr>\n",
1246 " <tr>\n",
1247 " <th>86</th>\n",
1248 " <td>Copper</td>\n",
1249 " <td>10555.3</td>\n",
1250 " <td>3.43512</td>\n",
1251 " <td>680.927</td>\n",
1252 " <td>8.903040</td>\n",
1253 " <td>0</td>\n",
1254 " <td>6.0</td>\n",
1255 " <td>1.0</td>\n",
1256 " <td>0.417947</td>\n",
1257 " <td>1.257088</td>\n",
1258 " <td>7.542528</td>\n",
1259 " <td>76.399189</td>\n",
1260 " <td>98.662103</td>\n",
1261 " <td>0</td>\n",
1262 " <td>8.836277</td>\n",
1263 " <td>True</td>\n",
1264 " <td>2.0</td>\n",
1265 " </tr>\n",
1266 " <tr>\n",
1267 " <th>87</th>\n",
1268 " <td>WCu</td>\n",
1269 " <td>10556.8</td>\n",
1270 " <td>3.43512</td>\n",
1271 " <td>681.017</td>\n",
1272 " <td>2.587410</td>\n",
1273 " <td>0</td>\n",
1274 " <td>1.5</td>\n",
1275 " <td>1.0</td>\n",
1276 " <td>0.292826</td>\n",
1277 " <td>1.811710</td>\n",
1278 " <td>2.717565</td>\n",
1279 " <td>79.116754</td>\n",
1280 " <td>101.249513</td>\n",
1281 " <td>0</td>\n",
1282 " <td>9.129103</td>\n",
1283 " <td>True</td>\n",
1284 " <td>2.0</td>\n",
1285 " </tr>\n",
1286 " <tr>\n",
1287 " <th>88</th>\n",
1288 " <td>Silicon</td>\n",
1289 " <td>10556.9</td>\n",
1290 " <td>3.43512</td>\n",
1291 " <td>681.028</td>\n",
1292 " <td>0.214043</td>\n",
1293 " <td>0</td>\n",
1294 " <td>0.1</td>\n",
1295 " <td>1.0</td>\n",
1296 " <td>0.001068</td>\n",
1297 " <td>0.387546</td>\n",
1298 " <td>0.038755</td>\n",
1299 " <td>79.155508</td>\n",
1300 " <td>101.463556</td>\n",
1301 " <td>0</td>\n",
1302 " <td>9.130170</td>\n",
1303 " <td>False</td>\n",
1304 " <td>2.0</td>\n",
1305 " </tr>\n",
1306 " <tr>\n",
1307 " <th>89</th>\n",
1308 " <td>Silicon</td>\n",
1309 " <td>10557.0</td>\n",
1310 " <td>3.43512</td>\n",
1311 " <td>681.036</td>\n",
1312 " <td>0.142696</td>\n",
1313 " <td>0</td>\n",
1314 " <td>0.1</td>\n",
1315 " <td>1.0</td>\n",
1316 " <td>0.001068</td>\n",
1317 " <td>0.387546</td>\n",
1318 " <td>0.038755</td>\n",
1319 " <td>79.194263</td>\n",
1320 " <td>101.606252</td>\n",
1321 " <td>0</td>\n",
1322 " <td>9.131238</td>\n",
1323 " <td>False</td>\n",
1324 " <td>2.0</td>\n",
1325 " </tr>\n",
1326 " <tr>\n",
1327 " <th>90</th>\n",
1328 " <td>HGC_G10-FR4</td>\n",
1329 " <td>10558.6</td>\n",
1330 " <td>3.43512</td>\n",
1331 " <td>681.139</td>\n",
1332 " <td>0.522997</td>\n",
1333 " <td>0</td>\n",
1334 " <td>1.6</td>\n",
1335 " <td>1.0</td>\n",
1336 " <td>0.009140</td>\n",
1337 " <td>0.320491</td>\n",
1338 " <td>0.512786</td>\n",
1339 " <td>79.707049</td>\n",
1340 " <td>102.129249</td>\n",
1341 " <td>0</td>\n",
1342 " <td>9.140378</td>\n",
1343 " <td>True</td>\n",
1344 " <td>3.0</td>\n",
1345 " </tr>\n",
1346 " <tr>\n",
1347 " <th>91</th>\n",
1348 " <td>Air</td>\n",
1349 " <td>10562.0</td>\n",
1350 " <td>3.43512</td>\n",
1351 " <td>681.356</td>\n",
1352 " <td>0.000684</td>\n",
1353 " <td>0</td>\n",
1354 " <td>3.4</td>\n",
1355 " <td>1.0</td>\n",
1356 " <td>0.000011</td>\n",
1357 " <td>0.000000</td>\n",
1358 " <td>0.000000</td>\n",
1359 " <td>79.707049</td>\n",
1360 " <td>102.129933</td>\n",
1361 " <td>0</td>\n",
1362 " <td>9.140389</td>\n",
1363 " <td>True</td>\n",
1364 " <td>3.0</td>\n",
1365 " </tr>\n",
1366 " <tr>\n",
1367 " <th>92</th>\n",
1368 " <td>StainlessSteel</td>\n",
1369 " <td>10563.0</td>\n",
1370 " <td>3.43512</td>\n",
1371 " <td>681.420</td>\n",
1372 " <td>1.410940</td>\n",
1373 " <td>0</td>\n",
1374 " <td>1.0</td>\n",
1375 " <td>1.0</td>\n",
1376 " <td>0.057619</td>\n",
1377 " <td>1.165737</td>\n",
1378 " <td>1.165737</td>\n",
1379 " <td>80.872786</td>\n",
1380 " <td>103.540873</td>\n",
1381 " <td>0</td>\n",
1382 " <td>9.198008</td>\n",
1383 " <td>True</td>\n",
1384 " <td>3.0</td>\n",
1385 " </tr>\n",
1386 " <tr>\n",
1387 " <th>93</th>\n",
1388 " <td>Lead</td>\n",
1389 " <td>10584.0</td>\n",
1390 " <td>3.43512</td>\n",
1391 " <td>682.775</td>\n",
1392 " <td>30.668000</td>\n",
1393 " <td>0</td>\n",
1394 " <td>21.0</td>\n",
1395 " <td>1.0</td>\n",
1396 " <td>3.742115</td>\n",
1397 " <td>1.273470</td>\n",
1398 " <td>26.742870</td>\n",
1399 " <td>107.615656</td>\n",
1400 " <td>134.208873</td>\n",
1401 " <td>0</td>\n",
1402 " <td>12.940123</td>\n",
1403 " <td>True</td>\n",
1404 " <td>3.0</td>\n",
1405 " </tr>\n",
1406 " <tr>\n",
1407 " <th>94</th>\n",
1408 " <td>StainlessSteel</td>\n",
1409 " <td>10585.0</td>\n",
1410 " <td>3.43512</td>\n",
1411 " <td>682.839</td>\n",
1412 " <td>1.410080</td>\n",
1413 " <td>0</td>\n",
1414 " <td>1.0</td>\n",
1415 " <td>1.0</td>\n",
1416 " <td>0.057619</td>\n",
1417 " <td>1.165737</td>\n",
1418 " <td>1.165737</td>\n",
1419 " <td>108.781393</td>\n",
1420 " <td>135.618953</td>\n",
1421 " <td>0</td>\n",
1422 " <td>12.997741</td>\n",
1423 " <td>True</td>\n",
1424 " <td>3.0</td>\n",
1425 " </tr>\n",
1426 " <tr>\n",
1427 " <th>95</th>\n",
1428 " <td>Air</td>\n",
1429 " <td>10588.3</td>\n",
1430 " <td>3.43512</td>\n",
1431 " <td>683.055</td>\n",
1432 " <td>0.000683</td>\n",
1433 " <td>0</td>\n",
1434 " <td>3.3</td>\n",
1435 " <td>1.0</td>\n",
1436 " <td>0.000011</td>\n",
1437 " <td>0.000000</td>\n",
1438 " <td>0.000000</td>\n",
1439 " <td>108.781393</td>\n",
1440 " <td>135.619636</td>\n",
1441 " <td>0</td>\n",
1442 " <td>12.997752</td>\n",
1443 " <td>True</td>\n",
1444 " <td>3.0</td>\n",
1445 " </tr>\n",
1446 " <tr>\n",
1447 " <th>96</th>\n",
1448 " <td>HGC_G10-FR4</td>\n",
1449 " <td>10589.9</td>\n",
1450 " <td>3.43512</td>\n",
1451 " <td>683.159</td>\n",
1452 " <td>0.522460</td>\n",
1453 " <td>0</td>\n",
1454 " <td>1.6</td>\n",
1455 " <td>1.0</td>\n",
1456 " <td>0.009140</td>\n",
1457 " <td>0.320491</td>\n",
1458 " <td>0.512786</td>\n",
1459 " <td>109.294179</td>\n",
1460 " <td>136.142096</td>\n",
1461 " <td>0</td>\n",
1462 " <td>13.006892</td>\n",
1463 " <td>True</td>\n",
1464 " <td>3.0</td>\n",
1465 " </tr>\n",
1466 " <tr>\n",
1467 " <th>97</th>\n",
1468 " <td>Silicon</td>\n",
1469 " <td>10590.1</td>\n",
1470 " <td>3.43512</td>\n",
1471 " <td>683.166</td>\n",
1472 " <td>0.142680</td>\n",
1473 " <td>0</td>\n",
1474 " <td>0.2</td>\n",
1475 " <td>1.0</td>\n",
1476 " <td>0.002135</td>\n",
1477 " <td>0.387546</td>\n",
1478 " <td>0.077509</td>\n",
1479 " <td>109.371689</td>\n",
1480 " <td>136.284776</td>\n",
1481 " <td>0</td>\n",
1482 " <td>13.009027</td>\n",
1483 " <td>False</td>\n",
1484 " <td>3.0</td>\n",
1485 " </tr>\n",
1486 " <tr>\n",
1487 " <th>98</th>\n",
1488 " <td>Silicon</td>\n",
1489 " <td>10590.2</td>\n",
1490 " <td>3.43512</td>\n",
1491 " <td>683.178</td>\n",
1492 " <td>0.214020</td>\n",
1493 " <td>0</td>\n",
1494 " <td>0.1</td>\n",
1495 " <td>1.0</td>\n",
1496 " <td>0.001068</td>\n",
1497 " <td>0.387546</td>\n",
1498 " <td>0.038755</td>\n",
1499 " <td>109.410443</td>\n",
1500 " <td>136.498796</td>\n",
1501 " <td>0</td>\n",
1502 " <td>13.010095</td>\n",
1503 " <td>False</td>\n",
1504 " <td>3.0</td>\n",
1505 " </tr>\n",
1506 " <tr>\n",
1507 " <th>99</th>\n",
1508 " <td>WCu</td>\n",
1509 " <td>10591.6</td>\n",
1510 " <td>3.43512</td>\n",
1511 " <td>683.268</td>\n",
1512 " <td>2.584590</td>\n",
1513 " <td>0</td>\n",
1514 " <td>1.4</td>\n",
1515 " <td>1.0</td>\n",
1516 " <td>0.273304</td>\n",
1517 " <td>1.811710</td>\n",
1518 " <td>2.536394</td>\n",
1519 " <td>111.946837</td>\n",
1520 " <td>139.083386</td>\n",
1521 " <td>0</td>\n",
1522 " <td>13.283399</td>\n",
1523 " <td>True</td>\n",
1524 " <td>4.0</td>\n",
1525 " </tr>\n",
1526 " <tr>\n",
1527 " <th>100</th>\n",
1528 " <td>Copper</td>\n",
1529 " <td>10597.6</td>\n",
1530 " <td>3.43512</td>\n",
1531 " <td>683.655</td>\n",
1532 " <td>8.890630</td>\n",
1533 " <td>0</td>\n",
1534 " <td>6.0</td>\n",
1535 " <td>1.0</td>\n",
1536 " <td>0.417947</td>\n",
1537 " <td>1.257088</td>\n",
1538 " <td>7.542528</td>\n",
1539 " <td>119.489365</td>\n",
1540 " <td>147.974016</td>\n",
1541 " <td>0</td>\n",
1542 " <td>13.701345</td>\n",
1543 " <td>True</td>\n",
1544 " <td>4.0</td>\n",
1545 " </tr>\n",
1546 " <tr>\n",
1547 " <th>101</th>\n",
1548 " <td>Air</td>\n",
1549 " <td>10601.0</td>\n",
1550 " <td>3.43512</td>\n",
1551 " <td>683.871</td>\n",
1552 " <td>0.000683</td>\n",
1553 " <td>0</td>\n",
1554 " <td>3.4</td>\n",
1555 " <td>1.0</td>\n",
1556 " <td>0.000011</td>\n",
1557 " <td>0.000000</td>\n",
1558 " <td>0.000000</td>\n",
1559 " <td>119.489365</td>\n",
1560 " <td>147.974699</td>\n",
1561 " <td>0</td>\n",
1562 " <td>13.701356</td>\n",
1563 " <td>True</td>\n",
1564 " <td>4.0</td>\n",
1565 " </tr>\n",
1566 " <tr>\n",
1567 " <th>102</th>\n",
1568 " <td>StainlessSteel</td>\n",
1569 " <td>10602.0</td>\n",
1570 " <td>3.43512</td>\n",
1571 " <td>683.936</td>\n",
1572 " <td>1.409420</td>\n",
1573 " <td>0</td>\n",
1574 " <td>1.0</td>\n",
1575 " <td>1.0</td>\n",
1576 " <td>0.057619</td>\n",
1577 " <td>1.165737</td>\n",
1578 " <td>1.165737</td>\n",
1579 " <td>120.655102</td>\n",
1580 " <td>149.384119</td>\n",
1581 " <td>0</td>\n",
1582 " <td>13.758975</td>\n",
1583 " <td>True</td>\n",
1584 " <td>4.0</td>\n",
1585 " </tr>\n",
1586 " <tr>\n",
1587 " <th>103</th>\n",
1588 " <td>Lead</td>\n",
1589 " <td>10623.0</td>\n",
1590 " <td>3.43512</td>\n",
1591 " <td>685.291</td>\n",
1592 " <td>30.637800</td>\n",
1593 " <td>0</td>\n",
1594 " <td>21.0</td>\n",
1595 " <td>1.0</td>\n",
1596 " <td>3.742115</td>\n",
1597 " <td>1.273470</td>\n",
1598 " <td>26.742870</td>\n",
1599 " <td>147.397972</td>\n",
1600 " <td>180.021919</td>\n",
1601 " <td>0</td>\n",
1602 " <td>17.501090</td>\n",
1603 " <td>True</td>\n",
1604 " <td>4.0</td>\n",
1605 " </tr>\n",
1606 " <tr>\n",
1607 " <th>104</th>\n",
1608 " <td>StainlessSteel</td>\n",
1609 " <td>10624.0</td>\n",
1610 " <td>3.43512</td>\n",
1611 " <td>685.355</td>\n",
1612 " <td>1.408570</td>\n",
1613 " <td>0</td>\n",
1614 " <td>1.0</td>\n",
1615 " <td>1.0</td>\n",
1616 " <td>0.057619</td>\n",
1617 " <td>1.165737</td>\n",
1618 " <td>1.165737</td>\n",
1619 " <td>148.563710</td>\n",
1620 " <td>181.430489</td>\n",
1621 " <td>0</td>\n",
1622 " <td>17.558709</td>\n",
1623 " <td>True</td>\n",
1624 " <td>4.0</td>\n",
1625 " </tr>\n",
1626 " <tr>\n",
1627 " <th>105</th>\n",
1628 " <td>Air</td>\n",
1629 " <td>10627.3</td>\n",
1630 " <td>3.43512</td>\n",
1631 " <td>685.571</td>\n",
1632 " <td>0.000682</td>\n",
1633 " <td>0</td>\n",
1634 " <td>3.3</td>\n",
1635 " <td>1.0</td>\n",
1636 " <td>0.000011</td>\n",
1637 " <td>0.000000</td>\n",
1638 " <td>0.000000</td>\n",
1639 " <td>148.563710</td>\n",
1640 " <td>181.431172</td>\n",
1641 " <td>0</td>\n",
1642 " <td>17.558720</td>\n",
1643 " <td>True</td>\n",
1644 " <td>4.0</td>\n",
1645 " </tr>\n",
1646 " <tr>\n",
1647 " <th>106</th>\n",
1648 " <td>Copper</td>\n",
1649 " <td>10633.3</td>\n",
1650 " <td>3.43512</td>\n",
1651 " <td>685.958</td>\n",
1652 " <td>8.880230</td>\n",
1653 " <td>0</td>\n",
1654 " <td>6.0</td>\n",
1655 " <td>1.0</td>\n",
1656 " <td>0.417947</td>\n",
1657 " <td>1.257088</td>\n",
1658 " <td>7.542528</td>\n",
1659 " <td>156.106238</td>\n",
1660 " <td>190.311402</td>\n",
1661 " <td>0</td>\n",
1662 " <td>17.976666</td>\n",
1663 " <td>True</td>\n",
1664 " <td>4.0</td>\n",
1665 " </tr>\n",
1666 " <tr>\n",
1667 " <th>107</th>\n",
1668 " <td>WCu</td>\n",
1669 " <td>10634.8</td>\n",
1670 " <td>3.43512</td>\n",
1671 " <td>686.049</td>\n",
1672 " <td>2.581140</td>\n",
1673 " <td>0</td>\n",
1674 " <td>1.5</td>\n",
1675 " <td>1.0</td>\n",
1676 " <td>0.292826</td>\n",
1677 " <td>1.811710</td>\n",
1678 " <td>2.717565</td>\n",
1679 " <td>158.823803</td>\n",
1680 " <td>192.892542</td>\n",
1681 " <td>0</td>\n",
1682 " <td>18.269492</td>\n",
1683 " <td>True</td>\n",
1684 " <td>4.0</td>\n",
1685 " </tr>\n",
1686 " <tr>\n",
1687 " <th>108</th>\n",
1688 " <td>Silicon</td>\n",
1689 " <td>10634.9</td>\n",
1690 " <td>3.43512</td>\n",
1691 " <td>686.060</td>\n",
1692 " <td>0.214078</td>\n",
1693 " <td>0</td>\n",
1694 " <td>0.1</td>\n",
1695 " <td>1.0</td>\n",
1696 " <td>0.001068</td>\n",
1697 " <td>0.387546</td>\n",
1698 " <td>0.038755</td>\n",
1699 " <td>158.862557</td>\n",
1700 " <td>193.106620</td>\n",
1701 " <td>0</td>\n",
1702 " <td>18.270559</td>\n",
1703 " <td>False</td>\n",
1704 " <td>4.0</td>\n",
1705 " </tr>\n",
1706 " <tr>\n",
1707 " <th>109</th>\n",
1708 " <td>Silicon</td>\n",
1709 " <td>10635.0</td>\n",
1710 " <td>3.43512</td>\n",
1711 " <td>686.068</td>\n",
1712 " <td>0.142719</td>\n",
1713 " <td>0</td>\n",
1714 " <td>0.1</td>\n",
1715 " <td>1.0</td>\n",
1716 " <td>0.001068</td>\n",
1717 " <td>0.387546</td>\n",
1718 " <td>0.038755</td>\n",
1719 " <td>158.901312</td>\n",
1720 " <td>193.249339</td>\n",
1721 " <td>0</td>\n",
1722 " <td>18.271627</td>\n",
1723 " <td>False</td>\n",
1724 " <td>4.0</td>\n",
1725 " </tr>\n",
1726 " <tr>\n",
1727 " <th>110</th>\n",
1728 " <td>HGC_G10-FR4</td>\n",
1729 " <td>10636.6</td>\n",
1730 " <td>3.43512</td>\n",
1731 " <td>686.171</td>\n",
1732 " <td>0.521664</td>\n",
1733 " <td>0</td>\n",
1734 " <td>1.6</td>\n",
1735 " <td>1.0</td>\n",
1736 " <td>0.009140</td>\n",
1737 " <td>0.320491</td>\n",
1738 " <td>0.512786</td>\n",
1739 " <td>159.414098</td>\n",
1740 " <td>193.771003</td>\n",
1741 " <td>0</td>\n",
1742 " <td>18.280767</td>\n",
1743 " <td>True</td>\n",
1744 " <td>5.0</td>\n",
1745 " </tr>\n",
1746 " <tr>\n",
1747 " <th>111</th>\n",
1748 " <td>Air</td>\n",
1749 " <td>10640.0</td>\n",
1750 " <td>3.43512</td>\n",
1751 " <td>686.387</td>\n",
1752 " <td>0.000682</td>\n",
1753 " <td>0</td>\n",
1754 " <td>3.4</td>\n",
1755 " <td>1.0</td>\n",
1756 " <td>0.000011</td>\n",
1757 " <td>0.000000</td>\n",
1758 " <td>0.000000</td>\n",
1759 " <td>159.414098</td>\n",
1760 " <td>193.771685</td>\n",
1761 " <td>0</td>\n",
1762 " <td>18.280778</td>\n",
1763 " <td>True</td>\n",
1764 " <td>5.0</td>\n",
1765 " </tr>\n",
1766 " <tr>\n",
1767 " <th>112</th>\n",
1768 " <td>StainlessSteel</td>\n",
1769 " <td>10641.0</td>\n",
1770 " <td>3.43512</td>\n",
1771 " <td>686.452</td>\n",
1772 " <td>1.407910</td>\n",
1773 " <td>0</td>\n",
1774 " <td>1.0</td>\n",
1775 " <td>1.0</td>\n",
1776 " <td>0.057619</td>\n",
1777 " <td>1.165737</td>\n",
1778 " <td>1.165737</td>\n",
1779 " <td>160.579835</td>\n",
1780 " <td>195.179595</td>\n",
1781 " <td>0</td>\n",
1782 " <td>18.338397</td>\n",
1783 " <td>True</td>\n",
1784 " <td>5.0</td>\n",
1785 " </tr>\n",
1786 " <tr>\n",
1787 " <th>113</th>\n",
1788 " <td>Lead</td>\n",
1789 " <td>10662.0</td>\n",
1790 " <td>3.43512</td>\n",
1791 " <td>687.807</td>\n",
1792 " <td>30.607900</td>\n",
1793 " <td>0</td>\n",
1794 " <td>21.0</td>\n",
1795 " <td>1.0</td>\n",
1796 " <td>3.742115</td>\n",
1797 " <td>1.273470</td>\n",
1798 " <td>26.742870</td>\n",
1799 " <td>187.322705</td>\n",
1800 " <td>225.787495</td>\n",
1801 " <td>0</td>\n",
1802 " <td>22.080512</td>\n",
1803 " <td>True</td>\n",
1804 " <td>5.0</td>\n",
1805 " </tr>\n",
1806 " <tr>\n",
1807 " <th>114</th>\n",
1808 " <td>StainlessSteel</td>\n",
1809 " <td>10663.0</td>\n",
1810 " <td>3.43512</td>\n",
1811 " <td>687.871</td>\n",
1812 " <td>1.407070</td>\n",
1813 " <td>0</td>\n",
1814 " <td>1.0</td>\n",
1815 " <td>1.0</td>\n",
1816 " <td>0.057619</td>\n",
1817 " <td>1.165737</td>\n",
1818 " <td>1.165737</td>\n",
1819 " <td>188.488442</td>\n",
1820 " <td>227.194565</td>\n",
1821 " <td>0</td>\n",
1822 " <td>22.138130</td>\n",
1823 " <td>True</td>\n",
1824 " <td>5.0</td>\n",
1825 " </tr>\n",
1826 " <tr>\n",
1827 " <th>115</th>\n",
1828 " <td>Air</td>\n",
1829 " <td>10666.3</td>\n",
1830 " <td>3.43512</td>\n",
1831 " <td>688.087</td>\n",
1832 " <td>0.000682</td>\n",
1833 " <td>0</td>\n",
1834 " <td>3.3</td>\n",
1835 " <td>1.0</td>\n",
1836 " <td>0.000011</td>\n",
1837 " <td>0.000000</td>\n",
1838 " <td>0.000000</td>\n",
1839 " <td>188.488442</td>\n",
1840 " <td>227.195247</td>\n",
1841 " <td>0</td>\n",
1842 " <td>22.138141</td>\n",
1843 " <td>True</td>\n",
1844 " <td>5.0</td>\n",
1845 " </tr>\n",
1846 " <tr>\n",
1847 " <th>116</th>\n",
1848 " <td>HGC_G10-FR4</td>\n",
1849 " <td>10667.9</td>\n",
1850 " <td>3.43512</td>\n",
1851 " <td>688.190</td>\n",
1852 " <td>0.521134</td>\n",
1853 " <td>0</td>\n",
1854 " <td>1.6</td>\n",
1855 " <td>1.0</td>\n",
1856 " <td>0.009140</td>\n",
1857 " <td>0.320491</td>\n",
1858 " <td>0.512786</td>\n",
1859 " <td>189.001228</td>\n",
1860 " <td>227.716381</td>\n",
1861 " <td>0</td>\n",
1862 " <td>22.147281</td>\n",
1863 " <td>True</td>\n",
1864 " <td>5.0</td>\n",
1865 " </tr>\n",
1866 " <tr>\n",
1867 " <th>117</th>\n",
1868 " <td>Silicon</td>\n",
1869 " <td>10668.1</td>\n",
1870 " <td>3.43512</td>\n",
1871 " <td>688.198</td>\n",
1872 " <td>0.142755</td>\n",
1873 " <td>0</td>\n",
1874 " <td>0.2</td>\n",
1875 " <td>1.0</td>\n",
1876 " <td>0.002135</td>\n",
1877 " <td>0.387546</td>\n",
1878 " <td>0.077509</td>\n",
1879 " <td>189.078738</td>\n",
1880 " <td>227.859136</td>\n",
1881 " <td>0</td>\n",
1882 " <td>22.149416</td>\n",
1883 " <td>False</td>\n",
1884 " <td>5.0</td>\n",
1885 " </tr>\n",
1886 " <tr>\n",
1887 " <th>118</th>\n",
1888 " <td>Silicon</td>\n",
1889 " <td>10668.2</td>\n",
1890 " <td>3.43512</td>\n",
1891 " <td>688.210</td>\n",
1892 " <td>0.214132</td>\n",
1893 " <td>0</td>\n",
1894 " <td>0.1</td>\n",
1895 " <td>1.0</td>\n",
1896 " <td>0.001068</td>\n",
1897 " <td>0.387546</td>\n",
1898 " <td>0.038755</td>\n",
1899 " <td>189.117492</td>\n",
1900 " <td>228.073268</td>\n",
1901 " <td>0</td>\n",
1902 " <td>22.150484</td>\n",
1903 " <td>False</td>\n",
1904 " <td>5.0</td>\n",
1905 " </tr>\n",
1906 " <tr>\n",
1907 " <th>119</th>\n",
1908 " <td>WCu</td>\n",
1909 " <td>10669.6</td>\n",
1910 " <td>3.43512</td>\n",
1911 " <td>688.300</td>\n",
1912 " <td>2.578370</td>\n",
1913 " <td>0</td>\n",
1914 " <td>1.4</td>\n",
1915 " <td>1.0</td>\n",
1916 " <td>0.273304</td>\n",
1917 " <td>1.811710</td>\n",
1918 " <td>2.536394</td>\n",
1919 " <td>191.653886</td>\n",
1920 " <td>230.651638</td>\n",
1921 " <td>0</td>\n",
1922 " <td>22.423788</td>\n",
1923 " <td>True</td>\n",
1924 " <td>6.0</td>\n",
1925 " </tr>\n",
1926 " <tr>\n",
1927 " <th>120</th>\n",
1928 " <td>Copper</td>\n",
1929 " <td>10675.6</td>\n",
1930 " <td>3.43512</td>\n",
1931 " <td>688.687</td>\n",
1932 " <td>8.868010</td>\n",
1933 " <td>0</td>\n",
1934 " <td>6.0</td>\n",
1935 " <td>1.0</td>\n",
1936 " <td>0.417947</td>\n",
1937 " <td>1.257088</td>\n",
1938 " <td>7.542528</td>\n",
1939 " <td>199.196414</td>\n",
1940 " <td>239.519648</td>\n",
1941 " <td>0</td>\n",
1942 " <td>22.841734</td>\n",
1943 " <td>True</td>\n",
1944 " <td>6.0</td>\n",
1945 " </tr>\n",
1946 " <tr>\n",
1947 " <th>121</th>\n",
1948 " <td>Air</td>\n",
1949 " <td>10679.0</td>\n",
1950 " <td>3.43512</td>\n",
1951 " <td>688.903</td>\n",
1952 " <td>0.000681</td>\n",
1953 " <td>0</td>\n",
1954 " <td>3.4</td>\n",
1955 " <td>1.0</td>\n",
1956 " <td>0.000011</td>\n",
1957 " <td>0.000000</td>\n",
1958 " <td>0.000000</td>\n",
1959 " <td>199.196414</td>\n",
1960 " <td>239.520329</td>\n",
1961 " <td>0</td>\n",
1962 " <td>22.841746</td>\n",
1963 " <td>True</td>\n",
1964 " <td>6.0</td>\n",
1965 " </tr>\n",
1966 " <tr>\n",
1967 " <th>122</th>\n",
1968 " <td>StainlessSteel</td>\n",
1969 " <td>10680.0</td>\n",
1970 " <td>3.43512</td>\n",
1971 " <td>688.968</td>\n",
1972 " <td>1.406410</td>\n",
1973 " <td>0</td>\n",
1974 " <td>1.0</td>\n",
1975 " <td>1.0</td>\n",
1976 " <td>0.057619</td>\n",
1977 " <td>1.165737</td>\n",
1978 " <td>1.165737</td>\n",
1979 " <td>200.362151</td>\n",
1980 " <td>240.926739</td>\n",
1981 " <td>0</td>\n",
1982 " <td>22.899364</td>\n",
1983 " <td>True</td>\n",
1984 " <td>6.0</td>\n",
1985 " </tr>\n",
1986 " <tr>\n",
1987 " <th>123</th>\n",
1988 " <td>Lead</td>\n",
1989 " <td>10701.0</td>\n",
1990 " <td>3.43512</td>\n",
1991 " <td>690.322</td>\n",
1992 " <td>30.578200</td>\n",
1993 " <td>0</td>\n",
1994 " <td>21.0</td>\n",
1995 " <td>1.0</td>\n",
1996 " <td>3.742115</td>\n",
1997 " <td>1.273470</td>\n",
1998 " <td>26.742870</td>\n",
1999 " <td>227.105021</td>\n",
2000 " <td>271.504939</td>\n",
2001 " <td>0</td>\n",
2002 " <td>26.641479</td>\n",
2003 " <td>True</td>\n",
2004 " <td>6.0</td>\n",
2005 " </tr>\n",
2006 " <tr>\n",
2007 " <th>124</th>\n",
2008 " <td>StainlessSteel</td>\n",
2009 " <td>10702.0</td>\n",
2010 " <td>3.43512</td>\n",
2011 " <td>690.387</td>\n",
2012 " <td>1.405580</td>\n",
2013 " <td>0</td>\n",
2014 " <td>1.0</td>\n",
2015 " <td>1.0</td>\n",
2016 " <td>0.057619</td>\n",
2017 " <td>1.165737</td>\n",
2018 " <td>1.165737</td>\n",
2019 " <td>228.270759</td>\n",
2020 " <td>272.910519</td>\n",
2021 " <td>0</td>\n",
2022 " <td>26.699098</td>\n",
2023 " <td>True</td>\n",
2024 " <td>6.0</td>\n",
2025 " </tr>\n",
2026 " <tr>\n",
2027 " <th>125</th>\n",
2028 " <td>Air</td>\n",
2029 " <td>10705.3</td>\n",
2030 " <td>3.43512</td>\n",
2031 " <td>690.603</td>\n",
2032 " <td>0.000681</td>\n",
2033 " <td>0</td>\n",
2034 " <td>3.3</td>\n",
2035 " <td>1.0</td>\n",
2036 " <td>0.000011</td>\n",
2037 " <td>0.000000</td>\n",
2038 " <td>0.000000</td>\n",
2039 " <td>228.270759</td>\n",
2040 " <td>272.911200</td>\n",
2041 " <td>0</td>\n",
2042 " <td>26.699109</td>\n",
2043 " <td>True</td>\n",
2044 " <td>6.0</td>\n",
2045 " </tr>\n",
2046 " <tr>\n",
2047 " <th>126</th>\n",
2048 " <td>Copper</td>\n",
2049 " <td>10711.3</td>\n",
2050 " <td>3.43512</td>\n",
2051 " <td>690.990</td>\n",
2052 " <td>8.857770</td>\n",
2053 " <td>0</td>\n",
2054 " <td>6.0</td>\n",
2055 " <td>1.0</td>\n",
2056 " <td>0.417947</td>\n",
2057 " <td>1.257088</td>\n",
2058 " <td>7.542528</td>\n",
2059 " <td>235.813287</td>\n",
2060 " <td>281.768970</td>\n",
2061 " <td>0</td>\n",
2062 " <td>27.117055</td>\n",
2063 " <td>True</td>\n",
2064 " <td>6.0</td>\n",
2065 " </tr>\n",
2066 " <tr>\n",
2067 " <th>127</th>\n",
2068 " <td>WCu</td>\n",
2069 " <td>10712.8</td>\n",
2070 " <td>3.43512</td>\n",
2071 " <td>691.080</td>\n",
2072 " <td>2.574970</td>\n",
2073 " <td>0</td>\n",
2074 " <td>1.5</td>\n",
2075 " <td>1.0</td>\n",
2076 " <td>0.292826</td>\n",
2077 " <td>1.811710</td>\n",
2078 " <td>2.717565</td>\n",
2079 " <td>238.530852</td>\n",
2080 " <td>284.343940</td>\n",
2081 " <td>0</td>\n",
2082 " <td>27.409881</td>\n",
2083 " <td>True</td>\n",
2084 " <td>6.0</td>\n",
2085 " </tr>\n",
2086 " <tr>\n",
2087 " <th>128</th>\n",
2088 " <td>Silicon</td>\n",
2089 " <td>10712.9</td>\n",
2090 " <td>3.43512</td>\n",
2091 " <td>691.092</td>\n",
2092 " <td>0.214095</td>\n",
2093 " <td>0</td>\n",
2094 " <td>0.1</td>\n",
2095 " <td>1.0</td>\n",
2096 " <td>0.001068</td>\n",
2097 " <td>0.387546</td>\n",
2098 " <td>0.038755</td>\n",
2099 " <td>238.569606</td>\n",
2100 " <td>284.558035</td>\n",
2101 " <td>0</td>\n",
2102 " <td>27.410948</td>\n",
2103 " <td>False</td>\n",
2104 " <td>6.0</td>\n",
2105 " </tr>\n",
2106 " <tr>\n",
2107 " <th>129</th>\n",
2108 " <td>Silicon</td>\n",
2109 " <td>10713.0</td>\n",
2110 " <td>3.43512</td>\n",
2111 " <td>691.100</td>\n",
2112 " <td>0.142730</td>\n",
2113 " <td>0</td>\n",
2114 " <td>0.1</td>\n",
2115 " <td>1.0</td>\n",
2116 " <td>0.001068</td>\n",
2117 " <td>0.387546</td>\n",
2118 " <td>0.038755</td>\n",
2119 " <td>238.608361</td>\n",
2120 " <td>284.700765</td>\n",
2121 " <td>0</td>\n",
2122 " <td>27.412016</td>\n",
2123 " <td>False</td>\n",
2124 " <td>6.0</td>\n",
2125 " </tr>\n",
2126 " <tr>\n",
2127 " <th>130</th>\n",
2128 " <td>HGC_G10-FR4</td>\n",
2129 " <td>10714.6</td>\n",
2130 " <td>3.43512</td>\n",
2131 " <td>691.203</td>\n",
2132 " <td>0.520350</td>\n",
2133 " <td>0</td>\n",
2134 " <td>1.6</td>\n",
2135 " <td>1.0</td>\n",
2136 " <td>0.009140</td>\n",
2137 " <td>0.320491</td>\n",
2138 " <td>0.512786</td>\n",
2139 " <td>239.121147</td>\n",
2140 " <td>285.221115</td>\n",
2141 " <td>0</td>\n",
2142 " <td>27.421156</td>\n",
2143 " <td>True</td>\n",
2144 " <td>7.0</td>\n",
2145 " </tr>\n",
2146 " <tr>\n",
2147 " <th>131</th>\n",
2148 " <td>Air</td>\n",
2149 " <td>10718.0</td>\n",
2150 " <td>3.43512</td>\n",
2151 " <td>691.419</td>\n",
2152 " <td>0.000680</td>\n",
2153 " <td>0</td>\n",
2154 " <td>3.4</td>\n",
2155 " <td>1.0</td>\n",
2156 " <td>0.000011</td>\n",
2157 " <td>0.000000</td>\n",
2158 " <td>0.000000</td>\n",
2159 " <td>239.121147</td>\n",
2160 " <td>285.221795</td>\n",
2161 " <td>0</td>\n",
2162 " <td>27.421167</td>\n",
2163 " <td>True</td>\n",
2164 " <td>7.0</td>\n",
2165 " </tr>\n",
2166 " <tr>\n",
2167 " <th>132</th>\n",
2168 " <td>StainlessSteel</td>\n",
2169 " <td>10788.0</td>\n",
2170 " <td>3.43512</td>\n",
2171 " <td>695.935</td>\n",
2172 " <td>98.308600</td>\n",
2173 " <td>0</td>\n",
2174 " <td>70.0</td>\n",
2175 " <td>1.0</td>\n",
2176 " <td>4.033304</td>\n",
2177 " <td>1.165737</td>\n",
2178 " <td>81.601601</td>\n",
2179 " <td>320.722748</td>\n",
2180 " <td>383.530395</td>\n",
2181 " <td>0</td>\n",
2182 " <td>31.454471</td>\n",
2183 " <td>True</td>\n",
2184 " <td>7.0</td>\n",
2185 " </tr>\n",
2186 " <tr>\n",
2187 " <th>133</th>\n",
2188 " <td>Air</td>\n",
2189 " <td>10791.3</td>\n",
2190 " <td>3.43512</td>\n",
2191 " <td>696.151</td>\n",
2192 " <td>0.000679</td>\n",
2193 " <td>0</td>\n",
2194 " <td>3.3</td>\n",
2195 " <td>1.0</td>\n",
2196 " <td>0.000011</td>\n",
2197 " <td>0.000000</td>\n",
2198 " <td>0.000000</td>\n",
2199 " <td>320.722748</td>\n",
2200 " <td>383.531074</td>\n",
2201 " <td>0</td>\n",
2202 " <td>31.454482</td>\n",
2203 " <td>True</td>\n",
2204 " <td>7.0</td>\n",
2205 " </tr>\n",
2206 " <tr>\n",
2207 " <th>134</th>\n",
2208 " <td>HGC_G10-FR4</td>\n",
2209 " <td>10792.9</td>\n",
2210 " <td>3.43512</td>\n",
2211 " <td>696.254</td>\n",
2212 " <td>0.519051</td>\n",
2213 " <td>0</td>\n",
2214 " <td>1.6</td>\n",
2215 " <td>1.0</td>\n",
2216 " <td>0.009140</td>\n",
2217 " <td>0.320491</td>\n",
2218 " <td>0.512786</td>\n",
2219 " <td>321.235534</td>\n",
2220 " <td>384.050125</td>\n",
2221 " <td>0</td>\n",
2222 " <td>31.463622</td>\n",
2223 " <td>True</td>\n",
2224 " <td>7.0</td>\n",
2225 " </tr>\n",
2226 " <tr>\n",
2227 " <th>135</th>\n",
2228 " <td>Silicon</td>\n",
2229 " <td>10793.1</td>\n",
2230 " <td>3.43512</td>\n",
2231 " <td>696.262</td>\n",
2232 " <td>0.142696</td>\n",
2233 " <td>0</td>\n",
2234 " <td>0.2</td>\n",
2235 " <td>1.0</td>\n",
2236 " <td>0.002135</td>\n",
2237 " <td>0.387546</td>\n",
2238 " <td>0.077509</td>\n",
2239 " <td>321.313043</td>\n",
2240 " <td>384.192821</td>\n",
2241 " <td>0</td>\n",
2242 " <td>31.465757</td>\n",
2243 " <td>False</td>\n",
2244 " <td>7.0</td>\n",
2245 " </tr>\n",
2246 " <tr>\n",
2247 " <th>136</th>\n",
2248 " <td>Silicon</td>\n",
2249 " <td>10793.2</td>\n",
2250 " <td>3.43512</td>\n",
2251 " <td>696.273</td>\n",
2252 " <td>0.214043</td>\n",
2253 " <td>0</td>\n",
2254 " <td>0.1</td>\n",
2255 " <td>1.0</td>\n",
2256 " <td>0.001068</td>\n",
2257 " <td>0.387546</td>\n",
2258 " <td>0.038755</td>\n",
2259 " <td>321.351798</td>\n",
2260 " <td>384.406864</td>\n",
2261 " <td>0</td>\n",
2262 " <td>31.466824</td>\n",
2263 " <td>False</td>\n",
2264 " <td>7.0</td>\n",
2265 " </tr>\n",
2266 " <tr>\n",
2267 " <th>137</th>\n",
2268 " <td>WCu</td>\n",
2269 " <td>10794.6</td>\n",
2270 " <td>3.43512</td>\n",
2271 " <td>696.364</td>\n",
2272 " <td>2.568580</td>\n",
2273 " <td>0</td>\n",
2274 " <td>1.4</td>\n",
2275 " <td>1.0</td>\n",
2276 " <td>0.273304</td>\n",
2277 " <td>1.811710</td>\n",
2278 " <td>2.536394</td>\n",
2279 " <td>323.888192</td>\n",
2280 " <td>386.975444</td>\n",
2281 " <td>0</td>\n",
2282 " <td>31.740128</td>\n",
2283 " <td>True</td>\n",
2284 " <td>8.0</td>\n",
2285 " </tr>\n",
2286 " <tr>\n",
2287 " <th>138</th>\n",
2288 " <td>Copper</td>\n",
2289 " <td>10800.6</td>\n",
2290 " <td>3.43512</td>\n",
2291 " <td>696.751</td>\n",
2292 " <td>8.832450</td>\n",
2293 " <td>0</td>\n",
2294 " <td>6.0</td>\n",
2295 " <td>1.0</td>\n",
2296 " <td>0.417947</td>\n",
2297 " <td>1.257088</td>\n",
2298 " <td>7.542528</td>\n",
2299 " <td>331.430720</td>\n",
2300 " <td>395.807894</td>\n",
2301 " <td>0</td>\n",
2302 " <td>32.158075</td>\n",
2303 " <td>True</td>\n",
2304 " <td>8.0</td>\n",
2305 " </tr>\n",
2306 " <tr>\n",
2307 " <th>139</th>\n",
2308 " <td>Air</td>\n",
2309 " <td>10804.0</td>\n",
2310 " <td>3.43512</td>\n",
2311 " <td>696.967</td>\n",
2312 " <td>0.000679</td>\n",
2313 " <td>0</td>\n",
2314 " <td>3.4</td>\n",
2315 " <td>1.0</td>\n",
2316 " <td>0.000011</td>\n",
2317 " <td>0.000000</td>\n",
2318 " <td>0.000000</td>\n",
2319 " <td>331.430720</td>\n",
2320 " <td>395.808573</td>\n",
2321 " <td>0</td>\n",
2322 " <td>32.158086</td>\n",
2323 " <td>True</td>\n",
2324 " <td>8.0</td>\n",
2325 " </tr>\n",
2326 " <tr>\n",
2327 " <th>140</th>\n",
2328 " <td>StainlessSteel</td>\n",
2329 " <td>10874.0</td>\n",
2330 " <td>3.43512</td>\n",
2331 " <td>701.483</td>\n",
2332 " <td>98.083000</td>\n",
2333 " <td>0</td>\n",
2334 " <td>70.0</td>\n",
2335 " <td>1.0</td>\n",
2336 " <td>4.033304</td>\n",
2337 " <td>1.165737</td>\n",
2338 " <td>81.601601</td>\n",
2339 " <td>413.032321</td>\n",
2340 " <td>493.891573</td>\n",
2341 " <td>0</td>\n",
2342 " <td>36.191390</td>\n",
2343 " <td>True</td>\n",
2344 " <td>8.0</td>\n",
2345 " </tr>\n",
2346 " <tr>\n",
2347 " <th>141</th>\n",
2348 " <td>Air</td>\n",
2349 " <td>10877.3</td>\n",
2350 " <td>3.43512</td>\n",
2351 " <td>701.699</td>\n",
2352 " <td>0.000677</td>\n",
2353 " <td>0</td>\n",
2354 " <td>3.3</td>\n",
2355 " <td>1.0</td>\n",
2356 " <td>0.000011</td>\n",
2357 " <td>0.000000</td>\n",
2358 " <td>0.000000</td>\n",
2359 " <td>413.032321</td>\n",
2360 " <td>493.892250</td>\n",
2361 " <td>0</td>\n",
2362 " <td>36.191401</td>\n",
2363 " <td>True</td>\n",
2364 " <td>8.0</td>\n",
2365 " </tr>\n",
2366 " <tr>\n",
2367 " <th>142</th>\n",
2368 " <td>Copper</td>\n",
2369 " <td>10883.3</td>\n",
2370 " <td>3.43512</td>\n",
2371 " <td>702.086</td>\n",
2372 " <td>8.809390</td>\n",
2373 " <td>0</td>\n",
2374 " <td>6.0</td>\n",
2375 " <td>1.0</td>\n",
2376 " <td>0.417947</td>\n",
2377 " <td>1.257088</td>\n",
2378 " <td>7.542528</td>\n",
2379 " <td>420.574849</td>\n",
2380 " <td>502.701640</td>\n",
2381 " <td>0</td>\n",
2382 " <td>36.609347</td>\n",
2383 " <td>True</td>\n",
2384 " <td>8.0</td>\n",
2385 " </tr>\n",
2386 " <tr>\n",
2387 " <th>143</th>\n",
2388 " <td>WCu</td>\n",
2389 " <td>10884.8</td>\n",
2390 " <td>3.43512</td>\n",
2391 " <td>702.176</td>\n",
2392 " <td>2.561680</td>\n",
2393 " <td>0</td>\n",
2394 " <td>1.5</td>\n",
2395 " <td>1.0</td>\n",
2396 " <td>0.292826</td>\n",
2397 " <td>1.811710</td>\n",
2398 " <td>2.717565</td>\n",
2399 " <td>423.292414</td>\n",
2400 " <td>505.263320</td>\n",
2401 " <td>0</td>\n",
2402 " <td>36.902173</td>\n",
2403 " <td>True</td>\n",
2404 " <td>8.0</td>\n",
2405 " </tr>\n",
2406 " <tr>\n",
2407 " <th>144</th>\n",
2408 " <td>Silicon</td>\n",
2409 " <td>10884.9</td>\n",
2410 " <td>3.43512</td>\n",
2411 " <td>702.188</td>\n",
2412 " <td>0.214112</td>\n",
2413 " <td>0</td>\n",
2414 " <td>0.1</td>\n",
2415 " <td>1.0</td>\n",
2416 " <td>0.001068</td>\n",
2417 " <td>0.387546</td>\n",
2418 " <td>0.038755</td>\n",
2419 " <td>423.331168</td>\n",
2420 " <td>505.477432</td>\n",
2421 " <td>0</td>\n",
2422 " <td>36.903241</td>\n",
2423 " <td>False</td>\n",
2424 " <td>8.0</td>\n",
2425 " </tr>\n",
2426 " <tr>\n",
2427 " <th>145</th>\n",
2428 " <td>Silicon</td>\n",
2429 " <td>10885.0</td>\n",
2430 " <td>3.43512</td>\n",
2431 " <td>702.196</td>\n",
2432 " <td>0.142742</td>\n",
2433 " <td>0</td>\n",
2434 " <td>0.1</td>\n",
2435 " <td>1.0</td>\n",
2436 " <td>0.001068</td>\n",
2437 " <td>0.387546</td>\n",
2438 " <td>0.038755</td>\n",
2439 " <td>423.369923</td>\n",
2440 " <td>505.620174</td>\n",
2441 " <td>0</td>\n",
2442 " <td>36.904308</td>\n",
2443 " <td>False</td>\n",
2444 " <td>8.0</td>\n",
2445 " </tr>\n",
2446 " <tr>\n",
2447 " <th>146</th>\n",
2448 " <td>HGC_G10-FR4</td>\n",
2449 " <td>10886.6</td>\n",
2450 " <td>3.43512</td>\n",
2451 " <td>702.299</td>\n",
2452 " <td>0.517521</td>\n",
2453 " <td>0</td>\n",
2454 " <td>1.6</td>\n",
2455 " <td>1.0</td>\n",
2456 " <td>0.009140</td>\n",
2457 " <td>0.320491</td>\n",
2458 " <td>0.512786</td>\n",
2459 " <td>423.882709</td>\n",
2460 " <td>506.137695</td>\n",
2461 " <td>0</td>\n",
2462 " <td>36.913448</td>\n",
2463 " <td>True</td>\n",
2464 " <td>9.0</td>\n",
2465 " </tr>\n",
2466 " <tr>\n",
2467 " <th>147</th>\n",
2468 " <td>Air</td>\n",
2469 " <td>10890.0</td>\n",
2470 " <td>3.43512</td>\n",
2471 " <td>702.515</td>\n",
2472 " <td>0.000677</td>\n",
2473 " <td>0</td>\n",
2474 " <td>3.4</td>\n",
2475 " <td>1.0</td>\n",
2476 " <td>0.000011</td>\n",
2477 " <td>0.000000</td>\n",
2478 " <td>0.000000</td>\n",
2479 " <td>423.882709</td>\n",
2480 " <td>506.138372</td>\n",
2481 " <td>0</td>\n",
2482 " <td>36.913459</td>\n",
2483 " <td>True</td>\n",
2484 " <td>9.0</td>\n",
2485 " </tr>\n",
2486 " <tr>\n",
2487 " <th>148</th>\n",
2488 " <td>Polyethylene</td>\n",
2489 " <td>10484.0</td>\n",
2490 " <td>0.00000</td>\n",
2491 " <td>0.000</td>\n",
2492 " <td>0.000000</td>\n",
2493 " <td>0</td>\n",
2494 " <td>6.0</td>\n",
2495 " <td>1.0</td>\n",
2496 " <td>0.011926</td>\n",
2497 " <td>0.197027</td>\n",
2498 " <td>1.182160</td>\n",
2499 " <td>425.064869</td>\n",
2500 " <td>506.138372</td>\n",
2501 " <td>0</td>\n",
2502 " <td>36.925385</td>\n",
2503 " <td>True</td>\n",
2504 " <td>9.0</td>\n",
2505 " </tr>\n",
2506 " </tbody>\n",
2507 "</table>\n",
2508 "</div>"
2509 ],
2510 "text/plain": [
2511 " Mat Z Eta R Etable Efull \\\n",
2512 "71 Polyethylene 10484.0 3.43512 676.324 10.293100 0 \n",
2513 "72 StainlessSteel 10485.0 3.43512 676.388 1.414020 0 \n",
2514 "73 Lead 10506.0 3.43512 677.743 30.729100 0 \n",
2515 "74 StainlessSteel 10507.0 3.43512 677.807 1.413140 0 \n",
2516 "75 Air 10510.3 3.43512 678.024 0.000685 0 \n",
2517 "76 HGC_G10-FR4 10511.9 3.43512 678.127 0.523805 0 \n",
2518 "77 Silicon 10512.1 3.43512 678.135 0.142698 0 \n",
2519 "78 Silicon 10512.2 3.43512 678.146 0.214046 0 \n",
2520 "79 WCu 10513.6 3.43512 678.236 2.590920 0 \n",
2521 "80 Copper 10519.6 3.43512 678.624 8.913600 0 \n",
2522 "81 Air 10523.0 3.43512 678.840 0.000685 0 \n",
2523 "82 StainlessSteel 10524.0 3.43512 678.904 1.412470 0 \n",
2524 "83 Lead 10545.0 3.43512 680.259 30.698500 0 \n",
2525 "84 StainlessSteel 10546.0 3.43512 680.323 1.411610 0 \n",
2526 "85 Air 10549.3 3.43512 680.539 0.000684 0 \n",
2527 "86 Copper 10555.3 3.43512 680.927 8.903040 0 \n",
2528 "87 WCu 10556.8 3.43512 681.017 2.587410 0 \n",
2529 "88 Silicon 10556.9 3.43512 681.028 0.214043 0 \n",
2530 "89 Silicon 10557.0 3.43512 681.036 0.142696 0 \n",
2531 "90 HGC_G10-FR4 10558.6 3.43512 681.139 0.522997 0 \n",
2532 "91 Air 10562.0 3.43512 681.356 0.000684 0 \n",
2533 "92 StainlessSteel 10563.0 3.43512 681.420 1.410940 0 \n",
2534 "93 Lead 10584.0 3.43512 682.775 30.668000 0 \n",
2535 "94 StainlessSteel 10585.0 3.43512 682.839 1.410080 0 \n",
2536 "95 Air 10588.3 3.43512 683.055 0.000683 0 \n",
2537 "96 HGC_G10-FR4 10589.9 3.43512 683.159 0.522460 0 \n",
2538 "97 Silicon 10590.1 3.43512 683.166 0.142680 0 \n",
2539 "98 Silicon 10590.2 3.43512 683.178 0.214020 0 \n",
2540 "99 WCu 10591.6 3.43512 683.268 2.584590 0 \n",
2541 "100 Copper 10597.6 3.43512 683.655 8.890630 0 \n",
2542 "101 Air 10601.0 3.43512 683.871 0.000683 0 \n",
2543 "102 StainlessSteel 10602.0 3.43512 683.936 1.409420 0 \n",
2544 "103 Lead 10623.0 3.43512 685.291 30.637800 0 \n",
2545 "104 StainlessSteel 10624.0 3.43512 685.355 1.408570 0 \n",
2546 "105 Air 10627.3 3.43512 685.571 0.000682 0 \n",
2547 "106 Copper 10633.3 3.43512 685.958 8.880230 0 \n",
2548 "107 WCu 10634.8 3.43512 686.049 2.581140 0 \n",
2549 "108 Silicon 10634.9 3.43512 686.060 0.214078 0 \n",
2550 "109 Silicon 10635.0 3.43512 686.068 0.142719 0 \n",
2551 "110 HGC_G10-FR4 10636.6 3.43512 686.171 0.521664 0 \n",
2552 "111 Air 10640.0 3.43512 686.387 0.000682 0 \n",
2553 "112 StainlessSteel 10641.0 3.43512 686.452 1.407910 0 \n",
2554 "113 Lead 10662.0 3.43512 687.807 30.607900 0 \n",
2555 "114 StainlessSteel 10663.0 3.43512 687.871 1.407070 0 \n",
2556 "115 Air 10666.3 3.43512 688.087 0.000682 0 \n",
2557 "116 HGC_G10-FR4 10667.9 3.43512 688.190 0.521134 0 \n",
2558 "117 Silicon 10668.1 3.43512 688.198 0.142755 0 \n",
2559 "118 Silicon 10668.2 3.43512 688.210 0.214132 0 \n",
2560 "119 WCu 10669.6 3.43512 688.300 2.578370 0 \n",
2561 "120 Copper 10675.6 3.43512 688.687 8.868010 0 \n",
2562 "121 Air 10679.0 3.43512 688.903 0.000681 0 \n",
2563 "122 StainlessSteel 10680.0 3.43512 688.968 1.406410 0 \n",
2564 "123 Lead 10701.0 3.43512 690.322 30.578200 0 \n",
2565 "124 StainlessSteel 10702.0 3.43512 690.387 1.405580 0 \n",
2566 "125 Air 10705.3 3.43512 690.603 0.000681 0 \n",
2567 "126 Copper 10711.3 3.43512 690.990 8.857770 0 \n",
2568 "127 WCu 10712.8 3.43512 691.080 2.574970 0 \n",
2569 "128 Silicon 10712.9 3.43512 691.092 0.214095 0 \n",
2570 "129 Silicon 10713.0 3.43512 691.100 0.142730 0 \n",
2571 "130 HGC_G10-FR4 10714.6 3.43512 691.203 0.520350 0 \n",
2572 "131 Air 10718.0 3.43512 691.419 0.000680 0 \n",
2573 "132 StainlessSteel 10788.0 3.43512 695.935 98.308600 0 \n",
2574 "133 Air 10791.3 3.43512 696.151 0.000679 0 \n",
2575 "134 HGC_G10-FR4 10792.9 3.43512 696.254 0.519051 0 \n",
2576 "135 Silicon 10793.1 3.43512 696.262 0.142696 0 \n",
2577 "136 Silicon 10793.2 3.43512 696.273 0.214043 0 \n",
2578 "137 WCu 10794.6 3.43512 696.364 2.568580 0 \n",
2579 "138 Copper 10800.6 3.43512 696.751 8.832450 0 \n",
2580 "139 Air 10804.0 3.43512 696.967 0.000679 0 \n",
2581 "140 StainlessSteel 10874.0 3.43512 701.483 98.083000 0 \n",
2582 "141 Air 10877.3 3.43512 701.699 0.000677 0 \n",
2583 "142 Copper 10883.3 3.43512 702.086 8.809390 0 \n",
2584 "143 WCu 10884.8 3.43512 702.176 2.561680 0 \n",
2585 "144 Silicon 10884.9 3.43512 702.188 0.214112 0 \n",
2586 "145 Silicon 10885.0 3.43512 702.196 0.142742 0 \n",
2587 "146 HGC_G10-FR4 10886.6 3.43512 702.299 0.517521 0 \n",
2588 "147 Air 10890.0 3.43512 702.515 0.000677 0 \n",
2589 "148 Polyethylene 10484.0 0.00000 0.000 0.000000 0 \n",
2590 "\n",
2591 " PhysThickInmm tracknum PhysThickInXo dEdx dEdxtimesdx \\\n",
2592 "71 0.0 1.0 0.000000 0.197027 0.000000 \n",
2593 "72 1.0 1.0 0.057619 1.165737 1.165737 \n",
2594 "73 21.0 1.0 3.742115 1.273470 26.742870 \n",
2595 "74 1.0 1.0 0.057619 1.165737 1.165737 \n",
2596 "75 3.3 1.0 0.000011 0.000000 0.000000 \n",
2597 "76 1.6 1.0 0.009140 0.320491 0.512786 \n",
2598 "77 0.2 1.0 0.002135 0.387546 0.077509 \n",
2599 "78 0.1 1.0 0.001068 0.387546 0.038755 \n",
2600 "79 1.4 1.0 0.273304 1.811710 2.536394 \n",
2601 "80 6.0 1.0 0.417947 1.257088 7.542528 \n",
2602 "81 3.4 1.0 0.000011 0.000000 0.000000 \n",
2603 "82 1.0 1.0 0.057619 1.165737 1.165737 \n",
2604 "83 21.0 1.0 3.742115 1.273470 26.742870 \n",
2605 "84 1.0 1.0 0.057619 1.165737 1.165737 \n",
2606 "85 3.3 1.0 0.000011 0.000000 0.000000 \n",
2607 "86 6.0 1.0 0.417947 1.257088 7.542528 \n",
2608 "87 1.5 1.0 0.292826 1.811710 2.717565 \n",
2609 "88 0.1 1.0 0.001068 0.387546 0.038755 \n",
2610 "89 0.1 1.0 0.001068 0.387546 0.038755 \n",
2611 "90 1.6 1.0 0.009140 0.320491 0.512786 \n",
2612 "91 3.4 1.0 0.000011 0.000000 0.000000 \n",
2613 "92 1.0 1.0 0.057619 1.165737 1.165737 \n",
2614 "93 21.0 1.0 3.742115 1.273470 26.742870 \n",
2615 "94 1.0 1.0 0.057619 1.165737 1.165737 \n",
2616 "95 3.3 1.0 0.000011 0.000000 0.000000 \n",
2617 "96 1.6 1.0 0.009140 0.320491 0.512786 \n",
2618 "97 0.2 1.0 0.002135 0.387546 0.077509 \n",
2619 "98 0.1 1.0 0.001068 0.387546 0.038755 \n",
2620 "99 1.4 1.0 0.273304 1.811710 2.536394 \n",
2621 "100 6.0 1.0 0.417947 1.257088 7.542528 \n",
2622 "101 3.4 1.0 0.000011 0.000000 0.000000 \n",
2623 "102 1.0 1.0 0.057619 1.165737 1.165737 \n",
2624 "103 21.0 1.0 3.742115 1.273470 26.742870 \n",
2625 "104 1.0 1.0 0.057619 1.165737 1.165737 \n",
2626 "105 3.3 1.0 0.000011 0.000000 0.000000 \n",
2627 "106 6.0 1.0 0.417947 1.257088 7.542528 \n",
2628 "107 1.5 1.0 0.292826 1.811710 2.717565 \n",
2629 "108 0.1 1.0 0.001068 0.387546 0.038755 \n",
2630 "109 0.1 1.0 0.001068 0.387546 0.038755 \n",
2631 "110 1.6 1.0 0.009140 0.320491 0.512786 \n",
2632 "111 3.4 1.0 0.000011 0.000000 0.000000 \n",
2633 "112 1.0 1.0 0.057619 1.165737 1.165737 \n",
2634 "113 21.0 1.0 3.742115 1.273470 26.742870 \n",
2635 "114 1.0 1.0 0.057619 1.165737 1.165737 \n",
2636 "115 3.3 1.0 0.000011 0.000000 0.000000 \n",
2637 "116 1.6 1.0 0.009140 0.320491 0.512786 \n",
2638 "117 0.2 1.0 0.002135 0.387546 0.077509 \n",
2639 "118 0.1 1.0 0.001068 0.387546 0.038755 \n",
2640 "119 1.4 1.0 0.273304 1.811710 2.536394 \n",
2641 "120 6.0 1.0 0.417947 1.257088 7.542528 \n",
2642 "121 3.4 1.0 0.000011 0.000000 0.000000 \n",
2643 "122 1.0 1.0 0.057619 1.165737 1.165737 \n",
2644 "123 21.0 1.0 3.742115 1.273470 26.742870 \n",
2645 "124 1.0 1.0 0.057619 1.165737 1.165737 \n",
2646 "125 3.3 1.0 0.000011 0.000000 0.000000 \n",
2647 "126 6.0 1.0 0.417947 1.257088 7.542528 \n",
2648 "127 1.5 1.0 0.292826 1.811710 2.717565 \n",
2649 "128 0.1 1.0 0.001068 0.387546 0.038755 \n",
2650 "129 0.1 1.0 0.001068 0.387546 0.038755 \n",
2651 "130 1.6 1.0 0.009140 0.320491 0.512786 \n",
2652 "131 3.4 1.0 0.000011 0.000000 0.000000 \n",
2653 "132 70.0 1.0 4.033304 1.165737 81.601601 \n",
2654 "133 3.3 1.0 0.000011 0.000000 0.000000 \n",
2655 "134 1.6 1.0 0.009140 0.320491 0.512786 \n",
2656 "135 0.2 1.0 0.002135 0.387546 0.077509 \n",
2657 "136 0.1 1.0 0.001068 0.387546 0.038755 \n",
2658 "137 1.4 1.0 0.273304 1.811710 2.536394 \n",
2659 "138 6.0 1.0 0.417947 1.257088 7.542528 \n",
2660 "139 3.4 1.0 0.000011 0.000000 0.000000 \n",
2661 "140 70.0 1.0 4.033304 1.165737 81.601601 \n",
2662 "141 3.3 1.0 0.000011 0.000000 0.000000 \n",
2663 "142 6.0 1.0 0.417947 1.257088 7.542528 \n",
2664 "143 1.5 1.0 0.292826 1.811710 2.717565 \n",
2665 "144 0.1 1.0 0.001068 0.387546 0.038755 \n",
2666 "145 0.1 1.0 0.001068 0.387546 0.038755 \n",
2667 "146 1.6 1.0 0.009140 0.320491 0.512786 \n",
2668 "147 3.4 1.0 0.000011 0.000000 0.000000 \n",
2669 "148 6.0 1.0 0.011926 0.197027 1.182160 \n",
2670 "\n",
2671 " dEdxtimesdxCum EtableCum EfullCum PhysThickInXoCum layerflag layer \n",
2672 "71 0.000000 10.293100 0 0.000000 True 1.0 \n",
2673 "72 1.165737 11.707120 0 0.057619 True 1.0 \n",
2674 "73 27.908607 42.436220 0 3.799733 True 1.0 \n",
2675 "74 29.074344 43.849360 0 3.857352 True 1.0 \n",
2676 "75 29.074344 43.850045 0 3.857363 True 1.0 \n",
2677 "76 29.587131 44.373850 0 3.866503 True 1.0 \n",
2678 "77 29.664640 44.516548 0 3.868638 False 1.0 \n",
2679 "78 29.703394 44.730594 0 3.869705 False 1.0 \n",
2680 "79 32.239788 47.321514 0 4.143010 True 2.0 \n",
2681 "80 39.782316 56.235114 0 4.560956 True 2.0 \n",
2682 "81 39.782316 56.235799 0 4.560967 True 2.0 \n",
2683 "82 40.948053 57.648269 0 4.618586 True 2.0 \n",
2684 "83 67.690923 88.346769 0 8.360701 True 2.0 \n",
2685 "84 68.856661 89.758379 0 8.418320 True 2.0 \n",
2686 "85 68.856661 89.759063 0 8.418330 True 2.0 \n",
2687 "86 76.399189 98.662103 0 8.836277 True 2.0 \n",
2688 "87 79.116754 101.249513 0 9.129103 True 2.0 \n",
2689 "88 79.155508 101.463556 0 9.130170 False 2.0 \n",
2690 "89 79.194263 101.606252 0 9.131238 False 2.0 \n",
2691 "90 79.707049 102.129249 0 9.140378 True 3.0 \n",
2692 "91 79.707049 102.129933 0 9.140389 True 3.0 \n",
2693 "92 80.872786 103.540873 0 9.198008 True 3.0 \n",
2694 "93 107.615656 134.208873 0 12.940123 True 3.0 \n",
2695 "94 108.781393 135.618953 0 12.997741 True 3.0 \n",
2696 "95 108.781393 135.619636 0 12.997752 True 3.0 \n",
2697 "96 109.294179 136.142096 0 13.006892 True 3.0 \n",
2698 "97 109.371689 136.284776 0 13.009027 False 3.0 \n",
2699 "98 109.410443 136.498796 0 13.010095 False 3.0 \n",
2700 "99 111.946837 139.083386 0 13.283399 True 4.0 \n",
2701 "100 119.489365 147.974016 0 13.701345 True 4.0 \n",
2702 "101 119.489365 147.974699 0 13.701356 True 4.0 \n",
2703 "102 120.655102 149.384119 0 13.758975 True 4.0 \n",
2704 "103 147.397972 180.021919 0 17.501090 True 4.0 \n",
2705 "104 148.563710 181.430489 0 17.558709 True 4.0 \n",
2706 "105 148.563710 181.431172 0 17.558720 True 4.0 \n",
2707 "106 156.106238 190.311402 0 17.976666 True 4.0 \n",
2708 "107 158.823803 192.892542 0 18.269492 True 4.0 \n",
2709 "108 158.862557 193.106620 0 18.270559 False 4.0 \n",
2710 "109 158.901312 193.249339 0 18.271627 False 4.0 \n",
2711 "110 159.414098 193.771003 0 18.280767 True 5.0 \n",
2712 "111 159.414098 193.771685 0 18.280778 True 5.0 \n",
2713 "112 160.579835 195.179595 0 18.338397 True 5.0 \n",
2714 "113 187.322705 225.787495 0 22.080512 True 5.0 \n",
2715 "114 188.488442 227.194565 0 22.138130 True 5.0 \n",
2716 "115 188.488442 227.195247 0 22.138141 True 5.0 \n",
2717 "116 189.001228 227.716381 0 22.147281 True 5.0 \n",
2718 "117 189.078738 227.859136 0 22.149416 False 5.0 \n",
2719 "118 189.117492 228.073268 0 22.150484 False 5.0 \n",
2720 "119 191.653886 230.651638 0 22.423788 True 6.0 \n",
2721 "120 199.196414 239.519648 0 22.841734 True 6.0 \n",
2722 "121 199.196414 239.520329 0 22.841746 True 6.0 \n",
2723 "122 200.362151 240.926739 0 22.899364 True 6.0 \n",
2724 "123 227.105021 271.504939 0 26.641479 True 6.0 \n",
2725 "124 228.270759 272.910519 0 26.699098 True 6.0 \n",
2726 "125 228.270759 272.911200 0 26.699109 True 6.0 \n",
2727 "126 235.813287 281.768970 0 27.117055 True 6.0 \n",
2728 "127 238.530852 284.343940 0 27.409881 True 6.0 \n",
2729 "128 238.569606 284.558035 0 27.410948 False 6.0 \n",
2730 "129 238.608361 284.700765 0 27.412016 False 6.0 \n",
2731 "130 239.121147 285.221115 0 27.421156 True 7.0 \n",
2732 "131 239.121147 285.221795 0 27.421167 True 7.0 \n",
2733 "132 320.722748 383.530395 0 31.454471 True 7.0 \n",
2734 "133 320.722748 383.531074 0 31.454482 True 7.0 \n",
2735 "134 321.235534 384.050125 0 31.463622 True 7.0 \n",
2736 "135 321.313043 384.192821 0 31.465757 False 7.0 \n",
2737 "136 321.351798 384.406864 0 31.466824 False 7.0 \n",
2738 "137 323.888192 386.975444 0 31.740128 True 8.0 \n",
2739 "138 331.430720 395.807894 0 32.158075 True 8.0 \n",
2740 "139 331.430720 395.808573 0 32.158086 True 8.0 \n",
2741 "140 413.032321 493.891573 0 36.191390 True 8.0 \n",
2742 "141 413.032321 493.892250 0 36.191401 True 8.0 \n",
2743 "142 420.574849 502.701640 0 36.609347 True 8.0 \n",
2744 "143 423.292414 505.263320 0 36.902173 True 8.0 \n",
2745 "144 423.331168 505.477432 0 36.903241 False 8.0 \n",
2746 "145 423.369923 505.620174 0 36.904308 False 8.0 \n",
2747 "146 423.882709 506.137695 0 36.913448 True 9.0 \n",
2748 "147 423.882709 506.138372 0 36.913459 True 9.0 \n",
2749 "148 425.064869 506.138372 0 36.925385 True 9.0 "
2750 ]
2751 },
2752 "execution_count": 21,
2753 "metadata": {},
2754 "output_type": "execute_result"
2755 }
2756 ],
2757 "source": [
2758 "#display(matZ)\n",
2759 "#matZ[(matZ[\"layer\"] == 1) & matZ[\"tracknum\"] == 20]\n",
2760 "#pd.set_option('display.max_columns', None)\n",
2761 "pd.set_option('display.max_rows', None)\n",
2762 "#matZ[(matZ[\"tracknum\"] == 6.0)]\n",
2763 "matZ = matZ.dropna()\n",
2764 "matZ[ matZ['tracknum'] == 1.0]\n",
2765 "#matZ[(matZ[\"Mat\"] == \"Air\")]\n",
2766 "#matZ[ (matZ[\"Z\"] > 3950) & (matZ[\"Z\"] < 4000) & (matZ[\"Mat\"] == \"Silicon\") & (matZ[\"R\"] < 900)]\n",
2767 "#matZ[ (matZ[\"Mat\"] == \"Silicon\") & (matZ[\"tracknum\"] == 5) ]\n",
2768 "#matZ[ (matZ[\"Mat\"] == \"Silicon\") & (matZ[\"Z\"] < 3900) & (matZ[\"Z\"] > 0) & (matZ[\"PhysThickInmm\"] > 0.22)]\n",
2769 "#matZ[ (matZ[\"Mat\"] == \"Silicon\") & (matZ[\"layer\"] > 28) & (matZ[\"layer\"] < 36)& (matZ[\"PhysThickInmm\"] < 0.22)]\n",
2770 "#matZ[matZ[\"layer\"] >= 36]\n",
2771 "#matZ[(matZ[\"Mat\"] == \"WCu\") & (matZ[\"Z\"] > 3950) & (matZ[\"Z\"] < 4000)]"
2772 ]
2773 },
2774 {
2775 "cell_type": "code",
2776 "execution_count": 22,
2777 "metadata": {},
2778 "outputs": [],
2779 "source": [
2780 "#matZ[matZ[\"Z\"]>0]\n",
2781 "#matZ[matZ[\"PhysThickInmm\"]> 2000.]\n",
2782 "#matZ[matZ[\"tracknum\"] == 2 ]\n",
2783 "#matZ[ (matZ[\"tracknum\"] == 20) & (matZ[\"layer\"] > 10) ]\n",
2784 "#matZ[450:500]\n",
2785 "#For the D28 cumulative material in front of layers. \n",
2786 "newmatZ = matZ[ (matZ['tracknum'] == 1.0) & (matZ['layerflag'] == True) & (matZ['layerflag'].shift(-1) == False)]\n",
2787 "newmatZ[['layer', 'PhysThickInXoCum']]\n",
2788 "newmatZ[['layer', 'PhysThickInXoCum']].round(3).to_csv(r'/eos/home-d/dalfonso/D44.cumulative.xo', header=None, index=None, sep=' ', mode='a')"
2789 ]
2790 },
2791 {
2792 "cell_type": "code",
2793 "execution_count": 23,
2794 "metadata": {},
2795 "outputs": [],
2796 "source": [
2797 "# The dEdx weights calculation doesn't include the sensitive material. \n",
2798 "matZ = matZ.drop(matZ[ (matZ.Mat == \"Silicon\") | (matZ.Mat == \"H_Scintillator\") ].index)"
2799 ]
2800 },
2801 {
2802 "cell_type": "code",
2803 "execution_count": 24,
2804 "metadata": {},
2805 "outputs": [],
2806 "source": [
2807 "# We will add a column with the dedx weights. First dedxtimesdx sum per layer \n",
2808 "# This is with the theoretical dedx\n",
2809 "mdEdxtimesdxsumperlayer = matZ.groupby(['tracknum','layer'])[\"dEdxtimesdx\"].sum()\n",
2810 "# And now with the detailed simulation\n",
2811 "mdEdxtimesdxsumperlayer_detailedtable = matZ.groupby(['tracknum','layer'])[\"Etable\"].sum()\n",
2812 "mdEdxtimesdxsumperlayer_detailedfull = matZ.groupby(['tracknum','layer'])[\"Efull\"].sum()\n",
2813 "\n",
2814 "#type(mdEdxtimesdxsumperlayer)\n",
2815 "mdEdxtimesdxsumperlayer = mdEdxtimesdxsumperlayer.to_frame().reset_index()\n",
2816 "mdEdxtimesdxsumperlayer_detailedtable = mdEdxtimesdxsumperlayer_detailedtable.to_frame().reset_index()\n",
2817 "mdEdxtimesdxsumperlayer_detailedfull = mdEdxtimesdxsumperlayer_detailedfull.to_frame().reset_index()\n",
2818 "\n",
2819 "#mdEdxtimesdxsumperlayer\n",
2820 "#type(mdEdxtimesdxsumperlayer)"
2821 ]
2822 },
2823 {
2824 "cell_type": "code",
2825 "execution_count": 25,
2826 "metadata": {},
2827 "outputs": [],
2828 "source": [
2829 "#Let's put also the accumulated energy loss\n",
2830 "mdEdxtimesdxsumperlayer[\"dEdxtimesdxCum\"] = mdEdxtimesdxsumperlayer.groupby('tracknum')[\"dEdxtimesdx\"].cumsum()\n",
2831 "mdEdxtimesdxsumperlayer_detailedtable[\"EtableCum\"] = mdEdxtimesdxsumperlayer_detailedtable.groupby('tracknum')[\"Etable\"].cumsum()\n",
2832 "mdEdxtimesdxsumperlayer_detailedfull[\"EfullCum\"] = mdEdxtimesdxsumperlayer_detailedfull.groupby('tracknum')[\"Efull\"].cumsum()"
2833 ]
2834 },
2835 {
2836 "cell_type": "code",
2837 "execution_count": 26,
2838 "metadata": {},
2839 "outputs": [],
2840 "source": [
2841 "#Final weights\n",
2842 "mdEdxtimesdxsumperlayer[\"dedxweights\"] = (mdEdxtimesdxsumperlayer[\"dEdxtimesdx\"] + mdEdxtimesdxsumperlayer[\"dEdxtimesdx\"].shift(-1))/2\n",
2843 "mdEdxtimesdxsumperlayer_detailedtable[\"dedxweights_detailedsimulationtable\"] = (mdEdxtimesdxsumperlayer_detailedtable[\"Etable\"] + mdEdxtimesdxsumperlayer_detailedtable[\"Etable\"].shift(-1))/2\n",
2844 "mdEdxtimesdxsumperlayer_detailedfull[\"dedxweights_detailedsimulationfull\"] = (mdEdxtimesdxsumperlayer_detailedfull[\"Efull\"] + mdEdxtimesdxsumperlayer_detailedfull[\"Efull\"].shift(-1))/2"
2845 ]
2846 },
2847 {
2848 "cell_type": "code",
2849 "execution_count": 27,
2850 "metadata": {},
2851 "outputs": [
2852 {
2853 "data": {
2854 "text/html": [
2855 "<div>\n",
2856 "<style scoped>\n",
2857 " .dataframe tbody tr th:only-of-type {\n",
2858 " vertical-align: middle;\n",
2859 " }\n",
2860 "\n",
2861 " .dataframe tbody tr th {\n",
2862 " vertical-align: top;\n",
2863 " }\n",
2864 "\n",
2865 " .dataframe thead th {\n",
2866 " text-align: right;\n",
2867 " }\n",
2868 "</style>\n",
2869 "<table border=\"1\" class=\"dataframe\">\n",
2870 " <thead>\n",
2871 " <tr style=\"text-align: right;\">\n",
2872 " <th></th>\n",
2873 " <th>tracknum</th>\n",
2874 " <th>layer</th>\n",
2875 " <th>dEdxtimesdx</th>\n",
2876 " <th>dEdxtimesdxCum</th>\n",
2877 " <th>dedxweights</th>\n",
2878 " <th>Etable</th>\n",
2879 " <th>EtableCum</th>\n",
2880 " <th>dedxweights_detailedsimulationtable</th>\n",
2881 " <th>Efull</th>\n",
2882 " <th>EfullCum</th>\n",
2883 " <th>dedxweights_detailedsimulationfull</th>\n",
2884 " <th>EtableoverEtheory</th>\n",
2885 " <th>EfulloverEtheory</th>\n",
2886 " </tr>\n",
2887 " </thead>\n",
2888 " <tbody>\n",
2889 " <tr>\n",
2890 " <th>9</th>\n",
2891 " <td>1.0</td>\n",
2892 " <td>1.0</td>\n",
2893 " <td>29.587131</td>\n",
2894 " <td>29.587131</td>\n",
2895 " <td>39.500245</td>\n",
2896 " <td>44.373850</td>\n",
2897 " <td>44.373850</td>\n",
2898 " <td>50.446385</td>\n",
2899 " <td>0</td>\n",
2900 " <td>0</td>\n",
2901 " <td>0.000000</td>\n",
2902 " <td>1.499769</td>\n",
2903 " <td>0.0</td>\n",
2904 " </tr>\n",
2905 " <tr>\n",
2906 " <th>10</th>\n",
2907 " <td>1.0</td>\n",
2908 " <td>2.0</td>\n",
2909 " <td>49.413359</td>\n",
2910 " <td>79.000490</td>\n",
2911 " <td>39.756638</td>\n",
2912 " <td>56.518919</td>\n",
2913 " <td>100.892769</td>\n",
2914 " <td>45.527382</td>\n",
2915 " <td>0</td>\n",
2916 " <td>0</td>\n",
2917 " <td>0.000000</td>\n",
2918 " <td>1.143798</td>\n",
2919 " <td>0.0</td>\n",
2920 " </tr>\n",
2921 " <tr>\n",
2922 " <th>11</th>\n",
2923 " <td>1.0</td>\n",
2924 " <td>3.0</td>\n",
2925 " <td>30.099917</td>\n",
2926 " <td>109.100407</td>\n",
2927 " <td>39.756638</td>\n",
2928 " <td>34.535844</td>\n",
2929 " <td>135.428613</td>\n",
2930 " <td>45.464795</td>\n",
2931 " <td>0</td>\n",
2932 " <td>0</td>\n",
2933 " <td>0.000000</td>\n",
2934 " <td>1.147373</td>\n",
2935 " <td>0.0</td>\n",
2936 " </tr>\n",
2937 " <tr>\n",
2938 " <th>12</th>\n",
2939 " <td>1.0</td>\n",
2940 " <td>4.0</td>\n",
2941 " <td>49.413359</td>\n",
2942 " <td>158.513766</td>\n",
2943 " <td>39.756638</td>\n",
2944 " <td>56.393746</td>\n",
2945 " <td>191.822359</td>\n",
2946 " <td>45.430394</td>\n",
2947 " <td>0</td>\n",
2948 " <td>0</td>\n",
2949 " <td>0.000000</td>\n",
2950 " <td>1.141265</td>\n",
2951 " <td>0.0</td>\n",
2952 " </tr>\n",
2953 " <tr>\n",
2954 " <th>13</th>\n",
2955 " <td>1.0</td>\n",
2956 " <td>5.0</td>\n",
2957 " <td>30.099917</td>\n",
2958 " <td>188.613683</td>\n",
2959 " <td>39.756638</td>\n",
2960 " <td>34.467042</td>\n",
2961 " <td>226.289401</td>\n",
2962 " <td>45.368857</td>\n",
2963 " <td>0</td>\n",
2964 " <td>0</td>\n",
2965 " <td>0.000000</td>\n",
2966 " <td>1.145088</td>\n",
2967 " <td>0.0</td>\n",
2968 " </tr>\n",
2969 " <tr>\n",
2970 " <th>14</th>\n",
2971 " <td>1.0</td>\n",
2972 " <td>6.0</td>\n",
2973 " <td>49.413359</td>\n",
2974 " <td>238.027042</td>\n",
2975 " <td>66.020266</td>\n",
2976 " <td>56.270672</td>\n",
2977 " <td>282.560073</td>\n",
2978 " <td>77.810016</td>\n",
2979 " <td>0</td>\n",
2980 " <td>0</td>\n",
2981 " <td>0.000000</td>\n",
2982 " <td>1.138774</td>\n",
2983 " <td>0.0</td>\n",
2984 " </tr>\n",
2985 " <tr>\n",
2986 " <th>15</th>\n",
2987 " <td>1.0</td>\n",
2988 " <td>7.0</td>\n",
2989 " <td>82.627173</td>\n",
2990 " <td>320.654216</td>\n",
2991 " <td>92.283895</td>\n",
2992 " <td>99.349360</td>\n",
2993 " <td>381.909433</td>\n",
2994 " <td>110.102908</td>\n",
2995 " <td>0</td>\n",
2996 " <td>0</td>\n",
2997 " <td>0.000000</td>\n",
2998 " <td>1.202381</td>\n",
2999 " <td>0.0</td>\n",
3000 " </tr>\n",
3001 " <tr>\n",
3002 " <th>16</th>\n",
3003 " <td>1.0</td>\n",
3004 " <td>8.0</td>\n",
3005 " <td>101.940616</td>\n",
3006 " <td>422.594832</td>\n",
3007 " <td>92.283895</td>\n",
3008 " <td>120.856456</td>\n",
3009 " <td>502.765889</td>\n",
3010 " <td>91.800143</td>\n",
3011 " <td>0</td>\n",
3012 " <td>0</td>\n",
3013 " <td>98.971647</td>\n",
3014 " <td>1.185557</td>\n",
3015 " <td>0.0</td>\n",
3016 " </tr>\n",
3017 " </tbody>\n",
3018 "</table>\n",
3019 "</div>"
3020 ],
3021 "text/plain": [
3022 " tracknum layer dEdxtimesdx dEdxtimesdxCum dedxweights Etable \\\n",
3023 "9 1.0 1.0 29.587131 29.587131 39.500245 44.373850 \n",
3024 "10 1.0 2.0 49.413359 79.000490 39.756638 56.518919 \n",
3025 "11 1.0 3.0 30.099917 109.100407 39.756638 34.535844 \n",
3026 "12 1.0 4.0 49.413359 158.513766 39.756638 56.393746 \n",
3027 "13 1.0 5.0 30.099917 188.613683 39.756638 34.467042 \n",
3028 "14 1.0 6.0 49.413359 238.027042 66.020266 56.270672 \n",
3029 "15 1.0 7.0 82.627173 320.654216 92.283895 99.349360 \n",
3030 "16 1.0 8.0 101.940616 422.594832 92.283895 120.856456 \n",
3031 "\n",
3032 " EtableCum dedxweights_detailedsimulationtable Efull EfullCum \\\n",
3033 "9 44.373850 50.446385 0 0 \n",
3034 "10 100.892769 45.527382 0 0 \n",
3035 "11 135.428613 45.464795 0 0 \n",
3036 "12 191.822359 45.430394 0 0 \n",
3037 "13 226.289401 45.368857 0 0 \n",
3038 "14 282.560073 77.810016 0 0 \n",
3039 "15 381.909433 110.102908 0 0 \n",
3040 "16 502.765889 91.800143 0 0 \n",
3041 "\n",
3042 " dedxweights_detailedsimulationfull EtableoverEtheory EfulloverEtheory \n",
3043 "9 0.000000 1.499769 0.0 \n",
3044 "10 0.000000 1.143798 0.0 \n",
3045 "11 0.000000 1.147373 0.0 \n",
3046 "12 0.000000 1.141265 0.0 \n",
3047 "13 0.000000 1.145088 0.0 \n",
3048 "14 0.000000 1.138774 0.0 \n",
3049 "15 0.000000 1.202381 0.0 \n",
3050 "16 98.971647 1.185557 0.0 "
3051 ]
3052 },
3053 "execution_count": 27,
3054 "metadata": {},
3055 "output_type": "execute_result"
3056 }
3057 ],
3058 "source": [
3059 "#Hack for the last layer\n",
3060 "mdEdxtimesdxsumperlayer.loc[mdEdxtimesdxsumperlayer[\"layer\"] == 8.0, \"dedxweights\"] = 92.283895\n",
3061 "mdEdxtimesdxsumperlayer_detailedtable.loc[mdEdxtimesdxsumperlayer_detailedtable[\"layer\"] == 8.0, \"dedxweights_detailedsimulationtable\"] = 91.800143\n",
3062 "mdEdxtimesdxsumperlayer_detailedfull.loc[mdEdxtimesdxsumperlayer_detailedfull[\"layer\"] == 8.0, \"dedxweights_detailedsimulationfull\"] = 98.971647\n",
3063 "\n",
3064 "#Drop duplicates not needed columns\n",
3065 "mdEdxtimesdxsumperlayer_detailedtable = mdEdxtimesdxsumperlayer_detailedtable.drop(['tracknum', 'layer'], axis=1)\n",
3066 "mdEdxtimesdxsumperlayer_detailedfull = mdEdxtimesdxsumperlayer_detailedfull.drop(['tracknum', 'layer'], axis=1)\n",
3067 "\n",
3068 "#mdEdxtimesdxsumperlayer[ mdEdxtimesdxsumperlayer['tracknum'] == 6.0 ]\n",
3069 "#mdEdxtimesdxsumperlayer_detailed[ mdEdxtimesdxsumperlayer_detailed['tracknum'] == 6.0 ]\n",
3070 "mdEdxtimesdxsumperlayer = pd.concat([mdEdxtimesdxsumperlayer, mdEdxtimesdxsumperlayer_detailedtable, mdEdxtimesdxsumperlayer_detailedfull], axis=1)\n",
3071 "mdEdxtimesdxsumperlayer = mdEdxtimesdxsumperlayer.dropna()\n",
3072 "\n",
3073 "#Adding two columns dedxtable/dedxtheory and dedxfull/dedxtheory\n",
3074 "mdEdxtimesdxsumperlayer[\"EtableoverEtheory\"] = mdEdxtimesdxsumperlayer[\"Etable\"]/mdEdxtimesdxsumperlayer[\"dEdxtimesdx\"]\n",
3075 "mdEdxtimesdxsumperlayer[\"EfulloverEtheory\"] = mdEdxtimesdxsumperlayer[\"Efull\"]/mdEdxtimesdxsumperlayer[\"dEdxtimesdx\"]\n",
3076 "mdEdxtimesdxsumperlayer[ mdEdxtimesdxsumperlayer['tracknum'] == 1.0 ][:-1].dropna()"
3077 ]
3078 },
3079 {
3080 "cell_type": "code",
3081 "execution_count": 28,
3082 "metadata": {},
3083 "outputs": [
3084 {
3085 "data": {
3086 "text/html": [
3087 "<div>\n",
3088 "<style scoped>\n",
3089 " .dataframe tbody tr th:only-of-type {\n",
3090 " vertical-align: middle;\n",
3091 " }\n",
3092 "\n",
3093 " .dataframe tbody tr th {\n",
3094 " vertical-align: top;\n",
3095 " }\n",
3096 "\n",
3097 " .dataframe thead th {\n",
3098 " text-align: right;\n",
3099 " }\n",
3100 "</style>\n",
3101 "<table border=\"1\" class=\"dataframe\">\n",
3102 " <thead>\n",
3103 " <tr style=\"text-align: right;\">\n",
3104 " <th></th>\n",
3105 " <th>layer</th>\n",
3106 " <th>EtableoverEtheory</th>\n",
3107 " <th>EfulloverEtheory</th>\n",
3108 " </tr>\n",
3109 " </thead>\n",
3110 " <tbody>\n",
3111 " <tr>\n",
3112 " <th>9</th>\n",
3113 " <td>1.0</td>\n",
3114 " <td>1.499769</td>\n",
3115 " <td>0.0</td>\n",
3116 " </tr>\n",
3117 " <tr>\n",
3118 " <th>10</th>\n",
3119 " <td>2.0</td>\n",
3120 " <td>1.143798</td>\n",
3121 " <td>0.0</td>\n",
3122 " </tr>\n",
3123 " <tr>\n",
3124 " <th>11</th>\n",
3125 " <td>3.0</td>\n",
3126 " <td>1.147373</td>\n",
3127 " <td>0.0</td>\n",
3128 " </tr>\n",
3129 " <tr>\n",
3130 " <th>12</th>\n",
3131 " <td>4.0</td>\n",
3132 " <td>1.141265</td>\n",
3133 " <td>0.0</td>\n",
3134 " </tr>\n",
3135 " <tr>\n",
3136 " <th>13</th>\n",
3137 " <td>5.0</td>\n",
3138 " <td>1.145088</td>\n",
3139 " <td>0.0</td>\n",
3140 " </tr>\n",
3141 " <tr>\n",
3142 " <th>14</th>\n",
3143 " <td>6.0</td>\n",
3144 " <td>1.138774</td>\n",
3145 " <td>0.0</td>\n",
3146 " </tr>\n",
3147 " <tr>\n",
3148 " <th>15</th>\n",
3149 " <td>7.0</td>\n",
3150 " <td>1.202381</td>\n",
3151 " <td>0.0</td>\n",
3152 " </tr>\n",
3153 " <tr>\n",
3154 " <th>16</th>\n",
3155 " <td>8.0</td>\n",
3156 " <td>1.185557</td>\n",
3157 " <td>0.0</td>\n",
3158 " </tr>\n",
3159 " </tbody>\n",
3160 "</table>\n",
3161 "</div>"
3162 ],
3163 "text/plain": [
3164 " layer EtableoverEtheory EfulloverEtheory\n",
3165 "9 1.0 1.499769 0.0\n",
3166 "10 2.0 1.143798 0.0\n",
3167 "11 3.0 1.147373 0.0\n",
3168 "12 4.0 1.141265 0.0\n",
3169 "13 5.0 1.145088 0.0\n",
3170 "14 6.0 1.138774 0.0\n",
3171 "15 7.0 1.202381 0.0\n",
3172 "16 8.0 1.185557 0.0"
3173 ]
3174 },
3175 "execution_count": 28,
3176 "metadata": {},
3177 "output_type": "execute_result"
3178 }
3179 ],
3180 "source": [
3181 "mdEdxtimesdxsumperlayer[[\"layer\",\"EtableoverEtheory\",\"EfulloverEtheory\"]][ mdEdxtimesdxsumperlayer['tracknum'] == 1.0 ][:-1] "
3182 ]
3183 },
3184 {
3185 "cell_type": "code",
3186 "execution_count": 29,
3187 "metadata": {},
3188 "outputs": [
3189 {
3190 "data": {
3191 "text/html": [
3192 "<div>\n",
3193 "<style scoped>\n",
3194 " .dataframe tbody tr th:only-of-type {\n",
3195 " vertical-align: middle;\n",
3196 " }\n",
3197 "\n",
3198 " .dataframe tbody tr th {\n",
3199 " vertical-align: top;\n",
3200 " }\n",
3201 "\n",
3202 " .dataframe thead th {\n",
3203 " text-align: right;\n",
3204 " }\n",
3205 "</style>\n",
3206 "<table border=\"1\" class=\"dataframe\">\n",
3207 " <thead>\n",
3208 " <tr style=\"text-align: right;\">\n",
3209 " <th></th>\n",
3210 " <th>layer</th>\n",
3211 " <th>dedxweights</th>\n",
3212 " <th>dedxweights_detailedsimulationtable</th>\n",
3213 " <th>dedxweights_detailedsimulationfull</th>\n",
3214 " </tr>\n",
3215 " </thead>\n",
3216 " <tbody>\n",
3217 " <tr>\n",
3218 " <th>9</th>\n",
3219 " <td>1.0</td>\n",
3220 " <td>39.500245</td>\n",
3221 " <td>50.446385</td>\n",
3222 " <td>0.000000</td>\n",
3223 " </tr>\n",
3224 " <tr>\n",
3225 " <th>10</th>\n",
3226 " <td>2.0</td>\n",
3227 " <td>39.756638</td>\n",
3228 " <td>45.527382</td>\n",
3229 " <td>0.000000</td>\n",
3230 " </tr>\n",
3231 " <tr>\n",
3232 " <th>11</th>\n",
3233 " <td>3.0</td>\n",
3234 " <td>39.756638</td>\n",
3235 " <td>45.464795</td>\n",
3236 " <td>0.000000</td>\n",
3237 " </tr>\n",
3238 " <tr>\n",
3239 " <th>12</th>\n",
3240 " <td>4.0</td>\n",
3241 " <td>39.756638</td>\n",
3242 " <td>45.430394</td>\n",
3243 " <td>0.000000</td>\n",
3244 " </tr>\n",
3245 " <tr>\n",
3246 " <th>13</th>\n",
3247 " <td>5.0</td>\n",
3248 " <td>39.756638</td>\n",
3249 " <td>45.368857</td>\n",
3250 " <td>0.000000</td>\n",
3251 " </tr>\n",
3252 " <tr>\n",
3253 " <th>14</th>\n",
3254 " <td>6.0</td>\n",
3255 " <td>66.020266</td>\n",
3256 " <td>77.810016</td>\n",
3257 " <td>0.000000</td>\n",
3258 " </tr>\n",
3259 " <tr>\n",
3260 " <th>15</th>\n",
3261 " <td>7.0</td>\n",
3262 " <td>92.283895</td>\n",
3263 " <td>110.102908</td>\n",
3264 " <td>0.000000</td>\n",
3265 " </tr>\n",
3266 " <tr>\n",
3267 " <th>16</th>\n",
3268 " <td>8.0</td>\n",
3269 " <td>92.283895</td>\n",
3270 " <td>91.800143</td>\n",
3271 " <td>98.971647</td>\n",
3272 " </tr>\n",
3273 " </tbody>\n",
3274 "</table>\n",
3275 "</div>"
3276 ],
3277 "text/plain": [
3278 " layer dedxweights dedxweights_detailedsimulationtable \\\n",
3279 "9 1.0 39.500245 50.446385 \n",
3280 "10 2.0 39.756638 45.527382 \n",
3281 "11 3.0 39.756638 45.464795 \n",
3282 "12 4.0 39.756638 45.430394 \n",
3283 "13 5.0 39.756638 45.368857 \n",
3284 "14 6.0 66.020266 77.810016 \n",
3285 "15 7.0 92.283895 110.102908 \n",
3286 "16 8.0 92.283895 91.800143 \n",
3287 "\n",
3288 " dedxweights_detailedsimulationfull \n",
3289 "9 0.000000 \n",
3290 "10 0.000000 \n",
3291 "11 0.000000 \n",
3292 "12 0.000000 \n",
3293 "13 0.000000 \n",
3294 "14 0.000000 \n",
3295 "15 0.000000 \n",
3296 "16 98.971647 "
3297 ]
3298 },
3299 "execution_count": 29,
3300 "metadata": {},
3301 "output_type": "execute_result"
3302 }
3303 ],
3304 "source": [
3305 "mdEdxtimesdxsumperlayer[[\"layer\",\"dedxweights\",\"dedxweights_detailedsimulationtable\",\"dedxweights_detailedsimulationfull\"]][ mdEdxtimesdxsumperlayer['tracknum'] == 1.0 ][:-1] "
3306 ]
3307 },
3308 {
3309 "cell_type": "code",
3310 "execution_count": 30,
3311 "metadata": {},
3312 "outputs": [
3313 {
3314 "data": {
3315 "text/html": [
3316 "<div>\n",
3317 "<style scoped>\n",
3318 " .dataframe tbody tr th:only-of-type {\n",
3319 " vertical-align: middle;\n",
3320 " }\n",
3321 "\n",
3322 " .dataframe tbody tr th {\n",
3323 " vertical-align: top;\n",
3324 " }\n",
3325 "\n",
3326 " .dataframe thead th {\n",
3327 " text-align: right;\n",
3328 " }\n",
3329 "</style>\n",
3330 "<table border=\"1\" class=\"dataframe\">\n",
3331 " <thead>\n",
3332 " <tr style=\"text-align: right;\">\n",
3333 " <th></th>\n",
3334 " <th>layer</th>\n",
3335 " <th>dEdxtimesdxCum</th>\n",
3336 " <th>EtableCum</th>\n",
3337 " <th>EfullCum</th>\n",
3338 " </tr>\n",
3339 " </thead>\n",
3340 " <tbody>\n",
3341 " <tr>\n",
3342 " <th>9</th>\n",
3343 " <td>1.0</td>\n",
3344 " <td>29.587131</td>\n",
3345 " <td>44.373850</td>\n",
3346 " <td>0</td>\n",
3347 " </tr>\n",
3348 " <tr>\n",
3349 " <th>10</th>\n",
3350 " <td>2.0</td>\n",
3351 " <td>79.000490</td>\n",
3352 " <td>100.892769</td>\n",
3353 " <td>0</td>\n",
3354 " </tr>\n",
3355 " <tr>\n",
3356 " <th>11</th>\n",
3357 " <td>3.0</td>\n",
3358 " <td>109.100407</td>\n",
3359 " <td>135.428613</td>\n",
3360 " <td>0</td>\n",
3361 " </tr>\n",
3362 " <tr>\n",
3363 " <th>12</th>\n",
3364 " <td>4.0</td>\n",
3365 " <td>158.513766</td>\n",
3366 " <td>191.822359</td>\n",
3367 " <td>0</td>\n",
3368 " </tr>\n",
3369 " <tr>\n",
3370 " <th>13</th>\n",
3371 " <td>5.0</td>\n",
3372 " <td>188.613683</td>\n",
3373 " <td>226.289401</td>\n",
3374 " <td>0</td>\n",
3375 " </tr>\n",
3376 " <tr>\n",
3377 " <th>14</th>\n",
3378 " <td>6.0</td>\n",
3379 " <td>238.027042</td>\n",
3380 " <td>282.560073</td>\n",
3381 " <td>0</td>\n",
3382 " </tr>\n",
3383 " <tr>\n",
3384 " <th>15</th>\n",
3385 " <td>7.0</td>\n",
3386 " <td>320.654216</td>\n",
3387 " <td>381.909433</td>\n",
3388 " <td>0</td>\n",
3389 " </tr>\n",
3390 " <tr>\n",
3391 " <th>16</th>\n",
3392 " <td>8.0</td>\n",
3393 " <td>422.594832</td>\n",
3394 " <td>502.765889</td>\n",
3395 " <td>0</td>\n",
3396 " </tr>\n",
3397 " </tbody>\n",
3398 "</table>\n",
3399 "</div>"
3400 ],
3401 "text/plain": [
3402 " layer dEdxtimesdxCum EtableCum EfullCum\n",
3403 "9 1.0 29.587131 44.373850 0\n",
3404 "10 2.0 79.000490 100.892769 0\n",
3405 "11 3.0 109.100407 135.428613 0\n",
3406 "12 4.0 158.513766 191.822359 0\n",
3407 "13 5.0 188.613683 226.289401 0\n",
3408 "14 6.0 238.027042 282.560073 0\n",
3409 "15 7.0 320.654216 381.909433 0\n",
3410 "16 8.0 422.594832 502.765889 0"
3411 ]
3412 },
3413 "execution_count": 30,
3414 "metadata": {},
3415 "output_type": "execute_result"
3416 }
3417 ],
3418 "source": [
3419 "mdEdxtimesdxsumperlayer[[\"layer\",\"dEdxtimesdxCum\",\"EtableCum\",\"EfullCum\"]][ mdEdxtimesdxsumperlayer['tracknum'] == 1.0 ][:-1] "
3420 ]
3421 },
3422 {
3423 "cell_type": "code",
3424 "execution_count": null,
3425 "metadata": {},
3426 "outputs": [],
3427 "source": []
3428 }
3429 ],
3430 "metadata": {
3431 "kernelspec": {
3432 "display_name": "Python 3",
3433 "language": "python",
3434 "name": "python3"
3435 },
3436 "language_info": {
3437 "codemirror_mode": {
3438 "name": "ipython",
3439 "version": 3
3440 },
3441 "file_extension": ".py",
3442 "mimetype": "text/x-python",
3443 "name": "python",
3444 "nbconvert_exporter": "python",
3445 "pygments_lexer": "ipythond3",
3446 "version": "3.7.1"
3447 }
3448 },
3449 "nbformat": 4,
3450 "nbformat_minor": 2
3451 }
