Fork me on GitHub

source: git/cards/delphes_card_CLICdet_380GeV.tcl@ 7a1e27d

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

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