From f2ff6b8d13eedf33e8d73a76faea825de350a18e Mon Sep 17 00:00:00 2001
From: sutt <sutt@cern.ch>
Date: Tue, 16 Jun 2020 15:58:12 +0100
Subject: [PATCH] Allow rejection based on pixel or sct holes individually
Previous the standard selection would select only on combined silicon holes,
and speciceal selections were implemented that also selected individually on
pixel or sct holes. These changes allow both to be set also for the standard
analyses
---
.../TrigInDetAnalysisUser/Analysis/src/rmain.cxx | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/rmain.cxx b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/rmain.cxx
index e3305778040..6f9eec0452c 100644
--- a/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/rmain.cxx
+++ b/Trigger/TrigAnalysis/TrigInDetAnalysisUser/Analysis/src/rmain.cxx
@@ -535,6 +535,8 @@ int main(int argc, char** argv)
int nbl = -1;
int nsiholes = 2;
+ int npixholes = 20; /// essentially no limit
+ int nsctholes = 20; /// essentially no limit
bool expectBL = false;
@@ -590,6 +592,8 @@ int main(int argc, char** argv)
if ( inputdata.isTagDefined("zed") ) zed = inputdata.GetValue("zed");
if ( inputdata.isTagDefined("npix") ) npix = inputdata.GetValue("npix");
if ( inputdata.isTagDefined("nsiholes") ) nsiholes = inputdata.GetValue("nsiholes");
+ if ( inputdata.isTagDefined("npixholes") ) npixholes = inputdata.GetValue("npixholes");
+ if ( inputdata.isTagDefined("nsctholes") ) nsctholes = inputdata.GetValue("nsctholes");
if ( inputdata.isTagDefined("expectBL") ) expectBL = ( inputdata.GetValue("expectBL") > 0.5 ? true : false );
if ( inputdata.isTagDefined("nsct") ) nsct = inputdata.GetValue("nsct");
if ( inputdata.isTagDefined("nbl") ) nbl = inputdata.GetValue("nbl");
@@ -930,7 +934,7 @@ int main(int argc, char** argv)
Filter_Track filter_offline( eta, 1000, zed, pT,
npix, nsct, -1, nbl,
-2, -2, chi2prob,
- 20, 20, nsiholes, expectBL ); /// include chi2 probability cut
+ npixholes, nsctholes, nsiholes, expectBL ); /// include chi2 probability cut
if ( selectcharge!=0 ) filter_offline.chargeSelection( selectcharge );
if ( pTMax>pT ) filter_offline.maxpT( pTMax );
--
GitLab