Physics/python/topSingleLeptonDQM_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 #Primary vertex selection
0004 PVCut = "abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake"
0005
0006 #Jet selection
0007 looseJetCut = "(chargedHadronEnergyFraction()>0 && chargedMultiplicity()>0 && chargedEmEnergyFraction()<0.99 && neutralHadronEnergyFraction()<0.99 && neutralEmEnergyFraction()<0.99 && (chargedMultiplicity()+neutralMultiplicity())>1) && abs(eta)<=2.4 "
0008
0009 tightJetCut = "(chargedHadronEnergyFraction()>0 && chargedMultiplicity()>0 && chargedEmEnergyFraction()<0.99 && neutralHadronEnergyFraction()<0.90 && neutralEmEnergyFraction()<0.90 && (chargedMultiplicity()+neutralMultiplicity())>1) && abs(eta)<=2.4 "
0010
0011 #Loose muon selection
0012 looseMuonCut  = "((isGlobalMuon || isTrackerMuon) && isPFMuon)"
0013 looseIsoCut   = "((pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.25)"
0014
0015 #Medium muon selection. Also requires either good global muon or tight segment compatibility
0016 mediumMuonCut = looseMuonCut + " innerTrack.validFraction > 0.8"
0017
0018 #Tight muon selection. Lacks distance to primary vertex variables, dz<0.5, dxy < 0.2. Now done at .cc
0019 tightMuonCut  = "isGlobalMuon && isPFMuon && globalTrack.normalizedChi2 < 10. && globalTrack.hitPattern.numberOfValidMuonHits > 0 && " + \
0020                "numberOfMatchedStations > 1 && innerTrack.hitPattern.numberOfValidPixelHits > 0 && innerTrack.hitPattern.trackerLayersWithMeasurement > 5 "
0021 tightIsoCut   = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.15"
0022
0023 #Electron selections
0024 looseEleCut = "((full5x5_sigmaIetaIeta() < 0.011 && superCluster().isNonnull() && superCluster().seed().isNonnull() && (deltaEtaSuperClusterTrackAtVtx() - superCluster().eta() + superCluster().seed().eta()) < 0.00477 && abs(deltaPhiSuperClusterTrackAtVtx()) < 0.222 && hadronicOverEm() < 0.298 && abs(1.0 - eSuperClusterOverP())*1.0/ecalEnergy() < 0.241 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster().eta()) < 1.479) || (full5x5_sigmaIetaIeta() < 0.0314 && superCluster().isNonnull() && superCluster().seed().isNonnull() && (deltaEtaSuperClusterTrackAtVtx() - superCluster().eta() + superCluster().seed().eta()) < 0.00868 && abs(deltaPhiSuperClusterTrackAtVtx()) < 0.213 && hadronicOverEm() < 0.101  && abs(1.0 - eSuperClusterOverP())*1.0/ecalEnergy() < 0.14 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster().eta()) > 1.479))"
0025
0026 tightEleCut = "((full5x5_sigmaIetaIeta() < 0.00998 && superCluster().isNonnull() && superCluster().seed().isNonnull() && (deltaEtaSuperClusterTrackAtVtx() - superCluster().eta() + superCluster().seed().eta()) < 0.00308  && abs(deltaPhiSuperClusterTrackAtVtx()) < 0.0816 && hadronicOverEm() < 0.0414 && abs(1.0 - eSuperClusterOverP())*1.0/ecalEnergy() < 0.0129 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster().eta()) < 1.479) ||  (full5x5_sigmaIetaIeta() < 0.0292 && superCluster().isNonnull() && superCluster().seed().isNonnull() && (deltaEtaSuperClusterTrackAtVtx() - superCluster().eta() + superCluster().seed().eta()) < 0.00605 && abs(deltaPhiSuperClusterTrackAtVtx()) < 0.0394  && hadronicOverEm() < 0.0641  && abs(1.0 - eSuperClusterOverP())*1.0/ecalEnergy() < 0.0129 && gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 1 && abs(superCluster().eta()) > 1.479))"
0027
0028
0029 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0030
0031 topSingleMuonMediumDQM = DQMEDAnalyzer('TopSingleLeptonDQM',
0032   ## ------------------------------------------------------
0033   ## SETUP
0034   ##
0035   ## configuration of the MonitoringEnsemble(s)
0036   ## [mandatory] : optional PSets may be omitted
0037   ##
0038   setup = cms.PSet(
0039     directory = cms.string("Physics/Top/TopSingleMuonMediumDQM/"),
0040     sources = cms.PSet(
0041       muons = cms.InputTag("muons"),
0042       elecs = cms.InputTag("gedGsfElectrons"),
0043       jets  = cms.InputTag("ak4PFJetsCHS"),
0044       mets  = cms.VInputTag("pfMet"),
0045       pvs   = cms.InputTag("offlinePrimaryVertices")
0046
0047     ),
0048     monitoring = cms.PSet(
0049       verbosity = cms.string("DEBUG")
0050     ),
0051     pvExtras = cms.PSet(
0052       select = cms.string(PVCut)
0053     ),
0054     elecExtras = cms.PSet(
0055       select     = cms.string(tightEleCut + "& pt>20 & abs(eta)<2.5 & (abs(superCluster().eta()) <= 1.4442 || abs(superCluster().eta()) >= 1.5660)"),
0056             rho = cms.InputTag("fixedGridRhoFastjetAll"),
0057     ),
0058     muonExtras = cms.PSet(
0059       select    = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4"),
0060       #select    = cms.string(looseMuonCut + " && pt>20 & abs(eta)<2.4"),
0061       isolation = cms.string(looseIsoCut)
0062     ),
0063     jetExtras = cms.PSet(
0064       jetCorrector = cms.InputTag("dqmAk4PFCHSL1FastL2L3Corrector"), #Use pak4PFCHSL1FastL2L3Residual for data!!!
0065       select = cms.string("pt>30 & abs(eta)<2.4"),
0066       jetBTaggers = cms.PSet(
0067                 cvsVertex = cms.PSet(
0068             label = cms.InputTag("pfDeepCSVJetTags:probb"),
0069               workingPoint = cms.double(0.4168)
0070               # deepCSV Medium from https://btv-wiki.docs.cern.ch/ScaleFactors/UL2018/
0071         )
0072       ),
0073     ),
0074     massExtras = cms.PSet(
0075       lowerEdge = cms.double( 70.),
0076       upperEdge = cms.double(110.)
0077     ),
0078   ),
0079
0080   preselection = cms.PSet(
0081     vertex = cms.PSet(
0082       src    = cms.InputTag("offlinePrimaryVertices"),
0083       select = cms.string(PVCut)
0084     )
0085   ),
0086
0087   selection = cms.VPSet(
0088     cms.PSet(
0089       label  = cms.string("muons:step0"),
0090       src    = cms.InputTag("muons"),
0091       select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4"),
0092       min    = cms.int32(1),
0093     ),
0094     cms.PSet(
0095       label  = cms.string("jets/pf:step1"),
0096       src    = cms.InputTag("ak4PFJetsCHS"),
0097       select = cms.string("pt>30 & abs(eta)<2.4"),
0098       min = cms.int32(4),
0099     ),
0100     cms.PSet(
0101       label  = cms.string("met:step2"),
0102       src    = cms.InputTag("pfMet"),
0103       select = cms.string("pt>30"),
0104     ),
0105   )
0106 )
0107
0108
0109 topSingleElectronMediumDQM = DQMEDAnalyzer('TopSingleLeptonDQM',
0110   ## ------------------------------------------------------
0111   ## SETUP
0112   ##
0113   ## configuration of the MonitoringEnsemble(s)
0114   ## [mandatory] : optional PSets may be omitted
0115   ##
0116   setup = cms.PSet(
0117     directory = cms.string("Physics/Top/TopSingleElectronMediumDQM/"),
0118     sources = cms.PSet(
0119       muons = cms.InputTag("muons"),
0120       elecs = cms.InputTag("gedGsfElectrons"),
0121       jets  = cms.InputTag("ak4PFJetsCHS"),
0122       mets  = cms.VInputTag("pfMet"),
0123       pvs   = cms.InputTag("offlinePrimaryVertices")
0124     ),
0125     monitoring = cms.PSet(
0126       verbosity = cms.string("DEBUG")
0127     ),
0128     pvExtras = cms.PSet(
0129       select = cms.string(PVCut)
0130     ),
0131     elecExtras = cms.PSet(
0132       select   = cms.string(tightEleCut + "& pt>20 & abs(eta)<2.5 & (abs(superCluster().eta()) <= 1.4442 || abs(superCluster().eta()) >= 1.5660)"),
0133             rho = cms.InputTag("fixedGridRhoFastjetAll"),
0134     ),
0135     muonExtras = cms.PSet(
0136       select     = cms.string(tightMuonCut + " & pt>20 & abs(eta)<2.4"),
0137       isolation  = cms.string(looseIsoCut),
0138     ),
0139     jetExtras = cms.PSet(
0140             jetCorrector = cms.InputTag("dqmAk4PFCHSL1FastL2L3Corrector"), #Use pak4PFCHSL1FastL2L3Residual for data!!!
0141       select = cms.string("pt>30 & abs(eta)<2.4"),
0142       jetBTaggers = cms.PSet(
0143               cvsVertex = cms.PSet(
0144           label = cms.InputTag("pfDeepCSVJetTags:probb"),
0145             workingPoint = cms.double(0.4168)
0146               # deepCSV Medium from https://btv-wiki.docs.cern.ch/ScaleFactors/UL2018/
0147         )
0148       )
0149     ),
0150     massExtras = cms.PSet(
0151       lowerEdge = cms.double( 70.),
0152       upperEdge = cms.double(110.)
0153     ),
0154   ),
0155   preselection = cms.PSet(
0156     vertex = cms.PSet(
0157       src    = cms.InputTag("offlinePrimaryVertices"),
0158       select = cms.string(PVCut)
0159     )
0160   ),
0161   selection = cms.VPSet(
0162     cms.PSet(
0163       label = cms.string("elecs:step0"),
0164       src   = cms.InputTag("gedGsfElectrons"),
0165       select = cms.string("pt>20 & abs(eta)<2.5 & (abs(superCluster().eta()) <= 1.4442 || abs(superCluster().eta()) >= 1.5660) &&" + tightEleCut),
0166  #    select = cms.string("pt>30 & abs(eta)<2.5 & abs(gsfTrack.d0)<0.02 & gsfTrack.hitPattern().numberOfLostHits('MISSING_INNER_HITS') <= 0 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) & " + EletightIsoCut),
0167       min = cms.int32(1),
0168     ),
0169     cms.PSet(
0170       label = cms.string("jets/pf:step1"),
0171       src   = cms.InputTag("ak4PFJetsCHS"),
0172       select = cms.string("pt>30 & abs(eta)<2.4"),
0173       min = cms.int32(4),
0174     ),
0175     cms.PSet(
0176       label  = cms.string("met:step2"),
0177       src    = cms.InputTag("pfMet"),
0178       select = cms.string("pt>30"),
0179     ),
0180   )
0181 )
0182