photon-diagnostics/wiki
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Update PBPS data

Browse Source
This commit is contained in:
arrell_c 2025-01-27 14:44:29 +01:00
parent d325547bc4
commit 93998bac21
1 changed files with 11 additions and 46 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

57
Aramis-Intensity-Diagnostics/PBPS-data.md
View File

@ -1,46 +1,11 @@
## Channel list ## IOC naming
All PSSS data channels are listed in the following sheet: | **PBPS Device** | **IOC Name** |
|--------------------|---------------------------|
[Sheet of PSSS data channels](https://docs.google.com/spreadsheets/d/16Hz7MZR3E7-OmLHfakkT2K6yPMFUzp4eIFZxrzXWjFk/edit?usp=sharing) | SARFE10-PBPS053 | SARFE10-CVME-PHO6212 |
| SAROP11-PBPS110 | SAROP11-CVME-PBPS2 |
## Nomenclature | SAROP11-PBPS122 | SAROP11-CVME-PBPS1 |
| SAROP21-PBPS103 | SAROP21-CVME-PBPS1 |
| Style | Meaning| | SAROP21-PBPS133 | SAROP21-CVME-PBPS3 |
|----------|----------| | SARES21-PBPS141 | SAROP21-CVME-PBPS2 |
| SARFE10-PSSS059:**SPECTRUM***|Raw spectral data | | SAROP31-PBPS113 | SAROP31-CVME-PBPS113 |
| SARFE10-PSSS059:**FIT*** | Results of Gaussian fit | | SAROP31-PBPS149 | SAROP31-CVME-PBPS149 |
| SARFE10-PSSS059:**SPECT*** | Statistics from the raw (without fit) spectrum |
| SARFE10-PSSS059:**AVG*** | Data calculated from average (not single shot) spectra |
## Concepts
**Width and center**
Spectral width and central energy are calculated by both a Gaussian fit and directly from the spectra. Depending on the spectral intensity profile, the Gaussian may return a poor representation of the spectra, for example for a spectral with a long tail. Before both fitting and calculation the spiky single shot spectra are smoothed with a savgol filter:
```
# smooth the spectrum with savgol filter with 51 window size and 3rd order polynomial
smoothed_spectrum = scipy.signal.savgol_filter(spectrum, 51, 3)
```
The standard deviation and center of mass for the spectra are calculated as:
```
smoothed_spectrum_normed = smoothed_spectrum / np.sum(smoothed_spectrum)
spectrum_com = np.sum(axis * smoothed_spectrum_normed)
spectrum_std = np.sqrt(np.sum((axis - spectrum_com) ** 2 * smoothed_spectrum_normed))
```
**Relative energy spread**
Channels ending in **RES** are a measure of the relative energy spread in units of ‰. For **FIT** data the relative energy spread is calculated as:
$`\frac{2.355*\sigma}{\mu}`$
and for directly calculated values (**SPECT**):
$`\frac{2.355*RMS}{\text{Center of mass}}`$
where the factor or 2.355 is included to allow these values to be comparable.
## Data flow schematic
<img src="https://docs.google.com/drawings/d/e/2PACX-1vRShOtDRMpQTAY1p5oZKvRcuyLToMvmgzTQHOS_JrdQArvqcul-RoMcIobmKFc8RPgJVjGzULKczS5-/pub?w=1440&amp;h=1080">