Modify Radius_StartNr_1_EndNr_1440.csv to report weight-averaged values of SigmaR and SigmaEta

[StopkaKris]

It would be ideal if Radius_StartNr_1_EndNr_1440.csv and Result_StartNr_1_EndNr_1440.csv contained the weighted averages of the Gaussian and Lorentzian fits for the radial (SigmaR) and azimuthal (SigmaEta) spot widths, instead of just the maximum that is reported right now. It is unclear how this data can be currently extracted for spots that span multiple frames in omega, since these are "aggregated" when moving from the *_PS.csv files in the Temp folder to Radius_StartNr_1_EndNr_1440.csv. Is this something we can implement please? It looks like maybe just a change to CalcRadiusAllZarr.c?

Furthermore, in the individual *_PS.csv files in the Temp folder, the values in the SigmaEta (and SigmaR) column are just the maximum between SigmaGEta and SigmaLEta. Since this file provides Mu (the ratio between the two fits), the values in the SigmaEta column might as well just be the actual computed width (weight averaged) using SigmaGEta, SigmaLEta, and Mu. And the user can just grab the maximum of the first two values if needed. What do you think?

[marinerhemant]

Kris, I have something in the works for this, stay tuned.

[jjsolano21]

Hey Hemant, I wanted to take a crack at this over the weekend. I edited MergeOverlappingPeaksAllZarr.c to track what spots are combined across diffraction frames - it's attached to this reply. My version creates a map 'Count_Combos.csv' to be used in the post-processing MATLAB script 'MergeWithMap.m' - also attached to this reply.

I've double checked that this map is capable of re-creating 'Result_StartNr_1_EndNr_1440.csv'. The post-processing MATLAB script also appends 8 values to 'Result_StartNr_1_EndNr_1440.csv':

1. BG for the spot that yielded the largest SigmaR

2. SigmaGR for the spot that yielded the largest SigmaR

3. SigmaLR for the spot that yielded the largest SigmaR

4. MU for the spot that yielded the largest SigmaR

5. BG for the spot that yielded the largest SigmaEta

6. SigmaGEta for the spot that yielded the largest SigmaEta

7. SigmaLEta for the spot that yielded the largest SigmaEta

8. MU for the spot that yielded the largest SigmaR

I figure you can probably do this more efficiently. But, this is a good example of the type of information we'd like to gather for a more detailed spot spreading analysis.

MergeOverlappingPeaksAllZarr.zip