 |
 |
Hi Bei, I don't have time to go through the paper you mention right now, but for a small fix like this, feel free to make a pull request against our GitHub mirror https://github.com/GENIEMC/GENIE <https://github.com/GENIEMC/GENIE> Normally I would recommend targeting the 'community' branch, but we're close to some big changes in GENIE, so right now I'd prefer you make changes on a copy of the 'trunk' branch and issue a pull request targeting that branch instead. Thanks for checking the numbers carefully! Gabriel Perdue Associate Scientist Scientific Computing Division Fermi National Accelerator Laboratory PO Box 500, MS 234, Batavia, IL 60510, USA [log in to unmask] Office: 630-840-6499 Cell: 630-605-8062 Connect with us! <http://www.fnal.gov/>Web <http://www.fnal.gov/> | Facebook <http://www.facebook.com/Fermilab> | Twitter <https://www.twitter.com/fermilab> > On Jul 6, 2018, at 12:13 PM, Bei Zhou <[log in to unmask]> wrote: > > 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 <https://urldefense.proofpoint.com/v2/url?u=http-3A__inspirehep.net_record_365249&d=DwMFaQ&c=gRgGjJ3BkIsb5y6s49QqsA&r=VuX71VzHBLywwZNXKJh6xw&m=qG6QkPxLFQKhpjvKLjdcl_2Jj1wzkWPoxfLcdD_PpAQ&s=4DlzG11pETOYhvM7XsQ13QZhjwWxJtTp49uKdlXL3p4&e=>) and a measurement (https://www.sciencedirect.com/science/article/pii/S0920563204007364 <https://urldefense.proofpoint.com/v2/url?u=https-3A__www.sciencedirect.com_science_article_pii_S0920563204007364&d=DwMFaQ&c=gRgGjJ3BkIsb5y6s49QqsA&r=VuX71VzHBLywwZNXKJh6xw&m=qG6QkPxLFQKhpjvKLjdcl_2Jj1wzkWPoxfLcdD_PpAQ&s=gpZ-CMdGcZSIreGcIu6PdmZmXtIYUzh-K8_L6ch-bKE&e=>). > > 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 <https://urldefense.proofpoint.com/v2/url?u=https-3A__www.jiscmail.ac.uk_cgi-2Dbin_webadmin-3FSUBED1-3DNEUTRINO-2DMC-2DSUPPORT-26A-3D1&d=DwMFaQ&c=gRgGjJ3BkIsb5y6s49QqsA&r=VuX71VzHBLywwZNXKJh6xw&m=qG6QkPxLFQKhpjvKLjdcl_2Jj1wzkWPoxfLcdD_PpAQ&s=WxWaNtT3F0YM7Gwrbre2dA_03rN8jwEY-rwOn5Vd4HY&e=> ######################################################################## 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