Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage3.tcl@ 27197df

Last change on this file since 27197df was 8d7319c, checked in by Aleksander Filip Zarnecki <zarnecki@…>, 4 years ago

Minor correction to description of tracking efficiency

  • Property mode set to 100644
File size: 41.4 KB
RevLine 
[7a1e27d]1#######################################
2# CLICdet model
[505a779]3# based on arXiv:1812.07337 and
4# CLICdp-Note-2017-001
5#
[7a1e27d]6# Ulrike Schnoor ulrike.schnoor@cern.ch
7#
[a15ac11]8# For the high energy stage of
9# CLIC: 3 TeV
10#
[7a1e27d]11# Jet finding with Valencia algorithm:
12# use exclusive clustering with njets
13# according to final state
[505a779]14#
15# c-tagging capabilities of CLICdet are
16# not yet implemented here. Please
17# contact us if you want to use it.
18#######################################
19
[7a1e27d]20#######################################
21# Order of execution of various modules
22#######################################
23
24set ExecutionPath {
25 ParticlePropagator
26
27 ChargedHadronTrackingEfficiency
28 ElectronTrackingEfficiency
29 MuonTrackingEfficiency
30
31 ChargedHadronMomentumSmearing
32 ElectronMomentumSmearing
33 MuonMomentumSmearing
34
35 TrackMerger
36
37 ECal
38 HCal
39
40 Calorimeter
41 EFlowMerger
42
43
44 PhotonEfficiency
45 PhotonIsolation
46
47 ElectronFilter
48 ElectronEfficiency
49 ElectronIsolation
50
51 ChargedHadronFilter
52
53 MuonEfficiency
54 MuonIsolation
55
56 EFlowFilter
57
58 NeutrinoFilter
59 GenJetFinder
60 FastJetFinderKt
61 FastJetFinderVLC_R05_N2
62 FastJetFinderVLC_R05_N3
63 FastJetFinderVLC_R05_N4
64 FastJetFinderVLC_R05_N5
65 FastJetFinderVLC_R05_N6
66 FastJetFinderVLC_R07_N2
67 FastJetFinderVLC_R07_N3
68 FastJetFinderVLC_R07_N4
69 FastJetFinderVLC_R07_N5
70 FastJetFinderVLC_R07_N6
71 FastJetFinderVLC_R10_N2
72 FastJetFinderVLC_R10_N3
73 FastJetFinderVLC_R10_N4
74 FastJetFinderVLC_R10_N5
75 FastJetFinderVLC_R10_N6
76 FastJetFinderVLC_R12_N2
77 FastJetFinderVLC_R12_N3
78 FastJetFinderVLC_R12_N4
79 FastJetFinderVLC_R12_N5
80 FastJetFinderVLC_R12_N6
81 FastJetFinderVLC_R15_N2
82 FastJetFinderVLC_R15_N3
83 FastJetFinderVLC_R15_N4
84 FastJetFinderVLC_R15_N5
85 FastJetFinderVLC_R15_N6
86
87 FastJetFinderVLC_R05_inclusive
88 FastJetFinderVLC_R07_inclusive
89 FastJetFinderVLC_R10_inclusive
90 FastJetFinderVLC_R12_inclusive
91 FastJetFinderVLC_R15_inclusive
92
93
94 MissingET
95 GenMissingET
96
97
98 JetMomentumSmearing_VLCR05N2
99 JetMomentumSmearing_VLCR05N3
100 JetMomentumSmearing_VLCR05N4
101 JetMomentumSmearing_VLCR05N5
102 JetMomentumSmearing_VLCR05N6
103 JetMomentumSmearing_VLCR05_inclusive
104 JetMomentumSmearing_VLCR07N2
105 JetMomentumSmearing_VLCR07N3
106 JetMomentumSmearing_VLCR07N4
107 JetMomentumSmearing_VLCR07N5
108 JetMomentumSmearing_VLCR07N6
109 JetMomentumSmearing_VLCR07_inclusive
110 JetMomentumSmearing_VLCR10N2
111 JetMomentumSmearing_VLCR10N3
112 JetMomentumSmearing_VLCR10N4
113 JetMomentumSmearing_VLCR10N5
114 JetMomentumSmearing_VLCR10N6
115 JetMomentumSmearing_VLCR10_inclusive
116 JetMomentumSmearing_VLCR12N2
117 JetMomentumSmearing_VLCR12N3
118 JetMomentumSmearing_VLCR12N4
119 JetMomentumSmearing_VLCR12N5
120 JetMomentumSmearing_VLCR12N6
121 JetMomentumSmearing_VLCR12_inclusive
122 JetMomentumSmearing_VLCR15N2
123 JetMomentumSmearing_VLCR15N3
124 JetMomentumSmearing_VLCR15N4
125 JetMomentumSmearing_VLCR15N5
126 JetMomentumSmearing_VLCR15N6
127 JetMomentumSmearing_VLCR15_inclusive
128
129
130
131
132 JetFlavorAssociation_R05N2
[0e6b14e]133 JetFlavorAssociation_R05N3
134 JetFlavorAssociation_R05N4
135 JetFlavorAssociation_R05N5
136 JetFlavorAssociation_R05N6
137
138 JetFlavorAssociation_R07N2
139 JetFlavorAssociation_R07N3
140 JetFlavorAssociation_R07N4
141 JetFlavorAssociation_R07N5
142 JetFlavorAssociation_R07N6
143
144 JetFlavorAssociation_R10N2
145 JetFlavorAssociation_R10N3
146 JetFlavorAssociation_R10N4
147 JetFlavorAssociation_R10N5
148 JetFlavorAssociation_R10N6
149
150 JetFlavorAssociation_R12N2
151 JetFlavorAssociation_R12N3
152 JetFlavorAssociation_R12N4
153 JetFlavorAssociation_R12N5
154 JetFlavorAssociation_R12N6
155
156 JetFlavorAssociation_R15N2
157 JetFlavorAssociation_R15N3
158 JetFlavorAssociation_R15N4
159 JetFlavorAssociation_R15N5
160 JetFlavorAssociation_R15N6
161
162 JetFlavorAssociation_R05_inclusive
163 JetFlavorAssociation_R07_inclusive
164 JetFlavorAssociation_R10_inclusive
165 JetFlavorAssociation_R12_inclusive
166 JetFlavorAssociation_R15_inclusive
167
168 JetFlavorAssociation_JER_R05N2
169 JetFlavorAssociation_JER_R05N3
170 JetFlavorAssociation_JER_R05N4
171 JetFlavorAssociation_JER_R05N5
172 JetFlavorAssociation_JER_R05N6
173
174 JetFlavorAssociation_JER_R07N2
175 JetFlavorAssociation_JER_R07N3
176 JetFlavorAssociation_JER_R07N4
177 JetFlavorAssociation_JER_R07N5
178 JetFlavorAssociation_JER_R07N6
179
180 JetFlavorAssociation_JER_R10N2
181 JetFlavorAssociation_JER_R10N3
182 JetFlavorAssociation_JER_R10N4
183 JetFlavorAssociation_JER_R10N5
184 JetFlavorAssociation_JER_R10N6
[7a1e27d]185
[0e6b14e]186 JetFlavorAssociation_JER_R12N2
187 JetFlavorAssociation_JER_R12N3
188 JetFlavorAssociation_JER_R12N4
189 JetFlavorAssociation_JER_R12N5
190 JetFlavorAssociation_JER_R12N6
191
192 JetFlavorAssociation_JER_R15N2
193 JetFlavorAssociation_JER_R15N3
194 JetFlavorAssociation_JER_R15N4
195 JetFlavorAssociation_JER_R15N5
196 JetFlavorAssociation_JER_R15N6
197
198 JetFlavorAssociation_JER_R05_inclusive
199 JetFlavorAssociation_JER_R07_inclusive
200 JetFlavorAssociation_JER_R10_inclusive
201 JetFlavorAssociation_JER_R12_inclusive
202 JetFlavorAssociation_JER_R15_inclusive
203
204 BTaggingWP50_R05N2
205 BTaggingWP70_R05N2
206 BTaggingWP90_R05N2
207 BTaggingWP50_R05N3
208 BTaggingWP70_R05N3
209 BTaggingWP90_R05N3
210 BTaggingWP50_R05N4
211 BTaggingWP70_R05N4
212 BTaggingWP90_R05N4
213 BTaggingWP50_R05N5
214 BTaggingWP70_R05N5
215 BTaggingWP90_R05N5
216 BTaggingWP50_R05N6
217 BTaggingWP70_R05N6
218 BTaggingWP90_R05N6
219 BTaggingWP50_R07N2
220 BTaggingWP70_R07N2
221 BTaggingWP90_R07N2
222 BTaggingWP50_R07N3
223 BTaggingWP70_R07N3
224 BTaggingWP90_R07N3
225 BTaggingWP50_R07N4
226 BTaggingWP70_R07N4
227 BTaggingWP90_R07N4
228 BTaggingWP50_R07N5
229 BTaggingWP70_R07N5
230 BTaggingWP90_R07N5
231 BTaggingWP50_R07N6
232 BTaggingWP70_R07N6
233 BTaggingWP90_R07N6
234 BTaggingWP50_R10N2
235 BTaggingWP70_R10N2
236 BTaggingWP90_R10N2
237 BTaggingWP50_R10N3
238 BTaggingWP70_R10N3
239 BTaggingWP90_R10N3
240 BTaggingWP50_R10N4
241 BTaggingWP70_R10N4
242 BTaggingWP90_R10N4
243 BTaggingWP50_R10N5
244 BTaggingWP70_R10N5
245 BTaggingWP90_R10N5
246 BTaggingWP50_R10N6
247 BTaggingWP70_R10N6
248 BTaggingWP90_R10N6
249 BTaggingWP50_R12N2
250 BTaggingWP70_R12N2
251 BTaggingWP90_R12N2
252 BTaggingWP50_R12N3
253 BTaggingWP70_R12N3
254 BTaggingWP90_R12N3
255 BTaggingWP50_R12N4
256 BTaggingWP70_R12N4
257 BTaggingWP90_R12N4
258 BTaggingWP50_R12N5
259 BTaggingWP70_R12N5
260 BTaggingWP90_R12N5
261 BTaggingWP50_R12N6
262 BTaggingWP70_R12N6
263 BTaggingWP90_R12N6
264 BTaggingWP50_R15N2
265 BTaggingWP70_R15N2
266 BTaggingWP90_R15N2
267 BTaggingWP50_R15N3
268 BTaggingWP70_R15N3
269 BTaggingWP90_R15N3
270 BTaggingWP50_R15N4
271 BTaggingWP70_R15N4
272 BTaggingWP90_R15N4
273 BTaggingWP50_R15N5
274 BTaggingWP70_R15N5
275 BTaggingWP90_R15N5
276 BTaggingWP50_R15N6
277 BTaggingWP70_R15N6
278 BTaggingWP90_R15N6
[7a1e27d]279 BTaggingWP50_R05_inclusive
280 BTaggingWP70_R05_inclusive
281 BTaggingWP90_R05_inclusive
282 BTaggingWP50_R07_inclusive
283 BTaggingWP70_R07_inclusive
284 BTaggingWP90_R07_inclusive
285 BTaggingWP50_R10_inclusive
286 BTaggingWP70_R10_inclusive
287 BTaggingWP90_R10_inclusive
288 BTaggingWP50_R12_inclusive
289 BTaggingWP70_R12_inclusive
290 BTaggingWP90_R12_inclusive
291 BTaggingWP50_R15_inclusive
292 BTaggingWP70_R15_inclusive
293 BTaggingWP90_R15_inclusive
294
[0e6b14e]295 BTagging_JER_WP50_R05N2
296 BTagging_JER_WP70_R05N2
297 BTagging_JER_WP90_R05N2
298 BTagging_JER_WP50_R05N3
299 BTagging_JER_WP70_R05N3
300 BTagging_JER_WP90_R05N3
301 BTagging_JER_WP50_R05N4
302 BTagging_JER_WP70_R05N4
303 BTagging_JER_WP90_R05N4
304 BTagging_JER_WP50_R05N5
305 BTagging_JER_WP70_R05N5
306 BTagging_JER_WP90_R05N5
307 BTagging_JER_WP50_R05N6
308 BTagging_JER_WP70_R05N6
309 BTagging_JER_WP90_R05N6
310 BTagging_JER_WP50_R07N2
311 BTagging_JER_WP70_R07N2
312 BTagging_JER_WP90_R07N2
313 BTagging_JER_WP50_R07N3
314 BTagging_JER_WP70_R07N3
315 BTagging_JER_WP90_R07N3
316 BTagging_JER_WP50_R07N4
317 BTagging_JER_WP70_R07N4
318 BTagging_JER_WP90_R07N4
319 BTagging_JER_WP50_R07N5
320 BTagging_JER_WP70_R07N5
321 BTagging_JER_WP90_R07N5
322 BTagging_JER_WP50_R07N6
323 BTagging_JER_WP70_R07N6
324 BTagging_JER_WP90_R07N6
325 BTagging_JER_WP50_R10N2
326 BTagging_JER_WP70_R10N2
327 BTagging_JER_WP90_R10N2
328 BTagging_JER_WP50_R10N3
329 BTagging_JER_WP70_R10N3
330 BTagging_JER_WP90_R10N3
331 BTagging_JER_WP50_R10N4
332 BTagging_JER_WP70_R10N4
333 BTagging_JER_WP90_R10N4
334 BTagging_JER_WP50_R10N5
335 BTagging_JER_WP70_R10N5
336 BTagging_JER_WP90_R10N5
337 BTagging_JER_WP50_R10N6
338 BTagging_JER_WP70_R10N6
339 BTagging_JER_WP90_R10N6
340 BTagging_JER_WP50_R12N2
341 BTagging_JER_WP70_R12N2
342 BTagging_JER_WP90_R12N2
343 BTagging_JER_WP50_R12N3
344 BTagging_JER_WP70_R12N3
345 BTagging_JER_WP90_R12N3
346 BTagging_JER_WP50_R12N4
347 BTagging_JER_WP70_R12N4
348 BTagging_JER_WP90_R12N4
349 BTagging_JER_WP50_R12N5
350 BTagging_JER_WP70_R12N5
351 BTagging_JER_WP90_R12N5
352 BTagging_JER_WP50_R12N6
353 BTagging_JER_WP70_R12N6
354 BTagging_JER_WP90_R12N6
355 BTagging_JER_WP50_R15N2
356 BTagging_JER_WP70_R15N2
357 BTagging_JER_WP90_R15N2
358 BTagging_JER_WP50_R15N3
359 BTagging_JER_WP70_R15N3
360 BTagging_JER_WP90_R15N3
361 BTagging_JER_WP50_R15N4
362 BTagging_JER_WP70_R15N4
363 BTagging_JER_WP90_R15N4
364 BTagging_JER_WP50_R15N5
365 BTagging_JER_WP70_R15N5
366 BTagging_JER_WP90_R15N5
367 BTagging_JER_WP50_R15N6
368 BTagging_JER_WP70_R15N6
369 BTagging_JER_WP90_R15N6
370 BTagging_JER_WP50_R05_inclusive
371 BTagging_JER_WP70_R05_inclusive
372 BTagging_JER_WP90_R05_inclusive
373 BTagging_JER_WP50_R07_inclusive
374 BTagging_JER_WP70_R07_inclusive
375 BTagging_JER_WP90_R07_inclusive
376 BTagging_JER_WP50_R10_inclusive
377 BTagging_JER_WP70_R10_inclusive
378 BTagging_JER_WP90_R10_inclusive
379 BTagging_JER_WP50_R12_inclusive
380 BTagging_JER_WP70_R12_inclusive
381 BTagging_JER_WP90_R12_inclusive
382 BTagging_JER_WP50_R15_inclusive
383 BTagging_JER_WP70_R15_inclusive
384 BTagging_JER_WP90_R15_inclusive
385
[7a1e27d]386 TauTagging_R05N2
387 TauTagging_R05N3
388 TauTagging_R05N4
389 TauTagging_R05N5
390 TauTagging_R05N6
391 TauTagging_R07N2
392 TauTagging_R07N3
393 TauTagging_R07N4
394 TauTagging_R07N5
395 TauTagging_R07N6
396 TauTagging_R10N2
397 TauTagging_R10N3
398 TauTagging_R10N4
399 TauTagging_R10N5
400 TauTagging_R10N6
401 TauTagging_R12N2
402 TauTagging_R12N3
403 TauTagging_R12N4
404 TauTagging_R12N5
405 TauTagging_R12N6
406 TauTagging_R15N2
407 TauTagging_R15N3
408 TauTagging_R15N4
409 TauTagging_R15N5
410 TauTagging_R15N6
411 TauTagging_R05_inclusive
412 TauTagging_R07_inclusive
413 TauTagging_R10_inclusive
414 TauTagging_R12_inclusive
415 TauTagging_R15_inclusive
416
[0e6b14e]417 TauTagging_JER_R05N2
418 TauTagging_JER_R05N3
419 TauTagging_JER_R05N4
420 TauTagging_JER_R05N5
421 TauTagging_JER_R05N6
422 TauTagging_JER_R07N2
423 TauTagging_JER_R07N3
424 TauTagging_JER_R07N4
425 TauTagging_JER_R07N5
426 TauTagging_JER_R07N6
427 TauTagging_JER_R10N2
428 TauTagging_JER_R10N3
429 TauTagging_JER_R10N4
430 TauTagging_JER_R10N5
431 TauTagging_JER_R10N6
432 TauTagging_JER_R12N2
433 TauTagging_JER_R12N3
434 TauTagging_JER_R12N4
435 TauTagging_JER_R12N5
436 TauTagging_JER_R12N6
437 TauTagging_JER_R15N2
438 TauTagging_JER_R15N3
439 TauTagging_JER_R15N4
440 TauTagging_JER_R15N5
441 TauTagging_JER_R15N6
442 TauTagging_JER_R05_inclusive
443 TauTagging_JER_R07_inclusive
444 TauTagging_JER_R10_inclusive
445 TauTagging_JER_R12_inclusive
446 TauTagging_JER_R15_inclusive
[7a1e27d]447
448
449 ScalarHT
450
451
452 TreeWriter
453}
454
455#################################
456# Propagate particles in cylinder
457#################################
458
459module ParticlePropagator ParticlePropagator {
460 set InputArray Delphes/stableParticles
461
462 set OutputArray stableParticles
463 set ChargedHadronOutputArray chargedHadrons
464 set ElectronOutputArray electrons
465 set MuonOutputArray muons
466
467 # radius of the magnetic field coverage in the calorimeter, in m
468 set Radius 1.5
469 # half-length of the magnetic field coverage in the calorimeter, in m
470 set HalfLength 2.31
471
472 # magnetic field, in T
473 set Bz 4.0
474}
475
476####################################
477# Charged hadron tracking efficiency
478####################################
479
480module Efficiency ChargedHadronTrackingEfficiency {
481 set InputArray ParticlePropagator/chargedHadrons
482 set OutputArray chargedHadrons
483 # Current full simulation with CLICdet provides for pions:
484
485 set EfficiencyFormula {
[8d7319c]486 (pt <= 0.1) * (0.000) +
487 (pt > 0.1) * (abs(eta) > 2.54) * (0.000) +
488 (pt > 0.1) * (energy >= 80) * (abs(eta) < 2.54) * (1.000) +
489 (pt > 0.1) * (energy < 80 && energy >= 3) * (abs(eta) <=2.54 && abs(eta) > 2.34) * (0.994) +
490 (pt > 0.1) * (energy < 80 && energy >= 3) * (abs(eta) <= 2.34) * (1.000) +
491 (pt > 0.1) * (energy < 3) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.990) +
492 (pt > 0.1) * (energy < 3) * (abs(eta) <= 0.55 ) * (1.000)
[7a1e27d]493 }
494}
495
496##############################
497# Electron tracking efficiency
498##############################
499
500module Efficiency ElectronTrackingEfficiency {
501 set InputArray ParticlePropagator/electrons
502 set OutputArray electrons
503
504
505 # Current full simulation with CLICdet provides for electrons:
506 set EfficiencyFormula {
[8d7319c]507 (pt <= 0.1) * (0.000) +
508 (pt > 0.1) * (abs(eta) > 2.54) * (0.000) +
509 (pt > 0.1) * (energy >= 80) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) +
510 (pt > 0.1) * (energy >= 80) * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) +
511 (pt > 0.1) * (energy >= 80) * (abs(eta) <= 2.34 ) * (1.000) +
512 (pt > 0.1) * (energy < 80 && energy >= 5) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) +
513 (pt > 0.1) * (energy < 80 && energy >= 5) * (abs(eta) <= 2.17) * (1.000) +
514 (pt > 0.1) * (energy < 5) * (abs(eta) <= 2.54 && abs(eta) > 2.34 ) * (1.000) +
515 (pt > 0.1) * (energy < 5) * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) +
516 (pt > 0.1) * (energy < 5) * (abs(eta) <= 0.76) * (0.999)
[7a1e27d]517 }
518}
519
520##########################
521# Muon tracking efficiency
522##########################
523
524module Efficiency MuonTrackingEfficiency {
525 set InputArray ParticlePropagator/muons
526 set OutputArray muons
527
528 # Current full simulation with CLICdet provides for muons:
529 set EfficiencyFormula {
[8d7319c]530 (pt <= 0.1) * (0.000) +
531 (pt > 0.1) * (abs(eta) > 2.54) * (0.000) +
532 (pt > 0.1) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 80) * (0.994) +
533 (pt > 0.1) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 5 && energy < 80) * (0.996) +
534 (pt > 0.1) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy < 5 ) * (0.996) +
535 (pt > 0.1) * (abs(eta) <= 2.44 ) * (energy >= 5 ) * (1.000) +
536 (pt > 0.1) * (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy < 5 ) * (0.999) +
537 (pt > 0.1) * (abs(eta) <= 2.25 ) * (energy < 5 ) * (1.000)
[7a1e27d]538 }
539}
540
541########################################
542# Momentum resolution for charged tracks
543########################################
544
545module MomentumSmearing ChargedHadronMomentumSmearing {
546 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
547 set OutputArray chargedHadrons
548
549
550 # Resolution given in dpT/pT.
551 # CLICdet internal studies
552 set ResolutionFormula {
553 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
[6f96f62]554 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
555 (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2 ) +
[7a1e27d]556 (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2 ) +
557 (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2 ) +
558 (abs(eta) < 0.18) * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2 )
559
560 }
561}
562
563###################################
564# Momentum resolution for electrons
565###################################
566
567module MomentumSmearing ElectronMomentumSmearing {
568 set InputArray ElectronTrackingEfficiency/electrons
569 set OutputArray electrons
570
571 # Resolution given in dpT/pT.
572 # CLICdet internal studies
573 set ResolutionFormula {
574 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
[6f96f62]575 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
576 (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.0915e-05 ^2 * pt^2 + 0.00663766^2 ) +
[7a1e27d]577 (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.15518e-05^2 * pt^2 + 0.00398644^2 ) +
578 (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2 + 0.00317807^2 ) +
[8c09a6b9]579 (abs(eta) < 0.18) * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2 )
[7a1e27d]580
581 }
582}
583
584###############################
585# Momentum resolution for muons
586###############################
587
588module MomentumSmearing MuonMomentumSmearing {
589 set InputArray MuonTrackingEfficiency/muons
590 set OutputArray muons
591
592 # Resolution given in dpT/pT.
593
594 # CLICdet internal studies
595 set ResolutionFormula {
596
[5ec2b23]597 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
[6f96f62]598 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
599 (abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt(9.81626e-06^2 * pt^2 + 0.00379895^2 ) +
[f9aa59d]600 (abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt(1.1959e-05^2 * pt^2 + 0.00242417^2 ) +
[7a1e27d]601 (abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2 + 0.00219291^2 ) +
602 (abs(eta) < 0.18) * sqrt(1.29686e-05^2 * pt^2 + 0.0020392^2 )
603
604 }
605}
606
607##############
608# Track merger
609##############
610
611module Merger TrackMerger {
612 # add InputArray InputArray
613 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
614 add InputArray ElectronMomentumSmearing/electrons
615 add InputArray MuonMomentumSmearing/muons
616 set OutputArray tracks
617}
618
619#############
620# ECAL
621#############
622
623module SimpleCalorimeter ECal {
624 set ParticleInputArray ParticlePropagator/stableParticles
625 set TrackInputArray TrackMerger/tracks
626
627 set TowerOutputArray ecalTowers
628 set EFlowTrackOutputArray eflowTracks
629 set EFlowTowerOutputArray eflowPhotons
630
631 set IsEcal true
632
633 set EnergyMin 0.5
634 set EnergySignificanceMin 1.0
635
636 set SmearTowerCenter true
637
638 set pi [expr {acos(-1)}]
639
640 # lists of the edges of each tower in eta and phi
641 # each list starts with the lower edge of the first tower
642 # the list ends with the higher edged of the last tower
643
644 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
645 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
646 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3
647 #ECAL cell sizes always 5x5 mm^2
648
649 #barrel:
650 #dphi = 0.2 degree towers up to eta <=1.2
651 set PhiBins {}
652 for {set i -900} {$i <= 900} {incr i} {
653 add PhiBins [expr {$i * $pi/900.0 }]
654 }
655 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
656 for {set i -400} {$i <=400} {incr i} {
657 set eta [expr {$i * 0.003}]
658 add EtaPhiBins $eta $PhiBins
659 }
660
661 #endcaps:
662 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
663 set PhiBins {}
664 for {set i -225} {$i <= 225} {incr i} {
665 add PhiBins [expr {$i * $pi/225.}]
666 }
667 #deta=0.02 units for 1.2 < |eta| <=2.5
668 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
669 for {set i 1} {$i <=66} {incr i} {
670 set eta [expr {-2.52 + $i * 0.02}]
671 add EtaPhiBins $eta $PhiBins
672 }
673 #same for 1.2 to 2.5
674 for {set i 1} {$i <=66} {incr i} {
675 set eta [expr {1.18 + $i*0.02}]
676 add EtaPhiBins $eta $PhiBins
677 }
678
679 #plug:
680 #dphi = 1 degree for 2.5 < eta <=3
681 set PhiBins {}
682 for {set i -180} {$i <= 180} {incr i} {
683 add PhiBins [expr {$i * $pi/180.}]
684 }
685 # deta = 0.02 for 2.5 < |eta| <=3
686 # from -3 to -2.5, there will be 25 segments
687 for {set i 1} {$i <= 26} {incr i} {
688 set eta [expr {-3.02 + $i * 0.02}]
689 add EtaPhiBins $eta $PhiBins
690 }
691 #same for 2.5 to 3
692 for {set i 1} {$i <= 26} {incr i} {
693 set eta [expr {2.48 + $i*0.02}]
694 add EtaPhiBins $eta $PhiBins
695 }
696
697
698
699 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
700
701 add EnergyFraction {0} {0.0}
702 # energy fractions for e, gamma and pi0
703 add EnergyFraction {11} {1.0}
704 add EnergyFraction {22} {1.0}
705 add EnergyFraction {111} {1.0}
706 # energy fractions for muon, neutrinos and neutralinos
707 add EnergyFraction {12} {0.0}
708 add EnergyFraction {13} {0.0}
709 add EnergyFraction {14} {0.0}
710 add EnergyFraction {16} {0.0}
711 add EnergyFraction {1000022} {0.0}
712 add EnergyFraction {1000023} {0.0}
713 add EnergyFraction {1000025} {0.0}
714 add EnergyFraction {1000035} {0.0}
715 add EnergyFraction {1000045} {0.0}
716 # energy fractions for K0short and Lambda
717 add EnergyFraction {310} {0.3}
718 add EnergyFraction {3122} {0.3}
719
720 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
721 set ResolutionFormula {
722 (abs(eta) <= 0.78 ) * sqrt(energy^2*0.01^2 + energy*0.156^2)+
723 (abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2 ) +
724 (abs(eta) <= 3 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2 )}
725}
726
727#############
728# HCAL
729#############
730
731module SimpleCalorimeter HCal {
732 set ParticleInputArray ParticlePropagator/stableParticles
733 set TrackInputArray ECal/eflowTracks
734
735 set TowerOutputArray hcalTowers
736 set EFlowTrackOutputArray eflowTracks
737 set EFlowTowerOutputArray eflowNeutralHadrons
738
739 set IsEcal false
740
741 set EnergyMin 1.0
742 set EnergySignificanceMin 1.0
743
744 set SmearTowerCenter true
745
746 set pi [expr {acos(-1)}]
747
748 # lists of the edges of each tower in eta and phi
749 # each list starts with the lower edge of the first tower
750 # the list ends with the higher edged of the last tower
751
752
753 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
754 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
755 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5
756 #HCAL cell sizes always 30x30 mm^2
757
758 #barrel and ring:
759 #dphi = 1 degree up to |eta| <=0.9
760 set PhiBins {}
761 for {set i -180} {$i <=180} {incr i} {
762 add PhiBins [expr {$i * $pi/180.0}]
763 }
764 #deta= 0.02 towers up to |eta| <=0.9
765 for {set i -45} {$i <=45} {incr i} {
766 set eta [expr {$i * 0.02}]
767 add EtaPhiBins $eta $PhiBins
768 }
769
770 #endcaps:
771 # dphi = 6 degree
772 set PhiBins {}
773 for {set i -30} {$i <=30} {incr i} {
774 add PhiBins [expr {$i * $pi/30.0}]
775 }
776 # deta =0.1 for 0.9 < |eta| <=3.0
777 #for -3.0 to -0.9, 21 segments
778 for {set i 1} {$i <=22} {incr i} {
779 set eta [expr {-3.1 + $i * 0.1}]
780 add EtaPhiBins $eta $PhiBins
781 }
782 #same for 0.9 to 3.0
783 for {set i 1} {$i <=22} {incr i} {
784 set eta [expr {0.8 + $i * 0.1 }]
785 add EtaPhiBins $eta $PhiBins
786 }
787
788 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
789 add EnergyFraction {0} {1.0}
790 # energy fractions for e, gamma and pi0
791 add EnergyFraction {11} {0.0}
792 add EnergyFraction {22} {0.0}
793 add EnergyFraction {111} {0.0}
794 # energy fractions for muon, neutrinos and neutralinos
795 add EnergyFraction {12} {0.0}
796 add EnergyFraction {13} {0.0}
797 add EnergyFraction {14} {0.0}
798 add EnergyFraction {16} {0.0}
799 add EnergyFraction {1000022} {0.0}
800 add EnergyFraction {1000023} {0.0}
801 add EnergyFraction {1000025} {0.0}
802 add EnergyFraction {1000035} {0.0}
803 add EnergyFraction {1000045} {0.0}
804 # energy fractions for K0short and Lambda
805 add EnergyFraction {310} {0.7}
806 add EnergyFraction {3122} {0.7}
807
808 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
809 #CLICdet internal studies
810 set ResolutionFormula {
811 (abs(eta)<= 0.3) * sqrt(1.38^2 + energy*0.308^2 + energy^2*0.050^2) +
812 (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2 + energy*0.322^2 + energy^2*0.048^2) +
813 (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt( 1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) +
814 (abs(eta)<=3 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2 )
815 }
816
817}
818
819#################
820# Electron filter
821#################
822
823module PdgCodeFilter ElectronFilter {
824 set InputArray HCal/eflowTracks
825 set OutputArray electrons
826 set Invert true
827 add PdgCode {11}
828 add PdgCode {-11}
829}
830
831######################
832# ChargedHadronFilter
833######################
834
835module PdgCodeFilter ChargedHadronFilter {
836 set InputArray HCal/eflowTracks
837 set OutputArray chargedHadrons
838
839 add PdgCode {11}
840 add PdgCode {-11}
841 add PdgCode {13}
842 add PdgCode {-13}
843}
844
845
846
847###################################################
848# Tower Merger (in case not using e-flow algorithm)
849###################################################
850
851module Merger Calorimeter {
852 # add InputArray InputArray
853 add InputArray ECal/ecalTowers
854 add InputArray HCal/hcalTowers
855 set OutputArray towers
856}
857
858
859####################
860# Energy flow merger
861####################
862
863module Merger EFlowMerger {
864 # add InputArray InputArray
865 add InputArray HCal/eflowTracks
866 add InputArray ECal/eflowPhotons
867 add InputArray HCal/eflowNeutralHadrons
868 set OutputArray eflow
869}
870
871
872###################
873# Photon efficiency
874###################
875
876module Efficiency PhotonEfficiency {
877 set InputArray ECal/eflowPhotons
878 set OutputArray photons
879
880 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
881
882 # efficiency formula for photons
883 # current full simulation of CLICdet yields:
884 set EfficiencyFormula {
885 (energy < 2.0 ) * (0.000) +
886 (energy >= 2.0) * (abs(eta) < 0.7)*(0.94) +
887 (energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=3.0) * (0.9) }
888
889}
890
[505a779]891
[7a1e27d]892##################
893# Photon isolation
894##################
895
896module Isolation PhotonIsolation {
897 set CandidateInputArray PhotonEfficiency/photons
898 set IsolationInputArray EFlowMerger/eflow
899
900 set OutputArray photons
901
[505a779]902 set DeltaRMax 0.1
[7a1e27d]903
904 set PTMin 0.5
905
[505a779]906 set PTRatioMax 0.2
[7a1e27d]907}
908
[505a779]909
[7a1e27d]910#####################
911# Electron efficiency
912#####################
913
914module Efficiency ElectronEfficiency {
915 set InputArray ElectronFilter/electrons
916 set OutputArray electrons
917
918 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
919
920 set EfficiencyFormula {
921 (energy < 3.0) * ( 0.00 ) +
922 ( energy >=3 && energy < 8 ) * (abs(eta) > 1.95) * (0.58 ) +
923 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) +
924 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) +
925 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) +
926 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) +
927 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.69) * (0.84 ) +
928 ( energy >=8 && energy < 13 ) * (abs(eta) > 1.95) * ( 0.6 ) +
929 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) +
930 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) +
931 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) +
932 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) +
933 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.69) * ( 0.88 ) +
934 ( energy >=13 && energy < 18 ) * (abs(eta) > 1.95) * ( 0.6 ) +
935 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.8 ) +
936 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) +
937 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) +
938 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) +
939 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.69) * ( 0.9 ) +
940 ( energy >=18 && energy < 23 ) * (abs(eta) > 1.95) * (0.64 ) +
941 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) +
942 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) +
943 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) +
944 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) +
945 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.69) * (0.92 ) +
946 ( energy >= 23 && energy < 28 ) * (abs(eta) > 1.95) * (0.64 ) +
947 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) +
948 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) +
949 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) +
950 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) +
951 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.69) * (0.94 ) +
952 ( energy >=28 && energy < 35 ) * (abs(eta) > 1.95) * (0.67 ) +
953 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) +
954 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) +
955 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) +
956 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
957 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.69) * (0.94 ) +
958 ( energy >=35 && energy < 45 ) * (abs(eta) > 1.95) * (0.68 ) +
959 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) +
960 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) +
961 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) +
962 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
963 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.69) * (0.96 ) +
964 ( energy >=45 && energy < 80 ) * (abs(eta) > 1.95) * ( 0.7 ) +
965 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) +
966 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.8 ) +
967 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) +
968 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) +
969 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.69) * ( 0.97 ) +
970 ( energy >=80 && energy < 200 ) * (abs(eta) > 1.95) * (0.68 ) +
971 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
972 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) +
973 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) +
974 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
975 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.69) * (0.98 ) +
976 ( energy >=200 && energy < 400 ) * (abs(eta) > 1.95) * ( 0.68 ) +
977 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) +
978 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) +
979 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) +
980 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) +
981 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.69) * ( 0.98 ) +
982 ( energy >=400 ) * (abs(eta) > 1.95) * (0.68 ) +
983 ( energy >=400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
984 ( energy >=400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) +
985 ( energy >=400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) +
986 ( energy >=400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
987 ( energy >=400 ) * (abs(eta) <= 0.69) * (0.98 )
988 }
989}
990
991####################
992# Electron isolation
993####################
994
995module Isolation ElectronIsolation {
996 set CandidateInputArray ElectronEfficiency/electrons
997 set IsolationInputArray EFlowMerger/eflow
998
999 set OutputArray electrons
1000
[505a779]1001 set DeltaRMax 0.1
[7a1e27d]1002
1003 set PTMin 0.5
1004
[505a779]1005 set PTRatioMax 0.2
[7a1e27d]1006}
1007
1008#################
1009# Muon efficiency
1010#################
1011
1012module Efficiency MuonEfficiency {
1013 set InputArray MuonMomentumSmearing/muons
1014 set OutputArray muons
1015
1016 # set EfficiencyFormula {efficiency as a function of eta and pt}
1017
1018 # efficiency formula for muons
1019 # current full simulation of CLICdet yields:
1020
1021
1022
1023 set EfficiencyFormula {
[8c09a6b9]1024 (energy < 2.0 ) * (0.00) +
1025 (energy>=2.0) * (0.999)
[7a1e27d]1026 }
1027}
1028
1029################
1030# Muon isolation
1031################
1032
1033module Isolation MuonIsolation {
1034 set CandidateInputArray MuonEfficiency/muons
1035 set IsolationInputArray EFlowMerger/eflow
1036
1037 set OutputArray muons
1038
[505a779]1039 set DeltaRMax 0.1
[7a1e27d]1040
1041 set PTMin 0.5
1042
[505a779]1043 set PTRatioMax 0.2
[7a1e27d]1044}
1045
1046
1047
1048###################
1049# Missing ET merger
1050###################
1051
1052module Merger MissingET {
1053 # add InputArray InputArray
1054 add InputArray EFlowMerger/eflow
1055 set MomentumOutputArray momentum
1056}
1057
1058
1059##################
1060# Scalar HT merger
1061##################
1062
1063module Merger ScalarHT {
1064 # add InputArray InputArray
1065 add InputArray EFlowMerger/eflow
1066 set EnergyOutputArray energy
1067}
1068######################
1069# EFlowFilter (UniqueObjectFinder)
1070######################
1071module UniqueObjectFinder EFlowFilter {
1072 add InputArray PhotonIsolation/photons photons
1073 add InputArray ElectronIsolation/electrons electrons
1074 add InputArray MuonIsolation/muons muons
1075 add InputArray EFlowMerger/eflow eflow
1076}
1077
1078#################
1079# Neutrino Filter
1080#################
1081
1082module PdgCodeFilter NeutrinoFilter {
1083
1084 set InputArray Delphes/stableParticles
1085 set OutputArray filteredParticles
1086
1087 set PTMin 0.0
1088
1089 add PdgCode {12}
1090 add PdgCode {14}
1091 add PdgCode {16}
1092 add PdgCode {-12}
1093 add PdgCode {-14}
1094 add PdgCode {-16}
1095
1096}
1097
1098
1099#####################
1100# MC truth jet finder
1101#####################
1102
1103module FastJetFinder GenJetFinder {
1104 set InputArray NeutrinoFilter/filteredParticles
1105
1106 set OutputArray jets
1107
1108 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
1109 set JetAlgorithm 9
1110 set ParameterR 0.5
1111
1112 set JetPTMin 20.0
1113}
1114
1115#########################
1116# Gen Missing ET merger
1117########################
1118
1119module Merger GenMissingET {
1120 # add InputArray InputArray
1121 add InputArray NeutrinoFilter/filteredParticles
1122 set MomentumOutputArray momentum
1123}
1124
1125
1126
1127############
1128# Jet finder
1129############
1130
1131module FastJetFinder FastJetFinderKt {
1132 # set InputArray Calorimeter/towers
[eee976c2]1133 set InputArray EFlowFilter/eflow
[7a1e27d]1134
1135 set OutputArray KTjets
1136
1137 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
1138 set JetAlgorithm 4
1139 set ParameterR 0.5
1140
1141 set JetPTMin 20.0
1142}
1143
1144################
1145# Jet finder VLC
1146################
1147
1148
[96fd15f]1149source CLIC/CLICdet_JetReco.tcl
[7a1e27d]1150
1151
1152#########################################
1153# Jet Momentum Smearing to mimick overlay
1154#########################################
1155
1156
[96fd15f]1157source CLIC/CLICdet_JetSmearing.tcl
[7a1e27d]1158
1159
1160
1161########################
1162# Jet Flavor Association
1163########################
1164
[96fd15f]1165source CLIC/CLICdet_JetFlavorAssociation.tcl
[0e6b14e]1166source CLIC/CLICdet_JetFlavorAssociation_JER.tcl
[7a1e27d]1167
1168###########
1169# b-tagging
1170###########
1171# based on CLICdp-Note-2014-002
1172
[96fd15f]1173source CLIC/CLICdet_BTagging.tcl
[0e6b14e]1174source CLIC/CLICdet_BTagging_JER.tcl
[7a1e27d]1175
1176
1177#############
1178# tau-tagging
1179#############
1180# based on LCD-2010-009
1181
[96fd15f]1182source CLIC/CLICdet_TauTagging.tcl
[0e6b14e]1183source CLIC/CLICdet_TauTagging_JER.tcl
[7a1e27d]1184
1185
1186
1187##################
1188# ROOT tree writer
1189##################
1190
1191module TreeWriter TreeWriter {
1192 # add Branch InputArray BranchName BranchClass
1193 add Branch Delphes/allParticles Particle GenParticle
1194
1195 add Branch GenJetFinder/jets GenJet Jet
1196
1197 add Branch FastJetFinderKt/KTjets KTjet Jet
1198 add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet
1199 add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet
1200 add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet
1201 add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet
1202 add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet
1203
1204 add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet
1205 add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet
1206 add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet
1207 add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet
1208 add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet
1209
1210 add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet
1211 add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet
1212 add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet
1213 add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet
1214 add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet
1215
1216 add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet
1217 add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet
1218 add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet
1219 add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet
1220 add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet
1221
1222 add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet
1223 add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet
1224 add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet
1225 add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet
1226 add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
1227
1228 add Branch FastJetFinderVLC_R05_inclusive/VLCjetsR05_inclusive VLCjetR05_inclusive Jet
1229 add Branch FastJetFinderVLC_R07_inclusive/VLCjetsR07_inclusive VLCjetR07_inclusive Jet
1230 add Branch FastJetFinderVLC_R10_inclusive/VLCjetsR10_inclusive VLCjetR10_inclusive Jet
1231 add Branch FastJetFinderVLC_R12_inclusive/VLCjetsR12_inclusive VLCjetR12_inclusive Jet
1232 add Branch FastJetFinderVLC_R15_inclusive/VLCjetsR15_inclusive VLCjetR15_inclusive Jet
1233
1234
1235
1236
[a89c86c]1237 ###with jet energy smearing
1238 add Branch JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2 JER_VLCjetR05N2 Jet
1239 add Branch JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3 JER_VLCjetR05N3 Jet
1240 add Branch JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4 JER_VLCjetR05N4 Jet
1241 add Branch JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5 JER_VLCjetR05N5 Jet
1242 add Branch JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6 JER_VLCjetR05N6 Jet
1243
1244 add Branch JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2 JER_VLCjetR07N2 Jet
1245 add Branch JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3 JER_VLCjetR07N3 Jet
1246 add Branch JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4 JER_VLCjetR07N4 Jet
1247 add Branch JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5 JER_VLCjetR07N5 Jet
1248 add Branch JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6 JER_VLCjetR07N6 Jet
1249
1250 add Branch JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2 JER_VLCjetR10N2 Jet
1251 add Branch JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3 JER_VLCjetR10N3 Jet
1252 add Branch JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4 JER_VLCjetR10N4 Jet
1253 add Branch JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5 JER_VLCjetR10N5 Jet
1254 add Branch JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6 JER_VLCjetR10N6 Jet
1255
1256 add Branch JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2 JER_VLCjetR12N2 Jet
1257 add Branch JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3 JER_VLCjetR12N3 Jet
1258 add Branch JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4 JER_VLCjetR12N4 Jet
1259 add Branch JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5 JER_VLCjetR12N5 Jet
1260 add Branch JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6 JER_VLCjetR12N6 Jet
1261
1262 add Branch JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2 JER_VLCjetR15N2 Jet
1263 add Branch JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3 JER_VLCjetR15N3 Jet
1264 add Branch JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4 JER_VLCjetR15N4 Jet
1265 add Branch JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5 JER_VLCjetR15N5 Jet
1266 add Branch JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6 JER_VLCjetR15N6 Jet
1267
1268 add Branch JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive JER_VLCjetR05_inclusive Jet
1269 add Branch JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive JER_VLCjetR07_inclusive Jet
1270 add Branch JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive JER_VLCjetR10_inclusive Jet
1271 add Branch JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive JER_VLCjetR12_inclusive Jet
1272 add Branch JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive JER_VLCjetR15_inclusive Jet
[7a1e27d]1273
1274 ####
1275
1276 add Branch GenMissingET/momentum GenMissingET MissingET
1277
1278 add Branch TrackMerger/tracks Track Track
1279 add Branch Calorimeter/towers Tower Tower
1280
1281 add Branch HCal/eflowTracks EFlowTrack Track
1282 add Branch ECal/eflowPhotons EFlowPhoton Tower
1283 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1284
1285 add Branch EFlowFilter/photons Photon Photon
1286 add Branch EFlowFilter/electrons Electron Electron
1287 add Branch EFlowFilter/muons Muon Muon
1288
1289 add Branch MissingET/momentum MissingET MissingET
1290 add Branch ScalarHT/energy ScalarHT ScalarHT
1291}
1292
Note: See TracBrowser for help on using the repository browser.