ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/L1RpcTriggerAnalysis/src/BestMuonFinder.cc
(Generate patch)

Comparing UserCode/L1RpcTriggerAnalysis/src/BestMuonFinder.cc (file contents):
Revision 1.1 by konec, Mon Nov 28 18:39:48 2011 UTC vs.
Revision 1.5 by konec, Mon May 21 10:50:08 2012 UTC

# Line 13 | Line 13
13   #include "DataFormats/TrackReco/interface/Track.h"
14  
15   #include "UserCode/L1RpcTriggerAnalysis/interface/Utilities.h"
16 + #include "DataFormats/Math/interface/deltaR.h"
17 + #include "DataFormats/Math/interface/deltaPhi.h"
18   #include "TObjArray.h"
19   #include "TH1D.h"
20   #include "TH2D.h"
# Line 46 | Line 48 | bool BestMuonFinder::run(const edm::Even
48    ev.getByLabel( theConfig.getParameter<std::string>("muonColl"), muons);
49    
50    for (reco::MuonCollection::const_iterator im = muons->begin(); im != muons->end(); ++im) {
51 +
52      if (!im->isTrackerMuon() || !im->innerTrack().isNonnull()) continue;
53      if (hMuPtVsEta) hMuPtVsEta->Fill(im->innerTrack()->eta(), im->innerTrack()->pt());
54      if (im->innerTrack()->pt() < theConfig.getParameter<double>("minPt")) continue;
# Line 55 | Line 58 | bool BestMuonFinder::run(const edm::Even
58      if (im->innerTrack()->normalizedChi2() >  theConfig.getParameter<double>("maxChi2Tk")) continue;
59      if (hMuChi2Gl && im->isGlobalMuon()) hMuChi2Gl->Fill(im->combinedMuon()->normalizedChi2());
60      if (    theConfig.getParameter<bool>("requireOuterTrack")){
61 <          if(!im->isStandAloneMuon() || !im->outerTrack().isNonnull())continue;
61 >        if(!im->isStandAloneMuon() || !im->outerTrack().isNonnull())continue;
62            if(im->standAloneMuon()->normalizedChi2() >  theConfig.getParameter<double>("maxChi2Sa")) continue;
63      }
64      if ( theConfig.getParameter<bool>("requireGlobalTrack")) {
# Line 63 | Line 66 | bool BestMuonFinder::run(const edm::Even
66            if(im->combinedMuon()->normalizedChi2() >  theConfig.getParameter<double>("maxChi2Mu")) continue;
67      }
68  
69 +    if (im->numberOfMatchedStations() <  theConfig.getParameter<int>("minNumberOfMatchedStations")) continue;
70 +    if (! (im->track()->hitPattern().numberOfValidPixelHits() > 0)) continue;
71 +
72 + //
73 + // TMP TIGHT SELECTION FOR IVAN
74 + //  if (! im->isGlobalMuon()) continue;
75 + //  if (! (im->globalTrack()->normalizedChi2() < 10)) continue;
76 + //  if (! (im->globalTrack()->hitPattern().numberOfValidMuonHits() > 0)) continue;
77 + //  if (! (im->numberOfMatchedStations() > 1)) continue;
78 + //  if (! (fabs(im->innerTrack()->dxy(reference)) < 0.2)) continue;
79 + //  if (! (im->track()->hitPattern().numberOfValidPixelHits() > 0)) continue;
80 + //  if (! (im->track()->hitPattern().numberOfValidTrackerHits() > 10)) continue;
81 + // TMP END OF TIGHT MUON SELECTION FROM IVAN
82 +
83      //remove muons without valid hits in tk and mu system
84      int nTrackerHits =0;
85      int nRPCHits =0;
# Line 79 | Line 96 | bool BestMuonFinder::run(const edm::Even
96        nDTHits = hp.numberOfValidMuonDTHits();
97        nCSCHits = hp.numberOfValidMuonCSCHits();
98      if (nDTHits==0 && hMuHitsRPCvsCSC) hMuHitsRPCvsCSC->Fill(nCSCHits,nRPCHits);
82 – //    if (nDTHits==0 && hMuHitsCSCvsEta) hMuHitsCSCvsEta->Fill(fabs(im->innerTrack()->eta()),nCSCHits);
99      if (nCSCHits==0 && hMuHitsRPCvsDT) hMuHitsRPCvsDT->Fill(nDTHits,nRPCHits);
100      }
101  
# Line 87 | Line 103 | bool BestMuonFinder::run(const edm::Even
103      if ( nRPCHits < theConfig.getParameter<int>("minNumberRpcHits")) continue;
104      if ( nDTHits + nCSCHits < theConfig.getParameter<int>("minNumberDtCscHits")  ) continue;
105  
90 –    if (theMuon) theUnique = false;
106      if (!theMuon || (im->track()->pt() > theMuon->track()->pt()) ) theMuon = &(*im);
107    }
108 +
109 +  //
110 +  // check if muon is unigue
111 +  //
112 +  if (theMuon) {
113 +    double muonEta = theMuon->innerTrack()->eta();
114 +    double muonPhi = theMuon->innerTrack()->phi();
115 +    for (reco::MuonCollection::const_iterator im = muons->begin(); im != muons->end(); ++im) {
116 +      if (!im->isTrackerMuon() || !im->innerTrack().isNonnull()) continue;
117 +      if (&(*im) == theMuon) continue;
118 +      if ( fabs(reco::deltaPhi(muonPhi, im->innerTrack()->phi())) > theConfig.getParameter<double>("deltaPhiUnique")) continue;
119 +      if ( fabs(muonEta-im->innerTrack()->eta()) > theConfig.getParameter<double>("deltaEtaUnique")) continue;
120 +      theUnique = false;
121 +    }
122 +  }
123    return true;
124   }
125  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines