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     IgnoreCompletely = cms.untracked.vstring(),
0036     Rethrow = cms.untracked.vstring(),
0037     TryToContinue = cms.untracked.vstring(),
0038     allowUnscheduled = cms.obsolete.untracked.bool,
0039     canDeleteEarly = cms.untracked.vstring(),
0040     emptyRunLumiMode = cms.obsolete.untracked.string,
0041     eventSetup = cms.untracked.PSet(
0042         forceNumberOfConcurrentIOVs = cms.untracked.PSet(
0043
0044         ),
0045         numberOfConcurrentIOVs = cms.untracked.uint32(1)
0046     ),
0047     fileMode = cms.untracked.string('FULLMERGE'),
0048     forceEventSetupCacheClearOnNewRun = cms.untracked.bool(False),
0049     makeTriggerResults = cms.obsolete.untracked.bool,
0050     numberOfConcurrentLuminosityBlocks = cms.untracked.uint32(1),
0051     numberOfConcurrentRuns = cms.untracked.uint32(1),
0052     numberOfStreams = cms.untracked.uint32(0),
0053     numberOfThreads = cms.untracked.uint32(1),
0054     printDependencies = cms.untracked.bool(False),
0055     sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0056     throwIfIllegalParameter = cms.untracked.bool(True),
0057     wantSummary = cms.untracked.bool(False)
0058 )
0059
0060 # Production Info
0061 process.configurationMetadata = cms.untracked.PSet(
0062     annotation = cms.untracked.string('step5 nevts:10'),
0063     name = cms.untracked.string('Applications'),
0064     version = cms.untracked.string('$Revision: 1.19 $')
0065 )
0066
0067 # Output definition
0068
0069 process.MINIAODSIMoutput = cms.OutputModule("PoolOutputModule",
0070     compressionAlgorithm = cms.untracked.string('LZMA'),
0071     compressionLevel = cms.untracked.int32(9),
0072     dataset = cms.untracked.PSet(
0073         dataTier = cms.untracked.string('MINIAODSIM'),
0074         filterName = cms.untracked.string('')
0075     ),
0076     fileName = cms.untracked.string('file:updatedMINIAODSIM.root'),
0077     outputCommands = process.MINIAODSIMEventContent.outputCommands
0078 )
0079
0080 # Additional output definition
0081
0082 # Other statements
0083 from Configuration.AlCa.GlobalTag import GlobalTag
0084 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0085
0086 # Rerun PUPPI MET and ak4 jets (but not ak8)
0087 from PhysicsTools.PatAlgos.tools.helpers import getPatAlgosToolsTask, addToProcessAndTask
0088 task = getPatAlgosToolsTask(process)
0089 from PhysicsTools.PatAlgos.slimming.puppiForMET_cff import makePuppiesFromMiniAOD
0090 makePuppiesFromMiniAOD(process,True)
0091 process.puppi.useExistingWeights = False
0092 process.puppiNoLep.useExistingWeights = False
0093 from PhysicsTools.PatUtils.tools.runMETCorrectionsAndUncertainties import runMetCorAndUncFromMiniAOD
0094 runMetCorAndUncFromMiniAOD(process,isData=False,metType="Puppi",postfix="Puppi",jetFlavor="AK4PFPuppi",recoMetFromPFCs=True)
0095 from PhysicsTools.PatAlgos.patPuppiJetSpecificProducer_cfi import patPuppiJetSpecificProducer
0096 addToProcessAndTask('patPuppiJetSpecificProducer', patPuppiJetSpecificProducer.clone(src=cms.InputTag("patJetsPuppi")), process, task)
0097 from PhysicsTools.PatAlgos.tools.jetTools import updateJetCollection
0098 updateJetCollection(
0099    process,
0100    labelName = 'PuppiJetSpecific',
0101    jetSource = cms.InputTag('patJetsPuppi'),
0102 )
0103 process.updatedPatJetsPuppiJetSpecific.userData.userFloats.src = ['patPuppiJetSpecificProducer:puppiMultiplicity', 'patPuppiJetSpecificProducer:neutralPuppiMultiplicity', 'patPuppiJetSpecificProducer:neutralHadronPuppiMultiplicity', 'patPuppiJetSpecificProducer:photonPuppiMultiplicity', 'patPuppiJetSpecificProducer:HFHadronPuppiMultiplicity', 'patPuppiJetSpecificProducer:HFEMPuppiMultiplicity' ]
0104 addToProcessAndTask('slimmedJetsPuppi', process.updatedPatJetsPuppiJetSpecific.clone(), process, task)
0105 del process.updatedPatJetsPuppiJetSpecific
0106 process.puppiSequence = cms.Sequence(process.puppiMETSequence+process.fullPatMetSequencePuppi+process.patPuppiJetSpecificProducer+process.slimmedJetsPuppi)
0107
0108 # Example how to not use vertex fit for track-vertex-association, but only dz
0109 #process.packedPrimaryVertexAssociationJME.assignment.useVertexFit = False
0110 #process.packedPrimaryVertexAssociationJME.assignment.maxDzSigForPrimaryAssignment = 5.0
0111 #process.packedPrimaryVertexAssociationJME.assignment.maxDzForPrimaryAssignment = 0.1
0112 #process.packedPrimaryVertexAssociationJME.assignment.maxDzErrorForPrimaryAssignment = 0.05
0113 #process.packedPrimaryVertexAssociationJME.assignment.OnlyUseFirstDz = False
0114 #process.pfCHS.vertexAssociationQuality = 6
0115
0116 # Path and EndPath definitions
0117 process.puppi_step = cms.Path(process.puppiSequence)
0118 process.endjob_step = cms.EndPath(process.endOfProcess)
0119 process.MINIAODSIMoutput_step = cms.EndPath(process.MINIAODSIMoutput)
0120
0121 # Schedule definition
0122 process.schedule = cms.Schedule(process.puppi_step,process.endjob_step,process.MINIAODSIMoutput_step)
0123 from PhysicsTools.PatAlgos.tools.helpers import associatePatAlgosToolsTask
0124 associatePatAlgosToolsTask(process)
0125
0126 #Setup FWK for multithreaded
0127 process.options.numberOfThreads=cms.untracked.uint32(2)
0128 process.options.numberOfStreams=cms.untracked.uint32(0)
0129 process.options.numberOfConcurrentLuminosityBlocks=cms.untracked.uint32(1)
0130
0131 # customisation of the process.
0132
0133 # End of customisation functions
0134
0135 # Customisation from command line
0136
0137 # Add early deletion of temporary data products to reduce peak memory need
0138 from Configuration.StandardSequences.earlyDeleteSettings_cff import customiseEarlyDelete
0139 process = customiseEarlyDelete(process)
0140 # End adding early deletion