Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage1.tcl@ 28077b2

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 28077b2 was 96fd15f, checked in by Ulrike Schnoor <ulrike.schnoor@…>, 6 years ago

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