TWiki> CMS Web>NewPhysics>CmsSusy>CMSSusyTWiki>SUS13007LeptonSelection (revision 2)EditAttach
Show... Hide

  • 1) SUS-13-007: covers the single lepton + btag using the Lepton Spectrum (LS) and Delta Phi (DP) methods. Lepton studies are documented in AN-2013-035. Targeting Moriond.

  • 2) AOD rereco is used for periods A, B, and part C; prompt reco is used for part C and D.

Summary tables:

Electron selections Hide

Quantity PAT Object and Member Function barrel cut endcap cut comment*
type pat::Electron      
p_T ele->pt() ≥ 20 GeV    
abs(eta) 
(ele->isEB() || ele->isEE()) && fabs(ele->superCluster()->eta()) <= 2.5
     
Rel. Isolation charged + max (0., photon + neutral - (*eleRho)*Aeff) <0.15 same  
hadOverEM 
ele->hadronicOverEm()
 <0.12 <0.10  
abs(deltaPhi) 
fabs(ele->deltaPhiSuperClusterTrackAtVtx())
 <0.06 <0.03  
abs(deltaEta) 
abs(ele->deltaEtaSuperClusterTrackAtVtx())
 <0.004 <0.007  
sigmaIetaIeta 
ele->scSigmaIEtaIEta()
 <0.01 <0.03  
Conversion rejection 
!ConversionTools::hasMatchedConversion(*gsfel,hConversions,beamSpotPosition)
 True True Use the gsfel obtained with originalObjectRef; don't use the pat::Electron
missingHits 
ele->gsfTrack()->trackerExpectedHitsInner().numberOfHits()
 <= 1 <= 1  
dxy 
fabs(ele->gsfTrack()->dxy(vertexPosition))
 < 0.02 same vertexPosition is the position of the first good vertex
dz 
fabs(ele->gsfTrack()->dz(vertexPosition))
 < 0.1 same  
abs(reco - PF) PT 
fabs(electron->pt() - (*ipfEl).pt())
< 10 GeV same 
(*ipfEl).particleId() == reco::PFCandidate::e && (*ipfEl) > 10.

veto electron selection

Quantity PAT Object and Member Function barrel cut endcap cut comment*
type pat::Electron      
p_T ele->pt() ≥ 15 GeV    
abs(eta) 
ele->isEB() || ele->isEE()
      
Rel. Isolation charged + max (0., photon + neutral - (*eleRho)*Aeff) <0.15 same  
hadOverEM 
ele->hadronicOverEm()
 <0.15 no cut  
abs(deltaPhi) 
fabs(ele->deltaPhiSuperClusterTrackAtVtx())
 <0.8 <0.7  
abs(deltaEta) 
abs(ele->deltaEtaSuperClusterTrackAtVtx())
 <0.007 <0.01  
sigmaIetaIeta 
ele->scSigmaIEtaIEta()
 <0.01 <0.03  
dxy 
fabs(ele->gsfTrack()->dxy(vertexPosition))
 < 0.04 same vertexPosition is the position of the first good vertex
dz 
fabs(ele->gsfTrack()->dz(vertexPosition))
 < 0.2 same  

Muon selections Hide

Muons are selected using the POG tight muon recommended selection at https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideMuonId

Quantity PAT Object and Member Function Cut Comment
muon type 
pat::Muon
    
ID 
muon -> isGlobalMuon()
True  
ID 
muon -> isPFMuon()
 True  
pT 
muon -> pt()
 ≥ 20 GeV  
abs(eta) 
muon-> eta()
≤ 2.4  
pf Rel. Isolation 
pfRelIso =  (muon->pfIsolationR03().sumChargedHadronPt 
+ max(0., muon->pfIsolationR03().sumNeutralHadronEt 
+ muon->pfIsolationR03().sumPhotonEt - 0.5*muon->pfIsolationR03().sumPUPt ) ) 
/ muon->pt()
 < 0.12 Don't use pat::Muon->isoDeposit()->depositWithin(0.3)
normChi2 
muon->globalTrack()->chi2() / muon->globalTrack()->ndof()
 ≤ 10  
nValMuonHits 
muon->globalTrack()->hitPattern().numberOfValidMuonHits()
 > 0  
numMatchedStations 
muon->numberOfMatchedStations()
 > 1  
pixelHits 
muon->innerTrack()->hitPattern().numberOfValidPixelHits()
> 0  
numTrackerLayersWithMeasurement 
muon->track()->hitPattern().trackerLayersWithMeasurement()
 > 5  
dxy 
fabs(muon->innerTrack()->dxy(vertexPosition))
 < 0.02 vertexPosition is the position of the first good vertex
dz 
fabs(muon->innerTrack()->dz(vertexPosition))
 < 0.5  
abs(reco - PF) 
fabs(muon->pt() - (*ipfMu).pt())
< 5 GeV 
(*ipfMu).particleId() == reco::PFCandidate::mu && (*ipfMu) > 10.

Muon Veto Selection

Quantity PAT Object and Member Function Cut Comment
Mu type 
pat::Muon
    
ID 
muon->isGlobalMuon() || muon->isTrackerMuon()
True  
ID 
muon -> isPFMuon()
True  
pT 
muon -> pt()
≥ 15 GeV  
abs(eta) 
fabs(muon -> eta())
 ≤ 2.5  
