Skip to content
Snippets Groups Projects
Commit f2ff6b8d authored by sutt's avatar sutt
Browse files

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
parent fafe2b8f
No related branches found
No related tags found
No related merge requests found
...@@ -535,6 +535,8 @@ int main(int argc, char** argv) ...@@ -535,6 +535,8 @@ int main(int argc, char** argv)
int nbl = -1; int nbl = -1;
int nsiholes = 2; int nsiholes = 2;
int npixholes = 20; /// essentially no limit
int nsctholes = 20; /// essentially no limit
bool expectBL = false; bool expectBL = false;
...@@ -590,6 +592,8 @@ int main(int argc, char** argv) ...@@ -590,6 +592,8 @@ int main(int argc, char** argv)
if ( inputdata.isTagDefined("zed") ) zed = inputdata.GetValue("zed"); if ( inputdata.isTagDefined("zed") ) zed = inputdata.GetValue("zed");
if ( inputdata.isTagDefined("npix") ) npix = inputdata.GetValue("npix"); if ( inputdata.isTagDefined("npix") ) npix = inputdata.GetValue("npix");
if ( inputdata.isTagDefined("nsiholes") ) nsiholes = inputdata.GetValue("nsiholes"); 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("expectBL") ) expectBL = ( inputdata.GetValue("expectBL") > 0.5 ? true : false );
if ( inputdata.isTagDefined("nsct") ) nsct = inputdata.GetValue("nsct"); if ( inputdata.isTagDefined("nsct") ) nsct = inputdata.GetValue("nsct");
if ( inputdata.isTagDefined("nbl") ) nbl = inputdata.GetValue("nbl"); if ( inputdata.isTagDefined("nbl") ) nbl = inputdata.GetValue("nbl");
...@@ -930,7 +934,7 @@ int main(int argc, char** argv) ...@@ -930,7 +934,7 @@ int main(int argc, char** argv)
Filter_Track filter_offline( eta, 1000, zed, pT, Filter_Track filter_offline( eta, 1000, zed, pT,
npix, nsct, -1, nbl, npix, nsct, -1, nbl,
-2, -2, chi2prob, -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 ( selectcharge!=0 ) filter_offline.chargeSelection( selectcharge );
if ( pTMax>pT ) filter_offline.maxpT( pTMax ); if ( pTMax>pT ) filter_offline.maxpT( pTMax );
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment