The neutrino - nucleon cross sections have been coded up and will be added to the VLE directory of GENIE soon (just need CVS permissions).
There are still some open issues. While waiting for responses on them, I thought it would be useful to show the current status. When the issues get resolved, I will plan to document the making of this class, in case it may be helpful to others.
Current cross section calculations:
or on log scale:
View cross sections on log scale.
The IBD x-sec agrees very well with the table displayed in the paper:
View ratio of IBD cross section: GENIE to values in the paper.
However, this agreement is only achievable with:
- The proton anomalous magnetic moment being 1.793, as given in the paper. PDG (and GENIE) have k_p = 2.793. I have emailed Strumia regarding this issue, but so far haven't received a response. Update: After talking with the authors, the parameter used is the anomalous magnetic moment, which is k_p-1 for the proton (and k_n for the neutron).
- The minimum q^2, used for the integral to calculate the total cross section, is limited in GENIE to not be smaller than 10^-4. For the plots above, this limit had been changed to 10^-10.
The limit of q^2 >= 10^-4 (default in GENIE) can cause some trouble for very low energy neutrinos:
View ratio (GENIE to paper) with default q^2 >= 10^-4
I've contacted Constantinos to talk about possible solutions.
Comments (2)
I don't know about any one else, but I would find your documentation on the creation of this class extremely helpful. I have played around with the code, but still only have a novices conception of how to go about adding VLE code to GENIE.
Posted by Cully Little | September 4, 2009 7:56 PM
Posted on September 4, 2009 19:56
Yes, sorry for the delay. I'm sending Constantinos my proposed changes outside of our VLE package (so that we can use a lower Q2) today. Once that's in place, I'll commit the class I wrote for the VLE package and post a macro which runs it.
The documentation will have to wait a little while longer, though, because I've got two meetings at the end of this month to prepare for...
Posted by Corey Reed | September 7, 2009 11:02 AM
Posted on September 7, 2009 11:02