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 91 forks source link

Failure to run gevgen with non-default tunes on Genie R-3_04_00 #386

Open CptBurtReynolds opened 5 months ago

CptBurtReynolds commented 5 months ago

When running the following gevgen command I run into a fatal error:

gevgen -n 20000 -p 14 -t 1000260560 -e 150 -r 1008 --seed 2989819 --cross-sections /home/archie/Documents/Simulation_Stuff/Genie_Cross_Sections/G2111a00000/data/gxspl-NUsmall.xml --tune G21_11a_00_000

With the console outputting the following after loading the splines:

1707909763 WARN Algorithm : [n] <Algorithm.cxx::FindConfig (157)> : No Configuration available for genie::EventGeneratorListAssembler/Default at the ConfigPool 1707909763 NOTICE Decay : [n] <Decayer.cxx::LoadConfig (221)> : Configured to decay: [-] |---o rho- (PDG code = -213) |---o rho0 (PDG code = 113) |---o rho+ (PDG code = 213) |---o eta (PDG code = 221) |---o omega (PDG code = 223) |---o eta' (PDG code = 331)

...

|---o Lambda_c+ (PDG code = 4122) |---o Sigma_c+ (PDG code = 4212) |---o Sigma_c++ (PDG code = 4222) 1707909763 FATAL Algorithm : [n] </home/archie/Documents/Software/Genie/src/Framework/Algorithm/Algorithm.icc::GetParam (102)> : Key: FermiConstant does not exist in pools from algorithm : genie::BergerSehgalRESPXSec2014/NoPauliBlock 1707909763 FATAL Algorithm : [n] </home/archie/Documents/Software/Genie/src/Framework/Algorithm/Algorithm.icc::GetParam (106)> : Current Configuration 1707909763 FATAL Algorithm : [n] </home/archie/Documents/Software/Genie/src/Framework/Algorithm/Algorithm.icc::GetParam (107)> : [-] Registry name: [genie::BergerSehgalRESPXSec2014/NoPauliBlock_summary] - Write Status: [unlocked] - Inhibited Item Locking: [on] - # entries: 65 |-----------------------------------------> CKM-Vcd [double] [unlocked] [l] : 0.22 |-----------------------------------------> CKM-Vcs [double] [unlocked] [l] : 0.995

The same occurs when I attempt to run gevgen with the tune G18_10a_02_11b. Running without specifying a tune (i.e using the default G18_02a_00_000) works perfectly fine, and generates the event file as expected.

nusense commented 5 months ago 
$ setup genie v3_04_00 -q e20:debug
$ setup genie_xsec v3_04_00 -q AR2320i00000:e1000:k250
$ echo $GENIEXSECFILE
/cvmfs/larsoft.opensciencegrid.org/products/genie_xsec/v3_04_00/NULL/AR2320i00000-k250-e1000/data/gxspl-NUsmall.xml
$ echo $GENIE_XSEC_TUNE
AR23_20i_00_000
# run fewer events
$ gevgen -n 200 -p 14 -t 1000260560 -e 150 -r 1008 --seed 2989819 --cross-sections \
     $GENIEXSECFILE --tune $GENIE_XSEC_TUNE

works fine, so it isn't just non-default tunes in general but must be specific to G21_11a_00_000.

Let me install the cross-section file I have for G21_11a_00_000 off of https://scisoft.fnal.gov/scisoft/packages/genie_xsec/v3_04_00/ 

$ setup genie_xsec v3_04_00 -q G2111a00000:k250:e1000
$ echo $GENIEXSECFILE
/grid/fermiapp/products/genie/local/genie_xsec/v3_04_00/NULL/G2111a00000-k250-e1000/data/gxspl-NUsmall.xml
$ echo $GENIE_XSEC_TUNE
G21_11a_00_000
$ gevgen -n 200 -p 14 -t 1000260560 -e 150 -r 1008 --seed 2989819 --cross-sections \
     $GENIEXSECFILE --tune $GENIE_XSEC_TUNE

This ran fine as well. Are you sure you're using the tagged version of R-3_04_00 without modifications? From whence did the spline file come? scisoft.fnal.gov?

CptBurtReynolds commented 5 months ago

Interesting. I should mention that G18_10a_02_11b also fails. What command/where can I see if it's the tagged version? The splines came from https://scisoft.fnal.gov/scisoft/packages/genie_xsec/v3_04_00/

karuboniru commented 5 months ago

What command/where can I see if it's the tagged version?

cd $GENIE ; git describe

CptBurtReynolds commented 5 months ago

Ah I think I figured it out. If I remember correctly there were issues with checking out the git repo so I might have grabbed it manually or something similar, but I'm pretty sure I'm running 3.6.0 given the large "in development banner" whenever I run anything.