enrico_sed - RuntimeError: Attempt to set bounds that exclude the existing value.

[Leonardo Supán]

Hi everyone!

I'm carrying out a binned analysis to get a SED, by running the enrico_sed script. After all the necessary initial steps and some iterations of gtlike, I was able to reach the point at which my energy interval is divided in NumEnergyBins bins of energy, but immediately the message I copy below pops up, killing the process.

As far as I could investigate, this error seems to be associated with trying to set/determine some values of the spectral index or its uncertainty, which are outside the interval you establish in your xml model (with min and max). Does anyone know if this is correct? I'm currently running enrico_sed, but setting a wider interval for gamma (that is, a lower min and a higher max), we'll see if this works...

Any suggestion will be very, very welcome and appreciated!! Thank you very much and all the best!

Leonardo

# ************************************************************

# *** Preparing submission of fit into energy bins ***

# ************************************************************

(' Emin = ', 500.0, ' Emax = ', 400000.0, ' Nbins = ', 8)

[Message]: Submition # 0 at energy 759

Traceback (most recent call last):

  File "/home/aips/Fermi/installation/enrico/bin/enrico_sed", line 58, in <module>

    sed(config,infile)

  File "/home/aips/Fermi/installation/enrico/bin/enrico_sed", line 13, in sed

    run(infile)

  File "/home/aips/Fermi/installation/enrico/enrico/RunGTlike.py", line 228, in run

    energybin.RunEbin(folder,Nbin,Fit,FitRunner)

  File "/home/aips/Fermi/installation/enrico/enrico/energybin.py", line 122, in RunEbin

    configfiles = PrepareEbin(Fit, FitRunner)

  File "/home/aips/Fermi/installation/enrico/enrico/energybin.py", line 98, in PrepareEbin

    NewFitObject = ChangeModel(Fit, ener[ibin], ener[ibin + 1], srcname, Pref[ibin] ,Gamma[ibin])

  File "/home/aips/Fermi/installation/enrico/enrico/energybin.py", line 29, in ChangeModel

    comp.logLike.getSource(name).getSrcFuncs()['Spectrum'].getParam('Index').setBounds(Gamma_min,Gamma_max)

  File "/home/aips/Fermi/installation/ScienceTools-v10r0p5-fssc-20150518-x86_64-unknown-linux-gnu-libc2.19-0/x86_64-unknown-linux-gnu-libc2.19-0/lib/python/pyLikelihood.py", line 1595, in setBounds

    return lib_pyLikelihood.Parameter_setBounds(*args)

RuntimeError: Attempt to set bounds that exclude the existing value.