Skip to content

Commit

Permalink
implemented maxIPSig and maxIPValue to reject seed tracks from PU
Browse files Browse the repository at this point in the history
  • Loading branch information
pvmulder committed Mar 20, 2014
1 parent 3799d7b commit f8c799f
Show file tree
Hide file tree
Showing 4 changed files with 9 additions and 2 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -44,6 +44,8 @@ class TracksClusteringFromDisplacedSeed {
std::pair<std::vector<reco::TransientTrack>,GlobalPoint> nearTracks(const reco::TransientTrack &seed, const std::vector<reco::TransientTrack> & tracks, const reco::Vertex & primaryVertex) const;

// unsigned int maxNTracks;
double max3DIPSignificance;
double max3DIPValue;
double min3DIPSignificance;
double min3DIPValue;
double clusterMaxDistance;
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -4,6 +4,8 @@

TracksClusteringFromDisplacedSeed::TracksClusteringFromDisplacedSeed(const edm::ParameterSet &params) :
// maxNTracks(params.getParameter<unsigned int>("maxNTracks")),
max3DIPSignificance(params.getParameter<double>("seedMax3DIPSignificance")),
max3DIPValue(params.getParameter<double>("seedMax3DIPValue")),
min3DIPSignificance(params.getParameter<double>("seedMin3DIPSignificance")),
min3DIPValue(params.getParameter<double>("seedMin3DIPValue")),
clusterMaxDistance(params.getParameter<double>("clusterMaxDistance")),
Expand Down Expand Up @@ -92,7 +94,7 @@ std::vector<TracksClusteringFromDisplacedSeed::Cluster> TracksClusteringFromDisp
std::vector<TransientTrack> seeds;
for(std::vector<TransientTrack>::const_iterator it = selectedTracks.begin(); it != selectedTracks.end(); it++){
std::pair<bool,Measurement1D> 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;
Expand Down
3 changes: 2 additions & 1 deletion RecoVertex/AdaptiveVertexFinder/test/sequenceTest.py
Original file line number Diff line number Diff line change
Expand Up @@ -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'


Expand Down

0 comments on commit f8c799f

Please sign in to comment.