Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:10:52

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 looseMuonCut = " (isGlobalMuon || isTrackerMuon) && isPFMuon"
0004 looseIsoCut  = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.25"
0005 
0006 tightMuonCut = " isGlobalMuon && isPFMuon && globalTrack.normalizedChi2 < 10. && globalTrack.hitPattern.numberOfValidMuonHits > 0 && " + \
0007                "numberOfMatchedStations > 1 && innerTrack.hitPattern.numberOfValidPixelHits > 0 && innerTrack.hitPattern.trackerLayersWithMeasurement > 5"
0008                # CB PV cut!
0009 tightIsoCut  = "(pfIsolationR04.sumChargedHadronPt + max(0., pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - 0.5 * pfIsolationR04.sumPUPt) ) / pt < 0.15"
0010 
0011 
0012 EletightIsoCut  = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.1"
0013 ElelooseIsoCut  = "(pfIsolationVariables.sumChargedHadronPt + max(0., pfIsolationVariables.sumNeutralHadronEt + pfIsolationVariables.sumPhotonEt - 0.5 * pfIsolationVariables.sumPUPt) ) / pt < 0.15"
0014 
0015 looseElecCut = "((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.numberOfHits('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.numberOfHits('MISSING_INNER_HITS') <= 1 && abs(superCluster.eta) > 1.479))"
0016 
0017 elecIPcut = "(abs(gsfTrack.d0)<0.05 & abs(gsfTrack.dz)<0.1 & abs(superCluster.eta) < 1.479)||(abs(gsfTrack.d0)<0.1 && abs(gsfTrack.dz)<0.2 && abs(superCluster.eta) > 1.479)"
0018 
0019 
0020 tightElecCut = "((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))"
0021 
0022 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0023 
0024 
0025 topSingleMuonMediumDQM_miniAOD = DQMEDAnalyzer('TopSingleLeptonDQM_miniAOD',
0026   setup = cms.PSet(
0027     directory = cms.string("Physics/Top/TopSingleMuonMediumDQM_miniAOD/"),
0028     sources = cms.PSet(
0029       muons = cms.InputTag("slimmedMuons"),
0030       elecs = cms.InputTag("slimmedElectrons"),
0031       jets  = cms.InputTag("slimmedJets"),
0032       mets  = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
0033       pvs   = cms.InputTag("offlineSlimmedPrimaryVertices")
0034 
0035     ),
0036     monitoring = cms.PSet(
0037       verbosity = cms.string("DEBUG")
0038     ),
0039     pvExtras = cms.PSet(
0040       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0041     ),
0042     elecExtras = cms.PSet(
0043       select = cms.string(tightElecCut + "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
0044       rho = cms.InputTag("fixedGridRhoFastjetAll"),
0045     #isolation = cms.string(ElelooseIsoCut),
0046     ),
0047                
0048     muonExtras = cms.PSet(
0049       select    = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
0050       isolation = cms.string(looseIsoCut)
0051     ),
0052     jetExtras = cms.PSet(                       
0053       select = cms.string("pt>30 & abs(eta)<2.4 "),
0054       jetBTaggers  = cms.PSet(
0055         trackCountingEff = cms.PSet(
0056           label = cms.InputTag("trackCountingHighEffBJetTags" ),
0057           workingPoint = cms.double(1.25)
0058         ),
0059         trackCountingPur = cms.PSet(
0060           label = cms.InputTag("trackCountingHighPurBJetTags" ),
0061           workingPoint = cms.double(3.00)
0062         ),
0063         secondaryVertex  = cms.PSet(
0064           label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
0065           workingPoint = cms.double(2.05)
0066         ),
0067         cvsVertex = cms.PSet(
0068           label = cms.InputTag("combinedSecondaryVertexBJetTags"),
0069           workingPoint = cms.double(0.898) 
0070         )
0071       ),                                                
0072     ),
0073     massExtras = cms.PSet(
0074       lowerEdge = cms.double( 70.),
0075       upperEdge = cms.double(110.)
0076     ),
0077   ),
0078   preselection = cms.PSet(
0079     vertex = cms.PSet(
0080       src    = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
0081       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0082     )
0083   ),
0084   selection = cms.VPSet(
0085     cms.PSet(
0086       label  = cms.string("muons:step0"),
0087       src    = cms.InputTag("slimmedMuons"),
0088       select = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut), #tightMuonCut +"&&"+ tightIsoCut + " && pt>20 & abs(eta)<2.1"), # CB what about iso? CD Added tightIso
0089       min    = cms.int32(1),
0090       #max    = cms.int32(1),
0091     ),
0092     cms.PSet(
0093       label  = cms.string("jets:step1"),
0094       src    = cms.InputTag("slimmedJets"),
0095       select = cms.string("pt>30 & abs(eta)<2.4 "),
0096       min = cms.int32(4),
0097     ), 
0098     cms.PSet(
0099       label  = cms.string("met:step2"),
0100       src    = cms.InputTag("slimmedMETs"),
0101       select = cms.string("pt>30"),
0102       #min = cms.int32(2),
0103     ), 
0104 #    cms.PSet(
0105 #      label  = cms.string("jets:step3"),
0106 #      src    = cms.InputTag("slimmedJets"),
0107 #      select = cms.string("pt>30 & abs(eta)<2.5 "),
0108 #      min = cms.int32(3),                                                
0109 #    ), 
0110 #    cms.PSet(
0111 #      label  = cms.string("jets:step4"),
0112 #      src    = cms.InputTag("slimmedJets"),
0113 #      select = cms.string("pt>30 & abs(eta)<2.5 "),
0114 #      min = cms.int32(4),                                                
0115 #    ),
0116   )
0117 )
0118 
0119 
0120 
0121 topSingleElectronMediumDQM_miniAOD = DQMEDAnalyzer('TopSingleLeptonDQM_miniAOD',
0122   setup = cms.PSet(
0123 
0124     directory = cms.string("Physics/Top/TopSingleElectronMediumDQM_miniAOD/"),
0125     sources = cms.PSet(
0126       muons = cms.InputTag("slimmedMuons"),
0127       elecs = cms.InputTag("slimmedElectrons"),
0128       jets  = cms.InputTag("slimmedJets"),
0129       mets  = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
0130       pvs   = cms.InputTag("offlineSlimmedPrimaryVertices")
0131 
0132     ),
0133     monitoring = cms.PSet(
0134       verbosity = cms.string("DEBUG")
0135     ),
0136     pvExtras = cms.PSet(
0137       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0138     ),
0139     elecExtras = cms.PSet(
0140       select = cms.string(tightElecCut + " && pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
0141                           #select     = cms.string(looseElecCut+ "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
0142       rho = cms.InputTag("fixedGridRhoFastjetAll"),
0143                           #isolation  = cms.string(ElelooseIsoCut),
0144     ),
0145     muonExtras = cms.PSet(
0146       select    = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
0147       isolation = cms.string(looseIsoCut)
0148     ),
0149     jetExtras = cms.PSet(
0150       select = cms.string("pt>30 & abs(eta)<2.4 "),
0151       jetBTaggers  = cms.PSet(
0152         trackCountingEff = cms.PSet(
0153           label = cms.InputTag("trackCountingHighEffBJetTags" ),
0154           workingPoint = cms.double(1.25)
0155         ),
0156         trackCountingPur = cms.PSet(
0157           label = cms.InputTag("trackCountingHighPurBJetTags" ),
0158           workingPoint = cms.double(3.00)
0159         ),
0160         secondaryVertex  = cms.PSet(
0161           label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
0162           workingPoint = cms.double(2.05)
0163         ),
0164         cvsVertex = cms.PSet(
0165           label = cms.InputTag("combinedSecondaryVertexBJetTags"),
0166           workingPoint = cms.double(0.898) 
0167         )
0168       ),
0169     ),
0170     massExtras = cms.PSet(
0171       lowerEdge = cms.double( 70.),
0172       upperEdge = cms.double(110.)
0173     ),
0174   ),
0175   preselection = cms.PSet(
0176     vertex = cms.PSet(
0177       src    = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
0178       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0179     )
0180   ),
0181   selection = cms.VPSet(
0182     cms.PSet(
0183       label  = cms.string("elecs:step0"),
0184       src    = cms.InputTag("slimmedElectrons"),
0185       select = cms.string("pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) &" + tightElecCut),
0186       min    = cms.int32(1),
0187       max    = cms.int32(1),
0188     ),
0189     cms.PSet(
0190       label  = cms.string("jets:step1"),
0191       src    = cms.InputTag("slimmedJets"),
0192       select = cms.string("pt>30 & abs(eta)<2.4 "),
0193       min = cms.int32(4),
0194     ),
0195     cms.PSet(
0196       label  = cms.string("met:step2"),
0197       src    = cms.InputTag("slimmedMETs"),
0198       select = cms.string("pt>30"),
0199       #min = cms.int32(2),
0200     ),
0201                         #cms.PSet(
0202                         #label  = cms.string("jets:step3"),
0203                         #src    = cms.InputTag("slimmedJets"),
0204                         #select = cms.string("pt>30 & abs(eta)<2.5 "),
0205                         #min = cms.int32(3),
0206                         #),
0207                         #cms.PSet(
0208                         #label  = cms.string("jets:step4"),
0209                         #src    = cms.InputTag("slimmedJets"),
0210                         #select = cms.string("pt>30 & abs(eta)<2.5 "),
0211                         #min = cms.int32(4),
0212                         #),
0213   )
0214 )