process.PATCMGSequence += process.selectedPatJetsCA8CHSprunedForBoostedTaus
process.PATCMGSequence += process.selectedPatJetsCA8CHSwithQJetsForBoostedTaus
patEventContentCMG+=['drop patJets_selectedPatJetsCA8CHSpruned_*_*']
patEventContentCMG+=['drop patJets_selectedPatJetsCA8CHSwithQjets_*_*']
######ADD PU JET ID
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
#process.puJetIdAK7CHS = puJetId.clone(
# jets ='selectedPatJetsAK7CHSwithQjets',
# jec = 'AK7chs'
# )
#process.PATCMGSequence += process.puJetIdAK7CHS
#patEventContentCMG+=['keep *_puJetIdAK7CHS_*_*']
process.puJetIdCA8CHS = puJetId.clone(
jets ='selectedPatJetsCA8CHSwithQjets',
jec = 'AK7chs'
)
process.puJetIdCA8CHSwithQJetsForBoostedTaus = puJetId.clone(
jets ='selectedPatJetsCA8CHSwithQJetsForBoostedTaus',
jec = 'AK7chs'
)
from CMGTools.External.pujetidproducer_cfi import stdalgos_4x, stdalgos_5x, stdalgos, cutbased, chsalgos_4x, chsalgos_5x, chsalgos
process.puJetMvaAK5CHS= puJetMva.clone(
jetids = cms.InputTag("puJetIdCHS"),
#### jets ='selectedPatJetsCHS',
jets ='patJetsWithVarCHS',
algos = chsalgos
)
cutoff=cms.double(100.0),
jetRad=cms.double(0.7),
jetAlgo=cms.string("AK"),
preclustering=cms.int32(30),
)
if not runQJets:
process.selectedPatJetsAK7CHSwithQjets.cutoff = 100000.0
process.PATCMGSequence += process.selectedPatJetsAK7CHSwithQjets
patEventContentCMG += ['drop patJets_selectedPatJetsAK7CHSwithNsub_*_*']
process.PATCMGSequence += process.selectedPatJetsCA8CHSwithQjets
patEventContentCMG += ['drop patJets_selectedPatJetsCA8CHSwithNsub_*_*']
######ADD PU JET ID
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
process.puJetIdAK7CHS = puJetId.clone(jets='selectedPatJetsAK7CHSwithQjets',
jec='AK7chs')
process.PATCMGSequence += process.puJetIdAK7CHS
patEventContentCMG += ['keep *_puJetIdAK7CHS_*_*']
process.puJetIdCA8CHS = puJetId.clone(jets='selectedPatJetsCA8CHSwithQjets',
jec='AK7chs')
from CMGTools.External.pujetidproducer_cfi import stdalgos_4x, stdalgos_5x, stdalgos, cutbased, chsalgos_4x, chsalgos_5x, chsalgos
process.puJetMvaAK5CHS = puJetMva.clone(
jetids=cms.InputTag("puJetIdCHS"),
#### jets ='selectedPatJetsCHS',
jets='patJetsWithVarCHS',
algos=chsalgos)
process.puJetMvaAK7CHS = puJetMva.clone(jetids=cms.InputTag("puJetIdAK7CHS"),
jets='selectedPatJetsAK7CHSwithQjets',
process.load("ExoDiBosonResonances.PATtupleProduction.PAT_boostedtaus_cff")
process.PATCMGSequence += process.selectedPatJetsCA8CHSprunedForBoostedTaus
process.PATCMGSequence += process.selectedPatJetsCA8CHSwithQJetsForBoostedTaus
patEventContentCMG += ['drop patJets_selectedPatJetsCA8CHSpruned_*_*']
patEventContentCMG += ['drop patJets_selectedPatJetsCA8CHSwithQjets_*_*']
######ADD PU JET ID
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
#process.puJetIdAK7CHS = puJetId.clone(
# jets ='selectedPatJetsAK7CHSwithQjets',
# jec = 'AK7chs'
# )
#process.PATCMGSequence += process.puJetIdAK7CHS
#patEventContentCMG+=['keep *_puJetIdAK7CHS_*_*']
process.puJetIdCA8CHS = puJetId.clone(jets='selectedPatJetsCA8CHSwithQjets',
jec='AK7chs')
process.puJetIdCA8CHSwithQJetsForBoostedTaus = puJetId.clone(
jets='selectedPatJetsCA8CHSwithQJetsForBoostedTaus', jec='AK7chs')
from CMGTools.External.pujetidproducer_cfi import stdalgos_4x, stdalgos_5x, stdalgos, cutbased, chsalgos_4x, chsalgos_5x, chsalgos
process.puJetMvaAK5CHS = puJetMva.clone(
jetids=cms.InputTag("puJetIdCHS"),
#### jets ='selectedPatJetsCHS',
jets='patJetsWithVarCHS',
algos=chsalgos)
#process.puJetMvaAK7CHS= puJetMva.clone(
# jetids = cms.InputTag("puJetIdAK7CHS"),
# jets ='selectedPatJetsAK7CHSwithQjets',
#process.PATCMGSequence += process.selectedPatJetsAK7CHSwithQjets
#patEventContentCMG+=['drop patJets_selectedPatJetsAK7CHSwithNsub_*_*']
process.PATCMGSequence += process.selectedPatJetsCA8CHSwithQjets
patEventContentCMG+=['drop patJets_selectedPatJetsCA8CHSwithNsub_*_*']
######ADD PU JET ID
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
#process.puJetIdAK7CHS = puJetId.clone(
# jets ='selectedPatJetsAK7CHSwithQjets',
# jec = 'AK7chs'
# )
#process.PATCMGSequence += process.puJetIdAK7CHS
#patEventContentCMG+=['keep *_puJetIdAK7CHS_*_*']
process.puJetIdCA8CHS = puJetId.clone(
jets ='selectedPatJetsCA8CHSwithQjets',
jec = 'AK7chs'
)
from CMGTools.External.pujetidproducer_cfi import stdalgos_4x, stdalgos_5x, stdalgos, cutbased, chsalgos_4x, chsalgos_5x, chsalgos
process.puJetMvaAK5CHS= puJetMva.clone(
jetids = cms.InputTag("puJetIdCHS"),
#### jets ='selectedPatJetsCHS',
jets ='patJetsWithVarCHS',
algos = chsalgos
)
#process.puJetMvaAK7CHS= puJetMva.clone(
# jetids = cms.InputTag("puJetIdAK7CHS"),
# jets ='selectedPatJetsAK7CHSwithQjets',
# algos = chsalgos
jets = cms.InputTag("selectedPatJetsAK5"),
rho = cms.InputTag("kt6PFJetsForIso:rho"),
jec = cms.string('ak5PFL1FastL2L3'),
isPatJet = cms.bool(True),
)
process.customPFJets = cms.EDProducer(
'PFJetUserData',
JetInputCollection=cms.untracked.InputTag("selectedPatJetsAK5"),
is2012Data=cms.untracked.bool(True),
qgMap=cms.untracked.InputTag("qglAK5PFCHS"),
Verbosity=cms.untracked.bool(False)
)
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
process.puJetIdAK5 = puJetId.clone( jets = 'customPFJets')
process.puJetMvaAK5= puJetMva.clone(
jetids = cms.InputTag("puJetIdAK5"),
jets = 'customPFJets',
)
process.puJetIdSequenceAK5 = cms.Sequence(process.puJetIdAK5*process.puJetMvaAK5)
# central jets for filtering and Z->jj candidates
process.customPFJetsCentral = cms.EDFilter(
"PATJetSelector",
src = cms.InputTag("customPFJets"),
cut = cms.string("abs(eta) < 2.4")
)
############## "Classic" PAT Muons and Electrons ########################
process.ak5PFJetsL1FastL2L3 = cms.EDProducer('PFJetCorrectionProducer',
src = cms.InputTag('ak5PFJets'),
correctors = cms.vstring(jetcorrection)
)
process.ak5PFJetsL1FastL2L3NoV = cms.EDProducer('PFJetCorrectionProducer',
src = cms.InputTag('ak5PFJetsNoV'),
correctors = cms.vstring(jetcorrection)
)
#PU jet ID
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
process.jetPUId = puJetId.clone(
jets = cms.InputTag("ak5PFJetsL1FastL2L3NoV"),
applyJec = cms.bool(False),
inputIsCorrected = cms.bool(True),
)
process.jetPUMVA = puJetMva.clone(
jets = cms.InputTag("ak5PFJetsL1FastL2L3NoV"),
jetids = cms.InputTag("jetPUId"),
applyJec = cms.bool(False),
inputIsCorrected = cms.bool(True),
)
process.jetPUIdSequence = cms.Sequence(process.jetPUId * process.jetPUMVA)
# MET corrections (type I)
process.load("EWKV.Analyzer.pfMETCorrections_cff")
# jetCorrections=('AK5PFchs', ['L1FastJet', 'L2Relative', 'L3Absolute']),
# pvCollection=cms.InputTag('offlinePrimaryVertices')
# )
#process.pfPileUpPFlow.checkClosestZVertex = False
#getattr(process, "patPF2PATSequence"+postfix).remove(getattr(process, "cleanPatCandidates"+postfix))
#getattr(process, "patPF2PATSequence"+postfix).remove(getattr(process, "countPatCandidates"+postfix))
#process.patseq = cms.Sequence(getattr(process,"patPF2PATSequence"+postfix))
#process.pat_step = cms.Path(process.patseq)
######Pileup Jet ID (Hendrik)
#https://twiki.cern.ch/twiki/bin/view/CMS/PileupJetID#Running_the_algorithm_on_reco_je
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
process.recoPuJetId = puJetId.clone(
jets = cms.InputTag("ak5PFJets"),
applyJec = cms.bool(True),
inputIsCorrected = cms.bool(False),
)
process.recoPuJetMva = puJetMva.clone(
jets = cms.InputTag("ak5PFJets"),
jetids = cms.InputTag("recoPuJetId"),
applyJec = cms.bool(True),
inputIsCorrected = cms.bool(False),
)
process.recoPuJetIdSqeuence = cms.Sequence(process.recoPuJetId * process.recoPuJetMva)
process.jetpuid_step = cms.Path(process.recoPuJetIdSqeuence)
######ROOTMAKER
process.makeroottree = cms.EDAnalyzer("RootMaker",
process.kt6PFJetsForIsolation.Rho_EtaMax = cms.double(2.5)
#
# particle flow isolation
#
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso, setupPFMuonIso
process.eleIsoSequence = setupPFElectronIso(process, 'gsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence + process.eleIsoSequence)
###########################################################################################
# MVA PU jets ID
process.load("CMGTools.External.pujetidsequence_cff")
from CMGTools.External.pujetidsequence_cff import puJetId
process.recoPuJetId = puJetId.clone(
jets = cms.InputTag("ak5PFJetsL1FastL2L3Residual"),
applyJec = cms.bool(False),
inputIsCorrected = cms.bool(True)
)
from CMGTools.External.pujetidsequence_cff import puJetMva
process.recoPuJetMva = puJetMva.clone(
jets = cms.InputTag("ak5PFJetsL1FastL2L3Residual"),
jetids = cms.InputTag("recoPuJetId"),
applyJec = cms.bool(False),
inputIsCorrected = cms.bool(True)
)
################## filters ##########################################################
## The good primary vertex filter
process.primaryVertexFilter = cms.EDFilter(
from JetMETCorrections.Configuration.JetCorrectionServices_cff import *
process.load(
'JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff')
from JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff import *
process.ak5PFJetsCorrection = ak5PFJetsL1FastL2L3.clone()
process.ak5PFJetsCorrection.src = cms.InputTag('hTozzTo4leptonsPFJetSelector')
process.ak5PFJetsCorrectionData = ak5PFJetsL1FastL2L3Residual.clone()
process.ak5PFJetsCorrectionData.src = cms.InputTag(
'hTozzTo4leptonsPFJetSelector')
process.load('CMGTools.External.pujetidsequence_cff')
from CMGTools.External.pujetidsequence_cff import puJetId, puJetMva
process.recoPuJetIdMC = puJetId.clone(
jets=cms.InputTag("ak5PFJetsCorrection"), )
process.recoPuJetMvaMC = puJetMva.clone(
jets=cms.InputTag("ak5PFJetsCorrection"),
jetids=cms.InputTag("recoPuJetIdMC"),
)
process.recoPuJetIdMCsequence = cms.Sequence(process.recoPuJetIdMC *
process.recoPuJetMvaMC)
process.recoPuJetIdData = puJetId.clone(
jets=cms.InputTag("ak5PFJetsCorrectionData"), )
process.recoPuJetMvaData = puJetMva.clone(
jets=cms.InputTag("ak5PFJetsCorrectionData"),
jetids=cms.InputTag("recoPuJetIdData"),
)
process.recoPuJetIdDatasequence = cms.Sequence(process.recoPuJetIdData *
process.recoPuJetMvaData)