PileupAlgos/test/recomputePuppiOnMiniAOD.py

0001 # Auto generated configuration file
0002 # using:
0003 # Revision: 1.19
0004 # Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
0005 # with command line options: step5 --conditions auto:run2_mc --nThreads 2 -n 10 --era Run2_2016 --eventcontent MINIAODSIM --filein file:step3_inMINIAODSIM.root -s NANO --datatier NANOAODSIM --mc --fileout file:step5.root
0006 import FWCore.ParameterSet.Config as cms
0007
0008 from Configuration.Eras.Era_Run2_2016_cff import Run2_2016
0009
0010 process = cms.Process('NANO',Run2_2016)
0011
0012 # import of standard configurations
0013 process.load('Configuration.StandardSequences.Services_cff')
0014 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0015 process.load('FWCore.MessageService.MessageLogger_cfi')
0016 process.load('Configuration.EventContent.EventContent_cff')
0017 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0018 process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
0019 process.load('Configuration.StandardSequences.MagneticField_cff')
0020 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0021 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0022
0023 process.maxEvents = cms.untracked.PSet(
0024     input = cms.untracked.int32(10),
0025     output = cms.optional.untracked.allowed(cms.int32,cms.PSet)
0026 )
0027
0028 # Input source
0029 process.source = cms.Source("PoolSource",
0030     fileNames = cms.untracked.vstring('file:/afs/cern.ch/user/h/hinzmann/stable_13TeV/jetmet/mtd/CMSSW_12_0_0_pre3/src/25202.0_TTbar_13+TTbar_13+DIGIUP15_PU25+RECOUP15_PU25+HARVESTUP15_PU25+NANOUP15_PU25/step3_inMINIAODSIM.root'),
0031     secondaryFileNames = cms.untracked.vstring()
0032 )
0033
0034 process.options = cms.untracked.PSet(
0035     FailPath = cms.untracked.vstring(),
0036     IgnoreCompletely = cms.untracked.vstring(),
0037     Rethrow = cms.untracked.vstring(),
0038     SkipEvent = cms.untracked.vstring(),
0039     allowUnscheduled = cms.obsolete.untracked.bool,
0040     canDeleteEarly = cms.untracked.vstring(),
0041     emptyRunLumiMode = cms.obsolete.untracked.string,
0042     eventSetup = cms.untracked.PSet(
0043         forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0044
0045         ),
0046         numberOfConcurrentIOVs = cms.untracked.uint32(1)
0047     ),
0048     fileMode = cms.untracked.string('FULLMERGE'),
0049     forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0050     makeTriggerResults = cms.obsolete.untracked.bool,
0051     numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(1),
0052     numberOfConcurrentRuns = cms.untracked.uint32(1),
0053     numberOfStreams = cms.untracked.uint32(0),
0054     numberOfThreads = cms.untracked.uint32(1),
0055     printDependencies = cms.untracked.bool(False),
0056     sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0057     throwIfIllegalParameter = cms.untracked.bool(True),
0058     wantSummary = cms.untracked.bool(False)
0059 )
0060
0061 # Production Info
0062 process.configurationMetadata = cms.untracked.PSet(
0063     annotation = cms.untracked.string('step5 nevts:10'),
0064     name = cms.untracked.string('Applications'),
0065     version = cms.untracked.string('$Revision: 1.19 $')
0066 )
0067
0068 # Output definition
0069
0070 process.MINIAODSIMoutput = cms.OutputModule("PoolOutputModule",
0071     compressionAlgorithm = cms.untracked.string('LZMA'),
0072     compressionLevel = cms.untracked.int32(9),
0073     dataset = cms.untracked.PSet(
0074         dataTier = cms.untracked.string('MINIAODSIM'),
0075         filterName = cms.untracked.string('')
0076     ),
0077     fileName = cms.untracked.string('file:updatedMINIAODSIM.root'),
0078     outputCommands = process.MINIAODSIMEventContent.outputCommands
0079 )
0080
0081 # Additional output definition
0082
0083 # Other statements
0084 from Configuration.AlCa.GlobalTag import GlobalTag
0085 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0086
0087 # Rerun PUPPI MET and ak4 jets (but not ak8)
0088 from PhysicsTools.PatAlgos.tools.helpers import getPatAlgosToolsTask, addToProcessAndTask
0089 task = getPatAlgosToolsTask(process)
0090 from PhysicsTools.PatAlgos.slimming.puppiForMET_cff import makePuppiesFromMiniAOD
0091 makePuppiesFromMiniAOD(process,True)
0092 process.puppi.useExistingWeights = False
0093 process.puppiNoLep.useExistingWeights = False
0094 from PhysicsTools.PatUtils.tools.runMETCorrectionsAndUncertainties import runMetCorAndUncFromMiniAOD
0095 runMetCorAndUncFromMiniAOD(process,isData=False,metType="Puppi",postfix="Puppi",jetFlavor="AK4PFPuppi",recoMetFromPFCs=True)
0096 from PhysicsTools.PatAlgos.patPuppiJetSpecificProducer_cfi import patPuppiJetSpecificProducer
0097 addToProcessAndTask('patPuppiJetSpecificProducer', patPuppiJetSpecificProducer.clone(src=cms.InputTag("patJetsPuppi")), process, task)
0098 from PhysicsTools.PatAlgos.tools.jetTools import updateJetCollection
0099 updateJetCollection(
0100    process,
0101    labelName = 'PuppiJetSpecific',
0102    jetSource = cms.InputTag('patJetsPuppi'),
0103 )
0104 process.updatedPatJetsPuppiJetSpecific.userData.userFloats.src = ['patPuppiJetSpecificProducer:puppiMultiplicity', 'patPuppiJetSpecificProducer:neutralPuppiMultiplicity', 'patPuppiJetSpecificProducer:neutralHadronPuppiMultiplicity', 'patPuppiJetSpecificProducer:photonPuppiMultiplicity', 'patPuppiJetSpecificProducer:HFHadronPuppiMultiplicity', 'patPuppiJetSpecificProducer:HFEMPuppiMultiplicity' ]
0105 addToProcessAndTask('slimmedJetsPuppi', process.updatedPatJetsPuppiJetSpecific.clone(), process, task)
0106 del process.updatedPatJetsPuppiJetSpecific
0107 process.puppiSequence = cms.Sequence(process.puppiMETSequence+process.fullPatMetSequencePuppi+process.patPuppiJetSpecificProducer+process.slimmedJetsPuppi)
0108
0109 # Example how to not use vertex fit for track-vertex-association, but only dz
0110 #process.packedPrimaryVertexAssociationJME.assignment.useVertexFit = False
0111 #process.packedPrimaryVertexAssociationJME.assignment.maxDzSigForPrimaryAssignment = 5.0
0112 #process.packedPrimaryVertexAssociationJME.assignment.maxDzForPrimaryAssignment = 0.1
0113 #process.packedPrimaryVertexAssociationJME.assignment.maxDzErrorForPrimaryAssignment = 0.05
0114 #process.packedPrimaryVertexAssociationJME.assignment.OnlyUseFirstDz = False
0115 #process.pfCHS.vertexAssociationQuality = 6
0116
0117 # Path and EndPath definitions
0118 process.puppi_step = cms.Path(process.puppiSequence)
0119 process.endjob_step = cms.EndPath(process.endOfProcess)
0120 process.MINIAODSIMoutput_step = cms.EndPath(process.MINIAODSIMoutput)
0121
0122 # Schedule definition
0123 process.schedule = cms.Schedule(process.puppi_step,process.endjob_step,process.MINIAODSIMoutput_step)
0124 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0125 associatePatAlgosToolsTask(process)
0126
0127 #Setup FWK for multithreaded
0128 process.options.numberOfThreads=cms.untracked.uint32(2)
0129 process.options.numberOfStreams=cms.untracked.uint32(0)
0130 process.options.numberOfConcurrentLuminosityBlocks=cms.untracked.uint32(1)
0131
0132 # customisation of the process.
0133
0134 # End of customisation functions
0135
0136 # Customisation from command line
0137
0138 # Add early deletion of temporary data products to reduce peak memory need
0139 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0140 process = customiseEarlyDelete(process)
0141 # End adding early deletion