From 9d8a4650bcedc3e4dbe8a498be1fd25d1912f776 Mon Sep 17 00:00:00 2001 From: RealPolitiX Date: Tue, 20 Oct 2020 04:43:36 -0500 Subject: [PATCH] Added core-level notebook --- .../Tutorial_02_W4f_Core-level_WSe2.ipynb | 242 ++++++++++++++++++ 1 file changed, 242 insertions(+) create mode 100644 Tutorial/Tutorial_02_W4f_Core-level_WSe2.ipynb diff --git a/Tutorial/Tutorial_02_W4f_Core-level_WSe2.ipynb b/Tutorial/Tutorial_02_W4f_Core-level_WSe2.ipynb new file mode 100644 index 0000000..a4e4cde --- /dev/null +++ b/Tutorial/Tutorial_02_W4f_Core-level_WSe2.ipynb @@ -0,0 +1,242 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "import matplotlib.pyplot as plt\n", + "import sys, os\n", + "import numpy as np\n", + "import dask\n", + "import h5py\n", + "\n", + "sys.path.append(\".../hextof-processor\")\n", + "# sys.path.append('.../pah')\n", + "from processor.DldFlashDataframeCreator import DldFlashProcessor, utils" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "# Binning parameter settings\n", + "mb_start=100\n", + "mb_end=399\n", + "\n", + "ToFfrom = 119.5\n", + "ToFto = 124.5\n", + "ToFstep=0.01\n", + "\n", + "delayFrom=-516.15\n", + "delayTo=-512.3\n", + "delayStep=0.02\n", + "\n", + "# dld parameters\n", + "dldXFrom=30\n", + "dldXTo=2700\n", + "dldXStep=30\n", + "\n", + "# dld parameters\n", + "dldYFrom=0\n", + "dldYTo=2850\n", + "dldYStep=30" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "# Run no.22191\n", + "processor = DldFlashProcessor()\n", + "processor.runNumber = 22191\n", + "processor.readDataframes(format='parquet')\n", + "\n", + "processor.resetBins()\n", + "ToF = processor.addBinning('dldTime', ToFfrom,ToFto,ToFstep)\n", + "Delay = processor.addBinning('delayStage', delayFrom,delayTo,delayStep)\n", + "processor.dd = processor.dd[((processor.dd['dldPosX']<1290)|(processor.dd['dldPosX']>1410))]\n", + "processor.dd = processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd = processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd = processor.dd[((processor.dd['dldMicrobunchId']<=mb_end)&(processor.dd['dldMicrobunchId']>=mb_start))]\n", + "processor.dd['delayStage'] = processor.dd['delayStage'] - processor.dd['bam']\n", + "\n", + "processor.shiftQuadrants()\n", + "processor.filterCircleDLDPos()\n", + "processor.correctSpaceCharge()\n", + "processor.correctSpaceCharge2()\n", + "processor.dd.dldTime=(1567*1567)/((processor.dd.dldTime-338.216)*(processor.dd.dldTime-338.216))+100-0.01\n", + "result_low1 = processor.computeBinnedData()\n", + "\n", + "\n", + "processor = DldFlashProcessor()\n", + "processor.runNumber = 22191\n", + "processor.readDataframes(format='parquet')\n", + "\n", + "processor.resetBins()\n", + "Delay = processor.addBinning('delayStage', delayFrom,delayTo,delayStep)\n", + "processor.dd=processor.dd[((processor.dd['dldPosX']<1290)|(processor.dd['dldPosX']>1410))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldMicrobunchId']<=mb_end)&(processor.dd['dldMicrobunchId']>=mb_start))]\n", + "processor.dd['delayStage'] = processor.dd['delayStage'] - processor.dd['bam']\n", + "\n", + "processor.shiftQuadrants()\n", + "processor.filterCircleDLDPos()\n", + "processor.correctSpaceCharge()\n", + "processor.correctSpaceCharge2()\n", + "processor.dd.dldTime=(1567*1567)/((processor.dd.dldTime-338.216)*(processor.dd.dldTime-338.216))+100-0.01\n", + "processor.dd=processor.dd[((processor.dd['dldTime']<124.5)&(processor.dd['dldTime']>119.5))]\n", + "hist_low1 = processor.computeBinnedData()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "# Run no.22192\n", + "processor = DldFlashProcessor()\n", + "processor.runNumber = 22192\n", + "processor.readDataframes(format='parquet')\n", + "\n", + "processor.resetBins()\n", + "ToF = processor.addBinning('dldTime', ToFfrom,ToFto,ToFstep)\n", + "Delay = processor.addBinning('delayStage', delayFrom,delayTo,delayStep)\n", + "processor.dd=processor.dd[((processor.dd['dldPosX']<1290)|(processor.dd['dldPosX']>1410))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldMicrobunchId']<=mb_end)&(processor.dd['dldMicrobunchId']>=mb_start))]\n", + "processor.dd['delayStage'] = processor.dd['delayStage'] - processor.dd['bam']\n", + "\n", + "processor.shiftQuadrants()\n", + "processor.filterCircleDLDPos()\n", + "processor.correctSpaceCharge()\n", + "processor.correctSpaceCharge2()\n", + "processor.dd.dldTime=(1567*1567)/((processor.dd.dldTime-338.216)*(processor.dd.dldTime-338.216))+100\n", + "result_low2 = processor.computeBinnedData()\n", + "\n", + "\n", + "processor = DldFlashProcessor()\n", + "processor.runNumber = 22192\n", + "processor.readDataframes(format='parquet')\n", + "processor.resetBins()\n", + "\n", + "Delay = processor.addBinning('delayStage', delayFrom,delayTo,delayStep)\n", + "processor.dd=processor.dd[((processor.dd['dldPosX']<1290)|(processor.dd['dldPosX']>1410))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldMicrobunchId']<=mb_end)&(processor.dd['dldMicrobunchId']>=mb_start))]\n", + "processor.dd['delayStage'] = processor.dd['delayStage'] - processor.dd['bam']\n", + "\n", + "processor.shiftQuadrants()\n", + "processor.filterCircleDLDPos()\n", + "processor.correctSpaceCharge()\n", + "processor.correctSpaceCharge2()\n", + "processor.dd.dldTime=(1567*1567)/((processor.dd.dldTime-338.216)*(processor.dd.dldTime-338.216))+100\n", + "processor.dd=processor.dd[((processor.dd['dldTime']<124.5)&(processor.dd['dldTime']>119.5))]\n", + "hist_low2 = processor.computeBinnedData()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "# Run no.22193\n", + "processor = DldFlashProcessor()\n", + "processor.runNumber = 22193\n", + "processor.readDataframes(format='parquet')\n", + "\n", + "processor.resetBins()\n", + "ToF = processor.addBinning('dldTime', ToFfrom,ToFto,ToFstep)\n", + "Delay = processor.addBinning('delayStage', delayFrom,delayTo,delayStep)\n", + "processor.dd=processor.dd[((processor.dd['dldPosX']<1290)|(processor.dd['dldPosX']>1410))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldMicrobunchId']<=mb_end)&(processor.dd['dldMicrobunchId']>=mb_start))]\n", + "processor.dd['delayStage'] = processor.dd['delayStage'] - processor.dd['bam']\n", + "\n", + "processor.shiftQuadrants()\n", + "processor.filterCircleDLDPos()\n", + "processor.correctSpaceCharge()\n", + "processor.correctSpaceCharge2()\n", + "processor.dd.dldTime=(1567*1567)/((processor.dd.dldTime-338.216)*(processor.dd.dldTime-338.216))+100\n", + "result_low3 = processor.computeBinnedData()\n", + "\n", + "\n", + "processor = DldFlashProcessor()\n", + "processor.runNumber = 22193\n", + "processor.readDataframes(format='parquet')\n", + "\n", + "processor.resetBins()\n", + "Delay = processor.addBinning('delayStage', delayFrom,delayTo,delayStep)\n", + "processor.dd=processor.dd[((processor.dd['dldPosX']<1290)|(processor.dd['dldPosX']>1410))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldPosY']<1380)|(processor.dd['dldPosY']>1500))]\n", + "processor.dd=processor.dd[((processor.dd['dldMicrobunchId']<=mb_end)&(processor.dd['dldMicrobunchId']>=mb_start))]\n", + "processor.dd['delayStage'] = processor.dd['delayStage'] - processor.dd['bam']\n", + "\n", + "processor.shiftQuadrants()\n", + "processor.filterCircleDLDPos()\n", + "processor.correctSpaceCharge()\n", + "processor.correctSpaceCharge2()\n", + "processor.dd.dldTime=(1567*1567)/((processor.dd.dldTime-338.216)*(processor.dd.dldTime-338.216))+100\n", + "processor.dd = processor.dd[((processor.dd['dldTime']<124.5)&(processor.dd['dldTime']>119.5))]\n", + "hist_low3 = processor.computeBinnedData()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "a = hist_low1[50:150].mean()\n", + "b = hist_low2[50:150].mean()\n", + "c = hist_low3[50:150].mean()\n", + "result_norm_low1 = []\n", + "result_norm_low2 = []\n", + "result_norm_low3 = []\n", + "result_norm_low1 = (result_low1/hist_low1)*a\n", + "result_norm_low2 = (result_low2/hist_low2)*b\n", + "result_norm_low3 = (result_low3/hist_low3)*c" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.6.4" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +}