-
Notifications
You must be signed in to change notification settings - Fork 1
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Replace scintillation time clipping by truncation #220
Conversation
Pull Request Test Coverage Report for Build 9193566598Details
💛 - Coveralls |
…nT/fuse into no_s1_photon_time_clipping
A note discussing the proposed changes can be found here. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This looks good! Like previously discussed, I prefer this re-iterative loop than assigning times for many more photons than needed, based on a safety factor. Is there a speed improvement this way?
Hi @ramirezdiego, I think there is a slight performance increase using the code in this PR compared to the main branch. I did not run detailed performance tests but for the simulations done in the note above the raw records simulation time decreased by about 10%. Maybe this is not relevant to the overall simulation time when S2s are also included. |
Thanks, I was actually just curious about the improved speed with respect to the previous approach in this PR. |
What does the code in this PR do / what does it improve?
This PR will replace the s1 photon timing clipping by truncation.
Can you briefly describe how it works?
Right now we are clipping the sampled S1 photon scintillation times in the S1PhotonPropagation plugin to the maximum value of
maximum_recombination_time
. To minimize the effect on the S1 waveform, this PR replaces the clipping by a truncation of the timing distribution at he samemaximum_recombination_time
. Two safety measures are added to make sure that we always have enough photon times to sample from. As a bonus theS1PhotonPropagation.photon_timings
function was cleaned up and the scintillation timing part moved to a dedicated function.Please include the following if applicable: