import FWCore.ParameterSet.Config as cms
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector
from PhysicsTools.SelectorUtils.pfJetIDSelector_cfi import pfJetIDSelector
PATElectronNtupleMaker = cms.EDFilter(
'PATElectronNtupleMaker',
hltTag = cms.InputTag("TriggerResults::REDIGI"),
MuonTag = cms.InputTag("selectedPatMuons"),
ElectronTag = cms.InputTag("selectedPatElectrons"),
PFElectronTag = cms.InputTag("patElectronsPFlowLoose"),
caloJetTag= cms.InputTag("selectedPatJets"),
JPTJetTag= cms.InputTag("selectedPatJetsAK5JPT"),
PFJetTag= cms.InputTag("selectedPatJetsAK5PF"),
PFlowJetTag= cms.InputTag("selectedPatJetsPFlow"),
caloMETTag= cms.InputTag("patMETs"),
tcMETTag= cms.InputTag("patMETsTC"),
PFMETTag= cms.InputTag("patMETsPF"),
PFlowMETTag = cms.InputTag("patMETsPFlow"),
jetIdLoose = jetIDSelector.clone(),
pfjetIdLoose = pfJetIDSelector.clone(),
ntupleFile = cms.string("NtupleMaker.root"),
inputType = cms.string("DATA")
)
## Maximal Number of Events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(100) )
from Analysis.JetAnalysis.multijetStudies_cfi import jetStudies as inputJetStudies
from Analysis.JetAnalysis.multijetStudies_cfi import pfJetStudies as inputPFJetStudies
from PhysicsTools.SelectorUtils.pfJetIDSelector_cfi import pfJetIDSelector as inputPFJetIDSelector
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector as inputJetIDSelector
process.TFileService = cms.Service("TFileService",
fileName = cms.string("multijetStudies.root")
)
process.multijetAna = cms.EDAnalyzer('EDMultijetAnalysis',
pfJetIDSelector = inputPFJetIDSelector.clone(),
jetIDSelector = inputJetIDSelector.clone(),
jetStudies = inputJetStudies.clone(ptMin = cms.double(50.0)),
pfJetStudies = inputPFJetStudies.clone(ptMin = cms.double(50.0)),
plotParameters = cms.PSet ( doTracks = cms.bool(False),
useMC = cms.bool(False)
)
)
process.p = cms.Path(
process.multijetAna
)
process.MessageLogger.cerr.FwkReport.reportEvery = 10000
import FWCore.ParameterSet.Config as cms
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector
from PhysicsTools.SelectorUtils.pfJetIDSelector_cfi import pfJetIDSelector
from Analysis.BoostedTopAnalysis.BoostedTopWTagParams_cfi import boostedTopWTagParams as inputBoostedTopWTagParams
MistagMakerParams = cms.PSet(
dijetSelectorParams=cms.PSet(
jetSrc=cms.InputTag('selectedPatJetsCA8PrunedCalo'),
pfJetSrc=cms.InputTag('goodPatJetsCA8PrunedPF'),
metSrc=cms.InputTag('patMETs'),
pfMetSrc=cms.InputTag('patMETsPF'),
useCalo=cms.bool(False),
minNJets=cms.uint32(2),
ptMin=cms.double(200),
etaMax=cms.double(2.4),
htMin=cms.double(0.0),
ht2Min=cms.double(0.0)),
jetIDParams=cms.PSet(jetIDSelector.clone()),
pfJetIDParams=cms.PSet(pfJetIDSelector.clone()),
boostedTopWTagParams=cms.PSet(inputBoostedTopWTagParams),
wMassMin=cms.double(60),
wMassMax=cms.double(100),
jetPtMin=cms.double(200))
'JEC/MCRUN2_74_V9::All_L1FastJet_AK4PFchs.txt',
'JEC/MCRUN2_74_V9::All_L2Relative_AK4PFchs.txt',
'JEC/MCRUN2_74_V9::All_L3Absolute_AK4PFchs.txt'
]
if corrMETonTheFly:
jecLevelsAK4 = [
'JEC/MCRUN2_74_V9::All_L1FastJet_AK4PF.txt',
'JEC/MCRUN2_74_V9::All_L2Relative_AK4PF.txt',
'JEC/MCRUN2_74_V9::All_L3Absolute_AK4PF.txt'
]
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector
process.goodSlimmedJets = cms.EDFilter(
"JetIDSelectionFunctorFilter",
filterParams=jetIDSelector.clone(quality=cms.string('LOOSE')),
src=cms.InputTag("slimmedJets"))
process.goodFatJets = cms.EDFilter(
"JetIDSelectionFunctorFilter",
filterParams=jetIDSelector.clone(quality=cms.string('LOOSE')),
src=cms.InputTag(jetsAK8))
################## Ntuplizer ###################
process.ntuplizer = cms.EDAnalyzer(
"Ntuplizer",
runOnMC=cms.bool(runOnMC),
doPruning=cms.bool(doAK8prunedReclustering),
doTausBoosted=cms.bool(doSemileptonicTausBoosted),
vertices=cms.InputTag("offlineSlimmedPrimaryVertices"),
muons=cms.InputTag("slimmedMuons"),
HCalVeto=cms.double(999.0),
RelIso=cms.double(0.2),
cutsToIgnore=cms.vstring('Chi2', 'D0', 'ED0', 'SD0',
'NHits', 'NValMuHits',
'ECalVeto', 'HCalVeto'),
RecalcFromBeamSpot=cms.bool(False),
beamLineSrc=cms.InputTag("offlineBeamSpot")),
# loose electrons
electronIdLoose=cms.PSet(version=cms.string('FIRSTDATA'),
D0=cms.double(999.0),
ED0=cms.double(999.0),
SD0=cms.double(999.0),
RelIso=cms.double(0.2),
cutsToIgnore=cms.vstring('D0', 'ED0', 'SD0')),
# loose jets
jetIdLoose=jetIDSelector.clone(),
pfjetIdLoose=pfJetIDSelector.clone(),
# kinematic cuts
minJets=cms.int32(1),
muPlusJets=cms.bool(True),
ePlusJets=cms.bool(False),
muPtMin=cms.double(20.0),
muEtaMax=cms.double(2.1),
eleEtMin=cms.double(20.0),
eleEtaMax=cms.double(2.4),
muPtMinLoose=cms.double(10.0),
muEtaMaxLoose=cms.double(2.5),
eleEtMinLoose=cms.double(15.0),
eleEtaMaxLoose=cms.double(2.5),
jetPtMin=cms.double(30.0),
jetEtaMax=cms.double(2.4),
jecLevelsAK4chs = [
'JEC/MCRUN2_74_V9::All_L1FastJet_AK4PFchs.txt',
'JEC/MCRUN2_74_V9::All_L2Relative_AK4PFchs.txt',
'JEC/MCRUN2_74_V9::All_L3Absolute_AK4PFchs.txt'
]
if corrMETonTheFly:
jecLevelsAK4 = [
'JEC/MCRUN2_74_V9::All_L1FastJet_AK4PF.txt',
'JEC/MCRUN2_74_V9::All_L2Relative_AK4PF.txt',
'JEC/MCRUN2_74_V9::All_L3Absolute_AK4PF.txt'
]
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector
process.goodSlimmedJets = cms.EDFilter("JetIDSelectionFunctorFilter",
filterParams=jetIDSelector.clone(quality = cms.string('LOOSE')),
src=cms.InputTag("slimmedJets")
)
process.goodFatJets = cms.EDFilter("JetIDSelectionFunctorFilter",
filterParams=jetIDSelector.clone(quality = cms.string('LOOSE')),
src=cms.InputTag(jetsAK8)
)
################## Ntuplizer ###################
process.ntuplizer = cms.EDAnalyzer("Ntuplizer",
runOnMC = cms.bool(runOnMC),
doPruning = cms.bool(doAK8prunedReclustering),
doTausBoosted = cms.bool(doSemileptonicTausBoosted),
vertices = cms.InputTag("offlineSlimmedPrimaryVertices"),
muons = cms.InputTag("slimmedMuons"),
"PoolSource",
fileNames=cms.untracked.vstring(
"dcap:///pnfs/cms/WAX/11/store/user/rappocc/QCD_Pt100to250-madgraph/shyft_35x_v1/7d6b33b2c5d85d3512a9f559ea910755/multijets_1_1.root"
),
)
## Maximal Number of Events
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(100))
from Analysis.JetAnalysis.multijetStudies_cfi import jetStudies as inputJetStudies
from Analysis.JetAnalysis.multijetStudies_cfi import pfJetStudies as inputPFJetStudies
from PhysicsTools.SelectorUtils.pfJetIDSelector_cfi import pfJetIDSelector as inputPFJetIDSelector
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector as inputJetIDSelector
process.TFileService = cms.Service("TFileService", fileName=cms.string("multijetStudies.root"))
process.multijetAna = cms.EDAnalyzer(
"EDMultijetAnalysis",
pfJetIDSelector=inputPFJetIDSelector.clone(),
jetIDSelector=inputJetIDSelector.clone(),
jetStudies=inputJetStudies.clone(ptMin=cms.double(50.0)),
pfJetStudies=inputPFJetStudies.clone(ptMin=cms.double(50.0)),
plotParameters=cms.PSet(doTracks=cms.bool(False), useMC=cms.bool(False)),
)
process.p = cms.Path(process.multijetAna)
process.MessageLogger.cerr.FwkReport.reportEvery = 10000
def selectObjects(process, cms):
print '=' * 60
print "Setting up selected objects"
print '=' * 60
###############################
#### Selections Setup #########
###############################
### AK5 Jets
#calo jets
process.selectedPatJets.cut = cms.string("pt > 20 & abs(eta) < 2.5")
#PF jets
process.selectedPatJetsPFlow.cut = cms.string("pt > 20 & abs(eta) < 2.5")
process.patJetsPFlow.addTagInfos = True
process.patJetsPFlow.tagInfoSources = cms.VInputTag(
cms.InputTag("secondaryVertexTagInfosAODPFlow"))
process.patJetsPFlow.userData.userFunctions = cms.vstring(
"? hasTagInfo('secondaryVertex') && tagInfoSecondaryVertex('secondaryVertex').nVertices() > 0 ? "
"tagInfoSecondaryVertex('secondaryVertex').secondaryVertex(0).p4().mass() : 0"
)
process.patJetsPFlow.userData.userFunctionLabels = cms.vstring(
'secvtxMass')
# electrons
process.selectedPatElectrons.cut = cms.string('pt > 10.0 & abs(eta) < 2.5')
process.patElectrons.embedTrack = cms.bool(True)
process.selectedPatElectronsPFlow.cut = cms.string(
'pt > 10.0 & abs(eta) < 2.5')
process.patElectronsPFlow.embedTrack = cms.bool(True)
process.selectedPatElectronsLoosePFlow.cut = cms.string(
'pt > 10.0 & abs(eta) < 2.5')
process.patElectronsLoosePFlow.embedTrack = cms.bool(True)
# muons
process.selectedPatMuons.cut = cms.string('pt > 10.0 & abs(eta) < 2.5')
process.patMuons.embedTrack = cms.bool(True)
process.selectedPatMuonsPFlow.cut = cms.string(
"pt > 10.0 & abs(eta) < 2.5")
process.patMuonsPFlow.embedTrack = cms.bool(True)
process.selectedPatMuonsLoosePFlow.cut = cms.string(
"pt > 10.0 & abs(eta) < 2.5")
process.patMuonsLoosePFlow.embedTrack = cms.bool(True)
# taus
process.selectedPatTausPFlow.cut = cms.string("pt > 10.0 & abs(eta) < 3")
process.selectedPatTaus.cut = cms.string("pt > 10.0 & abs(eta) < 3")
process.patTausPFlow.isoDeposits = cms.PSet()
process.patTaus.isoDeposits = cms.PSet()
# photons
process.patPhotonsPFlow.isoDeposits = cms.PSet()
process.patPhotons.isoDeposits = cms.PSet()
# Apply jet ID to all of the jets upstream. We aren't going to screw around
# with this, most likely. So, we don't really to waste time with it
# at the analysis level.
#default is version=PURE09, quality=loose ID
#https://twiki.cern.ch/twiki/bin/view/CMS/JetID#Documentation
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector
process.goodPatJets = cms.EDFilter("JetIDSelectionFunctorFilter",
filterParams=jetIDSelector.clone(),
src=cms.InputTag("selectedPatJets"))
from PhysicsTools.SelectorUtils.pfJetIDSelector_cfi import pfJetIDSelector
process.goodPatJetsPFlow = cms.EDFilter(
"PFJetIDSelectionFunctorFilter",
filterParams=pfJetIDSelector.clone(),
src=cms.InputTag("selectedPatJetsPFlow"))
process.selectedPatJetsForMETtype1p2Corr.src = cms.InputTag(
'selectedPatJetsPFlow')
process.selectedPatJetsForMETtype2Corr.src = cms.InputTag(
'selectedPatJetsPFlow')
import FWCore.ParameterSet.Config as cms
from PhysicsTools.SelectorUtils.jetIDSelector_cfi import jetIDSelector
from PhysicsTools.SelectorUtils.pfJetIDSelector_cfi import pfJetIDSelector
from Analysis.BoostedTopAnalysis.BoostedTopWTagParams_cfi import boostedTopWTagParams as inputBoostedTopWTagParams
MistagMakerParams = cms.PSet(
dijetSelectorParams=cms.PSet(
jetSrc=cms.InputTag("selectedPatJetsCA8PrunedCalo"),
pfJetSrc=cms.InputTag("goodPatJetsCA8PrunedPF"),
metSrc=cms.InputTag("patMETs"),
pfMetSrc=cms.InputTag("patMETsPF"),
useCalo=cms.bool(False),
minNJets=cms.uint32(2),
ptMin=cms.double(200),
etaMax=cms.double(2.4),
htMin=cms.double(0.0),
ht2Min=cms.double(0.0),
),
jetIDParams=cms.PSet(jetIDSelector.clone()),
pfJetIDParams=cms.PSet(pfJetIDSelector.clone()),
boostedTopWTagParams=cms.PSet(inputBoostedTopWTagParams),
wMassMin=cms.double(60),
wMassMax=cms.double(100),
jetPtMin=cms.double(200),
)