SUS13007LeptonSelection < CMS

CMS Web>NewPhysics>CmsSusy>CMSSusyTWiki>SUS13007LeptonSelection (2013-02-15, ViktorVeszpremi) (raw view)
%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. Analysis notes: AN-12-390 and AN-12-416. Pre-approval date: Fri, Feb 1, 2013. 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: 
         * <verbatim>HLT CleanPFHT350 Ele5 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50] </verbatim>
         * <verbatim>HLT CleanPFNoPUHT350 Ele5 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50] </verbatim>
         * <verbatim>HLT CleanPFHT300 Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50] </verbatim>
         * <verbatim>HLT CleanPFNoPUHT300 Ele15 CaloIdT CaloIsoVL TrkIdT TrkIsoVL PFMET[45,50] </verbatim>
         * <verbatim>HLT CleanPFHT300 Ele[40,60] CaloIdVT CaloIsoVL TrkIdT </verbatim>(control trigger)
         * <verbatim>HLT CleanPFNoPUHT300 Ele[40,60] CaloIdVT CaloIsoVL TrkIdT </verbatim>(control trigger)
      * and for muons: 
         * <verbatim>HLT PFHT350 Mu15 PFMET[45,50] </verbatim>
         * <verbatim>HLT PFNoPUHT350 Mu15 PFMET[45,50] </verbatim>
         * <verbatim>HLT PFHT400 Mu5 PFMET[45,50] </verbatim>
         * <verbatim>HLT PFNoPUHT400 Mu5 PFMET[45,50] </verbatim>
         * <verbatim>HLT Mu[40,60] PFHT350 </verbatim>(control trigger)
         * <verbatim>HLT Mu[40,60] PFNoPUHT350 </verbatim>(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 (and simulation) 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.)&plusmn;0.02(sys.) | 0.9613 &plusmn; 0.0014(stat.)&plusmn;0.08(sys.) |
|  Barrel | | 0.9900 &plusmn; 0.0007(stat.)&plusmn;0.03(sys.) |
|  Endcap | | 0.8725 &plusmn; 0.0014(stat.)&plusmn;0.07(sys.) |
| Electron GSF tracking | 0.9944 &plusmn; 0.0003(stat.)&plusmn;0.02(sys.) | 0.9902 &plusmn; 0.0004(stat.)&plusmn;0.01(sys.) |
| Muon ID1 | 0.9981 &plusmn; 0.0007(stat.)&plusmn;0.02(sys.) | 0.9531 &plusmn; 0.0007(stat.)&plusmn;0.03(sys.) |
| Muon ID2 | 1.0015 &plusmn; 0.0006(stat.)&plusmn;0.01(sys.) | 0.9902 &plusmn; 0.0005(stat.)&plusmn;0.01(sys.) |
| Muon track | 0.9953 &plusmn; 0.0001(stat.)&plusmn;0.005(sys.) | 0.9951 &plusmn; 0.0001(stat.)&plusmn;0.005(sys.) |

