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