Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage3_fcal.tcl@ 9cc5aeb

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