Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage1.tcl@ 4c9835d

Last change on this file since 4c9835d 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: 32.6 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#
[5a1f24a]8# For the low energy stage of
9# CLIC: 380 (350) GeV
10#
[7a1e27d]11# Jet finding with Valencia algorithm:
12# use exclusive clustering with njets
13# according to final state
14#
[505a779]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
99
100 JetFlavorAssociation_R05N2
101 JetFlavorAssociation_R05N3
102 JetFlavorAssociation_R05N4
103 JetFlavorAssociation_R05N5
104 JetFlavorAssociation_R05N6
105
106 JetFlavorAssociation_R07N2
107 JetFlavorAssociation_R07N3
108 JetFlavorAssociation_R07N4
109 JetFlavorAssociation_R07N5
110 JetFlavorAssociation_R07N6
111
112 JetFlavorAssociation_R10N2
113 JetFlavorAssociation_R10N3
114 JetFlavorAssociation_R10N4
115 JetFlavorAssociation_R10N5
116 JetFlavorAssociation_R10N6
117
118 JetFlavorAssociation_R12N2
119 JetFlavorAssociation_R12N3
120 JetFlavorAssociation_R12N4
121 JetFlavorAssociation_R12N5
122 JetFlavorAssociation_R12N6
123
124 JetFlavorAssociation_R15N2
125 JetFlavorAssociation_R15N3
126 JetFlavorAssociation_R15N4
127 JetFlavorAssociation_R15N5
128 JetFlavorAssociation_R15N6
129
130 JetFlavorAssociation_R05_inclusive
131 JetFlavorAssociation_R07_inclusive
132 JetFlavorAssociation_R10_inclusive
133 JetFlavorAssociation_R12_inclusive
134 JetFlavorAssociation_R15_inclusive
135
136
137 BTaggingWP50_R05N2
138 BTaggingWP70_R05N2
139 BTaggingWP90_R05N2
140 BTaggingWP50_R05N3
141 BTaggingWP70_R05N3
142 BTaggingWP90_R05N3
143 BTaggingWP50_R05N4
144 BTaggingWP70_R05N4
145 BTaggingWP90_R05N4
146 BTaggingWP50_R05N5
147 BTaggingWP70_R05N5
148 BTaggingWP90_R05N5
149 BTaggingWP50_R05N6
150 BTaggingWP70_R05N6
151 BTaggingWP90_R05N6
152 BTaggingWP50_R07N2
153 BTaggingWP70_R07N2
154 BTaggingWP90_R07N2
155 BTaggingWP50_R07N3
156 BTaggingWP70_R07N3
157 BTaggingWP90_R07N3
158 BTaggingWP50_R07N4
159 BTaggingWP70_R07N4
160 BTaggingWP90_R07N4
161 BTaggingWP50_R07N5
162 BTaggingWP70_R07N5
163 BTaggingWP90_R07N5
164 BTaggingWP50_R07N6
165 BTaggingWP70_R07N6
166 BTaggingWP90_R07N6
167 BTaggingWP50_R10N2
168 BTaggingWP70_R10N2
169 BTaggingWP90_R10N2
170 BTaggingWP50_R10N3
171 BTaggingWP70_R10N3
172 BTaggingWP90_R10N3
173 BTaggingWP50_R10N4
174 BTaggingWP70_R10N4
175 BTaggingWP90_R10N4
176 BTaggingWP50_R10N5
177 BTaggingWP70_R10N5
178 BTaggingWP90_R10N5
179 BTaggingWP50_R10N6
180 BTaggingWP70_R10N6
181 BTaggingWP90_R10N6
182 BTaggingWP50_R12N2
183 BTaggingWP70_R12N2
184 BTaggingWP90_R12N2
185 BTaggingWP50_R12N3
186 BTaggingWP70_R12N3
187 BTaggingWP90_R12N3
188 BTaggingWP50_R12N4
189 BTaggingWP70_R12N4
190 BTaggingWP90_R12N4
191 BTaggingWP50_R12N5
192 BTaggingWP70_R12N5
193 BTaggingWP90_R12N5
194 BTaggingWP50_R12N6
195 BTaggingWP70_R12N6
196 BTaggingWP90_R12N6
197 BTaggingWP50_R15N2
198 BTaggingWP70_R15N2
199 BTaggingWP90_R15N2
200 BTaggingWP50_R15N3
201 BTaggingWP70_R15N3
202 BTaggingWP90_R15N3
203 BTaggingWP50_R15N4
204 BTaggingWP70_R15N4
205 BTaggingWP90_R15N4
206 BTaggingWP50_R15N5
207 BTaggingWP70_R15N5
208 BTaggingWP90_R15N5
209 BTaggingWP50_R15N6
210 BTaggingWP70_R15N6
211 BTaggingWP90_R15N6
212 BTaggingWP50_R05_inclusive
213 BTaggingWP70_R05_inclusive
214 BTaggingWP90_R05_inclusive
215 BTaggingWP50_R07_inclusive
216 BTaggingWP70_R07_inclusive
217 BTaggingWP90_R07_inclusive
218 BTaggingWP50_R10_inclusive
219 BTaggingWP70_R10_inclusive
220 BTaggingWP90_R10_inclusive
221 BTaggingWP50_R12_inclusive
222 BTaggingWP70_R12_inclusive
223 BTaggingWP90_R12_inclusive
224 BTaggingWP50_R15_inclusive
225 BTaggingWP70_R15_inclusive
226 BTaggingWP90_R15_inclusive
227
228
229 TauTagging_R05N2
230 TauTagging_R05N3
231 TauTagging_R05N4
232 TauTagging_R05N5
233 TauTagging_R05N6
234 TauTagging_R07N2
235 TauTagging_R07N3
236 TauTagging_R07N4
237 TauTagging_R07N5
238 TauTagging_R07N6
239 TauTagging_R10N2
240 TauTagging_R10N3
241 TauTagging_R10N4
242 TauTagging_R10N5
243 TauTagging_R10N6
244 TauTagging_R12N2
245 TauTagging_R12N3
246 TauTagging_R12N4
247 TauTagging_R12N5
248 TauTagging_R12N6
249 TauTagging_R15N2
250 TauTagging_R15N3
251 TauTagging_R15N4
252 TauTagging_R15N5
253 TauTagging_R15N6
254
255 TauTagging_R05_inclusive
256 TauTagging_R07_inclusive
257 TauTagging_R10_inclusive
258 TauTagging_R12_inclusive
259 TauTagging_R15_inclusive
260
261
262
263 ScalarHT
264
265
266 TreeWriter
267}
268
269#################################
270# Propagate particles in cylinder
271#################################
272
273module ParticlePropagator ParticlePropagator {
274 set InputArray Delphes/stableParticles
275
276 set OutputArray stableParticles
277 set ChargedHadronOutputArray chargedHadrons
278 set ElectronOutputArray electrons
279 set MuonOutputArray muons
280
281 # radius of the magnetic field coverage in the calorimeter, in m
282 set Radius 1.5
283 # half-length of the magnetic field coverage in the calorimeter, in m
284 set HalfLength 2.31
285
286 # magnetic field, in T
287 set Bz 4.0
288}
289
290####################################
291# Charged hadron tracking efficiency
292####################################
293
294module Efficiency ChargedHadronTrackingEfficiency {
295 set InputArray ParticlePropagator/chargedHadrons
296 set OutputArray chargedHadrons
297 # Current full simulation with CLICdet provides for pions:
298
299 set EfficiencyFormula {
[8d7319c]300 (pt <= 0.1) * (0.000) +
301 (pt > 0.1) * (abs(eta) > 2.54) * (0.000) +
302 (pt > 0.1) * (energy >= 80) * (abs(eta) < 2.54) * (1.000) +
303 (pt > 0.1) * (energy < 80 && energy >= 3) * (abs(eta) <=2.54 && abs(eta) > 2.34) * (0.994) +
304 (pt > 0.1) * (energy < 80 && energy >= 3) * (abs(eta) <= 2.34) * (1.000) +
305 (pt > 0.1) * (energy < 3) * (abs(eta) <= 2.54 && abs(eta) > 0.55 ) * (0.990) +
306 (pt > 0.1) * (energy < 3) * (abs(eta) <= 0.55 ) * (1.000)
[7a1e27d]307 }
308}
309
310##############################
311# Electron tracking efficiency
312##############################
313
314module Efficiency ElectronTrackingEfficiency {
315 set InputArray ParticlePropagator/electrons
316 set OutputArray electrons
317
318
319 # Current full simulation with CLICdet provides for electrons:
320 set EfficiencyFormula {
[8d7319c]321 (pt <= 0.1) * (0.000) +
322 (pt > 0.1) * (abs(eta) > 2.54) * (0.000) +
323 (pt > 0.1) * (energy >= 80) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (0.993) +
324 (pt > 0.1) * (energy >= 80) * (abs(eta) <= 2.44 && abs(eta) > 2.34 ) * (0.997) +
325 (pt > 0.1) * (energy >= 80) * (abs(eta) <= 2.34 ) * (1.000) +
326 (pt > 0.1) * (energy < 80 && energy >= 5) * (abs(eta) <= 2.54 && abs(eta) > 2.17 ) * (0.998) +
327 (pt > 0.1) * (energy < 80 && energy >= 5) * (abs(eta) <= 2.17) * (1.000) +
328 (pt > 0.1) * (energy < 5) * (abs(eta) <= 2.54 && abs(eta) > 2.34 ) * (1.000) +
329 (pt > 0.1) * (energy < 5) * (abs(eta) <= 2.34 && abs(eta) > 0.76 ) * (0.997) +
330 (pt > 0.1) * (energy < 5) * (abs(eta) <= 0.76) * (0.999)
[7a1e27d]331 }
332}
333
334##########################
335# Muon tracking efficiency
336##########################
337
338module Efficiency MuonTrackingEfficiency {
339 set InputArray ParticlePropagator/muons
340 set OutputArray muons
341
342 # Current full simulation with CLICdet provides for muons:
343 set EfficiencyFormula {
[8d7319c]344 (pt <= 0.1) * (0.000) +
345 (pt > 0.1) * (abs(eta) > 2.54) * (0.000) +
346 (pt > 0.1) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 80) * (0.994) +
347 (pt > 0.1) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy >= 5 && energy < 80) * (0.996) +
348 (pt > 0.1) * (abs(eta) <= 2.54 && abs(eta) > 2.44 ) * (energy < 5 ) * (0.996) +
349 (pt > 0.1) * (abs(eta) <= 2.44 ) * (energy >= 5 ) * (1.000) +
350 (pt > 0.1) * (abs(eta) <= 2.44 && abs(eta) > 2.25 ) * (energy < 5 ) * (0.999) +
351 (pt > 0.1) * (abs(eta) <= 2.25 ) * (energy < 5 ) * (1.000)
[7a1e27d]352 }
353}
354
355########################################
356# Momentum resolution for charged tracks
357########################################
358
359module MomentumSmearing ChargedHadronMomentumSmearing {
360 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
361 set OutputArray chargedHadrons
362
363
364 # Resolution given in dpT/pT.
365 # CLICdet internal studies
366 set ResolutionFormula {
367 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
[6f96f62]368 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2 ) +
369 (abs(eta) < 1.32 && abs(eta) >= 0.88 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2 ) +
370 (abs(eta) < 0.88 && abs(eta) >= 0.45 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2 ) +
371 (abs(eta) < 0.45 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2 ) +
[7a1e27d]372 (abs(eta) < 0.18) * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2 )
373
374 }
375}
376
377###################################
378# Momentum resolution for electrons
379###################################
380
381module MomentumSmearing ElectronMomentumSmearing {
382 set InputArray ElectronTrackingEfficiency/electrons
383 set OutputArray electrons
384
385 # Resolution given in dpT/pT.
386 # CLICdet internal studies
387 set ResolutionFormula {
388 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
[6f96f62]389 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * sqrt( 8.62283e-05^2 * pt^2 + 0.0177556^2 ) +
390 (abs(eta) < 1.32 && abs(eta) >= 0.88 ) * sqrt( 1.0915e-05 ^2 * pt^2 + 0.00663766^2 ) +
391 (abs(eta) < 0.88 && abs(eta) >= 0.45 ) * sqrt( 1.15518e-05^2 * pt^2 + 0.00398644^2 ) +
392 (abs(eta) < 0.45 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2 + 0.00317807^2 ) +
[7a1e27d]393 (abs(eta) < 0.18) * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2 )
394
395 }
396}
397
398###############################
399# Momentum resolution for muons
400###############################
401
402module MomentumSmearing MuonMomentumSmearing {
403 set InputArray MuonTrackingEfficiency/muons
404 set OutputArray muons
405
406 # Resolution given in dpT/pT.
407
408 # CLICdet internal studies
409 set ResolutionFormula {
410
[5ec2b23]411 (abs(eta) < 2.66 && abs(eta) >= 1.74 ) * 2 * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
[6f96f62]412 (abs(eta) < 1.74 && abs(eta) >= 1.32 ) * sqrt(4.57439e-05^2 * pt^2 + 0.0149328^2 ) +
413 (abs(eta) < 1.32 && abs(eta) >= 0.88 ) * sqrt(9.81626e-06^2 * pt^2 + 0.00379895^2 ) +
414 (abs(eta) < 0.88 && abs(eta) >= 0.45 ) * sqrt(1.1959e-05^2 * pt^2 + 0.00242417^2 ) +
415 (abs(eta) < 0.45 && abs(eta) >= 0.18 ) * sqrt(1.20149e-05^2 * pt^2 + 0.00219291^2 ) +
[7a1e27d]416 (abs(eta) < 0.18) * sqrt(1.29686e-05^2 * pt^2 + 0.0020392^2 )
417
418 }
419}
420
421##############
422# Track merger
423##############
424
425module Merger TrackMerger {
426 # add InputArray InputArray
427 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
428 add InputArray ElectronMomentumSmearing/electrons
429 add InputArray MuonMomentumSmearing/muons
430 set OutputArray tracks
431}
432
433#############
434# ECAL
435#############
436
437module SimpleCalorimeter ECal {
438 set ParticleInputArray ParticlePropagator/stableParticles
439 set TrackInputArray TrackMerger/tracks
440
441 set TowerOutputArray ecalTowers
442 set EFlowTrackOutputArray eflowTracks
443 set EFlowTowerOutputArray eflowPhotons
444
445 set IsEcal true
446
447 set EnergyMin 0.5
448 set EnergySignificanceMin 1.0
449
450 set SmearTowerCenter true
451
452 set pi [expr {acos(-1)}]
453
454 # lists of the edges of each tower in eta and phi
455 # each list starts with the lower edge of the first tower
456 # the list ends with the higher edged of the last tower
457
458 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
459 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
460 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3
461 #ECAL cell sizes always 5x5 mm^2
462
463 #barrel:
464 #dphi = 0.2 degree towers up to eta <=1.2
465 set PhiBins {}
466 for {set i -900} {$i <= 900} {incr i} {
467 add PhiBins [expr {$i * $pi/900.0 }]
468 }
469 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
470 for {set i -400} {$i <=400} {incr i} {
471 set eta [expr {$i * 0.003}]
472 add EtaPhiBins $eta $PhiBins
473 }
474
475 #endcaps:
476 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
477 set PhiBins {}
478 for {set i -225} {$i <= 225} {incr i} {
479 add PhiBins [expr {$i * $pi/225.}]
480 }
481 #deta=0.02 units for 1.2 < |eta| <=2.5
482 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
483 for {set i 1} {$i <=66} {incr i} {
484 set eta [expr {-2.52 + $i * 0.02}]
485 add EtaPhiBins $eta $PhiBins
486 }
487 #same for 1.2 to 2.5
488 for {set i 1} {$i <=66} {incr i} {
489 set eta [expr {1.18 + $i*0.02}]
490 add EtaPhiBins $eta $PhiBins
491 }
492
493 #plug:
494 #dphi = 1 degree for 2.5 < eta <=3
495 set PhiBins {}
496 for {set i -180} {$i <= 180} {incr i} {
497 add PhiBins [expr {$i * $pi/180.}]
498 }
499 # deta = 0.02 for 2.5 < |eta| <=3
500 # from -3 to -2.5, there will be 25 segments
501 for {set i 1} {$i <= 26} {incr i} {
502 set eta [expr {-3.02 + $i * 0.02}]
503 add EtaPhiBins $eta $PhiBins
504 }
505 #same for 2.5 to 3
506 for {set i 1} {$i <= 26} {incr i} {
507 set eta [expr {2.48 + $i*0.02}]
508 add EtaPhiBins $eta $PhiBins
509 }
510
511
512
513 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
514
515 add EnergyFraction {0} {0.0}
516 # energy fractions for e, gamma and pi0
517 add EnergyFraction {11} {1.0}
518 add EnergyFraction {22} {1.0}
519 add EnergyFraction {111} {1.0}
520 # energy fractions for muon, neutrinos and neutralinos
521 add EnergyFraction {12} {0.0}
522 add EnergyFraction {13} {0.0}
523 add EnergyFraction {14} {0.0}
524 add EnergyFraction {16} {0.0}
525 add EnergyFraction {1000022} {0.0}
526 add EnergyFraction {1000023} {0.0}
527 add EnergyFraction {1000025} {0.0}
528 add EnergyFraction {1000035} {0.0}
529 add EnergyFraction {1000045} {0.0}
530 # energy fractions for K0short and Lambda
531 add EnergyFraction {310} {0.3}
532 add EnergyFraction {3122} {0.3}
533
534 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
535 set ResolutionFormula {
536 (abs(eta) <= 0.78 ) * sqrt(energy^2*0.01^2 + energy*0.156^2)+
537 (abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2 ) +
538 (abs(eta) <= 3 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2 )}
539}
540
541#############
542# HCAL
543#############
544
545module SimpleCalorimeter HCal {
546 set ParticleInputArray ParticlePropagator/stableParticles
547 set TrackInputArray ECal/eflowTracks
548
549 set TowerOutputArray hcalTowers
550 set EFlowTrackOutputArray eflowTracks
551 set EFlowTowerOutputArray eflowNeutralHadrons
552
553 set IsEcal false
554
555 set EnergyMin 1.0
556 set EnergySignificanceMin 1.0
557
558 set SmearTowerCenter true
559
560 set pi [expr {acos(-1)}]
561
562 # lists of the edges of each tower in eta and phi
563 # each list starts with the lower edge of the first tower
564 # the list ends with the higher edged of the last tower
565
566
567 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
568 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
569 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5
570 #HCAL cell sizes always 30x30 mm^2
571
572 #barrel and ring:
573 #dphi = 1 degree up to |eta| <=0.9
574 set PhiBins {}
575 for {set i -180} {$i <=180} {incr i} {
576 add PhiBins [expr {$i * $pi/180.0}]
577 }
578 #deta= 0.02 towers up to |eta| <=0.9
579 for {set i -45} {$i <=45} {incr i} {
580 set eta [expr {$i * 0.02}]
581 add EtaPhiBins $eta $PhiBins
582 }
583
584 #endcaps:
585 # dphi = 6 degree
586 set PhiBins {}
587 for {set i -30} {$i <=30} {incr i} {
588 add PhiBins [expr {$i * $pi/30.0}]
589 }
590 # deta =0.1 for 0.9 < |eta| <=3.0
591 #for -3.0 to -0.9, 21 segments
592 for {set i 1} {$i <=22} {incr i} {
593 set eta [expr {-3.1 + $i * 0.1}]
594 add EtaPhiBins $eta $PhiBins
595 }
596 #same for 0.9 to 3.0
597 for {set i 1} {$i <=22} {incr i} {
598 set eta [expr {0.8 + $i * 0.1 }]
599 add EtaPhiBins $eta $PhiBins
600 }
601
602 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
603 add EnergyFraction {0} {1.0}
604 # energy fractions for e, gamma and pi0
605 add EnergyFraction {11} {0.0}
606 add EnergyFraction {22} {0.0}
607 add EnergyFraction {111} {0.0}
608 # energy fractions for muon, neutrinos and neutralinos
609 add EnergyFraction {12} {0.0}
610 add EnergyFraction {13} {0.0}
611 add EnergyFraction {14} {0.0}
612 add EnergyFraction {16} {0.0}
613 add EnergyFraction {1000022} {0.0}
614 add EnergyFraction {1000023} {0.0}
615 add EnergyFraction {1000025} {0.0}
616 add EnergyFraction {1000035} {0.0}
617 add EnergyFraction {1000045} {0.0}
618 # energy fractions for K0short and Lambda
619 add EnergyFraction {310} {0.7}
620 add EnergyFraction {3122} {0.7}
621
622 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
623 #CLICdet internal studies
624 set ResolutionFormula {
625 (abs(eta)<= 0.3) * sqrt(1.38^2 + energy*0.308^2 + energy^2*0.050^2) +
626 (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2 + energy*0.322^2 + energy^2*0.048^2) +
627 (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt( 1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) +
628 (abs(eta)<=3 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2 )
629 }
630
631}
632
633#################
634# Electron filter
635#################
636
637module PdgCodeFilter ElectronFilter {
638 set InputArray HCal/eflowTracks
639 set OutputArray electrons
640 set Invert true
641 add PdgCode {11}
642 add PdgCode {-11}
643}
644
645######################
646# ChargedHadronFilter
647######################
648
649module PdgCodeFilter ChargedHadronFilter {
650 set InputArray HCal/eflowTracks
651 set OutputArray chargedHadrons
652
653 add PdgCode {11}
654 add PdgCode {-11}
655 add PdgCode {13}
656 add PdgCode {-13}
657}
658
659
660
661###################################################
662# Tower Merger (in case not using e-flow algorithm)
663###################################################
664
665module Merger Calorimeter {
666 # add InputArray InputArray
667 add InputArray ECal/ecalTowers
668 add InputArray HCal/hcalTowers
669 set OutputArray towers
670}
671
672
673####################
674# Energy flow merger
675####################
676
677module Merger EFlowMerger {
678 # add InputArray InputArray
679 add InputArray HCal/eflowTracks
680 add InputArray ECal/eflowPhotons
681 add InputArray HCal/eflowNeutralHadrons
682 set OutputArray eflow
683}
684
685
686###################
687# Photon efficiency
688###################
689
690module Efficiency PhotonEfficiency {
691 set InputArray ECal/eflowPhotons
692 set OutputArray photons
693
694 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
695
696 # efficiency formula for photons
697 # current full simulation of CLICdet yields:
698 set EfficiencyFormula {
699 (energy < 2.0 ) * (0.000) +
700 (energy >= 2.0) * (abs(eta) < 0.7)*(0.94) +
701 (energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=3.0) * (0.9) }
702
703}
704
705##################
706# Photon isolation
707##################
708
709module Isolation PhotonIsolation {
710 set CandidateInputArray PhotonEfficiency/photons
711 set IsolationInputArray EFlowMerger/eflow
712
713 set OutputArray photons
714
[505a779]715 set DeltaRMax 0.1
[7a1e27d]716
717 set PTMin 0.5
718
[505a779]719 set PTRatioMax 0.2
[7a1e27d]720}
721
722#####################
723# Electron efficiency
724#####################
725
726module Efficiency ElectronEfficiency {
727 set InputArray ElectronFilter/electrons
728 set OutputArray electrons
729
730 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
731
732 set EfficiencyFormula {
733 (energy < 3.0) * ( 0.00 ) +
734 ( energy >=3 && energy < 8 ) * (abs(eta) > 1.95) * (0.58 ) +
735 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) +
736 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) +
737 ( energy >=3 && energy < 8 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) +
738 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) +
739 ( energy >=3 && energy < 8 ) * (abs(eta) <= 0.69) * (0.84 ) +
740 ( energy >=8 && energy < 13 ) * (abs(eta) > 1.95) * ( 0.6 ) +
741 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) +
742 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) +
743 ( energy >=8 && energy < 13 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) +
744 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) +
745 ( energy >=8 && energy < 13 ) * (abs(eta) <= 0.69) * ( 0.88 ) +
746 ( energy >=13 && energy < 18 ) * (abs(eta) > 1.95) * ( 0.6 ) +
747 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.8 ) +
748 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) +
749 ( energy >=13 && energy < 18 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) +
750 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) +
751 ( energy >=13 && energy < 18 ) * (abs(eta) <= 0.69) * ( 0.9 ) +
752 ( energy >=18 && energy < 23 ) * (abs(eta) > 1.95) * (0.64 ) +
753 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) +
754 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) +
755 ( energy >=18 && energy < 23 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) +
756 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) +
757 ( energy >=18 && energy < 23 ) * (abs(eta) <= 0.69) * (0.92 ) +
758 ( energy >= 23 && energy < 28 ) * (abs(eta) > 1.95) * (0.64 ) +
759 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) +
760 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) +
761 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) +
762 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) +
763 ( energy >= 23 && energy < 28 ) * (abs(eta) <= 0.69) * (0.94 ) +
764 ( energy >=28 && energy < 35 ) * (abs(eta) > 1.95) * (0.67 ) +
765 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) +
766 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) +
767 ( energy >=28 && energy < 35 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) +
768 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
769 ( energy >=28 && energy < 35 ) * (abs(eta) <= 0.69) * (0.94 ) +
770 ( energy >=35 && energy < 45 ) * (abs(eta) > 1.95) * (0.68 ) +
771 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) +
772 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) +
773 ( energy >=35 && energy < 45 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) +
774 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
775 ( energy >=35 && energy < 45 ) * (abs(eta) <= 0.69) * (0.96 ) +
776 ( energy >=45 && energy < 80 ) * (abs(eta) > 1.95) * ( 0.7 ) +
777 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) +
778 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.8 ) +
779 ( energy >=45 && energy < 80 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) +
780 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) +
781 ( energy >=45 && energy < 80 ) * (abs(eta) <= 0.69) * ( 0.97 ) +
782 ( energy >=80 && energy < 200 ) * (abs(eta) > 1.95) * (0.68 ) +
783 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
784 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) +
785 ( energy >=80 && energy < 200 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) +
786 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
787 ( energy >=80 && energy < 200 ) * (abs(eta) <= 0.69) * (0.98 ) +
788 ( energy >=200 && energy < 400 ) * (abs(eta) > 1.95) * ( 0.68 ) +
789 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) +
790 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) +
791 ( energy >=200 && energy < 400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) +
792 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) +
793 ( energy >=200 && energy < 400 ) * (abs(eta) <= 0.69) * ( 0.98 ) +
794 ( energy >=400 ) * (abs(eta) > 1.95) * (0.68 ) +
795 ( energy >=400 ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
796 ( energy >=400 ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) +
797 ( energy >=400 ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) +
798 ( energy >=400 ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
799 ( energy >=400 ) * (abs(eta) <= 0.69) * (0.98 )
800 }
801}
802
803####################
804# Electron isolation
805####################
806
807module Isolation ElectronIsolation {
808 set CandidateInputArray ElectronEfficiency/electrons
809 set IsolationInputArray EFlowMerger/eflow
810
811 set OutputArray electrons
812
[505a779]813 set DeltaRMax 0.1
[7a1e27d]814
815 set PTMin 0.5
816
[505a779]817 set PTRatioMax 0.2
[7a1e27d]818}
819
820#################
821# Muon efficiency
822#################
823
824module Efficiency MuonEfficiency {
825 set InputArray MuonMomentumSmearing/muons
826 set OutputArray muons
827
828 # set EfficiencyFormula {efficiency as a function of eta and pt}
829
830 # efficiency formula for muons
831 # current full simulation of CLICdet yields:
832
833
834
835 set EfficiencyFormula {
[8c09a6b9]836 (energy < 2.0 ) * (0.00) +
837 (energy>=2.0) * (0.999)
[7a1e27d]838
839 }
840}
841
842################
843# Muon isolation
844################
845
846module Isolation MuonIsolation {
847 set CandidateInputArray MuonEfficiency/muons
848 set IsolationInputArray EFlowMerger/eflow
849
850 set OutputArray muons
851
[505a779]852 set DeltaRMax 0.1
[7a1e27d]853
854 set PTMin 0.5
855
[505a779]856 set PTRatioMax 0.2
[7a1e27d]857}
858
859
860
861###################
862# Missing ET merger
863###################
864
865module Merger MissingET {
866 # add InputArray InputArray
867 add InputArray EFlowMerger/eflow
868 set MomentumOutputArray momentum
869}
870
871
872##################
873# Scalar HT merger
874##################
875
876module Merger ScalarHT {
877 # add InputArray InputArray
878 add InputArray EFlowMerger/eflow
879 set EnergyOutputArray energy
880}
881######################
882# EFlowFilter (UniqueObjectFinder)
883######################
884module UniqueObjectFinder EFlowFilter {
885 add InputArray PhotonIsolation/photons photons
886 add InputArray ElectronIsolation/electrons electrons
887 add InputArray MuonIsolation/muons muons
888 add InputArray EFlowMerger/eflow eflow
889}
890
891#################
892# Neutrino Filter
893#################
894
895module PdgCodeFilter NeutrinoFilter {
896
897 set InputArray Delphes/stableParticles
898 set OutputArray filteredParticles
899
900 set PTMin 0.0
901
902 add PdgCode {12}
903 add PdgCode {14}
904 add PdgCode {16}
905 add PdgCode {-12}
906 add PdgCode {-14}
907 add PdgCode {-16}
908
909}
910
911
912#####################
913# MC truth jet finder
914#####################
915
916module FastJetFinder GenJetFinder {
917 set InputArray NeutrinoFilter/filteredParticles
918
919 set OutputArray jets
920
921 # 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
922 set JetAlgorithm 9
923 set ParameterR 0.5
924
925 set JetPTMin 20.0
926}
927
928#########################
929# Gen Missing ET merger
930########################
931
932module Merger GenMissingET {
933 # add InputArray InputArray
934 add InputArray NeutrinoFilter/filteredParticles
935 set MomentumOutputArray momentum
936}
937
938
939
940############
941# Jet finder
942############
943
944module FastJetFinder FastJetFinderKt {
945 # set InputArray Calorimeter/towers
[eee976c2]946 set InputArray EFlowFilter/eflow
[7a1e27d]947
948 set OutputArray KTjets
949
950 # 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
951 set JetAlgorithm 4
952 set ParameterR 0.5
953
954 set JetPTMin 20.0
955}
956
957################
958# Jet finder VLC
959################
960
961
[96fd15f]962source CLIC/CLICdet_JetReco.tcl
[7a1e27d]963
964
965########################
966# Jet Flavor Association
967########################
968
[96fd15f]969source CLIC/CLICdet_JetFlavorAssociation.tcl
[7a1e27d]970
971###########
972# b-tagging
973###########
974# based on CLICdp-Note-2014-002
975
[96fd15f]976source CLIC/CLICdet_BTagging.tcl
[7a1e27d]977
978
979#############
980# tau-tagging
981#############
982# based on LCD-2010-009
983
[96fd15f]984source CLIC/CLICdet_TauTagging.tcl
[7a1e27d]985
986
987
988##################
989# ROOT tree writer
990##################
991
992module TreeWriter TreeWriter {
993 # add Branch InputArray BranchName BranchClass
994 add Branch Delphes/allParticles Particle GenParticle
995
996 add Branch GenJetFinder/jets GenJet Jet
997
998 add Branch FastJetFinderKt/KTjets KTjet Jet
999 add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet
1000 add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet
1001 add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet
1002 add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet
1003 add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet
1004
1005 add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet
1006 add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet
1007 add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet
1008 add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet
1009 add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet
1010
1011 add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet
1012 add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet
1013 add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet
1014 add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet
1015 add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet
1016
1017 add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet
1018 add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet
1019 add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet
1020 add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet
1021 add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet
1022
1023 add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet
1024 add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet
1025 add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet
1026 add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet
1027 add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
1028
1029 add Branch FastJetFinderVLC_R05_inclusive/VLCjetsR05_inclusive VLCjetR05_inclusive Jet
1030 add Branch FastJetFinderVLC_R07_inclusive/VLCjetsR07_inclusive VLCjetR07_inclusive Jet
1031 add Branch FastJetFinderVLC_R10_inclusive/VLCjetsR10_inclusive VLCjetR10_inclusive Jet
1032 add Branch FastJetFinderVLC_R12_inclusive/VLCjetsR12_inclusive VLCjetR12_inclusive Jet
1033 add Branch FastJetFinderVLC_R15_inclusive/VLCjetsR15_inclusive VLCjetR15_inclusive Jet
1034
1035
1036
1037 ####
1038
1039 add Branch GenMissingET/momentum GenMissingET MissingET
1040
1041 add Branch TrackMerger/tracks Track Track
1042 add Branch Calorimeter/towers Tower Tower
1043
1044 add Branch HCal/eflowTracks EFlowTrack Track
1045 add Branch ECal/eflowPhotons EFlowPhoton Tower
1046 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1047
1048 add Branch EFlowFilter/photons Photon Photon
1049 add Branch EFlowFilter/electrons Electron Electron
1050 add Branch EFlowFilter/muons Muon Muon
1051
1052 add Branch MissingET/momentum MissingET MissingET
1053 add Branch ScalarHT/energy ScalarHT ScalarHT
1054}
1055
Note: See TracBrowser for help on using the repository browser.