Fork me on GitHub

source: git/cards/delphes_card_CEPC.tcl@ a0f2226

ImprovedOutputFile Timing dual_readout llp
Last change on this file since a0f2226 was 25ac711, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

first commit of CEPC card

  • Property mode set to 100644
File size: 16.9 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17
18 ECal
19 HCal
20
21 Calorimeter
22 EFlowMerger
23 EFlowFilter
24
25 PhotonEfficiency
26 PhotonIsolation
27
28 ElectronFilter
29 ElectronEfficiency
30 ElectronIsolation
31
32 MuonEfficiency
33 MuonIsolation
34
35 NeutrinoFilter
36 GenJetFinder
37 FastJetFinder
38
39 MissingET
40 GenMissingET
41
42 JetEnergyScale
43
44 JetFlavorAssociation
45
46 BTagging
47 CTagging
48
49 TauTagging
50
51 ScalarHT
52
53 UniqueObjectFinder
54
55 TreeWriter
56}
57
58#################################
59# Propagate particles in cylinder
60#################################
61
62module ParticlePropagator ParticlePropagator {
63 set InputArray Delphes/stableParticles
64
65 set OutputArray stableParticles
66 set ChargedHadronOutputArray chargedHadrons
67 set ElectronOutputArray electrons
68 set MuonOutputArray muons
69
70 # radius of the magnetic field coverage, in m
71 set Radius 1.81
72 # half-length of the magnetic field coverage, in m
73 set HalfLength 2.35
74
75 # magnetic field
76 set Bz 3.5
77}
78
79####################################
80# Charged hadron tracking efficiency
81####################################
82
83module Efficiency ChargedHadronTrackingEfficiency {
84 set InputArray ParticlePropagator/chargedHadrons
85 set OutputArray chargedHadrons
86
87 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
88
89 # tracking efficiency formula for charged hadrons
90 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
91 (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
92 (abs(eta) > 3.0) * (0.00)}
93}
94
95##############################
96# Electron tracking efficiency
97##############################
98
99module Efficiency ElectronTrackingEfficiency {
100 set InputArray ParticlePropagator/electrons
101 set OutputArray electrons
102
103 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
104
105 # tracking efficiency formula for electrons
106 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
107 (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
108 (abs(eta) > 3.0) * (0.00)}
109}
110
111##########################
112# Muon tracking efficiency
113##########################
114
115module Efficiency MuonTrackingEfficiency {
116 set InputArray ParticlePropagator/muons
117 set OutputArray muons
118
119 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
120
121 # tracking efficiency formula for muons
122 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
123 (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
124 (abs(eta) > 3.0) * (0.00)}
125}
126
127########################################
128# Momentum resolution for charged tracks
129########################################
130
131module MomentumSmearing ChargedHadronMomentumSmearing {
132 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
133 set OutputArray chargedHadrons
134
135 # set ResolutionFormula {resolution formula as a function of eta and pt}
136
137 # resolution formula for charged hadrons
138 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(0.001^2 + pt^2*1.e-5^2) +
139 (abs(eta) > 1.0 && abs(eta) <= 3.0) * sqrt(0.01^2 + pt^2*1.e-4^2)}
140
141
142}
143
144###################################
145# Momentum resolution for electrons
146###################################
147
148module MomentumSmearing ElectronMomentumSmearing {
149 set InputArray ElectronTrackingEfficiency/electrons
150 set OutputArray electrons
151
152 # set ResolutionFormula {resolution formula as a function of eta and energy}
153
154 # resolution formula for charged hadrons
155 set ResolutionFormula { (abs(eta) <= 1.0) * sqrt(0.001^2 + pt^2*1.e-5^2) +
156 (abs(eta) > 1.0 && abs(eta) <= 3.0) * sqrt(0.01^2 + pt^2*1.e-4^2)}
157}
158
159###############################
160# Momentum resolution for muons
161###############################
162
163module MomentumSmearing MuonMomentumSmearing {
164 set InputArray MuonTrackingEfficiency/muons
165 set OutputArray muons
166
167 # set ResolutionFormula {resolution formula as a function of eta and pt}
168
169 # resolution formula for charged hadrons
170 set ResolutionFormula { (abs(eta) <= 1.0) * sqrt(0.001^2 + pt^2*1.e-5^2) +
171 (abs(eta) > 1.0 && abs(eta) <= 3.0) * sqrt(0.01^2 + pt^2*1.e-4^2)}
172
173}
174
175##############
176# Track merger
177##############
178
179module Merger TrackMerger {
180# add InputArray InputArray
181 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
182 add InputArray ElectronMomentumSmearing/electrons
183 add InputArray MuonMomentumSmearing/muons
184 set OutputArray tracks
185}
186
187#############
188# ECAL
189#############
190
191module SimpleCalorimeter ECal {
192 set ParticleInputArray ParticlePropagator/stableParticles
193 set TrackInputArray TrackMerger/tracks
194
195 set TowerOutputArray ecalTowers
196 set EFlowTrackOutputArray eflowTracks
197 set EFlowTowerOutputArray eflowPhotons
198
199 set IsEcal true
200
201 set EnergyMin 0.5
202 set EnergySignificanceMin 1.0
203
204 set SmearTowerCenter true
205
206 set pi [expr {acos(-1)}]
207
208 # lists of the edges of each tower in eta and phi
209 # each list starts with the lower edge of the first tower
210 # the list ends with the higher edged of the last tower
211
212 # 1.0 degree towers (3 cm x 3 cm)
213 set PhiBins {}
214 for {set i -180} {$i <= 180} {incr i} {
215 add PhiBins [expr {$i * $pi/180.0}]
216 }
217
218 # 0.02 unit in eta up to eta = 3.0
219 for {set i -150} {$i <= 150} {incr i} {
220 set eta [expr {$i * 0.02}]
221 add EtaPhiBins $eta $PhiBins
222 }
223
224 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
225
226 add EnergyFraction {0} {0.0}
227 # energy fractions for e, gamma and pi0
228 add EnergyFraction {11} {1.0}
229 add EnergyFraction {22} {1.0}
230 add EnergyFraction {111} {1.0}
231 # energy fractions for muon, neutrinos and neutralinos
232 add EnergyFraction {12} {0.0}
233 add EnergyFraction {13} {0.0}
234 add EnergyFraction {14} {0.0}
235 add EnergyFraction {16} {0.0}
236 add EnergyFraction {1000022} {0.0}
237 add EnergyFraction {1000023} {0.0}
238 add EnergyFraction {1000025} {0.0}
239 add EnergyFraction {1000035} {0.0}
240 add EnergyFraction {1000045} {0.0}
241 # energy fractions for K0short and Lambda
242 add EnergyFraction {310} {0.3}
243 add EnergyFraction {3122} {0.3}
244
245 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
246
247 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.20^2) }
248
249}
250
251#############
252# HCAL
253#############
254
255module SimpleCalorimeter HCal {
256 set ParticleInputArray ParticlePropagator/stableParticles
257 set TrackInputArray ECal/eflowTracks
258
259 set TowerOutputArray hcalTowers
260 set EFlowTrackOutputArray eflowTracks
261 set EFlowTowerOutputArray eflowNeutralHadrons
262
263 set IsEcal false
264
265 set EnergyMin 1.0
266 set EnergySignificanceMin 1.0
267
268 set SmearTowerCenter true
269
270 set pi [expr {acos(-1)}]
271
272 # lists of the edges of each tower in eta and phi
273 # each list starts with the lower edge of the first tower
274 # the list ends with the higher edged of the last tower
275
276
277 # 2.0 degree towers (6 cm x 6 cm)
278 set PhiBins {}
279 for {set i -180} {$i <= 180} {incr i} {
280 add PhiBins [expr {$i * $pi/180.0}]
281 }
282
283 # 0.04 unit in eta up to eta = 3.0
284 for {set i -75} {$i <= 75} {incr i} {
285 set eta [expr {$i * 0.04}]
286 add EtaPhiBins $eta $PhiBins
287 }
288
289
290 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
291 add EnergyFraction {0} {1.0}
292 # energy fractions for e, gamma and pi0
293 add EnergyFraction {11} {0.0}
294 add EnergyFraction {22} {0.0}
295 add EnergyFraction {111} {0.0}
296 # energy fractions for muon, neutrinos and neutralinos
297 add EnergyFraction {12} {0.0}
298 add EnergyFraction {13} {0.0}
299 add EnergyFraction {14} {0.0}
300 add EnergyFraction {16} {0.0}
301 add EnergyFraction {1000022} {0.0}
302 add EnergyFraction {1000023} {0.0}
303 add EnergyFraction {1000025} {0.0}
304 add EnergyFraction {1000035} {0.0}
305 add EnergyFraction {1000045} {0.0}
306 # energy fractions for K0short and Lambda
307 add EnergyFraction {310} {0.7}
308 add EnergyFraction {3122} {0.7}
309
310 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
311
312 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.60^2)}
313
314}
315
316#################
317# Electron filter
318#################
319
320module PdgCodeFilter ElectronFilter {
321 set InputArray HCal/eflowTracks
322 set OutputArray electrons
323 set Invert true
324 add PdgCode {11}
325 add PdgCode {-11}
326}
327
328###################################################
329# Tower Merger (in case not using e-flow algorithm)
330###################################################
331
332module Merger Calorimeter {
333# add InputArray InputArray
334 add InputArray ECal/ecalTowers
335 add InputArray HCal/hcalTowers
336 set OutputArray towers
337}
338
339
340####################
341# Energy flow merger
342####################
343
344module Merger EFlowMerger {
345# add InputArray InputArray
346 add InputArray HCal/eflowTracks
347 add InputArray ECal/eflowPhotons
348 add InputArray HCal/eflowNeutralHadrons
349 set OutputArray eflow
350}
351
352######################
353# EFlowFilter
354######################
355
356module PdgCodeFilter EFlowFilter {
357 set InputArray EFlowMerger/eflow
358 set OutputArray eflow
359
360 add PdgCode {11}
361 add PdgCode {-11}
362 add PdgCode {13}
363 add PdgCode {-13}
364}
365
366
367###################
368# Missing ET merger
369###################
370
371module Merger MissingET {
372# add InputArray InputArray
373 add InputArray EFlowMerger/eflow
374 set MomentumOutputArray momentum
375}
376
377
378##################
379# Scalar HT merger
380##################
381
382module Merger ScalarHT {
383# add InputArray InputArray
384 add InputArray EFlowMerger/eflow
385 set EnergyOutputArray energy
386}
387
388#################
389# Neutrino Filter
390#################
391
392module PdgCodeFilter NeutrinoFilter {
393
394 set InputArray Delphes/stableParticles
395 set OutputArray filteredParticles
396
397 set PTMin 0.0
398
399 add PdgCode {12}
400 add PdgCode {14}
401 add PdgCode {16}
402 add PdgCode {-12}
403 add PdgCode {-14}
404 add PdgCode {-16}
405
406}
407
408
409#####################
410# MC truth jet finder
411#####################
412
413module FastJetFinder GenJetFinder {
414 set InputArray NeutrinoFilter/filteredParticles
415
416 set OutputArray jets
417
418 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
419 set JetAlgorithm 6
420 set ParameterR 0.5
421
422 set JetPTMin 20.0
423}
424
425#########################
426# Gen Missing ET merger
427########################
428
429module Merger GenMissingET {
430# add InputArray InputArray
431 add InputArray NeutrinoFilter/filteredParticles
432 set MomentumOutputArray momentum
433}
434
435############
436# Jet finder
437############
438
439module FastJetFinder FastJetFinder {
440# set InputArray Calorimeter/towers
441 set InputArray EFlowMerger/eflow
442
443 set OutputArray jets
444
445 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
446 set JetAlgorithm 6
447 set ParameterR 0.5
448
449 set JetPTMin 20.0
450}
451
452##################
453# Jet Energy Scale
454##################
455
456module EnergyScale JetEnergyScale {
457 set InputArray FastJetFinder/jets
458 set OutputArray jets
459
460 # scale formula for jets
461 set ScaleFormula {1.00}
462}
463
464
465########################
466# Jet Flavor Association
467########################
468
469module JetFlavorAssociation JetFlavorAssociation {
470
471 set PartonInputArray Delphes/partons
472 set ParticleInputArray Delphes/allParticles
473 set ParticleLHEFInputArray Delphes/allParticlesLHEF
474 set JetInputArray JetEnergyScale/jets
475
476 set DeltaR 0.5
477 set PartonPTMin 1.0
478 set PartonEtaMax 3.0
479
480}
481
482###################
483# Photon efficiency
484###################
485
486module Efficiency PhotonEfficiency {
487 set InputArray ECal/eflowPhotons
488 set OutputArray photons
489
490 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
491
492 # efficiency formula for photons
493 set EfficiencyFormula { (energy <= 2.0) * (0.00) +
494 (abs(eta) <= 1.5) * (energy > 2.0) * (0.99) +
495 (abs(eta) > 1.5 && abs(eta) <= 3.0) * (energy > 2.0) * (0.99) +
496 (abs(eta) > 3.0) * (0.00)}
497}
498
499##################
500# Photon isolation
501##################
502
503module Isolation PhotonIsolation {
504 set CandidateInputArray PhotonEfficiency/photons
505 set IsolationInputArray EFlowFilter/eflow
506
507 set OutputArray photons
508
509 set DeltaRMax 0.5
510
511 set PTMin 0.5
512
513 set PTRatioMax 0.12
514}
515
516#####################
517# Electron efficiency
518#####################
519
520module Efficiency ElectronEfficiency {
521 set InputArray ElectronFilter/electrons
522 set OutputArray electrons
523
524 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
525
526 # efficiency formula for electrons
527 set EfficiencyFormula { (energy <= 2.0) * (0.00) +
528 (abs(eta) <= 1.5) * (energy > 2.0) * (0.99) +
529 (abs(eta) > 1.5 && abs(eta) <= 3.0) * (energy > 2.0) * (0.99) +
530 (abs(eta) > 3.0) * (0.00)}
531}
532
533####################
534# Electron isolation
535####################
536
537module Isolation ElectronIsolation {
538 set CandidateInputArray ElectronEfficiency/electrons
539 set IsolationInputArray EFlowFilter/eflow
540
541 set OutputArray electrons
542
543 set DeltaRMax 0.5
544
545 set PTMin 0.5
546
547 set PTRatioMax 0.12
548}
549
550#################
551# Muon efficiency
552#################
553
554module Efficiency MuonEfficiency {
555 set InputArray MuonMomentumSmearing/muons
556 set OutputArray muons
557
558 # set EfficiencyFormula {efficiency as a function of eta and pt}
559 set EfficiencyFormula { (energy <= 2.0) * (0.00) +
560 (abs(eta) <= 1.5) * (energy > 2.0) * (0.99) +
561 (abs(eta) > 1.5 && abs(eta) <= 3.0) * (energy > 2.0) * (0.99) +
562 (abs(eta) > 3.0) * (0.00)}
563}
564
565################
566# Muon isolation
567################
568
569module Isolation MuonIsolation {
570 set CandidateInputArray MuonEfficiency/muons
571 set IsolationInputArray EFlowFilter/eflow
572
573 set OutputArray muons
574
575 set DeltaRMax 0.5
576
577 set PTMin 0.5
578
579 set PTRatioMax 0.25
580}
581
582
583###########
584# b-tagging
585###########
586
587module BTagging BTagging {
588 set JetInputArray JetEnergyScale/jets
589
590 set BitNumber 0
591
592 # default efficiency formula (misidentification rate)
593 add EfficiencyFormula {0} {0.001}
594
595 # efficiency formula for c-jets (misidentification rate)
596 add EfficiencyFormula {4} {0.10}
597
598 # efficiency formula for b-jets
599 add EfficiencyFormula {5} {0.80}
600}
601
602###########
603# c-tagging
604###########
605
606module BTagging CTagging {
607 set JetInputArray JetEnergyScale/jets
608
609 set BitNumber 1
610
611 # default efficiency formula (misidentification rate)
612 add EfficiencyFormula {0} {0.12}
613
614 # efficiency formula for c-jets (misidentification rate)
615 add EfficiencyFormula {4} {0.70}
616
617 # efficiency formula for b-jets
618 add EfficiencyFormula {5} {0.20}
619}
620
621
622#############
623# tau-tagging
624#############
625
626
627module TauTagging TauTagging {
628 set ParticleInputArray Delphes/allParticles
629 set PartonInputArray Delphes/partons
630 set JetInputArray JetEnergyScale/jets
631
632 set DeltaR 0.5
633
634 set TauPTMin 1.0
635
636 set TauEtaMax 3.0
637
638 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
639
640 # default efficiency formula (misidentification rate)
641 add EfficiencyFormula {0} {0.001}
642 # efficiency formula for tau-jets
643 add EfficiencyFormula {15} {0.4}
644}
645
646#####################################################
647# Find uniquely identified photons/electrons/tau/jets
648#####################################################
649
650module UniqueObjectFinder UniqueObjectFinder {
651# earlier arrays take precedence over later ones
652# add InputArray InputArray OutputArray
653 add InputArray PhotonIsolation/photons photons
654 add InputArray ElectronIsolation/electrons electrons
655 add InputArray MuonIsolation/muons muons
656 add InputArray JetEnergyScale/jets jets
657}
658
659
660##################
661# ROOT tree writer
662##################
663
664module TreeWriter TreeWriter {
665# add Branch InputArray BranchName BranchClass
666 add Branch Delphes/allParticles Particle GenParticle
667
668 add Branch GenJetFinder/jets GenJet Jet
669 add Branch GenMissingET/momentum GenMissingET MissingET
670
671 add Branch TrackMerger/tracks Track Track
672 add Branch Calorimeter/towers Tower Tower
673
674 add Branch HCal/eflowTracks EFlowTrack Track
675 add Branch ECal/eflowPhotons EFlowPhoton Tower
676 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
677
678 add Branch UniqueObjectFinder/photons Photon Photon
679 add Branch UniqueObjectFinder/electrons Electron Electron
680 add Branch UniqueObjectFinder/muons Muon Muon
681 add Branch UniqueObjectFinder/jets Jet Jet
682
683 add Branch MissingET/momentum MissingET MissingET
684 add Branch ScalarHT/energy ScalarHT ScalarHT
685}
686
Note: See TracBrowser for help on using the repository browser.