search

GENIE-MC / Generator

The popular GENIE Generator product is used by nearly all accelerator neutrino experiments and it plays a key role in the exploitation of neutrino data. The Generator implements a modern software framework and it includes state-of-the-art physics modules. It captures the latest results of the GENIE global analysis of neutrino scattering data and includes several tunes that were produced using the proprietary Comparisons and Tuning products. The GENIE physics model is universal and comprehensive: It handles all neutrinos and targets, and all processes relevant from MeV to PeV energy scales. The Generator includes several tools (flux drivers, detector geometry navigators, specialized event generation apps, event reweighting engines) to simulate complex experimental setups in full detail and to support generator-related analysis tasks.
http://www.genie-mc.org
42 stars 92 forks source link

Calculation of EM RES XSec #347

Closed jtenavidal closed 10 months ago

jtenavidal commented 10 months ago

The calculation of EM RES XSec on free nucleon takes very long time. This issue will look into possible ways to optimise it

jtenavidal commented 10 months ago

For neutrinos, we use the ReinSehgalRESXSecFast.cxx class which inherits from ReinSehgalRESXSecWithCacheFast.cxx for integration. A counter has been added to compare with electrons.

gmkspl -p 14 -t 1000000010 -n 30 -e 30 --tune G18_02a_00_000 -o e_on_n_EMRES.xml --event-generator-list CCRES

The output shows that at max it takes 1s to calculate a spline value:

 ** Creating cache branch - key = genie::BergerSehgalRESPXSec2014/NoPauliBlock/ResExcitationXSec/R:P33(1232);nu:14;int:Weak[CC]n
