process.printList = cms.EDAnalyzer( "ParticleListDrawer", src = cms.InputTag( "genParticles" ), maxEventsToPrint = cms.untracked.int32( 1 ) ) process.load( "SimGeneral.HepPDTESSource.pythiapdt_cfi") process.load("PhysicsTools.RecoAlgos.allTrackCandidates_cfi") # input pat sequences process.load("PhysicsTools.PatAlgos.patLayer0_cff") process.load("PhysicsTools.PatAlgos.patLayer1_cff") ## Necessary fixes to run 2.2.X on 2.1.X data from PhysicsTools.PatAlgos.tools.cmsswVersionTools import run22XonSummer08AODSIM run22XonSummer08AODSIM(process) # input pat analyzer sequence process.load("PhysicsTools.StarterKit.PatAnalyzerKit_cfi") # talk to TFileService for output histograms process.TFileService = cms.Service("TFileService", fileName = cms.string('PatAnalyzerKitHistos.root') ) # define path 'p': PAT Layer 0, PAT Layer 1, and the analyzer process.p = cms.Path(process.printList* process.allTrackCandidates* process.patLayer0* process.patLayer1* process.patAnalyzerKit)
# turn off sum-et dependent stuff. process.caTopJetsProducer.ptBins = cms.vdouble(0, 10e9) process.caTopJetsProducer.rBins = cms.vdouble(0.8, 0.8) process.caTopJetsProducer.ptFracBins = cms.vdouble(0.05, 0.05) process.caTopJetsProducer.nCellBins = cms.vint32(1, 1) print "About to input pat sequences" # input pat sequences process.load("PhysicsTools.PatAlgos.patLayer0_cff") process.load("PhysicsTools.PatAlgos.patLayer1_cff") ## Necessary fixes to run 2.2.X on 2.1.X data from PhysicsTools.PatAlgos.tools.cmsswVersionTools import run22XonSummer08AODSIM run22XonSummer08AODSIM(process) # switch jet collection to our juets from PhysicsTools.PatAlgos.tools.jetTools import * print "About to switch jet collection" ## ==== Example with CaloJets switchJetCollection( process, 'caTopJetsProducer', # Jet collection; must be already in the event when patLayer0 sequence is executed layers=[0, 1], # If you're not runnint patLayer1, set 'layers=[0]' runCleaner="BasicJet", # =None if not to clean doJTA=True, # Run Jet-Track association & JetCharge doBTagging=True, # Run b-tagging jetCorrLabel=(