Ejemplo n.º 1
0
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)
Ejemplo n.º 2
0
# 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=(