(r2) SUS13007LeptonSelection < CMS

CMS Web>NewPhysics>CmsSusy>CMSSusyTWiki>SUS13007LeptonSelection (revision 2) (raw view)~~EditAttach~~
%TWISTY{
mode="div"
showlink="Show..."
hidelink="Hide"
showimgleft="%ICONURLPATH{toggleopen-small}%"
hideimgleft="%ICONURLPATH{toggleclose-small}%"
}%

   * 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.

   * 3) Efficiencies of leptonic and hadronic legs are shown to be not strongly correlated, thus they are measured separately. The lepton legs are mesured via Z T&P in events that are acquired by the analysis triggers (therefore with correctly biased kinematics). 
      * The triggers used for electrons: 
         * HLT CleanPFHT350 Ele5 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50]
         * HLT CleanPFNoPUHT350 Ele5 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50]
         * HLT CleanPFHT300 Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50]
         * HLT CleanPFNoPUHT300 Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50]
         * HLT CleanPFHT300 Ele[40,60] CaloIdVT CaloIsoVL TrkIdT (control trigger)
         * HLT CleanPFNoPUHT300 Ele[40,60] CaloIdVT CaloIsoVL TrkIdT (control trigger)
      * and for muons: 
         * HLT PFHT350 Mu15 PFMET[45,50]
         * HLT PFNoPUHT350 Mu15 PFMET[45,50]
         * HLT PFHT400 Mu5 PFMET[45,50]
         * HLT PFNoPUHT400 Mu5 PFMET[45,50]
         * HLT Mu[40,60] PFHT350 (control trigger)
         * HLT Mu[40,60] PFNoPUHT350 (control trigger)

   * 4) All analysis selections are fully documented here: https://twiki.cern.ch/twiki/bin/view/CMS/SingleLepton2012#Electron_Selection https://twiki.cern.ch/twiki/bin/view/CMS/SingleLepton2012#Muon_Selection

Summary tables:

%TWISTY{mode="div" showlink="Electron selections" hidelink="Hide" showimgleft="%ICONURLPATH{toggleopen-small}%" hideimgleft="%ICONURLPATH{toggleclose-small}%" }%

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

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

%ENDTWISTY%

%TWISTY{mode="div" showlink="Muon selections" hidelink="Hide" showimgleft="%ICONURLPATH{toggleopen-small}%" hideimgleft="%ICONURLPATH{toggleclose-small}%" }%

Muons are selected using the POG tight muon recommended selection at <a href="https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideMuonId">https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideMuonId</a>

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

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

%ENDTWISTY%

   * 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 &plusmn; 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)

%ENDTWISTY%

* Link to RA4 Analysis: [[http://cms.cern.ch/iCMS/analysisadmin/cadi?ancode=SUS-11-015][paper draft]]<br /> * 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<br /> * 2) PromptReco until reReco are validated.<br /> * 3) For ElectronHad datasets:<br /> * HLT Ele10 CaloIdL CaloIsoVL TrkIdVL TrkIsoVL HT200<br /> * Ele10 CaloIdL CaloIsoVL TrkIdT TrkIsoVL HT200<br /> * HLT Ele10 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT200<br /> * HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT200<br /> * HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250<br /> * HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250 PFMHT25<br /> * HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250 PFMHT40<br /> * HLT Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL HT250 PFMHT50<br /> * HLT HT350 Ele30 CaloIdT TrkIdT (pT(e) background control trigger)<br /> * HLT HT400 Ele60 CaloIdT TrkIdT (pT(e) background control trigger)<br /> * HLT HT450 Ele60 CaloIdT TrkIdT (pT(e) background control trigger)<br /> * 4a) Signal electron requirements: Base on VBTF 80 CIso<br />%TWISTY{mode="div" showlink="Table" hidelink="Hide" showimgleft="%ICONURLPATH{toggleopen-small}%" hideimgleft="%ICONURLPATH{toggleclose-small}%"<br />}%<br />| *Quantity* | *PAT Object and Member Function* | *Cut* | *Comment* |<br />| p_T | GsfElectron =&gt; p4().Pt() | &ge; 20 GeV | ElecHad Trigger at 15 GeV is not 100% efficient at 20 GeV |<br />| abs(eta) | SuperCluster =&gt; Eta() | &le; 2.5 , &gt;1.566 or &lt;1.4442 | Veto EBEE transition region |<br />| Rel. Isolation | Barrel = ( electron-&gt;dr03TkSumPt() + max(0., electron-&gt;dr03EcalRecHitSumEt() - 1.) + electron-&gt;dr03HcalTowerSumEt() ) / electron-&gt;p4().Pt(), Endcaps = ( electron-&gt;dr03TkSumPt() + electron-&gt;dr03EcalRecHitSumEt() + electron-&gt;dr03HcalTowerSumEt() ) / electron-&gt;p4().Pt() ) | VBTF80: barrel&lt;0.07, endcap&lt;0.06 | |<br />| Identification | GsfElectron =&gt; hadronicOverEm(), deltaPhiSuperClusterTrackAtVtx(), deltaEtaSuperClusterTrackAtVtx(), sigmaIetaIeta() | VBTF80 |<br />| hadOverEM | GsfElectron =&gt; hadronicOverEm() | VBTF80: barrel&lt;0.04, endcap&lt;0.025 | |<br />| abs(deltaPhi) | GsfElectron =&gt; deltaPhiSuperClusterTrackAtVtx() | VBTF80: barrel&lt;0.06, endcap&lt;0.03 | |<br />| abs(deltaEta) | GsfElectron =&gt; deltaEtaSuperClusterTrackAtVtx() | VBTF80: barrel&lt;0.004, endcap&lt;0.007 | In 3_6_x this is not applied in endcap |<br />| sigmaIetaIeta | GsfElectron =&gt; sigmaIetaIeta() | VBTF80: barrel&lt;0.01, endcap&lt;0.03 | |<br />| Conversion rejection | gsfTrack() =&gt; trackerExpectedHitsInner() = 0, abs(dist) &ge; 0.02 && abs(delta cot(theta)) &ge; 0.02 | VBTF80 |<br />| Primary Vertex | abs(PV.z() - lepton.vertex.z()) | &lt;1.0 | |<br />| Impact Parameter | abs(e.gsfTrack().dxy(BS)) | &lt; 0.02 | from Muon |

-- Main.ViktorVeszpremi - 2013-02-15
Topic revision: r2 - 2013-02-15 - ViktorVeszpremi
CMS (CERN)
Copyright &© by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback