Skip to content

Commit

Permalink
added cut that each tau has at least one pixel hits (for at least one…
Browse files Browse the repository at this point in the history
… of the signal tracks),

in order to reduce increase of jet -> tau fake-rate observed for high pT jets in the Phys14 samples
  • Loading branch information
veelken committed Mar 5, 2015
1 parent eddeddf commit 895aafd
Show file tree
Hide file tree
Showing 2 changed files with 26 additions and 1 deletion.
Original file line number Diff line number Diff line change
Expand Up @@ -54,6 +54,8 @@ class PFRecoTauDiscriminationByHPSSelection : public PFTauDiscriminationProducer

bool requireTauChargedHadronsToBeChargedPFCands_;

int minPixelHits_;

int verbosity_;
};

Expand Down Expand Up @@ -103,6 +105,7 @@ PFRecoTauDiscriminationByHPSSelection::PFRecoTauDiscriminationByHPSSelection(con
));
}
requireTauChargedHadronsToBeChargedPFCands_ = pset.getParameter<bool>("requireTauChargedHadronsToBeChargedPFCands");
minPixelHits_ = pset.getParameter<int>("minPixelHits");
verbosity_ = pset.exists("verbosity") ?
pset.getParameter<int>("verbosity") : 0;
}
Expand Down Expand Up @@ -284,6 +287,26 @@ PFRecoTauDiscriminationByHPSSelection::discriminate(const reco::PFTauRef& tau) c
}
}
}

if ( minPixelHits_ > 0 ) {
int numPixelHits = 0;
const std::vector<reco::PFCandidatePtr>& chargedHadrCands = tau->signalPFChargedHadrCands();
for ( std::vector<reco::PFCandidatePtr>::const_iterator chargedHadrCand = chargedHadrCands.begin();
chargedHadrCand != chargedHadrCands.end(); ++chargedHadrCand ) {
const reco::Track* track = 0;
if ( (*chargedHadrCand)->trackRef().isNonnull() ) track = (*chargedHadrCand)->trackRef().get();
else if ( (*chargedHadrCand)->gsfTrackRef().isNonnull() ) track = (*chargedHadrCand)->gsfTrackRef().get();
if ( track ) {
numPixelHits += track->hitPattern().numberOfValidPixelHits();
}
}
if ( !(numPixelHits >= minPixelHits_) ) {
if ( verbosity_ ) {
edm::LogPrint("PFTauByHPSSelect") << " fails cut on sum of pixel hits." ;
}
return 0.0;
}
}

// Otherwise, we pass!
if ( verbosity_ ) {
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -75,7 +75,9 @@
decayMode_2Prong1Pi0,
decayMode_3Prong0Pi0
),
requireTauChargedHadronsToBeChargedPFCands = cms.bool(False)
requireTauChargedHadronsToBeChargedPFCands = cms.bool(False),
# CV: require at least one pixel hit for the sum of all tracks
minPixelHits = cms.int32(1)
)


Expand Down

0 comments on commit 895aafd

Please sign in to comment.