test/miniAOD/example_jet_and_met.py

0001 # THIS CONFIGURATION IS BROKEN. SINCE 2015 Geometry_cff has been deleted
0002 # and it is a fatal error to load it. And because of this I did not bother
0003 # to convert it to use tasks in 2017 when Tasks where implemented for unscheduled
0004 # mode (or remove the allowUnscheduled flag which no longer does anything).
0005 # Modules which are supposed to run unscheduled will not run.  Someone should
0006 # probably either fix or delete this ...
0007
0008 import FWCore.ParameterSet.Config as cms
0009
0010 process = cms.Process("S2")
0011 process.task = cms.Task()
0012
0013 process.source = cms.Source("PoolSource",
0014     fileNames = cms.untracked.vstring("file:patTuple_mini.root")
0015 )
0016 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0017
0018 from RecoJets.JetProducers.ak5PFJets_cfi import ak5PFJets
0019 from RecoJets.JetProducers.ak5GenJets_cfi import ak5GenJets
0020 from RecoMET.METProducers.pfMet_cfi import pfMet
0021
0022 process.load("CommonTools.ParticleFlow.pfCHS_cff")
0023 task.add(process.pfCHS)
0024
0025 process.ak5PFJets = ak5PFJets.clone(src = 'packedPFCandidates', doAreaFastjet = True) # no idea while doArea is false by default, but it's True in RECO so we have to set it
0026 process.ak5PFJetsCHS = ak5PFJets.clone(src = 'pfCHS', doAreaFastjet = True) # no idea while doArea is false by default, but it's True in RECO so we have to set it
0027 process.ak5GenJets = ak5GenJets.clone(src = 'packedGenParticles')
0028 process.pfMet = pfMet.clone(src = "packedPFCandidates")
0029 process.pfMet.calculateSignificance = False # this can't be easily implemented on packed PF candidates at the moment
0030
0031 process.load("PhysicsTools.PatAlgos.producersLayer1.patCandidates_cff")
0032 process.load("Configuration.EventContent.EventContent_cff")
0033 process.load('Configuration.StandardSequences.GeometryDB_cff')
0034 process.load('Configuration.StandardSequences.MagneticField_38T_cff')
0035 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0036 process.GlobalTag.globaltag = cms.string( autoCond[ 'phase1_2022_realistic' ] )
0037
0038
0039 process.patJets.addJetCharge   = False
0040 process.patJets.addBTagInfo    = True
0041 process.patJets.getJetMCFlavour = False
0042 process.patJets.addAssociatedTracks = False
0043 process.patJetPartonMatch.matched = "prunedGenParticles"
0044 process.patJetCorrFactors.primaryVertices = "offlineSlimmedPrimaryVertices"
0045 process.patMETs.addGenMET = False # There's no point in recalculating this, and we can't remake it since we don't have genParticles beyond |eta|=5
0046
0047 process.load('RecoBTag.Configuration.RecoBTag_cff')
0048 process.load('RecoJets.Configuration.RecoJetAssociations_cff')
0049
0050 process.load('PhysicsTools.PatAlgos.slimming.unpackedTracksAndVertices_cfi')
0051
0052 process.ak5JetTracksAssociatorAtVertexPF.jets = cms.InputTag("ak5PFJetsCHS")
0053 process.ak5JetTracksAssociatorAtVertexPF.tracks = cms.InputTag("unpackedTracksAndVertices")
0054 process.impactParameterTagInfos.primaryVertex = cms.InputTag("unpackedTracksAndVertices")
0055 process.inclusiveSecondaryVertexFinderTagInfos.extSVCollection = cms.InputTag("unpackedTracksAndVertices","secondary","")
0056
0057
0058 process.p = cms.Path(
0059     process.patJets + process.patMETs + process.inclusiveSecondaryVertexFinderTagInfos
0060 )
0061
0062 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0063 process.MessageLogger.cerr.FwkReport.reportEvery = 10
0064 process.MessageLogger.suppressWarning = cms.untracked.vstring('ecalLaserCorrFilter','manystripclus53X','toomanystripclus53X')
0065 process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0066 process.options.allowUnscheduled = cms.untracked.bool(True)
0067
0068 process.OUT = cms.OutputModule("PoolOutputModule",
0069     fileName = cms.untracked.string('test.root'),
0070     outputCommands = cms.untracked.vstring(['drop *','keep patJets_patJets_*_*','keep *_*_*_PAT','keep recoTracks_unp*_*_*','keep recoVertexs_unp*_*_*'])
0071 )
0072 process.endpath= cms.EndPath(process.OUT, process.task)
0073