Trigger efficiencies: scale factor is to be applied on Monte Carlo. Individual triggers are listed with scale factors and systematic errors. "Overall" is a value that is compatible with the result of every trigger in the table. In case of the muons, additionally, a correction is to be applied on the scale factor different for barrel and endcap to account for eta dependence.

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

    | Electron trigger path | Scale factor |
    | <verbatim>CleanPFHT350_Ele5_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET45     </verbatim>  | 0.96 &plusmn; 0.04 |
    | <verbatim>CleanPFHT350_Ele5_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET50     </verbatim>  | 0.96 &plusmn; 0.04 |
    | <verbatim>CleanPFNoPUHT350_Ele5_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET45 </verbatim>  | 0.94 &plusmn; 0.07 |
    | <verbatim>CleanPFNoPUHT350_Ele5_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET50 </verbatim>  | 0.95 &plusmn; 0.08 |
    | <verbatim>CleanPFHT300_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET45    </verbatim>  | 0.97 &plusmn; 0.07 |
    | <verbatim>CleanPFHT300_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET50    </verbatim>  | 0.97 &plusmn; 0.07 |
    | <verbatim>CleanPFNoPUHT300_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET45</verbatim>  | 0.98 &plusmn; 0.03 |
    | <verbatim>CleanPFNoPUHT300_Ele15_CaloIdT_CaloIsoVL_TrkIdT_TrkIsoVL_PFMET50</verbatim>  | 0.98 &plusmn; 0.04 |
    | <verbatim>CleanPFHT300_Ele40_CaloIdVT_TrkIdT                              </verbatim>  | 0.96 &plusmn; 0.02 |
    | <verbatim>CleanPFNoPUHT300_Ele40_CaloIdVT_TrkIdT                          </verbatim>  | 0.96 &plusmn; 0.02 |
    | <verbatim>CleanPFHT300_Ele60_CaloIdVT_TrkIdT                              </verbatim>  | 0.96 &plusmn; 0.02 |
    | <verbatim>CleanPFNoPUHT300_Ele60_CaloIdVT_TrkIdT                          </verbatim>  | 0.96 &plusmn; 0.02 |
    | <verbatim>Ele30_CaloIdVT_TrkIdT                                           </verbatim>  | 0.96 &plusmn; 0.02 |
    | <verbatim>Ele25_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30  </verbatim>  | 0.97 &plusmn; 0.02 |
    | Overall                                                          | 0.96 &plusmn; 0.05 |

    | Muon trigger path | Scale factor |
    | <verbatim>IsoMu24                  </verbatim>  | 0.88 &plusmn; 0.02 |
    | <verbatim>IsoMu24_eta2p1           </verbatim>  | 0.90 &plusmn; 0.01 |
    | <verbatim>Mu40_PFHT350             </verbatim>  | 0.91 &plusmn; 0.02 |
    | <verbatim>Mu40_PFNoPUHT350         </verbatim>  | 0.91 &plusmn; 0.02 |
    | <verbatim>Mu60_PFHT350             </verbatim>  | 0.91 &plusmn; 0.02 |
    | <verbatim>Mu60_PFNoPUHT350         </verbatim>  | 0.92 &plusmn; 0.01 |
    | <verbatim>PFHT400_Mu5_PFMET45      </verbatim>  | 0.92 &plusmn; 0.05 |
    | <verbatim>PFHT400_Mu5_PFMET50      </verbatim>  | 0.92 &plusmn; 0.05 |
    | <verbatim>PFNoPUHT400_Mu5_PFMET45  </verbatim>  | 0.93 &plusmn; 0.04 |
    | <verbatim>PFNoPUHT400_Mu5_PFMET50  </verbatim>  | 0.93 &plusmn; 0.04 |
    | <verbatim>PFHT350_Mu15_PFMET45     </verbatim>  | 0.93 &plusmn; 0.05 |
    | <verbatim>PFHT350_Mu15_PFMET50     </verbatim>  | 0.93 &plusmn; 0.05 |
    | <verbatim>PFNoPUHT350_Mu15_PFMET45 </verbatim>  | 0.92 &plusmn; 0.04 |
    | <verbatim>PFNoPUHT350_Mu15_PFMET50 </verbatim>  | 0.92 &plusmn; 0.05 |
    | Overall                  | 0.91 &plusmn; 0.04 |
    | Scale factor for abs(eta) < 0.9 | 1.076 |
    | Scale factor for abs(eta) &gt 0.9 | 0.928 |


%ENDTWISTY%

   * 6) Fakes are attributed to QCD and are derived from data. (eg Section 4.2 in AN-12-390)
   * 7) Both answers are no.
   * 8) Systematic errors on trigger and selection efficiencies are given in 5), more details are in AN-13-035.

%ENDTWISTY%

-- Main.ViktorVeszpremi - 2013-02-15
Topic revision: r3 - 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