CMS 3D CMS Logo
Classes | Enumerations | Functions | Variables
muon Namespace Reference

Classes

struct  SelectionTypeStringToEnum
 a lightweight "map" for selection type string label and enum value More...
 
struct  SelectorStringToEnum
 

Enumerations

enum  AlgorithmType {
  TMLastStation, TM2DCompatibility, TMOneStation, RPCMu,
  ME0Mu, GEMMu
}
 
enum  SelectionType {
  All = 0, AllGlobalMuons = 1, AllStandAloneMuons = 2, AllTrackerMuons = 3,
  TrackerMuonArbitrated = 4, AllArbitrated = 5, GlobalMuonPromptTight = 6, TMLastStationLoose = 7,
  TMLastStationTight = 8, TM2DCompatibilityLoose = 9, TM2DCompatibilityTight = 10, TMOneStationLoose = 11,
  TMOneStationTight = 12, TMLastStationOptimizedLowPtLoose = 13, TMLastStationOptimizedLowPtTight = 14, GMTkChiCompatibility = 15,
  GMStaChiCompatibility = 16, GMTkKinkTight = 17, TMLastStationAngLoose = 18, TMLastStationAngTight = 19,
  TMOneStationAngLoose = 20, TMOneStationAngTight = 21, TMLastStationOptimizedBarrelLowPtLoose = 22, TMLastStationOptimizedBarrelLowPtTight = 23,
  RPCMuLoose = 24, AllME0Muons = 25, ME0MuonArbitrated = 26, AllGEMMuons = 27,
  GEMMuonArbitrated = 28, TriggerIdLoose = 29
}
 Selector type. More...
 

Functions

float caloCompatibility (const reco::Muon &muon)
 
reco::TrackRef getTevRefitTrack (const reco::TrackRef &combinedTrack, const reco::TrackToTrackMap &map)
 
