Skip to content
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

physics: add option to disable long-lived daughter nuclei in decay chains #70

Merged
merged 1 commit into from
May 8, 2024

Conversation

ManuelHu
Copy link
Collaborator

@ManuelHu ManuelHu commented Apr 11, 2024

this code is somewhat inspired by MaGe, but reworked and modernized. The functionality can be enabled and configured via a macro command

Alternatives (we should document this somewhere...)

/process/had/rdm/thresholdForVeryLongDecayTime <threshold> ms

this has two problems:

  1. the time limit is applied to the sampled global time of the decay of an unstable nucleus. This is not deterministic (i.e. it is dependent on the random seed), and also not really helpful to limit decays to specific isotopes
  2. This also applies to primaries created with GPS, so one cannot simulate any decay to create secondaries.
/process/had/rdm/nucleusLimits <AMin> <AMax> <ZMin> <ZMax>

This method works and is reliable/deterministic, but you have to define the limits carefully, especially if you have multiple possible primary decays.

fixes #39

@ManuelHu ManuelHu merged commit ddd8fc2 into legend-exp:main May 8, 2024
5 checks passed
@ManuelHu ManuelHu deleted the lifetime branch May 8, 2024 09:29
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Kill daughter nuclei after radioactive decay if half-life is longer than 1 microsec
1 participant