pf Rel. Isolation 
pfRelIso =  (muon->pfIsolationR03().sumChargedHadronPt 
+ max(0., muon->pfIsolationR03().sumNeutralHadronEt 
+ muon->pfIsolationR03().sumPhotonEt - 0.5*muon->pfIsolationR03().sumPUPt ) ) 
/ muon->pt()
 < 0.2  
dxy 
fabs(muon->innerTrack()->dxy(vertexPosition))
< 0.2 vertexPosition is the position of the first good vertex
dz 
fabs(muon->innerTrack()->dz(vertexPosition))
< 0.5 vertexPosition is the position of the first good vertex

  • 5) Lepton reconstruction and ID efficiency scale factors between data and MC were determined with the T&P method. Backrgound subtraction was applied in the fit. The scale factors are computed as the ratio of constant fits on the efficiencies over the threshold value.

Selection data/MC (full simulation) data/MC (fast simulation)
Electron ID + Isolation (average) 0.9987 ± 0.0015(stat.)?0.02(sys.) 0.9613 ? 0.0014(stat.)?0.08(sys.)
Barrel   0.9900 ? 0.0007(stat.)?0.03(sys.)
Endcap   0.8725 ? 0.0014(stat.)?0.07(sys.)
Electron GSF tracking 0.9944 ? 0.0003(stat.)?0.02(sys.) 0.9902 ? 0.0004(stat.)?0.01(sys.)
Muon ID1 0.9981 ? 0.0007(stat.)?0.02(sys.) 0.9531 ? 0.0007(stat.)?0.03(sys.)
Muon ID2 1.0015 ? 0.0006(stat.)?0.01(sys.) 0.9902 ? 0.0005(stat.)?0.01(sys.)
Muon track 0.9953 ? 0.0001(stat.)?0.005(sys.) 0.9951 ? 0.0001(stat.)?0.005(sys.)

  • 6)
  • 7)
  • 8)

* Link to RA4 Analysis: paper draft
* 1) Search for supersymmetry in events with a lepton and missing energy, AN Notes: 2011/173, 2011/255, 2011/205, 2011/424, 2011/410, 2011/315
* 2) PromptReco until reReco are validated.
* 3) For ElectronHad datasets:
* HLT Ele10 CaloIdL CaloIsoVL TrkIdVL TrkIsoVL HT200
* Ele10 CaloIdL CaloIsoVL TrkIdT TrkIsoVL HT200
* HLT Ele10 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT200
* HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT200
* HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250
* HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250 PFMHT25
* HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250 PFMHT40
* HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250 PFMHT50
* HLT HT350 Ele30 CaloIdT TrkIdT (pT(e) background control trigger)
* HLT HT400 Ele60 CaloIdT TrkIdT (pT(e) background control trigger)
* HLT HT450 Ele60 CaloIdT TrkIdT (pT(e) background control trigger)
* 4a) Signal electron requirements: Base on VBTF 80 CIso

Table Hide

| Quantity | PAT Object and Member Function | Cut | Comment |
| p_T | GsfElectron => p4().Pt() | ≥ 20 GeV | ElecHad Trigger at 15 GeV is not 100% efficient at 20 GeV |
| abs(eta) | SuperCluster => Eta() | ≤ 2.5 , >1.566 or <1.4442 | Veto EBEE transition region |
| Rel. Isolation | Barrel = ( electron->dr03TkSumPt() + max(0., electron->dr03EcalRecHitSumEt() - 1.) + electron->dr03HcalTowerSumEt() ) / electron->p4().Pt(), Endcaps = ( electron->dr03TkSumPt() + electron->dr03EcalRecHitSumEt() + electron->dr03HcalTowerSumEt() ) / electron->p4().Pt() ) | VBTF80: barrel<0.07, endcap<0.06 | |
| Identification | GsfElectron => hadronicOverEm(), deltaPhiSuperClusterTrackAtVtx(), deltaEtaSuperClusterTrackAtVtx(), sigmaIetaIeta() | VBTF80 |
| hadOverEM | GsfElectron => hadronicOverEm() | VBTF80: barrel<0.04, endcap<0.025 | |
| abs(deltaPhi) | GsfElectron => deltaPhiSuperClusterTrackAtVtx() | VBTF80: barrel<0.06, endcap<0.03 | |
| abs(deltaEta) | GsfElectron => deltaEtaSuperClusterTrackAtVtx() | VBTF80: barrel<0.004, endcap<0.007 | In 3_6_x this is not applied in endcap |
| sigmaIetaIeta | GsfElectron => sigmaIetaIeta() | VBTF80: barrel<0.01, endcap<0.03 | |
| Conversion rejection | gsfTrack() => trackerExpectedHitsInner() = 0, abs(dist) ≥ 0.02 && abs(delta cot(theta)) ≥ 0.02 | VBTF80 |
| Primary Vertex | abs(PV.z() - lepton.vertex.z()) | <1.0 | |
| Impact Parameter | abs(e.gsfTrack().dxy(BS)) | < 0.02 | from Muon |

-- ViktorVeszpremi - 2013-02-15

Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r2 - 2013-02-15 - ViktorVeszpremi
 
This site is powered by the TWiki collaboration platform Powered by Perl This site is powered by the TWiki collaboration platformCopyright &© by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback