From 370ee4fac44fde8cfb3f801f1e792aae67fd2591 Mon Sep 17 00:00:00 2001 From: Jürgen Knödlseder Date: Thu, 2 Jul 2015 17:01:25 +0200 Subject: [PATCH] Track the number of background events that fall outside the ROI. --- src/ctobssim/ctobssim.cpp | 16 +++++++++++++++- 1 file changed, 15 insertions(+), 1 deletion(-) diff --git a/src/ctobssim/ctobssim.cpp b/src/ctobssim/ctobssim.cpp index 3baf666..40ee0f0 100644 --- a/src/ctobssim/ctobssim.cpp +++ b/src/ctobssim/ctobssim.cpp @@ -1084,6 +1084,10 @@ void ctobssim::simulate_background(GCTAObservation* obs, // Reserves space for events events->reserve(list->size()+events->size()); + // Initialise statistics + int n_appended = 0; + int n_outside_roi = 0; + // Append events for (int k = 0; k < list->size(); k++) { @@ -1100,14 +1104,24 @@ void ctobssim::simulate_background(GCTAObservation* obs, // Append event events->append(*event); + // Increment number of appended events + n_appended++; + } // endif: event was within ROI + // ... otherwise increment outside ROI counter + else { + n_outside_roi++; + } + } // endfor: looped over all events // Dump simulation results if (logNormal()) { + *wrklog << gammalib::parformat("MC events outside ROI"); + *wrklog << n_outside_roi << std::endl; *wrklog << gammalib::parformat("MC background events"); - *wrklog << list->size() << std::endl; + *wrklog << n_appended << std::endl; } // Free event list -- libgit2 0.21.4