1698243191 NOTICE ReinSehgalResCF : [n] <ReinSehgalRESXSecWithCacheFast.cxx::CacheResExcitationXSec (133)> : E threshold = 0.110881
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.01) = 0 x 1E-38 cm^2, evaluated in 7.59e-07 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.0301762) = 0 x 1E-38 cm^2, evaluated in 2.61e-07 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.0503525) = 0 x 1E-38 cm^2, evaluated in 2.3e-07 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.0705287) = 0 x 1E-38 cm^2, evaluated in 2.17e-07 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.090705) = 0 x 1E-38 cm^2, evaluated in 1.86e-07 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.110881) = 0 x 1E-38 cm^2, evaluated in 1.86e-07 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.120416) = 0 x 1E-38 cm^2, evaluated in 0.00218546 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.13077) = 0 x 1E-38 cm^2, evaluated in 0.000303582 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.142014) = 0 x 1E-38 cm^2, evaluated in 0.000287658 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.154226) = 0 x 1E-38 cm^2, evaluated in 0.000287645 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.167487) = 0 x 1E-38 cm^2, evaluated in 0.000330438 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.181889) = 0 x 1E-38 cm^2, evaluated in 0.000291976 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.197529) = 0 x 1E-38 cm^2, evaluated in 0.000292627 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.214514) = 0 x 1E-38 cm^2, evaluated in 0.000297099 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.232959) = 0 x 1E-38 cm^2, evaluated in 0.000324255 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.25299) = 0 x 1E-38 cm^2, evaluated in 0.000295287 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.274744) = 0 x 1E-38 cm^2, evaluated in 0.000298219 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.298369) = 2.98498e-06 x 1E-38 cm^2, evaluated in 0.0058318 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.324025) = 2.89092e-05 x 1E-38 cm^2, evaluated in 0.0536578 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.351887) = 0.000142363 x 1E-38 cm^2, evaluated in 0.0697802 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.382144) = 0.000537289 x 1E-38 cm^2, evaluated in 0.0751991 s 
1698243191 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.415004) = 0.00180187 x 1E-38 cm^2, evaluated in 0.0823352 s 
1698243192 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.450689) = 0.00548223 x 1E-38 cm^2, evaluated in 0.0874819 s 
1698243192 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.489442) = 0.0133552 x 1E-38 cm^2, evaluated in 0.119797 s 
1698243192 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.531528) = 0.0251629 x 1E-38 cm^2, evaluated in 0.151814 s 
1698243192 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.577232) = 0.040043 x 1E-38 cm^2, evaluated in 0.191684 s 
1698243192 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.626867) = 0.0573684 x 1E-38 cm^2, evaluated in 0.241587 s 
1698243192 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.680769) = 0.0765227 x 1E-38 cm^2, evaluated in 0.254264 s 
1698243193 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.739306) = 0.096846 x 1E-38 cm^2, evaluated in 0.295548 s 
1698243193 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.802877) = 0.117671 x 1E-38 cm^2, evaluated in 0.321787 s 
1698243193 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.871914) = 0.138372 x 1E-38 cm^2, evaluated in 0.353852 s 
1698243194 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.946887) = 0.158404 x 1E-38 cm^2, evaluated in 0.376094 s 
1698243194 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.02831) = 0.177332 x 1E-38 cm^2, evaluated in 0.412884 s 
1698243195 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.11673) = 0.194837 x 1E-38 cm^2, evaluated in 0.451079 s 
1698243195 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.21275) = 0.210719 x 1E-38 cm^2, evaluated in 0.4882 s 
1698243196 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.31703) = 0.224764 x 1E-38 cm^2, evaluated in 0.478178 s 
1698243196 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.43028) = 0.236621 x 1E-38 cm^2, evaluated in 0.502887 s 
1698243197 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.55327) = 0.246491 x 1E-38 cm^2, evaluated in 0.538583 s 
1698243197 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.68683) = 0.254582 x 1E-38 cm^2, evaluated in 0.542524 s 
1698243198 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.83187) = 0.261124 x 1E-38 cm^2, evaluated in 0.572472 s 
1698243198 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=1.98939) = 0.266345 x 1E-38 cm^2, evaluated in 0.575308 s 
1698243199 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=2.16045) = 0.270459 x 1E-38 cm^2, evaluated in 0.600859 s 
1698243200 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=2.34622) = 0.273659 x 1E-38 cm^2, evaluated in 0.630278 s 
1698243200 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=2.54797) = 0.276113 x 1E-38 cm^2, evaluated in 0.664909 s 
1698243201 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=2.76706) = 0.277964 x 1E-38 cm^2, evaluated in 0.67582 s 
1698243202 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=3.00499) = 0.279333 x 1E-38 cm^2, evaluated in 0.707687 s 
1698243202 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=3.26338) = 0.280319 x 1E-38 cm^2, evaluated in 0.71414 s 
1698243203 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=3.54399) = 0.281005 x 1E-38 cm^2, evaluated in 0.745496 s 
1698243204 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=3.84873) = 0.281455 x 1E-38 cm^2, evaluated in 0.765968 s 
1698243205 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=4.17967) = 0.281724 x 1E-38 cm^2, evaluated in 0.802026 s 
1698243205 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=4.53906) = 0.281854 x 1E-38 cm^2, evaluated in 0.806288 s 
1698243206 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=4.92936) = 0.281879 x 1E-38 cm^2, evaluated in 0.820928 s 
1698243207 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=5.35322) = 0.281827 x 1E-38 cm^2, evaluated in 0.832262 s 
1698243208 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=5.81353) = 0.281717 x 1E-38 cm^2, evaluated in 0.844114 s 
1698243209 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=6.31342) = 0.281568 x 1E-38 cm^2, evaluated in 0.864619 s 
1698243210 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=6.85629) = 0.281392 x 1E-38 cm^2, evaluated in 0.877327 s 
1698243211 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=7.44584) = 0.281198 x 1E-38 cm^2, evaluated in 0.924361 s 
1698243212 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=8.08609) = 0.280995 x 1E-38 cm^2, evaluated in 0.925233 s 
1698243213 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=8.78139) = 0.280788 x 1E-38 cm^2, evaluated in 0.92912 s 
1698243213 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=9.53647) = 0.280582 x 1E-38 cm^2, evaluated in 0.94142 s 
1698243214 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=10.3565) = 0.28038 x 1E-38 cm^2, evaluated in 0.945864 s 
1698243215 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=11.247) = 0.280185 x 1E-38 cm^2, evaluated in 0.969264 s 
1698243216 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=12.2141) = 0.279997 x 1E-38 cm^2, evaluated in 0.96165 s 
1698243217 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=13.2644) = 0.279818 x 1E-38 cm^2, evaluated in 0.987086 s 
1698243218 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=14.4049) = 0.279649 x 1E-38 cm^2, evaluated in 1.00671 s 
1698243219 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=15.6436) = 0.279491 x 1E-38 cm^2, evaluated in 1.02463 s 
1698243220 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=16.9887) = 0.279342 x 1E-38 cm^2, evaluated in 1.0315 s 
1698243221 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=18.4495) = 0.279203 x 1E-38 cm^2, evaluated in 1.02667 s 
1698243222 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=20.0359) = 0.279074 x 1E-38 cm^2, evaluated in 1.03718 s 
1698243223 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=21.7588) = 0.278955 x 1E-38 cm^2, evaluated in 1.0506 s 
1698243225 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=23.6297) = 0.278844 x 1E-38 cm^2, evaluated in 1.04072 s 
1698243226 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=25.6616) = 0.278742 x 1E-38 cm^2, evaluated in 1.07247 s 
1698243227 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=27.8681) = 0.278649 x 1E-38 cm^2, evaluated in 1.09897 s 
1698243228 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=30.2644) = 0.278563 x 1E-38 cm^2, evaluated in 1.10574 s 
1698243229 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=32.8668) = 0.278484 x 1E-38 cm^2, evaluated in 1.09874 s 
1698243230 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=35.6929) = 0.278411 x 1E-38 cm^2, evaluated in 1.10668 s 
1698243231 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=38.762) = 0.278344 x 1E-38 cm^2, evaluated in 1.10271 s 
1698243232 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=42.0951) = 0.278283 x 1E-38 cm^2, evaluated in 1.15429 s 
1698243233 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=45.7147) = 0.278227 x 1E-38 cm^2, evaluated in 1.18503 s 
1698243235 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=49.6456) = 0.278176 x 1E-38 cm^2, evaluated in 1.11436 s 
1698243236 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=53.9144) = 0.27813 x 1E-38 cm^2, evaluated in 1.17803 s 
1698243237 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=58.5504) = 0.278087 x 1E-38 cm^2, evaluated in 1.17971 s 
1698243238 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=63.5849) = 0.278049 x 1E-38 cm^2, evaluated in 1.18591 s 
1698243239 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=69.0524) = 0.278014 x 1E-38 cm^2, evaluated in 1.14291 s 
1698243240 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=74.99) = 0.277982 x 1E-38 cm^2, evaluated in 1.18382 s 
1698243242 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=81.4382) = 0.277954 x 1E-38 cm^2, evaluated in 1.21148 s 
1698243243 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=88.4408) = 0.277928 x 1E-38 cm^2, evaluated in 1.20471 s 
1698243244 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=96.0456) = 0.277904 x 1E-38 cm^2, evaluated in 1.18037 s 
1698243245 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=104.304) = 0.277883 x 1E-38 cm^2, evaluated in 1.24815 s 
1698243246 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=113.273) = 0.277864 x 1E-38 cm^2, evaluated in 1.2118 s 
1698243248 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=123.013) = 0.277847 x 1E-38 cm^2, evaluated in 1.17622 s 
1698243249 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=133.591) = 0.277832 x 1E-38 cm^2, evaluated in 1.18973 s 
1698243250 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=145.078) = 0.277818 x 1E-38 cm^2, evaluated in 1.20754 s 
1698243251 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=157.552) = 0.277806 x 1E-38 cm^2, evaluated in 1.18528 s 
1698243252 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=171.1) = 0.277795 x 1E-38 cm^2, evaluated in 1.17943 s 
1698243254 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=185.812) = 0.277785 x 1E-38 cm^2, evaluated in 1.19315 s 
1698243255 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=201.79) = 0.277777 x 1E-38 cm^2, evaluated in 1.21562 s 
1698243256 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=219.141) = 0.277769 x 1E-38 cm^2, evaluated in 1.22297 s 
1698243257 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=237.984) = 0.277762 x 1E-38 cm^2, evaluated in 1.2239 s 
1698243259 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=258.448) = 0.277756 x 1E-38 cm^2, evaluated in 1.23973 s 
1698243260 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=280.671) = 0.277751 x 1E-38 cm^2, evaluated in 1.22138 s 
1698243261 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=304.805) = 0.277746 x 1E-38 cm^2, evaluated in 1.24185 s 
1698243262 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=331.014) = 0.277742 x 1E-38 cm^2, evaluated in 1.21714 s 
1698243263 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=359.477) = 0.277738 x 1E-38 cm^2, evaluated in 1.21694 s 
1698243265 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=390.388) = 0.277735 x 1E-38 cm^2, evaluated in 1.23957 s 
1698243266 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=423.956) = 0.277732 x 1E-38 cm^2, evaluated in 1.24432 s 
1698243267 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=460.411) = 0.27773 x 1E-38 cm^2, evaluated in 1.26089 s 
1698243268 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=500) = 0.277728 x 1E-38 cm^2, evaluated in 1.27813 s
jtenavidal commented 10 months ago

