MEIntegrator crashes on initPDFSet

jpata commented 10 years ago

initPDFSet segfaults with

#5  0x00007f32b0fd9e79 in _wordcopy_fwd_aligned () from /lib64/libc.so.6
#6  0x00007f32b0fd3c72 in memmove () from /lib64/libc.so.6
#7  0x0000000000462793 in master.1.evolvepdfm_ ()
#8  0x0000000000466b04 in readevolve_ ()
#9  0x000000000045a5a9 in master.1.initevolve_.part.0 ()
#10 0x0000000000454eda in master.0.initpdfsetm_ ()
#11 0x00000000004557f2 in initpdfsetbynamem_ ()
#12 0x000000000044ff7d in LHAPDF::initPDFSet(int, std::string const&, int) ()
#13 0x00000000004277f4 in MEIntegratorNew::MEIntegratorNew(std::string, int, int) ()
#14 0x0000000000418a36 in main ()

when called as

LHAPDF::initPDFSet(0, "cteq65.LHgrid");

from https://github.com/bianchini/TTHStudies/blob/b5fe96ee63a432a99d9421ced326b5c6294341ec/interface/MEIntegratorNew.h#L339

but runs when called as in https://github.com/jpata/tthbb13/commit/77b2fa30e887756749b06add5e71fe886892425b#L340

LHAPDF documentation is unclear as to the indexing of the sets (first argument): http://lhapdf.hepforge.org/lhapdf5/manual#tth_sEcA

Also happens when using the LHAPATH from CMSSW5:

LHAPATH=/cvmfs/cms.cern.ch/slc5_amd64_gcc462/external/lhapdf/5.8.5-cms/share/lhapdf/PDFsets bin/slc6_amd64_gcc481/MEAnalysis src/TTH/MEAnalysis/python/test_cfg.py

jpata commented 10 years ago

old gfortran [cms:CMSSW_5_3_3_patch2@]$ gfortran --version GNU Fortran (GCC) 4.6.2

new gfortran [thebe:TTH@master]$ gfortran --version GNU Fortran (GCC) 4.8.1

jpata commented 9 years ago

Does not crash with lhapdf (as opposed to lhapdffull). In any case, this index was only a map index and should not have a physical impact.

jpata / tthbb13

MEIntegrator crashes on initPDFSet #1