search

e4nu / e4nuanalysiscode

Placeholder for all e4nu analyses
1 stars 2 forks source link

RES splines are empty #68

Closed jtenavidal closed 1 year ago

jtenavidal commented 1 year ago

The job logs show that it was computed correctly, but the output xml file has only 0s

jtenavidal commented 1 year ago

TEST:

  1. Running jobs with master + config dir
  2. Manually added RES splines from G18 into susa directory for Q2 min of 0.4. Can I compute the events now?
  3. Also running it for neutrinos. Xml files do not match the output of the jobs. Understand this - this is ok!
jtenavidal commented 1 year ago

Units

SLOG("XSecSplLst", pNOTICE)
                       << "xsec(E = " << E[i] << ") =  "
                       << (1E+38/units::cm2)*xsec[i] << " x 1E-38 cm^2";|

Where:

cm2 = 0.01/( hbarc*GeV) = 0.01/hbarc = 1.973269804e+14
hbarc         = 1.973269804e-16;
jtenavidal commented 1 year ago

Conversion works fine:

root [9] double hbarc= 1.973269804e-16;
root [10] double meter         = 1. / (hbarc * 1);
root [11] double centimeter  = 0.01*meter;
root [12] centimeter
(double) 5.0677307e+13
root [13]  double centimeter2 = centimeter*centimeter;
root [14] double cm2 = centimeter2;
root [15] 3.8579589407e-04*(1E+38/cm2)
(double) 15022096.
jtenavidal commented 1 year ago

Master version fails to produce the GEM21_11a EMRES splines. The values are non-null but too small - so it looks like 0 G18_01a seems ok

jtenavidal commented 1 year ago

GEM21a :

1690358756 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.188246) = 5.77137e-313 x 1E-38 cm^2
1690358786 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.208296) = 2.08924e-310 x 1E-38 cm^2
1690358819 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.230482) = 2.11752e-309 x 1E-38 cm^2
1690358850 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.255031) = 1.0537e-308 x 1E-38 cm^2
1690358883 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.282195) = 3.995e-308 x 1E-38 cm^2
1690358915 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.312252) = 1.35102e-307 x 1E-38 cm^2
1690358948 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.345511) = 4.04247e-307 x 1E-38 cm^2
1690358977 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.382312) = 8.96743e-307 x 1E-38 cm^2
1690359007 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.423032) = 1.46798e-306 x 1E-38 cm^2
1690359037 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.46809) = 2.00335e-306 x 1E-38 cm^2
1690359069 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.517947) = 2.47211e-306 x 1E-38 cm^2
1690359100 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.573115) = 2.87474e-306 x 1E-38 cm^2
1690359132 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.634158) = 3.22127e-306 x 1E-38 cm^2
1690359164 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.701704) = 3.52236e-306 x 1E-38 cm^2
jtenavidal commented 1 year ago

G18_01a

1690359052 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0757044) = 0 x 1E-38 cm^2
1690359072 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0837678) = 0.00262371 x 1E-38 cm^2
1690359092 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.09269) = 0.113033 x 1E-38 cm^2
1690359112 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.102563) = 0.958492 x 1E-38 cm^2
1690359132 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.113487) = 4.39894 x 1E-38 cm^2
1690359152 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.125574) = 14.6974 x 1E-38 cm^2
jtenavidal commented 1 year ago

Error seems to affect both G18_02a and the 10a series -> BS model

jtenavidal commented 1 year ago

Using v3.0.0:

 |---------------------------------> UseDRJoinScheme   [bool] [unlocked] [l] : 1
 |--------------------------------------------> Wcut [double] [unlocked] [l] : 1.7
 |-----------------------------------> WeinbergAngle [double] [unlocked] [l] : 0.501717
 |---------------------------------> XSec-Integrator    [alg] [unlocked] [l] : genie::ReinSehgalRESXSecFast/NoPauliBlock
 |-------------------------------------> minibooneGA   [bool] [unlocked] [l] : 1
 |-------------------------------------> minibooneGV   [bool] [unlocked] [l] : 1

But now the XSec-Integrator is : XSec-Integrator [alg] [unlocked] [l] : genie::ReinSehgalRESXSec/NoPauliBlock

Xsec != 0 in this version

jtenavidal commented 1 year ago

v3.2.2 ?

jtenavidal commented 1 year ago 
1690361275 WARN XSecSplLst : [s] <SplineExists (109)> : No splines for tune G18_02a_00_000 were found!
1690361275 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (128)> : No cached RES v-production data for input neutrino (pdgc: 11)
1690361275 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (131)> : Wait while computing/caching RES production xsec first...
1690361275 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (110)> : 
 ** Creating cache branch - key = genie::BergerSehgalRESPXSec2014/EM-NoPauliBlock/ResExcitationXSec/R:P33(1232);nu:11;int:EMp