For electrons: gmkspl -p 11 -t 1000000010 -n 30 -e 30 --tune G18_02a_00_000 -o e_on_n_EMRES.xml --event-generator-list EMRES

When using ReinSehgalRESXSecWithCacheFast: It holds after: 1698243378 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.187738) = 0 x 1E-38 cm^2, evaluated in 0.000163174 s

When using ReinSehgalRESXSecWithCache (default in GENIE for electrons):

1698243504 NOTICE XSecSplLst : [s] <CreateSpline (200)> : Energy threshold for current interaction = 0 GeV
1698243504 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (128)> : No cached RES v-production data for input neutrino (pdgc: 11)
1698243504 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (131)> : Wait while computing/caching RES production xsec first...
1698243504 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (110)> : 
 ** Creating cache branch - key = genie::BergerSehgalRESPXSec2014/EM-NoPauliBlock/ResExcitationXSec/R:P33(1232);nu:11;int:EMn
1698243504 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (118)> : E threshold = 0
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.01) = 0 x 1E-38 cm^2, evaluated in 5.2845e-05s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0110651) = 0 x 1E-38 cm^2, evaluated in 3.482e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0122437) = 0 x 1E-38 cm^2, evaluated in 2.867e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0135478) = 0 x 1E-38 cm^2, evaluated in 2.793e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0149908) = 0 x 1E-38 cm^2, evaluated in 2.809e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0165875) = 0 x 1E-38 cm^2, evaluated in 2.672e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0183542) = 0 x 1E-38 cm^2, evaluated in 2.748e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0203092) = 0 x 1E-38 cm^2, evaluated in 2.675e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0224723) = 0 x 1E-38 cm^2, evaluated in 2.697e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0248659) = 0 x 1E-38 cm^2, evaluated in 2.569e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0275144) = 0 x 1E-38 cm^2, evaluated in 2.721e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.030445) = 0 x 1E-38 cm^2, evaluated in 2.638e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0336878) = 0 x 1E-38 cm^2, evaluated in 2.713e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0372759) = 0 x 1E-38 cm^2, evaluated in 2.649e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0412463) = 0 x 1E-38 cm^2, evaluated in 2.663e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0456395) = 0 x 1E-38 cm^2, evaluated in 2.636e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0505006) = 0 x 1E-38 cm^2, evaluated in 2.779e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0558795) = 0 x 1E-38 cm^2, evaluated in 2.684e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0618314) = 0 x 1E-38 cm^2, evaluated in 2.688e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0684171) = 0 x 1E-38 cm^2, evaluated in 2.705e-06s 
1698243504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0757044) = 0 x 1E-38 cm^2, evaluated in 2.76e-06s 
1698243509 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0837678) = 0 x 1E-38 cm^2, evaluated in 5.31494s 
1698243515 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.09269) = 0 x 1E-38 cm^2, evaluated in 5.68847s 
1698243521 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.102563) = 0 x 1E-38 cm^2, evaluated in 5.99191s 
1698243527 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.113487) = 0 x 1E-38 cm^2, evaluated in 6.26513s 
1698243533 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.125574) = 0 x 1E-38 cm^2, evaluated in 6.48051s 
1698243540 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.13895) = 0 x 1E-38 cm^2, evaluated in 6.66477s 
1698243547 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.153749) = 0 x 1E-38 cm^2, evaluated in 6.80111s 
1698243554 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.170125) = 0 x 1E-38 cm^2, evaluated in 6.93199s 
1698243561 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.188246) = 15.579 x 1E-38 cm^2, evaluated in 7.47479s 
1698243569 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.208296) = 2552.78 x 1E-38 cm^2, evaluated in 7.45167s 
1698243576 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.230482) = 23817.4 x 1E-38 cm^2, evaluated in 7.40734s 
1698243584 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.255031) = 115189 x 1E-38 cm^2, evaluated in 7.49362s 
1698243591 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.282195) = 430522 x 1E-38 cm^2, evaluated in 7.48574s 
1698243599 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.312252) = 1.4437e+06 x 1E-38 cm^2, evaluated in 7.49681s 
1698243606 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.345511) = 4.26039e+06 x 1E-38 cm^2, evaluated in 7.51464s 
1698243614 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.382312) = 9.28245e+06 x 1E-38 cm^2, evaluated in 7.51461s 
1698243621 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.423032) = 1.50228e+07 x 1E-38 cm^2, evaluated in 7.66426s 
1698243629 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.46809) = 2.03853e+07 x 1E-38 cm^2, evaluated in 7.80291s 
1698243637 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.517947) = 2.50768e+07 x 1E-38 cm^2, evaluated in 7.98413s 
1698243645 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.573115) = 2.91084e+07 x 1E-38 cm^2, evaluated in 8.20845s 
1698243654 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.634158) = 3.25787e+07 x 1E-38 cm^2, evaluated in 8.44144s 
1698243662 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.701704) = 3.5596e+07 x 1E-38 cm^2, evaluated in 8.59592s 
1698243671 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.776444) = 3.82599e+07 x 1E-38 cm^2, evaluated in 8.77757s 
1698243680 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.859144) = 4.0652e+07 x 1E-38 cm^2, evaluated in 8.94694s 
1698243689 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.950653) = 4.28333e+07 x 1E-38 cm^2, evaluated in 9.00158s 
1698243698 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=1.05191) = 4.4849e+07 x 1E-38 cm^2, evaluated in 9.13706s 
1698243707 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=1.16395) = 4.65587e+07 x 1E-38 cm^2, evaluated in 9.20944s 
1698243717 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=1.28792) = 4.7948e+07 x 1E-38 cm^2, evaluated in 9.30677s 
1698243726 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=1.4251) = 4.92394e+07 x 1E-38 cm^2, evaluated in 9.27564s 
1698243735 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=1.57689) = 5.04703e+07 x 1E-38 cm^2, evaluated in 9.28437s 
1698243744 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=1.74485) = 5.16467e+07 x 1E-38 cm^2, evaluated in 9.21263s 
1698243754 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=1.9307) = 5.2768e+07 x 1E-38 cm^2, evaluated in 9.3302s 
1698243763 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=2.13634) = 5.38308e+07 x 1E-38 cm^2, evaluated in 9.15457s 
1698243772 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=2.36389) = 5.48344e+07 x 1E-38 cm^2, evaluated in 9.07668s 
1698243781 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=2.61567) = 5.57759e+07 x 1E-38 cm^2, evaluated in 9.13001s 
1698243790 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=2.89427) = 5.66557e+07 x 1E-38 cm^2, evaluated in 9.28033s 
1698243799 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=3.20254) = 5.74756e+07 x 1E-38 cm^2, evaluated in 9.05951s 
1698243809 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=3.54365) = 5.82357e+07 x 1E-38 cm^2, evaluated in 9.24975s 
1698243818 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=3.92109) = 5.89392e+07 x 1E-38 cm^2, evaluated in 9.07545s 
1698243827 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=4.33873) = 5.95876e+07 x 1E-38 cm^2, evaluated in 9.26801s 
1698243836 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=4.80086) = 6.01848e+07 x 1E-38 cm^2, evaluated in 9.21062s 
1698243845 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=5.3122) = 6.07336e+07 x 1E-38 cm^2, evaluated in 9.15156s 
1698243855 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=5.87802) = 6.12369e+07 x 1E-38 cm^2, evaluated in 9.11365s 
1698243864 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=6.50409) = 6.16971e+07 x 1E-38 cm^2, evaluated in 9.12155s 
1698243873 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=7.19686) = 6.21182e+07 x 1E-38 cm^2, evaluated in 9.14525s 
1698243882 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=7.96341) = 6.2502e+07 x 1E-38 cm^2, evaluated in 9.07401s 
1698243891 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=8.8116) = 6.28532e+07 x 1E-38 cm^2, evaluated in 9.13889s 
1698243900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=9.75014) = 6.31729e+07 x 1E-38 cm^2, evaluated in 9.07094s 
1698243909 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=10.7886) = 6.34639e+07 x 1E-38 cm^2, evaluated in 8.98067s 
1698243918 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=11.9378) = 6.37281e+07 x 1E-38 cm^2, evaluated in 9.06008s 
1698243927 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=13.2093) = 6.39688e+07 x 1E-38 cm^2, evaluated in 9.01209s 
1698243936 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=14.6162) = 6.4188e+07 x 1E-38 cm^2, evaluated in 8.89121s 
1698243945 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=16.173) = 6.43859e+07 x 1E-38 cm^2, evaluated in 9.06964s 
1698243954 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=17.8956) = 6.45658e+07 x 1E-38 cm^2, evaluated in 9.02427s

