|
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
90
91
92
93
94
|
# 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_data']
process.load('FWCore.MessageService.MessageLogger_cfi')
process.MessageLogger.GammaJetAnalysis=dict()
process.MessageLogger.cerr.FwkReport.reportEvery=cms.untracked.int32(1000)
#load the analyzer
process.load('Calibration.HcalCalibAlgos.gammaJetAnalysis_cfi')
# load energy corrector
process.load('JetMETCorrections.Configuration.CorrectedJetProducers_cff')
# run over files
process.GammaJetAnalysis.rootHistFilename = cms.string('PhoJet_tree_CHS_data2012_noGJetProd.root')
process.GammaJetAnalysis.doPFJets = cms.bool(True)
process.GammaJetAnalysis.doGenJets = cms.bool(False)
process.TFileService = cms.Service("TFileService",
fileName = cms.string('PhoJet_tree_CHS_data2012_noGJetProd.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()
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(
'file:/tmp/andriusj/Run2012A_Photon_22Jan2013-002618943913.root'
)
)
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(False) )
# Load jets and pfNoPileUP
process.load('RecoJets.Configuration.RecoPFJets_cff')
process.load('CommonTools.ParticleFlow.PF2PAT_cff')
process.load("CommonTools.ParticleFlow.pfNoPileUp_cff")
process.seq_ak4PFCHS= cms.Sequence( process.particleFlowPtrs *
process.pfNoPileUpJMESequence *
process.ak4PFJetsCHS )
# adapt input collections
process.GammaJetAnalysis.photonCollName= cms.string("photons")
process.GammaJetAnalysis.electronCollName= cms.string("gsfElectrons")
process.GammaJetAnalysis.photonIdLooseName= cms.InputTag("PhotonIDProd","PhotonCutBasedIDLoose")
process.GammaJetAnalysis.photonIdTightName= cms.InputTag("PhotonIDProd","PhotonCutBasedIDTight")
# 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(0)
process.GammaJetAnalysis.doGenJets = cms.bool(False)
process.GammaJetAnalysis.debug = cms.untracked.int32(0)
process.p = cms.Path(
process.seq_ak4PFCHS *
process.GammaJetAnalysis, process.ak4PFCHSL2L3CorrectorTask
)
|