Changes for page The FLASH HDF5 structure

Last modified by sndueste on 2025/02/06 10:55

From 39.1 to 38.1 From 41.1 to 40.1

From version 40.1

edited by sndueste
on 2021/02/09 14:35

Change comment: There is no comment for this version

To version 39.1

edited by sndueste
on 2021/02/09 11:53

Change comment: There is no comment for this version

Raw
Rendered

Summary

Page properties (1 modified, 0 added, 0 removed)

Details

Page properties

Content

@@ -18,56 +18,8 @@
  [[Contents>>doc:||anchor="Contents"]]
--== The new (starting 2021) HDF5 format ==
++== The current FLASH HDF5 structure ==
--\\
--
--FLASH provides a conversion of its data acquisition (DAQ) to the commonly used [[HDF5>>url:https://www.hdfgroup.org/solutions/hdf5/||shape="rect"]] format. Correlated  data are mapped by a primary index called **train ID**. Every data set has an individual index of train IDs to identify the data even when data are missing or deviate in update rate.If the data set of choice contains gaps,  users have to decide how to treat missing values. DAQ channels are spread across various files with one file name pattern for each DAQ. This means users have to assemble data from different files if necessary.
--
--The hierarchy is realized using a human readable named HDF tree with each DAQ channel containing the data sets "//value//" and "//index//". Additionally, the HDF group //zraw// contains a tree of the original DAQ channel names.
--
--Reference implementation, which follows the concept of Python libraries like Pandas, Xarray, or Dask, is given below.
--
--\\
--
--{{expand title="Discontinued HDF formats"}}
--== Comparison to FLASH's deprecated HDF formats ==
--
--Before 2021, FLASH provided two different HDF formats formally known as //near-online// and //offline// HDF files.
--
--All data sets were aligned by the same global index by filling in missing data. Fast DAQ channels  (pulse synchronous data - update 10 Hz) are filled in by NaN or 0, slow channels (not pulse synchronous data - update e.g. 1 Hz ) are interpolated by the DAQs libraries, by keeping the value constant up to the next change.
--
--While the near-online HDF files were converted live during the beamtime, the offline HDF files were manually compiled weeks/months later. Near-online HDF files were used by a provided API [[BeamtimeDaqAccess>>url:https://confluence.desy.de/display/FLASHUSER/Near-Online+data+analysis#Near-Onlinedataanalysis-BriefUserGuidetoBeamtimeDAQAccess||shape="rect"]], which also searches for DAQ channels in files. With assembled files "by run", all DAQ channels were existing in the same HDF file. While shorter runs usually fitted into one file, longer runs still had to be aggregated over several files. The creation of the assembled type of files still involves the use of fairly unstable DAQ Mex-functions As the environment required for using the Mex-functions is already deprecated, its continued existence is uncertain.
--
--The HDF hierarchy is almost identical in all formats. While in the the recent format each DAQ channel contains the data sets "value" and "index", the deprecated format had one data set for each DAQ channel. The deprecated formats supplies no time axis parameters for spectra data types.
--{{/expand}}
--
--== HDF excerpt: ==
--
--\\
--
--ADC data as example for **fast** **data** (10 Hz):
--
--\\
--
--[[image:attach:image2020-11-16_15-26-28.png||height="250"]]
--
--\\
--
--{{info title="Sample scripts in python"}}
--== Reference implementation (Python) ==
--
--[[~[~[image:attach:binder_badge.png~|~|thumbnail="true" width="120"~]~]>>url:https://mybinder.org/v2/git/https%3A%2F%2Fgitlab.desy.de%2Fchristopher.passow%2Fflash-daq-hdf/master||shape="rect"]]
--
--(% class="Object" %)[[https:~~/~~/gitlab.desy.de/christopher.passow/flash-daq-hdf>>url:https://gitlab.desy.de/christopher.passow/flash-daq-hdf||shape="rect"]]
--{{/info}}
--
--\\
--
--\\
--
--== The discontinued (till 2021) FLASH HDF5 structure ==
--
  The photon diagnostic, electron diagnostic and beamline information as well as the information about the pump-probe laser and the infrastructure offered for users (GHz/MHz ADCs) can be included in one HDF5 file which is organized according to train IDs. The general structure is:
  * Electron Diagnostic
@@ -102,6 +102,14 @@
  \\
++== The new (starting 2021) HDF5 format ==
++
++Here is [[some documentation on the changes of the HDF5 format>>doc:FLASHUSER.HDF5 format]] that well be available in 2021 (work in progress)
++
++\\
++
++\\
++
  == Most popular FLASH parameters and their names in HDF5, DOOCS and (raw) DAQ  ==
  {{id name="DOOCSparameters"/}}
@@ -198,7 +198,7 @@
  (% style="color: rgb(0,0,0);" %)**Discontinued GMD data recording / evaluation  (VME + PhotonFlux ML server)**
  {{code language="none"}}/FL1/Photon Diagnostic/GMD/Average energy/energy tunnel{{/code}}
--//always saved (PBD)//
++//always saved (PBD)//\\
  (% style="color: rgb(0,0,0);" %)DOOCS prop : {{code language="none"}}TTF2.DAQ/PHFLUX/OUT04/VAL{{/code}} (%%)
  (% style="color: rgb(0,0,0);" %)DAQ channel: {{code language="none"}}PBD.PHFLUX/TUNNEL.ENPULSEIC{{/code}}(%%)
@@ -456,7 +456,7 @@
  ===== electron bunch energy =====
  {{code language="none"}}/FL1/Electron Diagnostic/Electron energy/average electron energy{{/code}}
--//always saved (PBD)//
++_always saved (PBD)
  DOOCS prop : {{code language="none"}}TTF2.DAQ/ENERGY.DOGLEG/E_INTRA_MEAN/VAL{{/code}}
  DAQ channel: {{code language="none"}}PBD.ENERGY.DOGLEG/E_MEAN{{/code}}
  desc: electron bunch energy (average over the bunch train)
@@ -464,7 +464,7 @@
  \\
  {{code language="none"}}/FL1/Electron Diagnostic/Electron energy/pulse resolved energy{{/code}}
--//always saved (PBD)//
++_always saved (PBD)
  DOOCS prop : {{code language="none"}}TTF2.DAQ/ENERGY.DOGLEG/E_SPECT/VAL.TD{{/code}}
  DAQ channel: {{code language="none"}}PBD.ENERGY.DOGLEG/E_SPECT{{/code}}
  desc: electron bunch energy bunch resolved
@@ -473,7 +473,7 @@
  \\
  {{code language="none"}}/FL1/Electron Diagnostic/Electron energy/wavelength bunch train average{{/code}}
--//always saved (PBD)//
++_always saved (PBD)
  DOOCS prop : {{code language="none"}}TTF2.DAQ/ENERGY.DOGLEG/LAMBDA_MEAN/VAL{{/code}}
  DAQ channel: {{code language="none"}}PBD.ENERGY.DOGLEG/LAMBDA_MEAN{{/code}}
  desc: Wavelength calculated by the electron bunch energy (average over the bunch train) (FLASH1)
@@ -667,12 +667,12 @@
  ==== Beamline info (FLASH2) ====
--{{code language="none"}}/FL2/Beamlines/Attenuator/pressure  {{/code}}
--(% style="color: rgb(0,0,0);" %)//always saved (PBD2)// (%%)
--(% style="color: rgb(0,0,0);" %)DOOCS prop : {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
--(% style="color: rgb(0,0,0);" %)DAQ channel:  {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
--(% style="color: rgb(0,0,0);" %)desc: set pressure in the gas attenuator (%%)
--(% style="color: rgb(0,0,0);" %)units: mbar(%%)
++(% style="color: rgb(255,153,0);" %)(not yet available){{code language="none"}}/FL2/Beamlines/Attenuator/pressure  {{/code}}(%%)
++(% style="color: rgb(193,199,208);" %)//always saved (PBD2)// (%%)
++(% style="color: rgb(193,199,208);" %)DOOCS prop : {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
++(% style="color: rgb(193,199,208);" %)DAQ channel:  {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
++(% style="color: rgb(193,199,208);" %)desc: set pressure in the gas attenuator (%%)
++(% style="color: rgb(193,199,208);" %)units: mbar(%%)
  \\
  {{code language="none"}}/FL2/Beamlines/FL20/Shutter/open{{/code}}
@@ -691,13 +691,11 @@
  units: degree
  \\
--/FL2/Beamlines/Filter wheel/position wheel 2
--always saved (PBD2)
--DOOCS prop : {{code language="none"}}FLASH.FEL/FL20H.PH.MOTOR/MOTOR2.MOT3/FPOS{{/code}}
--DAQ channel: {{code language="none"}}FLASH.FEL/FL20H.PH.MOTOR/MOTOR2.MOT3/FPOS{{/code}}
--desc: Position of the BL filter wheel 2 - to correlate with the filter material please look [[here>>doc:FLASHUSER.jddd-linked help pages.Filter-Units.Filter wheels in FLASH1 and FLASH2.WebHome]]
--units: degree
++//always saved (PBD2)//
++\\\\**[[here>>doc:FLASHUSER.jddd-linked help pages.Filter-Units.Filter wheels in FLASH1 and FLASH2.WebHome]]**
++{{code language="none"}}/FL2/Beamlines/Filter wheel/position wheel 2  DOOCS prop : FLASH.FEL/FL20H.PH.MOTOR/MOTOR2.MOT3/FPOS DAQ channel: FLASH.FEL/FL20H.PH.MOTOR/MOTOR2.MOT3/FPOS desc: Position of the BL filter wheel 2 - to correlate with the filter material please look units: degree{{/code}}
++
  \\
  NOTE: Aperture positions in the beamline as well as the positions of the beam steering mirrors are also saved. for more Info contact your local contact
@@ -776,15 +776,13 @@
  \\
--/FL2/Photon Diagnostic/GMD/Pulse resolved beam position/position tunnel x
--always saved (PBD2)
--DOOCS prop : FLASH.FEL/XGM.BPM/FL2.HALL/X.TD
--DAQ channel: FLASH.FEL/XGM.BPM/FL2.HALL:2
--desc: Besides the well calibrated averaged beam position information there is also the option to measure the beam position on a single bunch level. HOWEVER this methide needs a perfectly adjusted signal level (talk to your local contact !!) and also then the signal to noise is rather small and one needs some averaging ... BUT with this option one can determine if there was a spatial slope on a burst (say forst bunches were lower than the last ones or so ...)
--units : mm
--(x=horizontal, y = vertial)
++//always saved (PBD2)//
++\\\\{{code language="none"}}/FL2/Photon Diagnostic/GMD/Pulse resolved beam position/position tunnel x DOOCS prop : FLASH.FEL/XGM.BPM/FL2.HALL/X.TDDAQ channel: FLASH.FEL/XGM.BPM/FL2.HALL:2 desc: {{/code}}Besides the well calibrated averaged beam position information there is also the option
++{{code language="none"}}to measure the beam position on a single bunch level. HOWEVER this methide needs a perfectly adjusted signal level (talk to your local contact !!) and also then the signal to noise is rather small and one needs some averaging ... BUT with this option one can determine if there was a spatial slope on a burst (say forst bunches were lower than the last ones or so ...)units : mm    {{/code}}
++{{code language="none"}}(x=horizontal, y = vertial){{/code}}
++
  again the same parameter set is available for the **HALL GMD**
  \\
@@ -1185,6 +1185,61 @@
  \\
++=== Samples how to read HDF5 with Matlab ===
++
++The examples apply to HDF files with HDF tree version before vers. 0.3.0.
++
++(% style="color: rgb(0,128,0);" %)%% read in the needed data p=path; path(p,'D:\mess-daten\DAQ-data\hdf5\Gotthard') % add the actual folder to the path hdf5file='FLASH1_EXP-2016-03-16T1420.h5'; % data with Gotthard and VLS
++
++(% style="color: rgb(0,128,0);" %)% h5disp(hdf5file,'/','min') % to get an idea what is in the file % h5info(hdf5file,'/Experiment/Gotthard1/BL.0') % get info about the individual channel
++
++(% style="color: rgb(0,128,0);" %)% read a Number per 10 Hz pulse train:
++
++(% class="code" %)
++(((
++FEL_Wavelength_energy_server=h5read(hdf5file,'/Photon Diagnostic/Wavelength/Calculated by energy/wavelength' );
++)))
++
++(% style="color: rgb(0,128,0);" %)% read in a 1D array (spectrum) %GMD data GMD_Spectrum=h5read(hdf5file,'/Photon Diagnostic/GMD/Pulse resolved energy/energy BDA');
++
++(% style="color: rgb(0,128,0);" %)% This reads in all the data . One can also limit the amount of data read to a subset of the data stored in the file. e.g. for image Data (Gotthard)
++
++\\
++
++(% class="code" %)
++(((
++Start_event= 500 ;   % define the first 10 Hz event
++ Number_of_events = 20; % how many 10 Hz events to load
++)))
++
++(% style="color: rgb(0,128,0);" %)% read only part of the data:
++
++(% class="code" %)
++(((
++Gotthard_data=h5read(hdf5file,'/Experiment/Gotthard1/BL.0',[2   650 Start_event],[50 85 Number_of_events]);
++)))
++
++(% style="color: rgb(0,128,0);" %)% start stop increment "manual" for 2 D data : [start bunch number in Gotthard data start pixelin spectrum start sample in 10 Hz trains], % [number of bunches in the Gotthard data number of points in the spectrum number of spectra ]
++
++(% class="code" %)
++(((
++VLS_Spectrum=h5read(hdf5file,'/Photon Diagnostic/Wavelength/VLS online spectrometer/PCO.ROI.X',[400 Start_event],[200 Number_of_events] );
++
++)))
++
++(% style="color: rgb(0,128,0);" %)% start stop increment "manual" for 1D data : [start sample in the spectrum start sample in 10 Hz trains],[number of points in the spectrum number of spectra ]
++
++(% style="color: rgb(0,128,0);" %)%GMD data
++
++(% class="code" %)
++(((
++GMD_Spectrum=h5read(hdf5file,'/Photon Diagnostic/GMD/Pulse resolved energy/energy BDA',[1 Start_event],[40 Number_of_events] );
++)))
++
++[[Contents>>doc:||anchor="Contents"]]
++
++\\
++
  == HDF5 and DOOCS ==
  Here is an outdated [[list with the available properties that are always saved (PBD) for FLASH1 as>>attach:FLASH1__DaqChannel2HdfNamePbd.xlsx]] HDF5 names and the corresponding DOOCS names

Changes for page The FLASH HDF5 structure

Summary

Details

Applications

Navigation