Commit e17cd5ffa722469fcc51903c85e6504a13234aa6
1 parent
def46ac6
Fix show_spectrum.py script in case that no upper limits exist
Showing
1 changed file
with
6 additions
and
3 deletions
examples/show_spectrum.py
| @@ -240,11 +240,14 @@ def plot_spectrum(filename, plotfile): | @@ -240,11 +240,14 @@ def plot_spectrum(filename, plotfile): | ||
| 240 | plt.errorbar(spec['energies'], spec['flux'], | 240 | plt.errorbar(spec['energies'], spec['flux'], |
| 241 | yerr=spec['e_flux'], xerr=[spec['ed_engs'], spec['eu_engs']], | 241 | yerr=spec['e_flux'], xerr=[spec['ed_engs'], spec['eu_engs']], |
| 242 | fmt='ro') | 242 | fmt='ro') |
| 243 | + | ||
| 243 | # Plot upper limits | 244 | # Plot upper limits |
| 244 | - plt.errorbar(spec['ul_energies'], spec['ul_flux'], | ||
| 245 | - yerr=spec['yerr'], xerr=[spec['ul_ed_engs'], spec['ul_eu_engs']], | ||
| 246 | - uplims=True, fmt='ro') | 245 | + if len(spec['ul_energies']) > 0: |
| 246 | + plt.errorbar(spec['ul_energies'], spec['ul_flux'], yerr=spec['yerr'], | ||
| 247 | + xerr=[spec['ul_ed_engs'], spec['ul_eu_engs']], | ||
| 248 | + uplims=True, fmt='ro') | ||
| 247 | 249 | ||
| 250 | + # Plot log-likelihood profiles | ||
| 248 | if len(spec['dll_scan']) > 0: | 251 | if len(spec['dll_scan']) > 0: |
| 249 | plot_dloglike(spec) | 252 | plot_dloglike(spec) |
| 250 | 253 |