Fork me on GitHub

source: git/cards/delphes_card_CLICdet.tcl@ a5f4828

ImprovedOutputFile Timing dual_readout llp
Last change on this file since a5f4828 was a5f4828, checked in by Ulrike Schnoor <schnooru@…>, 7 years ago

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