File indexing completed on 2023-03-17 11:16:10
0001 import FWCore.ParameterSet.Config as cms
0002 from PhysicsTools.NanoAOD.nano_eras_cff import *
0003 from PhysicsTools.NanoAOD.common_cff import *
0004 from PhysicsTools.NanoAOD.simpleCandidateFlatTableProducer_cfi import simpleCandidateFlatTableProducer
0005 from math import ceil,log
0006
0007
0008 photon_id_modules_WorkingPoints_nanoAOD = cms.PSet(
0009 modules = cms.vstring(
0010 'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_RunIIIWinter22_122X_V1_cff',
0011 'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Winter22_122X_V1_cff',
0012
0013 'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Fall17_94X_V2_cff',
0014 'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Fall17_94X_V2_cff',
0015 ),
0016 WorkingPoints = cms.vstring(
0017 "egmPhotonIDs:cutBasedPhotonID-RunIIIWinter22-122X-V1-loose",
0018 "egmPhotonIDs:cutBasedPhotonID-RunIIIWinter22-122X-V1-medium",
0019 "egmPhotonIDs:cutBasedPhotonID-RunIIIWinter22-122X-V1-tight",
0020 )
0021 )
0022
0023 photon_id_modules_WorkingPoints_nanoAOD_Run2 = cms.PSet(
0024 modules = cms.vstring(
0025 'RecoEgamma.PhotonIdentification.Identification.cutBasedPhotonID_Fall17_94X_V2_cff',
0026 'RecoEgamma.PhotonIdentification.Identification.mvaPhotonID_Fall17_94X_V2_cff',
0027 ),
0028 WorkingPoints = cms.vstring(
0029 "egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V2-loose",
0030 "egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V2-medium",
0031 "egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V2-tight",
0032 )
0033
0034 )
0035
0036
0037 run2_egamma.toModify(photon_id_modules_WorkingPoints_nanoAOD,
0038 modules=photon_id_modules_WorkingPoints_nanoAOD_Run2.modules).\
0039 toModify(photon_id_modules_WorkingPoints_nanoAOD,
0040 WorkingPoints=photon_id_modules_WorkingPoints_nanoAOD_Run2.WorkingPoints)
0041
0042 def make_bitmapVID_docstring(id_modules_working_points_pset):
0043 pset = id_modules_working_points_pset
0044
0045 for modname in pset.modules:
0046 ids = __import__(modname, globals(), locals(), ['idName','cutFlow'])
0047 for name in dir(ids):
0048 _id = getattr(ids,name)
0049 if hasattr(_id,'idName') and hasattr(_id,'cutFlow'):
0050 if (len(pset.WorkingPoints)>0 and _id.idName == pset.WorkingPoints[0].split(':')[-1]):
0051 cut_names = ','.join([cut.cutName.value() for cut in _id.cutFlow])
0052 n_bits_per_cut = int(ceil(log(len(pset.WorkingPoints)+1,2)))
0053 return 'VID compressed bitmap (%s), %d bits per cut'%(cut_names, n_bits_per_cut)
0054 raise ValueError("Something is wrong in the photon ID modules parameter set!")
0055
0056
0057 bitmapVIDForPho = cms.EDProducer("PhoVIDNestedWPBitmapProducer",
0058 src = cms.InputTag("slimmedPhotons"),
0059 srcForID = cms.InputTag("reducedEgamma","reducedGedPhotons"),
0060 WorkingPoints = photon_id_modules_WorkingPoints_nanoAOD.WorkingPoints,
0061 )
0062 _bitmapVIDForPho_docstring = make_bitmapVID_docstring(photon_id_modules_WorkingPoints_nanoAOD)
0063
0064 bitmapVIDForPhoRun2 = bitmapVIDForPho.clone(
0065 WorkingPoints = photon_id_modules_WorkingPoints_nanoAOD_Run2.WorkingPoints,
0066 )
0067 _bitmapVIDForPhoRun2_docstring = make_bitmapVID_docstring(photon_id_modules_WorkingPoints_nanoAOD_Run2)
0068
0069 isoForPho = cms.EDProducer("PhoIsoValueMapProducer",
0070 src = cms.InputTag("slimmedPhotons"),
0071 relative = cms.bool(False),
0072 doQuadratic = cms.bool(True),
0073 rho_PFIso = cms.InputTag("fixedGridRhoFastjetAll"),
0074 QuadraticEAFile_PFIso_Chg = cms.FileInPath("RecoEgamma/PhotonIdentification/data/RunIII_Winter22/effectiveArea_ChgHadronIso_95percentBased.txt"),
0075 QuadraticEAFile_PFIso_ECal = cms.FileInPath("RecoEgamma/PhotonIdentification/data/RunIII_Winter22/effectiveArea_ECalClusterIso_95percentBased.txt"),
0076 QuadraticEAFile_PFIso_HCal = cms.FileInPath("RecoEgamma/PhotonIdentification/data/RunIII_Winter22/effectiveArea_HCalClusterIso_95percentBased.txt"),
0077 )
0078
0079 hOverEForPho = cms.EDProducer("PhoHoverEValueMapProducer",
0080 src = cms.InputTag("slimmedPhotons"),
0081 relative = cms.bool(False),
0082 rho = cms.InputTag("fixedGridRhoFastjetAll"),
0083 QuadraticEAFile_HoverE = cms.FileInPath("RecoEgamma/PhotonIdentification/data/RunIII_Winter22/effectiveArea_coneBasedHoverE_95percentBased.txt"),
0084 )
0085
0086 isoForPhoFall17V2 = isoForPho.clone(
0087 doQuadratic = cms.bool(False),
0088 EAFile_PFIso_Chg = cms.FileInPath("RecoEgamma/PhotonIdentification/data/Fall17/effAreaPhotons_cone03_pfChargedHadrons_90percentBased_V2.txt"),
0089 EAFile_PFIso_Neu = cms.FileInPath("RecoEgamma/PhotonIdentification/data/Fall17/effAreaPhotons_cone03_pfNeutralHadrons_90percentBased_V2.txt"),
0090 EAFile_PFIso_Pho = cms.FileInPath("RecoEgamma/PhotonIdentification/data/Fall17/effAreaPhotons_cone03_pfPhotons_90percentBased_V2.txt"),
0091 )
0092
0093
0094 seedGainPho = cms.EDProducer("PhotonSeedGainProducer", src = cms.InputTag("slimmedPhotons"))
0095
0096 import RecoEgamma.EgammaTools.calibratedEgammas_cff
0097
0098 calibratedPatPhotonsNano = RecoEgamma.EgammaTools.calibratedEgammas_cff.calibratedPatPhotons.clone(
0099 produceCalibratedObjs = False,
0100 correctionFile = cms.string("EgammaAnalysis/ElectronTools/data/ScalesSmearings/Run2016_UltraLegacy_preVFP_RunFineEtaR9Gain"),
0101 )
0102
0103 (run2_egamma_2016 & tracker_apv_vfp30_2016).toModify(
0104 calibratedPatPhotonsNano,
0105 correctionFile = cms.string("EgammaAnalysis/ElectronTools/data/ScalesSmearings/Run2016_UltraLegacy_preVFP_RunFineEtaR9Gain")
0106 )
0107
0108 (run2_egamma_2016 & ~tracker_apv_vfp30_2016).toModify(
0109 calibratedPatPhotonsNano,
0110 correctionFile = cms.string("EgammaAnalysis/ElectronTools/data/ScalesSmearings/Run2016_UltraLegacy_postVFP_RunFineEtaR9Gain"),
0111 )
0112
0113 run2_egamma_2017.toModify(
0114 calibratedPatPhotonsNano,
0115 correctionFile = cms.string("EgammaAnalysis/ElectronTools/data/ScalesSmearings/Run2017_24Feb2020_runEtaR9Gain_v2")
0116 )
0117
0118 run2_egamma_2018.toModify(
0119 calibratedPatPhotonsNano,
0120 correctionFile = cms.string("EgammaAnalysis/ElectronTools/data/ScalesSmearings/Run2018_29Sep2020_RunFineEtaR9Gain")
0121 )
0122
0123 slimmedPhotonsWithUserData = cms.EDProducer("PATPhotonUserDataEmbedder",
0124 src = cms.InputTag("slimmedPhotons"),
0125 parentSrcs = cms.VInputTag("reducedEgamma:reducedGedPhotons"),
0126 userFloats = cms.PSet(
0127 mvaID = cms.InputTag("photonMVAValueMapProducer:PhotonMVAEstimatorRunIIIWinter22v1Values"),
0128 PFIsoChgQuadratic = cms.InputTag("isoForPho:PFIsoChgQuadratic"),
0129 PFIsoAllQuadratic = cms.InputTag("isoForPho:PFIsoAllQuadratic"),
0130 HoverEQuadratic = cms.InputTag("hOverEForPho:HoEForPhoEACorr"),
0131 mvaID_Fall17V2 = cms.InputTag("photonMVAValueMapProducer:PhotonMVAEstimatorRunIIFall17v2Values"),
0132 PFIsoChgFall17V2 = cms.InputTag("isoForPhoFall17V2:PFIsoChg"),
0133 PFIsoAllFall17V2 = cms.InputTag("isoForPhoFall17V2:PFIsoAll"),
0134 ),
0135 userIntFromBools = cms.PSet(
0136 cutBasedID_loose = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-RunIIIWinter22-122X-V1-loose"),
0137 cutBasedID_medium = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-RunIIIWinter22-122X-V1-medium"),
0138 cutBasedID_tight = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-RunIIIWinter22-122X-V1-tight"),
0139 cutBasedID_Fall17V2_loose = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V2-loose"),
0140 cutBasedID_Fall17V2_medium = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V2-medium"),
0141 cutBasedID_Fall17V2_tight = cms.InputTag("egmPhotonIDs:cutBasedPhotonID-Fall17-94X-V2-tight"),
0142 mvaID_WP90 = cms.InputTag("egmPhotonIDs:mvaPhoID-RunIIIWinter22-v1-wp90"),
0143 mvaID_WP80 = cms.InputTag("egmPhotonIDs:mvaPhoID-RunIIIWinter22-v1-wp80"),
0144 mvaID_Fall17V2_WP90 = cms.InputTag("egmPhotonIDs:mvaPhoID-RunIIFall17-v2-wp90"),
0145 mvaID_Fall17V2_WP80 = cms.InputTag("egmPhotonIDs:mvaPhoID-RunIIFall17-v2-wp80"),
0146 ),
0147 userInts = cms.PSet(
0148 VIDNestedWPBitmap = cms.InputTag("bitmapVIDForPho"),
0149 VIDNestedWPBitmapFall17V2 = cms.InputTag("bitmapVIDForPhoRun2"),
0150 seedGain = cms.InputTag("seedGainPho"),
0151
0152 )
0153 )
0154
0155
0156 run2_egamma.toModify(slimmedPhotonsWithUserData.userFloats,
0157 mvaID = None,
0158 PFIsoChgQuadratic = None,
0159 PFIsoAllQuadratic = None,
0160 HoverEQuadratic = None).\
0161 toModify(slimmedPhotonsWithUserData.userIntFromBools,
0162 cutBasedID_loose = None,
0163 cutBasedID_medium = None,
0164 cutBasedID_tight = None,
0165 mvaID_WP90 = None,
0166 mvaID_WP80 = None).\
0167 toModify(slimmedPhotonsWithUserData.userInts,
0168 VIDNestedWPBitmap = None)
0169
0170 run2_egamma.toModify(
0171 slimmedPhotonsWithUserData.userFloats,
0172 ecalEnergyErrPostCorrNew = cms.InputTag("calibratedPatPhotonsNano","ecalEnergyErrPostCorr"),
0173 ecalEnergyPreCorrNew = cms.InputTag("calibratedPatPhotonsNano","ecalEnergyPreCorr"),
0174 ecalEnergyPostCorrNew = cms.InputTag("calibratedPatPhotonsNano","ecalEnergyPostCorr"),
0175 energyScaleUpNew = cms.InputTag("calibratedPatPhotonsNano","energyScaleUp"),
0176 energyScaleDownNew = cms.InputTag("calibratedPatPhotonsNano","energyScaleDown"),
0177 energySigmaUpNew = cms.InputTag("calibratedPatPhotonsNano","energySigmaUp"),
0178 energySigmaDownNew = cms.InputTag("calibratedPatPhotonsNano","energySigmaDown"),
0179 )
0180
0181
0182 finalPhotons = cms.EDFilter("PATPhotonRefSelector",
0183 src = cms.InputTag("slimmedPhotonsWithUserData"),
0184 cut = cms.string("pt > 5 ")
0185 )
0186
0187 photonTable = simpleCandidateFlatTableProducer.clone(
0188 src = cms.InputTag("linkedObjects","photons"),
0189 name= cms.string("Photon"),
0190 doc = cms.string("slimmedPhotons after basic selection (" + finalPhotons.cut.value()+")"),
0191 variables = cms.PSet(P3Vars,
0192 jetIdx = Var("?hasUserCand('jet')?userCand('jet').key():-1", "int16", doc="index of the associated jet (-1 if none)"),
0193 electronIdx = Var("?hasUserCand('electron')?userCand('electron').key():-1", "int16", doc="index of the associated electron (-1 if none)"),
0194 energyErr = Var("getCorrectedEnergyError('regression2')",float,doc="energy error of the cluster from regression",precision=6),
0195 energyRaw = Var("superCluster().rawEnergy()",float,doc="raw energy of photon supercluster", precision=10),
0196 r9 = Var("full5x5_r9()",float,doc="R9 of the supercluster, calculated with full 5x5 region",precision=8),
0197 sieie = Var("full5x5_sigmaIetaIeta()",float,doc="sigma_IetaIeta of the supercluster, calculated with full 5x5 region",precision=8),
0198 sipip = Var("showerShapeVariables().sigmaIphiIphi", float, doc="sigmaIphiIphi of the supercluster", precision=8),
0199 sieip = Var("full5x5_showerShapeVariables().sigmaIetaIphi",float,doc="sigma_IetaIphi of the supercluster, calculated with full 5x5 region",precision=8),
0200 s4 = Var("full5x5_showerShapeVariables().e2x2/full5x5_showerShapeVariables().e5x5",float,doc="e2x2/e5x5 of the supercluster, calculated with full 5x5 region",precision=8),
0201 etaWidth = Var("superCluster().etaWidth()",float,doc="Width of the photon supercluster in eta", precision=8),
0202 phiWidth = Var("superCluster().phiWidth()",float,doc="Width of the photon supercluster in phi", precision=8),
0203 cutBased = Var(
0204 "userInt('cutBasedID_loose')+userInt('cutBasedID_medium')+userInt('cutBasedID_tight')",
0205 "uint8",
0206 doc="cut-based ID bitmap, RunIIIWinter22V1, (0:fail, 1:loose, 2:medium, 3:tight)",
0207 ),
0208 cutBased_Fall17V2 = Var(
0209 "userInt('cutBasedID_Fall17V2_loose')+userInt('cutBasedID_Fall17V2_medium')+userInt('cutBasedID_Fall17V2_tight')",
0210 "uint8",
0211 doc="cut-based ID bitmap, Fall17V2, (0:fail, 1:loose, 2:medium, 3:tight)",
0212 ),
0213 vidNestedWPBitmap = Var(
0214 "userInt('VIDNestedWPBitmap')",
0215 int,
0216 doc="RunIIIWinter22V1 " + _bitmapVIDForPho_docstring
0217 ),
0218 vidNestedWPBitmap_Fall17V2 = Var(
0219 "userInt('VIDNestedWPBitmapFall17V2')",
0220 int,
0221 doc="Fall17V2 " + _bitmapVIDForPhoRun2_docstring
0222 ),
0223 electronVeto = Var("passElectronVeto()",bool,doc="pass electron veto"),
0224 pixelSeed = Var("hasPixelSeed()",bool,doc="has pixel seed"),
0225 mvaID = Var("userFloat('mvaID')",float,doc="MVA ID score, Winter22V1",precision=10),
0226 mvaID_WP90 = Var("userInt('mvaID_WP90')",bool,doc="MVA ID WP90, Winter22V1"),
0227 mvaID_WP80 = Var("userInt('mvaID_WP80')",bool,doc="MVA ID WP80, Winter22V1"),
0228 mvaID_Fall17V2 = Var("userFloat('mvaID_Fall17V2')",float,doc="MVA ID score, Fall17V2",precision=10),
0229 mvaID_Fall17V2_WP90 = Var("userInt('mvaID_Fall17V2_WP90')",bool,doc="MVA ID WP90, Fall17V2"),
0230 mvaID_Fall17V2_WP80 = Var("userInt('mvaID_Fall17V2_WP80')",bool,doc="MVA ID WP80, Fall17V2"),
0231 trkSumPtHollowConeDR03 = Var("trkSumPtHollowConeDR03()",float,doc="Sum of track pT in a hollow cone of outer radius, inner radius", precision=8),
0232 pfPhoIso03 = Var("photonIso()",float,doc="PF absolute isolation dR=0.3, photon component (uncorrected)"),
0233 pfChargedIsoPFPV = Var("chargedHadronPFPVIso()",float,doc="PF absolute isolation dR=0.3, charged component (PF PV only)"),
0234 pfChargedIsoWorstVtx = Var("chargedHadronWorstVtxIso()",float,doc="PF absolute isolation dR=0.3, charged component (Vertex with largest isolation)"),
0235 pfRelIso03_chg_quadratic = Var("userFloat('PFIsoChgQuadratic')/pt",float,doc="PF relative isolation dR=0.3, charged hadron component (with quadraticEA*rho*rho + linearEA*rho Winter22V1 corrections)"),
0236 pfRelIso03_all_quadratic = Var("userFloat('PFIsoAllQuadratic')/pt",float,doc="PF relative isolation dR=0.3, total (with quadraticEA*rho*rho + linearEA*rho Winter22V1 corrections)"),
0237 pfRelIso03_chg_Fall17V2 = Var("userFloat('PFIsoChgFall17V2')/pt",float,doc="PF relative isolation dR=0.3, charged component (with Fall17V2 rho*EA PU corrections)"),
0238 pfRelIso03_all_Fall17V2 = Var("userFloat('PFIsoAllFall17V2')/pt",float,doc="PF relative isolation dR=0.3, total (with Fall17V2 rho*EA PU corrections)"),
0239 hoe = Var("hadronicOverEm()",float,doc="H over E",precision=8),
0240 hoe_PUcorr = Var("userFloat('HoverEQuadratic')",float,doc="PU corrected H/E (cone-based with quadraticEA*rho*rho + linearEA*rho Winter22V1 corrections)",precision=8),
0241 isScEtaEB = Var("abs(superCluster().eta()) < 1.4442",bool,doc="is supercluster eta within barrel acceptance"),
0242 isScEtaEE = Var("abs(superCluster().eta()) > 1.566 && abs(superCluster().eta()) < 2.5",bool,doc="is supercluster eta within endcap acceptance"),
0243 seedGain = Var("userInt('seedGain')","uint8",doc="Gain of the seed crystal"),
0244 seediEtaOriX = Var("superCluster().seedCrysIEtaOrIx","int8",doc="iEta or iX of seed crystal. iEta is barrel-only, iX is endcap-only. iEta runs from -85 to +85, with no crystal at iEta=0. iX runs from 1 to 100."),
0245 seediPhiOriY = Var("superCluster().seedCrysIPhiOrIy",int,doc="iPhi or iY of seed crystal. iPhi is barrel-only, iY is endcap-only. iPhi runs from 1 to 360. iY runs from 1 to 100."),
0246
0247 x_calo = Var("superCluster().seed().position().x()",float,doc="photon supercluster position on calorimeter, x coordinate (cm)",precision=10),
0248 y_calo = Var("superCluster().seed().position().y()",float,doc="photon supercluster position on calorimeter, y coordinate (cm)",precision=10),
0249 z_calo = Var("superCluster().seed().position().z()",float,doc="photon supercluster position on calorimeter, z coordinate (cm)",precision=10),
0250
0251 esEffSigmaRR = Var("full5x5_showerShapeVariables().effSigmaRR()", float, doc="preshower sigmaRR"),
0252 esEnergyOverRawE = Var("superCluster().preshowerEnergy()/superCluster().rawEnergy()", float, doc="ratio of preshower energy to raw supercluster energy"),
0253 haloTaggerMVAVal = Var("haloTaggerMVAVal()",float,doc="Value of MVA based BDT based beam halo tagger in the Ecal endcap (valid for pT > 200 GeV)",precision=8),
0254 )
0255 )
0256
0257
0258 run2_egamma.toModify(photonTable.variables,
0259 cutBased = photonTable.variables.cutBased_Fall17V2,
0260 cutBased_Fall17V2 = None,
0261 vidNestedWPBitmap = photonTable.variables.vidNestedWPBitmap_Fall17V2,
0262 vidNestedWPBitmap_Fall17V2 = None,
0263 mvaID = photonTable.variables.mvaID_Fall17V2,
0264 mvaID_Fall17V2 = None,
0265 mvaID_WP90 = photonTable.variables.mvaID_Fall17V2_WP90,
0266 mvaID_Fall17V2_WP90 = None,
0267 mvaID_WP80 = photonTable.variables.mvaID_Fall17V2_WP80,
0268 mvaID_Fall17V2_WP80 = None,
0269 pfRelIso03_chg = photonTable.variables.pfRelIso03_chg_Fall17V2,
0270 pfRelIso03_chg_Fall17V2 = None,
0271 pfRelIso03_all = photonTable.variables.pfRelIso03_all_Fall17V2,
0272 pfRelIso03_all_Fall17V2 = None,
0273 pfRelIso03_chg_quadratic=None,
0274 pfRelIso03_all_quadratic=None,
0275 hoe_PUcorr=None)
0276
0277
0278 run2_egamma.toModify(
0279 photonTable.variables,
0280 pt = Var("pt*userFloat('ecalEnergyPostCorrNew')/userFloat('ecalEnergyPreCorrNew')", float, precision=-1, doc="p_{T}"),
0281 energyErr = Var("userFloat('ecalEnergyErrPostCorrNew')",float,doc="energy error of the cluster from regression",precision=6),
0282 eCorr = Var("userFloat('ecalEnergyPostCorrNew')/userFloat('ecalEnergyPreCorrNew')",float,doc="ratio of the calibrated energy/miniaod energy"),
0283 hoe = Var("hadTowOverEm()",float,doc="H over E (Run2)",precision=8),
0284 )
0285
0286 photonsMCMatchForTable = cms.EDProducer("MCMatcher",
0287 src = photonTable.src,
0288 matched = cms.InputTag("finalGenParticles"),
0289 mcPdgId = cms.vint32(11,22),
0290 checkCharge = cms.bool(False),
0291 mcStatus = cms.vint32(1),
0292 maxDeltaR = cms.double(0.3),
0293 maxDPtRel = cms.double(0.5),
0294 resolveAmbiguities = cms.bool(True),
0295 resolveByMatchQuality = cms.bool(True),
0296 )
0297
0298 photonMCTable = cms.EDProducer("CandMCMatchTableProducer",
0299 src = photonTable.src,
0300 mcMap = cms.InputTag("photonsMCMatchForTable"),
0301 objName = photonTable.name,
0302 objType = photonTable.name,
0303 branchName = cms.string("genPart"),
0304 docString = cms.string("MC matching to status==1 photons or electrons"),
0305 )
0306
0307
0308 run2_egamma.toModify(
0309 photonTable.variables,
0310 dEscaleUp=Var("userFloat('ecalEnergyPostCorrNew') - userFloat('energyScaleUpNew')", float, doc="ecal energy scale shifted 1 sigma up (adding gain/stat/syst in quadrature)", precision=8),
0311 dEscaleDown=Var("userFloat('ecalEnergyPostCorrNew') - userFloat('energyScaleDownNew')", float, doc="ecal energy scale shifted 1 sigma down (adding gain/stat/syst in quadrature)", precision=8),
0312 dEsigmaUp=Var("userFloat('ecalEnergyPostCorrNew') - userFloat('energySigmaUpNew')", float, doc="ecal energy smearing value shifted 1 sigma up", precision=8),
0313 dEsigmaDown=Var("userFloat('ecalEnergyPostCorrNew') - userFloat('energySigmaDownNew')", float, doc="ecal energy smearing value shifted 1 sigma up", precision=8),
0314 )
0315
0316
0317 photonTask = cms.Task(bitmapVIDForPho, bitmapVIDForPhoRun2, isoForPho, hOverEForPho, isoForPhoFall17V2, seedGainPho, slimmedPhotonsWithUserData, finalPhotons)
0318
0319 photonTablesTask = cms.Task(photonTable)
0320 photonMCTask = cms.Task(photonsMCMatchForTable, photonMCTable)
0321
0322 _photonTask_Run2 = photonTask.copy()
0323 _photonTask_Run2.remove(bitmapVIDForPho)
0324 _photonTask_Run2.remove(isoForPho)
0325 _photonTask_Run2.remove(hOverEForPho)
0326 _photonTask_Run2.add(calibratedPatPhotonsNano)
0327 run2_egamma.toReplaceWith(photonTask, _photonTask_Run2)