bool isGoodMuon (const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
 main GoodMuon wrapper call More...
 
bool isGoodMuon (const reco::Muon &muon, AlgorithmType type, double minCompatibility, reco::Muon::ArbitrationType arbitrationType)
 
bool isGoodMuon (const reco::Muon &muon, AlgorithmType type, int minNumberOfMatches, double maxAbsDx, double maxAbsPullX, double maxAbsDy, double maxAbsPullY, double maxChamberDist, double maxChamberDistPull, reco::Muon::ArbitrationType arbitrationType, bool syncMinNMatchesNRequiredStationsInBarrelOnly=true, bool applyAlsoAngularCuts=false)
 
bool isHighPtMuon (const reco::Muon &, const reco::Vertex &)
 
bool isLooseMuon (const reco::Muon &)
 
bool isLooseTriggerMuon (const reco::Muon &)
 
bool isMediumMuon (const reco::Muon &, bool run2016_hip_mitigation=false)
 
bool isSoftMuon (const reco::Muon &, const reco::Vertex &, bool run2016_hip_mitigation=false)
 
bool isTightMuon (const reco::Muon &, const reco::Vertex &)
 
bool isTrackerHighPtMuon (const reco::Muon &, const reco::Vertex &)
 
reco::Muon::Selector makeSelectorBitset (reco::Muon const &muon, reco::Vertex const *vertex=nullptr, bool run2016_hip_mitigation=false)
 
bool overlap (const reco::Muon &muon1, const reco::Muon &muon2, double pullX=1.0, double pullY=1.0, bool checkAdjacentChambers=false)
 
unsigned int RequiredStationMask (const reco::Muon &muon, double maxChamberDist, double maxChamberDistPull, reco::Muon::ArbitrationType arbitrationType)
 
float segmentCompatibility (const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
 
SelectionType selectionTypeFromString (const std::string &label)
 
reco::Muon::Selector selectorFromString (const std::string &label)
 
int sharedSegments (const reco::Muon &muon1, const reco::Muon &muon2, unsigned int segmentArbitrationMask=reco::MuonSegmentMatch::BestInChamberByDR)
 
reco::Muon::MuonTrackTypePair sigmaSwitch (const reco::TrackRef &combinedTrack, const reco::TrackRef &trackerTrack, const double nSigma=2., const double ptThreshold=200.)
 
reco::Muon::MuonTrackTypePair sigmaSwitch (const reco::Muon &muon, const double nSigma=2., const double ptThreshold=200.)
 
reco::Muon::MuonTrackTypePair tevOptimized (const reco::TrackRef &combinedTrack, const reco::TrackRef &trackerTrack, const reco::TrackRef &tpfmsTrack, const reco::TrackRef &pickyTrack, const reco::TrackRef &dytTrack, const double ptThreshold=200., const double tune1=17., const double tune2=40., const double dptcut=0.25)
 
reco::Muon::MuonTrackTypePair tevOptimized (const reco::Muon &muon, const double ptThreshold=200., const double tune1=17., const double tune2=40., const double dptcut=0.25)
 
reco::Muon::MuonTrackTypePair TMR (const reco::TrackRef &trackerTrack, const reco::TrackRef &fmsTrack, const double tune=4.)
 
double trackProbability (const reco::TrackRef track)
 

Variables

static const SelectionTypeStringToEnum selectionTypeStringToEnumMap []
 
static const SelectorStringToEnum selectorStringToEnumMap []
 

Enumeration Type Documentation

◆ AlgorithmType

enum muon::AlgorithmType
Enumerator
TMLastStation 
TM2DCompatibility 
TMOneStation 
RPCMu 
ME0Mu 
GEMMu 

Definition at line 156 of file MuonSelectors.h.

156 { TMLastStation, TM2DCompatibility, TMOneStation, RPCMu, ME0Mu, GEMMu };

◆ SelectionType

enum muon::SelectionType

Selector type.

Enumerator
All 
AllGlobalMuons 
AllStandAloneMuons 
AllTrackerMuons 
TrackerMuonArbitrated 
AllArbitrated 
GlobalMuonPromptTight 
TMLastStationLoose 
TMLastStationTight 
TM2DCompatibilityLoose 
TM2DCompatibilityTight 
TMOneStationLoose 
TMOneStationTight 
TMLastStationOptimizedLowPtLoose 
TMLastStationOptimizedLowPtTight 
GMTkChiCompatibility 
GMStaChiCompatibility 
GMTkKinkTight 
TMLastStationAngLoose 
TMLastStationAngTight 
TMOneStationAngLoose 
TMOneStationAngTight 
TMLastStationOptimizedBarrelLowPtLoose 
TMLastStationOptimizedBarrelLowPtTight 
RPCMuLoose 
AllME0Muons 
ME0MuonArbitrated 
AllGEMMuons 
GEMMuonArbitrated 
TriggerIdLoose 

Definition at line 18 of file MuonSelectors.h.

18  {
19  All = 0, // dummy options - always true
20  AllGlobalMuons = 1, // checks isGlobalMuon flag
21  AllStandAloneMuons = 2, // checks isStandAloneMuon flag
22  AllTrackerMuons = 3, // checks isTrackerMuon flag
23  TrackerMuonArbitrated = 4, // resolve ambiguity of sharing segments
24  AllArbitrated = 5, // all muons with the tracker muon arbitrated
25  GlobalMuonPromptTight = 6, // global muons with tighter fit requirements
26  TMLastStationLoose = 7, // penetration depth loose selector
27  TMLastStationTight = 8, // penetration depth tight selector
28  TM2DCompatibilityLoose = 9, // likelihood based loose selector
29  TM2DCompatibilityTight = 10, // likelihood based tight selector
30  TMOneStationLoose = 11, // require one well matched segment
31  TMOneStationTight = 12, // require one well matched segment
32  TMLastStationOptimizedLowPtLoose = 13, // combination of TMLastStation and TMOneStation
33  TMLastStationOptimizedLowPtTight = 14, // combination of TMLastStation and TMOneStation
34  GMTkChiCompatibility = 15, // require tk stub have good chi2 relative to glb track
35  GMStaChiCompatibility = 16, // require sta stub have good chi2 compatibility relative to glb track
36  GMTkKinkTight = 17, // require a small kink value in the tracker stub
37  TMLastStationAngLoose = 18, // TMLastStationLoose with additional angular cuts
38  TMLastStationAngTight = 19, // TMLastStationTight with additional angular cuts
39  TMOneStationAngLoose = 20, // TMOneStationLoose with additional angular cuts
40  TMOneStationAngTight = 21, // TMOneStationTight with additional angular cuts
41  // The two algorithms that follow are identical to what were known as
42  // TMLastStationOptimizedLowPt* (sans the Barrel) as late as revision
43  // 1.7 of this file. The names were changed because indeed the low pt
44  // optimization applies only to the barrel region, whereas the sel-
45  // ectors above are more efficient at low pt in the endcaps, which is
46  // what we feel is more suggestive of the algorithm name. This will be
47  // less confusing for future generations of CMS members, I hope...
48  // combination of TMLastStation and TMOneStation but with low pT optimization in barrel only
49  TMLastStationOptimizedBarrelLowPtLoose = 22,
50  // combination of TMLastStation and TMOneStation but with low pT optimization in barrel only
51  TMLastStationOptimizedBarrelLowPtTight = 23,
52  RPCMuLoose = 24, // checks isRPCMuon flag (require two well matched hits in different RPC layers)
53  AllME0Muons = 25,
54  ME0MuonArbitrated = 26,
55  AllGEMMuons = 27,
56  GEMMuonArbitrated = 28,
57  TriggerIdLoose = 29
58  };

Function Documentation

◆ caloCompatibility()

float muon::caloCompatibility ( const reco::Muon &  muon)

Definition at line 58 of file MuonSelectors.cc.

Referenced by CSCEfficiency::filter(), isGoodMuon(), and pat::MuonSelector::muIdSelection_().

58 { return muon.caloCompatibility(); }

◆ getTevRefitTrack()

reco::TrackRef muon::getTevRefitTrack ( const reco::TrackRef &  combinedTrack,
const reco::TrackToTrackMap &  map 
)

Definition at line 131 of file MuonCocktails.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and genParticles_cff::map.

Referenced by MuonIdProducer::makeMuon().

131  {
132  reco::TrackToTrackMap::const_iterator it = map.find(combinedTrack);
133  return it == map.end() ? reco::TrackRef() : it->val;
134 }
edm::AssociationMap::const_iterator
const iterator
Definition: AssociationMap.h:78
genParticles_cff.map
map
Definition: genParticles_cff.py:11
reco::TrackRef
edm::Ref< TrackCollection > TrackRef
persistent reference to a Track
Definition: TrackFwd.h:20
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48
heppy_batch.val
val
Definition: heppy_batch.py:351

◆ isGoodMuon() [1/3]

bool muon::isGoodMuon ( const reco::Muon &  muon,
SelectionType  type,
reco::Muon::ArbitrationType  arbitrationType = reco::Muon::SegmentAndTrackArbitration 
)

main GoodMuon wrapper call

Definition at line 646 of file MuonSelectors.cc.

References funct::abs(), All, AllArbitrated, AllGEMMuons, AllGlobalMuons, AllME0Muons, AllStandAloneMuons, AllTrackerMuons, arbitrationType, GEMMu, GEMMuonArbitrated, GlobalMuonPromptTight, GMStaChiCompatibility, GMTkChiCompatibility, GMTkKinkTight, isLooseTriggerMuon(), ME0Mu, ME0MuonArbitrated, RPCMu, RPCMuLoose, TM2DCompatibility, TM2DCompatibilityLoose, TM2DCompatibilityTight, TMLastStation, TMLastStationAngLoose, TMLastStationAngTight, TMLastStationLoose, TMLastStationOptimizedBarrelLowPtLoose, TMLastStationOptimizedBarrelLowPtTight, TMLastStationOptimizedLowPtLoose, TMLastStationOptimizedLowPtTight, TMLastStationTight, TMOneStation, TMOneStationAngLoose, TMOneStationAngTight, TMOneStationLoose, TMOneStationTight, TrackerMuonArbitrated, and TriggerIdLoose.

Referenced by BPHSoftMuonSelect::accept(), ExampleMuonAnalyzer::analyze(), BPhysicsOniaDQM::analyze(), TrackEfficiencyMonitor::analyze(), HiggsDQM::analyze(), MuonIdVal::analyze(), MuonCosmicCompatibilityFiller::checkMuonID(), PFRecoTauDiscriminationAgainstMuon::discriminate(), ZtoMMEventSelector::filter(), ttbarEventSelector::filter(), HLTMuonTrkL1TFilter::hltFilter(), HLTMuonTrkFilter::hltFilter(), HLTDiMuonGlbTrkFilter::hltFilter(), HLTMuonL3PreFilter::hltFilter(), PFMuonAlgo::isGlobalLooseMuon(), PFMuonAlgo::isGlobalTightMuon(), isLooseTriggerMuon(), isSoftMuon(), isTightMuon(), PFMuonAlgo::isTightMuonPOG(), PFMuonAlgo::isTrackerLooseMuon(), PFMuonAlgo::isTrackerTightMuon(), PFB::match(), JetPlusTrackCorrector::matchMuons(), pat::MuonSelector::muIdSelection_(), MuonSelectorVIDWrapper< selectionType, arbitrationType >::operator()(), MuonBadTrackFilter::printMuonProperties(), PFMuonAlgo::printMuonProperties(), MuonIDFilterProducerForHLT::produce(), MuonSelectionTypeValueMapProducer::produce(), MuonRefProducer::produce(), MuonTrackProducer::produce(), ZtoMMMuonTrackProducer::produce(), pat::PackedCandidateMuonSelectorProducer::produce(), SoftLepton::produce(), MuonSpecialVariables::setArbitratedTracker(), and MuonSelectorVIDWrapper< selectionType, arbitrationType >::value().

646  {
647  switch (type) {
648  case muon::All:
649  return true;
650  break;
651  case muon::AllGlobalMuons:
652  return muon.isGlobalMuon();
653  break;
654  case muon::AllTrackerMuons:
655  return muon.isTrackerMuon();
656  break;
657  case muon::AllStandAloneMuons:
658  return muon.isStandAloneMuon();
659  break;
660  case muon::TrackerMuonArbitrated:
661  return muon.isTrackerMuon() && muon.numberOfMatches(arbitrationType) > 0;
662  break;
663  case muon::AllArbitrated:
664  return !muon.isTrackerMuon() || muon.numberOfMatches(arbitrationType) > 0;
665  break;
666  case muon::GlobalMuonPromptTight:
667  return muon.isGlobalMuon() && muon.globalTrack()->normalizedChi2() < 10. &&
668  muon.globalTrack()->hitPattern().numberOfValidMuonHits() > 0;
669  break;
670  // For "Loose" algorithms we choose maximum y quantity cuts of 1E9 instead of
671  // 9999 as before. We do this because the muon methods return 999999 (note
672  // there are six 9's) when the requested information is not available. For
673  // example, if a muon fails to traverse the z measuring superlayer in a station
674  // in the DT, then all methods involving segmentY in this station return
675  // 999999 to demonstrate that the information is missing. In order to not
676  // penalize muons for missing y information in Loose algorithms where we do
677  // not care at all about y information, we raise these limits. In the
678  // TMLastStation and TMOneStation algorithms we actually use this huge number
679  // to determine whether to consider y information at all.
680  case muon::TMLastStationLoose:
681  return muon.isTrackerMuon() &&
682  isGoodMuon(muon, TMLastStation, 2, 3, 3, 1E9, 1E9, -3, -3, arbitrationType, true, false);
683  break;
684  case muon::TMLastStationTight:
685  return muon.isTrackerMuon() &&
686  isGoodMuon(muon, TMLastStation, 2, 3, 3, 3, 3, -3, -3, arbitrationType, true, false);
687  break;
688  case muon::TMOneStationLoose:
689  return muon.isTrackerMuon() &&
690  isGoodMuon(muon, TMOneStation, 1, 3, 3, 1E9, 1E9, 1E9, 1E9, arbitrationType, false, false);
691  break;
692  case muon::TMOneStationTight:
693  return muon.isTrackerMuon() &&
694  isGoodMuon(muon, TMOneStation, 1, 3, 3, 3, 3, 1E9, 1E9, arbitrationType, false, false);
695  break;
696  case muon::TMLastStationOptimizedLowPtLoose:
697  if (muon.pt() < 8. && std::abs(muon.eta()) < 1.2)
698  return muon.isTrackerMuon() &&
699  isGoodMuon(muon, TMOneStation, 1, 3, 3, 1E9, 1E9, 1E9, 1E9, arbitrationType, false, false);
700  else
701  return muon.isTrackerMuon() &&
702  isGoodMuon(muon, TMLastStation, 2, 3, 3, 1E9, 1E9, -3, -3, arbitrationType, false, false);
703  break;
704  case muon::TMLastStationOptimizedLowPtTight:
705  if (muon.pt() < 8. && std::abs(muon.eta()) < 1.2)
706  return muon.isTrackerMuon() &&
707  isGoodMuon(muon, TMOneStation, 1, 3, 3, 3, 3, 1E9, 1E9, arbitrationType, false, false);
708  else
709  return muon.isTrackerMuon() &&
710  isGoodMuon(muon, TMLastStation, 2, 3, 3, 3, 3, -3, -3, arbitrationType, false, false);
711  break;
712  //compatibility loose
713  case muon::TM2DCompatibilityLoose:
714  return muon.isTrackerMuon() && isGoodMuon(muon, TM2DCompatibility, 0.7, arbitrationType);
715  break;
716  //compatibility tight
717  case muon::TM2DCompatibilityTight:
718  return muon.isTrackerMuon() && isGoodMuon(muon, TM2DCompatibility, 1.0, arbitrationType);
719  break;
720  case muon::GMTkChiCompatibility:
721  return muon.isGlobalMuon() && muon.isQualityValid() &&
722  std::abs(muon.combinedQuality().trkRelChi2 - muon.innerTrack()->normalizedChi2()) < 2.0;
723  break;
724  case muon::GMStaChiCompatibility:
725  return muon.isGlobalMuon() && muon.isQualityValid() &&
726  std::abs(muon.combinedQuality().staRelChi2 - muon.outerTrack()->normalizedChi2()) < 2.0;
727  break;
728  case muon::GMTkKinkTight:
729  return muon.isGlobalMuon() && muon.isQualityValid() && muon.combinedQuality().trkKink < 100.0;
730  break;
731  case muon::TMLastStationAngLoose:
732  return muon.isTrackerMuon() &&
733  isGoodMuon(muon, TMLastStation, 2, 3, 3, 1E9, 1E9, -3, -3, arbitrationType, false, true);
734  break;
735  case muon::TMLastStationAngTight:
736  return muon.isTrackerMuon() &&
737  isGoodMuon(muon, TMLastStation, 2, 3, 3, 3, 3, -3, -3, arbitrationType, false, true);
738  break;
739  case muon::TMOneStationAngLoose:
740  return muon.isTrackerMuon() &&
741  isGoodMuon(muon, TMOneStation, 1, 3, 3, 1E9, 1E9, 1E9, 1E9, arbitrationType, false, true);
742  break;
743  case muon::TMOneStationAngTight:
744  return muon.isTrackerMuon() &&
745  isGoodMuon(muon, TMOneStation, 1, 3, 3, 3, 3, 1E9, 1E9, arbitrationType, false, true);
746  break;
747  case muon::TMLastStationOptimizedBarrelLowPtLoose:
748  if (muon.pt() < 8. && std::abs(muon.eta()) < 1.2)
749  return muon.isTrackerMuon() &&
750  isGoodMuon(muon, TMOneStation, 1, 3, 3, 1E9, 1E9, 1E9, 1E9, arbitrationType, false, false);
751  else
752  return muon.isTrackerMuon() &&
753  isGoodMuon(muon, TMLastStation, 2, 3, 3, 1E9, 1E9, -3, -3, arbitrationType, true, false);
754  break;
755  case muon::TMLastStationOptimizedBarrelLowPtTight:
756  if (muon.pt() < 8. && std::abs(muon.eta()) < 1.2)
757  return muon.isTrackerMuon() &&
758  isGoodMuon(muon, TMOneStation, 1, 3, 3, 3, 3, 1E9, 1E9, arbitrationType, false, false);
759  else
760  return muon.isTrackerMuon() &&
761  isGoodMuon(muon, TMLastStation, 2, 3, 3, 3, 3, -3, -3, arbitrationType, true, false);
762  break;
763  case muon::RPCMuLoose:
764  return muon.isRPCMuon() && isGoodMuon(muon, RPCMu, 2, 20, 4, 1e9, 1e9, 1e9, 1e9, arbitrationType, false, false);
765  break;
766  case muon::AllME0Muons:
767  return muon.isME0Muon();
768  break;
769  case muon::ME0MuonArbitrated:
770  return muon.isME0Muon() &&
771  isGoodMuon(muon, ME0Mu, 1, 1e9, 1e9, 1e9, 1e9, 1e9, 1e9, arbitrationType, false, false);
772  break;
773  case muon::AllGEMMuons:
774  return muon.isGEMMuon();
775  break;
776  case muon::GEMMuonArbitrated:
777  return muon.isGEMMuon() &&
778  isGoodMuon(muon, GEMMu, 1, 1e9, 1e9, 1e9, 1e9, 1e9, 1e9, arbitrationType, false, false);
779  break;
780  case muon::TriggerIdLoose:
781  return isLooseTriggerMuon(muon);
782  break;
783  default:
784  return false;
785  }
786 }
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
muon::isGoodMuon
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
Definition: MuonSelectors.cc:646
arbitrationType
const reco::Muon::ArbitrationType arbitrationType
Definition: MuonMvaIDEstimator.cc:31
muon::isLooseTriggerMuon
bool isLooseTriggerMuon(const reco::Muon &)
Definition: MuonSelectors.cc:863

◆ isGoodMuon() [2/3]

bool muon::isGoodMuon ( const reco::Muon &  muon,
AlgorithmType  type,
double  minCompatibility,
reco::Muon::ArbitrationType  arbitrationType 
)

Definition at line 290 of file MuonSelectors.cc.

References arbitrationType, caloCompatibility(), segmentCompatibility(), and TM2DCompatibility.

293  {
294  if (!muon.isMatchesValid())
295  return false;
296  bool goodMuon = false;
297 
298  switch (type) {
299  case TM2DCompatibility:
300  // Simplistic first cut in the 2D segment- vs calo-compatibility plane. Will have to be refined!
301  if (((0.8 * caloCompatibility(muon)) + (1.2 * segmentCompatibility(muon, arbitrationType))) > minCompatibility)
302  goodMuon = true;
303  else
304  goodMuon = false;
305  return goodMuon;
306  break;
307  default:
308  // LogTrace("MuonIdentification")<<" // Invalid Algorithm Type called!";
309  goodMuon = false;
310  return goodMuon;
311  break;
312  }
313 }
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
muon::caloCompatibility
float caloCompatibility(const reco::Muon &muon)
Definition: MuonSelectors.cc:58
muon::segmentCompatibility
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
Definition: MuonSelectors.cc:61
arbitrationType
const reco::Muon::ArbitrationType arbitrationType
Definition: MuonMvaIDEstimator.cc:31

◆ isGoodMuon() [3/3]

bool muon::isGoodMuon ( const reco::Muon &  muon,
AlgorithmType  type,
int  minNumberOfMatches,
double  maxAbsDx,
double  maxAbsPullX,
double  maxAbsDy,
double  maxAbsPullY,
double  maxChamberDist,
double  maxChamberDistPull,
reco::Muon::ArbitrationType  arbitrationType,
bool  syncMinNMatchesNRequiredStationsInBarrelOnly = true,
bool  applyAlsoAngularCuts = false 
)

Definition at line 315 of file MuonSelectors.cc.

References funct::abs(), arbitrationType, hgcalTestNeighbor_cfi::detector, dX(), dY(), fftjetvertexadder_cfi::errX, MuonSubdetId::GEM, GEMMu, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, HLT_2024v12_cff::maxAbsDx, HLT_2024v12_cff::maxAbsDy, HLT_2024v12_cff::maxAbsPullX, HLT_2024v12_cff::maxAbsPullY, MuonSubdetId::ME0, ME0Mu, HLT_2024v12_cff::minNumberOfMatches, or, RequiredStationMask(), MuonSubdetId::RPC, RPCMu, mathSSE::sqrt(), relativeConstraints::station, TMLastStation, and TMOneStation.

326  {
327  if (!muon.isMatchesValid())
328  return false;
329  bool goodMuon = false;
330 
331  if (type == TMLastStation) {
332  // To satisfy my own paranoia, if the user specifies that the
333  // minimum number of matches is zero, then return true.
334  if (minNumberOfMatches == 0)
335  return true;
336 
337  unsigned int theStationMask = muon.stationMask(arbitrationType);
338  unsigned int theRequiredStationMask =
339  RequiredStationMask(muon, maxChamberDist, maxChamberDistPull, arbitrationType);
340 
341  // Require that there be at least a minimum number of segments
342  int numSegs = 0;
343  int numRequiredStations = 0;
344  for (int it = 0; it < 8; ++it) {
345  if (theStationMask & 1 << it)
346  ++numSegs;
347  if (theRequiredStationMask & 1 << it)
348  ++numRequiredStations;
349  }
350 
351  // Make sure the minimum number of matches is not greater than
352  // the number of required stations but still greater than zero
353  if (syncMinNMatchesNRequiredStationsInBarrelOnly) {
354  // Note that we only do this in the barrel region!
355  if (std::abs(muon.eta()) < 1.2) {
356  if (minNumberOfMatches > numRequiredStations)
357  minNumberOfMatches = numRequiredStations;
358  if (minNumberOfMatches < 1) //SK: this only happens for negative values
359  minNumberOfMatches = 1;
360  }
361  } else {
362  if (minNumberOfMatches > numRequiredStations)
363  minNumberOfMatches = numRequiredStations;
364  if (minNumberOfMatches < 1) //SK: this only happens for negative values
365  minNumberOfMatches = 1;
366  }
367 
368  if (numSegs >= minNumberOfMatches)
369  goodMuon = true;
370 
371  // Require that last required station have segment
372  // If there are zero required stations keep track
373  // of the last station with a segment so that we may
374  // apply the quality cuts below to it instead
375  int lastSegBit = 0;
376  if (theRequiredStationMask) {
377  for (int stationIdx = 7; stationIdx >= 0; --stationIdx)
378  if (theRequiredStationMask & 1 << stationIdx) {
379  if (theStationMask & 1 << stationIdx) {
380  lastSegBit = stationIdx;
381  goodMuon &= 1;
382  break;
383  } else {
384  goodMuon = false;
385  break;
386  }
387  }
388  } else {
389  for (int stationIdx = 7; stationIdx >= 0; --stationIdx)
390  if (theStationMask & 1 << stationIdx) {
391  lastSegBit = stationIdx;
392  break;
393  }
394  }
395 
396  if (!goodMuon)
397  return false;
398 
399  // Impose pull cuts on last segment
400  int station = 0, detector = 0;
401  station = lastSegBit < 4 ? lastSegBit + 1 : lastSegBit - 3;
402  detector = lastSegBit < 4 ? 1 : 2;
403 
404  // Check x information
405  if (std::abs(muon.pullX(station, detector, arbitrationType, true)) > maxAbsPullX &&
406  std::abs(muon.dX(station, detector, arbitrationType)) > maxAbsDx)
407  return false;
408 
409  if (applyAlsoAngularCuts && std::abs(muon.pullDxDz(station, detector, arbitrationType, true)) > maxAbsPullX)
410  return false;
411 
412  // Is this a tight algorithm, i.e. do we bother to check y information?
413  if (maxAbsDy < 999999) { // really if maxAbsDy < 1E9 as currently defined
414 
415  // Check y information
416  if (detector == 2) { // CSC
417  if (std::abs(muon.pullY(station, 2, arbitrationType, true)) > maxAbsPullY &&
418  std::abs(muon.dY(station, 2, arbitrationType)) > maxAbsDy)
419  return false;
420 
421  if (applyAlsoAngularCuts && std::abs(muon.pullDyDz(station, 2, arbitrationType, true)) > maxAbsPullY)
422  return false;
423  } else {
424  //
425  // In DT, if this is a "Tight" algorithm and the last segment is
426  // missing y information (always the case in station 4!!!), impose
427  // respective cuts on the next station in the stationMask that has
428  // a segment with y information. If there are no segments with y
429  // information then there is nothing to penalize. Should we
430  // penalize in Tight for having zero segments with y information?
431  // That is the fundamental question. Of course I am being uber
432  // paranoid; if this is a good muon then there will probably be at
433  // least one segment with y information but not always. Suppose
434  // somehow a muon only creates segments in station 4, then we
435  // definitely do not want to require that there be at least one
436  // segment with y information because we will lose it completely.
437  //
438 
439  for (int stationIdx = station; stationIdx > 0; --stationIdx) {
440  if (!(theStationMask & 1 << (stationIdx - 1))) // don't bother if the station is not in the stationMask
441  continue;
442 
443  if (muon.dY(stationIdx, 1, arbitrationType) > 999998) // no y-information
444  continue;
445 
446  if (std::abs(muon.pullY(stationIdx, 1, arbitrationType, true)) > maxAbsPullY &&
447  std::abs(muon.dY(stationIdx, 1, arbitrationType)) > maxAbsDy) {
448  return false;
449  }
450 
451  if (applyAlsoAngularCuts && std::abs(muon.pullDyDz(stationIdx, 1, arbitrationType, true)) > maxAbsPullY)
452  return false;
453 
454  // If we get this far then great this is a good muon
455  return true;
456  }
457  }
458  }
459 
460  return goodMuon;
461  } // TMLastStation
462 
463  // TMOneStation requires only that there be one "good" segment, regardless
464  // of the required stations. We do not penalize if there are absolutely zero
465  // segments with y information in the Tight algorithm. Maybe I'm being
466  // paranoid but so be it. If it's really a good muon then we will probably
467  // find at least one segment with both x and y information but you never
468  // know, and I don't want to deal with a potential inefficiency in the DT
469  // like we did with the original TMLastStation. Incidentally, not penalizing
470  // for total lack of y information in the Tight algorithm is what is done in
471  // the new TMLastStation
472  //
473  if (type == TMOneStation) {
474  unsigned int theStationMask = muon.stationMask(arbitrationType);
475 
476  // Of course there must be at least one segment
477  if (!theStationMask)
478  return false;
479 
480  int station = 0, detector = 0;
481  // Keep track of whether or not there is a DT segment with y information.
482  // In the end, if it turns out there are absolutely zero DT segments with
483  // y information, require only that there was a segment with good x info.
484  // This of course only applies to the Tight algorithms.
485  bool existsGoodDTSegX = false;
486  bool existsDTSegY = false;
487 
488  // Impose cuts on the segments in the station mask until we find a good one
489  // Might as well start with the lowest bit to speed things up.
490  for (int stationIdx = 0; stationIdx <= 7; ++stationIdx)
491  if (theStationMask & 1 << stationIdx) {
492  station = stationIdx < 4 ? stationIdx + 1 : stationIdx - 3;
493  detector = stationIdx < 4 ? 1 : 2;
494 
495  if ((std::abs(muon.pullX(station, detector, arbitrationType, true)) > maxAbsPullX &&
496  std::abs(muon.dX(station, detector, arbitrationType)) > maxAbsDx) ||
497  (applyAlsoAngularCuts && std::abs(muon.pullDxDz(station, detector, arbitrationType, true)) > maxAbsPullX))
498  continue;
499  else if (detector == 1)
500  existsGoodDTSegX = true;
501 
502  // Is this a tight algorithm? If yes, use y information
503  if (maxAbsDy < 999999) {
504  if (detector == 2) { // CSC
505  if ((std::abs(muon.pullY(station, 2, arbitrationType, true)) > maxAbsPullY &&
506  std::abs(muon.dY(station, 2, arbitrationType)) > maxAbsDy) ||
507  (applyAlsoAngularCuts && std::abs(muon.pullDyDz(station, 2, arbitrationType, true)) > maxAbsPullY))
508  continue;
509  } else {
510  if (muon.dY(station, 1, arbitrationType) > 999998) // no y-information
511  continue;
512  else
513  existsDTSegY = true;
514 
515  if ((std::abs(muon.pullY(station, 1, arbitrationType, true)) > maxAbsPullY &&
516  std::abs(muon.dY(station, 1, arbitrationType)) > maxAbsDy) ||
517  (applyAlsoAngularCuts && std::abs(muon.pullDyDz(station, 1, arbitrationType, true)) > maxAbsPullY)) {
518  continue;
519  }
520  }
521  }
522 
523  // If we get this far then great this is a good muon
524  return true;
525  }
526 
527  // If we get this far then for sure there are no "good" CSC segments. For
528  // DT, check if there were any segments with y information. If there
529  // were none, but there was a segment with good x, then we're happy. If
530  // there WERE segments with y information, then they must have been shit
531  // since we are here so fail it. Of course, if this is a Loose algorithm
532  // then fail immediately since if we had good x we would already have
533  // returned true
534  if (maxAbsDy < 999999) {
535  if (existsDTSegY)
536  return false;
537  else if (existsGoodDTSegX)
538  return true;
539  } else
540  return false;
541  } // TMOneStation
542 
543  if (type == RPCMu) {
544  if (minNumberOfMatches == 0)
545  return true;
546 
547  int nMatch = 0;
548  for (const auto& chamberMatch : muon.matches()) {
549  if (chamberMatch.detector() != MuonSubdetId::RPC)
550  continue;
551 
552  const float trkX = chamberMatch.x;
553  const float errX = chamberMatch.xErr;
554 
555  for (const auto& rpcMatch : chamberMatch.rpcMatches) {
556  const float rpcX = rpcMatch.x;
557  const float dX = std::abs(rpcX - trkX);
558  if (dX < maxAbsDx or dX < maxAbsPullX * errX) {
559  ++nMatch;
560  break;
561  }
562  }
563  }
564 
565  if (nMatch >= minNumberOfMatches)
566  return true;
567  else
568  return false;
569  } // RPCMu
570 
571  if (type == ME0Mu) {
572  if (minNumberOfMatches == 0)
573  return true;
574 
575  int nMatch = 0;
576  for (const auto& chamberMatch : muon.matches()) {
577  if (chamberMatch.detector() != MuonSubdetId::ME0)
578  continue;
579 
580  const float trkX = chamberMatch.x;
581  const float errX2 = chamberMatch.xErr * chamberMatch.xErr;
582  const float trkY = chamberMatch.y;
583  const float errY2 = chamberMatch.yErr * chamberMatch.yErr;
584 
585  for (const auto& segment : chamberMatch.me0Matches) {
586  const float me0X = segment.x;
587  const float me0ErrX2 = segment.xErr * segment.xErr;
588  const float me0Y = segment.y;
589  const float me0ErrY2 = segment.yErr * segment.yErr;
590 
591  const float dX = std::abs(me0X - trkX);
592  const float dY = std::abs(me0Y - trkY);
593  const float invPullX2 = errX2 + me0ErrX2;
594  const float invPullY2 = errY2 + me0ErrY2;
595 
596  if ((dX < maxAbsDx or dX < maxAbsPullX * std::sqrt(invPullX2)) and
597  (dY < maxAbsDy or dY < maxAbsPullY * std::sqrt(invPullY2))) {
598  ++nMatch;
599  break;
600  }
601  }
602  }
603 
604  return (nMatch >= minNumberOfMatches);
605  } // ME0Mu
606 
607  if (type == GEMMu) {
608  if (minNumberOfMatches == 0)
609  return true;
610 
611  int nMatch = 0;
612  for (const auto& chamberMatch : muon.matches()) {
613  if (chamberMatch.detector() != MuonSubdetId::GEM)
614  continue;
615 
616  const float trkX = chamberMatch.x;
617  const float errX2 = chamberMatch.xErr * chamberMatch.xErr;
618  const float trkY = chamberMatch.y;
619  const float errY2 = chamberMatch.yErr * chamberMatch.yErr;
620 
621  for (const auto& segment : chamberMatch.gemMatches) {
622  const float gemX = segment.x;
623  const float gemErrX2 = segment.xErr * segment.xErr;
624  const float gemY = segment.y;
625  const float gemErrY2 = segment.yErr * segment.yErr;
626 
627  const float dX = std::abs(gemX - trkX);
628  const float dY = std::abs(gemY - trkY);
629  const float invPullX2 = errX2 + gemErrX2;
630  const float invPullY2 = errY2 + gemErrY2;
631 
632  if ((dX < maxAbsDx or dX < maxAbsPullX * std::sqrt(invPullX2)) and
633  (dY < maxAbsDy or dY < maxAbsPullY * std::sqrt(invPullY2))) {
634  ++nMatch;
635  break;
636  }
637  }
638  }
639 
640  return (nMatch >= minNumberOfMatches);
641  } // GEMMu
642 
643  return goodMuon;
644 }
MuonSubdetId::GEM
static constexpr int GEM
Definition: MuonSubdetId.h:14
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
fftjetvertexadder_cfi.errX
errX
Definition: fftjetvertexadder_cfi.py:37
HLT_2024v12_cff.maxAbsDy
maxAbsDy
Definition: HLT_2024v12_cff.py:9508
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
MuonSubdetId::ME0
static constexpr int ME0
Definition: MuonSubdetId.h:15
HLT_2024v12_cff.maxAbsDx
maxAbsDx
Definition: HLT_2024v12_cff.py:9511
dX
float dX(const MatchPair &match)
Definition: SoftMuonMvaRun3Estimator.cc:23
arbitrationType
const reco::Muon::ArbitrationType arbitrationType
Definition: MuonMvaIDEstimator.cc:31
MuonSubdetId::RPC
static constexpr int RPC
Definition: MuonSubdetId.h:13
muon::RequiredStationMask
unsigned int RequiredStationMask(const reco::Muon &muon, double maxChamberDist, double maxChamberDistPull, reco::Muon::ArbitrationType arbitrationType)
Definition: MuonSelectors.cc:40
hgcalTestNeighbor_cfi.detector
detector
Definition: hgcalTestNeighbor_cfi.py:6
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48
HLT_2024v12_cff.maxAbsPullX
maxAbsPullX
Definition: HLT_2024v12_cff.py:9496
HLT_2024v12_cff.minNumberOfMatches
minNumberOfMatches
Definition: HLT_2024v12_cff.py:9514
HLT_2024v12_cff.maxAbsPullY
maxAbsPullY
Definition: HLT_2024v12_cff.py:9507
dY
float dY(const MatchPair &match)
Definition: SoftMuonMvaRun3Estimator.cc:45

◆ isHighPtMuon()

bool muon::isHighPtMuon ( const reco::Muon &  muon,
const reco::Vertex &  vtx 
)

Definition at line 929 of file MuonSelectors.cc.

References funct::abs(), hfClusterShapes_cfi::hits, and L1BJetProducer_cff::vtx.

Referenced by MuonKinVsEtaAnalyzer::analyze(), makeSelectorBitset(), Muon.Muon::muonID(), MuonPOGStandardCut::operator()(), MuonMiniAOD::PassesCut_A(), HLTDQMMuonSelector::passMuonSel(), and MuonPOGStandardCut::value().

929  {
930  if (!muon.isGlobalMuon())
931  return false;
932 
933  bool muValHits = (muon.globalTrack()->hitPattern().numberOfValidMuonHits() > 0 ||
934  muon.tunePMuonBestTrack()->hitPattern().numberOfValidMuonHits() > 0);
935 
936  bool muMatchedSt = muon.numberOfMatchedStations() > 1;
937  if (!muMatchedSt) {
938  if (muon.isTrackerMuon() && muon.numberOfMatchedStations() == 1) {
939  if (muon.expectedNnumberOfMatchedStations() < 2 || !(muon.stationMask() == 1 || muon.stationMask() == 16) ||
940  muon.numberOfMatchedRPCLayers() > 2)
941  muMatchedSt = true;
942  }
943  }
944 
945  bool muID = muValHits && muMatchedSt;
946 
947  bool hits = muon.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
948  muon.innerTrack()->hitPattern().numberOfValidPixelHits() > 0;
949 
950  bool momQuality = muon.tunePMuonBestTrack()->ptError() / muon.tunePMuonBestTrack()->pt() < 0.3;
951 
952  bool ip =
953  std::abs(muon.innerTrack()->dxy(vtx.position())) < 0.2 && std::abs(muon.innerTrack()->dz(vtx.position())) < 0.5;
954 
955  return muID && hits && momQuality && ip;
956 }
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
L1BJetProducer_cff.vtx
vtx
Definition: L1BJetProducer_cff.py:10
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22

◆ isLooseMuon()

bool muon::isLooseMuon ( const reco::Muon &  muon)

Definition at line 891 of file MuonSelectors.cc.

Referenced by DiMuonHistograms::analyze(), EfficiencyAnalyzer::analyze(), MuonKinVsEtaAnalyzer::analyze(), L1TTauOffline::getTightMuons(), isMediumMuon(), Muon.Muon::looseId(), makeSelectorBitset(), Muon.Muon::muonID(), MuonPOGStandardCut::operator()(), MuonMiniAOD::PassesCut_A(), ZCounting::passMuonID(), HLTDQMMuonSelector::passMuonSel(), SoftPFMuonTagInfoProducer::produce(), and MuonPOGStandardCut::value().

891  {
892  return muon.isPFMuon() && (muon.isGlobalMuon() || muon.isTrackerMuon());
893 }

◆ isLooseTriggerMuon()

bool muon::isLooseTriggerMuon ( const reco::Muon &  muon)

Definition at line 863 of file MuonSelectors.cc.

References isGoodMuon(), or, and TMOneStationTight.

Referenced by isGoodMuon(), makeSelectorBitset(), and MuonIDFilterProducerForHLT::produce().

863  {
864  // Requirements:
865  // - no depencence on information not availabe in the muon object
866  // - use only robust inputs
867  bool tk_id = muon::isGoodMuon(muon, TMOneStationTight);
868  if (not tk_id)
869  return false;
870  bool layer_requirements = muon.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5;
871  bool match_requirements =
872  (muon.expectedNnumberOfMatchedStations() < 2) or (muon.numberOfMatchedStations() > 1) or (muon.pt() < 8);
873  return layer_requirements and match_requirements;
874 }
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
muon::isGoodMuon
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
Definition: MuonSelectors.cc:646

◆ isMediumMuon()

bool muon::isMediumMuon ( const reco::Muon &  muon,
bool  run2016_hip_mitigation = false 
)

Definition at line 895 of file MuonSelectors.cc.

References isLooseMuon(), and segmentCompatibility().

Referenced by DiMuonHistograms::analyze(), EfficiencyAnalyzer::analyze(), MuonKinVsEtaAnalyzer::analyze(), HcalHBHEMuonAnalyzer::analyze(), makeSelectorBitset(), MuonPOGStandardCut::operator()(), MuonMiniAOD::PassesCut_A(), ZCounting::passMuonID(), HLTDQMMuonSelector::passMuonSel(), AlCaHcalHBHEMuonProducer::produce(), BadGlobalMuonTagger::tighterId(), and MuonPOGStandardCut::value().

895  {
896  if (not isLooseMuon(muon))
897  return false;
898  if (run2016_hip_mitigation) {
899  if (muon.innerTrack()->validFraction() < 0.49)
900  return false;
901  } else {
902  if (muon.innerTrack()->validFraction() < 0.8)
903  return false;
904  }
905 
906  bool goodGlb = muon.isGlobalMuon() && muon.globalTrack()->normalizedChi2() < 3. &&
907  muon.combinedQuality().chi2LocalPosition < 12. && muon.combinedQuality().trkKink < 20.;
908 
909  return (segmentCompatibility(muon) > (goodGlb ? 0.303 : 0.451));
910 }
muon::isLooseMuon
bool isLooseMuon(const reco::Muon &)
Definition: MuonSelectors.cc:891
muon::segmentCompatibility
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
Definition: MuonSelectors.cc:61

◆ isSoftMuon()

bool muon::isSoftMuon ( const reco::Muon &  muon,
const reco::Vertex &  vtx,
bool  run2016_hip_mitigation = false 
)

Definition at line 912 of file MuonSelectors.cc.

References funct::abs(), reco::TrackBase::highPurity, isGoodMuon(), hgcalTBTopologyTester_cfi::layers, TMOneStationTight, and L1BJetProducer_cff::vtx.

Referenced by DiMuonHistograms::analyze(), MuonKinVsEtaAnalyzer::analyze(), LeptonSkimming::filter(), makeSelectorBitset(), Muon.Muon::muonID(), MuonPOGStandardCut::operator()(), MuonMiniAOD::PassesCut_A(), HLTDQMMuonSelector::passMuonSel(), and MuonPOGStandardCut::value().

912  {
913  bool muID = muon::isGoodMuon(muon, TMOneStationTight);
914 
915  if (!muID)
916  return false;
917 
918  bool layers = muon.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
919  muon.innerTrack()->hitPattern().pixelLayersWithMeasurement() > 0;
920 
921  bool ishighq = muon.innerTrack()->quality(reco::Track::highPurity);
922 
923  bool ip =
924  std::abs(muon.innerTrack()->dxy(vtx.position())) < 0.3 && std::abs(muon.innerTrack()->dz(vtx.position())) < 20.;
925 
926  return layers && ip && (ishighq | run2016_hip_mitigation);
927 }
hgcalTBTopologyTester_cfi.layers
layers
Definition: hgcalTBTopologyTester_cfi.py:8
L1BJetProducer_cff.vtx
vtx
Definition: L1BJetProducer_cff.py:10
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
muon::isGoodMuon
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
Definition: MuonSelectors.cc:646

◆ isTightMuon()

bool muon::isTightMuon ( const reco::Muon &  muon,
const reco::Vertex &  vtx 
)

Definition at line 876 of file MuonSelectors.cc.

References funct::abs(), GlobalMuonPromptTight, hfClusterShapes_cfi::hits, isGoodMuon(), and L1BJetProducer_cff::vtx.

Referenced by METplusTrackMonitor::analyze(), RecoDisplacedMuonValidator::analyze(), RecoMuonValidator::analyze(), MuonKinVsEtaAnalyzer::analyze(), DiMuonHistograms::analyze(), EfficiencyAnalyzer::analyze(), MuonPFAnalyzer::analyze(), MuonRecoOneHLT::analyze(), HcalHBHEMuonAnalyzer::analyze(), L1TMuonDQMOffline::getTightMuons(), makeSelectorBitset(), Muon.Muon::muonID(), MuonPOGStandardCut::operator()(), MuonMiniAOD::PassesCut_A(), ZCounting::passMuonID(), HLTDQMMuonSelector::passMuonSel(), AlCaHcalHBHEMuonProducer::produce(), and MuonPOGStandardCut::value().

876  {
877  if (!muon.isPFMuon() || !muon.isGlobalMuon())
878  return false;
879 
880  bool muID = isGoodMuon(muon, GlobalMuonPromptTight) && (muon.numberOfMatchedStations() > 1);
881 
882  bool hits = muon.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
883  muon.innerTrack()->hitPattern().numberOfValidPixelHits() > 0;
884 
885  bool ip = std::abs(muon.muonBestTrack()->dxy(vtx.position())) < 0.2 &&
886  std::abs(muon.muonBestTrack()->dz(vtx.position())) < 0.5;
887 
888  return muID && hits && ip;
889 }
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
L1BJetProducer_cff.vtx
vtx
Definition: L1BJetProducer_cff.py:10
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
muon::isGoodMuon
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
Definition: MuonSelectors.cc:646

◆ isTrackerHighPtMuon()

bool muon::isTrackerHighPtMuon ( const reco::Muon &  muon,
const reco::Vertex &  vtx 
)

Definition at line 958 of file MuonSelectors.cc.

References funct::abs(), hfClusterShapes_cfi::hits, and L1BJetProducer_cff::vtx.

Referenced by makeSelectorBitset().

958  {
959  bool muID = muon.isTrackerMuon() && muon.track().isNonnull() && (muon.numberOfMatchedStations() > 1);
960  if (!muID)
961  return false;
962 
963  bool hits = muon.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
964  muon.innerTrack()->hitPattern().numberOfValidPixelHits() > 0;
965 
966  bool momQuality = muon.tunePMuonBestTrack()->ptError() < 0.3 * muon.tunePMuonBestTrack()->pt();
967 
968  bool ip =
969  std::abs(muon.innerTrack()->dxy(vtx.position())) < 0.2 && std::abs(muon.innerTrack()->dz(vtx.position())) < 0.5;
970 
971  return muID && hits && momQuality && ip;
972 }
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
L1BJetProducer_cff.vtx
vtx
Definition: L1BJetProducer_cff.py:10
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22

◆ makeSelectorBitset()

reco::Muon::Selector muon::makeSelectorBitset ( reco::Muon const &  muon,
reco::Vertex const *  vertex = nullptr,
bool  run2016_hip_mitigation = false 
)

Definition at line 1029 of file MuonSelectors.cc.

References funct::abs(), reco::Muon::CutBasedIdGlobalHighPt, reco::Muon::CutBasedIdLoose, reco::Muon::CutBasedIdMedium, reco::Muon::CutBasedIdMediumPrompt, reco::Muon::CutBasedIdTight, reco::Muon::CutBasedIdTrkHighPt, reco::Muon::InTimeMuon, isHighPtMuon(), isLooseMuon(), isLooseTriggerMuon(), isMediumMuon(), isSoftMuon(), isTightMuon(), isTrackerHighPtMuon(), SiStripPI::max, outOfTimeMuon(), reco::Muon::PFIsoLoose, reco::Muon::PFIsoMedium, reco::Muon::PFIsoTight, reco::Muon::PFIsoVeryLoose, reco::Muon::PFIsoVeryTight, reco::Muon::PFIsoVeryVeryTight, TrackValidation_cff::selectors, reco::Muon::SoftCutBasedId, reco::Muon::TkIsoLoose, reco::Muon::TkIsoTight, muons_cff::tkRelIso, reco::Muon::TriggerIdLoose, and bphysicsOniaDQM_cfi::vertex.

Referenced by MuonProducer::produce(), pat::PATMuonProducer::produce(), and pat::LeptonUpdater< T >::recomputeMuonBasicSelectors().

1031  {
1032  // https://twiki.cern.ch/twiki/bin/viewauth/CMS/SWGuideMuonIdRun2
1033  unsigned int selectors = muon.selectors();
1034  // Compute Id and Isolation variables
1035  double chIso = muon.pfIsolationR04().sumChargedHadronPt;
1036  double nIso = muon.pfIsolationR04().sumNeutralHadronEt;
1037  double phoIso = muon.pfIsolationR04().sumPhotonEt;
1038  double puIso = muon.pfIsolationR04().sumPUPt;
1039  double dbCorrectedIsolation = chIso + std::max(nIso + phoIso - .5 * puIso, 0.);
1040  double dbCorrectedRelIso = dbCorrectedIsolation / muon.pt();
1041  double tkRelIso = muon.isolationR03().sumPt / muon.pt();
1042 
1043  // Base selectors
1044  if (muon::isLooseMuon(muon))
1045  selectors |= (1UL << reco::Muon::CutBasedIdLoose);
1046  if (vertex) {
1047  if (muon::isTightMuon(muon, *vertex))
1048  selectors |= (1UL << reco::Muon::CutBasedIdTight);
1049  if (muon::isSoftMuon(muon, *vertex, run2016_hip_mitigation))
1050  selectors |= (1UL << reco::Muon::SoftCutBasedId);
1051  if (muon::isHighPtMuon(muon, *vertex))
1052  selectors |= (1UL << reco::Muon::CutBasedIdGlobalHighPt);
1053  if (muon::isTrackerHighPtMuon(muon, *vertex))
1054  selectors |= (1UL << reco::Muon::CutBasedIdTrkHighPt);
1055  }
1056  if (muon::isMediumMuon(muon, run2016_hip_mitigation)) {
1057  selectors |= (1UL << reco::Muon::CutBasedIdMedium);
1058  if (vertex and std::abs(muon.muonBestTrack()->dz(vertex->position())) < 0.1 and
1059  std::abs(muon.muonBestTrack()->dxy(vertex->position())) < 0.02)
1060  selectors |= (1UL << reco::Muon::CutBasedIdMediumPrompt);
1061  }
1062 
1063  // PF isolation
1064  if (dbCorrectedRelIso < 0.40)
1065  selectors |= (1UL << reco::Muon::PFIsoVeryLoose);
1066  if (dbCorrectedRelIso < 0.25)
1067  selectors |= (1UL << reco::Muon::PFIsoLoose);
1068  if (dbCorrectedRelIso < 0.20)
1069  selectors |= (1UL << reco::Muon::PFIsoMedium);
1070  if (dbCorrectedRelIso < 0.15)
1071  selectors |= (1UL << reco::Muon::PFIsoTight);
1072  if (dbCorrectedRelIso < 0.10)
1073  selectors |= (1UL << reco::Muon::PFIsoVeryTight);
1074  if (dbCorrectedRelIso < 0.05)
1075  selectors |= (1UL << reco::Muon::PFIsoVeryVeryTight);
1076 
1077  // Tracker isolation
1078  if (tkRelIso < 0.10)
1079  selectors |= (1UL << reco::Muon::TkIsoLoose);
1080  if (tkRelIso < 0.05)
1081  selectors |= (1UL << reco::Muon::TkIsoTight);
1082 
1083  // Trigger selectors
1084  if (isLooseTriggerMuon(muon))
1085  selectors |= (1UL << reco::Muon::TriggerIdLoose);
1086 
1087  // Timing
1088  if (!outOfTimeMuon(muon))
1089  selectors |= (1UL << reco::Muon::InTimeMuon);
1090 
1091  return static_cast<reco::Muon::Selector>(selectors);
1092 }
bphysicsOniaDQM_cfi.vertex
vertex
Definition: bphysicsOniaDQM_cfi.py:7
muon::isMediumMuon
bool isMediumMuon(const reco::Muon &, bool run2016_hip_mitigation=false)
Definition: MuonSelectors.cc:895
muon::isLooseMuon
bool isLooseMuon(const reco::Muon &)
Definition: MuonSelectors.cc:891
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
muon::isSoftMuon
bool isSoftMuon(const reco::Muon &, const reco::Vertex &, bool run2016_hip_mitigation=false)
Definition: MuonSelectors.cc:912
muon::isHighPtMuon
bool isHighPtMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:929
TrackValidation_cff.selectors
selectors
Definition: TrackValidation_cff.py:271
outOfTimeMuon
bool outOfTimeMuon(const reco::Muon &muon)
Definition: MuonSelectors.cc:1013
muon::isTrackerHighPtMuon
bool isTrackerHighPtMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:958
muon::isTightMuon
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:876
muon::isLooseTriggerMuon
bool isLooseTriggerMuon(const reco::Muon &)
Definition: MuonSelectors.cc:863
reco::Muon::Selector
Selector
Definition: Muon.h:202
muons_cff.tkRelIso
tkRelIso
Definition: muons_cff.py:173

◆ overlap()

bool muon::overlap ( const reco::Muon &  muon1,
const reco::Muon &  muon2,
double  pullX = 1.0,
double  pullY = 1.0,
bool  checkAdjacentChambers = false 
)

Definition at line 788 of file MuonSelectors.cc.

References funct::abs(), MuonSubdetId::CSC, globals_cff::id1, globals_cff::id2, reco::Muon::matches(), reco::Muon::numberOfMatches(), reco::LeafCandidate::pt(), pullX(), pullY(), reco::Muon::SegmentAndTrackArbitration, and mathSSE::sqrt().

789  {
790  unsigned int nMatches1 = muon1.numberOfMatches(reco::Muon::SegmentAndTrackArbitration);
791  unsigned int nMatches2 = muon2.numberOfMatches(reco::Muon::SegmentAndTrackArbitration);
792  unsigned int betterMuon = (muon1.pt() > muon2.pt() ? 1 : 2);
793  for (std::vector<reco::MuonChamberMatch>::const_iterator chamber1 = muon1.matches().begin();
794  chamber1 != muon1.matches().end();
795  ++chamber1)
796  for (std::vector<reco::MuonChamberMatch>::const_iterator chamber2 = muon2.matches().begin();
797  chamber2 != muon2.matches().end();
798  ++chamber2) {
799  // if ( (chamber1->segmentMatches.empty() || chamber2->segmentMatches.empty()) ) continue;
800 
801  // handle case where both muons have information about the same chamber
802  // here we know how close they are
803  if (chamber1->id == chamber2->id) {
804  // found the same chamber
805  if (std::abs(chamber1->x - chamber2->x) <
806  pullX * sqrt(chamber1->xErr * chamber1->xErr + chamber2->xErr * chamber2->xErr)) {
807  if (betterMuon == 1)
808  nMatches2--;
809  else
810  nMatches1--;
811  if (nMatches1 == 0 || nMatches2 == 0)
812  return true;
813  continue;
814  }
815  if (std::abs(chamber1->y - chamber2->y) <
816  pullY * sqrt(chamber1->yErr * chamber1->yErr + chamber2->yErr * chamber2->yErr)) {
817  if (betterMuon == 1)
818  nMatches2--;
819  else
820  nMatches1--;
821  if (nMatches1 == 0 || nMatches2 == 0)
822  return true;
823  }
824  } else {
825  if (!checkAdjacentChambers)
826  continue;
827  // check if tracks are pointing into overlaping region of the CSC detector
828  if (chamber1->id.subdetId() != MuonSubdetId::CSC || chamber2->id.subdetId() != MuonSubdetId::CSC)
829  continue;
830  CSCDetId id1(chamber1->id);
831  CSCDetId id2(chamber2->id);
832  if (id1.endcap() != id2.endcap())
833  continue;
834  if (id1.station() != id2.station())
835  continue;
836  if (id1.ring() != id2.ring())
837  continue;
838  if (std::abs(id1.chamber() - id2.chamber()) > 1)
839  continue;
840  // FIXME: we don't handle 18->1; 36->1 transitions since
841  // I don't know how to check for sure how many chambers
842  // are there. Probably need to hard code some checks.
843 
844  // Now we have to make sure that both tracks are close to an edge
845  // FIXME: ignored Y coordinate for now
846  if (std::abs(chamber1->edgeX) > chamber1->xErr * pullX)
847  continue;
848  if (std::abs(chamber2->edgeX) > chamber2->xErr * pullX)
849  continue;
850  if (chamber1->x * chamber2->x < 0) { // check if the same edge
851  if (betterMuon == 1)
852  nMatches2--;
853  else
854  nMatches1--;
855  if (nMatches1 == 0 || nMatches2 == 0)
856  return true;
857  }
858  }
859  }
860  return false;
861 }
globals_cff.id1
id1
Definition: globals_cff.py:53
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::Muon::numberOfMatches
int numberOfMatches(ArbitrationType type=SegmentAndTrackArbitration) const
get number of chambers with matched segments
reco::Muon::matches
std::vector< MuonChamberMatch > & matches()
get muon matching information
Definition: Muon.h:145
globals_cff.id2
id2
Definition: globals_cff.py:54
pullY
float pullY(const MatchPair &match)
Definition: SoftMuonMvaRun3Estimator.cc:52
MuonSubdetId::CSC
static constexpr int CSC
Definition: MuonSubdetId.h:12
pullX
float pullX(const MatchPair &match)
Definition: SoftMuonMvaRun3Estimator.cc:30

◆ RequiredStationMask()

unsigned int muon::RequiredStationMask ( const reco::Muon &  muon,
double  maxChamberDist,
double  maxChamberDistPull,
reco::Muon::ArbitrationType  arbitrationType 
)

Definition at line 40 of file MuonSelectors.cc.

References arbitrationType.

Referenced by isGoodMuon().

43  {
44  unsigned int theMask = 0;
45 
46  for (int stationIdx = 1; stationIdx < 5; ++stationIdx) {
47  for (int detectorIdx = 1; detectorIdx < 3; ++detectorIdx) {
48  float dist = muon.trackDist(stationIdx, detectorIdx, arbitrationType);
49  if (dist < maxChamberDist &&
50  dist < maxChamberDistPull * muon.trackDistErr(stationIdx, detectorIdx, arbitrationType))
51  theMask += 1 << ((stationIdx - 1) + 4 * (detectorIdx - 1));
52  }
53  }
54  return theMask;
55 }
arbitrationType
const reco::Muon::ArbitrationType arbitrationType
Definition: MuonMvaIDEstimator.cc:31

◆ segmentCompatibility()

float muon::segmentCompatibility ( const reco::Muon &  muon,
reco::Muon::ArbitrationType  arbitrationType = reco::Muon::SegmentAndTrackArbitration 
)

Definition at line 61 of file MuonSelectors.cc.

References arbitrationType, f, mps_fire::i, and funct::pow().

Referenced by ExampleMuonAnalyzer::analyze(), HiggsDQM::analyze(), StudyTriggerHLT::analyze(), MuonIdVal::analyze(), HcalHBHEMuonAnalyzer::analyze(), StudyCaloResponse::analyze(), HcalHBHEMuonHighEtaAnalyzer::analyzeMuon(), MuonCosmicCompatibilityFiller::checkMuonSegments(), PFRecoTauDiscriminationAgainstMuon::discriminate(), ParticleNetFeatureEvaluator::fillParticleFeatures(), BadParticleFilter::filter(), isGoodMuon(), isMediumMuon(), isMediumMuonCustom(), pat::MuonSelector::muIdSelection_(), Muon.Muon::muonID(), HcalIsoTrkAnalyzer::notaMuon(), AlCaHcalIsotrkProducer::notaMuon(), MuonSegmentCompatibilityCut::operator()(), AlCaHcalHBHEMuonProducer::produce(), and MuonSegmentCompatibilityCut::value().

61  {
62  bool use_weight_regain_at_chamber_boundary = true;
63  bool use_match_dist_penalty = true;
64 
65  int nr_of_stations_crossed = 0;
66  std::vector<int> stations_w_track(8);
67  std::vector<int> station_has_segmentmatch(8);
68  std::vector<int> station_was_crossed(8);
69  std::vector<float> stations_w_track_at_boundary(8);
70  std::vector<float> station_weight(8);
71  int position_in_stations = 0;
72  float full_weight = 0.;
73 
74  for (int i = 1; i <= 8; ++i) {
75  // ********************************************************;
76  // *** fill local info for this muon (do some counting) ***;
77  // ************** begin ***********************************;
78  if (i <= 4) { // this is the section for the DTs
79  float thisTrackDist = muon.trackDist(i, 1, arbitrationType);
80  if (thisTrackDist < 999999) { //current "raw" info that a track is close to a chamber
81  ++nr_of_stations_crossed;
82  station_was_crossed[i - 1] = 1;
83  if (thisTrackDist > -10.)
84  stations_w_track_at_boundary[i - 1] = thisTrackDist;
85  else
86  stations_w_track_at_boundary[i - 1] = 0.;
87  }
88  //current "raw" info that a segment is matched to the current track
89  if (muon.segmentX(i, 1, arbitrationType) < 999999) {
90  station_has_segmentmatch[i - 1] = 1;
91  }
92  } else { // this is the section for the CSCs
93  float thisTrackDist = muon.trackDist(i - 4, 2, arbitrationType);
94  if (thisTrackDist < 999999) { //current "raw" info that a track is close to a chamber
95  ++nr_of_stations_crossed;
96  station_was_crossed[i - 1] = 1;
97  if (thisTrackDist > -10.)
98  stations_w_track_at_boundary[i - 1] = thisTrackDist;
99  else
100  stations_w_track_at_boundary[i - 1] = 0.;
101  }
102  //current "raw" info that a segment is matched to the current track
103  if (muon.segmentX(i - 4, 2, arbitrationType) < 999999) {
104  station_has_segmentmatch[i - 1] = 1;
105  }
106  }
107  // rough estimation of chamber border efficiency (should be parametrized better, this is just a quick guess):
108  // TF1 * merf = new TF1("merf","-0.5*(TMath::Erf(x/6.)-1)",-100,100);
109  // use above value to "unpunish" missing segment if close to border, i.e. rather than not adding any weight, add
110  // the one from the function. Only for dist ~> -10 cm, else full punish!.
111 
112  // ********************************************************;
113  // *** fill local info for this muon (do some counting) ***;
114  // ************** end *************************************;
115  }
116 
117  // ********************************************************;
118  // *** calculate weights for each station *****************;
119  // ************** begin ***********************************;
120  // const float slope = 0.5;
121  // const float attenuate_weight_regain = 1.;
122  // if attenuate_weight_regain < 1., additional punishment if track is close to boundary and no segment
123  const float attenuate_weight_regain = 0.5;
124 
125  for (int i = 1; i <= 8; ++i) { // loop over all possible stations
126 
127  // first set all weights if a station has been crossed
128  // later penalize if a station did not have a matching segment
129 
130  //old logic if(station_has_segmentmatch[i-1] > 0 ) { // the track has an associated segment at the current station
131  if (station_was_crossed[i - 1] > 0) { // the track crossed this chamber (or was nearby)
132  // - Apply a weight depending on the "depth" of the muon passage.
133  // - The station_weight is later reduced for stations with badly matched segments.
134  // - Even if there is no segment but the track passes close to a chamber boundary, the
135  // weight is set non zero and can go up to 0.5 of the full weight if the track is quite
136  // far from any station.
137  ++position_in_stations;
138 
139  switch (nr_of_stations_crossed) { // define different weights depending on how many stations were crossed
140  case 1:
141  station_weight[i - 1] = 1.f;
142  break;
143  case 2:
144  if (position_in_stations == 1)
145  station_weight[i - 1] = 0.33f;
146  else
147  station_weight[i - 1] = 0.67f;
148  break;
149  case 3:
150  if (position_in_stations == 1)
151  station_weight[i - 1] = 0.23f;
152  else if (position_in_stations == 2)
153  station_weight[i - 1] = 0.33f;
154  else
155  station_weight[i - 1] = 0.44f;
156  break;
157  case 4:
158  if (position_in_stations == 1)
159  station_weight[i - 1] = 0.10f;
160  else if (position_in_stations == 2)
161  station_weight[i - 1] = 0.20f;
162  else if (position_in_stations == 3)
163  station_weight[i - 1] = 0.30f;
164  else
165  station_weight[i - 1] = 0.40f;
166  break;
167 
168  default:
169  // LogTrace("MuonIdentification")<<" // Message: A muon candidate track has more than 4 stations with matching segments.";
170  // LogTrace("MuonIdentification")<<" // Did not expect this - please let me know: ibloch@fnal.gov";
171  // for all other cases
172  station_weight[i - 1] = 1.f / nr_of_stations_crossed;
173  }
174 
175  if (use_weight_regain_at_chamber_boundary) { // reconstitute some weight if there is no match but the segment is close to a boundary:
176  if (station_has_segmentmatch[i - 1] <= 0 && stations_w_track_at_boundary[i - 1] != 0.) {
177  // if segment is not present but track in inefficient region, do not count as "missing match" but add some reduced weight.
178  // original "match weight" is currently reduced by at least attenuate_weight_regain, variing with an error function down to 0 if the track is
179  // inside the chamber.
180  // remark: the additional scale of 0.5 normalizes Err to run from 0 to 1 in y
181  station_weight[i - 1] = station_weight[i - 1] * attenuate_weight_regain * 0.5f *
182  (std::erf(stations_w_track_at_boundary[i - 1] / 6.f) + 1.f);
183  } else if (station_has_segmentmatch[i - 1] <= 0 &&
184  stations_w_track_at_boundary[i - 1] == 0.f) { // no segment match and track well inside chamber
185  // full penalization
186  station_weight[i - 1] = 0.f;
187  }
188  } else { // always fully penalize tracks with no matching segment, whether the segment is close to the boundary or not.
189  if (station_has_segmentmatch[i - 1] <= 0)
190  station_weight[i - 1] = 0.f;
191  }
192 
193  // if track has matching segment, but the matching is not high quality, penalize
194  if (station_has_segmentmatch[i - 1] > 0) {
195  if (i <= 4) { // we are in the DTs
196  if (muon.dY(i, 1, arbitrationType) < 999999.f &&
197  muon.dX(i, 1, arbitrationType) < 999999.f) { // have both X and Y match
198  const float pullTot2 =
199  std::pow(muon.pullX(i, 1, arbitrationType), 2.) + std::pow(muon.pullY(i, 1, arbitrationType), 2.);
200  if (pullTot2 > 1.f) {
201  const float dxy2 =
202  std::pow(muon.dX(i, 1, arbitrationType), 2.) + std::pow(muon.dY(i, 1, arbitrationType), 2.);
203  // reduce weight
204  if (use_match_dist_penalty) {
205  // only use pull if 3 sigma is not smaller than 3 cm
206  if (dxy2 < 9.f && pullTot2 > 9.f) {
207  if (dxy2 > 1.f)
208  station_weight[i - 1] *= 1.f / std::pow(dxy2, .125);
209  } else {
210  station_weight[i - 1] *= 1.f / std::pow(pullTot2, .125);
211  }
212  }
213  }
214  } else if (muon.dY(i, 1, arbitrationType) >= 999999.f) { // has no match in Y
215  // has a match in X. Pull larger that 1 to avoid increasing the weight (just penalize, don't anti-penalize)
216  if (muon.pullX(i, 1, arbitrationType) > 1.f) {
217  // reduce weight
218  if (use_match_dist_penalty) {
219  // only use pull if 3 sigma is not smaller than 3 cm
220  if (muon.dX(i, 1, arbitrationType) < 3.f && muon.pullX(i, 1, arbitrationType) > 3.f) {
221  if (muon.dX(i, 1, arbitrationType) > 1.f)
222  station_weight[i - 1] *= 1.f / std::pow(muon.dX(i, 1, arbitrationType), .25);
223  } else {
224  station_weight[i - 1] *= 1.f / std::pow(muon.pullX(i, 1, arbitrationType), .25);
225  }
226  }
227  }
228  } else { // has no match in X
229  // has a match in Y. Pull larger that 1 to avoid increasing the weight (just penalize, don't anti-penalize)
230  if (muon.pullY(i, 1, arbitrationType) > 1.f) {
231  // reduce weight
232  if (use_match_dist_penalty) {
233  // only use pull if 3 sigma is not smaller than 3 cm
234  if (muon.dY(i, 1, arbitrationType) < 3. && muon.pullY(i, 1, arbitrationType) > 3.) {
235  if (muon.dY(i, 1, arbitrationType) > 1.f)
236  station_weight[i - 1] *= 1.f / std::pow(muon.dY(i, 1, arbitrationType), .25);
237  } else {
238  station_weight[i - 1] *= 1.f / std::pow(muon.pullY(i, 1, arbitrationType), .25);
239  }
240  }
241  }
242  }
243  } else { // We are in the CSCs
244  const float pullTot2 =
245  std::pow(muon.pullX(i - 4, 2, arbitrationType), 2.) + std::pow(muon.pullY(i - 4, 2, arbitrationType), 2.);
246  if (pullTot2 > 1.f) {
247  // reduce weight
248  if (use_match_dist_penalty) {
249  const float dxy2 =
250  std::pow(muon.dX(i - 4, 2, arbitrationType), 2.) + std::pow(muon.dY(i - 4, 2, arbitrationType), 2.);
251  // only use pull if 3 sigma is not smaller than 3 cm
252  if (dxy2 < 9.f && pullTot2 > 9.f) {
253  if (dxy2 > 1.f)
254  station_weight[i - 1] *= 1.f / std::pow(dxy2, .125);
255  } else {
256  station_weight[i - 1] *= 1.f / std::pow(pullTot2, .125);
257  }
258  }
259  }
260  }
261  }
262 
263  // Thoughts:
264  // - should penalize if the segment has only x OR y info
265  // - should also use the segment direction, as it now works!
266 
267  } else { // track did not pass a chamber in this station - just reset weight
268  station_weight[i - 1] = 0.f;
269  }
270 
271  //increment final weight for muon:
272  full_weight += station_weight[i - 1];
273  }
274 
275  // if we don't expect any matches, we set the compatibility to
276  // 0.5 as the track is as compatible with a muon as it is with
277  // background - we should maybe rather set it to -0.5!
278  if (nr_of_stations_crossed == 0) {
279  // full_weight = attenuate_weight_regain*0.5;
280  full_weight = 0.5f;
281  }
282 
283  // ********************************************************;
284  // *** calculate weights for each station *****************;
285  // ************** end *************************************;
286 
287  return full_weight;
288 }
mps_fire.i
i
Definition: mps_fire.py:429
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
arbitrationType
const reco::Muon::ArbitrationType arbitrationType
Definition: MuonMvaIDEstimator.cc:31
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29

◆ selectionTypeFromString()

SelectionType muon::selectionTypeFromString ( const std::string &  label)

Definition at line 9 of file MuonSelectors.cc.

References Exception, newFWLiteAna::found, mps_fire::i, label, muon::SelectionTypeStringToEnum::label, selectionTypeStringToEnumMap, relativeConstraints::value, and muon::SelectionTypeStringToEnum::value.

Referenced by MuonSelectionTypeValueMapProducer::MuonSelectionTypeValueMapProducer(), MuonTrackProducer::produce(), and pat::PackedCandidateMuonSelectorProducer::produce().

9  {
10  SelectionType value = (SelectionType)-1;
11  bool found = false;
12  for (int i = 0; selectionTypeStringToEnumMap[i].label && (!found); ++i)
13  if (!strcmp(label.c_str(), selectionTypeStringToEnumMap[i].label)) {
14  found = true;
15  value = selectionTypeStringToEnumMap[i].value;
16  }
17 
18  // in case of unrecognized selection type
19  if (!found)
20  throw cms::Exception("MuonSelectorError") << label << " is not a recognized SelectionType";
21  return value;
22  }
mps_fire.i
i
Definition: mps_fire.py:429
muon::SelectionType
SelectionType
Selector type.
Definition: MuonSelectors.h:18
label
char const * label
Definition: PFTauDecayModeTools.cc:11
value
Definition: value.py:1
muon::selectionTypeStringToEnumMap
static const SelectionTypeStringToEnum selectionTypeStringToEnumMap[]
Definition: MuonSelectors.h:66

◆ selectorFromString()

reco::Muon::Selector muon::selectorFromString ( const std::string &  label)

Definition at line 24 of file MuonSelectors.cc.

References Exception, newFWLiteAna::found, mps_fire::i, label, muon::SelectorStringToEnum::label, selectorStringToEnumMap, relativeConstraints::value, and muon::SelectorStringToEnum::value.

24  {
25  reco::Muon::Selector value = (reco::Muon::Selector)-1;
26  bool found = false;
27  for (int i = 0; selectorStringToEnumMap[i].label && (!found); ++i)
28  if (!strcmp(label.c_str(), selectorStringToEnumMap[i].label)) {
29  found = true;
30  value = selectorStringToEnumMap[i].value;
31  }
32 
33  // in case of unrecognized selection type
34  if (!found)
35  throw cms::Exception("MuonSelectorError") << label << " is not a recognized reco::Muon::Selector";
36  return value;
37  }
mps_fire.i
i
Definition: mps_fire.py:429
muon::SelectorStringToEnum::value
reco::Muon::Selector value
Definition: MuonSelectors.h:104
label
char const * label
Definition: PFTauDecayModeTools.cc:11
value
Definition: value.py:1
muon::selectorStringToEnumMap
static const SelectorStringToEnum selectorStringToEnumMap[]
Definition: MuonSelectors.h:107
reco::Muon::Selector
Selector
Definition: Muon.h:202

◆ sharedSegments()

int muon::sharedSegments ( const reco::Muon &  muon1,
const reco::Muon &  muon2,
unsigned int  segmentArbitrationMask = reco::MuonSegmentMatch::BestInChamberByDR 
)

Determine the number of shared segments between two muons. Comparison is done using the segment references in the reco::Muon object.

Definition at line 974 of file MuonSelectors.cc.

References reco::Muon::matches(), amptDefaultParameters_cff::mu, and runTheMatrix::ret.

Referenced by heppy::CMGMuonCleanerBySegmentsAlgo::clean(), BadGlobalMuonTagger::filter(), and modules::MuonCleanerBySegmentsT< T >::produce().

974  {
975  int ret = 0;
976 
977  // Will do with a stupid double loop, since creating and filling a map is probably _more_ inefficient for a single lookup.
978  for (std::vector<reco::MuonChamberMatch>::const_iterator chamberMatch = mu.matches().begin();
979  chamberMatch != mu.matches().end();
980  ++chamberMatch) {
981  if (chamberMatch->segmentMatches.empty())
982  continue;
983  for (std::vector<reco::MuonChamberMatch>::const_iterator chamberMatch2 = mu2.matches().begin();
984  chamberMatch2 != mu2.matches().end();
985  ++chamberMatch2) {
986  if (chamberMatch2->segmentMatches.empty())
987  continue;
988  if (chamberMatch2->id() != chamberMatch->id())
989  continue;
990  for (std::vector<reco::MuonSegmentMatch>::const_iterator segmentMatch = chamberMatch->segmentMatches.begin();
991  segmentMatch != chamberMatch->segmentMatches.end();
992  ++segmentMatch) {
993  if (!segmentMatch->isMask(segmentArbitrationMask))
994  continue;
995  for (std::vector<reco::MuonSegmentMatch>::const_iterator segmentMatch2 = chamberMatch2->segmentMatches.begin();
996  segmentMatch2 != chamberMatch2->segmentMatches.end();
997  ++segmentMatch2) {
998  if (!segmentMatch2->isMask(segmentArbitrationMask))
999  continue;
1000  if ((segmentMatch->cscSegmentRef.isNonnull() &&
1001  segmentMatch->cscSegmentRef == segmentMatch2->cscSegmentRef) ||
1002  (segmentMatch->dtSegmentRef.isNonnull() && segmentMatch->dtSegmentRef == segmentMatch2->dtSegmentRef)) {
1003  ++ret;
1004  } // is the same
1005  } // segment of mu2 in chamber
1006  } // segment of mu1 in chamber
1007  } // chamber of mu2
1008  } // chamber of mu1
1009 
1010  return ret;
1011 }
runTheMatrix.ret
ret
prodAgent to be discontinued
Definition: runTheMatrix.py:759
amptDefaultParameters_cff.mu
mu
Definition: amptDefaultParameters_cff.py:16

◆ sigmaSwitch() [1/2]

reco::Muon::MuonTrackTypePair muon::sigmaSwitch ( const reco::TrackRef &  combinedTrack,
const reco::TrackRef &  trackerTrack,
const double  nSigma = 2.,
const double  ptThreshold = 200. 
)

Definition at line 139 of file MuonCocktails.cc.

References reco::Muon::CombinedTrack, dumpMFGeometry_cfg::delta, reco::Muon::InnerTrack, HLTSiStripMonitoring_cff::nSigma, ewkMuLumiMonitorDQM_cfi::ptThreshold, and DiMuonV_cfg::threshold.

Referenced by MuonIdProducer::makeMuon(), and sigmaSwitch().

142  {
143  // If either the global or tracker-only fits have pT below threshold
144  // (default 200 GeV), return the tracker-only fit.
145  if (combinedTrack->pt() < ptThreshold || trackerTrack->pt() < ptThreshold)
146  return make_pair(trackerTrack, reco::Muon::InnerTrack);
147 
148  // If both are above the pT threshold, compare the difference in
149  // q/p: if less than two sigma of the tracker-only track, switch to
150  // global. Otherwise, use tracker-only.
151  const double delta = fabs(trackerTrack->qoverp() - combinedTrack->qoverp());
152  const double threshold = nSigma * trackerTrack->qoverpError();
153  return delta > threshold ? make_pair(trackerTrack, reco::Muon::InnerTrack)
154  : make_pair(combinedTrack, reco::Muon::CombinedTrack);
155 }
HLTSiStripMonitoring_cff.nSigma
nSigma
Definition: HLTSiStripMonitoring_cff.py:151
dumpMFGeometry_cfg.delta
delta
Definition: dumpMFGeometry_cfg.py:25
Muon
Definition: Muon.py:1
ewkMuLumiMonitorDQM_cfi.ptThreshold
ptThreshold
Definition: ewkMuLumiMonitorDQM_cfi.py:13
DiMuonV_cfg.threshold
threshold
Definition: DiMuonV_cfg.py:101
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:46

◆ sigmaSwitch() [2/2]

reco::Muon::MuonTrackTypePair muon::sigmaSwitch ( const reco::Muon &  muon,
const double  nSigma = 2.,
const double  ptThreshold = 200. 
)
inline

Definition at line 57 of file MuonCocktails.h.

References HLTSiStripMonitoring_cff::nSigma, ewkMuLumiMonitorDQM_cfi::ptThreshold, and sigmaSwitch().

59  {
60  return muon::sigmaSwitch(muon.globalTrack(), muon.innerTrack(), nSigma, ptThreshold);
61  }
HLTSiStripMonitoring_cff.nSigma
nSigma
Definition: HLTSiStripMonitoring_cff.py:151
muon::sigmaSwitch
reco::Muon::MuonTrackTypePair sigmaSwitch(const reco::TrackRef &combinedTrack, const reco::TrackRef &trackerTrack, const double nSigma=2., const double ptThreshold=200.)
Definition: MuonCocktails.cc:139
ewkMuLumiMonitorDQM_cfi.ptThreshold
ptThreshold
Definition: ewkMuLumiMonitorDQM_cfi.py:13

◆ tevOptimized() [1/2]

reco::Muon::MuonTrackTypePair muon::tevOptimized ( const reco::TrackRef &  combinedTrack,
const reco::TrackRef &  trackerTrack,
const reco::TrackRef &  tpfmsTrack,
const reco::TrackRef &  pickyTrack,
const reco::TrackRef &  dytTrack,
const double  ptThreshold = 200.,
const double  tune1 = 17.,
const double  tune2 = 40.,
const double  dptcut = 0.25 
)

Definition at line 9 of file MuonCocktails.cc.

References reco::Muon::CombinedTrack, reco::Muon::DYT, dqmdumpme::first, mps_fire::i, reco::Muon::InnerTrack, reco::Muon::Picky, TtFullHadEvtBuilder_cfi::prob, ewkMuLumiMonitorDQM_cfi::ptThreshold, reco::Muon::TPFMS, trackProbability(), and validateGeometry_cfg::valid.

Referenced by MuonIdProducer::makeMuon(), MuonsFromRefitTracksProducer::produce(), PFMuonAlgo::reconstructMuon(), and tevOptimized().

17  {
18  const unsigned int nAlgo = 5;
19 
20  // Array for convenience below.
21  const reco::Muon::MuonTrackTypePair refit[nAlgo] = {make_pair(trackerTrack, reco::Muon::InnerTrack),
22  make_pair(combinedTrack, reco::Muon::CombinedTrack),
23  make_pair(tpfmsTrack, reco::Muon::TPFMS),
24  make_pair(pickyTrack, reco::Muon::Picky),
25  make_pair(dytTrack, reco::Muon::DYT)};
26 
27  // Calculate the log(tail probabilities). If there's a problem,
28  // signify this with prob == 0. The current problems recognized are:
29  // the track being not available, whether the (re)fit failed or it's
30  // just not in the event, or if the (re)fit ended up with no valid
31  // hits.
32  double prob[nAlgo] = {0., 0., 0., 0., 0.};
33  bool valid[nAlgo] = {false, false, false, false, false};
34 
35  double dptmin = 1.;
36 
37  if (dptcut > 0) {
38  for (unsigned int i = 0; i < nAlgo; ++i)
39  if (refit[i].first.isNonnull())
40  if (refit[i].first->ptError() / refit[i].first->pt() < dptmin)
41  dptmin = refit[i].first->ptError() / refit[i].first->pt();
42 
43  if (dptmin > dptcut)
44  dptcut = dptmin + 0.15;
45  }
46 
47  for (unsigned int i = 0; i < nAlgo; ++i)
48  if (refit[i].first.isNonnull()) {
49  valid[i] = true;
50  if (refit[i].first->numberOfValidHits() &&
51  (refit[i].first->ptError() / refit[i].first->pt() < dptcut || dptcut < 0))
52  prob[i] = muon::trackProbability(refit[i].first);
53  }
54 
55  // Start with picky.
56  int chosen = 3;
57 
58  // If there's a problem with picky, make the default one of the
59  // other tracks. Try TPFMS first, then global, then tracker-only.
60  if (prob[3] == 0.) {
61  // split so that passing dptcut<0 recreates EXACTLY the old tuneP behavior
62  if (dptcut > 0) {
63  if (prob[4] > 0.)
64  chosen = 4;
65  else if (prob[0] > 0.)
66  chosen = 0;
67  else if (prob[2] > 0.)
68  chosen = 2;
69  else if (prob[1] > 0.)
70  chosen = 1;
71  } else {
72  if (prob[2] > 0.)
73  chosen = 2;
74  else if (prob[1] > 0.)
75  chosen = 1;
76  else if (prob[0] > 0.)
77  chosen = 0;
78  }
79  }
80 
81  // Now the algorithm: switch from picky to dyt if the difference is lower than a tuned value. Then
82  // switch from picky to tracker-only if the
83  // difference, log(tail prob(picky)) - log(tail prob(tracker-only))
84  // is greater than a tuned value. Then compare the
85  // so-picked track to TPFMS in the same manner using another tuned
86  // value.
87  if (prob[4] > 0. && prob[3] > 0.) {
88  if (refit[3].first->pt() > 0 && refit[4].first->pt() > 0 &&
89  (refit[4].first->ptError() / refit[4].first->pt() - refit[3].first->ptError() / refit[3].first->pt()) <= 0)
90  chosen = 4; // dyt
91  }
92 
93  if (prob[0] > 0. && prob[chosen] > 0. && (prob[chosen] - prob[0]) > tune1)
94  chosen = 0;
95  if (prob[2] > 0. && (prob[chosen] - prob[2]) > tune2)
96  chosen = 2;
97 
98  // Sanity checks
99  if (chosen == 4 && !valid[4])
100  chosen = 3;
101  if (chosen == 3 && !valid[3])
102  chosen = 2;
103  if (chosen == 2 && !valid[2])
104  chosen = 1;
105  if (chosen == 1 && !valid[1])
106  chosen = 0;
107 
108  // Done. If pT of the chosen track (or pT of the tracker track) is below the threshold value, return the tracker track.
109  if (valid[chosen] && refit[chosen].first->pt() < ptThreshold && prob[0] > 0.)
110  return make_pair(trackerTrack, reco::Muon::InnerTrack);
111  if (trackerTrack->pt() < ptThreshold && prob[0] > 0.)
112  return make_pair(trackerTrack, reco::Muon::InnerTrack);
113 
114  // Return the chosen track (which can be the global track in
115  // very rare cases).
116  return refit[chosen];
117 }
mps_fire.i
i
Definition: mps_fire.py:429
TtFullHadEvtBuilder_cfi.prob
prob
Definition: TtFullHadEvtBuilder_cfi.py:33
validateGeometry_cfg.valid
valid
Definition: validateGeometry_cfg.py:21
ewkMuLumiMonitorDQM_cfi.ptThreshold
ptThreshold
Definition: ewkMuLumiMonitorDQM_cfi.py:13
muon::trackProbability
double trackProbability(const reco::TrackRef track)
Definition: MuonCocktails.cc:122
reco::Muon::MuonTrackTypePair
std::pair< TrackRef, Muon::MuonTrackType > MuonTrackTypePair
Definition: Muon.h:38

◆ tevOptimized() [2/2]

reco::Muon::MuonTrackTypePair muon::tevOptimized ( const reco::Muon &  muon,
const double  ptThreshold = 200.,
const double  tune1 = 17.,
const double  tune2 = 40.,
const double  dptcut = 0.25 
)
inline

Definition at line 31 of file MuonCocktails.h.

References ewkMuLumiMonitorDQM_cfi::ptThreshold, and tevOptimized().

35  {
36  return tevOptimized(muon.globalTrack(),
37  muon.innerTrack(),
38  muon.tpfmsTrack(),
39  muon.pickyTrack(),
40  muon.dytTrack(),
41  ptThreshold,
42  tune1,
43  tune2,
44  dptcut);
45  }
muon::tevOptimized
reco::Muon::MuonTrackTypePair tevOptimized(const reco::Muon &muon, const double ptThreshold=200., const double tune1=17., const double tune2=40., const double dptcut=0.25)
Definition: MuonCocktails.h:31
ewkMuLumiMonitorDQM_cfi.ptThreshold
ptThreshold
Definition: ewkMuLumiMonitorDQM_cfi.py:13

◆ TMR()

reco::Muon::MuonTrackTypePair muon::TMR ( const reco::TrackRef &  trackerTrack,
const reco::TrackRef &  fmsTrack,
const double  tune = 4. 
)

Definition at line 160 of file MuonCocktails.cc.

References reco::Muon::InnerTrack, edm::Ref< C, T, F >::isNonnull(), reco::Muon::None, reco::Muon::TPFMS, and trackProbability().

162  {
163  double probTK = 0;
164  double probFMS = 0;
165 
166  if (trackerTrack.isNonnull() && trackerTrack->numberOfValidHits())
167  probTK = muon::trackProbability(trackerTrack);
168  if (fmsTrack.isNonnull() && fmsTrack->numberOfValidHits())
169  probFMS = muon::trackProbability(fmsTrack);
170 
171  bool TKok = probTK > 0;
172  bool FMSok = probFMS > 0;
173 
174  if (TKok && FMSok) {
175  if (probFMS - probTK > tune)
176  return make_pair(trackerTrack, reco::Muon::InnerTrack);
177  else
178  return make_pair(fmsTrack, reco::Muon::TPFMS);
179  } else if (FMSok)
180  return make_pair(fmsTrack, reco::Muon::TPFMS);
181  else if (TKok)
182  return make_pair(trackerTrack, reco::Muon::InnerTrack);
183  else
184  return make_pair(reco::TrackRef(), reco::Muon::None);
185 }
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:238
muon::trackProbability
double trackProbability(const reco::TrackRef track)
Definition: MuonCocktails.cc:122

◆ trackProbability()

double muon::trackProbability ( const reco::TrackRef  track)

Definition at line 122 of file MuonCocktails.cc.

References createfilelist::int, dqm-mbProfile::log, and HLT_2024v12_cff::track.

Referenced by HistogramProbabilityEstimator::probability(), tevOptimized(), tevOptimizedTMR(), and TMR().

122  {
123  int nDOF = (int)track->ndof();
124  if (nDOF > 0 && track->chi2() > 0) {
125  return -log(TMath::Prob(track->chi2(), nDOF));
126  } else {
127  return 0.0;
128  }
129 }
HLT_2024v12_cff.track
track
Definition: HLT_2024v12_cff.py:9556

Variable Documentation

◆ selectionTypeStringToEnumMap

const SelectionTypeStringToEnum muon::selectionTypeStringToEnumMap[]
static

Definition at line 66 of file MuonSelectors.h.

Referenced by selectionTypeFromString().

◆ selectorStringToEnumMap

const SelectorStringToEnum muon::selectorStringToEnumMap[]
static

Definition at line 107 of file MuonSelectors.h.

Referenced by selectorFromString().

Generated for CMSSW Reference Manual by doxygen 1.8.14