Fork me on GitHub

source: git/cards/delphes_card_CLICdet_Stage3.tcl@ d180044

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

fixed muon efficiencies

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