Luigi Tibaldo / gammalib

Commit 1f0f05574fd174b7ba3a8aa9845508d3e327bf5a

Authored by Jürgen Knödlseder
1 parent e235d725

Make sure that COMPTEL background model is not negative

Inline Side-by-side
Showing 1 changed file with 27 additions and 0 deletions
inst/com/src/GCOMObservation.cpp
  View file @1f0f055
@@ -1271,6 +1271,15 @@ void GCOMObservation::compute_drb_phinor(const GCOMDri& drm) @@ -1271,6 +1271,15 @@ void GCOMObservation::compute_drb_phinor(const GCOMDri& drm)
1271 } 1271 }
1272 } 1272 }
1273 1273
  1274 + // Make sure that DRB is non-negative
  1275 + for (int iphibar = 0; iphibar < nphibar; ++iphibar) {
  1276 + for (int ipix = 0; ipix < npix; ++ipix) {
  1277 + if (map_drb(ipix, iphibar) < 0.0) {
  1278 + map_drb(ipix, iphibar) = 0.0;
  1279 + }
  1280 + }
  1281 + }
  1282 +
1274 // Return 1283 // Return
1275 return; 1284 return;
1276 } 1285 }
@@ -1617,6 +1626,15 @@ void GCOMObservation::compute_drb_bgdlixa(const GCOMDri&amp; drm, @@ -1617,6 +1626,15 @@ void GCOMObservation::compute_drb_bgdlixa(const GCOMDri&amp; drm,
1617 } 1626 }
1618 } 1627 }
1619 1628
  1629 + // Make sure that DRB is non-negative
  1630 + for (int iphibar = 0; iphibar < nphibar; ++iphibar) {
  1631 + for (int ipix = 0; ipix < npix; ++ipix) {
  1632 + if (map_drb(ipix, iphibar) < 0.0) {
  1633 + map_drb(ipix, iphibar) = 0.0;
  1634 + }
  1635 + }
  1636 + }
  1637 +
1620 // Return 1638 // Return
1621 return; 1639 return;
1622 } 1640 }
@@ -1821,6 +1839,15 @@ void GCOMObservation::compute_drb_bgdlixe(const GCOMDri&amp; drm, @@ -1821,6 +1839,15 @@ void GCOMObservation::compute_drb_bgdlixe(const GCOMDri&amp; drm,
1821 } 1839 }
1822 } 1840 }
1823 1841
  1842 + // Make sure that DRB is non-negative
  1843 + for (int iphibar = 0; iphibar < nphibar; ++iphibar) {
  1844 + for (int ipix = 0; ipix < npix; ++ipix) {
  1845 + if (map_drb(ipix, iphibar) < 0.0) {
  1846 + map_drb(ipix, iphibar) = 0.0;
  1847 + }
  1848 + }
  1849 + }
  1850 +
1824 // Return 1851 // Return
1825 return; 1852 return;
1826 } 1853 }