which is 5 times more than for neutrinos.

jtenavidal commented 10 months ago

Test: change xsec shape by commenting out 1/Q4 dependence:

TICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0372759) = 0 x 1E-38 cm^2, evaluated in 4.291e-06s 
1698245236 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0412463) = 0 x 1E-38 cm^2, evaluated in 4.24e-06s 
1698245236 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0456395) = 0 x 1E-38 cm^2, evaluated in 4.256e-06s 
1698245236 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0505006) = 0 x 1E-38 cm^2, evaluated in 4.31e-06s 
1698245236 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0558795) = 0 x 1E-38 cm^2, evaluated in 4.324e-06s 
1698245236 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0618314) = 0 x 1E-38 cm^2, evaluated in 4.248e-06s 
1698245236 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0684171) = 0 x 1E-38 cm^2, evaluated in 4.345e-06s 
1698245236 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0757044) = 0 x 1E-38 cm^2, evaluated in 4.349e-06s 
1698245241 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0837678) = 0 x 1E-38 cm^2, evaluated in 5.47223s 
1698245247 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.09269) = 0 x 1E-38 cm^2, evaluated in 6.02697s 
1698245253 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.102563) = 0 x 1E-38 cm^2, evaluated in 6.13451s 
1698245260 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.113487) = 0 x 1E-38 cm^2, evaluated in 6.62378s 
1698245266 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.125574) = 0 x 1E-38 cm^2, evaluated in 6.61589s 
1698245273 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.13895) = 0 x 1E-38 cm^2, evaluated in 6.73813s 
1698245280 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.153749) = 0 x 1E-38 cm^2, evaluated in 6.88129s 
1698245287 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.170125) = 0 x 1E-38 cm^2, evaluated in 6.96975s 
1698245295 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.188246) = 2.18732e-10 x 1E-38 cm^2, evaluated in 7.50659s 
1698245302 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.208296) = 4.57181e-08 x 1E-38 cm^2, evaluated in 7.43542s 
1698245309 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.230482) = 5.38935e-07 x 1E-38 cm^2, evaluated in 7.41365s 
1698245317 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.255031) = 3.2389e-06 x 1E-38 cm^2, evaluated in 7.40036s 
1698245324 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.282195) = 1.46838e-05 x 1E-38 cm^2, evaluated in 7.44022s 
1698245332 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.312252) = 5.82326e-05 x 1E-38 cm^2, evaluated in 7.48901s 
1698245339 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.345511) = 0.000207954 x 1E-38 cm^2, evaluated in 7.50988s 
1698245347 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.382312) = 0.000613476 x 1E-38 cm^2, evaluated in 7.42539s 
1698245354 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.423032) = 0.00145019 x 1E-38 cm^2, evaluated in 7.52762s 
1698245362 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.46809) = 0.00287587 x 1E-38 cm^2, evaluated in 7.57736s 
1698245369 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.517947) = 0.00500026 x 1E-38 cm^2, evaluated in 7.61892s 
1698245377 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.573115) = 0.00786001 x 1E-38 cm^2, evaluated in 7.73012s 
1698245385 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.634158) = 0.011407 x 1E-38 cm^2, evaluated in 7.73342s

