|Message: Re: Has Geant4 code to simulate neutino simulation?||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
There is mention of neutrino interactions in section 25.8 of the Physics Manual (Geant4 9.4).
There is a class G4QNeutrinoPhysics that attempts to implement neutrino interactions. My playing with it indicates it does not behave reasonably at all, and frequently crashes ("cannot recover GSM", whatever that means). The code is opaque and not easily fixed or modified.
I have been working on implementing neutrino interactions in G4beamline (http://g4beamline.muonsinc.com), which is based on the Geant4 toolkit. My approach is to interface to the Genie neutrino monte-carlo program (http://www.genie-mc.org/), via a new physics process which uses pipes to a new interface program that calls Genie code -- this way the two huge codes do not need to be linked as a single executable. The user specifies a list of materials in which neutrinos have a specified interaction length, and particle weights are adjusted accordingly; you can set it up so they traverse 1,000 km of rock and then interact in the last few meters near the detector and in the detector.
At present it is working and close to release (a week or two of testing remain), but note it uses G4beamline, not just Geant4, and depends heavily on G4beamline infrastructure. I suppose the basic process code could be adapted to other programs, or you could use G4beamline :-). It is intended for studying neutrino backgrounds and radiation from a muon collider, and does not do oscillations. The physics-process code is only about 150 lines long (not counting the millions of lines in Genie (:-)), and is not opaque, so oscillations could probably be added (but that is a non-trivial physics problem).