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 |
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); |