| 46 |
|
ev.getByLabel( theConfig.getParameter<std::string>("muonColl"), muons); |
| 47 |
|
|
| 48 |
|
for (reco::MuonCollection::const_iterator im = muons->begin(); im != muons->end(); ++im) { |
| 49 |
+ |
|
| 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; |
| 53 |
– |
if (im->innerTrack()->dxy(reference) > theConfig.getParameter<double>("maxTIP")) continue; |
| 54 |
|
if (hMuChi2Tk) hMuChi2Tk->Fill(im->innerTrack()->normalizedChi2()); |
| 55 |
|
if (im->innerTrack()->normalizedChi2() > theConfig.getParameter<double>("maxChi2Tk")) continue; |
| 56 |
|
if (hMuChi2Gl && im->isGlobalMuon()) hMuChi2Gl->Fill(im->combinedMuon()->normalizedChi2()); |
| 57 |
< |
if ( theConfig.getParameter<bool>("requireOuterTrack")){ |
| 58 |
< |
if(!im->isStandAloneMuon() || !im->outerTrack().isNonnull())continue; |
| 59 |
< |
if(im->standAloneMuon()->normalizedChi2() > theConfig.getParameter<double>("maxChi2Sa")) continue; |
| 60 |
< |
} |
| 61 |
< |
if ( theConfig.getParameter<bool>("requireGlobalTrack")) { |
| 62 |
< |
if(!im->isGlobalMuon() || !im->globalTrack().isNonnull()) continue; |
| 63 |
< |
if(im->combinedMuon()->normalizedChi2() > theConfig.getParameter<double>("maxChi2Mu")) continue; |
| 64 |
< |
} |
| 57 |
> |
if (im->numberOfMatchedStations() < theConfig.getParameter<int>("minNumberOfMatchedStations")) continue; |
| 58 |
> |
|
| 59 |
> |
// |
| 60 |
> |
// TMP TIGHT SELECTION FOR IVAN |
| 61 |
> |
// if (! im->isGlobalMuon()) continue; |
| 62 |
> |
// if (! (im->globalTrack()->normalizedChi2() < 10)) continue; |
| 63 |
> |
// if (! (im->globalTrack()->hitPattern().numberOfValidMuonHits() > 0)) continue; |
| 64 |
> |
// if (! (im->numberOfMatchedStations() > 1)) continue; |
| 65 |
> |
// if (! (fabs(im->innerTrack()->dxy(reference)) < 0.2)) continue; |
| 66 |
> |
// if (! (im->track()->hitPattern().numberOfValidPixelHits() > 0)) continue; |
| 67 |
> |
// if (! (im->track()->hitPattern().numberOfValidTrackerHits() > 10)) continue; |
| 68 |
> |
// TMP END OF TIGHT MUON SELECTION FROM IVAN |
| 69 |
|
|
| 70 |
|
//remove muons without valid hits in tk and mu system |
| 71 |
|
int nTrackerHits =0; |
| 83 |
|
nDTHits = hp.numberOfValidMuonDTHits(); |
| 84 |
|
nCSCHits = hp.numberOfValidMuonCSCHits(); |
| 85 |
|
if (nDTHits==0 && hMuHitsRPCvsCSC) hMuHitsRPCvsCSC->Fill(nCSCHits,nRPCHits); |
| 82 |
– |
// if (nDTHits==0 && hMuHitsCSCvsEta) hMuHitsCSCvsEta->Fill(fabs(im->innerTrack()->eta()),nCSCHits); |
| 86 |
|
if (nCSCHits==0 && hMuHitsRPCvsDT) hMuHitsRPCvsDT->Fill(nDTHits,nRPCHits); |
| 87 |
|
} |
| 88 |
|
|
