|
|
|
#include <PFMuonAlgo.h>
Public Member Functions | |
| PFMuonAlgo () | |
| constructor | |
| virtual | ~PFMuonAlgo () |
| destructor | |
Static Public Member Functions | |
| static bool | isGlobalLooseMuon (const reco::PFBlockElement &elt) |
| static bool | isGlobalLooseMuon (const reco::MuonRef &muonRef) |
| static bool | isGlobalTightMuon (const reco::PFBlockElement &elt) |
| static bool | isGlobalTightMuon (const reco::MuonRef &muonRef) |
| static bool | isIsolatedMuon (const reco::PFBlockElement &elt) |
| static bool | isIsolatedMuon (const reco::MuonRef &muonRef) |
| static bool | isLooseMuon (const reco::PFBlockElement &elt) |
| static bool | isLooseMuon (const reco::MuonRef &muonRef) |
| static bool | isMuon (const reco::MuonRef &muonRef) |
| static bool | isMuon (const reco::PFBlockElement &elt) |
| Check if a block element is a muon. | |
| static bool | isTightMuonPOG (const reco::MuonRef &muonRef) |
| static bool | isTrackerLooseMuon (const reco::MuonRef &muonRef) |
| static bool | isTrackerLooseMuon (const reco::PFBlockElement &elt) |
| static bool | isTrackerTightMuon (const reco::PFBlockElement &elt) |
| static bool | isTrackerTightMuon (const reco::MuonRef &muonRef) |
| static void | printMuonProperties (const reco::MuonRef &muonRef) |
Definition at line 7 of file PFMuonAlgo.h.
| PFMuonAlgo::PFMuonAlgo | ( | ) | [inline] |
| virtual PFMuonAlgo::~PFMuonAlgo | ( | ) | [inline, virtual] |
| bool PFMuonAlgo::isGlobalLooseMuon | ( | const reco::PFBlockElement & | elt | ) | [static] |
Definition at line 50 of file PFMuonAlgo.cc.
References reco::PFBlockElementTrack::muonRef().
Referenced by isLooseMuon(), and PFAlgoTestBenchElectrons::processBlock().
{
const reco::PFBlockElementTrack* eltTrack
= dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
assert ( eltTrack );
reco::MuonRef muonRef = eltTrack->muonRef();
return isGlobalLooseMuon(muonRef);
}
| bool PFMuonAlgo::isGlobalLooseMuon | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 194 of file PFMuonAlgo.cc.
References muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), min, query::result, and muon::TMLastStationAngTight.
{
if ( !muonRef.isNonnull() ) return false;
if ( !muonRef->isGlobalMuon() ) return false;
if ( !muonRef->isStandAloneMuon() ) return false;
reco::TrackRef standAloneMu = muonRef->standAloneMuon();
reco::TrackRef combinedMu = muonRef->combinedMuon();
reco::TrackRef trackerMu = muonRef->track();
unsigned nMuonHits =
standAloneMu->hitPattern().numberOfValidMuonDTHits() +
2*standAloneMu->hitPattern().numberOfValidMuonCSCHits();
bool quality = false;
if ( muonRef->isTrackerMuon() ){
bool result = combinedMu->normalizedChi2() < 100.;
bool laststation =
muon::isGoodMuon(*muonRef,muon::TMLastStationAngTight);
int nMatches = muonRef->numberOfMatches();
quality = laststation && nMuonHits > 12 && nMatches > 1;
return result && quality;
}
else{
// Check the quality of the stand-alone muon :
// good chi**2 and large number of hits and good pt error
if ( nMuonHits <=15 ||
standAloneMu->normalizedChi2() > 10. ||
standAloneMu->ptError()/standAloneMu->pt() > 0.20 ) {
quality = false;
}
else {
// If the stand-alone muon is good, check the global muon
if ( combinedMu->normalizedChi2() > standAloneMu->normalizedChi2() ) {
// If the combined muon is worse than the stand-alone, it
// means that either the corresponding tracker track was not
// reconstructed, or that the sta muon comes from a late
// pion decay (hence with a momentum smaller than the track)
// Take the stand-alone muon only if its momentum is larger
// than that of the track
// Note that here we even take the standAlone if it has a smaller pT, in contrast to GlobalTight
if(standAloneMu->pt() > trackerMu->pt() || combinedMu->normalizedChi2()<5.) quality = true;
}
else {
// If the combined muon is better (and good enough), take the
// global muon
if(combinedMu->ptError()/combinedMu->pt() < std::min(0.20,standAloneMu->ptError()/standAloneMu->pt()))
quality = true;
}
}
}
return quality;
}
| bool PFMuonAlgo::isGlobalTightMuon | ( | const reco::PFBlockElement & | elt | ) | [static] |
Definition at line 37 of file PFMuonAlgo.cc.
References reco::PFBlockElementTrack::muonRef().
Referenced by isMuon(), and PFAlgo::reconstructTrack().
{
const reco::PFBlockElementTrack* eltTrack
= dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
assert ( eltTrack );
reco::MuonRef muonRef = eltTrack->muonRef();
return isGlobalTightMuon(muonRef);
}
| bool PFMuonAlgo::isGlobalTightMuon | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 103 of file PFMuonAlgo.cc.
References muon::GlobalMuonPromptTight, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), min, query::result, and muon::TM2DCompatibilityTight.
{
if ( !muonRef.isNonnull() ) return false;
if ( !muonRef->isGlobalMuon() ) return false;
if ( !muonRef->isStandAloneMuon() ) return false;
if ( muonRef->isTrackerMuon() ) {
bool result = muon::isGoodMuon(*muonRef,muon::GlobalMuonPromptTight);
bool isTM2DCompatibilityTight = muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityTight);
int nMatches = muonRef->numberOfMatches();
bool quality = nMatches > 2 || isTM2DCompatibilityTight;
return result && quality;
} else {
reco::TrackRef standAloneMu = muonRef->standAloneMuon();
// No tracker muon -> Request a perfect stand-alone muon, or an even better global muon
bool result = false;
// Check the quality of the stand-alone muon :
// good chi**2 and large number of hits and good pt error
if ( ( standAloneMu->hitPattern().numberOfValidMuonDTHits() < 22 &&
standAloneMu->hitPattern().numberOfValidMuonCSCHits() < 15 ) ||
standAloneMu->normalizedChi2() > 10. ||
standAloneMu->ptError()/standAloneMu->pt() > 0.20 ) {
result = false;
} else {
reco::TrackRef combinedMu = muonRef->combinedMuon();
reco::TrackRef trackerMu = muonRef->track();
// If the stand-alone muon is good, check the global muon
if ( combinedMu->normalizedChi2() > standAloneMu->normalizedChi2() ) {
// If the combined muon is worse than the stand-alone, it
// means that either the corresponding tracker track was not
// reconstructed, or that the sta muon comes from a late
// pion decay (hence with a momentum smaller than the track)
// Take the stand-alone muon only if its momentum is larger
// than that of the track
result = standAloneMu->pt() > trackerMu->pt() ;
} else {
// If the combined muon is better (and good enough), take the
// global muon
result =
combinedMu->ptError()/combinedMu->pt() <
std::min(0.20,standAloneMu->ptError()/standAloneMu->pt());
}
}
return result;
}
return false;
}
| bool PFMuonAlgo::isIsolatedMuon | ( | const reco::PFBlockElement & | elt | ) | [static] |
Definition at line 76 of file PFMuonAlgo.cc.
References reco::PFBlockElementTrack::muonRef().
Referenced by GreedyMuonPFCandidateFilter::filter(), isMuon(), PFAlgo::postMuonCleaning(), PFAlgoTestBenchElectrons::processBlock(), and PFAlgo::reconstructTrack().
{
const reco::PFBlockElementTrack* eltTrack
= dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
assert ( eltTrack );
reco::MuonRef muonRef = eltTrack->muonRef();
return isIsolatedMuon(muonRef);
}
| bool PFMuonAlgo::isIsolatedMuon | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 285 of file PFMuonAlgo.cc.
References edm::Ref< C, T, F >::isNonnull().
{
if ( !muonRef.isNonnull() ) return false;
if ( !muonRef->isIsolationValid() ) return false;
// Isolated Muons which are missed by standard cuts are nearly always global+tracker
if ( !muonRef->isGlobalMuon() ) return false;
// If it's not a tracker muon, only take it if there are valid muon hits
reco::TrackRef standAloneMu = muonRef->standAloneMuon();
if ( !muonRef->isTrackerMuon() ){
if(standAloneMu->hitPattern().numberOfValidMuonDTHits() == 0 &&
standAloneMu->hitPattern().numberOfValidMuonCSCHits() ==0) return false;
}
// for isolation, take the smallest pt available to reject fakes
reco::TrackRef combinedMu = muonRef->combinedMuon();
double smallestMuPt = combinedMu->pt();
if(standAloneMu->pt()<smallestMuPt) smallestMuPt = standAloneMu->pt();
if(muonRef->isTrackerMuon())
{
reco::TrackRef trackerMu = muonRef->track();
if(trackerMu->pt() < smallestMuPt) smallestMuPt= trackerMu->pt();
}
double sumPtR03 = muonRef->isolationR03().sumPt;
double emEtR03 = muonRef->isolationR03().emEt;
double hadEtR03 = muonRef->isolationR03().hadEt;
double relIso = (sumPtR03 + emEtR03 + hadEtR03)/smallestMuPt;
if(relIso<0.1) return true;
else return false;
}
| bool PFMuonAlgo::isLooseMuon | ( | const reco::PFBlockElement & | elt | ) | [static] |
Definition at line 24 of file PFMuonAlgo.cc.
References reco::PFBlockElementTrack::muonRef().
{
const reco::PFBlockElementTrack* eltTrack
= dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
assert ( eltTrack );
reco::MuonRef muonRef = eltTrack->muonRef();
return isLooseMuon(muonRef);
}
| bool PFMuonAlgo::isLooseMuon | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 96 of file PFMuonAlgo.cc.
References isGlobalLooseMuon(), and isTrackerLooseMuon().
{
return isGlobalLooseMuon(muonRef) || isTrackerLooseMuon(muonRef);
}
| bool PFMuonAlgo::isMuon | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 90 of file PFMuonAlgo.cc.
References isGlobalTightMuon(), isIsolatedMuon(), and isTrackerTightMuon().
{
return isGlobalTightMuon(muonRef) || isTrackerTightMuon(muonRef) || isIsolatedMuon(muonRef);
}
| bool PFMuonAlgo::isMuon | ( | const reco::PFBlockElement & | elt | ) | [static] |
Check if a block element is a muon.
Definition at line 11 of file PFMuonAlgo.cc.
References reco::PFBlockElementTrack::muonRef().
{
const reco::PFBlockElementTrack* eltTrack
= dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
assert ( eltTrack );
reco::MuonRef muonRef = eltTrack->muonRef();
return isMuon(muonRef);
}
| bool PFMuonAlgo::isTightMuonPOG | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 327 of file PFMuonAlgo.cc.
References muon::GlobalMuonPromptTight, and muon::isGoodMuon().
{
if(!muon::isGoodMuon(*muonRef,muon::GlobalMuonPromptTight)) return false;
if(!muonRef->isTrackerMuon()) return false;
if(muonRef->numberOfMatches()<2) return false;
//const reco::TrackRef& combinedMuon = muonRef->combinedMuon();
const reco::TrackRef& combinedMuon = muonRef->globalTrack();
if(combinedMuon->hitPattern().numberOfValidTrackerHits()<11) return false;
if(combinedMuon->hitPattern().numberOfValidPixelHits()==0) return false;
if(combinedMuon->hitPattern().numberOfValidMuonHits()==0) return false;
return true;
}
| bool PFMuonAlgo::isTrackerLooseMuon | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 263 of file PFMuonAlgo.cc.
References muon::AllArbitrated, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), and muon::TMLastStationAngTight.
{
if ( !muonRef.isNonnull() ) return false;
if(!muonRef->isTrackerMuon()) return false;
reco::TrackRef trackerMu = muonRef->track();
if(trackerMu->ptError()/trackerMu->pt() > 0.20) return false;
// this doesn't seem to be necessary on the small samples looked at, but keep it around as insurance
if(trackerMu->pt()>20.) return false;
bool isAllArbitrated = muon::isGoodMuon(*muonRef,muon::AllArbitrated);
bool isTMLastStationAngTight = muon::isGoodMuon(*muonRef,muon::TMLastStationAngTight);
bool quality = isAllArbitrated && isTMLastStationAngTight;
return quality;
}
| static bool PFMuonAlgo::isTrackerLooseMuon | ( | const reco::PFBlockElement & | elt | ) | [static] |
Referenced by isLooseMuon().
| bool PFMuonAlgo::isTrackerTightMuon | ( | const reco::PFBlockElement & | elt | ) | [static] |
Definition at line 63 of file PFMuonAlgo.cc.
References reco::PFBlockElementTrack::muonRef().
Referenced by isMuon(), and PFAlgo::reconstructTrack().
{
const reco::PFBlockElementTrack* eltTrack
= dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
assert ( eltTrack );
reco::MuonRef muonRef = eltTrack->muonRef();
return isTrackerTightMuon(muonRef);
}
| bool PFMuonAlgo::isTrackerTightMuon | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 166 of file PFMuonAlgo.cc.
References muon::AllArbitrated, reco::TrackBase::hitPattern(), muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), and muon::TM2DCompatibilityTight.
{
if ( !muonRef.isNonnull() ) return false;
if(!muonRef->isTrackerMuon()) return false;
reco::TrackRef trackerMu = muonRef->track();
const reco::Track& track = *trackerMu;
unsigned nTrackerHits = track.hitPattern().numberOfValidTrackerHits();
if(nTrackerHits<=12) return false;
bool isAllArbitrated = muon::isGoodMuon(*muonRef,muon::AllArbitrated);
bool isTM2DCompatibilityTight = muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityTight);
if(!isAllArbitrated || !isTM2DCompatibilityTight) return false;
if((trackerMu->ptError()/trackerMu->pt() > 0.10)){
//std::cout<<" PT ERROR > 10 % "<< trackerMu->pt() <<std::endl;
return false;
}
return true;
}
| void PFMuonAlgo::printMuonProperties | ( | const reco::MuonRef & | muonRef | ) | [static] |
Definition at line 349 of file PFMuonAlgo.cc.
References runregparse::combined, gather_cfg::cout, delta, reco::TrackBase::hitPattern(), muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), max(), min, mathSSE::sqrt(), muon::TM2DCompatibilityLoose, muon::TM2DCompatibilityTight, muon::TMLastStationAngLoose, muon::TMLastStationAngTight, muon::TMLastStationLoose, muon::TMLastStationOptimizedBarrelLowPtLoose, muon::TMLastStationOptimizedBarrelLowPtTight, muon::TMLastStationOptimizedLowPtLoose, muon::TMLastStationOptimizedLowPtTight, muon::TMLastStationTight, muon::TMOneStationLoose, muon::TMOneStationTight, and patCandidatesForDimuonsSequences_cff::tracker.
Referenced by PFAlgoTestBenchElectrons::processBlock(), and PFAlgo::reconstructTrack().
{
if ( !muonRef.isNonnull() ) return;
bool isGL = muonRef->isGlobalMuon();
bool isTR = muonRef->isTrackerMuon();
bool isST = muonRef->isStandAloneMuon();
std::cout<<" GL: "<<isGL<<" TR: "<<isTR<<" ST: "<<isST<<std::endl;
std::cout<<" nMatches "<<muonRef->numberOfMatches()<<std::endl;
if ( muonRef->isGlobalMuon() ){
reco::TrackRef combinedMu = muonRef->combinedMuon();
std::cout<<" GL, pt: " << combinedMu->pt()
<< " +/- " << combinedMu->ptError()/combinedMu->pt()
<< " chi**2 GBL : " << combinedMu->normalizedChi2()<<std::endl;
std::cout<< " Total Muon Hits : " << combinedMu->hitPattern().numberOfValidMuonHits()
<< "/" << combinedMu->hitPattern().numberOfLostMuonHits()
<< " DT Hits : " << combinedMu->hitPattern().numberOfValidMuonDTHits()
<< "/" << combinedMu->hitPattern().numberOfLostMuonDTHits()
<< " CSC Hits : " << combinedMu->hitPattern().numberOfValidMuonCSCHits()
<< "/" << combinedMu->hitPattern().numberOfLostMuonCSCHits()
<< " RPC Hits : " << combinedMu->hitPattern().numberOfValidMuonRPCHits()
<< "/" << combinedMu->hitPattern().numberOfLostMuonRPCHits()<<std::endl;
std::cout<<" # of Valid Tracker Hits "<<combinedMu->hitPattern().numberOfValidTrackerHits()<<std::endl;
std::cout<<" # of Valid Pixel Hits "<<combinedMu->hitPattern().numberOfValidPixelHits()<<std::endl;
}
if ( muonRef->isStandAloneMuon() ){
reco::TrackRef standAloneMu = muonRef->standAloneMuon();
std::cout<<" ST, pt: " << standAloneMu->pt()
<< " +/- " << standAloneMu->ptError()/standAloneMu->pt()
<< " eta : " << standAloneMu->eta()
<< " DT Hits : " << standAloneMu->hitPattern().numberOfValidMuonDTHits()
<< "/" << standAloneMu->hitPattern().numberOfLostMuonDTHits()
<< " CSC Hits : " << standAloneMu->hitPattern().numberOfValidMuonCSCHits()
<< "/" << standAloneMu->hitPattern().numberOfLostMuonCSCHits()
<< " RPC Hits : " << standAloneMu->hitPattern().numberOfValidMuonRPCHits()
<< "/" << standAloneMu->hitPattern().numberOfLostMuonRPCHits()
<< " chi**2 STA : " << standAloneMu->normalizedChi2()<<std::endl;
}
if ( muonRef->isTrackerMuon() ){
reco::TrackRef trackerMu = muonRef->track();
const reco::Track& track = *trackerMu;
std::cout<<" TR, pt: " << trackerMu->pt()
<< " +/- " << trackerMu->ptError()/trackerMu->pt()
<< " chi**2 TR : " << trackerMu->normalizedChi2()<<std::endl;
std::cout<<" nTrackerHits "<<track.hitPattern().numberOfValidTrackerHits()<<std::endl;
std::cout<< "TMLastStationAngLoose "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationAngLoose) << std::endl
<< "TMLastStationAngTight "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationAngTight) << std::endl
<< "TMLastStationLoose "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationLoose) << std::endl
<< "TMLastStationTight "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationTight) << std::endl
<< "TMOneStationLoose "
<< muon::isGoodMuon(*muonRef,muon::TMOneStationLoose) << std::endl
<< "TMOneStationTight "
<< muon::isGoodMuon(*muonRef,muon::TMOneStationTight) << std::endl
<< "TMLastStationOptimizedLowPtLoose "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedLowPtLoose) << std::endl
<< "TMLastStationOptimizedLowPtTight "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedLowPtTight) << std::endl
<< "TMLastStationOptimizedBarrelLowPtLoose "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedBarrelLowPtLoose) << std::endl
<< "TMLastStationOptimizedBarrelLowPtTight "
<< muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedBarrelLowPtTight) << std::endl
<< std::endl;
}
std::cout<< "TM2DCompatibilityLoose "
<< muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityLoose) << std::endl
<< "TM2DCompatibilityTight "
<< muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityTight) << std::endl;
if ( muonRef->isGlobalMuon() && muonRef->isTrackerMuon() && muonRef->isStandAloneMuon() ){
reco::TrackRef combinedMu = muonRef->combinedMuon();
reco::TrackRef trackerMu = muonRef->track();
reco::TrackRef standAloneMu = muonRef->standAloneMuon();
double sigmaCombined = combinedMu->ptError()/(combinedMu->pt()*combinedMu->pt());
double sigmaTracker = trackerMu->ptError()/(trackerMu->pt()*trackerMu->pt());
double sigmaStandAlone = standAloneMu->ptError()/(standAloneMu->pt()*standAloneMu->pt());
bool combined = combinedMu->ptError()/combinedMu->pt() < 0.20;
bool tracker = trackerMu->ptError()/trackerMu->pt() < 0.20;
bool standAlone = standAloneMu->ptError()/standAloneMu->pt() < 0.20;
double delta1 = combined && tracker ?
fabs(1./combinedMu->pt() -1./trackerMu->pt())
/sqrt(sigmaCombined*sigmaCombined + sigmaTracker*sigmaTracker) : 100.;
double delta2 = combined && standAlone ?
fabs(1./combinedMu->pt() -1./standAloneMu->pt())
/sqrt(sigmaCombined*sigmaCombined + sigmaStandAlone*sigmaStandAlone) : 100.;
double delta3 = standAlone && tracker ?
fabs(1./standAloneMu->pt() -1./trackerMu->pt())
/sqrt(sigmaStandAlone*sigmaStandAlone + sigmaTracker*sigmaTracker) : 100.;
double delta =
standAloneMu->hitPattern().numberOfValidMuonDTHits()+
standAloneMu->hitPattern().numberOfValidMuonCSCHits() > 0 ?
std::min(delta3,std::min(delta1,delta2)) : std::max(delta3,std::max(delta1,delta2));
std::cout << "delta = " << delta << " delta1 "<<delta1<<" delta2 "<<delta2<<" delta3 "<<delta3<<std::endl;
double ratio =
combinedMu->ptError()/combinedMu->pt()
/ (trackerMu->ptError()/trackerMu->pt());
//if ( ratio > 2. && delta < 3. ) std::cout << "ALARM ! " << ratio << ", " << delta << std::endl;
std::cout<<" ratio "<<ratio<<" combined mu pt "<<combinedMu->pt()<<std::endl;
//bool quality3 = ( combinedMu->pt() < 50. || ratio < 2. ) && delta < 3.;
}
double sumPtR03 = muonRef->isolationR03().sumPt;
double emEtR03 = muonRef->isolationR03().emEt;
double hadEtR03 = muonRef->isolationR03().hadEt;
double relIsoR03 = (sumPtR03 + emEtR03 + hadEtR03)/muonRef->pt();
double sumPtR05 = muonRef->isolationR05().sumPt;
double emEtR05 = muonRef->isolationR05().emEt;
double hadEtR05 = muonRef->isolationR05().hadEt;
double relIsoR05 = (sumPtR05 + emEtR05 + hadEtR05)/muonRef->pt();
std::cout<<" 0.3 Radion Rel Iso: "<<relIsoR03<<" sumPt "<<sumPtR03<<" emEt "<<emEtR03<<" hadEt "<<hadEtR03<<std::endl;
std::cout<<" 0.5 Radion Rel Iso: "<<relIsoR05<<" sumPt "<<sumPtR05<<" emEt "<<emEtR05<<" hadEt "<<hadEtR05<<std::endl;
return;
}
1.7.3