Commit 72358246c9659ec5df0dcb1a334d931060ff8858

Authored by Jürgen Knödlseder
1 parent ed1c0c27

Use true energy instead of reconstructued energy in ctobssim area computation

Showing 1 changed file with 4 additions and 4 deletions
src/ctobssim/ctobssim.cpp
@@ -805,7 +805,7 @@ void ctobssim::simulate_source(GCTAObservation* obs, @@ -805,7 +805,7 @@ void ctobssim::simulate_source(GCTAObservation* obs,
805 } 805 }
806 806
807 // Determine simulation area 807 // Determine simulation area
808 - double area = get_area(obs, ereco_min, ereco_max); 808 + double area = get_area(obs, etrue_min, etrue_max);
809 809
810 // Log energy range and simulation area 810 // Log energy range and simulation area
811 if (logNormal()) { 811 if (logNormal()) {
@@ -1225,14 +1225,14 @@ GEbounds ctobssim::get_ebounds(const GEbounds& ebounds) const @@ -1225,14 +1225,14 @@ GEbounds ctobssim::get_ebounds(const GEbounds& ebounds) const
1225 * @brief Get simulation area (cm^2) 1225 * @brief Get simulation area (cm^2)
1226 * 1226 *
1227 * @param[in] obs Pointer on CTA observation. 1227 * @param[in] obs Pointer on CTA observation.
1228 - * @param[in] emin Minimum energy.  
1229 - * @param[in] emax Maximum energy. 1228 + * @param[in] emin Minimum true energy.
  1229 + * @param[in] emax Maximum true energy.
1230 * @return Simulation area (cm^2). 1230 * @return Simulation area (cm^2).
1231 * 1231 *
1232 * Get the simulation area for an energy interval in units of cm^2. This is 1232 * Get the simulation area for an energy interval in units of cm^2. This is
1233 * done by extracting the maximum effective area value within the energy 1233 * done by extracting the maximum effective area value within the energy
1234 * range [emin,emax] and by multiplying this value by 2 for security. The 1234 * range [emin,emax] and by multiplying this value by 2 for security. The
1235 - * effective area is sampled at 10 energy values within the energy interval. 1235 + * effective area is sampled at 100 energy values within the energy interval.
1236 ***************************************************************************/ 1236 ***************************************************************************/
1237 double ctobssim::get_area(GCTAObservation* obs, 1237 double ctobssim::get_area(GCTAObservation* obs,
1238 const GEnergy& emin, 1238 const GEnergy& emin,