EgammaPhotonProducers/python/reducedEgamma_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoEgamma.EgammaIsolationAlgos.interestingEgammaIsoDetIdsSequence_cff import interestingEgammaIsoHCALSel
0004 from RecoJets.Configuration.CaloTowersES_cfi import *
0005
0006 reducedEgamma = cms.EDProducer("ReducedEGProducer",
0007   keepPfSuperclusterPtMin = cms.double(5.),
0008   keepPfSuperclusterAbsetaMax = cms.double(2.5),
0009   relinkSuperclusterPtMin = cms.double(99999.), # no SC linking
0010   keepPhotons = cms.string("hadTowOverEm()<0.15 && pt>10 && (pt>14 || chargedHadronIso()<10)"), #keep in output
0011   slimRelinkPhotons = cms.string("hadTowOverEm()<0.15 && pt>10 && (pt>14 || chargedHadronIso()<10)"), #keep only slimmed SuperCluster plus seed cluster
0012   relinkPhotons = cms.string("(r9()>0.8 || chargedHadronIso()<20 || chargedHadronIso()<0.3*pt())"), #keep all associated clusters/rechits/conversions
0013   keepOOTPhotons = cms.string("pt>10"), #keep in output
0014   slimRelinkOOTPhotons = cms.string("pt>10"), #keep only slimmed SuperCluster plus seed cluster
0015   relinkOOTPhotons = cms.string("(r9()>0.8)"), #keep all associated clusters/rechits/conversions
0016   keepGsfElectrons = cms.string(""), #keep in output
0017   slimRelinkGsfElectrons = cms.string(""), #keep only slimmed SuperCluster plus seed cluster
0018   relinkGsfElectrons = cms.string("pt>5"), #keep all associated clusters/rechits/conversions
0019   pflowSuperclusters = cms.InputTag("superClusterMerger"),
0020   photons = cms.InputTag("gedPhotons"),
0021   ootPhotons = cms.InputTag("ootPhotons"),
0022   gsfElectrons = cms.InputTag("gedGsfElectrons"),
0023   conversions = cms.InputTag("allConversions"),
0024   gsfTracks = cms.InputTag("electronGsfTracks"),
0025   singleConversions = cms.InputTag("particleFlowEGamma"),
0026   barrelEcalHits = cms.InputTag("reducedEcalRecHitsEB"),
0027   endcapEcalHits = cms.InputTag("reducedEcalRecHitsEE"),
0028   preshowerEcalHits = cms.InputTag("reducedEcalRecHitsES"),
0029   hbheHits= cms.InputTag("reducedHcalRecHits","hbhereco"),
0030   photonsPFValMap = cms.InputTag("particleBasedIsolation","gedPhotons"),
0031   gsfElectronsPFValMap = cms.InputTag("particleBasedIsolation","gedGsfElectrons"),
0032   photonIDSources = cms.VInputTag(),
0033   photonIDOutput = cms.vstring(),
0034   gsfElectronIDSources = cms.VInputTag(),
0035   gsfElectronIDOutput = cms.vstring(),
0036   photonFloatValueMapSources = cms.VInputTag(),
0037   photonFloatValueMapOutput = cms.vstring(),
0038   ootPhotonFloatValueMapSources = cms.VInputTag(),
0039   ootPhotonFloatValueMapOutput = cms.vstring(),
0040   hiPhotonIsolationMapInput = cms.InputTag(""),
0041   hiPhotonIsolationMapOutput = cms.string(""),
0042   gsfElectronFloatValueMapSources = cms.VInputTag(),
0043   gsfElectronFloatValueMapOutput = cms.vstring(),
0044   applyPhotonCalibOnData = cms.bool(False),
0045   applyPhotonCalibOnMC = cms.bool(False),
0046   applyGsfElectronCalibOnData = cms.bool(False),
0047   applyGsfElectronCalibOnMC = cms.bool(False),
0048   photonCalibEnergySource = cms.InputTag(""),
0049   photonCalibEnergyErrSource = cms.InputTag(""),
0050   gsfElectronCalibEnergySource = cms.InputTag(""),
0051   gsfElectronCalibEnergyErrSource = cms.InputTag(""),
0052   gsfElectronCalibEcalEnergySource = cms.InputTag(""),
0053   gsfElectronCalibEcalEnergyErrSource = cms.InputTag(""),
0054   hcalHitSel = interestingEgammaIsoHCALSel
0055 )
0056
0057 superClusterMerger = cms.EDProducer("EgammaSuperClusterMerger",
0058   src = cms.VInputTag(
0059     cms.InputTag("particleFlowSuperClusterECAL:particleFlowSuperClusterECALBarrel"),
0060     cms.InputTag("particleFlowSuperClusterECAL:particleFlowSuperClusterECALEndcapWithPreshower"),
0061   )
0062 )
0063
0064 from Configuration.Eras.Modifier_phase2_common_cff import phase2_common
0065 phase2_common.toModify(reducedEgamma,
0066         preshowerEcalHits = "",
0067 )
0068
0069 from RecoEgamma.EgammaPhotonProducers.reducedEgamma_tools import calibrateReducedEgamma
0070 from Configuration.ProcessModifiers.run2_miniAOD_UL_cff import run2_miniAOD_UL
0071 modifyReducedEGammaRun2MiniAOD = (run2_miniAOD_UL).makeProcessModifier(calibrateReducedEgamma)
0072
0073 from Configuration.ProcessModifiers.pp_on_AA_cff import pp_on_AA
0074
0075 pp_on_AA.toModify(
0076     reducedEgamma,
0077     ootPhotons = "",
0078     keepPhotons = "pt>15 && abs(eta)<2.5",
0079     slimRelinkPhotons = "pt>15 && abs(eta)<2.5",
0080     relinkPhotons = "pt>15 && abs(eta)<2.5",
0081     hiPhotonIsolationMapInput = "photonIsolationHIProducerppGED",
0082     hiPhotonIsolationMapOutput = "photonIsolationHIProducerppGED"
0083 )
0084
0085 from Configuration.ProcessModifiers.run2_miniAOD_pp_on_AA_103X_cff import run2_miniAOD_pp_on_AA_103X
0086 run2_miniAOD_pp_on_AA_103X.toModify(
0087     reducedEgamma,
0088     photonsPFValMap = "pfEGammaToCandidateRemapperCleaned:photons",
0089     gsfElectronsPFValMap = "pfEGammaToCandidateRemapperCleaned:electrons"
0090 )
0091
0092 from Configuration.ProcessModifiers.egamma_lowPt_exclusive_cff import egamma_lowPt_exclusive
0093 egamma_lowPt_exclusive.toModify(reducedEgamma, keepPfSuperclusterPtMin = 1.0, keepPhotons = "", slimRelinkPhotons = "", relinkPhotons = "", relinkGsfElectrons = "")