Tom,
Multiple issues here.
On 5/27/2011 4:44 PM, Tom Roberts wrote:
> I am working on interfacing Genie into G4beamline, which is a
> single-particle-tracking simulation program based on Geant4. I want to
> be able to directly simulate neutrino backgrounds from a muon
> collider. My approach is to write a main program for Genie (based on
> gevgen) that reads neutrino tracks from a pipe, and writes the
> secondary tracks to a different pipe; the other ends of the pipes are
> handled by G4beamline. I will arrange to track neutrinos normally
> through most rock and dirt (with very large steps), but for the last
> few meters it will use an artificial cross-section that makes a few
> percent of them interact in the last few meters; it will then track
> the secondaries as usual.
>
I hope this is consistent with the copyright, Costas will know for sure.
> To compute the weights of the secondary tracks I need to know the
> total inelastic cross-section of this particular neutrino (type,
> energy) on this particular material, computing it for each neutrino
> interaction.
>
> Given the neutrino PDGid and energy, how can I compute its total
> inelastic cross-section?
>
You have the total CC and NC cross sections available in spline. One
definition of the total inel xs
is all NC except NC elastic. Another would add all the CC xs's.
> If there is a routine that can be used for a mixed target, that would
> be great, but I could also put together the individual cross-sections
> for individual atoms in the mixed material.
>
check geometry section of User Manual. There are various ways to do
this. You can use a list or
a root file.
regards,
Steve
>
> Tom Roberts
|