From e7e36ea362779ea876152febf93b2a27e484e8d9 Mon Sep 17 00:00:00 2001
From: Karim Ahmed <karim.ahmed@xfel.eu>
Date: Wed, 22 Jul 2020 17:15:59 +0200
Subject: [PATCH] read voltage, gain and temperature
---
.../pnCCD/Characterize_pnCCD_Dark_NBC.ipynb | 49 +++++++++++++------
1 file changed, 33 insertions(+), 16 deletions(-)
diff --git a/notebooks/pnCCD/Characterize_pnCCD_Dark_NBC.ipynb b/notebooks/pnCCD/Characterize_pnCCD_Dark_NBC.ipynb
index 10da19ed8..d9cbdd4e0 100644
--- a/notebooks/pnCCD/Characterize_pnCCD_Dark_NBC.ipynb
+++ b/notebooks/pnCCD/Characterize_pnCCD_Dark_NBC.ipynb
@@ -31,16 +31,19 @@
"outputs": [],
"source": [
"cluster_profile = \"noDB\" # ipcluster profile to use\n",
- "in_folder = \"/gpfs/exfel/exp/SQS/201930/p900075/raw\" # input folder, required\n",
+ "in_folder = \"/gpfs/exfel/exp/SQS/202031/p900166/raw\" # input folder, required\n",
"out_folder = '/gpfs/exfel/data/scratch/ahmedk/test/pnccd' # output folder, required\n",
"sequence = 0 # sequence file to use\n",
- "run = 364 # which run to read data from, required\n",
+ "run = 34 # which run to read data from, required\n",
"\n",
+ "db_module = \"PnCCD1\"\n",
"karabo_da = ['PNCCD01'] # data aggregators\n",
+ "karabo_da_control = \"PNCCD02\" # file inset for control data\n",
"karabo_id = \"SQS_NQS_PNCCD1MP\" # karabo prefix of PNCCD devices\n",
"receiver_id = \"PNCCD_FMT-0\" # inset for receiver devices\n",
"path_template = 'RAW-R{:04d}-{}-S{{:05d}}.h5' # the template to use to access data\n",
"h5path = '/INSTRUMENT/{}/CAL/{}:output/data/image/' # path in the HDF5 file the data is at\n",
+ "h5path_ctrl = '/CONTROL/{}/CTRL/TCTRL'\n",
"\n",
"# for database time derivation:\n",
"use_dir_creation_date = True # use dir creation date as data production reference date\n",
@@ -82,6 +85,7 @@
"import copy\n",
"import datetime\n",
"import warnings\n",
+ "import h5py\n",
"warnings.filterwarnings('ignore')\n",
"\n",
"import numpy as np\n",
@@ -107,17 +111,6 @@
"from cal_tools.enums import BadPixels"
]
},
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {},
- "outputs": [],
- "source": [
- "# Output Folder Creation:\n",
- "if not os.path.exists(out_folder):\n",
- " os.makedirs(out_folder)"
- ]
- },
{
"cell_type": "code",
"execution_count": null,
@@ -143,7 +136,7 @@
"outputs": [],
"source": [
"proposal = list(filter(None, in_folder.strip('/').split('/')))[-2]\n",
- "file_loc = 'Proposal: {}, Run: {}'.format(proposal, run)\n",
+ "file_loc = f'Proposal: {proposal}, Run: {run}'\n",
"print(\"File Location:\", file_loc)"
]
},
@@ -171,6 +164,9 @@
"filename = fp_path.format(sequence)\n",
"h5path = h5path.format(karabo_id, receiver_id)\n",
"\n",
+ "# Output Folder Creation:\n",
+ "os.makedirs(out_folder, exist_ok=True)\n",
+ "\n",
"# Run's creation time:\n",
"if creation_time:\n",
" try:\n",
@@ -196,6 +192,23 @@
"print(\"Run number: {}\".format(run))"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "# extract slow data 'RAW-R{:04d}-{}-S{{:05d}}.h5'\n",
+ "ctrl_fname = os.path.join(ped_dir, path_template.format(run, karabo_da_control)).format(sequence)\n",
+ "ctrl_path = h5path_ctrl.format(karabo_id)\n",
+ "mdl_ctrl_path = \"/CONTROL/{}/MDL/\".format(karabo_id)\n",
+ "\n",
+ "with h5py.File(ctrl_fname, \"r\") as f:\n",
+ " bias_voltage = abs(f[os.path.join(mdl_ctrl_path, \"DAQ_MPOD/u0voltage/value\")][0])\n",
+ " gain = f[os.path.join(mdl_ctrl_path, \"DAQ_GAIN/pNCCDGain/value\")][0]\n",
+ " fix_temperature = f[os.path.join(ctrl_path, \"inputA/krdg/value\")][0]"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
@@ -993,7 +1006,6 @@
" det = Constants.CCD(DetectorTypes.pnCCD)\n",
" const = getattr(det, const_name)()\n",
" const.data = constant_maps[const_name].data\n",
- "\n",
" metadata.calibration_constant = const\n",
"\n",
" # set the operating condition\n",
@@ -1004,7 +1016,12 @@
" pixels_x=pixels_x,\n",
" pixels_y=pixels_y)\n",
"\n",
- " device = Detectors.PnCCD1\n",
+ " for parm in condition.parameters:\n",
+ " if parm.name == \"Sensor Temperature\":\n",
+ " parm.lower_deviation = time_limits\n",
+ " parm.upper_deviation = time_limits\n",
+ "\n",
+ " device = device = getattr(Detectors, db_module)\n",
" metadata.detector_condition = condition\n",
"\n",
" # specify the a version for this constant\n",
--
GitLab