Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Failure of generators generated with gevgen (18O with nu_e_bar, G21_11b, latest version of genie) #383

Open
WentaiLuo opened this issue Feb 7, 2024 · 0 comments

Comments

@WentaiLuo
Copy link

I used the latest version of genie that is producing nu_e_bar with oxygen 18 generators, the energy spectrum used is Kam's Atmospheric Neutrino Spectrum, which produces 1,500,000 events with energies of 0.01-1 GeV. when running up to 468897 there is an error reported, with the error being

gevgen: Spline.cxx:366: double genie::Spline::Evaluate(double) const: Assertion `!TMath::IsNaN(x)' failed.

Aborted

make: *** [Makefile:18: run_-12_1000080180_G21_11_b] Error

The rest of the information for this event is as follows
[-] Generated flux neutrino:
|----o PDG-code : -12
|----o 4-momentum : (E = 0.156052, px = 0, py = 0, pz = 0.156052) / M^2 = 0 / P = 0.156052
|----o 4-position : (t = 0, x = 0, y = 0, z = 0)
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateEvent1Try (886)> : Computing interaction probabilities for max. path lengths
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::ComputeInteractionProbabilities (1112)> : Computing relative interaction probabilities for each material
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateEvent1Try (891)> : The no-interaction probability (max. path lengths) is: 95.2386 %
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateEvent1Try (902)> : ** Rejecting current flux neutrino
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateEvent (833)> : Flux neutrino didn't interact - Trying the next one...
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateFluxNeutrino (1042)> : Generating a flux neutrino
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateFluxNeutrino (1056)> :
[-] Generated flux neutrino:
|----o PDG-code : -12
|----o 4-momentum : (E = 0.0270707, px = 0, py = 0, pz = 0.0270707) / M^2 = 0 / P = 0.0270707
|----o 4-position : (t = 0, x = 0, y = 0, z = 0)
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateEvent1Try (886)> : Computing interaction probabilities for max. path lengths
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::ComputeInteractionProbabilities (1112)> : Computing relative interaction probabilities for each material
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateEvent1Try (891)> : The no-interaction probability (max. path lengths) is: 99.9449 %
1707242135 NOTICE GMCJDriver : [n] <GMCJDriver.cxx::GenerateEvent1Try (902)> : ** Rejecting current flux neutrino
1707242135 NOTICE GM

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant