[Prev][Next][Index][Thread]
Weak interaction generator
Hi. I've finally succeeded in testing the v3 implementation of
the weak decay code that I sent some time ago. There were a couple
of small errors to correct in hypdec for the code to run at all.
In addition, I found a very rare divide by zero exception in PROCOM1,
which is deep into Angel Ramos's code and put in some protection
from that. I was able to compare the results with earlier ones that
were given in her paper. I left in the code a couple of histos
created just for the code testing; these could be commented out, but
it will be convenient to leave them around if further testing is required.
I haven't looked at how well the fitting procedure for the
proton in everec works.
The changes will be sent to Gervasio as a correction cradle, together
with the one needed to load in the GEANT commons when they are included in
the fidamc.cmz.
In response to Aldo's question about the FLUKA call in GUHADR, that
call is made to test for normal hadronic interactions when tracking
hadrons. When a K- stops in the target, a hypernucleus is formed,
GUHADR executes HYPDEC and then returns. Otherwise FLUKA is called to
invoke possible strong interactions.
For the case of stopping K- the FLUKA code is not very good. I've
written a routine KMICAPT which identifies the nucleus on which
the K is captured, gets a nucleon Fermi momentum, interacts to
produce a quasi-free lambda or sigma. However I don't feel that this
code has been sufficiently tested to release, so it's commented out
at present.
Ciao, Art