GENIE-MC/Generator

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

CptBurtReynolds opened this issue · 4 comments

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.

$ 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?

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/

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

cd $GENIE ; git describe

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.