added tau eta option

PhysicsAnalysis/TopPhys/xAOD/TopCPTools/Root/TopTauCPTools.cxx

@@ -138,9 +138,10 @@ namespace top {
     int tauEleBDTWPLoose = tau_EleBDTWP_to_enum(m_config->tauEleBDTWPLoose());
     top::check(tauEleBDTWPLoose >= 0, m_config->tauEleBDTWPLoose() + " is not a valid tau WP");
 
-    const std::vector<float> absEtaRegion = {
-      0., 1.37, 1.37, 1.52, 1.52, 2.5
-    };
+//    const std::vector<float> absEtaRegion = {
+//      0., 1.37, 1.37, 1.52, 1.52, 2.5
+//    };
+
     const double absCharge = 1.;
     const std::vector<size_t> nTracks = {
       1, 3
@@ -218,7 +219,7 @@ namespace top {
                    "Failed to set tau SelectionCuts");
         top::check(asg::setProperty(tauSelectionTool, "PtMin", static_cast<double>(m_config->tauPtcut() / 1e3)),
                    "Failed to set tau pT cut");
-        top::check(asg::setProperty(tauSelectionTool, "AbsEtaRegion", absEtaRegion),
+        top::check(asg::setProperty(tauSelectionTool, "AbsEtaRegion", m_config->tauEtaRegions() ),
                    "Failed to set tau AbsEtaRegion");
         top::check(asg::setProperty(tauSelectionTool, "AbsCharge", absCharge),
                    "Failed to set tau AbsCharge");
@@ -287,7 +288,7 @@ namespace top {
                    "Failed to set tau SelectionCuts");
         top::check(asg::setProperty(tauSelectionTool, "PtMin", static_cast<double>(m_config->tauPtcut() / 1e3)),
                    "Failed to set tau pT cut");
-        top::check(asg::setProperty(tauSelectionTool, "AbsEtaRegion", absEtaRegion),
+        top::check(asg::setProperty(tauSelectionTool, "AbsEtaRegion", m_config->tauEtaRegions()),
                    "Failed to set loose tau AbsEtaRegion");
         top::check(asg::setProperty(tauSelectionTool, "AbsCharge", absCharge),
                    "Failed to set loose tau AbsCharge");

PhysicsAnalysis/TopPhys/xAOD/TopConfiguration/Root/ConfigurationSettings.cxx

@@ -265,6 +265,10 @@ namespace top {
                       "Pt cut applied to both tight and loose taus (in MeV)."
                       "Default 25 GeV.",
                       "25000");
+    registerParameter("TauEtaRegions",
+					  "Eta regions used for both tight and loose taus."
+					  "Default vetoing crack region [0., 1.37, 1.52, 2.5]",
+					  "[0., 1.37, 1.52, 2.5]");
     registerParameter("TauJetIDWP",
                       "Tau jet IDWP (None, Loose, Medium, Tight, LooseNotMedium, LooseNotTight, MediumNotTight, NotLoose, RNNLoose, RNNMedium, RNNTight)."
                       "Default RNNMedium.",

PhysicsAnalysis/TopPhys/xAOD/TopConfiguration/Root/TopConfig.cxx

@@ -1268,6 +1268,7 @@ namespace top {
 
     //tau configuration
     this->tauPtcut(std::stof(settings->value("TauPt")));
+    this->tauEtaRegions(settings->value("TauEtaRegions"));
     this->tauJetIDWP(settings->value("TauJetIDWP"));
     this->tauJetIDWPLoose(settings->value("TauJetIDWPLoose"));
     this->tauEleBDTWP(settings->value("TauEleBDTWP"));

PhysicsAnalysis/TopPhys/xAOD/TopConfiguration/TopConfiguration/TopConfig.h

@@ -1341,6 +1341,13 @@ namespace top {
       }
     }
 
+    inline virtual void tauEtaRegions(const std::string& etaRegions) {
+      if (!m_configFixed) {
+        m_tau_configuration.etaRegions = etaRegions;
+        m_tau_configuration_loose.etaRegions = etaRegions;
+      }
+    }
+
     inline virtual void tauJetIDWP(const std::string& s) {
       if (!m_configFixed) m_tau_configuration.jetIDWP = s;
     }
@@ -1394,6 +1401,10 @@ namespace top {
       return m_tau_configuration.pt;
     }
 
+    inline const std::string& tauEtaRegions() const {
+       return m_tau_configuration.etaRegions;
+    }
+
     inline const std::string& tauJetIDWP() const {
       return m_tau_configuration.jetIDWP;
     }
@@ -2245,6 +2256,7 @@ namespace top {
       bool eleOLR = false;
       // pT cut on taus
       float pt = 20000;
+      std::string etaRegions="0., 1.37, 1.52, 2.5";
     } m_tau_configuration, m_tau_configuration_loose;
 
     // photon configuration