|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
# This was migrated to use reco::JetCorrector in Nov 2022.
# This was not fully tested because this configuration
# before the migration already failed with multiple errors
# unrelated to JetCorrectors. At the least:
#
# Something unknown but unrelated to JetCorrectors in these
# three lines in gammaJetAnalysis_cfi.py causes this file to
# be unparseable by Python:
#
# from RecoJets.Configuration.RecoJets_cff import *
# from RecoJets.Configuration.RecoPFJets_cff import *
# from CommonTools.ParticleFlow.pfNoPileUp_cff import *
#
# The input file does not exist in a publicly available
# space. There may be other problems.
import FWCore.ParameterSet.Config as cms
process = cms.Process('ANALYSIS')
process.load('Configuration.StandardSequences.Services_cff')
# Specify IdealMagneticField ESSource (needed for CMSSW 730)
process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
from Configuration.AlCa.autoCond import autoCond
process.GlobalTag.globaltag=autoCond['run1_mc']
process.load('FWCore.MessageService.MessageLogger_cfi')
process.MessageLogger.GammaJetAnalysis = dict()
process.MessageLogger.cerr.FwkReport.reportEvery=cms.untracked.int32(1000)
#load the gammaJet analyzer
process.load('Calibration.HcalCalibAlgos.gammaJetAnalysis_cfi')
# needed for nonCHS
process.load('JetMETCorrections.Configuration.CorrectedJetProducers_cff')
# run over files
process.GammaJetAnalysis.rootHistFilename = cms.string('PhoJet_tree_CHS.root')
process.GammaJetAnalysis.doPFJets = cms.bool(True)
process.GammaJetAnalysis.doGenJets = cms.bool(False)
process.TFileService = cms.Service("TFileService",
fileName = cms.string('PhoJet_tree_CHS.root'))
# trigger names should not end with '_'
process.GammaJetAnalysis.photonTriggers = cms.vstring(
'HLT_Photon20_CaloIdVL_IsoL','HLT_Photon30_CaloIdVL_IsoL',
'HLT_Photon50_CaloIdVL_IsoL','HLT_Photon75_CaloIdVL_IsoL',
'HLT_Photon90_CaloIdVL_IsoL','HLT_Photon135',
'HLT_Photon150','HLT_Photon160')
# triggers for CMSSW 730
process.GammaJetAnalysis.photonTriggers += cms.vstring(
'HLT_Photon22', 'HLT_Photon30', 'HLT_Photon36',
'HLT_Photon50', 'HLT_Photon75',
'HLT_Photon90', 'HLT_Photon120', 'HLT_Photon175',
'HLT_Photon250_NoHE', 'HLT_Photon300_NoHE'
)
# to disable photonTriggers assign an empty vstring
#process.GammaJetAnalysis.photonTriggers = cms.vstring()
# a clone without CHS
process.GammaJetAnalysis_noCHS= process.GammaJetAnalysis.clone()
process.GammaJetAnalysis_noCHS.rootHistFilename = cms.string('PhoJet_tree_nonCHS.root')
# for 7XY use ak4* instead of ak5
process.GammaJetAnalysis_noCHS.pfJetCollName = cms.string('ak4PFJets')
process.GammaJetAnalysis_noCHS.JetCorrections = cms.InputTag("ak4PFL2L3Corrector")
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(
'file:../../HcalAlCaRecoProducers/test/gjet.root'
# '/store/relval/CMSSW_7_3_0/RelValPhotonJets_Pt_10_13/GEN-SIM-RECO/MCRUN2_73_V7-v1/00000/522CE329-7B81-E411-B6C3-0025905A6110.root',
# '/store/relval/CMSSW_7_3_0/RelValPhotonJets_Pt_10_13/GEN-SIM-RECO/MCRUN2_73_V7-v1/00000/5279D224-7B81-E411-BCAA-002618943930.root'
# '/store/relval/CMSSW_7_3_0/RelValPhotonJets_Pt_10_13/GEN-SIM-RECO/MCRUN2_73_V7-v1/00000/522CE329-7B81-E411-B6C3-0025905A6110.root'
# '/store/relval/CMSSW_7_4_0_pre6/RelValPhotonJets_Pt_10_13/GEN-SIM-RECO/MCRUN2_74_V1-v1/00000/6EC8FCC8-E2A8-E411-9506-002590596468.root'
)
)
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(False) )
# name of the process that used the GammaJetProd producer
process.GammaJetAnalysis.prodProcess = cms.untracked.string('MYGAMMAJET')
# specify 'workOnAOD=2' to apply tokens from GammaJetProd producer
process.GammaJetAnalysis.workOnAOD = cms.int32(2)
process.GammaJetAnalysis.doGenJets = cms.bool(False)
process.GammaJetAnalysis.debug = cms.untracked.int32(0)
process.p = cms.Path(
process.GammaJetAnalysis, process.ak4PFL2L3CorrectorTask, process.ak4PFCHSL2L3CorrectorTask
)
|