Fork me on GitHub

source: git/cards/delphes_card_IDEAtrkCov.tcl@ 8d7319c

Last change on this file since 8d7319c was 95e4e44, checked in by michele <michele.selvaggi@…>, 5 years ago

updated IDA card with geometry in delphes card

  • Property mode set to 100644
File size: 24.7 KB
Line 
1#################################################################### l
2# FCC-ee IDEA detector model
3#
4# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
5# email: efontane@bo.infn.it,
6# lorenzo.pezzotti01@universitadipavia.it,
7# m.antonello@uninsubria.it,
8# michele.selvaggi@cern.ch
9#####################################################################
10#
11#######################################
12# Order of execution of various modules
13#######################################
14
15set ExecutionPath {
16 ParticlePropagator
17
18 ChargedHadronTrackingEfficiency
19 ElectronTrackingEfficiency
20 MuonTrackingEfficiency
21
22 TrackMergerPre
23 TrackSmearing
24
25 TrackMerger
26 Calorimeter
27 EFlowMerger
28
29 PhotonEfficiency
30 PhotonIsolation
31
32 MuonFilter
33
34 ElectronFilter
35 ElectronEfficiency
36 ElectronIsolation
37
38 MuonEfficiency
39 MuonIsolation
40
41 MissingET
42
43 NeutrinoFilter
44 GenJetFinder
45 GenMissingET
46
47 FastJetFinder
48
49 JetEnergyScale
50
51 JetFlavorAssociation
52
53 BTagging
54 TauTagging
55
56 UniqueObjectFinder
57
58 ScalarHT
59 TreeWriter
60}
61
62
63#################################
64# Propagate particles in cylinder
65#################################
66
67module ParticlePropagator ParticlePropagator {
68 set InputArray Delphes/stableParticles
69
70 set OutputArray stableParticles
71 set ChargedHadronOutputArray chargedHadrons
72 set ElectronOutputArray electrons
73 set MuonOutputArray muons
74
75 # inner radius of the solenoid, in m
76 set Radius 2.25
77
78 # half-length: z of the solenoid, in m
79 set HalfLength 2.5
80
81 # magnetic field, in T
82 set Bz 2.0
83}
84
85####################################
86# Charged hadron tracking efficiency
87####################################
88
89module Efficiency ChargedHadronTrackingEfficiency {
90 set InputArray ParticlePropagator/chargedHadrons
91 set OutputArray chargedHadrons
92 # We use only one efficiency, we set only 0 effincency out of eta bounds:
93
94 set EfficiencyFormula {
95 (abs(eta) > 3.0) * (0.000) +
96 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
97 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
98 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
99 }
100}
101
102# (pt <= 0.1) * (0.00) +
103# (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
104# (abs(eta) > 3) * (0.00)
105
106
107
108##############################
109# Electron tracking efficiency
110##############################
111
112module Efficiency ElectronTrackingEfficiency {
113 set InputArray ParticlePropagator/electrons
114 set OutputArray electrons
115
116
117 # Current full simulation with CLICdet provides for electrons:
118 set EfficiencyFormula {
119 (abs(eta) > 3.0) * (0.000) +
120 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
121 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
122 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
123 }
124}
125
126
127##########################
128# Muon tracking efficiency
129##########################
130
131module Efficiency MuonTrackingEfficiency {
132 set InputArray ParticlePropagator/muons
133 set OutputArray muons
134
135 # Current full simulation with CLICdet provides for muons:
136 set EfficiencyFormula {
137 (abs(eta) > 3.0) * (0.000) +
138 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
139 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
140 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
141 }
142}
143
144##############
145# Track merger
146##############
147
148module Merger TrackMergerPre {
149# add InputArray InputArray
150 add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
151 add InputArray ElectronTrackingEfficiency/electrons
152 add InputArray MuonTrackingEfficiency/muons
153 set OutputArray tracks
154}
155
156
157########################################
158# Smearing for charged tracks
159########################################
160
161module TrackCovariance TrackSmearing {
162 set InputArray TrackMergerPre/tracks
163 set OutputArray tracks
164
165 ## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt
166 set DetectorGeometry {
167
168 1 PIPE -100 100 0.015 0.0012 0.35276 0 0 0 0 0 0
169 1 VTXLOW -0.12 0.12 0.017 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
170 1 VTXLOW -0.16 0.16 0.023 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
171 1 VTXLOW -0.16 0.16 0.031 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
172 1 VTXHIGH -1 1 0.32 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
173 1 VTXHIGH -1.05 1.05 0.34 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
174 1 DCHCANI -2.125 2.125 0.345 0.0002 0.237223 0 0 0 0 0 0
175 1 DCH -2 2 0.36 0.0147748 1400 1 0.0203738 0 0.0001 0 1
176 1 DCH -2 2 0.374775 0.0147748 1400 1 -0.0212097 0 0.0001 0 1
177 1 DCH -2 2 0.38955 0.0147748 1400 1 0.0220456 0 0.0001 0 1
178 1 DCH -2 2 0.404324 0.0147748 1400 1 -0.0228814 0 0.0001 0 1
179 1 DCH -2 2 0.419099 0.0147748 1400 1 0.0237172 0 0.0001 0 1
180 1 DCH -2 2 0.433874 0.0147748 1400 1 -0.024553 0 0.0001 0 1
181 1 DCH -2 2 0.448649 0.0147748 1400 1 0.0253888 0 0.0001 0 1
182 1 DCH -2 2 0.463423 0.0147748 1400 1 -0.0262245 0 0.0001 0 1
183 1 DCH -2 2 0.478198 0.0147748 1400 1 0.0270602 0 0.0001 0 1
184 1 DCH -2 2 0.492973 0.0147748 1400 1 -0.0278958 0 0.0001 0 1
185 1 DCH -2 2 0.507748 0.0147748 1400 1 0.0287314 0 0.0001 0 1
186 1 DCH -2 2 0.522523 0.0147748 1400 1 -0.029567 0 0.0001 0 1
187 1 DCH -2 2 0.537297 0.0147748 1400 1 0.0304025 0 0.0001 0 1
188 1 DCH -2 2 0.552072 0.0147748 1400 1 -0.031238 0 0.0001 0 1
189 1 DCH -2 2 0.566847 0.0147748 1400 1 0.0320734 0 0.0001 0 1
190 1 DCH -2 2 0.581622 0.0147748 1400 1 -0.0329088 0 0.0001 0 1
191 1 DCH -2 2 0.596396 0.0147748 1400 1 0.0337442 0 0.0001 0 1
192 1 DCH -2 2 0.611171 0.0147748 1400 1 -0.0345795 0 0.0001 0 1
193 1 DCH -2 2 0.625946 0.0147748 1400 1 0.0354147 0 0.0001 0 1
194 1 DCH -2 2 0.640721 0.0147748 1400 1 -0.0362499 0 0.0001 0 1
195 1 DCH -2 2 0.655495 0.0147748 1400 1 0.0370851 0 0.0001 0 1
196 1 DCH -2 2 0.67027 0.0147748 1400 1 -0.0379202 0 0.0001 0 1
197 1 DCH -2 2 0.685045 0.0147748 1400 1 0.0387552 0 0.0001 0 1
198 1 DCH -2 2 0.69982 0.0147748 1400 1 -0.0395902 0 0.0001 0 1
199 1 DCH -2 2 0.714595 0.0147748 1400 1 0.0404252 0 0.0001 0 1
200 1 DCH -2 2 0.729369 0.0147748 1400 1 -0.04126 0 0.0001 0 1
201 1 DCH -2 2 0.744144 0.0147748 1400 1 0.0420949 0 0.0001 0 1
202 1 DCH -2 2 0.758919 0.0147748 1400 1 -0.0429296 0 0.0001 0 1
203 1 DCH -2 2 0.773694 0.0147748 1400 1 0.0437643 0 0.0001 0 1
204 1 DCH -2 2 0.788468 0.0147748 1400 1 -0.044599 0 0.0001 0 1
205 1 DCH -2 2 0.803243 0.0147748 1400 1 0.0454336 0 0.0001 0 1
206 1 DCH -2 2 0.818018 0.0147748 1400 1 -0.0462681 0 0.0001 0 1
207 1 DCH -2 2 0.832793 0.0147748 1400 1 0.0471025 0 0.0001 0 1
208 1 DCH -2 2 0.847568 0.0147748 1400 1 -0.0479369 0 0.0001 0 1
209 1 DCH -2 2 0.862342 0.0147748 1400 1 0.0487713 0 0.0001 0 1
210 1 DCH -2 2 0.877117 0.0147748 1400 1 -0.0496055 0 0.0001 0 1
211 1 DCH -2 2 0.891892 0.0147748 1400 1 0.0504397 0 0.0001 0 1
212 1 DCH -2 2 0.906667 0.0147748 1400 1 -0.0512738 0 0.0001 0 1
213 1 DCH -2 2 0.921441 0.0147748 1400 1 0.0521079 0 0.0001 0 1
214 1 DCH -2 2 0.936216 0.0147748 1400 1 -0.0529418 0 0.0001 0 1
215 1 DCH -2 2 0.950991 0.0147748 1400 1 0.0537757 0 0.0001 0 1
216 1 DCH -2 2 0.965766 0.0147748 1400 1 -0.0546095 0 0.0001 0 1
217 1 DCH -2 2 0.980541 0.0147748 1400 1 0.0554433 0 0.0001 0 1
218 1 DCH -2 2 0.995315 0.0147748 1400 1 -0.056277 0 0.0001 0 1
219 1 DCH -2 2 1.01009 0.0147748 1400 1 0.0571106 0 0.0001 0 1
220 1 DCH -2 2 1.02486 0.0147748 1400 1 -0.0579441 0 0.0001 0 1
221 1 DCH -2 2 1.03964 0.0147748 1400 1 0.0587775 0 0.0001 0 1
222 1 DCH -2 2 1.05441 0.0147748 1400 1 -0.0596108 0 0.0001 0 1
223 1 DCH -2 2 1.06919 0.0147748 1400 1 0.0604441 0 0.0001 0 1
224 1 DCH -2 2 1.08396 0.0147748 1400 1 -0.0612773 0 0.0001 0 1
225 1 DCH -2 2 1.09874 0.0147748 1400 1 0.0621104 0 0.0001 0 1
226 1 DCH -2 2 1.11351 0.0147748 1400 1 -0.0629434 0 0.0001 0 1
227 1 DCH -2 2 1.12829 0.0147748 1400 1 0.0637763 0 0.0001 0 1
228 1 DCH -2 2 1.14306 0.0147748 1400 1 -0.0646092 0 0.0001 0 1
229 1 DCH -2 2 1.15784 0.0147748 1400 1 0.0654419 0 0.0001 0 1
230 1 DCH -2 2 1.17261 0.0147748 1400 1 -0.0662746 0 0.0001 0 1
231 1 DCH -2 2 1.18739 0.0147748 1400 1 0.0671071 0 0.0001 0 1
232 1 DCH -2 2 1.20216 0.0147748 1400 1 -0.0679396 0 0.0001 0 1
233 1 DCH -2 2 1.21694 0.0147748 1400 1 0.068772 0 0.0001 0 1
234 1 DCH -2 2 1.23171 0.0147748 1400 1 -0.0696042 0 0.0001 0 1
235 1 DCH -2 2 1.24649 0.0147748 1400 1 0.0704364 0 0.0001 0 1
236 1 DCH -2 2 1.26126 0.0147748 1400 1 -0.0712685 0 0.0001 0 1
237 1 DCH -2 2 1.27604 0.0147748 1400 1 0.0721005 0 0.0001 0 1
238 1 DCH -2 2 1.29081 0.0147748 1400 1 -0.0729324 0 0.0001 0 1
239 1 DCH -2 2 1.30559 0.0147748 1400 1 0.0737642 0 0.0001 0 1
240 1 DCH -2 2 1.32036 0.0147748 1400 1 -0.0745958 0 0.0001 0 1
241 1 DCH -2 2 1.33514 0.0147748 1400 1 0.0754274 0 0.0001 0 1
242 1 DCH -2 2 1.34991 0.0147748 1400 1 -0.0762589 0 0.0001 0 1
243 1 DCH -2 2 1.36468 0.0147748 1400 1 0.0770903 0 0.0001 0 1
244 1 DCH -2 2 1.37946 0.0147748 1400 1 -0.0779215 0 0.0001 0 1
245 1 DCH -2 2 1.39423 0.0147748 1400 1 0.0787527 0 0.0001 0 1
246 1 DCH -2 2 1.40901 0.0147748 1400 1 -0.0795837 0 0.0001 0 1
247 1 DCH -2 2 1.42378 0.0147748 1400 1 0.0804147 0 0.0001 0 1
248 1 DCH -2 2 1.43856 0.0147748 1400 1 -0.0812455 0 0.0001 0 1
249 1 DCH -2 2 1.45333 0.0147748 1400 1 0.0820762 0 0.0001 0 1
250 1 DCH -2 2 1.46811 0.0147748 1400 1 -0.0829068 0 0.0001 0 1
251 1 DCH -2 2 1.48288 0.0147748 1400 1 0.0837373 0 0.0001 0 1
252 1 DCH -2 2 1.49766 0.0147748 1400 1 -0.0845677 0 0.0001 0 1
253 1 DCH -2 2 1.51243 0.0147748 1400 1 0.0853979 0 0.0001 0 1
254 1 DCH -2 2 1.52721 0.0147748 1400 1 -0.086228 0 0.0001 0 1
255 1 DCH -2 2 1.54198 0.0147748 1400 1 0.087058 0 0.0001 0 1
256 1 DCH -2 2 1.55676 0.0147748 1400 1 -0.0878879 0 0.0001 0 1
257 1 DCH -2 2 1.57153 0.0147748 1400 1 0.0887177 0 0.0001 0 1
258 1 DCH -2 2 1.58631 0.0147748 1400 1 -0.0895474 0 0.0001 0 1
259 1 DCH -2 2 1.60108 0.0147748 1400 1 0.0903769 0 0.0001 0 1
260 1 DCH -2 2 1.61586 0.0147748 1400 1 -0.0912063 0 0.0001 0 1
261 1 DCH -2 2 1.63063 0.0147748 1400 1 0.0920356 0 0.0001 0 1
262 1 DCH -2 2 1.64541 0.0147748 1400 1 -0.0928647 0 0.0001 0 1
263 1 DCH -2 2 1.66018 0.0147748 1400 1 0.0936937 0 0.0001 0 1
264 1 DCH -2 2 1.67495 0.0147748 1400 1 -0.0945226 0 0.0001 0 1
265 1 DCH -2 2 1.68973 0.0147748 1400 1 0.0953514 0 0.0001 0 1
266 1 DCH -2 2 1.7045 0.0147748 1400 1 -0.09618 0 0.0001 0 1
267 1 DCH -2 2 1.71928 0.0147748 1400 1 0.0970085 0 0.0001 0 1
268 1 DCH -2 2 1.73405 0.0147748 1400 1 -0.0978369 0 0.0001 0 1
269 1 DCH -2 2 1.74883 0.0147748 1400 1 0.0986651 0 0.0001 0 1
270 1 DCH -2 2 1.7636 0.0147748 1400 1 -0.0994932 0 0.0001 0 1
271 1 DCH -2 2 1.77838 0.0147748 1400 1 0.100321 0 0.0001 0 1
272 1 DCH -2 2 1.79315 0.0147748 1400 1 -0.101149 0 0.0001 0 1
273 1 DCH -2 2 1.80793 0.0147748 1400 1 0.101977 0 0.0001 0 1
274 1 DCH -2 2 1.8227 0.0147748 1400 1 -0.102804 0 0.0001 0 1
275 1 DCH -2 2 1.83748 0.0147748 1400 1 0.103632 0 0.0001 0 1
276 1 DCH -2 2 1.85225 0.0147748 1400 1 -0.104459 0 0.0001 0 1
277 1 DCH -2 2 1.86703 0.0147748 1400 1 0.105286 0 0.0001 0 1
278 1 DCH -2 2 1.8818 0.0147748 1400 1 -0.106113 0 0.0001 0 1
279 1 DCH -2 2 1.89658 0.0147748 1400 1 0.10694 0 0.0001 0 1
280 1 DCH -2 2 1.91135 0.0147748 1400 1 -0.107766 0 0.0001 0 1
281 1 DCH -2 2 1.92613 0.0147748 1400 1 0.108593 0 0.0001 0 1
282 1 DCH -2 2 1.9409 0.0147748 1400 1 -0.109419 0 0.0001 0 1
283 1 DCH -2 2 1.95568 0.0147748 1400 1 0.110246 0 0.0001 0 1
284 1 DCH -2 2 1.97045 0.0147748 1400 1 -0.111072 0 0.0001 0 1
285 1 DCH -2 2 1.98523 0.0147748 1400 1 0.111898 0 0.0001 0 1
286 1 DCH -2 2 2 0.0147748 1400 1 -0.112723 0 0.0001 0 1
287 1 DCHCANO -2.125 2.125 2.02 0.02 1.667 0 0 0 0 0 0
288 1 BSILWRP -2.35 2.35 2.04 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
289 1 BSILWRP -2.35 2.35 2.06 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
290 1 MAG -2.5 2.5 2.25 0.05 0.0658 0 0 0 0 0 0
291 1 BPRESH -2.55 2.55 2.45 0.02 1 2 0 1.5708 7e-005 0.01 1
292 2 VTXDSK 0.141 0.3 -0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
293 2 VTXDSK 0.138 0.3 -0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
294 2 VTXDSK 0.065 0.3 -0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
295 2 VTXDSK 0.062 0.3 -0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
296 2 VTXDSK 0.062 0.3 0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
297 2 VTXDSK 0.065 0.3 0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
298 2 VTXDSK 0.138 0.3 0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
299 2 VTXDSK 0.141 0.3 0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
300 2 DCHWALL 0.345 2.02 2.125 0.25 5.55 0 0 0 0 0 0
301 2 DCHWALL 0.345 2.02 -2.125 0.25 5.55 0 0 0 0 0 0
302 2 FSILWRP 0.354 2.02 -2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
303 2 FSILWRP 0.35 2.02 -2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
304 2 FSILWRP 0.35 2.02 2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
305 2 FSILWRP 0.354 2.02 2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
306 2 FRAD 0.38 2.09 2.49 0.0043 0.005612 0 0 0 0 0 0
307 2 FRAD 0.38 2.09 -2.49 0.0043 0.005612 0 0 0 0 0 0
308 2 FPRESH 0.39 2.43 -2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
309 2 FPRESH 0.39 2.43 2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
310 }
311
312 set Bz 2.0
313}
314
315
316##############
317# Track merger
318##############
319
320module Merger TrackMerger {
321# add InputArray InputArray
322 add InputArray TrackSmearing/tracks
323 set OutputArray tracks
324}
325
326
327#############
328# Calorimeter
329#############
330module DualReadoutCalorimeter Calorimeter {
331 set ParticleInputArray ParticlePropagator/stableParticles
332 set TrackInputArray TrackMerger/tracks
333
334 set TowerOutputArray towers
335 set PhotonOutputArray photons
336
337 set EFlowTrackOutputArray eflowTracks
338 set EFlowPhotonOutputArray eflowPhotons
339 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
340
341 set ECalEnergyMin 0.5
342 set HCalEnergyMin 0.5
343 set EnergyMin 0.5
344 set ECalEnergySignificanceMin 1.0
345 set HCalEnergySignificanceMin 1.0
346 set EnergySignificanceMin 1.0
347
348 set SmearTowerCenter true
349 set pi [expr {acos(-1)}]
350
351 # Lists of the edges of each tower in eta and phi;
352 # each list starts with the lower edge of the first tower;
353 # the list ends with the higher edged of the last tower.
354 # Barrel: deta=0.02 towers up to |eta| <= 0.88 ( up to 45°)
355 # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad)
356 # Cell size: about 6 cm x 6 cm
357
358 #barrel:
359 set PhiBins {}
360 for {set i -120} {$i <= 120} {incr i} {
361 add PhiBins [expr {$i * $pi/120}]
362 }
363 #deta=0.02 units for |eta| <= 0.88
364 for {set i -44} {$i < 45} {incr i} {
365 set eta [expr {$i * 0.02}]
366 add EtaPhiBins $eta $PhiBins
367 }
368
369 #endcaps:
370 set PhiBins {}
371 for {set i -120} {$i <= 120} {incr i} {
372 add PhiBins [expr {$i* $pi/120}]
373 }
374 #deta=0.02 units for 0.88 < |eta| <= 3.0
375 #first, from -3.0 to -0.88
376 for {set i 1} {$i <=106} {incr i} {
377 set eta [expr {-3.00 + $i * 0.02}]
378 add EtaPhiBins $eta $PhiBins
379 }
380 #same for 0.88 to 3.0
381 for {set i 1} {$i <=106} {incr i} {
382 set eta [expr {0.88 + $i * 0.02}]
383 add EtaPhiBins $eta $PhiBins
384 }
385
386 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
387 add EnergyFraction {0} {0.0 1.0}
388 # energy fractions for e, gamma and pi0
389 add EnergyFraction {11} {1.0 0.0}
390 add EnergyFraction {22} {1.0 0.0}
391 add EnergyFraction {111} {1.0 0.0}
392 # energy fractions for muon, neutrinos and neutralinos
393 add EnergyFraction {12} {0.0 0.0}
394 add EnergyFraction {13} {0.0 0.0}
395 add EnergyFraction {14} {0.0 0.0}
396 add EnergyFraction {16} {0.0 0.0}
397 add EnergyFraction {1000022} {0.0 0.0}
398 add EnergyFraction {1000023} {0.0 0.0}
399 add EnergyFraction {1000025} {0.0 0.0}
400 add EnergyFraction {1000035} {0.0 0.0}
401 add EnergyFraction {1000045} {0.0 0.0}
402 # energy fractions for K0short and Lambda
403 add EnergyFraction {310} {0.3 0.7}
404 add EnergyFraction {3122} {0.3 0.7}
405
406
407 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
408 set ECalResolutionFormula {
409 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.11^2)+
410 (abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.11^2)
411 }
412
413 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
414 set HCalResolutionFormula {
415 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.30^2)+
416 (abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.30^2)
417 }
418}
419
420####################
421# Energy flow merger
422####################
423
424module Merger EFlowMerger {
425# add InputArray InputArray
426 add InputArray Calorimeter/eflowTracks
427 add InputArray Calorimeter/eflowPhotons
428 add InputArray Calorimeter/eflowNeutralHadrons
429 set OutputArray eflow
430}
431
432###################
433# Photon efficiency
434###################
435
436module Efficiency PhotonEfficiency {
437 set InputArray Calorimeter/eflowPhotons
438 set OutputArray photons
439
440 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
441 # efficiency formula for photons
442 set EfficiencyFormula {
443 (energy < 2.0) * (0.000)+
444 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
445 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
446 (abs(eta) > 3.0) * (0.000)
447 }
448}
449
450##################
451# Photon isolation
452##################
453
454module Isolation PhotonIsolation {
455 set CandidateInputArray PhotonEfficiency/photons
456 set IsolationInputArray EFlowMerger/eflow
457
458 set OutputArray photons
459
460 set DeltaRMax 0.5
461
462 set PTMin 0.5
463
464 set PTRatioMax 999.
465}
466
467#################
468# Electron filter
469#################
470
471module PdgCodeFilter ElectronFilter {
472 set InputArray Calorimeter/eflowTracks
473 set OutputArray electrons
474 set Invert true
475 add PdgCode {11}
476 add PdgCode {-11}
477}
478
479#################
480# Muon filter
481#################
482
483module PdgCodeFilter MuonFilter {
484 set InputArray Calorimeter/eflowTracks
485 set OutputArray muons
486 set Invert true
487 add PdgCode {13}
488 add PdgCode {-13}
489}
490
491
492#####################
493# Electron efficiency
494#####################
495
496module Efficiency ElectronEfficiency {
497 set InputArray ElectronFilter/electrons
498 set OutputArray electrons
499
500 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
501
502 # efficiency formula for electrons
503 set EfficiencyFormula {
504 (energy < 2.0) * (0.000)+
505 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
506 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
507 (abs(eta) > 3.0) * (0.000)
508 }
509}
510
511####################
512# Electron isolation
513####################
514
515module Isolation ElectronIsolation {
516 set CandidateInputArray ElectronEfficiency/electrons
517 set IsolationInputArray EFlowMerger/eflow
518
519 set OutputArray electrons
520
521 set DeltaRMax 0.5
522
523 set PTMin 0.5
524
525 set PTRatioMax 0.12
526}
527
528#################
529# Muon efficiency
530#################
531
532module Efficiency MuonEfficiency {
533 set InputArray MuonFilter/muons
534 set OutputArray muons
535
536 # set EfficiencyFormula {efficiency as a function of eta and pt}
537
538 # efficiency formula for muons
539 set EfficiencyFormula {
540 (energy < 2.0) * (0.000)+
541 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
542 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
543 (abs(eta) > 3.0) * (0.000)
544 }
545}
546
547################
548# Muon isolation
549################
550
551module Isolation MuonIsolation {
552 set CandidateInputArray MuonEfficiency/muons
553 set IsolationInputArray EFlowMerger/eflow
554
555 set OutputArray muons
556
557 set DeltaRMax 0.5
558
559 set PTMin 0.5
560
561 set PTRatioMax 0.25
562}
563
564###################
565# Missing ET merger
566###################
567
568module Merger MissingET {
569# add InputArray InputArray
570 add InputArray EFlowMerger/eflow
571 set MomentumOutputArray momentum
572}
573
574##################
575# Scalar HT merger
576##################
577
578module Merger ScalarHT {
579# add InputArray InputArray
580 add InputArray UniqueObjectFinder/jets
581 add InputArray UniqueObjectFinder/electrons
582 add InputArray UniqueObjectFinder/photons
583 add InputArray UniqueObjectFinder/muons
584 set EnergyOutputArray energy
585}
586
587#####################
588# Neutrino Filter
589#####################
590
591module PdgCodeFilter NeutrinoFilter {
592
593 set InputArray Delphes/stableParticles
594 set OutputArray filteredParticles
595
596 set PTMin 0.0
597
598 add PdgCode {12}
599 add PdgCode {14}
600 add PdgCode {16}
601 add PdgCode {-12}
602 add PdgCode {-14}
603 add PdgCode {-16}
604}
605
606
607#####################
608# MC truth jet finder
609#####################
610
611module FastJetFinder GenJetFinder {
612 set InputArray NeutrinoFilter/filteredParticles
613
614 set OutputArray jets
615
616 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
617 set JetAlgorithm 6
618 set ParameterR 0.4
619 set JetPTMin 1.0
620}
621
622
623#########################
624# Gen Missing ET merger
625########################
626
627module Merger GenMissingET {
628# add InputArray InputArray
629 add InputArray NeutrinoFilter/filteredParticles
630 set MomentumOutputArray momentum
631}
632
633############
634# Jet finder
635############
636
637module FastJetFinder FastJetFinder {
638# set InputArray Calorimeter/towers
639 set InputArray EFlowMerger/eflow
640
641 set OutputArray jets
642
643 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
644 set JetAlgorithm 6
645 set ParameterR 0.4
646 set JetPTMin 1.0
647}
648
649##################
650# Jet Energy Scale
651##################
652
653module EnergyScale JetEnergyScale {
654 set InputArray FastJetFinder/jets
655 set OutputArray jets
656
657 # scale formula for jets
658 set ScaleFormula {1.08}
659}
660
661########################
662# Jet Flavor Association
663########################
664
665module JetFlavorAssociation JetFlavorAssociation {
666
667 set PartonInputArray Delphes/partons
668 set ParticleInputArray Delphes/allParticles
669 set ParticleLHEFInputArray Delphes/allParticlesLHEF
670 set JetInputArray JetEnergyScale/jets
671
672 set DeltaR 0.5
673 set PartonPTMin 1.0
674 set PartonEtaMax 3.0
675}
676
677###########
678# b-tagging
679###########
680
681module BTagging BTagging {
682 set JetInputArray JetEnergyScale/jets
683
684 set BitNumber 0
685
686 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
687
688 # default efficiency formula (misidentification rate)
689 add EfficiencyFormula {0} {0.01}
690
691 # efficiency formula for c-jets (misidentification rate)
692 add EfficiencyFormula {4} {0.10}
693
694 # efficiency formula for b-jets
695 add EfficiencyFormula {5} {0.80}
696}
697
698#############
699# tau-tagging
700#############
701
702module TauTagging TauTagging {
703 set ParticleInputArray Delphes/allParticles
704 set PartonInputArray Delphes/partons
705 set JetInputArray JetEnergyScale/jets
706
707 set DeltaR 0.5
708 set TauPTMin 1.0
709 set TauEtaMax 3.0
710
711 # default efficiency formula (misidentification rate)
712 add EfficiencyFormula {0} {0.001}
713 # efficiency formula for tau-jets
714 add EfficiencyFormula {15} {0.6}
715}
716
717
718#####################################################
719# Find uniquely identified photons/electrons/tau/jets
720#####################################################
721
722module UniqueObjectFinder UniqueObjectFinder {
723# earlier arrays take precedence over later ones
724# add InputArray InputArray OutputArray
725 add InputArray PhotonIsolation/photons photons
726 add InputArray ElectronIsolation/electrons electrons
727 add InputArray MuonIsolation/muons muons
728 add InputArray JetEnergyScale/jets jets
729}
730
731
732
733##################
734# ROOT tree writer
735##################
736
737# Tracks, towers and eflow objects are not stored by default in the output.
738# If needed (for jet constituent or other studies), uncomment the relevant
739# "add Branch ..." lines.
740
741module TreeWriter TreeWriter {
742 # add Branch InputArray BranchName BranchClass
743
744 add Branch Delphes/allParticles Particle GenParticle
745
746 add Branch TrackMerger/tracks Track Track
747 add Branch Calorimeter/towers Tower Tower
748
749 add Branch Calorimeter/eflowTracks EFlowTrack Track
750 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
751 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
752
753 add Branch Calorimeter/photons CaloPhoton Photon
754 add Branch PhotonEfficiency/photons PhotonEff Photon
755 add Branch PhotonIsolation/photons PhotonIso Photon
756
757 add Branch GenJetFinder/jets GenJet Jet
758 add Branch GenMissingET/momentum GenMissingET MissingET
759
760 add Branch UniqueObjectFinder/jets Jet Jet
761 add Branch UniqueObjectFinder/electrons Electron Electron
762 add Branch UniqueObjectFinder/photons Photon Photon
763 add Branch UniqueObjectFinder/muons Muon Muon
764
765 add Branch JetEnergyScale/jets AntiKtJet Jet
766
767 add Branch MissingET/momentum MissingET MissingET
768 add Branch ScalarHT/energy ScalarHT ScalarHT
769}
Note: See TracBrowser for help on using the repository browser.