(r9) ElectronIDRun2 < Main

Tags: , view all tags

---+!! Efficiency and fake rate of the cut-based electron ID for Run 2.

---++ Samples

Signal samples from the CSA production can be found in [[https://cmsweb.cern.ch/das/][DAS]], by typing in the query window:
   * dataset = /RSGravToZZ*/*Spring14dr*/AODSIM

To implement the cut-based electron identification we produced private MINIAODSIM samples with few events (100K), using the common ID and Isolation Toolkit [[https://twiki.cern.ch/twiki/bin/viewauth/CMS/CommonIDAndIsolationFW][CITK]]. 

To access the private samples from lxplus, you have to obtain a valid grid proxy (=voms-proxy-init --voms cms=). The =process.source= in your configure file should look like the following: 

<verbatim>
process.source = cms.Source("PoolSource",
    fileNames = cms.untracked.vstring('root://cmsxrootd.fnal.gov//store/user/jruizvar/RSGravToZZ/M1000/RSGravToZZ_kMpl01_M-1000_Tune4C_13TeV-pythia8_MINIAODSIM_PU20bx25_1.root')
)
</verbatim>

---++ Matching categories

MINIAODSIM samples contain information about generated and reconstructed particles. Our matching algorithm identifies the closest GenElectron for every RecoElectron, as well as the particle ID of the  GenElectron's ancestor. This procedure allows to define four matching categories as described in the following table:

| *1* | ==Unmatched== | deltaR(closest gen electron, reco electron) > 0.02 |
| *2* | ==True electron from W== | deltaR(closest gen electron, reco electron) < 0.02 AND  (ancestor PID == 24) |
| *3* | ==True electron from Z== | deltaR(closest gen electron, reco electron) < 0.02 AND  (ancestor PID == 23) |
| *4* | ==Anything else== | What remains from above |

---++ Isolation

Particle flow isolation with delta beta correction is calculated as follows:

<verbatim>
    reco::GsfElectron::PflowIsolationVariables pfIso = el->pfIsolationVariables();
    absiso = pfIso.sumChargedHadronPt + std::max(0.0 , pfIso.sumNeutralHadronEt + pfIso.sumPhotonEt - 0.5 * pfIso.sumPUPt );
    relIsoWithDBeta = absiso/pt;
</verbatim>

Internally, the isolation variables use the standard recommended cone size =deltaR = 0.3=. Details [[https://twiki.cern.ch/twiki/bin/viewauth/CMS/EgammaCutBasedIdentification][here]].

  * Electron isolation after different cut-based selections 
<br /> <img width="1300" alt="iso.png" src="%ATTACHURLPATH%/iso.png" height="600" />

---++ Efficiency and fake rate

The efficiency is the number of signal electrons passing the ID divided by the total. The fake rate was estimated in a similar way but using the non-signal electrons. 

<br /> <img width="400" alt="formulas.png" src="%ATTACHURLPATH%/formulas.png" height="150" />

To access the [[https://cmssdt.cern.ch/SDT/doxygen/CMSSW_7_1_7/doc/html/d2/d1f/classpat_1_1Electron.html#aa0fbe7bf4bacf8035991013a67e1a40e][electronID]], which is a member function of the =pat::Electron= class, we used a line like the following:

=int passID = ( el->electronID("cutBasedElectronID-CSA14-PU20bx25-V0-standalone-veto") > 0.5 );=

---++ Working points

The cut-based electron identification recommended by EGM for Run 2 analyses, determines four working points: *veto*, *loose*, *medium*, and *tight*. A detailed description [[https://twiki.cern.ch/twiki/bin/viewauth/CMS/CutBasedElectronIdentificationRun2][here]].

---++ Preliminary results

To estimate the efficiency and fake rate of the electron identification, we classified the electrons according to its matching category and the *passID* associated to the *cutBasedElectronID-CSA14* working points.

The classification for each working point is displayed as a TH2 histogram with two bins in the Y axis associated to the passID variable, and four bins in the X axis corresponding to the matching categories. The number of entries inside each box of the histogram indicates the number of electrons in a given matching category with a given passID.   

   * Electron classification, global efficiency, and fake rate using sample RSGravToZZ 1 TeV
<br /> <img width="1000" alt="RSGravToZZ-M1000_matchId.png" src="%ATTACHURLPATH%/RSGravToZZ-M1000_matchId.png" height="750" />

   * Electron classification, global efficiency, and fake rate using sample RSGravToZZ 4.5 TeV
<br /> <img width="1000" alt="RSGravToZZ-M4500_matchId.png" src="%ATTACHURLPATH%/RSGravToZZ-M4500_matchId.png" height="750" />

Efficiency and fake rate as function of pT and etaSC are shown below:

   * Efficiency and fake rate as function of electron's pT
<br /> <img width="1300" alt="eff_pt.png" src="%ATTACHURLPATH%/eff_pt.png" height="600" />

   * Efficiency and fake rate as function of etaSC 
<br /> <img width="1300" alt="eff_eta.png" src="%ATTACHURLPATH%/eff_eta.png" height="600" />


---++ Comments

%COMMENT%

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
png	RSGravToZZ-M1000_matchId.png	r1	manage	34.6 K	2014-09-19 - 11:40	JoseRuiz	Classification of electrons in sample RSGravToZZ-M1000
png	RSGravToZZ-M1000_matchIdA.png	r1	manage	34.3 K	2014-09-18 - 15:23	JoseRuiz	Classification of electrons in sample RSGravToZZ-M1000
png	RSGravToZZ-M4500_matchId.png	r1	manage	34.3 K	2014-09-19 - 11:41	JoseRuiz	Classification of electrons in sample RSGravToZZ-M4500
png	RSGravToZZ-M4500_matchIdA.png	r1	manage	33.9 K	2014-09-18 - 15:24	JoseRuiz	Classification of electrons in sample RSGravToZZ-M4500
png	eff_eta.png	r2 r1	manage	77.4 K	2014-09-19 - 11:41	JoseRuiz	Efficiency and fake rate as function of super cluster eta
png	eff_pt.png	r2 r1	manage	76.6 K	2014-09-19 - 11:42	JoseRuiz	Efficiency and fake rate as function of electron's pT
png	formulas.png	r1	manage	117.6 K	2014-09-15 - 09:36	JoseRuiz	Efficiency and fake rate formulas
png	iso.png	r3 r2 r1	manage	50.1 K	2014-09-19 - 12:34	JoseRuiz	Electron isolation with delta beta correction

Topic revision: r9 - 2014-09-26 - JoseRuiz

antalya escort bursa escort eskisehir escort istanbul escort izmir escort