EgammaIsolationAlgos/python/interestingEgammaIsoDetIdsSequence_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoEgamma.EgammaIsolationAlgos.interestingEleIsoDetIdModule_cff import *
0004 from RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff import *
0005
0006
0007 import RecoEgamma.EgammaIsolationAlgos.interestingEleIsoDetIdModule_cff
0008 interestingGedEleIsoDetIdEB = RecoEgamma.EgammaIsolationAlgos.interestingEleIsoDetIdModule_cff.interestingEleIsoDetId.clone(
0009     recHitsLabel  = 'ecalRecHit:EcalRecHitsEB',
0010     emObjectLabel = 'gedGsfElectrons',
0011     etCandCut     = 0.0,
0012     energyCut     = 0.095,
0013     etCut         = 0.0,
0014     outerRadius   = 0.6,
0015     innerRadius   = 0.0
0016 )
0017
0018 import RecoEgamma.EgammaIsolationAlgos.interestingEleIsoDetIdModule_cff
0019 interestingGedEleIsoDetIdEE = RecoEgamma.EgammaIsolationAlgos.interestingEleIsoDetIdModule_cff.interestingEleIsoDetId.clone(
0020     recHitsLabel  = 'ecalRecHit:EcalRecHitsEE',
0021     emObjectLabel = 'gedGsfElectrons',
0022     etCandCut     = 0.0,
0023     energyCut     = 0.0,
0024     etCut         = 0.110,
0025     outerRadius   = 0.6,
0026     innerRadius   = 0.0
0027 )
0028
0029 import RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff
0030 interestingGedGamIsoDetIdEB = RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff.interestingGamIsoDetId.clone(
0031     recHitsLabel  = 'ecalRecHit:EcalRecHitsEB',
0032     emObjectLabel = 'gedPhotons',
0033     etCandCut     = 0.0,
0034     energyCut     = 0.095,
0035     etCut         = 0.0,
0036     outerRadius   = 0.6,
0037     innerRadius   = 0.0
0038 )
0039
0040 import RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff
0041 interestingGedGamIsoDetIdEE = RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff.interestingGamIsoDetId.clone(
0042     recHitsLabel  = 'ecalRecHit:EcalRecHitsEE',
0043     emObjectLabel = 'gedPhotons',
0044     etCandCut     = 0.0,
0045     energyCut     = 0.0,
0046     etCut         = 0.110,
0047     outerRadius   = 0.6,
0048     innerRadius   = 0.0
0049 )
0050 ## OOT photons
0051 interestingOotGamIsoDetIdEB = interestingGedGamIsoDetIdEB.clone(emObjectLabel = 'ootPhotons')
0052 interestingOotGamIsoDetIdEE = interestingGedGamIsoDetIdEE.clone(emObjectLabel = 'ootPhotons')
0053
0054 import RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff
0055 interestingGamIsoDetIdEB = RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff.interestingGamIsoDetId.clone(
0056     recHitsLabel  = 'ecalRecHit:EcalRecHitsEB',
0057     emObjectLabel = 'photons',
0058     etCandCut     = 0.0,
0059     energyCut     = 0.095,
0060     etCut         = 0.0,
0061     outerRadius   = 0.6,
0062     innerRadius   = 0.0
0063 )
0064
0065 import RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff
0066 interestingGamIsoDetIdEE = RecoEgamma.EgammaIsolationAlgos.interestingGamIsoDetIdModule_cff.interestingGamIsoDetId.clone(
0067     recHitsLabel  = 'ecalRecHit:EcalRecHitsEE',
0068     emObjectLabel = 'photons',
0069     etCandCut     = 0.0,
0070     energyCut     = 0.0,
0071     etCut         = 0.110,
0072     outerRadius   = 0.6,
0073     innerRadius   = 0.0
0074 )
0075
0076 import RecoEgamma.EgammaIsolationAlgos.interestingGedEgammaIsoHCALDetId_cfi
0077 interestingEgammaIsoHCALSel = cms.PSet(
0078     maxDIEta           = cms.int32(5),
0079     maxDIPhi           = cms.int32(5),
0080     minEnergyHB        = cms.double(0.8),
0081     minEnergyHEDepth1  = cms.double(0.1),
0082     minEnergyHEDefault = cms.double(0.2),
0083 )
0084 interestingGedEgammaIsoHCALDetId = RecoEgamma.EgammaIsolationAlgos.interestingGedEgammaIsoHCALDetId_cfi.interestingGedEgammaIsoHCALDetId.clone(
0085     recHitsLabel       = "hbhereco",
0086     elesLabel          = "gedGsfElectrons",
0087     phosLabel          = "gedPhotons",
0088     superClustersLabel = "particleFlowEGamma",
0089     minSCEt            = 20,
0090     minEleEt           = 20,
0091     minPhoEt           = 20,
0092     hitSelection       = interestingEgammaIsoHCALSel
0093 )
0094
0095 ## OOT Photons
0096 interestingOotEgammaIsoHCALDetId = interestingGedEgammaIsoHCALDetId.clone(
0097     phosLabel          = "ootPhotons",
0098     elesLabel          = "",
0099     superClustersLabel = ""
0100 )
0101
0102 import RecoEgamma.EgammaIsolationAlgos.interestingEgammaIsoESDetIdModule_cff
0103 interestingOotEgammaIsoESDetId = RecoEgamma.EgammaIsolationAlgos.interestingEgammaIsoESDetIdModule_cff.interestingEgammaIsoESDetId.clone(
0104     eeClusToESMapLabel = "particleFlowClusterOOTECAL",
0105     ecalPFClustersLabel= "particleFlowClusterOOTECAL",
0106     phosLabel          = "ootPhotons",
0107     elesLabel          = "",
0108     superClustersLabel = "",
0109     minSCEt            = 500,
0110     minEleEt           = 20,
0111     minPhoEt           = 20,
0112     maxDR              = 0.4
0113 )
0114
0115 interestingEgammaIsoDetIdsTask = cms.Task(
0116     interestingGedEleIsoDetIdEB ,
0117     interestingGedEleIsoDetIdEE ,
0118     interestingGedGamIsoDetIdEB ,
0119     interestingGedGamIsoDetIdEE ,
0120     interestingOotGamIsoDetIdEB ,
0121     interestingOotGamIsoDetIdEE ,
0122     interestingGamIsoDetIdEB ,
0123     interestingGamIsoDetIdEE ,
0124     interestingGedEgammaIsoHCALDetId,
0125     interestingOotEgammaIsoHCALDetId,
0126     interestingOotEgammaIsoESDetId
0127 )
0128 interestingEgammaIsoDetIds = cms.Sequence(interestingEgammaIsoDetIdsTask)
0129
0130 _pp_on_AA_interestingEgammaIsoDetIdsTask = interestingEgammaIsoDetIdsTask.copy()
0131 _pp_on_AA_interestingEgammaIsoDetIdsTask.remove(interestingOotGamIsoDetIdEB)
0132 _pp_on_AA_interestingEgammaIsoDetIdsTask.remove(interestingOotGamIsoDetIdEE)
0133 _pp_on_AA_interestingEgammaIsoDetIdsTask.remove(interestingOotEgammaIsoHCALDetId)
0134 _pp_on_AA_interestingEgammaIsoDetIdsTask.remove(interestingOotEgammaIsoESDetId)
0135
0136 from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
0137 pp_on_AA.toReplaceWith(interestingEgammaIsoDetIdsTask, _pp_on_AA_interestingEgammaIsoDetIdsTask)
0138
0139 from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
0140
0141 egamma_lowPt_exclusive.toModify(interestingGedEgammaIsoHCALDetId,
0142                    minSCEt = 1.0, #default 20
0143                    minEleEt= 1.0, #default 20
0144                minPhoEt= 1.0 #default 20
0145 )
0146
0147 from Configuration.Eras.Modifier_run3_HB_cff import run3_HB
0148 run3_HB.toModify(interestingEgammaIsoHCALSel,
0149                  minEnergyHB = 0.1)