PatExamples/test/analyzePatBJetTracks_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("Test")
0004
0005 # initialize MessageLogger and output report
0006 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0007 process.MessageLogger.cerr.threshold = 'INFO'
0008 process.MessageLogger.cerr.INFO = cms.untracked.PSet(
0009     default          = cms.untracked.PSet( limit = cms.untracked.int32(0)  ),
0010     PATSummaryTables = cms.untracked.PSet( limit = cms.untracked.int32(-1) )
0011 )
0012 process.options   = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
0013
0014 process.load("Configuration.StandardSequences.GeometryDB_cff")
0015 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0016 process.load("Configuration.StandardSequences.MagneticField_cff")
0017
0018 process.GlobalTag.globaltag = cms.string( autoCond[ 'phase1_2022_realistic' ] )
0019 # produce PAT Layer 1
0020 process.load("PhysicsTools.PatAlgos.patSequences_cff")
0021 # switch old trigger matching off
0022 from PhysicsTools.PatAlgos.tools.trigTools import switchOffTriggerMatchingOld
0023 switchOffTriggerMatchingOld( process )
0024
0025 # source
0026 process.source = cms.Source("PoolSource",
0027     fileNames = cms.untracked.vstring(
0028         '/store/mc/Fall08/TTJets-madgraph/GEN-SIM-RECO/IDEAL_V11_redigi_v10/0000/06FC3959-4DFC-DD11-B504-00E08178C091.root',
0029     )
0030 )
0031 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0032
0033 process.TFileService = cms.Service("TFileService",
0034     fileName = cms.string("analyzePatBJetTracks.root")
0035 )
0036
0037 process.analyzeBJetTracks = cms.EDAnalyzer("PatBJetTrackAnalyzer",
0038     # input collections
0039     primaryVertices = cms.InputTag("offlinePrimaryVertices"),
0040     beamSpot = cms.InputTag("offlineBeamSpot"),
0041     tracks = cms.InputTag("generalTracks"),
0042     jets = cms.InputTag("selectedLayer1Jets"),
0043
0044     #jet cuts
0045     jetPtCut = cms.double(30),
0046     jetEtaCut = cms.double(2.4),
0047
0048     #track cuts
0049     maxDeltaR = cms.double(0.5),
0050     minPt = cms.double(1.0),
0051     minPixelHits = cms.uint32(2),
0052     minTotalHits = cms.uint32(8),
0053
0054     # selecting the second-highest signed IP (i.e. "TrackCountingHighEff")
0055     nThTrack = cms.uint32(2)
0056 )
0057
0058 process.p = cms.Path(
0059     process.patDefaultSequence *
0060     process.analyzeBJetTracks
0061 )