It seems not to affect the computing time

jtenavidal commented 10 months ago

Neutrino spline using ReinSehgalRESXSec:

1698245459 NOTICE XSecSplLst : [s] <CreateSpline (200)> : Energy threshold for current interaction = 0.110881 GeV
1698245459 NOTICE XSecSplLst : [s] <CreateSpline (256)> : xsec(E = 0.01) =  0 x 1E-38 cm^2, evaluated in 5.4099e-05 s
1698245459 NOTICE XSecSplLst : [s] <CreateSpline (256)> : xsec(E = 0.0301762) =  0 x 1E-38 cm^2, evaluated in 2.306e-06 s
1698245459 NOTICE XSecSplLst : [s] <CreateSpline (256)> : xsec(E = 0.0503525) =  0 x 1E-38 cm^2, evaluated in 1.363e-06 s
1698245459 NOTICE XSecSplLst : [s] <CreateSpline (256)> : xsec(E = 0.0705287) =  0 x 1E-38 cm^2, evaluated in 1.199e-06 s
1698245459 NOTICE XSecSplLst : [s] <CreateSpline (256)> : xsec(E = 0.090705) =  0 x 1E-38 cm^2, evaluated in 1.054e-06 s
1698245459 NOTICE XSecSplLst : [s] <CreateSpline (256)> : xsec(E = 0.110881) =  0 x 1E-38 cm^2, evaluated in 1.009e-06 s
1698245459 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (128)> : No cached RES v-production data for input neutrino (pdgc: 14)
1698245459 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (131)> : Wait while computing/caching RES production xsec first...
1698245459 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (110)> : 
 ** Creating cache branch - key = genie::BergerSehgalRESPXSec2014/NoPauliBlock/ResExcitationXSec/R:P33(1232);nu:14;int:Weak[CC]n
