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 
0016 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))"
0017 
0018 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)"
0019 
0020 
0021 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))"
0022 
0023 from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
0024 
0025 singleTopMuonMediumDQM_miniAOD = DQMEDAnalyzer('SingleTopTChannelLeptonDQM_miniAOD',
0026   setup = cms.PSet(
0027 
0028     directory = cms.string("Physics/Top/SingleTopMuonMediumDQM_miniAOD/"),
0029     sources = cms.PSet(
0030       muons = cms.InputTag("slimmedMuons"),
0031       elecs = cms.InputTag("slimmedElectrons"),
0032       jets  = cms.InputTag("slimmedJets"),
0033       mets  = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
0034       pvs   = cms.InputTag("offlineSlimmedPrimaryVertices")
0035     ),
0036 
0037     monitoring = cms.PSet(
0038       verbosity = cms.string("DEBUG")
0039     ),
0040 
0041     pvExtras = cms.PSet(
0042       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0043     ),
0044     elecExtras = cms.PSet(
0045       select = cms.string(tightElecCut + "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
0046       rho = cms.InputTag("fixedGridRhoFastjetAll"),
0047       #isolation = cms.string(ElelooseIsoCut),
0048     ),
0049                      
0050     muonExtras = cms.PSet(
0051       select    = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
0052       isolation = cms.string(looseIsoCut)
0053     ),
0054     jetExtras = cms.PSet(
0055 
0056       select = cms.string("pt>30 & abs(eta)<2.4"), # & neutralEmEnergyFraction >0.01 & chargedEmEnergyFraction>0.01"),
0057       jetBTaggers  = cms.PSet(
0058         trackCountingEff = cms.PSet(
0059           label = cms.InputTag("trackCountingHighEffBJetTags" ),
0060           workingPoint = cms.double(1.25)
0061         ),
0062         trackCountingPur = cms.PSet(
0063          label = cms.InputTag("trackCountingHighPurBJetTags" ),
0064           workingPoint = cms.double(3.41)
0065         ),
0066         secondaryVertex  = cms.PSet(
0067           label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
0068           workingPoint = cms.double(2.05)
0069         ),
0070         combinedSecondaryVertex  = cms.PSet(
0071           label = cms.InputTag("combinedSecondaryVertexBJetTags"),
0072           workingPoint = cms.double(0.898)
0073         )
0074      ),
0075    ),
0076    massExtras = cms.PSet(
0077      lowerEdge = cms.double( 70.),
0078      upperEdge = cms.double(110.)
0079    ),
0080   ),
0081   preselection = cms.PSet(
0082     vertex = cms.PSet(
0083       src    = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
0084       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0085     )
0086   ),
0087 
0088   selection = cms.VPSet(
0089     cms.PSet(
0090       label  = cms.string("muons:step0"),
0091       src    = cms.InputTag("slimmedMuons"),
0092       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
0093       min    = cms.int32(1),
0094       #max    = cms.int32(1),
0095     ),
0096     cms.PSet(
0097       label  = cms.string("jets:step1"),
0098       src    = cms.InputTag("slimmedJets"),
0099       select = cms.string("pt>30 & abs(eta)<2.4 "),
0100       min = cms.int32(2),
0101     ),
0102     cms.PSet(
0103       label  = cms.string("met:step2"),
0104       src    = cms.InputTag("slimmedMETs"),
0105       select = cms.string("pt>30"),
0106       #min = cms.int32(2),
0107     ),
0108   )
0109 )
0110 
0111 singleTopElectronMediumDQM_miniAOD = DQMEDAnalyzer('SingleTopTChannelLeptonDQM_miniAOD',
0112 
0113   setup = cms.PSet(
0114  
0115     directory = cms.string("Physics/Top/SingleTopElectronMediumDQM_miniAOD/"),
0116 
0117     sources = cms.PSet(
0118       muons = cms.InputTag("slimmedMuons"),
0119       elecs = cms.InputTag("slimmedElectrons"),
0120       jets  = cms.InputTag("slimmedJets"),
0121       mets  = cms.VInputTag("slimmedMETs", "slimmedMETsNoHF", "slimmedMETsPuppi"),
0122       pvs   = cms.InputTag("offlineSlimmedPrimaryVertices")
0123 
0124     ),
0125 
0126     monitoring = cms.PSet(
0127       verbosity = cms.string("DEBUG")
0128     ),
0129 
0130     pvExtras = cms.PSet(
0131       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0132     ),
0133     elecExtras = cms.PSet(
0134       select = cms.string(tightElecCut + " && pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
0135       #select     = cms.string(looseElecCut+ "&& pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660)"),
0136       rho = cms.InputTag("fixedGridRhoFastjetAll"),
0137       #isolation  = cms.string(ElelooseIsoCut),
0138     ),
0139     muonExtras = cms.PSet(
0140       select    = cms.string(tightMuonCut + " && pt>20 & abs(eta)<2.4 && " + looseIsoCut),
0141       isolation = cms.string(looseIsoCut)
0142     ),
0143 
0144     jetExtras = cms.PSet(
0145 
0146       select = cms.string("pt>30 & abs(eta)<2.4"),
0147       jetBTaggers  = cms.PSet(
0148         trackCountingEff = cms.PSet(
0149           label = cms.InputTag("trackCountingHighEffBJetTags" ),
0150           workingPoint = cms.double(1.25)
0151         ),
0152         trackCountingPur = cms.PSet(
0153           label = cms.InputTag("trackCountingHighPurBJetTags" ),
0154           workingPoint = cms.double(3.41)
0155         ),
0156         secondaryVertex  = cms.PSet(
0157           label = cms.InputTag("simpleSecondaryVertexHighEffBJetTags"),
0158           workingPoint = cms.double(2.05)
0159         ),
0160         combinedSecondaryVertex  = cms.PSet(
0161           label = cms.InputTag("combinedSecondaryVertexBJetTags"),
0162           workingPoint = cms.double(0.898)
0163         ),
0164       ),
0165     ),
0166     massExtras = cms.PSet(
0167       lowerEdge = cms.double( 70.),
0168       upperEdge = cms.double(110.)
0169     ),
0170                
0171   ),
0172 
0173   preselection = cms.PSet(
0174     vertex = cms.PSet(
0175       src    = cms.InputTag("offlineSlimmedPrimaryVertices"),#,
0176       select = cms.string("abs(z) < 24. & position.rho < 2. & ndof > 4 & !isFake")
0177     )
0178   ),
0179 
0180   selection = cms.VPSet(
0181     cms.PSet(
0182       label  = cms.string("elecs:step0"),
0183       src    = cms.InputTag("slimmedElectrons"),
0184       select = cms.string("pt>20 & abs(eta)<2.5 & (abs(superCluster.eta) <= 1.4442 || abs(superCluster.eta) >= 1.5660) &&" + tightElecCut),
0185       min    = cms.int32(1),
0186       max    = cms.int32(1),
0187     ),
0188     cms.PSet(
0189       label  = cms.string("jets:step1"),
0190       src    = cms.InputTag("slimmedJets"),
0191       select = cms.string("pt>30 & abs(eta)<2.4 "),
0192       min = cms.int32(2),
0193     ),
0194     cms.PSet(
0195       label  = cms.string("met:step2"),
0196       src    = cms.InputTag("slimmedMETs"),
0197       select = cms.string("pt>30"),
0198       #min = cms.int32(2),
0199     ),
0200   )
0201 )