 |
 |
Dear GENIE developers, It seems there are a few bugs in "./src/EVGModules/NucDeExcitationSim.cxx", which is for the deexcitation of the residual nucleus from ν+O. The code is not consistent with the papers it is based on. The NucDeExcitationSim.cxx is based on Ejiri1993 ( http://inspirehep.net/record/365249) and a measurement ( https://www.sciencedirect.com/science/article/pii/S0920563204007364). For the p hole. The Ejiri1993 gives that the s1/2 p hole branch ratio (BR) is 28%, while the p1/2 25% and p3/2 47%. And the measurement gives more details about the γ BR from the s1/2 sate. The total γ (>3 MeV) fraction from 1/2s p hole is 8%. (That's to say there is a probability of 8/28=0.2857 that s1/2 p hole gives γ > 3 MeV). However, in line 114 and line 126-131, the NucDeExcitationSim.cxx seems wrongly set all the s1/2 p hole give γ > 3 MeV. I means, the probabilites for different final γ's are scaled by 28/8=3.5. So an overprediction of γ fraction is made. I think it can be easily fixed by changing some numbers in these lines. Also, in line 141-142, it seems there is typos of the probabilities. The 0.05, 0.033 0.017 should be 0.5, 0.33 and 0.17 which matches the measurement. Obviously the sun of former does not give 1. Thanks!! Best, Bei -- Bei Zhou Graduate student, Center for Cosmology and AstroParticle Physics (CCAPP), The Ohio State University ######################################################################## To unsubscribe from the NEUTRINO-MC-SUPPORT list, click the following link: https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=NEUTRINO-MC-SUPPORT&A=1