Implement values of exclusion regions for failed PMTs in D2 modules

Jürgen Knödlseder
1 parent b0570719
Showing 2 changed files with 72 additions and 3 deletions
inst/com/src/GCOMSupport.cpp
inst/com/src/GCOMSupport.hpp
@@ -131,3 +131,69 @@ double gammalib::com_energy2(const double&amp; energy, const double&amp; phigeo)
     // Return D2 energy deposit
     return e2;
 }
+
+
+/***********************************************************************//**
+ * @brief Return D2 module exclusion region X position
+ *
+ * @param[in] id2 Module identifier [0,...,13].
+ * @return Exclusion region X position (cm).
+ *
+ * Returns the D2 module exclusion region X position for a given module.
+ * The method does not check the validity of the module identifier.
+ *
+ * The values have been implemented from the MPE-FPM-7 file.
+ ***************************************************************************/
+const double& gammalib::com_exd2x(const int& id2)
+{
+    // Set D2 module exclusion region X positions
+    static const double exd2x[] = {  0.0,    0.0,    0.0,   0.0,   0.0,   0.0,   0.0,
+                                     0.0,    0.0,    0.0, -51.7,   0.0,   9.0, -34.7};
+
+    // Return value
+    return exd2x[id2];
+}
+
+
+/***********************************************************************//**
+ * @brief Return D2 module exclusion region Y position
+ *
+ * @param[in] id2 Module identifier [0,...,13].
+ * @return Exclusion region Y position (cm).
+ *
+ * Returns the D2 module exclusion region Y position for a given module.
+ * The method does not check the validity of the module identifier.
+ *
+ * The values have been implemented from the MPE-FPM-7 file.
+ ***************************************************************************/
+const double& gammalib::com_exd2y(const int& id2)
+{
+    // Set D2 module exclusion region Y positions
+    static const double exd2y[] = {  0.0,-41.254,    0.0,   0.0,   0.0,   0.0,   0.0,
+                                     0.0,    0.0,    0.0,   8.7,   0.0, +41.2, +49.0};
+
+    // Return value
+    return exd2y[id2];
+}
+
+
+/***********************************************************************//**
+ * @brief Return D2 module exclusion region radius
+ *
+ * @param[in] id2 Module identifier [0,...,13].
+ * @return Exclusion region radius (cm).
+ *
+ * Returns the D2 module exclusion region radius for a given module.
+ * The method does not check the validity of the module identifier.
+ *
+ * The values have been implemented from the MPE-FPM-7 file.
+ ***************************************************************************/
+const double& gammalib::com_exd2r(const int& id2)
+{
+    // Set D2 module exclusion region radii
+    static const double exd2r[] = {  0.0, 14.085,    0.0,   0.0,   0.0,   0.0,   0.0,
+                                     0.0,    0.0,    0.0,   9.0,   0.0,   9.0,   9.0};
+
+    // Return value
+    return exd2r[id2];
+}
@@ -38,9 +38,12 @@ class GSkyMap;
 /* __ Prototypes _________________________________________________________ */
 namespace gammalib {
-    void   com_wcs_mer2car(GSkyMap& map);
-    double com_energy1(const double& energy, const double& phigeo);
-    double com_energy2(const double& energy, const double& phigeo);
+    void          com_wcs_mer2car(GSkyMap& map);
+    double        com_energy1(const double& energy, const double& phigeo);
+    double        com_energy2(const double& energy, const double& phigeo);
+    const double& com_exd2x(const int& id2);
+    const double& com_exd2y(const int& id2);
+    const double& com_exd2r(const int& id2);
 }
 #endif /* GCOMSUPPORT_HPP */