decaychain generator segfaults when using chains with alpha decays

[metelko]

an example. If beta_decays.dat hacked to give the chain
CHAIN: 238U
4
238U 0001 1.0000 2 1.410e+17
234Th 0001 1.0000 1 2.080e+06
234Pa 0001 1.0000 1 2.412e+04
#234U 0001 1.0000 2 7.747e+12
230Th 0001 1.0000 2 2.379e+12
This works...
CHAIN: 238U
4
238U 0001 1.0000 2 1.410e+17
234Th 0001 1.0000 1 2.080e+06
234Pa 0001 1.0000 1 2.412e+04
234U 0001 1.0000 2 7.747e+12
#230Th 0001 1.0000 2 2.379e+12
seg faults.

NB Just using 234U as the source (not a chain) works too..

[MorganAskins]

Along these lines, we should consider merging in the Decay0 generator (which came from Snoman into SNO+ Rat and eventually Andy merged into the default rat-pac). It's here:
https://github.com/rat-pac/rat-pac/blob/master/src/gen/Decay0.cc
https://github.com/rat-pac/rat-pac/blob/master/src/gen/Decay0.hh

[metelko]

Hi Morgan, a question about how the decay chain works. If simulating an event, does the simulator pick a decay from a random decay from the decay chain assuming all the decay products are in equilibrium? So if there's currently a problem with the decay chain could this be simulated by eg 100KBq of 238U 100kBq of 234Th 100kBq of 234Pa etc for the rest of the chain. Thanks, Carl

…

On 2020-02-27 18:23, Morgan Askins wrote: Along these lines, we should consider merging in the Decay0 generator (which came from Snoman into SNO+ Rat and eventually Andy merged into the default rat-pac). It's here: https://github.com/rat-pac/rat-pac/blob/master/src/gen/Decay0.cc https://github.com/rat-pac/rat-pac/blob/master/src/gen/Decay0.hh -- You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub [1], or unsubscribe [2]. Links: ------ [1] #80?email_source=notifications&email_token=AMCGDXDYK6G7KAOIFEVA4XTRFAAJNA5CNFSM4K4XM5BKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOENFNB3Q#issuecomment-592105710 [2] https://github.com/notifications/unsubscribe-auth/AMCGDXBHAMWFL5NKEMX3NUTRFAAJNANCNFSM4K4XM5BA