process.logErrorTooManySeeds
)
######################
### NoPU & MVA MET ###
######################
####### SET THIS TO TRUE TO RUN THE SEQUENCE (CHECK THE README FILE FOR INSTRUCTIONS) ########
runNoPUMVAMetSequence = False
##############################################################################################
if runNoPUMVAMetSequence:
from RecoJets.JetProducers.PileupJetIDParams_cfi import JetIdParams
from JetMETCorrections.Configuration.JetCorrectionProducers_cff import ak5PFJetsL1
if isRealData:
process.ak5PFJetsL123Corrected = ak5PFJetsL1.clone(correctors = ['ak5PFL1FastL2L3Residual'])
else:
process.ak5PFJetsL123Corrected = ak5PFJetsL1.clone(correctors = ['ak5PFL1FastL2L3'])
process.recoPuJetIdCorrected = pileupJetIdProducer.clone(
jets = cms.InputTag("ak5PFJetsL123Corrected"),
algos = cms.VPSet(full_5x, cutbased),
residualsTxt = cms.FileInPath("RecoJets/JetProducers/data/mva_JetID_v1.weights.xml") # not used, but has to point to an existing file
)
# NoPU MET
from JetMETCorrections.METPUSubtraction.noPileUpPFMET_cff import noPileUpPFMEtData, noPileUpPFMEt
process.pfNoPileUpMetData = noPileUpPFMEtData.clone(
srcJets = cms.InputTag('ak5PFJetsL123Corrected'),
srcJetIds = cms.InputTag('recoPuJetIdCorrected', 'fullId'),
srcPFCandToVertexAssociations = cms.InputTag('pfCandidateToVertexAssoc'),
def configure( dataset, sourceNames=[], hltPaths=[], maxEvents = -1, outputName = '', runNoPUMVAMetSequence=False ):
""" This is the main routine used to configure the nTuplizer. Please check
the correctness of all relevant default parameters.
dataset: Please choose one of the entries in the lists.
sourceNames: list of AOD filenames
hltPaths: list of trigger hlt paths (will be ignored for MC)
maxEvents: maximal number of events
runNoPUMVAMetSequence: please read README for more details
"""
collisionDatasets = [
'52xPrompt', # Run2012[AB]{C}-PromptReco{_v1}
'52x23May2012', # Run2012A-23May2012
'53xPromptC', # Run2012C-PromptReco-v2
'53xPromptD', # Run2012D-PromptReco
'53x13July2012', # Run2012[AB]-13Jul2012
'53x06Aug2012', # Run2012A-06Aug2012
'53x24Aug2012', # Run2012C-24Aug2012
'53x11Dec2012', # Run2012C-EcalRecover_11Dec2012
'53x16Jan2013', # Run2012D-16Jan2013
'53x22Jan2013' # Run2012[ABCD]-22Jan2013
]
mcDatasets = [
'52xFullSim',
'52xFastSim',
'53xFullSim',
'53xFastSim'
]
isRealData = dataset in collisionDatasets
isMC = dataset in mcDatasets
isFastSim = isMC and 'FastSim' in dataset
is53x = '53x' in dataset
is52x = '52x' in dataset
if not isRealData and not isMC:
raise RuntimeError("Dataset " + dataset + " not defined")
import FWCore.ParameterSet.Config as cms
##########################
### Initialize process ###
##########################
process = cms.Process("RA3")
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(maxEvents))
process.options = cms.untracked.PSet(
wantSummary = cms.untracked.bool(False),
FailPath = cms.untracked.vstring(
'FatalRootError'
)
)
process.source = cms.Source("PoolSource",
noEventSort = cms.untracked.bool(True),
duplicateCheckMode = cms.untracked.string('noDuplicateCheck'),
fileNames = cms.untracked.vstring(sourceNames)
)
#####################
### MessageLogger ###
#####################
process.load('FWCore.MessageService.MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
process.MessageLogger.suppressWarning = cms.untracked.vstring(
'newSecondaryVertexTagInfos',
'chsSecondaryVertexTagInfos',
'pfCandidateToVertexAssoc',
'manystripclus53X',
'toomanystripclus53X'
)
process.MessageLogger.suppressError = cms.untracked.vstring('ecalLaserCorrFilter')
process.MessageLogger.categories.append('SusyNtuplizer')
process.MessageLogger.cerr.SusyNtuplizer = cms.untracked.PSet( limit = cms.untracked.int32(100) )
#############################
### Conditions & Services ###
#############################
process.load('Configuration.Geometry.GeometryIdeal_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
# The following global tags are the latest as of April 8th 2013. The user is strongly recommended to check the
# SWGuideFrontierConditions twiki before running the configuration.
if dataset == '52xPrompt':
process.GlobalTag.globaltag = 'GR_P_V39_AN3::All'
elif dataset == '52x23May2012':
process.GlobalTag.globaltag = 'FT_P_V32B_AN4::All'
elif dataset == '53xPromptC' or dataset == '53xPromptD':
process.GlobalTag.globaltag = 'GR_P_V42_AN4::All'
elif dataset == '53x13July2012':
process.GlobalTag.globaltag = 'FT_53_V6C_AN4::All'
elif dataset == '53x06Aug2012':
process.GlobalTag.globaltag = 'FT_53_V6C_AN4::All'
elif dataset == '53x24Aug2012':
process.GlobalTag.globaltag = 'FT53_V10A_AN4::All'
elif dataset == '53x11Dec2012':
process.GlobalTag.globaltag = 'FT_P_V42C_AN4::All'
elif dataset == '53x16Jan2013':
process.GlobalTag.globaltag = 'FT_P_V43E_AN4::All'
elif dataset == '53x22Jan2013':
process.GlobalTag.globaltag = 'FT_53_V21_AN4::All'
elif isMC and is52x:
process.GlobalTag.globaltag = 'START52_V16::All'
elif isMC and is53x:
process.GlobalTag.globaltag = 'START53_V25::All'
#####################
### SusyNtuplizer ###
#####################
process.load("SUSYPhotonAnalysis.SusyNtuplizer.susyNtuplizer_cfi")
process.susyNtuplizer.debugLevel = 0
process.susyNtuplizer.isFastSim = isFastSim
process.susyNtuplizer.caloJetCollectionTags = []
if outputName:
process.susyNtuplizer.outputFileName = outputName
#########################
### HLT result filter ###
#########################
if len(hltPaths) != 0:
process.load("HLTrigger.HLTfilters.hltHighLevel_cfi")
process.hltHighLevel.HLTPaths = hltPaths
else:
process.hltHighLevel = cms.Sequence() # placeholder
##########################################
### Good vertex collection (transient) ###
##########################################
process.goodVertices = cms.EDFilter("VertexSelector",
filter = cms.bool(False),
src = cms.InputTag("offlinePrimaryVertices"),
cut = cms.string("!isFake && ndof >= 4 && abs(z) <= 24 && position.rho < 2")
)
process.primaryVertex = cms.EDFilter("PATSingleVertexSelector",
mode = cms.string('firstVertex'),
vertices = cms.InputTag('goodVertices'),
filter = cms.bool(False)
)
process.vertexSelectionSequence = cms.Sequence(
process.goodVertices +
process.primaryVertex
)
###########################
### PU Rho calculations ###
###########################
from RecoJets.JetProducers.kt4PFJets_cfi import kt4PFJets
#Calculate rho restricted to barrel for photon pileup subtraction
process.kt6PFJetsRhoBarrelOnly = kt4PFJets.clone(
src = cms.InputTag('particleFlow'),
rParam = cms.double(0.6),
#Eta range of jets to be considered for Rho calculation
#Should be at most (jet acceptance - jet radius)
doRhoFastjet = cms.bool(True),
Rho_EtaMax=cms.double(1.4442),
#Eta range of ghost jets to be considered for Rho calculation - must be greater than Rho_EtaMax
Ghost_EtaMax=cms.double(2.5)
)
#Calculate rho restricted to barrel for photon pileup subtraction
process.kt6PFJetsRho25 = kt4PFJets.clone(
src = cms.InputTag('particleFlow'),
rParam = cms.double(0.6),
#Eta range of jets to be considered for Rho calculation
#Should be at most (jet acceptance - jet radius)
doRhoFastjet = cms.bool(True),
Rho_EtaMax=cms.double(2.5)
)
process.puRhoSequence = cms.Sequence(
process.kt6PFJetsRhoBarrelOnly +
process.kt6PFJetsRho25
)
###############################
### PF-based reconstruction ###
###############################
# Run the full event reconstruction starting from the full list of PFCandidates (particleFlow).
# Serves as inputs to isoDeposit, CHS jets, and NoPUMET.
process.load("CommonTools.ParticleFlow.PFBRECO_cff")
process.pfPileUp.Vertices = cms.InputTag("goodVertices")
process.pfPileUp.checkClosestZVertex = False
process.pfBasedRecoSequence = cms.Sequence(
process.pfNoPileUpSequence +
process.pfParticleSelectionSequence +
process.pfPhotonSequence +
process.pfMuonSequence +
process.pfNoMuon +
process.pfElectronSequence +
process.pfNoElectron
)
####################################
### Photon & electron isoDeposit ###
####################################
# Runs PU identificaton over particleFlow, then calculates isodeposits around
# the given particles using the NoPU collection.
# The end product is equivalent to running setupPFIso in
# CommonTools.ParticleFlow.Tools.pfIsolation, just with fewer modules to run.
process.photonPFIsoDepositCharged = process.phPFIsoDepositCharged.clone(src = 'photons')
process.photonPFIsoDepositNeutral = process.phPFIsoDepositNeutral.clone(src = 'photons')
process.photonPFIsoDepositGamma = process.phPFIsoDepositGamma.clone(src = 'photons')
process.gsfElectronPFIsoDepositCharged = process.elPFIsoDepositCharged.clone(src = 'gsfElectrons')
process.gsfElectronPFIsoDepositNeutral = process.elPFIsoDepositNeutral.clone(src = 'gsfElectrons')
process.gsfElectronPFIsoDepositGamma = process.elPFIsoDepositGamma.clone(src = 'gsfElectrons')
process.photonPFIsoValueCharged03 = process.phPFIsoValueCharged03PFId.clone()
process.photonPFIsoValueCharged03.deposits[0].src = cms.InputTag('photonPFIsoDepositCharged')
process.photonPFIsoValueNeutral03 = process.phPFIsoValueNeutral03PFId.clone()
process.photonPFIsoValueNeutral03.deposits[0].src = cms.InputTag('photonPFIsoDepositNeutral')
process.photonPFIsoValueGamma03 = process.phPFIsoValueGamma03PFId.clone()
process.photonPFIsoValueGamma03.deposits[0].src = cms.InputTag('photonPFIsoDepositGamma')
process.gsfElectronPFIsoValueCharged03 = process.elPFIsoValueCharged03PFId.clone()
process.gsfElectronPFIsoValueCharged03.deposits[0].src = cms.InputTag('gsfElectronPFIsoDepositCharged')
process.gsfElectronPFIsoValueNeutral03 = process.elPFIsoValueNeutral03PFId.clone()
process.gsfElectronPFIsoValueNeutral03.deposits[0].src = cms.InputTag('gsfElectronPFIsoDepositNeutral')
process.gsfElectronPFIsoValueGamma03 = process.elPFIsoValueGamma03PFId.clone()
process.gsfElectronPFIsoValueGamma03.deposits[0].src = cms.InputTag('gsfElectronPFIsoDepositGamma')
process.photonIsoDepositSequence = cms.Sequence(
process.photonPFIsoDepositCharged +
process.photonPFIsoDepositNeutral +
process.photonPFIsoDepositGamma +
process.photonPFIsoValueCharged03 +
process.photonPFIsoValueNeutral03 +
process.photonPFIsoValueGamma03
)
process.gsfElectronIsoDepositSequence = cms.Sequence(
process.gsfElectronPFIsoDepositCharged +
process.gsfElectronPFIsoDepositNeutral +
process.gsfElectronPFIsoDepositGamma +
process.gsfElectronPFIsoValueCharged03 +
process.gsfElectronPFIsoValueNeutral03 +
process.gsfElectronPFIsoValueGamma03
)
process.pfIsolationSequence = cms.Sequence(
process.photonIsoDepositSequence +
process.gsfElectronIsoDepositSequence
)
###############################################
### PFchs (charged hadron subtraction) jets ###
###############################################
# pfJets comes from PFBRECO.
process.ak5PFchsJets = process.pfJets.clone()
process.ak5PFchsJets.doAreaFastjet = True
process.pfCHSJetSequence = cms.Sequence(
process.ak5PFchsJets
)
###########
### JEC ###
###########
process.load("JetMETCorrections.Configuration.DefaultJEC_cff")
process.load("JetMETCorrections.Configuration.JetCorrectionServices_cff")
# Setting up JEC ESProducers for ak5PFchs. This block will be included in the JetCorrectionServices_cff
# in a future tag by JetMET (April 8, 2013)
process.ak5PFchsL1Fastjet = process.ak5PFL1Fastjet.clone(algorithm = cms.string('AK5PFchs'))
process.ak5PFchsL2Relative = process.ak5PFL2Relative.clone(algorithm = cms.string('AK5PFchs'))
process.ak5PFchsL3Absolute = process.ak5PFL3Absolute.clone(algorithm = cms.string('AK5PFchs'))
process.ak5PFchsResidual = process.ak5PFResidual.clone(algorithm = cms.string('AK5PFchs'))
process.ak5PFchsL2L3 = cms.ESProducer('JetCorrectionESChain',
correctors = cms.vstring('ak5PFchsL2Relative', 'ak5PFchsL3Absolute')
)
process.ak5PFchsL2L3Residual = process.ak5PFchsL2L3.clone()
process.ak5PFchsL2L3Residual.correctors.append('ak5PFchsResidual')
process.ak5PFchsL1FastL2L3 = process.ak5PFchsL2L3.clone()
process.ak5PFchsL1FastL2L3.correctors.insert(0, 'ak5PFchsL1Fastjet')
process.ak5PFchsL1FastL2L3Residual = process.ak5PFchsL1FastL2L3.clone()
process.ak5PFchsL1FastL2L3Residual.correctors.append('ak5PFchsResidual')
###############################
### Corrected MET producers ###
###############################
process.load("JetMETCorrections.Type1MET.caloMETCorrections_cff")
process.load("JetMETCorrections.Type1MET.pfMETCorrections_cff")
if isRealData:
process.caloJetMETcorr.jetCorrLabel = cms.string("ak5CaloL2L3Residual")
process.pfJetMETcorr.jetCorrLabel = cms.string("ak5PFL1FastL2L3Residual")
else:
process.caloJetMETcorr.jetCorrLabel = cms.string("ak5CaloL2L3")
process.pfJetMETcorr.jetCorrLabel = cms.string("ak5PFL1FastL2L3")
from JetMETCorrections.Type1MET.pfMETsysShiftCorrections_cfi import pfMEtSysShiftCorr
process.pfMEtSysShiftCorr = pfMEtSysShiftCorr.clone(
srcVertices = cms.InputTag('goodVertices')
)
# http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/JetMETCorrections/Type1MET/python/pfMETsysShiftCorrections_cfi.py?revision=1.6&view=markup
# pfMEtSysShiftCorrParameters_2012runABCvsNvtx_data (_mc)
if isRealData:
process.pfMEtSysShiftCorr.parameter = cms.PSet(
px = cms.string("+0.2661 + 0.3217*Nvtx"),
py = cms.string("-0.2251 - 0.1747*Nvtx")
)
else:
process.pfMEtSysShiftCorr.parameter = cms.PSet(
px = cms.string("+0.1166 + 0.0200*Nvtx"),
py = cms.string("+0.2764 - 0.1280*Nvtx")
)
# Remove the track-based type 0 correction
process.producePFMETCorrections.remove(process.pfchsMETcorr)
# Use PF-based type 0 correction. Reference: AN-2012/333
# Note that type 0 and type 1 are not truly orthogonal; the implicit assumption is that the PU
# contributions are fully captured in the high-pt jet offset, low-pt jet, and non-clustering
# energy terms in the type 1 MET correction formula.
from CommonTools.RecoUtils.pfcand_assomap_cfi import PFCandAssoMap
from JetMETCorrections.Type1MET.pfMETCorrectionType0_cfi import pfMETcorrType0
# Full implementation requires the reconstruction of displaced vertex, which we are skipping here.
# (Displaced vertex reconstruction requires TrackExtra information and is impossible from AOD)
process.pfCandidateToVertexAssoc = PFCandAssoMap.clone()
# Attention: Type0PFMETcorrInputProducer does not produce sumEt of the correction.
process.pfType0MetCorrection = pfMETcorrType0.clone(
srcPFCandidateToVertexAssociations = cms.InputTag('pfCandidateToVertexAssoc'),
srcHardScatterVertex = cms.InputTag('primaryVertex')
)
process.pfType01CorrectedMet = process.pfType1CorrectedMet.clone(
srcType1Corrections = cms.VInputTag(
cms.InputTag('pfType0MetCorrection'),
cms.InputTag('pfJetMETcorr', 'type1')
)
)
process.pfType01p2CorrectedMet = process.pfType1p2CorrectedMet.clone(
srcType1Corrections = cms.VInputTag(
cms.InputTag('pfType0MetCorrection'),
cms.InputTag('pfJetMETcorr', 'type1')
)
)
process.pfSysShiftCorrectedMet = process.pfType1CorrectedMet.clone(
srcType1Corrections = cms.VInputTag(
cms.InputTag('pfMEtSysShiftCorr')
)
)
process.pfType01SysShiftCorrectedMet = process.pfType1CorrectedMet.clone(
srcType1Corrections = cms.VInputTag(
cms.InputTag('pfType0MetCorrection'),
cms.InputTag('pfJetMETcorr', 'type1'),
cms.InputTag('pfMEtSysShiftCorr')
)
)
process.correctedMetSequence = cms.Sequence(
process.produceCaloMETCorrections +
process.producePFMETCorrections +
process.pfCandidateToVertexAssoc +
process.pfType0MetCorrection +
process.pfType01CorrectedMet +
process.pfType01p2CorrectedMet +
process.pfMEtSysShiftCorr +
process.pfSysShiftCorrectedMet +
process.pfType01SysShiftCorrectedMet
)
#############################
### MVA-based electron ID ###
#############################
process.load('EGamma.EGammaAnalysisTools.electronIdMVAProducer_cfi')
process.eidMVASequence = cms.Sequence(
process.mvaTrigV0 *
process.mvaNonTrigV0
)
#################
### PU jet ID ###
#################
from RecoJets.JetProducers.PileupJetIDParams_cfi import full_5x, full_5x_chs, cutbased
from RecoJets.JetProducers.PileupJetID_cfi import pileupJetIdProducer, pileupJetIdProducerChs
process.recoPuJetId = pileupJetIdProducer.clone(
jets = cms.InputTag("ak5PFJets"),
algos = cms.VPSet(full_5x, cutbased),
applyJec = cms.bool(True),
inputIsCorrected = cms.bool(False),
residualsTxt = cms.FileInPath("RecoJets/JetProducers/data/mva_JetID_v1.weights.xml") # not used, but has to point to an existing file
)
process.recoPuJetIdChs = pileupJetIdProducerChs.clone(
jets = cms.InputTag("ak5PFchsJets"),
algos = cms.VPSet(full_5x_chs, cutbased),
applyJec = cms.bool(True),
inputIsCorrected = cms.bool(False),
residualsTxt = cms.FileInPath("RecoJets/JetProducers/data/mva_JetID_v1.weights.xml") # not used, but has to point to an existing file
)
process.recoPuJetIdSequence = cms.Sequence(
process.recoPuJetId +
process.recoPuJetIdChs
)
##################################
### Jet-parton matching for MC ###
##################################
if isMC:
process.myPartons = cms.EDProducer("PartonSelector",
withLeptons = cms.bool(False),
src = cms.InputTag("genParticles")
)
process.flavourByRef = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("ak5PFJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("myPartons")
)
process.flavourAssociationAlg = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("flavourByRef"),
physicsDefinition = cms.bool(False)
)
process.flavourAssociationPhy = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("flavourByRef"),
physicsDefinition = cms.bool(True)
)
process.flavourByRefCHS = cms.EDProducer("JetPartonMatcher",
jets = cms.InputTag("ak5PFchsJets"),
coneSizeToAssociate = cms.double(0.3),
partons = cms.InputTag("myPartons")
)
process.flavourAssociationCHSAlg = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("flavourByRefCHS"),
physicsDefinition = cms.bool(False)
)
process.flavourAssociationCHSPhy = cms.EDProducer("JetFlavourIdentifier",
srcByReference = cms.InputTag("flavourByRefCHS"),
physicsDefinition = cms.bool(True)
)
process.JetFlavourMatchingSequence = cms.Sequence(
process.myPartons *
process.flavourByRef *
process.flavourByRefCHS *
process.flavourAssociationAlg *
process.flavourAssociationPhy *
process.flavourAssociationCHSAlg *
process.flavourAssociationCHSPhy
)
else:
process.JetFlavourMatchingSequence = cms.Sequence() # placeholder
#################
### b-tagging ###
#################
# Re-run b-tagging with PFJets as input
# b-tagging general configuration
process.load("RecoJets.JetAssociationProducers.ic5PFJetTracksAssociatorAtVertex_cfi")
process.load("RecoBTag.Configuration.RecoBTag_cff")
# create a new jets and tracks associaiton
process.newJetTracksAssociatorAtVertex = process.ic5PFJetTracksAssociatorAtVertex.clone(
jets = cms.InputTag("ak5PFJets"),
tracks = cms.InputTag("generalTracks")
)
process.chsJetTracksAssociatorAtVertex = process.ic5PFJetTracksAssociatorAtVertex.clone(
jets = cms.InputTag("ak5PFchsJets"),
tracks = cms.InputTag("generalTracks")
)
# impact parameter b-tag
process.newImpactParameterTagInfos = process.impactParameterTagInfos.clone(
jetTracks = cms.InputTag("newJetTracksAssociatorAtVertex")
)
process.chsImpactParameterTagInfos = process.impactParameterTagInfos.clone(
jetTracks = cms.InputTag("chsJetTracksAssociatorAtVertex")
)
# TCHE
process.newTrackCountingHighEffBJetTags = process.trackCountingHighEffBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newImpactParameterTagInfos") )
)
process.chsTrackCountingHighEffBJetTags = process.trackCountingHighEffBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsImpactParameterTagInfos") )
)
# TCHP
process.newTrackCountingHighPurBJetTags = process.trackCountingHighPurBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newImpactParameterTagInfos") )
)
process.chsTrackCountingHighPurBJetTags = process.trackCountingHighPurBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsImpactParameterTagInfos") )
)
# JP
process.newJetProbabilityBJetTags = process.jetProbabilityBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newImpactParameterTagInfos") )
)
process.chsJetProbabilityBJetTags = process.jetProbabilityBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsImpactParameterTagInfos") )
)
# JBP
process.newJetBProbabilityBJetTags = process.jetBProbabilityBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newImpactParameterTagInfos") )
)
process.chsJetBProbabilityBJetTags = process.jetBProbabilityBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsImpactParameterTagInfos") )
)
# secondary vertex b-tag
process.newSecondaryVertexTagInfos = process.secondaryVertexTagInfos.clone(
trackIPTagInfos = cms.InputTag("newImpactParameterTagInfos")
)
process.chsSecondaryVertexTagInfos = process.secondaryVertexTagInfos.clone(
trackIPTagInfos = cms.InputTag("chsImpactParameterTagInfos")
)
# SSV
process.newSimpleSecondaryVertexBJetTags = process.simpleSecondaryVertexBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newSecondaryVertexTagInfos") )
)
process.chsSimpleSecondaryVertexBJetTags = process.simpleSecondaryVertexBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsSecondaryVertexTagInfos") )
)
# CSV
process.newCombinedSecondaryVertexBJetTags = process.combinedSecondaryVertexBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newImpactParameterTagInfos"),
cms.InputTag("newSecondaryVertexTagInfos")
)
)
process.chsCombinedSecondaryVertexBJetTags = process.combinedSecondaryVertexBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsImpactParameterTagInfos"),
cms.InputTag("chsSecondaryVertexTagInfos")
)
)
# CSVMVA
process.newCombinedSecondaryVertexMVABJetTags = process.combinedSecondaryVertexMVABJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newImpactParameterTagInfos"),
cms.InputTag("newSecondaryVertexTagInfos")
)
)
process.chsCombinedSecondaryVertexMVABJetTags = process.combinedSecondaryVertexMVABJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsImpactParameterTagInfos"),
cms.InputTag("chsSecondaryVertexTagInfos")
)
)
# soft electron b-tag
process.newSoftElectronTagInfos = process.softElectronTagInfos.clone(
jets = "ak5PFJets"
)
process.newSoftElectronBJetTags = process.softElectronBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newSoftElectronTagInfos") )
)
process.chsSoftElectronTagInfos = process.softElectronTagInfos.clone(
jets = "ak5PFchsJets"
)
process.chsSoftElectronBJetTags = process.softElectronBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsSoftElectronTagInfos") )
)
# soft muon b-tag
process.newSoftMuonTagInfos = process.softMuonTagInfos.clone(
jets = "ak5PFJets"
)
process.newSoftMuonBJetTags = process.softMuonBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newSoftMuonTagInfos") )
)
process.newSoftMuonByIP3dBJetTags = process.softMuonByIP3dBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newSoftMuonTagInfos") )
)
process.newSoftMuonByPtBJetTags = process.softMuonByPtBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("newSoftMuonTagInfos") )
)
process.chsSoftMuonTagInfos = process.softMuonTagInfos.clone(
jets = "ak5PFchsJets"
)
process.chsSoftMuonBJetTags = process.softMuonBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsSoftMuonTagInfos") )
)
process.chsSoftMuonByIP3dBJetTags = process.softMuonByIP3dBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsSoftMuonTagInfos") )
)
process.chsSoftMuonByPtBJetTags = process.softMuonByPtBJetTags.clone(
tagInfos = cms.VInputTag( cms.InputTag("chsSoftMuonTagInfos") )
)
process.newJetTracksAssociator = cms.Sequence(
process.newJetTracksAssociatorAtVertex
)
process.chsJetTracksAssociator = cms.Sequence(
process.chsJetTracksAssociatorAtVertex
)
process.newJetBtaggingIP = cms.Sequence(
process.newImpactParameterTagInfos * (
process.newTrackCountingHighEffBJetTags +
process.newTrackCountingHighPurBJetTags +
process.newJetProbabilityBJetTags +
process.newJetBProbabilityBJetTags
)
)
process.chsJetBtaggingIP = cms.Sequence(
process.chsImpactParameterTagInfos * (
process.chsTrackCountingHighEffBJetTags +
process.chsTrackCountingHighPurBJetTags +
process.chsJetProbabilityBJetTags +
process.chsJetBProbabilityBJetTags
)
)
process.newJetBtaggingSV = cms.Sequence(
process.newImpactParameterTagInfos *
process.newSecondaryVertexTagInfos * (
process.newSimpleSecondaryVertexBJetTags +
process.newCombinedSecondaryVertexBJetTags +
process.newCombinedSecondaryVertexMVABJetTags
)
)
process.chsJetBtaggingSV = cms.Sequence(
process.chsImpactParameterTagInfos *
process.chsSecondaryVertexTagInfos * (
process.chsSimpleSecondaryVertexBJetTags +
process.chsCombinedSecondaryVertexBJetTags +
process.chsCombinedSecondaryVertexMVABJetTags
)
)
process.newJetBtaggingEle = cms.Sequence(
process.softElectronCands *
process.newSoftElectronTagInfos *
process.newSoftElectronBJetTags
)
process.chsJetBtaggingEle = cms.Sequence(
process.softElectronCands *
process.chsSoftElectronTagInfos *
process.chsSoftElectronBJetTags
)
process.newJetBtaggingMu = cms.Sequence(
process.newSoftMuonTagInfos * (
process.newSoftMuonBJetTags +
process.newSoftMuonByIP3dBJetTags +
process.newSoftMuonByPtBJetTags
)
)
process.chsJetBtaggingMu = cms.Sequence(
process.chsSoftMuonTagInfos * (
process.chsSoftMuonBJetTags +
process.chsSoftMuonByIP3dBJetTags +
process.chsSoftMuonByPtBJetTags
)
)
process.newJetBtagging = cms.Sequence(
process.newJetBtaggingIP +
process.newJetBtaggingSV +
process.newJetBtaggingEle +
process.newJetBtaggingMu
)
process.chsJetBtagging = cms.Sequence(
process.chsJetBtaggingIP +
process.chsJetBtaggingSV +
process.chsJetBtaggingEle +
process.chsJetBtaggingMu
)
# These are fixes for the JetProbability b-tagger calibrations as recommended by BTV.
# See https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideBTagJetProbabilityCalibration#Calibration_in_52x_Data_and_MC
useSpecialBTagCalibration = True
if isRealData and is52x:
btagTrackProbability2DTag = "TrackProbabilityCalibration_2D_2012DataTOT_v1_offline"
btagTrackProbability3DTag = "TrackProbabilityCalibration_3D_2012DataTOT_v1_offline"
elif isRealData and is53x:
btagTrackProbability2DTag = "TrackProbabilityCalibration_2D_Data53X_v2"
btagTrackProbability3DTag = "TrackProbabilityCalibration_3D_Data53X_v2"
elif isMC and is53x:
btagTrackProbability2DTag = "TrackProbabilityCalibration_2D_MC53X_v2"
btagTrackProbability3DTag = "TrackProbabilityCalibration_3D_MC53X_v2"
else:
useSpecialBTagCalibration = False
if useSpecialBTagCalibration:
process.GlobalTag.toGet = cms.VPSet(
cms.PSet(record = cms.string("BTagTrackProbability2DRcd"),
tag = cms.string(btagTrackProbability2DTag),
connect = cms.untracked.string("frontier://FrontierPrep/CMS_COND_BTAU")
),
cms.PSet(record = cms.string("BTagTrackProbability3DRcd"),
tag = cms.string(btagTrackProbability3DTag),
connect = cms.untracked.string("frontier://FrontierPrep/CMS_COND_BTAU")
)
)
###########################
### Quark-gluon tagging ###
###########################
from QuarkGluonTagger.EightTeV.QGTagger_RecoJets_cff import goodOfflinePrimaryVerticesQG, QGTagger
process.goodOfflinePrimaryVerticesQG = goodOfflinePrimaryVerticesQG
process.QGTaggerAK5 = QGTagger.clone(
srcJets = cms.InputTag('ak5PFJets'),
srcRho = cms.InputTag('kt6PFJets', 'rho'),
srcRhoIso = cms.InputTag('kt6PFJetsRho25', 'rho')
)
process.QGTaggerAK5chs = QGTagger.clone(
srcJets = cms.InputTag('ak5PFchsJets'),
useCHS = cms.untracked.bool(True),
srcRho = cms.InputTag('kt6PFJets', 'rho'),
srcRhoIso = cms.InputTag('kt6PFJetsRho25', 'rho')
)
if isRealData:
process.QGTaggerAK5.jec = cms.untracked.string('ak5PFL1FastL2L3Residual')
process.QGTaggerAK5chs.jec = cms.untracked.string('ak5PFchsL1FastL2L3Residual')
else:
process.QGTaggerAK5.jec = cms.untracked.string('ak5PFL1FastL2L3')
process.QGTaggerAK5chs.jec = cms.untracked.string('ak5PFchsL1FastL2L3')
process.QGTaggingSequence = cms.Sequence(
process.goodOfflinePrimaryVerticesQG +
process.QGTaggerAK5 +
process.QGTaggerAK5chs
)
###################
### MET filters ###
###################
# HBHENoiseFilterResultProducer
process.load('CommonTools.RecoAlgos.HBHENoiseFilterResultProducer_cfi')
# HCAL laser events filter
process.load("RecoMET.METFilters.hcalLaserEventFilter_cfi")
process.hcalLaserEventFilter.taggingMode = cms.bool(True)
# EcalDeadCellTriggerPrimitiveFilter
process.load('RecoMET.METFilters.EcalDeadCellTriggerPrimitiveFilter_cfi')
process.EcalDeadCellTriggerPrimitiveFilter.taggingMode = cms.bool(True)
# EcalDeadCellBoundaryEnergyFilter
process.load('RecoMET.METFilters.EcalDeadCellBoundaryEnergyFilter_cfi')
process.EcalDeadCellBoundaryEnergyFilter.taggingMode = cms.bool(True)
# Tracking failure filter
process.load('RecoMET.METFilters.trackingFailureFilter_cfi')
process.trackingFailureFilter.taggingMode = cms.bool(True)
# EE Bad SC Filter
process.load('RecoMET.METFilters.eeBadScFilter_cfi')
process.eeBadScFilter.taggingMode = cms.bool(True)
# EE ring of fire
process.load('RecoMET.METFilters.eeNoiseFilter_cfi')
process.eeNoiseFilter.taggingMode = cms.bool(True)
# Inconsistent muon pf candidate filter
process.load('RecoMET.METFilters.inconsistentMuonPFCandidateFilter_cfi')
process.inconsistentMuonPFCandidateFilter.taggingMode = cms.bool(True)
# Greedy muon pf candidate filter
process.load('RecoMET.METFilters.greedyMuonPFCandidateFilter_cfi')
process.greedyMuonPFCandidateFilter.taggingMode = cms.bool(True)
# The ECAL laser correction filter (needs correct GT to work)
process.load('RecoMET.METFilters.ecalLaserCorrFilter_cfi')
process.ecalLaserCorrFilter.taggingMode = cms.bool(True)
# The tracking POG filters
process.load('RecoMET.METFilters.trackingPOGFilters_cff')
process.manystripclus53X.taggedMode = cms.untracked.bool(True)
process.toomanystripclus53X.taggedMode = cms.untracked.bool(True)
process.logErrorTooManyClusters.taggedMode = cms.untracked.bool(True)
process.logErrorTooManyTripletsPairs.taggedMode = cms.untracked.bool(True)
process.logErrorTooManySeeds.taggedMode = cms.untracked.bool(True)
#Add up all MET filters
process.metFiltersSequence = cms.Sequence(
process.HBHENoiseFilterResultProducer +
process.hcalLaserEventFilter +
process.EcalDeadCellTriggerPrimitiveFilter +
process.EcalDeadCellBoundaryEnergyFilter +
process.trackingFailureFilter +
process.eeBadScFilter +
process.eeNoiseFilter +
process.inconsistentMuonPFCandidateFilter +
process.greedyMuonPFCandidateFilter +
process.ecalLaserCorrFilter +
process.manystripclus53X +
process.toomanystripclus53X +
process.logErrorTooManyClusters +
process.logErrorTooManyTripletsPairs +
process.logErrorTooManySeeds
)
######################
### NoPU & MVA MET ###
######################
if runNoPUMVAMetSequence:
from RecoJets.JetProducers.PileupJetIDParams_cfi import JetIdParams
from JetMETCorrections.Configuration.JetCorrectionProducers_cff import ak5PFJetsL1
if isRealData:
process.ak5PFJetsL123Corrected = ak5PFJetsL1.clone(correctors = ['ak5PFL1FastL2L3Residual'])
else:
process.ak5PFJetsL123Corrected = ak5PFJetsL1.clone(correctors = ['ak5PFL1FastL2L3'])
process.recoPuJetIdCorrected = pileupJetIdProducer.clone(
jets = cms.InputTag("ak5PFJetsL123Corrected"),
algos = cms.VPSet(full_5x, cutbased),
residualsTxt = cms.FileInPath("RecoJets/JetProducers/data/mva_JetID_v1.weights.xml") # not used, but has to point to an existing file
)
# NoPU MET
from JetMETCorrections.METPUSubtraction.noPileUpPFMET_cff import noPileUpPFMEtData, noPileUpPFMEt
process.pfNoPileUpMetData = noPileUpPFMEtData.clone(
srcJets = cms.InputTag('ak5PFJetsL123Corrected'),
srcJetIds = cms.InputTag('recoPuJetIdCorrected', 'fullId'),
srcPFCandToVertexAssociations = cms.InputTag('pfCandidateToVertexAssoc'),
srcHardScatterVertex = cms.InputTag('primaryVertex')
)
process.pfNoPileUpMet = noPileUpPFMEt.clone(
srcMVAMEtData = cms.InputTag('pfNoPileUpMetData'),
srcLeptons = cms.VInputTag(
'pfIsolatedPhotons',
'pfIsolatedMuons',
'pfIsolatedElectrons'
),
srcType0Correction = cms.InputTag('pfType0MetCorrection'),
saveInputs = cms.bool(False)
)
# MVA MET
from JetMETCorrections.METPUSubtraction.mvaPFMET_cff import pfMEtMVA
process.pfMVAMet = pfMEtMVA.clone(
srcCorrJets = cms.InputTag('ak5PFJetsL123Corrected'),
srcLeptons = cms.VInputTag(
'pfIsolatedPhotons',
'pfIsolatedMuons',
'pfIsolatedElectrons'
)
)
process.noPUMVAMetSequence = cms.Sequence(
process.ak5PFJetsL123Corrected +
process.recoPuJetIdCorrected +
process.pfNoPileUpMetData +
process.pfNoPileUpMet +
process.pfMVAMet
)
else:
process.noPUMVAMetSequence = cms.Sequence() # placeholder
###################################################################
### Dataset-dependent sequence and event content configurations ###
###################################################################
if isRealData:
process.susyNtuplizer.metCollectionTags.remove('genMetTrue')
process.susyNtuplizer.jetFlavourMatchingTags = cms.PSet()
process.susyNtuplizer.gridParams = cms.vstring()
if isFastSim:
process.susyNtuplizer.muonCollectionTags = cms.vstring("muons")
process.susyNtuplizer.muonIdTags = cms.PSet()
process.metFiltersSequence.remove(process.HBHENoiseFilterResultProducer)
process.metFiltersSequence.remove(process.logErrorTooManyClusters)
process.metFiltersSequence.remove(process.logErrorTooManyTripletsPairs)
process.metFiltersSequence.remove(process.logErrorTooManySeeds)
# met filters with run = False will automatically be default = False
process.susyNtuplizer.metFilters.CSCBeamHalo.run = False
process.susyNtuplizer.metFilters.HcalNoise.run = False
process.susyNtuplizer.metFilters.LogErrorTooManyClusters.run = False
process.susyNtuplizer.metFilters.LogErrorTooManyTripletsPairs.run = False
process.susyNtuplizer.metFilters.LogErrorTooManySeeds.run = False
if is52x or isFastSim:
process.metFiltersSequence.remove(process.manystripclus53X)
process.metFiltersSequence.remove(process.toomanystripclus53X)
process.susyNtuplizer.metFilters.ManyStripClus53X.run = False
process.susyNtuplizer.metFilters.TooManyStripClus53X.run = False
if not runNoPUMVAMetSequence:
process.susyNtuplizer.metCollectionTags.remove('pfNoPileUpMet')
process.susyNtuplizer.metCollectionTags.remove('pfMVAMet')
if dataset in ['52xPrompt', '52x23May2012', '53xPromptC', '53x13July2012', '53x06Aug2012', '53x24Aug2012', '53x11Dec2013']:
process.susyNtuplizer.metFilters.HcalLaserEventList.default = True
if dataset in ['53xPromptD', '53x16Jan2013']:
process.susyNtuplizer.metFilters.HcalLaserOccupancy.default = True
if dataset == '53x22Jan2013':
process.susyNtuplizer.metFilters.HcalLaserRECOUserStep.default = True
else:
process.susyNtuplizer.metFilters.HcalLaserRECOUserStep.run = False
if dataset in ['53x13July2012', '53x24Aug2012']:
process.susyNtuplizer.metFilters.EcalLaserCorr.default = True
process.susyNtuplizer.storeLumiInfo = cms.bool(False)
if dataset == '53x22Jan2013':
process.correctedMetSequence.remove(process.pfMEtSysShiftCorr)
process.correctedMetSequence.remove(process.pfSysShiftCorrectedMet)
process.correctedMetSequence.remove(process.pfType01SysShiftCorrectedMet)
process.susyNtuplizer.metCollectionTags.remove('pfSysShiftCorrectedMet')
process.susyNtuplizer.metCollectionTags.remove('pfType01SysShiftCorrectedMet')
#####################
### Finalize path ###
#####################
process.standard_step = cms.Path(
process.hltHighLevel +
process.vertexSelectionSequence +
process.puRhoSequence +
process.pfBasedRecoSequence +
process.pfIsolationSequence +
process.pfCHSJetSequence +
process.correctedMetSequence +
process.eidMVASequence +
process.recoPuJetIdSequence +
process.JetFlavourMatchingSequence +
process.newJetTracksAssociator +
process.chsJetTracksAssociator +
process.newJetBtagging +
process.chsJetBtagging +
process.QGTaggingSequence +
process.metFiltersSequence
)
process.optional_step = cms.Path(
process.hltHighLevel +
process.noPUMVAMetSequence
)
process.ntuplizer_step = cms.Path(
process.hltHighLevel +
process.susyNtuplizer
)
process.schedule = cms.Schedule(process.standard_step,process.optional_step,process.ntuplizer_step)
return process
process.logErrorTooManyClusters + process.logErrorTooManyTripletsPairs +
process.logErrorTooManySeeds)
######################
### NoPU & MVA MET ###
######################
####### SET THIS TO TRUE TO RUN THE SEQUENCE (CHECK THE README FILE FOR INSTRUCTIONS) ########
runNoPUMVAMetSequence = True
##############################################################################################
if runNoPUMVAMetSequence:
from RecoJets.JetProducers.PileupJetIDParams_cfi import JetIdParams
from JetMETCorrections.Configuration.JetCorrectionProducers_cff import ak5PFJetsL1
if isRealData:
process.ak5PFJetsL123Corrected = ak5PFJetsL1.clone(
correctors=['ak5PFL1FastL2L3Residual'])
else:
process.ak5PFJetsL123Corrected = ak5PFJetsL1.clone(
correctors=['ak5PFL1FastL2L3'])
process.recoPuJetIdCorrected = pileupJetIdProducer.clone(
jets=cms.InputTag("ak5PFJetsL123Corrected"),
algos=cms.VPSet(full_5x, cutbased),
residualsTxt=cms.FileInPath(
"RecoJets/JetProducers/data/mva_JetID_v1.weights.xml"
) # not used, but has to point to an existing file
)
# NoPU MET
from JetMETCorrections.METPUSubtraction.noPileUpPFMET_cff import noPileUpPFMEtData, noPileUpPFMEt
process.pfNoPileUpMetData = noPileUpPFMEtData.clone(