diff --git a/MuonSpectrometer/MuonReconstruction/MuonSegmentMakers/MuonSegmentMakerTools/DCMathSegmentMaker/src/MuonClusterSegmentFinderTool.cxx b/MuonSpectrometer/MuonReconstruction/MuonSegmentMakers/MuonSegmentMakerTools/DCMathSegmentMaker/src/MuonClusterSegmentFinderTool.cxx index 1a73965426b340636e6fe95ada8967c14e84d74a..7a01d416776c4c88df5e9249cd55f64d7f3d3284 100755 --- a/MuonSpectrometer/MuonReconstruction/MuonSegmentMakers/MuonSegmentMakerTools/DCMathSegmentMaker/src/MuonClusterSegmentFinderTool.cxx +++ b/MuonSpectrometer/MuonReconstruction/MuonSegmentMakers/MuonSegmentMakerTools/DCMathSegmentMaker/src/MuonClusterSegmentFinderTool.cxx @@ -27,6 +27,7 @@ MuonClusterSegmentFinderTool::MuonClusterSegmentFinderTool(const std::string& ty // declareProperty("IPConstraint", m_ipConstraint = true); declareProperty("ClusterDistance", m_maxClustDist = 5); + declareProperty("NOfSeedLayers", m_nOfSeedLayers=1); } StatusCode @@ -547,11 +548,11 @@ MuonClusterSegmentFinderTool::segmentSeed(std::vector<std::vector<const Muon::Mu { 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 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; 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 const Amg::Vector3D& gp1 = (*cit)->prepRawData()->globalPosition(); 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; for (std::vector<const Muon::MuonClusterOnTrack*>::const_iterator cit2 = orderedClusters[k].begin();