Fork me on GitHub

source: git/cards/delphes_card_FCC_basic.tcl@ 914fb04

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 914fb04 was fa7f333, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

fixed eflow at high pt
  • Property mode set to 100644
File size: 14.5 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 ElectronEnergySmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17 AngularSmearing
18 ImpactParameterSmearing
19
20 ECal
21 HCal
22
23 TowerMerger
24 EFlowMerger
25
26 MissingET
27
28 NeutrinoFilter
29 GenJetFinder
30 FastJetFinder
31
32 JetEnergyScale
33
34 TrackCountingBTagging
35 TauTagging
36
37 ScalarHT
38
39 TreeWriter
40}
41
42#################################
43# Propagate particles in cylinder
44#################################
45
46module ParticlePropagator ParticlePropagator {
47 set InputArray Delphes/stableParticles
48
49 set OutputArray stableParticles
50 set ChargedHadronOutputArray chargedHadrons
51 set ElectronOutputArray electrons
52 set MuonOutputArray muons
53
54 # radius of the magnetic field coverage, in m
55 set Radius 2.60
56 # half-length of the magnetic field coverage, in m
57 set HalfLength 6.00
58
59 # magnetic field
60 set Bz 6.0
61}
62
63####################################
64# Charged hadron tracking efficiency
65####################################
66
67module Efficiency ChargedHadronTrackingEfficiency {
68 set InputArray ParticlePropagator/chargedHadrons
69 set OutputArray chargedHadrons
70
71 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
72
73 # tracking efficiency formula for charged hadrons
74 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
75 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
76 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
77 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.60) +
78 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.90) +
79 (abs(eta) > 4.0) * (0.00)}
80
81}
82
83##############################
84# Electron tracking efficiency
85##############################
86
87module Efficiency ElectronTrackingEfficiency {
88 set InputArray ParticlePropagator/electrons
89 set OutputArray electrons
90
91 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
92
93 # tracking efficiency formula for electrons
94 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
95 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
96 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
97 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) +
98 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) +
99 (abs(eta) > 4.0) * (0.00)}
100}
101
102##########################
103# Muon tracking efficiency
104##########################
105
106module Efficiency MuonTrackingEfficiency {
107 set InputArray ParticlePropagator/muons
108 set OutputArray muons
109
110 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
111
112 # tracking efficiency formula for muons
113 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
114 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
115 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
116 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) +
117 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.98) +
118 (abs(eta) > 4.0) * (0.00)}
119}
120
121########################################
122# Momentum resolution for charged tracks
123########################################
124
125module MomentumSmearing ChargedHadronMomentumSmearing {
126 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
127 set OutputArray chargedHadrons
128
129 # set ResolutionFormula {resolution formula as a function of eta and pt}
130
131 # resolution formula for charged hadrons
132 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*2.e-5) +
133 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)}
134
135
136}
137
138#################################
139# Energy resolution for electrons
140#################################
141
142module EnergySmearing ElectronEnergySmearing {
143 set InputArray ElectronTrackingEfficiency/electrons
144 set OutputArray electrons
145
146 # set ResolutionFormula {resolution formula as a function of eta and energy}
147
148 # resolution formula for electrons
149 set ResolutionFormula { (abs(eta) <= 4.0) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.007) +
150 (abs(eta) <= 4.0) * (energy > 2.0e1) * sqrt(energy^2*0.005^2 + energy*0.02^2) +
151 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
152
153}
154
155###############################
156# Momentum resolution for muons
157###############################
158
159module MomentumSmearing MuonMomentumSmearing {
160 set InputArray MuonTrackingEfficiency/muons
161 set OutputArray muons
162
163 # set ResolutionFormula {resolution formula as a function of eta and pt}
164
165 # resolution formula for muons
166 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1) * (0.01 + pt*5.e-6) +
167 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.e-5)}
168
169
170}
171
172##############
173# Track merger
174##############
175
176module Merger TrackMerger {
177# add InputArray InputArray
178 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
179 add InputArray ElectronEnergySmearing/electrons
180 add InputArray MuonMomentumSmearing/muons
181 set OutputArray tracks
182}
183
184
185################################
186# Track angular smearing
187################################
188
189module AngularSmearing AngularSmearing {
190 set InputArray TrackMerger/tracks
191 set OutputArray tracks
192
193
194 # angular smearing in eta formula as a function of pt and eta
195 set EtaResolutionFormula { 0.001 }
196
197 # angular smearing in phi formula as a function of pt and eta
198 set PhiResolutionFormula { 0.001 }
199
200}
201
202################################
203# Track impact parameter smearing
204################################
205
206module ImpactParameterSmearing ImpactParameterSmearing {
207 set InputArray AngularSmearing/tracks
208 set OutputArray tracks
209
210
211 # absolute impact parameter smearing formula (in mm) as a function of pt and eta
212 set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) +
213 (pt > 5.0) * (0.005)}
214
215}
216
217#############
218# ECAL
219#############
220
221module SimpleCalorimeter ECal {
222 set ParticleInputArray ParticlePropagator/stableParticles
223 set TrackInputArray ImpactParameterSmearing/tracks
224
225 set TowerOutputArray ecalTowers
226 set EFlowTrackOutputArray eflowTracks
227 set EFlowTowerOutputArray eflowPhotons
228
229 set IsEcal true
230
231 set EnergyMin 0.5
232 set EnergySignificanceMin 1.0
233
234 set SmearTowerCenter true
235
236 set pi [expr {acos(-1)}]
237
238 # lists of the edges of each tower in eta and phi
239 # each list starts with the lower edge of the first tower
240 # the list ends with the higher edged of the last tower
241
242 # 0.5 degree towers
243 set PhiBins {}
244 for {set i -360} {$i <= 360} {incr i} {
245 add PhiBins [expr {$i * $pi/360.0}]
246 }
247
248 # 0.01 unit in eta up to eta = 3
249 for {set i -300} {$i <= 300} {incr i} {
250 set eta [expr {$i * 0.01}]
251 add EtaPhiBins $eta $PhiBins
252 }
253
254 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
255
256 add EnergyFraction {0} {0.0}
257 # energy fractions for e, gamma and pi0
258 add EnergyFraction {11} {1.0}
259 add EnergyFraction {22} {1.0}
260 add EnergyFraction {111} {1.0}
261 # energy fractions for muon, neutrinos and neutralinos
262 add EnergyFraction {12} {0.0}
263 add EnergyFraction {13} {0.0}
264 add EnergyFraction {14} {0.0}
265 add EnergyFraction {16} {0.0}
266 add EnergyFraction {1000022} {0.0}
267 add EnergyFraction {1000023} {0.0}
268 add EnergyFraction {1000025} {0.0}
269 add EnergyFraction {1000035} {0.0}
270 add EnergyFraction {1000045} {0.0}
271 # energy fractions for K0short and Lambda
272 add EnergyFraction {310} {0.3}
273 add EnergyFraction {3122} {0.3}
274
275 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
276
277 # This is the CMS ECAL resolution, extended up eta = 6.0
278 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2) +
279 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
280
281
282}
283
284#############
285# HCAL
286#############
287
288module SimpleCalorimeter HCal {
289 set ParticleInputArray ParticlePropagator/stableParticles
290 set TrackInputArray ECal/eflowTracks
291
292 set TowerOutputArray hcalTowers
293 set EFlowTrackOutputArray eflowTracks
294 set EFlowTowerOutputArray eflowNeutralHadrons
295
296 set IsEcal false
297
298 set EnergyMin 1.0
299 set EnergySignificanceMin 1.0
300
301 set SmearTowerCenter true
302
303 set pi [expr {acos(-1)}]
304
305 # lists of the edges of each tower in eta and phi
306 # each list starts with the lower edge of the first tower
307 # the list ends with the higher edged of the last tower
308
309
310 # 5 degree towers
311 set PhiBins {}
312 for {set i -72} {$i <= 72} {incr i} {
313 add PhiBins [expr {$i * $pi/75.0}]
314 }
315
316 # 0.05 unit in eta up to eta = 3
317 for {set i -60} {$i <= 60} {incr i} {
318 set eta [expr {$i * 0.05}]
319 add EtaPhiBins $eta $PhiBins
320 }
321
322
323 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
324 add EnergyFraction {0} {1.0}
325 # energy fractions for e, gamma and pi0
326 add EnergyFraction {11} {0.0}
327 add EnergyFraction {22} {0.0}
328 add EnergyFraction {111} {0.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.7}
341 add EnergyFraction {3122} {0.7}
342
343 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
344
345 # This is the ATLAS HCAL resolution, extended up eta = 6.0
346 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
347 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
348 (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
349
350}
351
352####################
353# Tower Merger (in case not using e-flow algorithm)
354####################
355
356module Merger TowerMerger {
357# add InputArray InputArray
358 add InputArray ECal/ecalTowers
359 add InputArray HCal/hcalTowers
360 set OutputArray towers
361}
362
363####################
364# Energy flow merger
365####################
366
367module Merger EFlowMerger {
368# add InputArray InputArray
369 add InputArray HCal/eflowTracks
370 add InputArray ECal/eflowPhotons
371 add InputArray HCal/eflowNeutralHadrons
372 set OutputArray eflow
373}
374
375
376###################
377# Missing ET merger
378###################
379
380module Merger MissingET {
381# add InputArray InputArray
382 add InputArray EFlowMerger/eflow
383 set MomentumOutputArray momentum
384}
385
386
387##################
388# Scalar HT merger
389##################
390
391module Merger ScalarHT {
392# add InputArray InputArray
393 add InputArray EFlowMerger/eflow
394 set EnergyOutputArray energy
395}
396
397#####################
398# Neutrino Filter
399#####################
400
401module PdgCodeFilter NeutrinoFilter {
402
403 set InputArray Delphes/stableParticles
404 set OutputArray filteredParticles
405
406 set PTMin 0.0
407
408 add PdgCode {12}
409 add PdgCode {14}
410 add PdgCode {16}
411 add PdgCode {-12}
412 add PdgCode {-14}
413 add PdgCode {-16}
414
415}
416
417
418#####################
419# MC truth jet finder
420#####################
421
422module FastJetFinder GenJetFinder {
423 set InputArray NeutrinoFilter/filteredParticles
424
425 set OutputArray jets
426
427 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
428 set JetAlgorithm 6
429 set ParameterR 0.5
430
431 set JetPTMin 100.0
432}
433
434############
435# Jet finder
436############
437
438module FastJetFinder FastJetFinder {
439# set InputArray TowerMerger/towers
440 set InputArray EFlowMerger/eflow
441
442 set OutputArray jets
443
444 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
445 set JetAlgorithm 6
446 set ParameterR 0.5
447
448 set JetPTMin 100.0
449}
450
451##################
452# Jet Energy Scale
453##################
454
455module EnergyScale JetEnergyScale {
456 set InputArray FastJetFinder/jets
457 set OutputArray jets
458
459 # scale formula for jets
460 set ScaleFormula {1.00}
461}
462
463##########################
464# Track Counting b-tagging
465##########################
466
467module TrackCountingBTagging TrackCountingBTagging {
468 set TrackInputArray ImpactParameterSmearing/tracks
469 set JetInputArray JetEnergyScale/jets
470
471 set BitNumber 0
472
473 # maximum distance between jet and track
474 set DeltaR 0.3
475
476 # minimum pt of tracks
477 set TrackPTMin 1.0
478
479 # minimum transverse impact parameter (in mm)
480 set TrackIPMax 2.0
481
482 # minimum ip significance for the track to be counted
483 set SigMin 6.5
484
485 # minimum number of tracks (high efficiency n=2, high purity n=3)
486 set Ntracks 3
487}
488
489
490##########################
491# tau-tagging
492##########################
493
494
495module TauTagging TauTagging {
496 set ParticleInputArray Delphes/allParticles
497 set PartonInputArray Delphes/partons
498 set JetInputArray JetEnergyScale/jets
499
500 set DeltaR 0.5
501
502 set TauPTMin 1.0
503
504 set TauEtaMax 4.0
505
506 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
507
508 # default efficiency formula (misidentification rate)
509 add EfficiencyFormula {0} {0.001}
510 # efficiency formula for tau-jets
511 add EfficiencyFormula {15} {0.4}
512}
513
514##################
515# ROOT tree writer
516##################
517
518module TreeWriter TreeWriter {
519# add Branch InputArray BranchName BranchClass
520 add Branch Delphes/allParticles Particle GenParticle
521 add Branch GenJetFinder/jets GenJet Jet
522
523 add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
524 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
525 add Branch ECal/eflowPhotons Photon Photon
526
527 add Branch ElectronEnergySmearing/electrons Electron Electron
528 add Branch MuonMomentumSmearing/muons Muon Muon
529 add Branch JetEnergyScale/jets Jet Jet
530 add Branch MissingET/momentum MissingET MissingET
531 add Branch ScalarHT/energy ScalarHT ScalarHT
532}
Note: See TracBrowser for help on using the repository browser.