1690361275 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (118)> : E threshold = 0.00115938
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.01) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0110651) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0122437) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0135478) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0149908) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0165875) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0183542) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0203092) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0224723) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0248659) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0275144) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.030445) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0336878) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0372759) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0412463) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0456395) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0505006) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0558795) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0618314) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0684171) = 0 x 1E-38 cm^2
1690361275 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0757044) = 0 x 1E-38 cm^2
1690361297 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0837678) = 0 x 1E-38 cm^2
1690361320 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.09269) = 0 x 1E-38 cm^2
1690361345 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.102563) = 0 x 1E-38 cm^2
1690361371 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.113487) = 0 x 1E-38 cm^2
1690361398 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.125574) = 0 x 1E-38 cm^2
1690361426 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.13895) = 0 x 1E-38 cm^2
1690361455 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.153749) = 0 x 1E-38 cm^2
1690361485 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.170125) = 0 x 1E-38 cm^2
1690361517 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.188246) = 4706.73 x 1E-38 cm^2
1690361551 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.208296) = 1.47797e+06 x 1E-38 cm^2
1690361581 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.230482) = 1.49569e+07 x 1E-38 cm^2
1690361611 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.255031) = 7.45273e+07 x 1E-38 cm^2
1690361640 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.282195) = 2.82431e+08 x 1E-38 cm^2
1690361671 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.312252) = 9.5494e+08 x 1E-38 cm^2
1690361702 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.345511) = 2.85694e+09 x 1E-38 cm^2
^C
-bash-4.2$ git branch
* (HEAD detached at 5421588d0)
  master

It worked in this version: https://github.com/GENIE-MC/Generator/commit/5421588d01882e6441dda5b06c5abd7e37437f4b

But values are too big?

Same version with tune G18_01: 

   |---o       Sigma_c++ (PDG code = 4222)

1690361900 WARN XSecSplLst : [s] <SplineExists (109)> : No splines for tune G18_01a_00_000 were found!
1690361900 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (128)> : No cached RES v-production data for input neutrino (pdgc: 11)
1690361900 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (131)> : Wait while computing/caching RES production xsec first...
1690361900 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (110)> : 
 ** Creating cache branch - key = genie::ReinSehgalRESPXSec/EM-NoPauliBlock/ResExcitationXSec/R:P33(1232);nu:11;int:EMp
1690361900 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (118)> : E threshold = 0.00115938
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.01) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0110651) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0122437) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0135478) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0149908) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0165875) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0183542) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0203092) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0224723) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0248659) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0275144) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.030445) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0336878) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0372759) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0412463) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0456395) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0505006) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0558795) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0618314) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0684171) = 0 x 1E-38 cm^2
1690361900 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0757044) = 0 x 1E-38 cm^2
1690361920 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0837678) = 0.00262371 x 1E-38 cm^2
1690361940 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.09269) = 0.113033 x 1E-38 cm^2
1690361959 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.102563) = 0.958492 x 1E-38 cm^2
1690361979 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.113487) = 4.39894 x 1E-38 cm^2
1690361999 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.125574) = 14.6974 x 1E-38 cm^2
1690362020 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.13895) = 40.6472 x 1E-38 cm^2
1690362040 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.153749) = 99.5613 x 1E-38 cm^2
1690362061 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.170125) = 224.962 x 1E-38 cm^2
1690362082 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.188246) = 489.949 x 1E-38 cm^2
1690362104 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.208296) = 3453.46 x 1E-38 cm^2
1690362126 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.230482) = 26873.3 x 1E-38 cm^2
1690362146 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.255031) = 127992 x 1E-38 cm^2
jtenavidal commented 1 year ago

EM QE with GEM21_11a:

1690272372 NOTICE XSecSplLst : [s] <CreateSpline (195)> : Energy threshold for current interaction = 0.000511138 GeV
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.01) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.0131796) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.0173701) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.022893) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.0301719) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.0397652) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.0524088) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.0690725) =  0 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.0910345) =  7.40934e+07 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.119979) =  2.74623e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.158127) =  4.76416e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.208405) =  6.37241e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.274668) =  7.51863e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.362001) =  8.28442e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.477101) =  8.77903e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.628798) =  9.09609e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 0.828727) =  9.30204e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 1.09223) =  9.43913e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 1.43951) =  9.53283e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 1.8972) =  9.59837e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 2.50043) =  9.64505e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 3.29546) =  9.67876e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 4.34326) =  9.70336e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 5.72423) =  9.72148e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 7.54428) =  9.73491e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 9.94302) =  9.74491e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 13.1045) =  9.7524e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 17.2711) =  9.75802e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 22.7625) =  9.76225e+08 x 1E-38 cm^2
1690272372 NOTICE XSecSplLst : [s] <CreateSpline (243)> : xsec(E = 30) =  9.76544e+08 x 1E-38 cm^2
jtenavidal commented 1 year ago

