Fork me on GitHub

source: git/cards/delphes_card_CLICdet.tcl@ 5a69ba1d

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

fixed input to jet clustering from EFlowFilter
  • Property mode set to 100644
File size: 31.2 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
78
79 MissingET
80 GenMissingET
81
82 JetEnergyScale
83
84 JetFlavorAssociation
85
86 BTagging
87
88 TauTagging
89
90 ScalarHT
91
92
93
94 TreeWriter
95}
96
97#################################
98# Propagate particles in cylinder
99#################################
100
101module ParticlePropagator ParticlePropagator {
102 set InputArray Delphes/stableParticles
103
104 set OutputArray stableParticles
105 set ChargedHadronOutputArray chargedHadrons
106 set ElectronOutputArray electrons
107 set MuonOutputArray muons
108
109 # radius of the magnetic field coverage in the calorimeter, in m
110 set Radius 1.5
111 # half-length of the magnetic field coverage in the calorimeter, in m
112 set HalfLength 2.31
113
114 # magnetic field, in T
115 set Bz 4.0
116}
117
118####################################
119# Charged hadron tracking efficiency
120####################################
121
122module Efficiency ChargedHadronTrackingEfficiency {
123 set InputArray ParticlePropagator/chargedHadrons
124 set OutputArray chargedHadrons
125
126 # FIXME currently uses tracking efficiency from muon simulation
127 set EfficiencyFormula {
128 (abs(eta) > 2.66) * (0.000) +
129 (abs(eta) < 2.66 && abs(eta) > 2.44) * (pt > 0.1) * (0.997) +
130 (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.1 && pt < 0.174) * (0.997) +
131 (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.174) * (1.000) +
132 (abs(eta) < 2.25) * (pt > 0.1) * (1.000) }
133}
134
135##############################
136# Electron tracking efficiency
137##############################
138
139module Efficiency ElectronTrackingEfficiency {
140 set InputArray ParticlePropagator/electrons
141 set OutputArray electrons
142
143 # FIXME currently uses tracking efficiency from muon simulation
144 set EfficiencyFormula {
145 (abs(eta) > 2.66) * (0.000) +
146 (abs(eta) < 2.66 && abs(eta) > 2.44) * (pt > 0.1) * (0.997) +
147 (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.1 && pt < 0.174) * (0.997) +
148 (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.174) * (1.000) +
149 (abs(eta) < 2.25) * (pt > 0.1) * (1.000)
150 }
151}
152
153##########################
154# Muon tracking efficiency
155##########################
156
157module Efficiency MuonTrackingEfficiency {
158 set InputArray ParticlePropagator/muons
159 set OutputArray muons
160
161 # Current full simulation with CLICdet provides for muons:
162 set EfficiencyFormula {
163 (abs(eta) > 2.66) * (0.000) +
164 (abs(eta) < 2.66 && abs(eta) > 2.44) * (pt > 0.1) * (0.997) +
165 (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.1 && pt < 0.174) * (0.997) +
166 (abs(eta) < 2.44 && abs(eta) > 2.25) * (pt > 0.174) * (1.000) +
167 (abs(eta) < 2.25) * (pt > 0.1) * (1.000)
168 }
169}
170
171########################################
172# Momentum resolution for charged tracks
173########################################
174
175module MomentumSmearing ChargedHadronMomentumSmearing {
176 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
177 set OutputArray chargedHadrons
178
179 # Using eta mid-points between evaluated resolutions from full simulation
180 # Resolution given in dpT/pT.
181 # FIXME: currently uses momentum resolution from single muon simulation!
182 set ResolutionFormula {
183 (abs(eta) <= 0.26) * (pt > 0.1) * sqrt(0.0162^2 + pt^2*5.863e-4^2) +
184 (abs(eta) > 0.26 && abs(eta) <= 0.45) * (pt > 0.1) * sqrt(0.0065^2 + pt^2*5.949e-5^2) +
185 (abs(eta) > 0.45 && abs(eta) <= 0.65) * (pt > 0.1) * sqrt(0.0041^2 + pt^2*3.014e-5^2) +
186 (abs(eta) > 0.65 && abs(eta) <= 2.75) * (pt > 0.1) * sqrt(0.0036^2 + pt^2*2.977e-5^2) +
187 (abs(eta) > 2.75) * (pt > 0.1) * sqrt(0.0021^2 + pt^2*2.189e-5^2)
188 }
189}
190
191###################################
192# Momentum resolution for electrons
193###################################
194
195module MomentumSmearing ElectronMomentumSmearing {
196 set InputArray ElectronTrackingEfficiency/electrons
197 set OutputArray electrons
198
199 # Using eta mid-points between evaluated resolutions from full simulation
200 # Resolution given in dpT/pT.
201 # FIXME: currently uses momentum resolution from single muon simulation!
202 set ResolutionFormula {
203 (abs(eta) <= 0.26) * (pt > 0.1) * sqrt(0.0162^2 + pt^2*5.863e-4^2) +
204 (abs(eta) > 0.26 && abs(eta) <= 0.45) * (pt > 0.1) * sqrt(0.0065^2 + pt^2*5.949e-5^2) +
205 (abs(eta) > 0.45 && abs(eta) <= 0.65) * (pt > 0.1) * sqrt(0.0041^2 + pt^2*3.014e-5^2) +
206 (abs(eta) > 0.65 && abs(eta) <= 2.75) * (pt > 0.1) * sqrt(0.0036^2 + pt^2*2.977e-5^2) +
207 (abs(eta) > 2.75) * (pt > 0.1) * sqrt(0.0021^2 + pt^2*2.189e-5^2)
208 }
209}
210
211###############################
212# Momentum resolution for muons
213###############################
214
215module MomentumSmearing MuonMomentumSmearing {
216 set InputArray MuonTrackingEfficiency/muons
217 set OutputArray muons
218
219 # Using eta mid-points between evaluated resolutions from full simulation
220 # Resolution given in dpT/pT.
221 set ResolutionFormula {
222 (abs(eta) <= 0.26) * (pt > 0.1) * sqrt(0.0162^2 + pt^2*5.863e-4^2) +
223 (abs(eta) > 0.26 && abs(eta) <= 0.45) * (pt > 0.1) * sqrt(0.0065^2 + pt^2*5.949e-5^2) +
224 (abs(eta) > 0.45 && abs(eta) <= 0.65) * (pt > 0.1) * sqrt(0.0041^2 + pt^2*3.014e-5^2) +
225 (abs(eta) > 0.65 && abs(eta) <= 2.75) * (pt > 0.1) * sqrt(0.0036^2 + pt^2*2.977e-5^2) +
226 (abs(eta) > 2.75) * (pt > 0.1) * sqrt(0.0021^2 + pt^2*2.189e-5^2)
227 }
228}
229
230##############
231# Track merger
232##############
233
234module Merger TrackMerger {
235 # add InputArray InputArray
236 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
237 add InputArray ElectronMomentumSmearing/electrons
238 add InputArray MuonMomentumSmearing/muons
239 set OutputArray tracks
240}
241
242#############
243# ECAL
244#############
245
246module SimpleCalorimeter ECal {
247 set ParticleInputArray ParticlePropagator/stableParticles
248 set TrackInputArray TrackMerger/tracks
249
250 set TowerOutputArray ecalTowers
251 set EFlowTrackOutputArray eflowTracks
252 set EFlowTowerOutputArray eflowPhotons
253
254 set IsEcal true
255
256 set EnergyMin 0.5
257 set EnergySignificanceMin 1.0
258
259 set SmearTowerCenter true
260
261 set pi [expr {acos(-1)}]
262
263 # lists of the edges of each tower in eta and phi
264 # each list starts with the lower edge of the first tower
265 # the list ends with the higher edged of the last tower
266
267 #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
268 #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
269 #ECAL plug: dphi = 1 degree, deta = 0.02 up to |eta| <=3
270 #ECAL cell sizes always 5x5 mm^2
271
272 #barrel:
273 #dphi = 0.2 degree towers up to eta <=1.2
274 set PhiBins {}
275 for {set i -900} {$i <= 900} {incr i} {
276 add PhiBins [expr {$i * $pi/900.0 }]
277 }
278 # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
279 for {set i -400} {$i <=400} {incr i} {
280 set eta [expr {$i * 0.003}]
281 add EtaPhiBins $eta $PhiBins
282 }
283
284 #endcaps:
285 #dphi = 0.8 degree towers for 1.2 < eta <=2.5
286 set PhiBins {}
287 for {set i -225} {$i <= 225} {incr i} {
288 add PhiBins [expr {$i * $pi/225.}]
289 }
290 #deta=0.02 units for 1.2 < |eta| <=2.5
291 #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
292 for {set i 1} {$i <=66} {incr i} {
293 set eta [expr {-2.52 + $i * 0.02}]
294 add EtaPhiBins $eta $PhiBins
295 }
296 #same for 1.2 to 2.5
297 for {set i 1} {$i <=66} {incr i} {
298 set eta [expr {1.18 + $i*0.02}]
299 add EtaPhiBins $eta $PhiBins
300 }
301
302 #plug:
303 #dphi = 1 degree for 2.5 < eta <=3
304 set PhiBins {}
305 for {set i -180} {$i <= 180} {incr i} {
306 add PhiBins [expr {$i * $pi/180.}]
307 }
308 # deta = 0.02 for 2.5 < |eta| <=3
309 # from -3 to -2.5, there will be 25 segments
310 for {set i 1} {$i <= 26} {incr i} {
311 set eta [expr {-3.02 + $i * 0.02}]
312 add EtaPhiBins $eta $PhiBins
313 }
314 #same for 2.5 to 3
315 for {set i 1} {$i <= 26} {incr i} {
316 set eta [expr {2.48 + $i*0.02}]
317 add EtaPhiBins $eta $PhiBins
318 }
319
320
321
322 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
323
324 add EnergyFraction {0} {0.0}
325 # energy fractions for e, gamma and pi0
326 add EnergyFraction {11} {1.0}
327 add EnergyFraction {22} {1.0}
328 add EnergyFraction {111} {1.0}
329 # energy fractions for muon, neutrinos and neutralinos
330 add EnergyFraction {12} {0.0}
331 add EnergyFraction {13} {0.0}
332 add EnergyFraction {14} {0.0}
333 add EnergyFraction {16} {0.0}
334 add EnergyFraction {1000022} {0.0}
335 add EnergyFraction {1000023} {0.0}
336 add EnergyFraction {1000025} {0.0}
337 add EnergyFraction {1000035} {0.0}
338 add EnergyFraction {1000045} {0.0}
339 # energy fractions for K0short and Lambda
340 add EnergyFraction {310} {0.3}
341 add EnergyFraction {3122} {0.3}
342
343 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
344
345 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.15^2) }
346
347}
348
349#############
350# HCAL
351#############
352
353module SimpleCalorimeter HCal {
354 set ParticleInputArray ParticlePropagator/stableParticles
355 set TrackInputArray ECal/eflowTracks
356
357 set TowerOutputArray hcalTowers
358 set EFlowTrackOutputArray eflowTracks
359 set EFlowTowerOutputArray eflowNeutralHadrons
360
361 set IsEcal false
362
363 set EnergyMin 1.0
364 set EnergySignificanceMin 1.0
365
366 set SmearTowerCenter true
367
368 set pi [expr {acos(-1)}]
369
370 # lists of the edges of each tower in eta and phi
371 # each list starts with the lower edge of the first tower
372 # the list ends with the higher edged of the last tower
373
374
375 #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
376 #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
377 #HCAL endcaps: dphi = 6 degree, deta = 0.1 up to |eta| <=3.5
378 #HCAL cell sizes always 30x30 mm^2
379
380 #barrel and ring:
381 #dphi = 1 degree up to |eta| <=0.9
382 set PhiBins {}
383 for {set i -180} {$i <=180} {incr i} {
384 add PhiBins [expr {$i * $pi/180.0}]
385 }
386 #deta= 0.02 towers up to |eta| <=0.9
387 for {set i -45} {$i <=45} {incr i} {
388 set eta [expr {$i * 0.02}]
389 add EtaPhiBins $eta $PhiBins
390 }
391
392 #endcaps:
393 # dphi = 6 degree
394 set PhiBins {}
395 for {set i -30} {$i <=30} {incr i} {
396 add PhiBins [expr {$i * $pi/30.0}]
397 }
398 # deta =0.1 for 0.9 < |eta| <=3.5
399 #for -3.5 to -0.9, 26 segments
400 for {set i 1} {$i <=27} {incr i} {
401 set eta [expr {-3.6 + $i * 0.1}]
402 add EtaPhiBins $eta $PhiBins
403 }
404 #same for 0.9 to 3.5
405 for {set i 1} {$i <=27} {incr i} {
406 set eta [expr {0.8 + $i * 0.1 }]
407 add EtaPhiBins $eta $PhiBins
408 }
409
410 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
411 add EnergyFraction {0} {1.0}
412 # energy fractions for e, gamma and pi0
413 add EnergyFraction {11} {0.0}
414 add EnergyFraction {22} {0.0}
415 add EnergyFraction {111} {0.0}
416 # energy fractions for muon, neutrinos and neutralinos
417 add EnergyFraction {12} {0.0}
418 add EnergyFraction {13} {0.0}
419 add EnergyFraction {14} {0.0}
420 add EnergyFraction {16} {0.0}
421 add EnergyFraction {1000022} {0.0}
422 add EnergyFraction {1000023} {0.0}
423 add EnergyFraction {1000025} {0.0}
424 add EnergyFraction {1000035} {0.0}
425 add EnergyFraction {1000045} {0.0}
426 # energy fractions for K0short and Lambda
427 add EnergyFraction {310} {0.7}
428 add EnergyFraction {3122} {0.7}
429
430 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
431
432 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.015^2 + energy*0.50^2)}
433
434}
435
436#################
437# Electron filter
438#################
439
440module PdgCodeFilter ElectronFilter {
441 set InputArray HCal/eflowTracks
442 set OutputArray electrons
443 set Invert true
444 add PdgCode {11}
445 add PdgCode {-11}
446}
447
448######################
449# ChargedHadronFilter
450######################
451
452module PdgCodeFilter ChargedHadronFilter {
453 set InputArray HCal/eflowTracks
454 set OutputArray chargedHadrons
455
456 add PdgCode {11}
457 add PdgCode {-11}
458 add PdgCode {13}
459 add PdgCode {-13}
460}
461
462
463
464###################################################
465# Tower Merger (in case not using e-flow algorithm)
466###################################################
467
468module Merger Calorimeter {
469 # add InputArray InputArray
470 add InputArray ECal/ecalTowers
471 add InputArray HCal/hcalTowers
472 set OutputArray towers
473}
474
475
476####################
477# Energy flow merger
478####################
479
480module Merger EFlowMerger {
481 # add InputArray InputArray
482 add InputArray HCal/eflowTracks
483 add InputArray ECal/eflowPhotons
484 add InputArray HCal/eflowNeutralHadrons
485 set OutputArray eflow
486}
487
488##here
489
490###################
491# Photon efficiency
492###################
493
494module Efficiency PhotonEfficiency {
495 set InputArray ECal/eflowPhotons
496 set OutputArray photons
497
498 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
499
500 # efficiency formula for photons
501 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
502 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
503 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
504 (abs(eta) > 2.5) * (0.00)}
505}
506
507##################
508# Photon isolation
509##################
510
511module Isolation PhotonIsolation {
512 set CandidateInputArray PhotonEfficiency/photons
513 set IsolationInputArray EFlowMerger/eflow
514
515 set OutputArray photons
516
517 set DeltaRMax 0.5
518
519 set PTMin 0.5
520
521 set PTRatioMax 0.12
522}
523
524#####################
525# Electron efficiency
526#####################
527
528module Efficiency ElectronEfficiency {
529 set InputArray ElectronFilter/electrons
530 set OutputArray electrons
531
532 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
533
534 # efficiency formula for electrons
535 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
536 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
537 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
538 (abs(eta) > 2.5) * (0.00)}
539}
540
541####################
542# Electron isolation
543####################
544
545module Isolation ElectronIsolation {
546 set CandidateInputArray ElectronEfficiency/electrons
547 set IsolationInputArray EFlowMerger/eflow
548
549 set OutputArray electrons
550
551 set DeltaRMax 0.5
552
553 set PTMin 0.5
554
555 set PTRatioMax 0.12
556}
557
558#################
559# Muon efficiency
560#################
561
562module Efficiency MuonEfficiency {
563 set InputArray MuonMomentumSmearing/muons
564 set OutputArray muons
565
566 # set EfficiencyFormula {efficiency as a function of eta and pt}
567
568 # efficiency formula for muons
569 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
570 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
571 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
572 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
573 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
574 (abs(eta) > 2.4) * (0.00)}
575}
576
577################
578# Muon isolation
579################
580
581module Isolation MuonIsolation {
582 set CandidateInputArray MuonEfficiency/muons
583 set IsolationInputArray EFlowMerger/eflow
584
585 set OutputArray muons
586
587 set DeltaRMax 0.5
588
589 set PTMin 0.5
590
591 set PTRatioMax 0.25
592}
593
594
595
596###################
597# Missing ET merger
598###################
599
600module Merger MissingET {
601 # add InputArray InputArray
602 add InputArray EFlowMerger/eflow
603 set MomentumOutputArray momentum
604}
605
606
607##################
608# Scalar HT merger
609##################
610
611module Merger ScalarHT {
612 # add InputArray InputArray
613 add InputArray EFlowMerger/eflow
614 set EnergyOutputArray energy
615}
616######################
617# EFlowFilter (UniqueObjectFinder)
618######################
619module UniqueObjectFinder EFlowFilter {
620 add InputArray PhotonIsolation/photons photons
621 add InputArray ElectronIsolation/electrons electrons
622 add InputArray MuonIsolation/muons muons
623 add InputArray EFlowMerger/eflow eflow
624}
625
626
627#
628#module PdgCodeFilter EFlowFilter {
629# set InputArray EFlowMerger/eflow
630# set OutputArray eflow
631#
632# add PdgCode {11}
633# add PdgCode {-11}
634# add PdgCode {13}
635# add PdgCode {-13}
636#}
637
638#################
639# Neutrino Filter
640#################
641
642module PdgCodeFilter NeutrinoFilter {
643
644 set InputArray Delphes/stableParticles
645 set OutputArray filteredParticles
646
647 set PTMin 0.0
648
649 add PdgCode {12}
650 add PdgCode {14}
651 add PdgCode {16}
652 add PdgCode {-12}
653 add PdgCode {-14}
654 add PdgCode {-16}
655
656}
657
658
659#####################
660# MC truth jet finder
661#####################
662
663module FastJetFinder GenJetFinder {
664 set InputArray NeutrinoFilter/filteredParticles
665
666 set OutputArray jets
667
668 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
669 set JetAlgorithm 9
670 set ParameterR 0.5
671
672 set JetPTMin 20.0
673}
674
675#########################
676# Gen Missing ET merger
677########################
678
679module Merger GenMissingET {
680 # add InputArray InputArray
681 add InputArray NeutrinoFilter/filteredParticles
682 set MomentumOutputArray momentum
683}
684
685
686
687############
688# Jet finder
689############
690
691module FastJetFinder FastJetFinderKt {
692 # set InputArray Calorimeter/towers
693 set InputArray EFlowMerger/eflow
694
695 set OutputArray KTjets
696
697 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
698 set JetAlgorithm 4
699 set ParameterR 0.5
700
701 set JetPTMin 20.0
702}
703
704############
705# Jet finder VLC
706############
707#R05 N2
708module FastJetFinder FastJetFinderVLC_R05_N2 {
709 # set InputArray Calorimeter/towers
710 set InputArray EFlowFilter/eflow
711
712 set OutputArray VLCjetsR05N2
713
714 # 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
715 set NJets 2
716 set ExclusiveClustering true
717 set JetAlgorithm 9
718 set ParameterR 0.5
719 set Beta 1.0
720 set Gamma 1.0
721
722 set JetPTMin 20.0
723}
724#R05 N3
725module FastJetFinder FastJetFinderVLC_R05_N3 {
726 set InputArray EFlowFilter/eflow
727 set OutputArray VLCjetsR05N3
728 set NJets 3
729 set ExclusiveClustering true
730 set JetAlgorithm 9
731 set ParameterR 0.5
732 set Beta 1.0
733 set Gamma 1.0
734 set JetPTMin 20.0
735}
736#R05 N4
737module FastJetFinder FastJetFinderVLC_R05_N4 {
738 set InputArray EFlowFilter/eflow
739 set OutputArray VLCjetsR05N4
740 set NJets 4
741 set ExclusiveClustering true
742 set JetAlgorithm 9
743 set ParameterR 0.5
744 set Beta 1.0
745 set Gamma 1.0
746 set JetPTMin 20.0
747}
748#R05 N5
749module FastJetFinder FastJetFinderVLC_R05_N5 {
750 set InputArray EFlowFilter/eflow
751 set OutputArray VLCjetsR05N5
752 set NJets 5
753 set ExclusiveClustering true
754 set JetAlgorithm 9
755 set ParameterR 0.5
756 set Beta 1.0
757 set Gamma 1.0
758 set JetPTMin 20.0
759}
760#R05 N6
761module FastJetFinder FastJetFinderVLC_R05_N6 {
762 set InputArray EFlowFilter/eflow
763 set OutputArray VLCjetsR05N6
764 set NJets 6
765 set ExclusiveClustering true
766 set JetAlgorithm 9
767 set ParameterR 0.5
768 set Beta 1.0
769 set Gamma 1.0
770 set JetPTMin 20.0
771}
772#R07 N2
773module FastJetFinder FastJetFinderVLC_R07_N2 {
774 set InputArray EFlowFilter/eflow
775 set OutputArray VLCjetsR07N2
776 set NJets 2
777 set ExclusiveClustering true
778 set JetAlgorithm 9
779 set ParameterR 0.7
780 set Beta 1.0
781 set Gamma 1.0
782 set JetPTMin 20.0
783}
784#R07 N3
785module FastJetFinder FastJetFinderVLC_R07_N3 {
786 set InputArray EFlowFilter/eflow
787 set OutputArray VLCjetsR07N3
788 set NJets 3
789 set ExclusiveClustering true
790 set JetAlgorithm 9
791 set ParameterR 0.7
792 set Beta 1.0
793 set Gamma 1.0
794 set JetPTMin 20.0
795}
796#R07 N4
797module FastJetFinder FastJetFinderVLC_R07_N4 {
798 set InputArray EFlowFilter/eflow
799 set OutputArray VLCjetsR07N4
800 set NJets 4
801 set ExclusiveClustering true
802 set JetAlgorithm 9
803 set ParameterR 0.7
804 set Beta 1.0
805 set Gamma 1.0
806 set JetPTMin 20.0
807}
808#R07 N5
809module FastJetFinder FastJetFinderVLC_R07_N5 {
810 set InputArray EFlowFilter/eflow
811 set OutputArray VLCjetsR07N5
812 set NJets 5
813 set ExclusiveClustering true
814 set JetAlgorithm 9
815 set ParameterR 0.7
816 set Beta 1.0
817 set Gamma 1.0
818 set JetPTMin 20.0
819}
820#R07 N6
821module FastJetFinder FastJetFinderVLC_R07_N6 {
822 set InputArray EFlowFilter/eflow
823 set OutputArray VLCjetsR07N6
824 set NJets 6
825 set ExclusiveClustering true
826 set JetAlgorithm 9
827 set ParameterR 0.7
828 set Beta 1.0
829 set Gamma 1.0
830 set JetPTMin 20.0
831}
832
833#R10N2
834module FastJetFinder FastJetFinderVLC_R10_N2 {
835 set InputArray EFlowFilter/eflow
836 set OutputArray VLCjetsR10N2
837 set NJets 2
838 set ExclusiveClustering true
839 set JetAlgorithm 9
840 set ParameterR 1.0
841 set Beta 1.0
842 set Gamma 1.0
843 set JetPTMin 20.0
844}
845#R10 N3
846module FastJetFinder FastJetFinderVLC_R10_N3 {
847 set InputArray EFlowFilter/eflow
848 set OutputArray VLCjetsR10N3
849 set NJets 3
850 set ExclusiveClustering true
851 set JetAlgorithm 9
852 set ParameterR 1.0
853 set Beta 1.0
854 set Gamma 1.0
855 set JetPTMin 20.0
856}
857#R10 N4
858module FastJetFinder FastJetFinderVLC_R10_N4 {
859 set InputArray EFlowFilter/eflow
860 set OutputArray VLCjetsR10N4
861 set NJets 4
862 set ExclusiveClustering true
863 set JetAlgorithm 9
864 set ParameterR 1.0
865 set Beta 1.0
866 set Gamma 1.0
867 set JetPTMin 20.0
868}
869#R10 N5
870module FastJetFinder FastJetFinderVLC_R10_N5 {
871 set InputArray EFlowFilter/eflow
872 set OutputArray VLCjetsR10N5
873 set NJets 5
874 set ExclusiveClustering true
875 set JetAlgorithm 9
876 set ParameterR 1.0
877 set Beta 1.0
878 set Gamma 1.0
879 set JetPTMin 20.0
880}
881#R10 N6
882module FastJetFinder FastJetFinderVLC_R10_N6 {
883 set InputArray EFlowFilter/eflow
884 set OutputArray VLCjetsR10N6
885 set NJets 6
886 set ExclusiveClustering true
887 set JetAlgorithm 9
888 set ParameterR 1.0
889 set Beta 1.0
890 set Gamma 1.0
891 set JetPTMin 20.0
892}
893
894#R12 N2
895module FastJetFinder FastJetFinderVLC_R12_N2 {
896 set InputArray EFlowFilter/eflow
897 set OutputArray VLCjetsR12N2
898 set NJets 2
899 set ExclusiveClustering true
900 set JetAlgorithm 9
901 set ParameterR 1.2
902 set Beta 1.0
903 set Gamma 1.0
904 set JetPTMin 20.0
905}
906#R12 N3
907module FastJetFinder FastJetFinderVLC_R12_N3 {
908 set InputArray EFlowFilter/eflow
909 set OutputArray VLCjetsR12N3
910 set NJets 3
911 set ExclusiveClustering true
912 set JetAlgorithm 9
913 set ParameterR 1.2
914 set Beta 1.0
915 set Gamma 1.0
916 set JetPTMin 20.0
917}
918#R12 N4
919module FastJetFinder FastJetFinderVLC_R12_N4 {
920 set InputArray EFlowFilter/eflow
921 set OutputArray VLCjetsR12N4
922 set NJets 4
923 set ExclusiveClustering true
924 set JetAlgorithm 9
925 set ParameterR 1.2
926 set Beta 1.0
927 set Gamma 1.0
928 set JetPTMin 20.0
929}
930#R12 N5
931module FastJetFinder FastJetFinderVLC_R12_N5 {
932 set InputArray EFlowFilter/eflow
933 set OutputArray VLCjetsR12N5
934 set NJets 5
935 set ExclusiveClustering true
936 set JetAlgorithm 9
937 set ParameterR 1.2
938 set Beta 1.0
939 set Gamma 1.0
940 set JetPTMin 20.0
941}
942#R12 N6
943module FastJetFinder FastJetFinderVLC_R12_N6 {
944 set InputArray EFlowFilter/eflow
945 set OutputArray VLCjetsR12N6
946 set NJets 6
947 set ExclusiveClustering true
948 set JetAlgorithm 9
949 set ParameterR 1.2
950 set Beta 1.0
951 set Gamma 1.0
952 set JetPTMin 20.0
953}
954
955
956#R15 N2
957module FastJetFinder FastJetFinderVLC_R15_N2 {
958 set InputArray EFlowFilter/eflow
959 set OutputArray VLCjetsR15N2
960 set NJets 2
961 set ExclusiveClustering true
962 set JetAlgorithm 9
963 set ParameterR 1.5
964 set Beta 1.0
965 set Gamma 1.0
966 set JetPTMin 20.0
967}
968#R15 N3
969module FastJetFinder FastJetFinderVLC_R15_N3 {
970 set InputArray EFlowFilter/eflow
971 set OutputArray VLCjetsR15N3
972 set NJets 3
973 set ExclusiveClustering true
974 set JetAlgorithm 9
975 set ParameterR 1.5
976 set Beta 1.0
977 set Gamma 1.0
978 set JetPTMin 20.0
979}
980#R15 N4
981module FastJetFinder FastJetFinderVLC_R15_N4 {
982 set InputArray EFlowFilter/eflow
983 set OutputArray VLCjetsR15N4
984 set NJets 4
985 set ExclusiveClustering true
986 set JetAlgorithm 9
987 set ParameterR 1.5
988 set Beta 1.0
989 set Gamma 1.0
990 set JetPTMin 20.0
991}
992#R15 N5
993module FastJetFinder FastJetFinderVLC_R15_N5 {
994 set InputArray EFlowFilter/eflow
995 set OutputArray VLCjetsR15N5
996 set NJets 5
997 set ExclusiveClustering true
998 set JetAlgorithm 9
999 set ParameterR 1.5
1000 set Beta 1.0
1001 set Gamma 1.0
1002 set JetPTMin 20.0
1003}
1004#R15 N6
1005module FastJetFinder FastJetFinderVLC_R15_N6 {
1006 set InputArray EFlowFilter/eflow
1007 set OutputArray VLCjetsR15N6
1008 set NJets 6
1009 set ExclusiveClustering true
1010 set JetAlgorithm 9
1011 set ParameterR 1.5
1012 set Beta 1.0
1013 set Gamma 1.0
1014 set JetPTMin 20.0
1015}
1016##################
1017# Jet Energy Scale
1018##################
1019
1020module EnergyScale JetEnergyScale {
1021 set InputArray FastJetFinderVLC_R10_N4/VLCjetsR10N4
1022 set OutputArray jets
1023
1024 # scale formula for jets
1025 set ScaleFormula {1.00}
1026}
1027
1028
1029########################
1030# Jet Flavor Association
1031########################
1032
1033module JetFlavorAssociation JetFlavorAssociation {
1034
1035 set PartonInputArray Delphes/partons
1036 set ParticleInputArray Delphes/allParticles
1037 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1038 set JetInputArray JetEnergyScale/jets
1039
1040 set DeltaR 0.5
1041 set PartonPTMin 1.0
1042 set PartonEtaMax 2.5
1043
1044}
1045
1046
1047###########
1048# b-tagging
1049###########
1050
1051module BTagging BTagging {
1052 set JetInputArray JetEnergyScale/jets
1053
1054 set BitNumber 0
1055
1056 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1057 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
1058 # gluon's PDG code has the lowest priority
1059
1060 # based on arXiv:1211.4462
1061
1062 # default efficiency formula (misidentification rate)
1063 add EfficiencyFormula {0} {0.01+0.000038*pt}
1064
1065 # efficiency formula for c-jets (misidentification rate)
1066 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
1067
1068 # efficiency formula for b-jets
1069 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
1070}
1071
1072#############
1073# tau-tagging
1074#############
1075
1076
1077module TauTagging TauTagging {
1078 set ParticleInputArray Delphes/allParticles
1079 set PartonInputArray Delphes/partons
1080 set JetInputArray JetEnergyScale/jets
1081
1082 set DeltaR 0.5
1083
1084 set TauPTMin 1.0
1085
1086 set TauEtaMax 4.0
1087
1088 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1089
1090 # default efficiency formula (misidentification rate)
1091 add EfficiencyFormula {0} {0.001}
1092 # efficiency formula for tau-jets
1093 add EfficiencyFormula {15} {0.4}
1094}
1095
1096#####################################################
1097# Find uniquely identified photons/electrons/tau/jets
1098#####################################################
1099
1100#module UniqueObjectFinder UniqueObjectFinder {
1101# # earlier arrays take precedence over later ones
1102# # add InputArray InputArray OutputArray
1103# add InputArray PhotonIsolation/photons photons
1104# add InputArray ElectronIsolation/electrons electrons
1105# add InputArray MuonIsolation/muons muons
1106# add InputArray JetEnergyScale/jets jets
1107####not necessary any more, since these objects have already been removed if necessary??
1108#}
1109#
1110
1111##################
1112# ROOT tree writer
1113##################
1114
1115module TreeWriter TreeWriter {
1116 # add Branch InputArray BranchName BranchClass
1117 add Branch Delphes/allParticles Particle GenParticle
1118
1119 add Branch GenJetFinder/jets GenJet Jet
1120
1121 add Branch FastJetFinderKt/KTjets KTjet Jet
1122 add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 VLCjetR05N2 Jet
1123 add Branch FastJetFinderVLC_R05_N3/VLCjetsR05N3 VLCjetR05N3 Jet
1124 add Branch FastJetFinderVLC_R05_N4/VLCjetsR05N4 VLCjetR05N4 Jet
1125 add Branch FastJetFinderVLC_R05_N5/VLCjetsR05N5 VLCjetR05N5 Jet
1126 add Branch FastJetFinderVLC_R05_N6/VLCjetsR05N6 VLCjetR05N6 Jet
1127
1128 add Branch FastJetFinderVLC_R07_N2/VLCjetsR07N2 VLCjetR07N2 Jet
1129 add Branch FastJetFinderVLC_R07_N3/VLCjetsR07N3 VLCjetR07N3 Jet
1130 add Branch FastJetFinderVLC_R07_N4/VLCjetsR07N4 VLCjetR07N4 Jet
1131 add Branch FastJetFinderVLC_R07_N5/VLCjetsR07N5 VLCjetR07N5 Jet
1132 add Branch FastJetFinderVLC_R07_N6/VLCjetsR07N6 VLCjetR07N6 Jet
1133
1134 add Branch FastJetFinderVLC_R10_N2/VLCjetsR10N2 VLCjetR10N2 Jet
1135 add Branch FastJetFinderVLC_R10_N3/VLCjetsR10N3 VLCjetR10N3 Jet
1136 add Branch FastJetFinderVLC_R10_N4/VLCjetsR10N4 VLCjetR10N4 Jet
1137 add Branch FastJetFinderVLC_R10_N5/VLCjetsR10N5 VLCjetR10N5 Jet
1138 add Branch FastJetFinderVLC_R10_N6/VLCjetsR10N6 VLCjetR10N6 Jet
1139
1140 add Branch FastJetFinderVLC_R12_N2/VLCjetsR12N2 VLCjetR12N2 Jet
1141 add Branch FastJetFinderVLC_R12_N3/VLCjetsR12N3 VLCjetR12N3 Jet
1142 add Branch FastJetFinderVLC_R12_N4/VLCjetsR12N4 VLCjetR12N4 Jet
1143 add Branch FastJetFinderVLC_R12_N5/VLCjetsR12N5 VLCjetR12N5 Jet
1144 add Branch FastJetFinderVLC_R12_N6/VLCjetsR12N6 VLCjetR12N6 Jet
1145
1146 add Branch FastJetFinderVLC_R15_N2/VLCjetsR15N2 VLCjetR15N2 Jet
1147 add Branch FastJetFinderVLC_R15_N3/VLCjetsR15N3 VLCjetR15N3 Jet
1148 add Branch FastJetFinderVLC_R15_N4/VLCjetsR15N4 VLCjetR15N4 Jet
1149 add Branch FastJetFinderVLC_R15_N5/VLCjetsR15N5 VLCjetR15N5 Jet
1150 add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
1151
1152
1153 add Branch GenMissingET/momentum GenMissingET MissingET
1154
1155 add Branch TrackMerger/tracks Track Track
1156 add Branch Calorimeter/towers Tower Tower
1157
1158 add Branch HCal/eflowTracks EFlowTrack Track
1159 add Branch ECal/eflowPhotons EFlowPhoton Tower
1160 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1161
1162# add Branch UniqueObjectFinder/photons Photon Photon
1163# add Branch UniqueObjectFinder/electrons Electron Electron
1164# add Branch UniqueObjectFinder/muons Muon Muon
1165# add Branch UniqueObjectFinder/jets Jet Jet
1166#
1167
1168# add Branch PhotonIsolation/photons Photon Photon
1169# add Branch ElectronIsolation/electrons Electron Electron
1170# add Branch MuonIsolation/muons Muon Muon
1171
1172 add Branch EFlowFilter/photons Photon Photon
1173 add Branch EFlowFilter/electrons Electron Electron
1174 add Branch EFlowFilter/muons Muon Muon
1175
1176
1177
1178 add Branch MissingET/momentum MissingET MissingET
1179 add Branch ScalarHT/energy ScalarHT ScalarHT
1180}
1181
Note: See TracBrowser for help on using the repository browser.