1698245459 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (118)> : E threshold = 0.110881
1698245459 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.01) = 0 x 1E-38 cm^2, evaluated in 6.01e-07s 
1698245459 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0301762) = 0 x 1E-38 cm^2, evaluated in 2.41e-07s 
1698245459 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0503525) = 0 x 1E-38 cm^2, evaluated in 2.18e-07s 
1698245459 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.0705287) = 0 x 1E-38 cm^2, evaluated in 2.07e-07s 
1698245459 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.090705) = 0 x 1E-38 cm^2, evaluated in 1.77e-07s 
1698245459 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.110881) = 0 x 1E-38 cm^2, evaluated in 1.78e-07s 

1698245475 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.120011) = 0 x 1E-38 cm^2, evaluated in 15.1933s 
1698245490 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.129892) = 0 x 1E-38 cm^2, evaluated in 15.1922s 
1698245505 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.140587) = 0 x 1E-38 cm^2, evaluated in 15.0089s 
1698245520 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.152162) = 0 x 1E-38 cm^2, evaluated in 14.9857s 
1698245535 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.16469) = 0 x 1E-38 cm^2, evaluated in 14.9505s 
1698245550 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.17825) = 0 x 1E-38 cm^2, evaluated in 14.928s 
1698245565 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.192927) = 0 x 1E-38 cm^2, evaluated in 14.9712s 
1698245580 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.208812) = 0 x 1E-38 cm^2, evaluated in 14.9592s 
1698245594 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.226004) = 0 x 1E-38 cm^2, evaluated in 14.7472s 
1698245609 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.244613) = 0 x 1E-38 cm^2, evaluated in 14.6613s 
1698245624 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.264753) = 0 x 1E-38 cm^2, evaluated in 14.7s 
1698245638 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.286552) = 4.89154e-07 x 1E-38 cm^2, evaluated in 14.4772s 
1698245652 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.310145) = 9.93385e-06 x 1E-38 cm^2, evaluated in 13.8802s 
1698245666 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.335682) = 6.00687e-05 x 1E-38 cm^2, evaluated in 13.71s 
1698245680 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.36332) = 0.000243198 x 1E-38 cm^2, evaluated in 13.6842s 
1698245693 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (181)> : RES XSec (R:P33(1232), E=0.393235) = 0.00082508 x 1E-38 cm^2, evaluated in 13.7725s

