Skip to content
Snippets Groups Projects
Commit 36f9ec8c authored by Stefano Rosati's avatar Stefano Rosati
Browse files

adding cuts for NSW hough

parent a0f337ad
6 merge requests!58791DataQualityConfigurations: Modify L1Calo config for web display,!46784MuonCondInterface: Enable thread-safety checking.,!46776Updated LArMonitoring config file for WD to match new files produced using MT,!45405updated ART test cron job,!42417Draft: DIRE and VINCIA Base Fragments for Pythia 8.3,!39435SWEEP of Fixes in NSW segment finding and in treatment of MM stereo clusters (Sweep of MR !38871, !39241 and !39423)
...@@ -1534,12 +1534,24 @@ namespace Muon { ...@@ -1534,12 +1534,24 @@ namespace Muon {
updateHits(hits,hough); updateHits(hits,hough);
} }
Identifier id_hit = hits.front()->tgc ? hits.front()->tgc->etaCluster.hitList.front()->identify() : hits.front()->prd->identify();
MuonHough::MuonLayerHoughSelector selectorLoose;
MuonHough::MuonLayerHoughSelector selector;
if ( m_idHelperSvc->issTgc(id_hit) || m_idHelperSvc->isMM(id_hit) ) {
selectorLoose = MuonHough::MuonLayerHoughSelector({std::make_pair(0,9.9)});
selector = MuonHough::MuonLayerHoughSelector({std::make_pair(0,13.9)});
}
else {
selectorLoose = m_selectorsLoose[hough.m_descriptor.chIndex];
selector = m_selectors[hough.m_descriptor.chIndex];
}
// Muon::MuonStationIndex::StIndex stIndex = Muon::MuonStationIndex::toStationIndex(hough.m_descriptor.chIndex); // Muon::MuonStationIndex::StIndex stIndex = Muon::MuonStationIndex::toStationIndex(hough.m_descriptor.chIndex);
unsigned int nmaxima = 0; unsigned int nmaxima = 0;
while( nmaxima < 5 ){ while( nmaxima < 5 ){
MuonHough::MuonLayerHough::Maximum maximum; MuonHough::MuonLayerHough::Maximum maximum;
if( hough.findMaximum( maximum, m_selectorsLoose[hough.m_descriptor.chIndex] ) ) { if( hough.findMaximum( maximum, selectorLoose ) ) {
hough.associateHitsToMaximum(maximum,hits); hough.associateHitsToMaximum(maximum,hits);
ATH_MSG_DEBUG("findMaxima: Found Eta Maximum " << nmaxima ATH_MSG_DEBUG("findMaxima: Found Eta Maximum " << nmaxima
<< " " << maximum.max << " " << maximum.max
...@@ -1579,7 +1591,7 @@ namespace Muon { ...@@ -1579,7 +1591,7 @@ namespace Muon {
if( nmdt > 0 || (nmm + nstgc) > 0) { if( nmdt > 0 || (nmm + nstgc) > 0) {
maxima.push_back( new MuonHough::MuonLayerHough::Maximum(maximum) ); maxima.push_back( new MuonHough::MuonLayerHough::Maximum(maximum) );
// add to seed list if // add to seed list if
if( maximum.max > m_selectors[hough.m_descriptor.chIndex].getCutValue(maximum.pos) ) seedMaxima.push_back(maxima.back()); if( maximum.max > selector.getCutValue(maximum.pos) ) seedMaxima.push_back(maxima.back());
++nmaxima; ++nmaxima;
} }
hough.fillLayer2(maximum.hits,true); hough.fillLayer2(maximum.hits,true);
......
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