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.2 by konec, Sun Dec 4 15:36:14 2011 UTC vs.
Revision 1.4 by konec, Mon Jan 23 23:47:03 2012 UTC

# Line 47 | Line 47 | bool BestMuonFinder::run(const edm::Even
47    
48    for (reco::MuonCollection::const_iterator im = muons->begin(); im != muons->end(); ++im) {
49  
50 //
51 // FIXME - UNCOMMENT SELECTED LINES
52 //
53
54
50      if (!im->isTrackerMuon() || !im->innerTrack().isNonnull()) continue;
51      if (hMuPtVsEta) hMuPtVsEta->Fill(im->innerTrack()->eta(), im->innerTrack()->pt());
52      if (im->innerTrack()->pt() < theConfig.getParameter<double>("minPt")) continue;
53      if (fabs(im->innerTrack()->eta()) >  theConfig.getParameter<double>("maxEta")) continue;
54 +    if (im->innerTrack()->dxy(reference) >  theConfig.getParameter<double>("maxTIP")) continue;
55      if (hMuChi2Tk) hMuChi2Tk->Fill(im->innerTrack()->normalizedChi2());
56      if (im->innerTrack()->normalizedChi2() >  theConfig.getParameter<double>("maxChi2Tk")) continue;
57      if (hMuChi2Gl && im->isGlobalMuon()) hMuChi2Gl->Fill(im->combinedMuon()->normalizedChi2());
62 /*
63    if (im->innerTrack()->dxy(reference) >  theConfig.getParameter<double>("maxTIP")) continue;
58      if (    theConfig.getParameter<bool>("requireOuterTrack")){
59 <          if(!im->isStandAloneMuon() || !im->outerTrack().isNonnull())continue;
59 >        if(!im->isStandAloneMuon() || !im->outerTrack().isNonnull())continue;
60            if(im->standAloneMuon()->normalizedChi2() >  theConfig.getParameter<double>("maxChi2Sa")) continue;
61      }
62      if ( theConfig.getParameter<bool>("requireGlobalTrack")) {
63            if(!im->isGlobalMuon() ||  !im->globalTrack().isNonnull()) continue;  
64            if(im->combinedMuon()->normalizedChi2() >  theConfig.getParameter<double>("maxChi2Mu")) continue;
65      }
66 < */
66 >
67 >    if (im->numberOfMatchedStations() <  theConfig.getParameter<int>("minNumberOfMatchedStations")) continue;
68 >    if (! (im->track()->hitPattern().numberOfValidPixelHits() > 0)) continue;
69 >
70   //
71   // TMP TIGHT SELECTION FOR IVAN
72 <  if (! im->isGlobalMuon()) continue;
73 <  if (! (im->globalTrack()->normalizedChi2() < 10)) continue;
74 <  if (! (im->globalTrack()->hitPattern().numberOfValidMuonHits() > 0)) continue;
75 <  if (! (im->numberOfMatchedStations() > 1)) continue;
76 <  if (! (fabs(im->innerTrack()->dxy(reference)) < 0.2)) continue;
77 <  if (! (im->track()->hitPattern().numberOfValidPixelHits() > 0)) continue;
78 <  if (! (im->track()->hitPattern().numberOfValidTrackerHits() > 10)) continue;
72 > //  if (! im->isGlobalMuon()) continue;
73 > //  if (! (im->globalTrack()->normalizedChi2() < 10)) continue;
74 > //  if (! (im->globalTrack()->hitPattern().numberOfValidMuonHits() > 0)) continue;
75 > //  if (! (im->numberOfMatchedStations() > 1)) continue;
76 > //  if (! (fabs(im->innerTrack()->dxy(reference)) < 0.2)) continue;
77 > //  if (! (im->track()->hitPattern().numberOfValidPixelHits() > 0)) continue;
78 > //  if (! (im->track()->hitPattern().numberOfValidTrackerHits() > 10)) continue;
79   // TMP END OF TIGHT MUON SELECTION FROM IVAN
80  
81      //remove muons without valid hits in tk and mu system
# Line 100 | Line 97 | bool BestMuonFinder::run(const edm::Even
97      if (nCSCHits==0 && hMuHitsRPCvsDT) hMuHitsRPCvsDT->Fill(nDTHits,nRPCHits);
98      }
99  
103 /*
100      if (nTrackerHits< theConfig.getParameter<int>("minNumberTrackerHits")) continue;
101      if ( nRPCHits < theConfig.getParameter<int>("minNumberRpcHits")) continue;
102      if ( nDTHits + nCSCHits < theConfig.getParameter<int>("minNumberDtCscHits")  ) continue;
107 */
103  
104      if (theMuon) theUnique = false;
105      if (!theMuon || (im->track()->pt() > theMuon->track()->pt()) ) theMuon = &(*im);

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines