Changes for page The FLASH HDF5 structure

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

From version 34.1
edited by sndueste
on 2020/10/28 16:14
Change comment: There is no comment for this version
To version 52.1
edited by sndueste
on 2021/04/15 12:08
Change comment: There is no comment for this version

Summary

Details

Page properties
Content
... ... @@ -1,3 +1,6 @@
1 +{{layout}}
2 +{{layout-section ac:type="single"}}
3 +{{layout-cell}}
1 1  == Contents ==
2 2  
3 3  
... ... @@ -8,35 +8,46 @@
8 8  
9 9  == General information about HDF5 ==
10 10  
11 -HDF5 is a data format maintained by the HDF group. For detailed general information see [[here~[~[image:url:http://hasfweb.desy.de/pub/TWiki/TWikiDocGraphics/external-link.gif~|~|width="13" height="12"~]~]>>url:https://www.hdfgroup.org/HDF5/||shape="rect"]]
14 +HDF5 is a data format maintained by the HDF group. For detailed general information see [[here~[~[image:url:http://hasfweb.desy.de/pub/TWiki/TWikiDocGraphics/external-link.gif~|~|width="13" height="12"~]~]>>url:https://portal.hdfgroup.org/display/support||shape="rect"]]
12 12  
13 -To download the HDF5 Viewer click [[here~[~[image:url:http://hasfweb.desy.de/pub/TWiki/TWikiDocGraphics/external-link.gif~|~|width="13" height="12"~]~]>>url:https://support.hdfgroup.org/products/java/hdfview/index.html||shape="rect"]]
16 +To download the HDF5 Viewer click [[here~[~[image:url:http://hasfweb.desy.de/pub/TWiki/TWikiDocGraphics/external-link.gif~|~|width="13" height="12"~]~]>>url:https://portal.hdfgroup.org/display/support/Download+HDFView||shape="rect"]]
14 14  
15 15  For use on the DESY Maxwell or PAL (max-fsc or pal) **hdfview** is available in the xray module :
16 16  {{code language="none"}}> module load xray{{/code}}
17 17  {{code language="none"}}> hdfview{{/code}}
18 18  
22 +or you can use
23 +
24 +{{code language="none"}}> silx view{{/code}}
25 +
26 +\\
27 +
19 19  [[Contents>>doc:||anchor="Contents"]]
20 20  
21 -== FLASH HDF5 structure ==
30 +== The FLASH HDF5 format ==
22 22  
23 -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 organised according to train IDs. The general structure is:
32 +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 gapsusers 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.
24 24  
25 -* Electron Diagnostic
26 -* Photon Diagnostics
27 -* Beamlines
28 -* Experiment
29 -* Timing
30 -[[~[~[image:url:http://hasfweb.desy.de/pub/Setup/HDF5main/HDF5_structure.jpg~|~|alt="HDF5_structure.jpg" width="700" height="500"~]~]>>attach:HDF5_structure.jpg]]
34 +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.
31 31  
32 -A detailled description of (most) channels can be found in the lower part of the hdf5 viewer:
33 -[[~[~[image:url:http://hasfweb.desy.de/pub/Setup/HDF5main/HDF5_structure_desc.jpg~|~|alt="HDF5_structure_desc.jpg" width="341" height="413"~]~]>>attach:HDF5_structure_desc.jpg]]
36 +Reference implementation, which follows the concept of Python libraries like Pandas, Xarray, or Dask, is given below.
34 34  
35 -[[Contents>>doc:||anchor="Contents"]]
38 +\\
36 36  
40 +{{expand title="Discontinued HDF formats"}}
41 +=== Comparison to FLASH's deprecated HDF formats ===
42 +
43 +Before 2021, FLASH provided two different HDF formats formally known as //near-online// and //offline// HDF files.
44 +
45 +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.
46 +
47 +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.
48 +
49 +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.
50 +
37 37  \\
38 38  
39 -== HDF5 example files ==
53 +**HDF5 example files (old format)**
40 40  
41 41  Here we have a few HDF5 samples (User data combined with Photon diagnostics data) from a few beamtimes showing the different kind options.
42 42  
... ... @@ -49,20 +49,67 @@
49 49  \\
50 50  
51 51  \\
66 +{{/expand}}
52 52  
53 -[[Contents>>doc:||anchor="Contents"]]
68 +=== HDF examples: ===
69 +{{/layout-cell}}
70 +{{/layout-section}}
54 54  
72 +{{layout-section ac:type="three_equal"}}
73 +{{layout-cell}}
74 +* ADC data as example for **fast** **data** (10 Hz):           
75 +
55 55  \\
56 56  
78 +[[image:attach:image2020-11-16_15-26-28.png||height="250"]]
79 +{{/layout-cell}}
80 +
81 +{{layout-cell}}
82 +* The //average// FEL pulse energy as example for the **slow** **data** (different e.g. 1 Hz):
83 +
84 + ( every 10th train ID is listed in the HDF group "index")
85 +
57 57  \\
58 58  
88 +[[image:attach:image2020-11-16_15-31-45.png||height="250"]]
89 +{{/layout-cell}}
90 +
91 +{{layout-cell}}
92 +* (((
93 +//zraw// group contains the **original DAQ (DOOCS) names**
94 +)))
95 +
96 + of the properties saved in the DESY internal raw format. (For experts)
97 +
98 +\\
99 +
100 +[[image:attach:image2020-11-16_16-26-3.png||height="400"]]
101 +{{/layout-cell}}
102 +{{/layout-section}}
103 +
104 +{{layout-section ac:type="single"}}
105 +{{layout-cell}}
106 +{{info title="Sample scripts in python"}}
107 +=== Sample scripts / Reference implementation (Python) ===
108 +
109 +[[~[~[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"]]
110 +
111 +(% class="Object" %)[[https:~~/~~/gitlab.desy.de/christopher.passow/flash-daq-hdf>>url:https://gitlab.desy.de/christopher.passow/flash-daq-hdf||shape="rect"]]
112 +{{/info}}
113 +{{/layout-cell}}
114 +{{/layout-section}}
115 +
116 +{{layout-section ac:type="single"}}
117 +{{layout-cell}}
118 +== Complete list of recordable parameters ==
119 +
120 +The complete list for the relation between DOOCS names and HDF5 names for the recordable parameters can be found in [[DESY's Repository~[~[image:url:http://hasfweb.desy.de/pub/TWiki/TWikiDocGraphics/external-link.gif~|~|width="13" height="12"~]~]>>url:https://stash.desy.de/projects/CS/repos/pah/browse/src/camp/data/channel2HdfName.dat||shape="rect"]].
121 +\\
122 +
59 59  == Most popular FLASH parameters and their names in HDF5, DOOCS and (raw) DAQ ==
60 60  
61 61  {{id name="DOOCSparameters"/}}
62 62  
63 -The complete list for the relation between DOOCS names and HDF5 names for the recordable parameters can be found in [[DESY's Repository~[~[image:url:http://hasfweb.desy.de/pub/TWiki/TWikiDocGraphics/external-link.gif~|~|width="13" height="12"~]~]>>url:https://stash.desy.de/projects/CS/repos/pah/browse/src/camp/data/channel2HdfName.dat||shape="rect"]].
64 -The most common and often used ones are summarized below:
65 -
66 66  Note, the HDF group and data set names apply to our HDF tree version since vers. 0.3.0.
67 67  
68 68  \\
... ... @@ -69,8 +69,6 @@
69 69  
70 70  === FLASH1 ===
71 71  
72 -\\
73 -
74 74  ==== Beamline info (FLASH1) ====
75 75  
76 76  {{code language="none"}}/FL1/Beamlines/Attenuator/pressure{{/code}}
... ... @@ -145,12 +145,14 @@
145 145  
146 146  \\
147 147  
148 -==== Photon Diagnostics SASE ([[GMD>>url:http://photon-science.desy.de/facilities/flash/photon_diagnostics/gmd_intensity_and_position/index_eng.html||shape="rect"]]) / Spectrometer (FLASH1) ====
207 +==== Photon Diagnostics SASE ([[GMD>>url:http://photon-science.desy.de/facilities/flash/photon_diagnostics/gmd_intensity_and_position/index_eng.html||shape="rect"]]) ====
149 149  
209 +\\
210 +
211 +(% style="color: rgb(0,0,0);" %)**Discontinued GMD data recording / evaluation  (VME + PhotonFlux ML server)**
212 +
150 150  {{code language="none"}}/FL1/Photon Diagnostic/GMD/Average energy/energy tunnel{{/code}}
151 -//always saved (PBD)//
152 -(% style="color: rgb(165,173,186);" %)New MTCA Version: DOOCS prop : {{code language="none"}}FLASH.FEL/XGM.PHOTONFLUX/FL1.TUNNEL/PHOTONFLUX.UJ{{/code}} (%%)
153 -(% style="color: rgb(165,173,186);" %)New MTCA Version: DAQ channel: {{code language="none"}}FLASH.FEL/XGM.PHOTONFLUX/FL1.TUNNEL/PHOTONFLUX.UJ{{/code}}
214 +//always saved (PBD)//
154 154  
155 155  (% style="color: rgb(0,0,0);" %)DOOCS prop : {{code language="none"}}TTF2.DAQ/PHFLUX/OUT04/VAL{{/code}} (%%)
156 156  (% style="color: rgb(0,0,0);" %)DAQ channel: {{code language="none"}}PBD.PHFLUX/TUNNEL.ENPULSEIC{{/code}}(%%)
... ... @@ -160,9 +160,6 @@
160 160  {{code language="none"}}/FL1/Photon Diagnostic/GMD/Pulse resolved energy/energy tunnel{{/code}}
161 161  //always saved (PBD)//
162 162  
163 -(% style="color: rgb(193,199,208);" %)New MTCA Version: DOOCS prop : 
164 -{{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL1.TUNNEL/INTENSITY.TD{{/code}}New MTCA Version: DAQ channel: {{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL1.TUNNEL/INTENSITY.TD{{/code}}
165 -
166 166  DOOCS prop : {{code language="none"}}TTF2.DAQ/PHFLUX/OUT34/VAL{{/code}}
167 167  DAQ channel: {{code language="none"}}PBD.PHFLUX/TUNNEL.ENERGYPULSE.USER{{/code}}
168 168  desc :Energy per pulse Tunnel (from e-) - the values are set to "0" if there was no SASE beam in the FEL
... ... @@ -178,9 +178,6 @@
178 178  {{code language="none"}}/FL1/Photon Diagnostic/GMD/Average energy/energy BDA{{/code}}
179 179  //always saved (PBD)//
180 180  
181 -(% style="color: rgb(165,173,186);" %)New MTCA Version: DOOCS prop : {{code language="none"}}FLASH.FEL/XGM.PHOTONFLUX/FL1.BDA/PHOTONFLUX.UJ{{/code}} (%%)
182 -(% style="color: rgb(165,173,186);" %)New MTCA Version: DAQ channel: {{code language="none"}}FLASH.FEL/XGM.PHOTONFLUX/FL1.BDA/PHOTONFLUX.UJ{{/code}}
183 -
184 184  DOOCS prop : {{code language="none"}}TTF2.DAQ/PHFLUX/OUT05/VAL{{/code}}
185 185  DAQ channel: {{code language="none"}}PBD.PHFLUX/BDA.ENPULSEIC{{/code}}
186 186  desc : calibrated average SASE Energy/pulse measured in the BDA (in the experimental hall) after the attenuator (ion current)
... ... @@ -189,8 +189,6 @@
189 189  {{code language="none"}}/FL1/Photon Diagnostic/GMD/Pulse resolved energy/energy BDA{{/code}}
190 190  //always saved (PBD)//
191 191  
192 -(% style="color: rgb(193,199,208);" %)New MTCA Version: DOOCS prop : 
193 -{{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL1.BDA/INTENSITY.TD{{/code}}New MTCA Version: DAQ channel: {{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL1.BDA/INTENSITY.TD{{/code}}(%%)
194 194  DOOCS prop : {{code language="none"}}TTF2.DAQ/PHFLUX/OUT35/VAL{{/code}}
195 195  DAQ channel: {{code language="none"}}PBD.PHFLUX/BDA.ENERGYPULSE.USER{{/code}}
196 196  desc :Energy per pulse BDA (from e-) - the values are set to "0" if there was no SASE beam in the FEL
... ... @@ -205,6 +205,55 @@
205 205  
206 206  \\
207 207  
261 +(% style="color: rgb(0,0,0);" %)**NEW (2021) GMD data recording / evaluation  (MTCA, analog to FLASH2 and XFEL)**
262 +
263 +{{code language="none"}}/FL1/Photon Diagnostic/GMD/Average energy/energy tunnel{{/code}}
264 +//always saved (PBD)//
265 +(% style="color: rgb(0,0,0);" %)DOOCS prop : {{code language="none"}}FLASH.FEL/XGM.PHOTONFLUX/FL1.TUNNEL/PHOTONFLUX.UJ{{/code}} (%%)
266 +(% style="color: rgb(0,0,0);" %)DAQ channel: {{code language="none"}}FLASH.FEL/XGM.PHOTONFLUX/FL1.TUNNEL/PHOTONFLUX.UJ{{/code}}
267 +
268 +desc : calibrated average SASE Energy/pulse measured in the TUNNEL upstream the gas attenuator
269 +units : microJ
270 +
271 +\\
272 +
273 +{{code language="none"}}/FL1/Photon Diagnostic/GMD/Pulse resolved energy/energy tunnel{{/code}}
274 +//always saved (PBD)//
275 +
276 +(% style="color: rgb(0,0,0);" %)DOOCS prop : 
277 +{{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL1.TUNNEL/INTENSITY.TD{{/code}}DAQ channel: {{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL1.TUNNEL/INTENSITY.TD{{/code}}
278 +
279 +desc : Energy per pulse measured in the Tunnel.  In addition measurement errors and beam position are included for EACH pulse in the pulse train  !! (see below)
280 +units : a.u. (more or less µJ but need to be calibrated with the "Average energy" for good precision)** [[see here for help>>doc:FLASHUSER.jddd-linked help pages.Calibrating the pulse resolved (electron) data from GMD.WebHome]]**
281 +
282 +{{info title="GMD pulse resolved data structure"}}
283 +For every pulse in the pulse train the information is saved:
284 +
285 +
286 +1. **Intensity per pulse** (a.u. (more or less µJ ))
287 +1. Intensity per pulse (auxillary GMD) - not used
288 +1. Position horizontal (mm, for a single pulse the position information may be very noisy  - talk to your local contact)
289 +1. Position vertical (mm, for a single pulse the position information may be very noisy  - talk to your local contact)
290 +1. **Intensity per pulse sigma** (a.u. (more or less µJ ), (% style="color: rgb(0,0,0);" %)This parameter gives an indication of the error of the measurement of the pulse energy. This takes signal to noise,  detector resolution, uncertainties in crossection etc into account. (it is NOT the measurement of the statistical fluctuation of the SASE pulses)(%%))
291 +1. Position horizontal sigma  (mm, indicates the error (RMS, sigma) of the measurement according to known uncertainties and signal to noise)
292 +1. Position vertical sigma  (mm, indicates the error (RMS, sigma) of the measurement according to known uncertainties and signal to noise)
293 +1. Combined warning and error flags
294 +
295 +The pulse energy and the error are plotted for the first bunch of the pulse trains saved in this 
296 +file
297 +
298 +
299 +[[image:attach:image2021-2-9_10-51-6.png||height="250"]]
300 +
301 +\\
302 +{{/info}}
303 +
304 +All values for the GMD are also available for the** BDA GMD **which is located in the experimental hall down stream the gas attenuator. If the attenuator is on the ratio between BDA and Tunnel signal shows the attenuation. **BUT NOTE** that the filter units  are downstream the GMD. So if filters  are used this influence is NOT measured by the GMD BDA !
305 +
306 +\\
307 +
308 +Besides  pulse energy the GMD also provides information about the beam position
309 +
208 208  {{code language="none"}}/FL1/Photon Diagnostic/GMD/Beam position/position BDA horizontal{{/code}}
209 209  //always saved (PBD)//
210 210  DOOCS prop : {{code language="none"}}FLASH.FEL/XGM.POSMON/FL1.BDA/IX.POS{{/code}}
... ... @@ -233,6 +233,10 @@
233 233  desc :Beam position of the photon Beam determined by the GMD (TUNNEL, y=vertical)
234 234  units : mm
235 235  
338 +\\
339 +
340 +==== Spectrometer (FLASH1) ====
341 +
236 236  {{code language="none"}}/FL1/Photon Diagnostic/Wavelength/Tunnelspectrometer/wavelength{{/code}}
237 237  _always saved (when Spectrum is measured !!) (PBD) _
238 238  DOOCS prop : {{code language="none"}}TTF2.EXP/PHOTONWL.ML/WAVE_LENGTH/VAL.TD{{/code}}
... ... @@ -363,7 +363,7 @@
363 363  ===== electron bunch energy =====
364 364  
365 365  {{code language="none"}}/FL1/Electron Diagnostic/Electron energy/average electron energy{{/code}}
366 -_always saved (PBD)
472 +//always saved (PBD)//
367 367  DOOCS prop : {{code language="none"}}TTF2.DAQ/ENERGY.DOGLEG/E_INTRA_MEAN/VAL{{/code}}
368 368  DAQ channel: {{code language="none"}}PBD.ENERGY.DOGLEG/E_MEAN{{/code}}
369 369  desc: electron bunch energy (average over the bunch train)
... ... @@ -371,7 +371,7 @@
371 371  \\
372 372  
373 373  {{code language="none"}}/FL1/Electron Diagnostic/Electron energy/pulse resolved energy{{/code}}
374 -_always saved (PBD)
480 +//always saved (PBD)//
375 375  DOOCS prop : {{code language="none"}}TTF2.DAQ/ENERGY.DOGLEG/E_SPECT/VAL.TD{{/code}}
376 376  DAQ channel: {{code language="none"}}PBD.ENERGY.DOGLEG/E_SPECT{{/code}}
377 377  desc: electron bunch energy bunch resolved
... ... @@ -380,7 +380,7 @@
380 380  \\
381 381  
382 382  {{code language="none"}}/FL1/Electron Diagnostic/Electron energy/wavelength bunch train average{{/code}}
383 -_always saved (PBD)
489 +//always saved (PBD)//
384 384  DOOCS prop : {{code language="none"}}TTF2.DAQ/ENERGY.DOGLEG/LAMBDA_MEAN/VAL{{/code}}
385 385  DAQ channel: {{code language="none"}}PBD.ENERGY.DOGLEG/LAMBDA_MEAN{{/code}}
386 386  desc: Wavelength calculated by the electron bunch energy (average over the bunch train) (FLASH1)
... ... @@ -428,7 +428,7 @@
428 428  //always saved (PBD)//
429 429  DOOCS prop : {{code language="none"}}TTF2.DIAG/PBD.TOROID.ML/12EXP/CHARGE.TD{{/code}}
430 430  DAQ channel: {{code language="none"}}TTF2.DIAG/PBD.TOROID.ML/12EXP{{/code}}
431 -desc: The bunch pattern as function of time in a burst recorded by toroide diagnostic BEHIND the undulator. (FLASH1)
537 +desc: The bunch pattern as function of time in a burst recorded by toroid diagnostic BEHIND the undulator. (FLASH1)
432 432  units:
433 433  \\
434 434  
... ... @@ -501,7 +501,7 @@
501 501  //always saved (PBD)//
502 502  DOOCS prop : {{code language="none"}}TTF2.FEL/TDOLFEL/TDOLFEL/STREAK.CAM.TIME{{/code}}
503 503  DOOCS prop : {{code language="none"}}TTF2.FEL/TDOLFEL/TDOLFEL/STREAK.CAM.TIME{{/code}}
504 -desc: delaytime between the optical laser and the FEL units: ps
610 +desc: delay time between the optical laser and the FEL units: ps
505 505  
506 506  \\
507 507  
... ... @@ -511,7 +511,7 @@
511 511  
512 512  ==== User Data (FLASH1) ====
513 513  
514 -The data saved specifically for detectors at an experiment will show up in /Experiment/ there is a large number of options for cameras or monitoring pslow properties (motor positons etc) for user experiments. For details please ask your local contact.
620 +The data saved specifically for detectors at an experiment will show up in /Experiment/ there is a large number of options for cameras or monitoring of slow properties (motor positions etc) for user experiments. For details please ask your local contact.
515 515  
516 516  NOTE: If parameters for an experiment are included on short notice the correct naming in the HDF5 may not be in time and the data will show up in /uncategorized/ with the DOOCS names
517 517  
... ... @@ -539,13 +539,13 @@
539 539  {{code language="none"}}/FL1/Experiment/BL3/ADQ412 GHz ADC/CH03/TD{{/code}}
540 540  
541 541  DOOCS prop : {{code language="none"}}FLASH.FEL/ADC.ADQ.PG/EXP1.CH00/CH00.TD  or CH00.DAQ.TD{{/code}}
542 -here the {{code language="none"}}CH00.TD{{/code}} is the full ADC trace as it is sampled ( typically several 100.000 samples per pulse train) while the {{code language="none"}}CH00.DAQ.TD{{/code}} trace only has the number of samples which are sent to the DAQ OR if //grouping// is activated the {{code language="none"}}CH00.DAQ.TD{{/code}} conatins only the grouped spectra. To read the ADC trace with an online analysis program the {{code language="none"}}CH00.DAQ.TD{{/code}} is used preferablly.
648 +here the {{code language="none"}}CH00.TD{{/code}} is the full ADC trace as it is sampled ( typically several 100.000 samples per pulse train) while the {{code language="none"}}CH00.DAQ.TD{{/code}} trace only has the number of samples which are sent to the DAQ OR if //grouping// is activated the {{code language="none"}}CH00.DAQ.TD{{/code}} contains only the grouped spectra. To read the ADC trace with an online analysis program the {{code language="none"}}CH00.DAQ.TD{{/code}} is used preferably.
543 543  DAQ channel: {{code language="none"}}FLASH.FEL/ADC.ADQ.PG/EXP1.CH00{{/code}}
544 544  
545 545  In addition there are also additional parameters saved like:
546 546  
547 547  * {{code language="none"}}sample frequency{{/code}}: it shows the sample frequency in MHz (number of samples per µs). NOTE: the clock of the ADC is NOT synchronized to the FLASH timing system. Thus the number of samples between bunches in the bunch train may be not integer numbers which will be show up for long bunch trains.
548 -* {{code language="none"}}number of samples{{/code}}: total number of samoles recorded for each 10 Hz trigger
654 +* {{code language="none"}}number of samples{{/code}}: total number of samples recorded for each 10 Hz trigger
549 549  * {{code language="none"}}error (ADC):{{/code}} 0 indicates that there was no error
550 550  
551 551  ===== MHz ADCs =====
... ... @@ -558,7 +558,7 @@
558 558  In addition there are also additional parameters saved like:
559 559  
560 560  * {{code language="none"}}sample frequency{{/code}}: it shows the sample frequency in MHz (number of samples per µs). NOTE: the clock of the ADC is NOT synchronized to the FLASH timing system. Thus the number of samples between bunches in the bunch train may be not integer numbers which will be show up for long bunch trains.
561 -* {{code language="none"}}number of samples{{/code}}: total number of samoles recorded for each 10 Hz trigger
667 +* {{code language="none"}}number of samples{{/code}}: total number of samples recorded for each 10 Hz trigger
562 562  
563 563  [[Contents>>doc:||anchor="Contents"]]
564 564  
... ... @@ -574,12 +574,12 @@
574 574  
575 575  ==== Beamline info (FLASH2) ====
576 576  
577 -(% style="color: rgb(255,153,0);" %)(not yet available){{code language="none"}}/FL2/Beamlines/Attenuator/pressure  {{/code}}(%%)
578 -(% style="color: rgb(193,199,208);" %)//always saved (PBD2)// (%%)
579 -(% style="color: rgb(193,199,208);" %)DOOCS prop : {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
580 -(% style="color: rgb(193,199,208);" %)DAQ channel:  {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
581 -(% style="color: rgb(193,199,208);" %)desc: set pressure in the gas attenuator (%%)
582 -(% style="color: rgb(193,199,208);" %)units: mbar(%%)
683 +{{code language="none"}}/FL2/Beamlines/Attenuator/pressure  {{/code}}
684 +(% style="color: rgb(0,0,0);" %)//always saved (PBD2)// (%%)
685 +(% style="color: rgb(0,0,0);" %)DOOCS prop : {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
686 +(% style="color: rgb(0,0,0);" %)DAQ channel:  {{code language="none"}}FLASH.FEL/ATT.GAS_DOSING/FL2.HALL/PRESSURE{{/code}} (%%)
687 +(% style="color: rgb(0,0,0);" %)desc: set pressure in the gas attenuator (%%)
688 +(% style="color: rgb(0,0,0);" %)units: mbar(%%)
583 583  \\
584 584  
585 585  {{code language="none"}}/FL2/Beamlines/FL20/Shutter/open{{/code}}
... ... @@ -598,11 +598,13 @@
598 598  units: degree
599 599  \\
600 600  
707 +/FL2/Beamlines/Filter wheel/position wheel 2
708 +always saved (PBD2)
709 +DOOCS prop : {{code language="none"}}FLASH.FEL/FL20H.PH.MOTOR/MOTOR2.MOT3/FPOS{{/code}}
710 +DAQ channel: {{code language="none"}}FLASH.FEL/FL20H.PH.MOTOR/MOTOR2.MOT3/FPOS{{/code}}
711 +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]]
712 +units: degree
601 601  
602 -//always saved (PBD2)//
603 -\\\\**[[here>>doc:FLASHUSER.jddd-linked help pages.Filter-Units.Filter wheels in FLASH1 and FLASH2.WebHome]]**
604 -{{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}}
605 -
606 606  \\
607 607  
608 608  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
... ... @@ -624,40 +624,37 @@
624 624  
625 625  \\
626 626  
627 -(% style="color: rgb(255,204,0);" %)average Sigma    to be included in
628 -
629 -(% style="color: rgb(193,199,208);" %){{code language="none"}}/FL2/Photon Diagnostic/GMD/Average energy/energy tunnel uncertainty{{/code}} (%%)
630 -(% style="color: rgb(193,199,208);" %)//always saved (PBD2)// (%%)
631 -(% style="color: rgb(193,199,208);" %)DOOCS prop : FLASH.FEL/XGM.PHOTONFLUX/FL2.TUNNEL/PHOTONFLUX.UJ.SIGMA (%%)
632 -(% style="color: rgb(193,199,208);" %)DAQ channel: FLASH.FEL/XGM.PHOTONFLUX/FL2.TUNNEL/PHOTONFLUX.UJ.SIGMA (%%)
633 -(% style="color: rgb(193,199,208);" %)desc :  This parameter gives an indication of the error of the measurement of the average pulse energy. This takes signal to noise,  detector resolution, uncertainties in crossection etc into account. (it is NOT the measurement of the statistical fluctuation of the SASE pulses)(%%)
634 -(% style="color: rgb(193,199,208);" %)units : microJ
635 -
636 -\\
637 -
638 -\\
639 -
640 640  {{code language="none"}}/FL2/Photon Diagnostic/GMD/Pulse resolved energy/energy tunnel{{/code}}
641 641  //always saved (PBD2)//
642 642  DOOCS prop : {{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL2.TUNNEL/INTENSITY.TD{{/code}}
643 643  DAQ channel: {{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL2.TUNNEL/INTENSITY.TD{{/code}}
644 -desc : Energy per pulse measured in the Tunnel  (in fromt of the gas attenuator and the apertures in the Hall)
739 +desc : Energy per pulse measured in the Tunnel  (in front of the gas attenuator and the apertures in the Hall) In addition measurement errors and beam position are included for EACH pulse in the pulse train  !! (see below)
645 645  units : a.u. (more or less µJ but need to be calibrated with the "Average energy" for good precision)** [[see here for help>>doc:FLASHUSER.jddd-linked help pages.Calibrating the pulse resolved (electron) data from GMD.WebHome]]**
646 646  
647 -\\
742 +{{info title="GMD pulse resolved data structure"}}
743 +For every pulse in the pulse train the information is saved:
744 +
648 648  
649 -\\
746 +1. **Intensity per pulse** (a.u. (more or less µJ ))
747 +1. Intensity per pulse (auxillary GMD) - not used
748 +1. Position horizontal (mm, for a single pulse the position information may be very noisy  - talk to your local contact)
749 +1. Position vertical (mm, for a single pulse the position information may be very noisy  - talk to your local contact)
750 +1. **Intensity per pulse sigma** (a.u. (more or less µJ ), (% style="color: rgb(0,0,0);" %)This parameter gives an indication of the error of the measurement of the pulse energy. This takes signal to noise,  detector resolution, uncertainties in crossection etc into account. (it is NOT the measurement of the statistical fluctuation of the SASE pulses)(%%))
751 +1. Position horizontal sigma  (mm, indicates the error (RMS, sigma) of the measurement according to known uncertainties and signal to noise)
752 +1. Position vertical sigma  (mm, indicates the error (RMS, sigma) of the measurement according to known uncertainties and signal to noise)
753 +1. Combined warning and error flags
650 650  
651 -{{code language="none"}}/FL2/Photon Diagnostic/GMD/Pulse resolved energy/energy tunnel uncertainty (sigma){{/code}}
652 -//always saved (PBD2)//
653 -DOOCS prop : {{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL2.TUNNEL/INTENSITY.SIGMA.TD{{/code}}
654 -DAQ channel: {{code language="none"}}FLASH.FEL/XGM.INTENSITY/FL2.TUNNEL:4{{/code}}
655 -desc :(% style="color: rgb(0,0,0);" %)This parameter gives an indication of the error of the measurement of the pulse energy. This takes signal to noise,  detector resolution, uncertainties in crossection etc into account. (it is NOT the measurement of the statistical fluctuation of the SASE pulses)(%%)
656 -units : a.u. (more or less µJ - related to the value of the pulse energy )
755 +The pulse energy and the error are plotted for the first bunch of the pulse trains saved in this 
756 +file
657 657  
758 +
759 +[[image:attach:image2021-2-9_10-51-6.png||height="250"]]
658 658  
659 -{{code language="none"}}{{/code}}
761 +\\
762 +{{/info}}
660 660  
764 +\\
765 +
661 661  All values for the GMD are also available for the** HALL GMD **which is located in the experimental hall down stream the gas attenuator. If the attenuator is on the ratio between Hall and Tunnel signal shows the attenuation. **BUT NOTE** that the filter unit and the Aperture 4 are downstream the GMD. So if filters and aperture are used this influence is NOT measured by the GMD hall !
662 662  
663 663  \\
... ... @@ -684,13 +684,15 @@
684 684  
685 685  \\
686 686  
792 +/FL2/Photon Diagnostic/GMD/Pulse resolved beam position/position tunnel x
793 +always saved (PBD2)
794 +DOOCS prop : FLASH.FEL/XGM.BPM/FL2.HALL/X.TD
795 +DAQ channel: FLASH.FEL/XGM.BPM/FL2.HALL:2
796 +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 ...)
797 +units : mm   
687 687  
688 -//always saved (PBD2)//
689 -\\\\{{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
690 -{{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}}
799 +(x=horizontal, y = vertial)
691 691  
692 -{{code language="none"}}(x=horizontal, y = vertial){{/code}}
693 -
694 694  again the same parameter set is available for the **HALL GMD**
695 695  
696 696  \\
... ... @@ -718,7 +718,7 @@
718 718  // saved opon request (PBD2)//
719 719  DOOCS prop : {{code language="none"}}FLASH.UTIL/STORE/FL2.TUNNEL.OPIS/VAL040{{/code}}
720 720  DAQ channel:{{code language="none"}} FLASH.UTIL/STORE/FL2.TUNNEL.OPIS/VAL040{{/code}}
721 -desc : meanwavelength ( ~~ 1 sec averaging time ) measured in the TUNNEL for a specific bunch out of the bunch train (via photoelectron spectroscopy)
828 +desc : mean wavelength ( ~~ 1 sec averaging time ) measured in the TUNNEL for a specific bunch out of the bunch train (via photoelectron spectroscopy)
722 722  units : nm
723 723  
724 724  \\
... ... @@ -732,11 +732,11 @@
732 732  
733 733  \\
734 734  
735 -If Opis is running typically on the the averaged data is saved. For several experiments it may make sense to save the information for each single bunch. This is up to now done by savng the compleate ADC trace of the TOF setup. This is a huge amount of data and needs processing. This has to be performed after the beamtime in close contact to [[Markus Braune>>mailto:markus.braune@desy.de||shape="rect"]] ( respobsible for [[OPIS>>url:http://photon-science.desy.de/facilities/flash/photon_diagnostics/opis_spectrometer/index_eng.html||shape="rect"]])
842 +If Opis is running typically on the the averaged data is saved. For several experiments it may make sense to save the information for each single bunch. This is up to now done by saving the complete ADC trace of the TOF setup. This is a huge amount of data and needs processing. This has to be performed after the beamtime in close contact to [[Markus Braune>>mailto:markus.braune@desy.de||shape="rect"]] ( responsible for [[OPIS>>url:http://photon-science.desy.de/facilities/flash/photon_diagnostics/opis_spectrometer/index_eng.html||shape="rect"]])
736 736  
737 737  \\
738 738  
739 -In case OPIS was not operating there is still informaton about the **set wavelength** for the undulators (see below) which may differ by up to 5 % from the actual wavelength due to different settings in the FEL ...
846 +In case OPIS was not operating there is still information about the **set wavelength** for the undulators (see below) which may differ by up to 5 % from the actual wavelength due to different settings in the FEL ...
740 740  
741 741  \\
742 742  
... ... @@ -755,8 +755,21 @@
755 755  
756 756  \\
757 757  
758 -===== undulator settings =====
865 +===== electron bunch energy =====
759 759  
867 +{{code language="none"}}/FL2/Electron Diagnostic/Electron energy/energy of first bunch/behind undulators{{/code}}
868 +//always saved (PBD2)//
869 +DOOCS prop : {{code language="none"}}FLASH.DIAG/BEAM_ENERGY_MEASUREMENT/FL2XTDS/ENERGY.FLASH2{{/code}}
870 +DAQ channel: (% style="color: rgb(94,108,132);" %)TTF2.DAQ/PBD2.BEAM.ENERGY.MEAS.ML.COPY/FL2XTDS.ENERGY.FLASH2{{code language="none"}}{{/code}}(%%)
871 +desc: electron bunch energy measured behind the undulator. Data is saved with 10 Hz - BUT (for computation reasons) only the energy of the FIRST bunch is recorded. The data is also available for (% style="color: rgb(94,108,132);" %)extraction and septum  in the beginning of FLASH2
872 +
873 +(% style="letter-spacing: 0.0px;" %)units: (% class="twikiNewLink" %)MeV
874 +
875 +\\
876 +
877 +(% style="color: rgb(94,108,132);font-weight: 600;letter-spacing: 0.0px;" %)
878 +undulator settings
879 +
760 760  {{code language="none"}}/FL2/Electron Diagnostic/Undulator setting/set wavelength{{/code}}
761 761  //always saved (PBD2)//
762 762  DOOCS prop : {{code language="none"}}TTF2.FEEDBACK/FL2.WAVELENGTHCONTROL/FLASH2/WAVELENGTH{{/code}}
... ... @@ -1086,72 +1086,6 @@
1086 1086  [[Contents>>doc:||anchor="Contents"]]
1087 1087  
1088 1088  \\
1089 -
1090 -== Example code showing how to access HDF5 files ==
1091 -
1092 -\\
1093 -
1094 -=== Samples how to read HDF5 with Matlab ===
1095 -
1096 -The examples apply to HDF files with HDF tree version before vers. 0.3.0.
1097 -
1098 -(% 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
1099 -
1100 -(% 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
1101 -
1102 -(% style="color: rgb(0,128,0);" %)% read a Number per 10 Hz pulse train:
1103 -
1104 -(% class="code" %)
1105 -(((
1106 -FEL_Wavelength_energy_server=h5read(hdf5file,'/Photon Diagnostic/Wavelength/Calculated by energy/wavelength' );
1107 -)))
1108 -
1109 -(% 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');
1110 -
1111 -(% 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)
1112 -
1113 -\\
1114 -
1115 -(% class="code" %)
1116 -(((
1117 -Start_event= 500 ; % define the first 10 Hz event
1118 - Number_of_events = 20; % how many 10 Hz events to load
1119 -)))
1120 -
1121 -(% style="color: rgb(0,128,0);" %)% read only part of the data:
1122 -
1123 -(% class="code" %)
1124 -(((
1125 -Gotthard_data=h5read(hdf5file,'/Experiment/Gotthard1/BL.0',[2 650 Start_event],[50 85 Number_of_events]);
1126 -)))
1127 -
1128 -(% 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 ]
1129 -
1130 -(% class="code" %)
1131 -(((
1132 -VLS_Spectrum=h5read(hdf5file,'/Photon Diagnostic/Wavelength/VLS online spectrometer/PCO.ROI.X',[400 Start_event],[200 Number_of_events] );
1133 -
1134 -)))
1135 -
1136 -(% 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 ]
1137 -
1138 -(% style="color: rgb(0,128,0);" %)%GMD data
1139 -
1140 -(% class="code" %)
1141 -(((
1142 -GMD_Spectrum=h5read(hdf5file,'/Photon Diagnostic/GMD/Pulse resolved energy/energy BDA',[1 Start_event],[40 Number_of_events] );
1143 -)))
1144 -
1145 -[[Contents>>doc:||anchor="Contents"]]
1146 -
1147 -=== How to read HDF5 with Python via FLASHH5 ===
1148 -
1149 -With [[FlashH5>>doc:FLASHUSER.FLASHH5 for easier access of FLASH's HDF5 files]] we provide h5py based examples to access and evaluate FLASH's HDF5 files.
1150 -
1151 -\\
1152 -
1153 -== HDF5 and DOOCS ==
1154 -
1155 -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
1156 -
1157 -[[Contents>>doc:||anchor="Contents"]]
1209 +{{/layout-cell}}
1210 +{{/layout-section}}
1211 +{{/layout}}