v R-3_02_02

1690361367 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (128)> : No cached RES v-production data for input neutrino (pdgc: 11)
1690361367 WARN ReinSehgalResT : [n] <ReinSehgalRESXSec.cxx::Integrate (131)> : Wait while computing/caching RES production xsec first...
1690361367 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (110)> : 
 ** Creating cache branch - key = genie::BergerSehgalRESPXSec2014/EM-NoPauliBlock/ResExcitationXSec/R:P33(1232);nu:11;int:EMp
1690361367 NOTICE ReinSehgalResC : [n] <ReinSehgalRESXSecWithCache.cxx::CacheResExcitationXSec (118)> : E threshold = 0.00115938
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.01) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0110651) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0122437) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0135478) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0149908) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0165875) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0183542) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0203092) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0224723) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0248659) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0275144) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.030445) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0336878) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0372759) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0412463) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0456395) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0505006) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0558795) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0618314) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0684171) = 0 x 1E-38 cm^2
1690361367 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0757044) = 0 x 1E-38 cm^2
1690361390 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.0837678) = 0 x 1E-38 cm^2
1690361415 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.09269) = 0 x 1E-38 cm^2
1690361441 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.102563) = 0 x 1E-38 cm^2
1690361469 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.113487) = 0 x 1E-38 cm^2
1690361497 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.125574) = 0 x 1E-38 cm^2
1690361527 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.13895) = 0 x 1E-38 cm^2
1690361558 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.153749) = 0 x 1E-38 cm^2
1690361589 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.170125) = 0 x 1E-38 cm^2
1690361620 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.188246) = 11.9386 x 1E-38 cm^2
1690361651 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.208296) = 3750.79 x 1E-38 cm^2
1690361683 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.230482) = 37986 x 1E-38 cm^2
1690361716 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.255031) = 189412 x 1E-38 cm^2
1690361748 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.282195) = 718384 x 1E-38 cm^2
1690361779 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.312252) = 2.43148e+06 x 1E-38 cm^2
1690361812 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.345511) = 7.28042e+06 x 1E-38 cm^2
1690361846 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.382312) = 1.61642e+07 x 1E-38 cm^2
1690361880 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.423032) = 2.64683e+07 x 1E-38 cm^2
1690361914 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.46809) = 3.61325e+07 x 1E-38 cm^2
1690361947 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.517947) = 4.45923e+07 x 1E-38 cm^2
1690361981 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.573115) = 5.18571e+07 x 1E-38 cm^2
1690362017 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.634158) = 5.81057e+07 x 1E-38 cm^2
1690362055 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.701704) = 6.35338e+07 x 1E-38 cm^2
1690362094 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.776444) = 6.83223e+07 x 1E-38 cm^2
1690362134 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.859144) = 7.26163e+07 x 1E-38 cm^2
1690362175 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=0.950653) = 7.65271e+07 x 1E-38 cm^2
1690362216 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=1.05191) = 8.01379e+07 x 1E-38 cm^2
1690362257 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=1.16395) = 8.32119e+07 x 1E-38 cm^2
1690362298 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=1.28792) = 8.57034e+07 x 1E-38 cm^2
1690362339 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=1.4251) = 8.80192e+07 x 1E-38 cm^2
1690362381 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=1.57689) = 9.02277e+07 x 1E-38 cm^2
1690362421 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=1.74485) = 9.23397e+07 x 1E-38 cm^2
1690362462 NOTICE ReinSehgalResC : [s] <CacheResExcitationXSec (179)> : RES XSec (R:P33(1232), E=1.9307) = 9.43515e+07 x 1E-38 cm^2
jtenavidal commented 1 year ago

From vR-3_2_2 to v3_4 :

v3_4 vs more recent ones:

Master version:

For neutrinos there's no change

QUESTIONS:

jtenavidal commented 1 year ago

v3.2.2 (original values) -> https://github.com/GENIE-MC/Generator/commit/2ce2ecef58ba0fd0814d16a21f3cb1e2e4bfce9b (original values) -> https://github.com/GENIE-MC/Generator/commit/ce771f2d5f1a5d542f21b73e189d1adec2795c5b (values with dipole, ok ) -> v3.4.0 (values from dipole) -> https://github.com/GENIE-MC/Generator/pull/297 (? as dipole values. Contains EM parameters separation) -> master (0 xsec values)

jtenavidal commented 1 year ago

Fixed. Pull request opened in GENIE-Generator https://github.com/GENIE-MC/Generator/pull/306

There was a loading missing