Same results but more than 10 times more to compute

jtenavidal commented 10 months ago

Test at constant E=2GeV: I did the test to compute it for a constant energy. In such case, the result is always the same and it takes the same amount of time:

1698251408 NOTICE XSecSplLst : [s] <CreateSpline (200)> : Energy threshold for current interaction = 0 GeV
1698251408 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (128)> : No cached RES v-production data for input neutrino (pdgc: 11)
1698251408 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (131)> : Wait while computing/caching RES production xsec first...
1698251408 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (110)> : 
 ** Creating cache branch - key = genie::BergerSehgalRESPXSec2014/EM-NoPauliBlock/ResExcitationXSec/R:P33(1232);nu:11;int:EMn
1698251408 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (118)> : E threshold = 0
1698251417 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.80129s 
1698251425 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.75868s 
1698251434 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.78768s 
1698251443 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.74729s 
1698251452 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.73148s 
1698251460 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.75798s 
1698251469 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.78451s 
1698251478 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.76787s 
1698251487 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.75775s 
1698251495 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.75706s 
1698251504 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.76795s 
1698251513 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.74977s 
1698251522 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.79459s 
1698251530 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.75448s 
1698251539 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.76685s 
1698251548 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.75748s 
1698251557 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.77577s 
1698251566 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.79081s 
1698251574 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (183)> : RES XSec (R:P33(1232), E=2) = 0.0542877 x 1E-38 cm^2, evaluated in 8.797s

Notice that for this calculation I didn't put back the 1/Q4 dependence hence the results are not the same. It takes more or less the same time as it did in the original calculation hence we can rule out it is due to a memory leak but the behaviour is as expected.

idkakorin commented 10 months ago

When using ReinSehgalRESXSecWithCacheFast: It holds after: 1698243378 NOTICE ReinSehgalResCF : [s] <CacheResExcitationXSec (192)> : RES XSec (R:P33(1232), E=0.187738) = 0 x 1E-38 cm^2, evaluated in 0.000163174 s

I suggest to play with parameters gsl-max-eval and gsl-relative-tolerance in the file ReinSehgalRESXSecFast.xml to find optimal ones. For example, the first attempt to set gsl-max-eval=10000000 gives the positive result: the calculation of the first non-zero node takes some time due to the peculiarities of the adaptive algorithm and way of calculation the cross section, but the remaining nodes are calculated quite quickly. Based on the results of these games, you can create parameter sets EM-NoPauliBlock and EM-WithPauliBlock with found optimal values for gsl-max-eval and gsl-relative-tolerance in the file ReinSehgalRESXSecFast.xml.