diff --git a/PhysicsTools/NanoAOD/python/lowPtElectrons_cff.py b/PhysicsTools/NanoAOD/python/lowPtElectrons_cff.py index 571bc72fcfd0a..ddec293ac8656 100644 --- a/PhysicsTools/NanoAOD/python/lowPtElectrons_cff.py +++ b/PhysicsTools/NanoAOD/python/lowPtElectrons_cff.py @@ -16,10 +16,22 @@ ) ) +import PhysicsTools.PatAlgos.producersLayer1.electronProducer_cfi +updatedLowPtElectrons = cms.EDProducer( + "PATElectronUpdater", + src = cms.InputTag("modifiedLowPtElectrons"), + vertices = cms.InputTag("offlineSlimmedPrimaryVertices"), + computeMiniIso = cms.bool(True), + fixDxySign = cms.bool(False), + pfCandsForMiniIso = cms.InputTag("packedPFCandidates"), + miniIsoParamsB = PhysicsTools.PatAlgos.producersLayer1.electronProducer_cfi.patElectrons.miniIsoParamsB, + miniIsoParamsE = PhysicsTools.PatAlgos.producersLayer1.electronProducer_cfi.patElectrons.miniIsoParamsE, +) + from RecoEgamma.EgammaElectronProducers.lowPtGsfElectronID_cfi import lowPtGsfElectronID lowPtPATElectronID = lowPtGsfElectronID.clone( usePAT = True, - electrons = "modifiedLowPtElectrons", + electrons = "updatedLowPtElectrons", unbiased = "", ModelWeights = [ 'RecoEgamma/ElectronIdentification/data/LowPtElectrons/LowPtElectrons_ID_2020Nov28.root', @@ -28,8 +40,8 @@ isoForLowPtEle = cms.EDProducer( "EleIsoValueMapProducer", - src = cms.InputTag("modifiedLowPtElectrons"), - relative = cms.bool(False), + src = cms.InputTag("updatedLowPtElectrons"), + relative = cms.bool(True), rho_MiniIso = cms.InputTag("fixedGridRhoFastjetAll"), rho_PFIso = cms.InputTag("fixedGridRhoFastjetAll"), EAFile_MiniIso = cms.FileInPath("RecoEgamma/ElectronIdentification/data/Fall17/effAreaElectrons_cone03_pfNeuHadronsAndPhotons_94X.txt"), @@ -38,7 +50,7 @@ updatedLowPtElectronsWithUserData = cms.EDProducer( "PATElectronUserDataEmbedder", - src = cms.InputTag("modifiedLowPtElectrons"), + src = cms.InputTag("updatedLowPtElectrons"), userFloats = cms.PSet( ID = cms.InputTag("lowPtPATElectronID"), miniIsoChg = cms.InputTag("isoForLowPtEle:miniIsoChg"), @@ -76,9 +88,9 @@ unbiased = Var("electronID('unbiased')",float,doc="ElectronSeed, pT- and dxy- agnostic BDT (raw) score"), ptbiased = Var("electronID('ptbiased')",float,doc="ElectronSeed, pT- and dxy- dependent BDT (raw) score"), # Isolation - miniPFRelIso_chg = Var("userFloat('miniIsoChg')/pt",float, + miniPFRelIso_chg = Var("userFloat('miniIsoChg')",float, doc="mini PF relative isolation, charged component"), - miniPFRelIso_all = Var("userFloat('miniIsoAll')/pt",float, + miniPFRelIso_all = Var("userFloat('miniIsoAll')",float, doc="mini PF relative isolation, total (with scaled rho*EA PU corrections)"), # Conversions convVeto = Var("passConversionVeto()",bool,doc="pass conversion veto"), @@ -100,8 +112,6 @@ dxyErr = Var("edB('PV2D')",float,doc="dxy uncertainty, in cm",precision=6), dz = Var("dB('PVDZ')",float,doc="dz (with sign) wrt first PV, in cm",precision=10), dzErr = Var("abs(edB('PVDZ'))",float,doc="dz uncertainty, in cm",precision=6), - ip3d = Var("abs(dB('PV3D'))",float,doc="3D impact parameter wrt first PV, in cm",precision=10), - sip3d = Var("abs(dB('PV3D')/edB('PV3D'))",float,doc="3D impact parameter significance wrt first PV, in cm",precision=10), # Cross-referencing #jetIdx #photonIdx @@ -175,6 +185,7 @@ ################################################################################ lowPtElectronSequence = cms.Sequence(modifiedLowPtElectrons + +updatedLowPtElectrons +lowPtPATElectronID +isoForLowPtEle +updatedLowPtElectronsWithUserData diff --git a/PhysicsTools/NanoAOD/python/nanoDQM_cfi.py b/PhysicsTools/NanoAOD/python/nanoDQM_cfi.py index 19d6f305b16f2..4ac38307e74f2 100644 --- a/PhysicsTools/NanoAOD/python/nanoDQM_cfi.py +++ b/PhysicsTools/NanoAOD/python/nanoDQM_cfi.py @@ -133,8 +133,8 @@ Plot1D('unbiased', 'unbiased', 40, -10., 10., 'ElectronSeed, pT- and dxy- agnostic BDT (raw) score'), Plot1D('ptbiased', 'ptbiased', 40, -10., 10., 'ElectronSeed, pT- and dxy- dependent BDT (raw) score'), # Isolation - Plot1D('miniPFRelIso_chg', 'miniPFRelIso_chg', 20, 0, 2, 'mini PF relative isolation, charged component'), - Plot1D('miniPFRelIso_all', 'miniPFRelIso_all', 20, 0, 2, 'mini PF relative isolation, total (with scaled rho*EA PU corrections)'), + Plot1D('miniPFRelIso_chg', 'miniPFRelIso_chg', 20, 0., 1., 'mini PF relative isolation, charged component'), + Plot1D('miniPFRelIso_all', 'miniPFRelIso_all', 20, 0., 1., 'mini PF relative isolation, total (with scaled rho*EA PU corrections)'), # Conversions Plot1D('convVeto', 'convVeto', 2, -0.5, 1.5, 'pass conversion veto'), Plot1D('convWP', 'convWP', 8, -0.5, 7.5, 'conversion flag bit map: 1=Veto, 2=Loose, 3=Tight'), @@ -154,8 +154,6 @@ Plot1D('dz', 'dz', 20, -0.3, 0.3, 'dz (with sign) wrt first PV, in cm'), Plot1D('dxyErr', 'dxyErr', 20, 0., 0.2, 'dxy uncertainty, in cm'), Plot1D('dzErr', 'dzErr', 20, 0., 0.2, 'dz uncertainty, in cm'), - Plot1D('ip3d', 'ip3d', 20, 0., 0.2, '3D impact parameter wrt first PV, in cm'), - Plot1D('sip3d', 'sip3d', 20, 0., 20., '3D impact parameter significance wrt first PV, in cm'), ), ),