PatExamples/test/analyzePatAnalysisTasks_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("Test")
0004
0005 process.source = cms.Source("PoolSource",
0006   fileNames = cms.untracked.vstring(
0007     "file:patTuple.root"
0008   )
0009 )
0010 process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
0011
0012 process.MessageLogger = cms.Service("MessageLogger")
0013
0014 #################
0015 #               #
0016 # EXERCISE 1    #
0017 #               #
0018 #################
0019
0020 process.jecAnalyzer = cms.EDAnalyzer("WrappedEDAnalysisTasksAnalyzerJEC",
0021   Jets = cms.InputTag("cleanPatJets"),
0022   jecLevel=cms.string("SomeWrongName"),
0023 #  jecLevel=cms.string("L3Absolute"),
0024   patJetCorrFactors= cms.string('SomeWrongName'),
0025 #  patJetCorrFactors= cms.string('patJetCorrFactors'),
0026   help=cms.bool(False)
0027 )
0028
0029 process.p = cms.Path(process.jecAnalyzer)
0030
0031
0032 #################
0033 #               #
0034 # EXERCISE 2    #
0035 #               #
0036 #################
0037
0038 process.TFileService = cms.Service("TFileService",
0039   fileName = cms.string('TFileServiceOutput.root')
0040 )
0041
0042
0043 process.jecAnalyzerRel=process.jecAnalyzer.clone(jecLevel="L2Relative")
0044 process.jecAnalyzerNon=process.jecAnalyzer.clone(jecLevel="Uncorrected")
0045 process.p.replace(process.jecAnalyzer, process.jecAnalyzer * process.jecAnalyzerRel * process.jecAnalyzerNon)
0046
0047
0048 #################
0049 #               #
0050 # EXERCISE 3.2  #
0051 #               #
0052 #################
0053 ## Geometry and Detector Conditions (needed for a Scaling up and down the JEC)
0054 process.load("Configuration.StandardSequences.GeometryDB_cff")
0055 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0056 from Configuration.AlCa.autoCond import autoCond
0057 process.GlobalTag.globaltag = cms.string( autoCond[ 'startup' ] )
0058 process.load("Configuration.StandardSequences.MagneticField_cff")
0059
0060 # load the standard PAT config
0061 process.load("PhysicsTools.PatAlgos.patSequences_cff")
0062 from PhysicsTools.PatUtils.tools.metUncertaintyTools import runMEtUncertainties
0063 runMEtUncertainties(process)
0064
0065 #process.shiftedPatJets=process.shiftedPatJetsEnUp.clone(src = cms.InputTag("cleanPatJets"),shiftBy=cms.double(1))
0066
0067 #process.jecAnalyzerEnUp=process.jecAnalyzer.clone(Jets = cms.InputTag("shiftedPatJets"))
0068 #process.p.replace(process.jecAnalyzer,  process.shiftedPatJets * process.jecAnalyzerEnUp)
0069
0070 #################
0071 #               #
0072 # EXERCISE 3.3  #
0073 #               #
0074 #################
0075
0076 #process.btagAnalyzer = cms.EDAnalyzer("WrappedEDAnalysisTasksAnalyzerBTag",
0077 #  Jets = cms.InputTag("cleanPatJets"),
0078 #   bTagAlgo=cms.string('trackCountingHighEffBJetTags'),
0079 #   bins=cms.uint32(100),
0080 #   lowerbin=cms.double(0.),
0081 #   upperbin=cms.double(10.)
0082 #)
0083
0084 #process.btagAnalyzerTCHP=process.btagAnalyzer.clone(bTagAlgo="trackCountingHighPurBJetTags")
0085
0086 #process.p.replace(process.jecAnalyzer, process.jecAnalyzer * process.btagAnalyzer * process.btagAnalyzerTCHP)
0087
0088
0089
0090