Note that there are some explanatory texts on larger screens.

plurals
  1. POFind nucleus involved in nCapture process in Geant4
    primarykey
    data
    text
    <p>I'm trying to determine the target nucleus in an (n,gamma) reaction in Geant4. I have been unable to extract this information. The only place I have found it to be stored is in G4IsoParticleChange which is created by the G4HadronicProcess if Isotope Counting is enabled.</p> <p>Unfortunately this information is lost (not in a nice way, either, but memory leak style) every time the process is invoked. I cannot find a user hook-in to ask for this information in between particle creations. G4UserStackingAction is not sufficient as far as I can tell because the secondary particles are all created for a step before they are stacked (* though writing that last sentence has given me an idea).</p> <p>Could anyone help me determine the nucleus which captures the neutron in the (n,gamma) reaction? Is there an easier way to get this information?</p> <p>Thanks</p> <p>P.S. Since a neutron can only be absorbed once within a step, would it be safe to just wait until the secondaries are stacked to obtain the IsoParticleChange info or am I risking a memory leak?</p> <p>Edit to be more clear:</p> <p>I'm asking if someone knows how to retrieve the nucleus which was the target in the nCapture process in Geant4. It is clear from the source that a memory leak will occur whenever the G4HadronicProcess is invoked if the G4IsoParticleChange information is not retrieved. There appears to be no user hook in point which is appropriate to grab this without missing some information and causing said memory leak, and yet it is possible to turn this information storage on. I am wondering what the correct way to grab this information is, or if there is a better way to obtain this type of info about the target nucleus.</p> <p>I have previously gotten some help on SO from people who are familiar with Geant. I have not been able to get a response from the Slac Geant4 forum to give me access to post there. The forum doesn't seem terribly active, anyway.</p> <p>Thanks</p> <p>Follow-up:</p> <p>In case someone stumbles upon this, the answer I got from one of the authors of the G4HadronicProcess class was "This part of the implementation has not been maintained for many years" and "Getting the nucleus involved in the interaction is very difficult if not impossible without editing the source code." So I'm in the process of setting up a new Geant workspace where I can do that. FYI.</p>
    singulars
    1. This table or related slice is empty.
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
 

Querying!

 
Guidance

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload