diff --git a/RecoVertex/AdaptiveVertexFinder/interface/TracksClusteringFromDisplacedSeed.h b/RecoVertex/AdaptiveVertexFinder/interface/TracksClusteringFromDisplacedSeed.h index 5e5705a038034..3b909a512add1 100644 --- a/RecoVertex/AdaptiveVertexFinder/interface/TracksClusteringFromDisplacedSeed.h +++ b/RecoVertex/AdaptiveVertexFinder/interface/TracksClusteringFromDisplacedSeed.h @@ -44,6 +44,8 @@ class TracksClusteringFromDisplacedSeed { std::pair,GlobalPoint> nearTracks(const reco::TransientTrack &seed, const std::vector & tracks, const reco::Vertex & primaryVertex) const; // unsigned int maxNTracks; + double max3DIPSignificance; + double max3DIPValue; double min3DIPSignificance; double min3DIPValue; double clusterMaxDistance; diff --git a/RecoVertex/AdaptiveVertexFinder/python/inclusiveVertexFinder_cfi.py b/RecoVertex/AdaptiveVertexFinder/python/inclusiveVertexFinder_cfi.py index fd9eace13e510..562b59d17f2c0 100644 --- a/RecoVertex/AdaptiveVertexFinder/python/inclusiveVertexFinder_cfi.py +++ b/RecoVertex/AdaptiveVertexFinder/python/inclusiveVertexFinder_cfi.py @@ -10,6 +10,8 @@ maxNTracks = cms.uint32(30), clusterizer = cms.PSet( + seedMax3DIPSignificance = cms.double(9999.), + seedMax3DIPValue = cms.double(5), seedMin3DIPSignificance = cms.double(1.2), seedMin3DIPValue = cms.double(0.005), clusterMaxDistance = cms.double(0.05), #500um diff --git a/RecoVertex/AdaptiveVertexFinder/src/TracksClusteringFromDisplacedSeed.cc b/RecoVertex/AdaptiveVertexFinder/src/TracksClusteringFromDisplacedSeed.cc index 7a607df67558c..8e3ca58afde10 100644 --- a/RecoVertex/AdaptiveVertexFinder/src/TracksClusteringFromDisplacedSeed.cc +++ b/RecoVertex/AdaptiveVertexFinder/src/TracksClusteringFromDisplacedSeed.cc @@ -4,6 +4,8 @@ TracksClusteringFromDisplacedSeed::TracksClusteringFromDisplacedSeed(const edm::ParameterSet ¶ms) : // maxNTracks(params.getParameter("maxNTracks")), + max3DIPSignificance(params.getParameter("seedMax3DIPSignificance")), + max3DIPValue(params.getParameter("seedMax3DIPValue")), min3DIPSignificance(params.getParameter("seedMin3DIPSignificance")), min3DIPValue(params.getParameter("seedMin3DIPValue")), clusterMaxDistance(params.getParameter("clusterMaxDistance")), @@ -92,7 +94,7 @@ std::vector TracksClusteringFromDisp std::vector seeds; for(std::vector::const_iterator it = selectedTracks.begin(); it != selectedTracks.end(); it++){ std::pair ip = IPTools::absoluteImpactParameter3D(*it,pv); - if(ip.first && ip.second.value() >= min3DIPValue && ip.second.significance() >= min3DIPSignificance) + if(ip.first && ip.second.value() >= min3DIPValue && ip.second.significance() >= min3DIPSignificance && ip.second.value() <= max3DIPValue && ip.second.significance() <= max3DIPSignificance) { #ifdef VTXDEBUG std::cout << "new seed " << it-selectedTracks.begin() << " ref " << it->trackBaseRef().key() << " " << ip.second.value() << " " << ip.second.significance() << " " << it->track().hitPattern().trackerLayersWithMeasurement() << " " << it->track().pt() << " " << it->track().eta() << std::endl; diff --git a/RecoVertex/AdaptiveVertexFinder/test/sequenceTest.py b/RecoVertex/AdaptiveVertexFinder/test/sequenceTest.py index 071eadbc92004..03d66da8cdb2a 100644 --- a/RecoVertex/AdaptiveVertexFinder/test/sequenceTest.py +++ b/RecoVertex/AdaptiveVertexFinder/test/sequenceTest.py @@ -44,7 +44,8 @@ ) # Other statements -process.GlobalTag.globaltag = 'POSTLS162_V1::All' +process.GlobalTag.globaltag = 'POSTLS161_V15::All' +#process.GlobalTag.globaltag = 'POSTLS162_V5::All' #process.GlobalTag.globaltag = 'START53_V26::All'