Skip to content

Commit

Permalink
Cutbased photonID Fall17 V2
Browse files Browse the repository at this point in the history
  • Loading branch information
@swagata87 @guitargeek
swagata87 authored and guitargeek committed Dec 3, 2018
1 parent 4090b10 commit daebbcd
Show file tree
Hide file tree
Showing 4 changed files with 181 additions and 0 deletions.
9 changes: 9 additions & 0 deletions ...onIdentification/data/Fall17/effAreaPhotons_cone03_pfChargedHadrons_90percentBased_V2.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,9 @@
# The constants are based on 90% contours of isolation.
# |eta| min |eta| max effective area
0.0000 1.0000 0.0404
1.0000 1.4790 0.0396
1.4790 2.0000 0.0364
2.0000 2.2000 0.0353
2.2000 2.3000 0.0351
2.3000 2.4000 0.0347
2.4000 5.0000 0.0337
9 changes: 9 additions & 0 deletions ...onIdentification/data/Fall17/effAreaPhotons_cone03_pfNeutralHadrons_90percentBased_V2.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,9 @@
# The constants are based on 90% contours of isolation.
# |eta| min |eta| max effective area
0.0000 1.0000 0.067
1.0000 1.4790 0.105
1.4790 2.0000 0.079
2.0000 2.2000 0.023
2.2000 2.3000 0.008
2.3000 2.4000 0.003
2.4000 5.0000 0.014
9 changes: 9 additions & 0 deletions ...ma/PhotonIdentification/data/Fall17/effAreaPhotons_cone03_pfPhotons_90percentBased_V2.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,9 @@
# The constants are based on 90% contours of isolation.
# |eta| min |eta| max effective area
0.0000 1.0000 0.111
1.0000 1.4790 0.095
1.4790 2.0000 0.062
2.0000 2.2000 0.084
2.2000 2.3000 0.107
2.3000 2.4000 0.121
2.4000 5.0000 0.147
154 changes: 154 additions & 0 deletions RecoEgamma/PhotonIdentification/python/Identification/cutBasedPhotonID_Fall17_100X_V2_cff.py
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,154 @@
from PhysicsTools.SelectorUtils.centralIDRegistry import central_id_registry

import FWCore.ParameterSet.Config as cms

# Common functions and classes for ID definition are imported here:
from RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_tools \
import ( WorkingPoint_V2,
IsolationCutInputs,
configureVIDCutBasedPhoID_V5 )

# The cut values are taken from the twiki:
# https://twiki.cern.ch/twiki/bin/viewauth/CMS/CutBasedPhotonIdentificationRun2
# (where they may not stay, if a newer version of cuts becomes available for these
# conditions)

# See also the presentation explaining these working points (this will not change):
# https://indico.cern.ch/event/732974/contributions/3072291/attachments/1685029/2709189/PhotonIDStudy.pdf


#
# First, define cut values
#

# Loose working point Barrel and Endcap
idName = "cutBasedPhotonID-Fall17-100X-V2-loose"
WP_Loose_EB = WorkingPoint_V2(
idName , # idName
0.0313 , # hOverECut
0.0108 , # full5x5_SigmaIEtaIEtaCut
# Isolation cuts are generally absIso < C1 + pt*C2, except for NeuHad is < C1 + pt*C2 + pt*pt*C3
1.272 , # absPFChaHadIsoWithEACut_C1
0 , # absPFChaHadIsoWithEACut_C2
8.3828 , # absPFNeuHadIsoWithEACut_C1
0.01512 , # absPFNeuHadIsoWithEACut_C2
0.00002259 , # absPFNeuHadIsoWithEACut_C3
4.2891 , # absPFPhoIsoWithEACut_C1
0.004017 # absPFPhoIsoWithEACut_C2
)
WP_Loose_EE = WorkingPoint_V2(
idName , #idName
0.0431 , # hOverECut
0.0281 , # full5x5_SigmaIEtaIEtaCut
# Isolation cuts are generally absIso < C1 + pt*C2, except for NeuHad is < C1 + pt*C2 + pt*pt*C3
1.4449 , # absPFChaHadIsoWithEACut_C1
0.00 , # absPFChaHadIsoWithEACut_C2
13.6894 , # absPFNeuHadIsoWithEACut_C1
0.0117 , # absPFNeuHadIsoWithEACut_C2
0.000023 , # absPFNeuHadIsoWithEACut_C3
4.2969 , # absPFPhoIsoWithEACut_C1
0.0037 # absPFPhoIsoWithEACut_C2
)

# Medium working point Barrel and Endcap
idName = "cutBasedPhotonID-Fall17-100X-V2-medium"
WP_Medium_EB = WorkingPoint_V2(
idName , # idName
0.02 , # hOverECut
0.01 , # full5x5_SigmaIEtaIEtaCut
# Isolation cuts are generally absIso < C1 + pt*C2, except for NeuHad is < C1 + pt*C2 + pt*pt*C3
0.6127 , # absPFChaHadIsoWithEACut_C1
0.00 , # absPFChaHadIsoWithEACut_C2
3.6265 , # absPFNeuHadIsoWithEACut_C1
0.01512 , # absPFNeuHadIsoWithEACut_C2
0.00002259 , # absPFNeuHadIsowithEACut_C3
2.3692 , # absPFPhoIsoWithEACut_C1
0.004017 # absPFPhoIsoWithEACut_C2
)

