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
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
|
import FWCore.ParameterSet.Config as cms
process = cms.Process("hcalCalib")
process.load("Configuration.StandardSequences.GeometryHCAL_cff")
# no CMS input files
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(0) )
process.source = cms.Source("EmptySource")
process.hcalCalib = cms.EDAnalyzer("HcalCalibrator",
#-------- File containing the list of the input root files
inputFileList = cms.untracked.string("test/inputFiles_IsoTrkCalib.list"),
#-------- calibration type: ISO_TRACK or DI_JET
calibType = cms.untracked.string("ISO_TRACK"),
#-------- calibration method: L3, MATRIX_INV_OF_ETA_AVE, L3_AND_MTRX_INV
calibMethod = cms.untracked.string("L3_AND_MTRX_INV"),
# --------min target energy: ref jet (had) energy, or isotrack momentum
minTargetE = cms.untracked.double(40.0),
# ------- max target energy: ref jet (had) energy, or isotrack momentum
maxTargetE = cms.untracked.double(60.0),
# ------- min energy in the cell
minCellE = cms.untracked.double(0.0),
# ------- minimum e/p to accept an IsoTrack for calibration (not used for DiJEts)
minEOverP = cms.untracked.double(0.2),
# ------- maximum e/p to accept an IsoTrack for calibration (not used for DiJEts)
maxEOverP = cms.untracked.double(999.0),
# ------- MIP cut: maximum EM energy associated with the track (not used for DiJEts)
maxTrkEmE = cms.untracked.double(1.0),
#--------Maximum Et allowed for third jet in the event (affects only DiJets)
maxEtThirdJet = cms.untracked.double(5.0),
# ------- Minimum deltaPhi for the dijets in degrees: 0-180 (affects only DiJets)
minDPhiDiJets = cms.untracked.double(150.0),
# ------- logical flag to sum depths in HCAL
sumDepths = cms.untracked.bool(True),
# ------- logical flag to sum depths 1,2 in HB for towers 15,16.
# --------If sumDepths="true" this flag has no effect
sumSmallDepths = cms.untracked.bool(True),
# ------- logical flag to combine phi in HCAL
combinePhi = cms.untracked.bool(True),
# ------- cluster size in HB for isotracks: 3 or 5 (means 3x3, 5x5)
hbClusterSize = cms.untracked.int32(3),
# ------- cluster size in HE for isotracks: 3 or 5 (means 3x3, 5x5)
heClusterSize = cms.untracked.int32(3),
# -------- flag to use cone clustering -> overrides the above cluster sizes
# -------- and uses cone size as specified below
useConeClustering = cms.untracked.bool(True),
# -------- size of the cone (when useConeClustering=True)
maxConeDist = cms.untracked.double(26.2),
# ------- max ABS(iEta) used in the calibration: for matrix inversion sets the range for performing inversion
# ------- For all methods: controls the range of correction factors are saved in the output file
calibAbsIEtaMax = cms.untracked.int32(24),
# ------- max ABS(iEta) used in the calibration: for matrix inversion sets the range for performing inversion
# ------- For all methods: controls the range of correction factors are saved in the output file
calibAbsIEtaMin = cms.untracked.int32(1),
# ------- max EmFraction of probe jet allowed (does not affect isotracks)
maxProbeJetEmFrac = cms.untracked.double(0.1),
# ------- max EmFraction of tag jet allowed (does not affect isotracks)
maxTagJetEmFrac = cms.untracked.double(1.0),
# ------- max abs(eta) of tag jet allowed (does not affect isotracks)
maxTagJetAbsEta = cms.untracked.double(1.0),
# ------- min Et of tag jet allowed (does not affect isotracks)
minTagJetEt = cms.untracked.double(30.0),
# ------- min abs(eta) of probe jet allowed (does not affect isotracks)
minProbeJetAbsEta = cms.untracked.double(1.4),
# ------- file containing the phi symmetry corrections (to apply on the fly in isoTrack calibration)
# ------- USER PROVIDED!
phiSymCorFileName = cms.untracked.string("phiSymCor.txt"),
# ------- Flag to read phi symmetry corrections and apply on the fly
# ------- Set to false if already applied in Analyzer
applyPhiSymCorFlag = cms.untracked.bool(False),
# ------- output file with correction factors
outputCorCoefFileName = cms.untracked.string("calibConst_IsoTrk.txt"),
# ------- output file with some hisograms
histoFileName = cms.untracked.string("histoFile.root")
)
process.p = cms.Path(process.hcalCalib)
|