From 4177a13a367dd8b1d05fffeb62a02f88ea42f30c Mon Sep 17 00:00:00 2001
From: matthieu <matthieu.carrere@lupm.in2p3.fr>
Date: Tue, 14 Sep 2021 10:44:10 +0200
Subject: [PATCH] added new custom feature for telescopes array

---
 ...A_PROD3-Paranal-baseline-template-20deg-nopipe | 15 +++++++++++++--
 1 file changed, 13 insertions(+), 2 deletions(-)

diff --git a/CHERENKOV/examples_input/INPUTS_CTA_PROD3-Paranal-baseline-template-20deg-nopipe b/CHERENKOV/examples_input/INPUTS_CTA_PROD3-Paranal-baseline-template-20deg-nopipe
index d3a86f6..8eea7ae 100644
--- a/CHERENKOV/examples_input/INPUTS_CTA_PROD3-Paranal-baseline-template-20deg-nopipe
+++ b/CHERENKOV/examples_input/INPUTS_CTA_PROD3-Paranal-baseline-template-20deg-nopipe
@@ -75,7 +75,12 @@ NSHOW   2000                          // number of showers to generate
 #define DIFFUSE 1
 #endif
 *
+#ifdef THETA
+THETAP  $(THETA) $(THETA)      // Range of zenith angles (degree)
+#else
 THETAP  20.  20.      // Range of zenith angles (degree)
+#endif
+
 #ifdef FROM_SOUTH
 PHIP    0. 0.         // Range of azimuth angles (degree): primaries coming from South
 #else
@@ -102,7 +107,7 @@ NSHOW    $(NSHOW)            // Requires NSHOW environment variable.
 #endif
 #ifdef EMIN
 # ifdef EMAX
-ERANGE   $(EMIN)E3 $(EMAX)E3 // Requires EMIN and EMAX in units of TeV.
+ERANGE   $(EMIN) $(EMAX) // Requires EMIN and EMAX in units of GeV.
 # endif
 #endif
 #ifdef ESLOPE
@@ -120,7 +125,11 @@ ESLOPE $(ESLOPE)             // Requires spectral slope ESLOPE (<0)
 *   Maximum Cherenkov angle at ground level: 1.2014 deg.
 OBSLEV 2150.E2      // Observation level (in cm) for CTA near Paranal
 ATMOSPHERE 26 Y     // Should be slightly better for Paranal than profiles 1 (tropical) or 10 (HESS)
+#ifdef MAGNET_BX
+MAGNET $(MAGNET_BX) $(MAGNET_BZ)
+#else
 MAGNET 21.325 -8.926   // Magnetic field at assumed site [H, Z] (muT) (about the same as for Armazones site)
+#endif
 *
 ARRANG  0.            // Rotation of array to north [D] (degree); use zero here for any site for now
 *
@@ -132,7 +141,9 @@ CSCAT  10  2000e2  0. // Use shower several times (gammas, point source only)
 CSCAT  20  2500e2  0. // Use shower several times (protons+electrons+..., larger area for diffuse origin)
 #endif
 #ifdef NSCAT
-# ifdef CSCAT
+# ifdef CSCAT2
+CSCAT    $(NSCAT) $(CSCAT)E2 $(CSCAT2)E2 // Requires NSCAT (10/20) and CSCAT and CSCAT2 (in units of meters).
+#else
 CSCAT    $(NSCAT) $(CSCAT)E2 0. // Requires NSCAT (10/20) and CSCAT (in units of meters).
 # endif
 #endif
-- 
GitLab