Electron Identification using Likelihood

Complete: 3

Introduction

The following variables are currently used to discriminate between real and fake electrons in the likelihood:

  • Energy of closest BasicCluster to track impact point at ECAL / outermost track momentum
  • DeltaEta between SuperCluster position and track direction at vertex extrapolated to ECAL assuming no radiation
  • DeltaPhi between SuperCluster position and track direction at vertex
  • Ratio of energy in HCAL behind SuperCluster to SuperCluster energy
  • Combination in a linear Fisher discriminant of cluster shape variables:
    • Energy in 3x3 crystals / energy in 5x5 crystals
    • SigmaEtaEta cluster shape covariance

Specialized pdfs are for:

  • ECAL barrel, ECAL endcap
  • electron pT < 15GeV / >15GeV
  • electron classes (only for signal). Two possibilities are allowed:
    • split by "fullclass": golden, bigbrem, narrow, showering
    • split by "class": non-showering (golden+bigbrem+narrow all together weighted by their fractions), showering

Species: only two species are currently considered:

  • Electrons: pdfs taken from pure di-electron sample, particle gun
  • Jets: QCD samples, pdfs are the weighted sum of pdf's for di-jets from different pT bins from enriched samples. The weight is the product of the L1 filter efficiency and the cross section.

Configuration

Configuration can be done changing the .cfi file:

  • EgammaAnalysis/ElectronIDESSources/data/likelihoodESetup.cfi

In the default configuration, signal species (electrons) pdf's are split by "fullclass" and background species is not split. For the signal species, the possibilities are:

  • split by "fullclass". The configuration is:
    string signalWeightSplitting = "fullclass"
    bool splitSignalPdfs = true
  • split by "class":
    string signalWeightSplitting = "class"
    bool splitSignalPdfs = true
  • do not split at all:
    string signalWeightSplitting = "fullclass"
    bool splitSignalPdfs = false

In principle also for background the same procedure is available. Due to the low MC statistics available in the Bigbrem/Narrow classes, right now only the unsplit pdf's are available in the DB. Configuration is:

  • do not split at all:
    string backgroundWeightSplitting = "fullclass"
    bool splitBackgroundPdfs = false

Setup of the code

Tags in 17X

cvs co -r V00-00-07-01   EgammaAnalysis/ElectronIDAlgos
cvs co -r V00-00-04          EgammaAnalysis/ElectronIDProducers
cvs co -r V00-00-03-01   EgammaAnalysis/ElectronIDESSources
cvs co -r V00-01-01          CondCore/EgammaPlugins                           
cvs co -r V00-02-00          CondFormats/EgammaObjects 

Tags in 16X

cvs co -r V00-00-07-01   EgammaAnalysis/ElectronIDAlgos
cvs co -r V00-00-04          EgammaAnalysis/ElectronIDProducers
cvs co -r V00-00-03-01   EgammaAnalysis/ElectronIDESSources
cvs co -r V00-01-01          CondCore/EgammaPlugins
cvs co -r V02-06-02          CondFormats/DataRecord
cvs co -r V00-02-00          CondFormats/EgammaObjects
cvs co -r V00-01-03          CondFormats/PhysicsToolsObjects
cvs co -r V00-05-00          DataFormats/EgammaReco
cvs co -r V00-05-07          RecoEcal/EgammaClusterProducers
cvs co -r V00-03-00          RecoEcal/EgammaCoreTools

for releases 1_6_8, 1_6_9 apply the following backportings:

cvs up -r 1.10 DataFormats/EgammaReco/interface/BasicCluster.h
cvs up -r 1.10 RecoEcal/EgammaCoreTools/interface/PositionCalc.h

Tags in 131HLT6 (obsolete!)

cvs co -r V00-00-04      EgammaAnalysis/ElectronIDAlgos
cvs co -r V00-00-01      EgammaAnalysis/ElectronIDESSources
cvs co -r V00-05-00      DataFormats/EgammaReco
cvs co -r V00-01-00      CondCore/EgammaPlugins
cvs co -r V00-01-10      CondFormats/BTauObjects
cvs co -r V02-06-02      CondFormats/DataRecord
cvs co -r V00-01-00      CondFormats/EgammaObjects
cvs co -r V00-03-04      RecoEcal/EgammaClusterAlgos
cvs co -r V00-05-07      RecoEcal/EgammaClusterProducers
cvs co -r V00-03-00      RecoEcal/EgammaCoreTools
cvs co -r shr-8aug07     PhysicsTools/StatPatternRecognition

Run with CRAB

Since the SQLite DB file is not yet distributed on sites, it has to be exported in the CRAB InputSandBox. To do it, change, in EgammaAnalysis/ElectronIDESSources/data/likelihoodPdfsDB.cfi:

replace CondDBCommon.connect = "sqlite_file:electronIdLikelihoodTkIsolated.db"    
replace CondDBCommon.catalog = "file:electronIdLikelihoodTkIsolatedCRAB.xml"

and, in the CRAB cfg file:

## files to be put in InputSandBox, full path or in the current directory
## (wildcard * are allowed): comma separated list
additional_input_files = /afs/cern.ch/user/e/emanuele/public/4Likelihood/PDFsSQLite/CMSSW_1_6_8/electronIdLikelihoodTkIsolated.db, 
/afs/cern.ch/user/e/emanuele/public/4Likelihood/PDFsSQLite/CMSSW_1_6_8/electronIdLikelihoodTkIsolatedCRAB.xml

Test the likelihood

Examples of analyzers tested in CMSSW_168 can be found here:

> cvs co -d emanuele/EgammaTools UserCode/emanuele/EgammaTools
> scramv1 b
> cmsRun emanuele/EgammaTools/test/testElectronIDWithLH.cfg


This topic: CMSPublic > CMSCommunicationsGroup > CMSCommunicationsProjects > WebHome > SWGuide > SWGuideEgamma > SWGuideElectronID > SWGuideLikelihoodElectronID
Topic revision: r11 - 2008-02-22 - EmanueleDiMarco
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback