Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage2.tcl@ a26a130

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