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

adding configurable number of seed layers

parent 8275197b
No related branches found
No related tags found
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)
...@@ -27,6 +27,7 @@ MuonClusterSegmentFinderTool::MuonClusterSegmentFinderTool(const std::string& ty ...@@ -27,6 +27,7 @@ MuonClusterSegmentFinderTool::MuonClusterSegmentFinderTool(const std::string& ty
// //
declareProperty("IPConstraint", m_ipConstraint = true); declareProperty("IPConstraint", m_ipConstraint = true);
declareProperty("ClusterDistance", m_maxClustDist = 5); declareProperty("ClusterDistance", m_maxClustDist = 5);
declareProperty("NOfSeedLayers", m_nOfSeedLayers=1);
} }
StatusCode StatusCode
...@@ -547,11 +548,11 @@ MuonClusterSegmentFinderTool::segmentSeed(std::vector<std::vector<const Muon::Mu ...@@ -547,11 +548,11 @@ MuonClusterSegmentFinderTool::segmentSeed(std::vector<std::vector<const Muon::Mu
{ {
std::vector<std::pair<Amg::Vector3D, Amg::Vector3D> > seeds; std::vector<std::pair<Amg::Vector3D, Amg::Vector3D> > seeds;
if (orderedClusters.size() < 2) return seeds; if (orderedClusters.size() < 4) return seeds;
// calculate the straight line between the two furthest points // calculate the straight line between the two furthest points
int seedlayers1 = 0; int seedlayers1 = 0;
for (unsigned int i = 0; (i < orderedClusters.size() && seedlayers1 < 2); ++i) { for (unsigned int i = 0; (i < orderedClusters.size() && seedlayers1 < m_nOfSeedLayers); ++i) {
bool usedLayer1 = false; bool usedLayer1 = false;
for (std::vector<const Muon::MuonClusterOnTrack*>::const_iterator cit = orderedClusters[i].begin(); for (std::vector<const Muon::MuonClusterOnTrack*>::const_iterator cit = orderedClusters[i].begin();
...@@ -563,7 +564,7 @@ MuonClusterSegmentFinderTool::segmentSeed(std::vector<std::vector<const Muon::Mu ...@@ -563,7 +564,7 @@ MuonClusterSegmentFinderTool::segmentSeed(std::vector<std::vector<const Muon::Mu
const Amg::Vector3D& gp1 = (*cit)->prepRawData()->globalPosition(); const Amg::Vector3D& gp1 = (*cit)->prepRawData()->globalPosition();
int seedlayers2 = 0; int seedlayers2 = 0;
for (unsigned int k = orderedClusters.size() - 1; (k > i && seedlayers2 < 2); --k) { for (unsigned int k = orderedClusters.size() - 1; (k > i && seedlayers2 < m_nOfSeedLayers); --k) {
bool usedLayer2 = false; bool usedLayer2 = false;
for (std::vector<const Muon::MuonClusterOnTrack*>::const_iterator cit2 = orderedClusters[k].begin(); for (std::vector<const Muon::MuonClusterOnTrack*>::const_iterator cit2 = orderedClusters[k].begin();
......
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