CMS 3D CMS Logo

/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/Validation/RecoJets/analysis/allJets_validation_cfg.py

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 #from JetMETCorrections.Configuration.JetCorrectionsRecord_cfi import *
00004 #from RecoJets.Configuration.RecoJetAssociations_cff import *
00005 
00006 process = cms.Process("JETVALIDATION")
00007 
00008 process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
00009 process.GlobalTag.globaltag = 'START42_V17::All'
00010 #process.GlobalTag.globaltag = 'MC_38Y_V14::All'
00011 
00012 #process.load("Configuration.StandardSequences.Services_cff")
00013 #process.load("Configuration.StandardSequences.Simulation_cff")
00014 #process.load("Configuration.StandardSequences.MixingNoPileUp_cff")
00015 #process.load("Configuration.StandardSequences.VtxSmearedGauss_cff")
00016 #process.load("Configuration.StandardSequences.Geometry_cff")
00017 #process.load("Configuration.StandardSequences.MagneticField_cff")
00018 process.load("FWCore.MessageLogger.MessageLogger_cfi")
00019 
00020 from JetMETCorrections.Configuration.JetCorrectionServicesAllAlgos_cff import *
00021 #
00022 #
00023 # DQM
00024 #
00025 process.load("DQMServices.Core.DQM_cfg")
00026 
00027 # check # of bins
00028 #process.load("DQMServices.Components.DQMStoreStats_cfi")
00029 
00030 #process.load("JetMETCorrections.Configuration.JetPlusTrackCorrections_cff")
00031 #process.load("JetMETCorrections.Configuration.ZSPJetCorrections332_cff")
00032 process.load("JetMETCorrections.Configuration.DefaultJEC_cff")
00033 #process.ak5L1JPTOffset.era = 'Jec11V0'
00034 #process.ak5L1JPTOffset.useCondDB = False
00035 
00036 process.ak5L1JPTOffset.offsetService = cms.string('')
00037 
00038 process.load('RecoJets.Configuration.RecoPFJets_cff')
00039 process.kt6PFJets.doRhoFastjet = True
00040 process.ak5PFJets.doAreaFastjet = True
00041 
00042 # Validation module
00043 process.load("Validation.RecoJets.JetValidation_cff")
00044 
00045 process.maxEvents = cms.untracked.PSet(
00046        input = cms.untracked.int32(-1)
00047 )
00048 
00049 process.source = cms.Source("PoolSource",
00050 #    debugFlag = cms.untracked.bool(True),
00051 #    debugVebosity = cms.untracked.uint32(0),
00052 
00053     fileNames = cms.untracked.vstring(
00054 
00055 'file:/afs/cern.ch/user/k/kovitang/scratch0/60A53D23-A1BB-E011-987B-001A92971B90.root'
00056     )
00057 
00058 )
00059 
00060 process.fileSaver = cms.EDAnalyzer("JetFileSaver",
00061                         OutputFile = cms.untracked.string('histo.root')
00062 )
00063 
00064 ## Test for corrected jets - available only for 
00065 #process.prefer("L2L3CorJetIC5Calo")
00066 
00067 #process.L2L3CorJetIcone5 = cms.EDProducer("CaloJetCorrectionProducer",
00068 #    src = cms.InputTag("iterativeCone5CaloJets"),
00069 #    correctors = cms.vstring('L2L3JetCorrectorIC5Calo')
00070 #)
00071 
00072 
00073 ## AK5 Corrected jets
00074 #process.JetAnalyzerAK5Cor = cms.EDAnalyzer("CaloJetTester",
00075 #    src = cms.InputTag('ak5CaloJetsL2L3'),
00076 #    JetCorrectionService = cms.string('ak5CaloJetsL2L3'),
00077     ##src = cms.InputTag("L2L3CorJetAK5Calo"),
00078 #    srcGen = cms.InputTag("ak5GenJets"),
00079 ##    TurnOnEverything = cms.untracked.string('yes'),
00080 ##    TurnOnEverything = cms.untracked.string('no'),
00081 ##    outputFile = cms.untracked.string('histo.root'),
00082 ##    outputFile = cms.untracked.string('test.root'),
00083 #    genEnergyFractionThreshold = cms.double(0.05),
00084 #    genPtThreshold = cms.double(1.0),
00085 #    RThreshold = cms.double(0.3),
00086 #    reverseEnergyFractionThreshold = cms.double(0.5)
00087 #)
00088 
00089 
00090 ### IC5 JPT jets
00091 #JetAnalyzerIC5JPT = cms.EDFilter("CaloJetTester",
00092 #    src = cms.InputTag("ic5JPTJetsL2L3"),
00093 #    srcGen = cms.InputTag("iterativeCone5GenJets"),
00094 ##    TurnOnEverything = cms.untracked.string('yes'),
00095 ##    TurnOnEverything = cms.untracked.string('no'),
00096 ##    outputFile = cms.untracked.string('histo.root'),
00097 ##    outputFile = cms.untracked.string('test.root'),
00098 #    genEnergyFractionThreshold = cms.double(0.05),
00099 #    genPtThreshold = cms.double(1.0),
00100 #    RThreshold = cms.double(0.3),
00101 #    reverseEnergyFractionThreshold = cms.double(0.5)
00102 #)
00103 
00104 ### AntiKt5 JPT jets
00105 #JetAnalyzerAk5JPT = cms.EDFilter("CaloJetTester",
00106 #    src = cms.InputTag("ak5JPTJetsL2L3"),
00107 #    srcGen = cms.InputTag("ak5GenJets"),
00108 ##    TurnOnEverything = cms.untracked.string('yes'),
00109 ##    TurnOnEverything = cms.untracked.string('no'),
00110 ##    outputFile = cms.untracked.string('histo.root'),
00111 ##    outputFile = cms.untracked.string('test.root'),
00112 #    genEnergyFractionThreshold = cms.double(0.05),
00113 #    genPtThreshold = cms.double(1.0),
00114 #    RThreshold = cms.double(0.3),
00115 #    reverseEnergyFractionThreshold = cms.double(0.5)
00116 #)
00117 
00118 process.p1 = cms.Path(process.fileSaver
00119                       *process.kt6PFJets 
00120                       #* process.ak5PFJets 
00121                       #--- Non-Standard sequence (that involve Producers)
00122                       #*process.ak5CaloJetsL2L3
00123  #                     *process.ZSPJetCorrectionsIcone5
00124  #                     *process.ZSPJetCorrectionsAntiKt5
00125  #                     *process.JetPlusTrackCorrectionsIcone5
00126  #                     *process.JetPlusTrackCorrectionsAntiKt5
00127 #                      *process.JetAnalyzerAK5Cor
00128 #                      *process.JetAnalyzerIC5JPT
00129 #                      *process.JetAnalyzerAk5JPT
00130                       #--- Standard sequence
00131                       *process.JetValidation
00132                       #--- DQM stats module
00133 #                      *process.dqmStoreStats
00134 )
00135 process.DQM.collectorHost = ''
00136