File indexing completed on 2024-04-25 02:14:04
0001 import FWCore.ParameterSet.Config as cms
0002 from PhysicsTools.NanoAOD.common_cff import *
0003 from PhysicsTools.NanoAOD.simplePATIsolatedTrackFlatTableProducer_cfi import simplePATIsolatedTrackFlatTableProducer
0004
0005 finalIsolatedTracks = cms.EDProducer("IsolatedTrackCleaner",
0006 tracks = cms.InputTag("isolatedTracks"),
0007 cut = cms.string("((pt>5 && (abs(pdgId) == 11 || abs(pdgId) == 13)) || pt > 10) && (abs(pdgId) < 15 || abs(eta) < 2.5) && ((abs(dxy) < 0.2 && abs(dz) < 0.1) || pt>15) && ((pfIsolationDR03().chargedHadronIso < 5 && pt < 25) || pfIsolationDR03().chargedHadronIso/pt < 0.2)"),
0008 finalLeptons = cms.VInputTag(
0009 cms.InputTag("finalElectrons"),
0010 cms.InputTag("finalLooseMuons"),
0011 ),
0012 )
0013
0014 isoForIsoTk = cms.EDProducer("IsoTrackIsoValueMapProducer",
0015 src = cms.InputTag("finalIsolatedTracks"),
0016 relative = cms.bool(True),
0017 rho_MiniIso = cms.InputTag("fixedGridRhoFastjetAll"),
0018 EAFile_MiniIso = cms.FileInPath("PhysicsTools/NanoAOD/data/effAreaMuons_cone03_pfNeuHadronsAndPhotons_80X.txt"),
0019 )
0020
0021 isFromLostTrackForIsoTk = cms.EDProducer("IsFromLostTrackMapProducer",
0022 srcIsoTracks = cms.InputTag("finalIsolatedTracks"),
0023 lostTracks = cms.InputTag("lostTracks"),
0024 )
0025
0026 isoTrackTable = simplePATIsolatedTrackFlatTableProducer.clone(
0027 src = cms.InputTag("finalIsolatedTracks"),
0028 name = cms.string("IsoTrack"),
0029 doc = cms.string("isolated tracks after basic selection (" + finalIsolatedTracks.cut.value() + ") and lepton veto"),
0030 variables = cms.PSet(P3Vars,
0031 dz = Var("dz",float,doc="dz (with sign) wrt first PV, in cm",precision=10),
0032 dxy = Var("dxy",float,doc="dxy (with sign) wrt first PV, in cm",precision=10),
0033 pfRelIso03_chg = Var("pfIsolationDR03().chargedHadronIso/pt",float,doc="PF relative isolation dR=0.3, charged component",precision=10),
0034 pfRelIso03_all = Var("(pfIsolationDR03().chargedHadronIso + max(pfIsolationDR03().neutralHadronIso + pfIsolationDR03().photonIso - pfIsolationDR03().puChargedHadronIso/2,0.0))/pt",float,doc="PF relative isolation dR=0.3, total (deltaBeta corrections)",precision=10),
0035 isPFcand = Var("packedCandRef().isNonnull()",bool,doc="if isolated track is a PF candidate"),
0036 fromPV = Var("fromPV", "int16", doc="isolated track comes from PV"),
0037 pdgId = Var("pdgId",int,doc="PDG id of PF cand"),
0038 isHighPurityTrack = Var("isHighPurityTrack",bool,doc="track is high purity"),
0039 charge = Var("charge", "int16", doc="electric charge"),
0040 ),
0041 externalVariables = cms.PSet(
0042 miniPFRelIso_chg = ExtVar("isoForIsoTk:miniIsoChg",float,doc="mini PF relative isolation, charged component",precision=10),
0043 miniPFRelIso_all = ExtVar("isoForIsoTk:miniIsoAll",float,doc="mini PF relative isolation, total (with scaled rho*EA PU corrections)",precision=10),
0044 isFromLostTrack = ExtVar("isFromLostTrackForIsoTk:isFromLostTrack",bool,doc="if isolated track comes from a lost track"),
0045 ),
0046 )
0047
0048 isoTrackTask = cms.Task(finalIsolatedTracks,isoForIsoTk,isFromLostTrackForIsoTk)
0049 isoTrackTablesTask = cms.Task(isoTrackTable)
0050