WP_Medium_EE = WorkingPoint_V2(
idName , #idName
0.0248 , # hOverECut
0.0271 , # full5x5_SigmaIEtaIEtaCut
# Isolation cuts are generally absIso < C1 + pt*C2, except for NeuHad is < C1 + pt*C2 + pt*pt*C3
0.7981 , # absPFChaHadIsoWithEACut_C1
0.00 , # absPFChaHadIsoWithEACut_C2
2.3737 , # absPFNeuHadIsoWithEACut_C1
0.0117 , # absPFNeuHadIsoWithEACut_C2
0.000023 , # absPFNeuHadIsowithEACut_C3
3.1760 , # absPFPhoIsoWithEACut_C1
0.0037 # absPFPhoIsoWithEACut_C2
)

# Tight working point Barrel and Endcap
idName = "cutBasedPhotonID-Fall17-100X-V2-tight"
WP_Tight_EB = WorkingPoint_V2(
idName , # idName
0.02 , # hOverECut
0.01 , # full5x5_SigmaIEtaIEtaCut
# Isolation cuts are generally absIso < C1 + pt*C2, except for NeuHad is < C1 + pt*C2 + pt*pt*C3
0.0939 , # absPFChaHadIsoWithEACut_C1
0.00 , # absPFChaHadIsoWithEACut_C2
3.6195 , # absPFNeuHadIsoWithEACut_C1
0.01512 , # absPFNeuHadIsoWithEACut_C2
0.00002259 , # absPFNeuHadIsowithEACut_C3
2.0018 , # absPFPhoIsoWithEACut_C1
0.004017 # absPFPhoIsoWithEACut_C2
)

WP_Tight_EE = WorkingPoint_V2(
idName , #idName
0.0219 , # hOverECut
0.0268 , # full5x5_SigmaIEtaIEtaCut
# Isolation cuts are generally absIso < C1 + pt*C2, except for NeuHad is < C1 + pt*C2 + pt*pt*C3
0.2629 , # absPFChaHadIsoWithEACut_C1
0.00 , # absPFChaHadIsoWithEACut_C2
2.1436 , # absPFNeuHadIsoWithEACut_C1
0.0117 , # absPFNeuHadIsoWithEACut_C2
0.000023 , # absPFNeuHadIsowithEACut_C3
3.0133 , # absPFPhoIsoWithEACut_C1
0.0037 # absPFPhoIsoWithEACut_C2
)


# Second, define where to find the precomputed isolations and what effective
# areas to use for pile-up correction
isoInputs = IsolationCutInputs(
# chHadIsolationMapName
'photonIDValueMapProducer:phoChargedIsolation' ,
# chHadIsolationEffAreas
"RecoEgamma/PhotonIdentification/data/Fall17/effAreaPhotons_cone03_pfChargedHadrons_90percentBased_V2.txt",
# neuHadIsolationMapName
'photonIDValueMapProducer:phoNeutralHadronIsolation' ,
# neuHadIsolationEffAreas
"RecoEgamma/PhotonIdentification/data/Fall17/effAreaPhotons_cone03_pfNeutralHadrons_90percentBased_V2.txt" ,
# phoIsolationMapName
"photonIDValueMapProducer:phoPhotonIsolation" ,
# phoIsolationEffAreas
"RecoEgamma/PhotonIdentification/data/Fall17/effAreaPhotons_cone03_pfPhotons_90percentBased_V2.txt"
)

#
# Finally, set up VID configuration for all cuts
#
cutBasedPhotonID_Fall17_100X_V2_loose = configureVIDCutBasedPhoID_V5 ( WP_Loose_EB, WP_Loose_EE, isoInputs)
cutBasedPhotonID_Fall17_100X_V2_medium = configureVIDCutBasedPhoID_V5 ( WP_Medium_EB, WP_Medium_EE, isoInputs)
cutBasedPhotonID_Fall17_100X_V2_tight = configureVIDCutBasedPhoID_V5 ( WP_Tight_EB, WP_Tight_EE, isoInputs)

## The MD5 sum numbers below reflect the exact set of cut variables
# and values above. If anything changes, one has to
# 1) comment out the lines below about the registry,
# 2) run "calculateMD5 <this file name> <one of the VID config names just above>
# 3) update the MD5 sum strings below and uncomment the lines again.
#

central_id_registry.register(cutBasedPhotonID_Fall17_100X_V2_loose.idName,
'4578dfcceb0bfd1ba5ac28973c843fd0')
central_id_registry.register(cutBasedPhotonID_Fall17_100X_V2_medium.idName,
'28b186c301061395f394a81266c8d7de')
central_id_registry.register(cutBasedPhotonID_Fall17_100X_V2_tight.idName,
'6f4f0ed6a8bf2de8dcf0bc3349b0546d')

cutBasedPhotonID_Fall17_100X_V2_loose.isPOGApproved = cms.untracked.bool(True)
cutBasedPhotonID_Fall17_100X_V2_medium.isPOGApproved = cms.untracked.bool(True)
cutBasedPhotonID_Fall17_100X_V2_tight.isPOGApproved = cms.untracked.bool(True)

0 comments on commit daebbcd

Please sign in to comment.