change S1 split threshold main contexts
Created by: jorana
For the data we are currently taking we see that too many peaks are split. This is due to the fact that we are detecting a lot of large S1s, the split threshold for these S1-like signals is to low. This PR solves this.
For more info see Raise peak splitting requirement for S1-like signals in GXe on the XENONnT wiki as this repo is public.
Implications of merging this PR If we merge this PR, that means that effectively all nT data we currently have is out of date and could not be able to be loaded using this context until we rebuild from records. Only the the datatypes in the figure below would be loadable of the data we currently have. This includes simulation data!
Everything above this we would have to reprocess (see https://straxen.readthedocs.io/en/latest/reference/datastructure.html#event-positions).
1T data / context / analyses would not be affected.
Alternative(s)
We could alternatively change only the xenonnt_online context by adding lines and inserting context
here:
config = straxen.contexts.xnt_common_config
config['peak_split_gof_threshold'] = (
None, # Reserved
((0.5, 1), (4, 0.4)),
((2, 1), (4.5, 0.4)))
I choose not to do this as the simulation- and data-contexts would diverge and we couldn't compare the two easily anymore.
Related PRs and notes
- https://github.com/AxFoundation/strax/pull/225
- https://github.com/XENONnT/straxen/pull/45
- Peak splitting note
- Natural breaks splitting
- Example wfs
Unrelated I fixed some error in the simulation context.