Fork me on GitHub

source: git/cards/delphes_card_FCC_basic.tcl@ d77b51d

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

fixed typo in Hcal resolution formula
  • 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 EFlowTowerOutputArray eflowPhotons
227
228 set IsEcal true
229
230 set EnergyMin 0.5
231 set EnergySignificanceMin 1.0
232
233 set SmearTowerCenter true
234
235 set pi [expr {acos(-1)}]
236
237 # lists of the edges of each tower in eta and phi
238 # each list starts with the lower edge of the first tower
239 # the list ends with the higher edged of the last tower
240
241 # 0.5 degree towers
242 set PhiBins {}
243 for {set i -360} {$i <= 360} {incr i} {
244 add PhiBins [expr {$i * $pi/360.0}]
245 }
246
247 # 0.01 unit in eta up to eta = 3
248 for {set i -300} {$i <= 300} {incr i} {
249 set eta [expr {$i * 0.01}]
250 add EtaPhiBins $eta $PhiBins
251 }
252
253 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
254
255 add EnergyFraction {0} {0.0}
256 # energy fractions for e, gamma and pi0
257 add EnergyFraction {11} {1.0}
258 add EnergyFraction {22} {1.0}
259 add EnergyFraction {111} {1.0}
260 # energy fractions for muon, neutrinos and neutralinos
261 add EnergyFraction {12} {0.0}
262 add EnergyFraction {13} {0.0}
263 add EnergyFraction {14} {0.0}
264 add EnergyFraction {16} {0.0}
265 add EnergyFraction {1000022} {0.0}
266 add EnergyFraction {1000023} {0.0}
267 add EnergyFraction {1000025} {0.0}
268 add EnergyFraction {1000035} {0.0}
269 add EnergyFraction {1000045} {0.0}
270 # energy fractions for K0short and Lambda
271 add EnergyFraction {310} {0.3}
272 add EnergyFraction {3122} {0.3}
273
274 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
275
276 # This is the CMS ECAL resolution, extended up eta = 6.0
277 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2) +
278 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
279
280
281}
282
283#############
284# HCAL
285#############
286
287module SimpleCalorimeter HCal {
288 set ParticleInputArray ParticlePropagator/stableParticles
289 set TrackInputArray ImpactParameterSmearing/tracks
290
291 set TowerOutputArray hcalTowers
292 set EFlowTowerOutputArray eflowNeutralHadrons
293
294 set IsEcal false
295
296 set EnergyMin 1.0
297 set EnergySignificanceMin 1.0
298
299 set SmearTowerCenter true
300
301 set pi [expr {acos(-1)}]
302
303 # lists of the edges of each tower in eta and phi
304 # each list starts with the lower edge of the first tower
305 # the list ends with the higher edged of the last tower
306
307
308 # 5 degree towers
309 set PhiBins {}
310 for {set i -72} {$i <= 72} {incr i} {
311 add PhiBins [expr {$i * $pi/75.0}]
312 }
313
314 # 0.05 unit in eta up to eta = 3
315 for {set i -60} {$i <= 60} {incr i} {
316 set eta [expr {$i * 0.05}]
317 add EtaPhiBins $eta $PhiBins
318 }
319
320
321 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
322 add EnergyFraction {0} {1.0}
323 # energy fractions for e, gamma and pi0
324 add EnergyFraction {11} {0.0}
325 add EnergyFraction {22} {0.0}
326 add EnergyFraction {111} {0.0}
327 # energy fractions for muon, neutrinos and neutralinos
328 add EnergyFraction {12} {0.0}
329 add EnergyFraction {13} {0.0}
330 add EnergyFraction {14} {0.0}
331 add EnergyFraction {16} {0.0}
332 add EnergyFraction {1000022} {0.0}
333 add EnergyFraction {1000023} {0.0}
334 add EnergyFraction {1000025} {0.0}
335 add EnergyFraction {1000035} {0.0}
336 add EnergyFraction {1000045} {0.0}
337 # energy fractions for K0short and Lambda
338 add EnergyFraction {310} {0.7}
339 add EnergyFraction {3122} {0.7}
340
341 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
342
343 # This is the ATLAS HCAL resolution, extended up eta = 6.0
344 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
345 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
346 (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
347
348}
349
350####################
351# Tower Merger (in case not using e-flow algorithm)
352####################
353
354module Merger TowerMerger {
355# add InputArray InputArray
356 add InputArray ECal/ecalTowers
357 add InputArray HCal/hcalTowers
358 set OutputArray towers
359}
360
361####################
362# Energy flow merger
363####################
364
365module Merger EFlowMerger {
366# add InputArray InputArray
367 add InputArray ImpactParameterSmearing/tracks
368 add InputArray ECal/eflowPhotons
369 add InputArray HCal/eflowNeutralHadrons
370 set OutputArray eflow
371}
372
373
374###################
375# Missing ET merger
376###################
377
378module Merger MissingET {
379# add InputArray InputArray
380 add InputArray EFlowMerger/eflow
381 set MomentumOutputArray momentum
382}
383
384
385##################
386# Scalar HT merger
387##################
388
389module Merger ScalarHT {
390# add InputArray InputArray
391 add InputArray EFlowMerger/eflow
392 set EnergyOutputArray energy
393}
394
395#####################
396# Neutrino Filter
397#####################
398
399module PdgCodeFilter NeutrinoFilter {
400
401 set InputArray Delphes/stableParticles
402 set OutputArray filteredParticles
403
404 set PTMin 0.0
405
406 add PdgCode {12}
407 add PdgCode {14}
408 add PdgCode {16}
409 add PdgCode {-12}
410 add PdgCode {-14}
411 add PdgCode {-16}
412
413}
414
415
416#####################
417# MC truth jet finder
418#####################
419
420module FastJetFinder GenJetFinder {
421 set InputArray NeutrinoFilter/filteredParticles
422
423 set OutputArray jets
424
425 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
426 set JetAlgorithm 6
427 set ParameterR 0.5
428
429 set JetPTMin 100.0
430}
431
432############
433# Jet finder
434############
435
436module FastJetFinder FastJetFinder {
437# set InputArray TowerMerger/towers
438 set InputArray EFlowMerger/eflow
439
440 set OutputArray jets
441
442 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
443 set JetAlgorithm 6
444 set ParameterR 0.5
445
446 set JetPTMin 100.0
447}
448
449##################
450# Jet Energy Scale
451##################
452
453module EnergyScale JetEnergyScale {
454 set InputArray FastJetFinder/jets
455 set OutputArray jets
456
457 # scale formula for jets
458 set ScaleFormula {1.00}
459}
460
461##########################
462# Track Counting b-tagging
463##########################
464
465module TrackCountingBTagging TrackCountingBTagging {
466 set TrackInputArray ImpactParameterSmearing/tracks
467 set JetInputArray JetEnergyScale/jets
468
469 set BitNumber 0
470
471 # maximum distance between jet and track
472 set DeltaR 0.3
473
474 # minimum pt of tracks
475 set TrackPTMin 1.0
476
477 # minimum transverse impact parameter (in mm)
478 set TrackIPMax 2.0
479
480 # minimum ip significance for the track to be counted
481 set SigMin 6.5
482
483 # minimum number of tracks (high efficiency n=2, high purity n=3)
484 set Ntracks 3
485}
486
487
488##########################
489# tau-tagging
490##########################
491
492
493module TauTagging TauTagging {
494 set ParticleInputArray Delphes/allParticles
495 set PartonInputArray Delphes/partons
496 set JetInputArray JetEnergyScale/jets
497
498 set DeltaR 0.5
499
500 set TauPTMin 1.0
501
502 set TauEtaMax 4.0
503
504 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
505
506 # default efficiency formula (misidentification rate)
507 add EfficiencyFormula {0} {0.001}
508 # efficiency formula for tau-jets
509 add EfficiencyFormula {15} {0.4}
510}
511
512##################
513# ROOT tree writer
514##################
515
516module TreeWriter TreeWriter {
517# add Branch InputArray BranchName BranchClass
518 add Branch Delphes/allParticles Particle GenParticle
519 add Branch GenJetFinder/jets GenJet Jet
520
521 add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
522 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
523 add Branch ECal/eflowPhotons Photon Photon
524
525 add Branch ElectronEnergySmearing/electrons Electron Electron
526 add Branch MuonMomentumSmearing/muons Muon Muon
527 add Branch JetEnergyScale/jets Jet Jet
528 add Branch MissingET/momentum MissingET MissingET
529 add Branch ScalarHT/energy ScalarHT ScalarHT
530}
Note: See TracBrowser for help on using the repository browser.