MA5SandBox: delphes_card_atlas_susy_2018_32.tcl

File delphes_card_atlas_susy_2018_32.tcl, 18.8 KB (added by Benjamin Fuks, 4 years ago)
Line 
1#######################################
2# MODIF OF STANDARD DELFES CARD #
3#######################################
4
5# => Medium electron 2013 efficiency implemented
6# => 2013 efficiency for muons implemented
7# => Efficiency for b-tag, c-jets, light jets taken
8# from ATLAS 0-leptons + 2b-jets (ATLAS SUSY 13 05)
9
10
11#######################################
12# Order of execution of various modules
13#######################################
14
15set ExecutionPath {
16 ParticlePropagator
17
18 ChargedHadronTrackingEfficiency
19 ElectronTrackingEfficiency
20 MuonTrackingEfficiency
21
22 ChargedHadronMomentumSmearing
23 ElectronEnergySmearing
24 MuonMomentumSmearing
25
26 TrackMerger
27 Calorimeter
28 EFlowMerger
29
30 PhotonEfficiency
31
32 ElectronEfficiency
33
34 MuonEfficiency
35
36 MissingET
37
38 GenJetFinder
39 FastJetFinder
40
41 JetEnergyScale
42
43 JetFlavorAssociation
44
45 BTagging
46 TauTagging
47
48 UniqueObjectFinder
49
50
51 TreeWriter
52}
53
54
55
56#################################
57# Propagate particles in cylinder
58#################################
59
60module ParticlePropagator ParticlePropagator {
61 set InputArray Delphes/stableParticles
62
63 set OutputArray stableParticles
64 set ChargedHadronOutputArray chargedHadrons
65 set ElectronOutputArray electrons
66 set MuonOutputArray muons
67
68 # radius of the magnetic field coverage, in m
69 set Radius 1.15
70 # half-length of the magnetic field coverage, in m
71 set HalfLength 3.51
72
73 # magnetic field
74 set Bz 2.0
75}
76
77####################################
78# Charged hadron tracking efficiency
79####################################
80
81module Efficiency ChargedHadronTrackingEfficiency {
82 set InputArray ParticlePropagator/chargedHadrons
83 set OutputArray chargedHadrons
84
85 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
86
87 # tracking efficiency formula for charged hadrons
88 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
89 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
90 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
91 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
92 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
93 (abs(eta) > 2.5) * (0.00)}
94}
95
96##############################
97# Electron tracking efficiency
98##############################
99
100module Efficiency ElectronTrackingEfficiency {
101 set InputArray ParticlePropagator/electrons
102 set OutputArray electrons
103
104 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
105
106 # tracking efficiency formula for electrons
107 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
108 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
109 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
110 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
111 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
112 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
113 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
114 (abs(eta) > 2.5) * (0.00)}
115}
116
117##########################
118# Muon tracking efficiency
119##########################
120
121module Efficiency MuonTrackingEfficiency {
122 set InputArray ParticlePropagator/muons
123 set OutputArray muons
124
125 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
126
127 # tracking efficiency formula for muons
128 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
129 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
130 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
131 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
132 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
133 (abs(eta) > 2.5) * (0.00)}
134}
135
136########################################
137# Momentum resolution for charged tracks
138########################################
139
140module MomentumSmearing ChargedHadronMomentumSmearing {
141 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
142 set OutputArray chargedHadrons
143
144 # set ResolutionFormula {resolution formula as a function of eta and pt}
145
146 # resolution formula for charged hadrons
147 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
148 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
149 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
150 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
151 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
152 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
153 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
154 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
155}
156
157#################################
158# Energy resolution for electrons
159#################################
160
161module EnergySmearing ElectronEnergySmearing {
162 set InputArray ElectronTrackingEfficiency/electrons
163 set OutputArray electrons
164
165 # set ResolutionFormula {resolution formula as a function of eta and energy}
166
167 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
168 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
169 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
170 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
171
172}
173
174###############################
175# Momentum resolution for muons
176###############################
177
178module MomentumSmearing MuonMomentumSmearing {
179 set InputArray MuonTrackingEfficiency/muons
180 set OutputArray muons
181
182 # set ResolutionFormula {resolution formula as a function of eta and pt}
183
184 # resolution formula for muons
185 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
186 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \
187 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
188 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \
189 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \
190 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \
191 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
192 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)}
193}
194
195##############
196# Track merger
197##############
198
199module Merger TrackMerger {
200# add InputArray InputArray
201 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
202 add InputArray ElectronEnergySmearing/electrons
203 add InputArray MuonMomentumSmearing/muons
204 set OutputArray tracks
205}
206
207#############
208# Calorimeter
209#############
210
211module Calorimeter Calorimeter {
212 set ParticleInputArray ParticlePropagator/stableParticles
213 set TrackInputArray TrackMerger/tracks
214
215 set TowerOutputArray towers
216 set PhotonOutputArray photons
217
218 set EFlowTrackOutputArray eflowTracks
219 set EFlowPhotonOutputArray eflowPhotons
220 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
221
222 set pi [expr {acos(-1)}]
223
224 # lists of the edges of each tower in eta and phi
225 # each list starts with the lower edge of the first tower
226 # the list ends with the higher edged of the last tower
227
228 # 10 degrees towers
229 set PhiBins {}
230 for {set i -18} {$i <= 18} {incr i} {
231 add PhiBins [expr {$i * $pi/18.0}]
232 }
233 foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
234 add EtaPhiBins $eta $PhiBins
235 }
236
237 # 20 degrees towers
238 set PhiBins {}
239 for {set i -9} {$i <= 9} {incr i} {
240 add PhiBins [expr {$i * $pi/9.0}]
241 }
242 foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
243 add EtaPhiBins $eta $PhiBins
244 }
245
246 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
247 add EnergyFraction {0} {0.0 1.0}
248 # energy fractions for e, gamma and pi0
249 add EnergyFraction {11} {1.0 0.0}
250 add EnergyFraction {22} {1.0 0.0}
251 add EnergyFraction {111} {1.0 0.0}
252 # energy fractions for muon, neutrinos and neutralinos
253 add EnergyFraction {12} {0.0 0.0}
254 add EnergyFraction {13} {0.0 0.0}
255 add EnergyFraction {14} {0.0 0.0}
256 add EnergyFraction {16} {0.0 0.0}
257 add EnergyFraction {1000022} {0.0 0.0}
258 add EnergyFraction {1000023} {0.0 0.0}
259 add EnergyFraction {1000025} {0.0 0.0}
260 add EnergyFraction {1000035} {0.0 0.0}
261 add EnergyFraction {1000045} {0.0 0.0}
262 # energy fractions for K0short and Lambda
263 add EnergyFraction {310} {0.3 0.7}
264 add EnergyFraction {3122} {0.3 0.7}
265
266 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
267 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
268 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
269 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
270 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
271 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
272
273 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
274 # http://arxiv.org/pdf/hep-ex/0004009v1
275 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
276 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
277 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
278 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
279}
280
281####################
282# Energy flow merger
283####################
284
285module Merger EFlowMerger {
286# add InputArray InputArray
287 add InputArray Calorimeter/eflowTracks
288 add InputArray Calorimeter/eflowPhotons
289 add InputArray Calorimeter/eflowNeutralHadrons
290 set OutputArray eflow
291}
292
293###################
294# Photon efficiency
295###################
296
297module Efficiency PhotonEfficiency {
298 set InputArray Calorimeter/photons
299 set OutputArray photons
300
301 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
302
303 # efficiency formula for photons
304 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
305 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
306 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
307 (abs(eta) > 2.5) * (0.00)}
308}
309
310##################
311# Photon isolation
312##################
313
314module Isolation PhotonIsolation {
315 set CandidateInputArray PhotonEfficiency/photons
316 set IsolationInputArray EFlowMerger/eflow
317
318 set OutputArray photons
319
320 set DeltaRMax 0.5
321
322 set PTMin 0.5
323
324 set PTRatioMax 0.1
325}
326
327#####################
328# Electron efficiency
329#####################
330
331module Efficiency ElectronEfficiency {
332 set InputArray ElectronEnergySmearing/electrons
333 set OutputArray electrons
334
335 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
336
337 # efficiency formula for "tight" electrons retreived from https://arxiv.org/abs/1908.00005
338 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
339 (abs(eta) <= 1.5) * (pt > 5.0) * (.95) + \
340 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0) * (.85) + \
341 (abs(eta) > 2.5)* (0.00)}
342# set EfficiencyFormula {
343# (pt <= 4.5) * (0.00) +\
344# (abs(eta) > 2.5 ) * (0.00) +\
345# (abs(eta) < 0.5 && pt > 4.5) * 0.778 +\
346# (abs(eta)>=0.5 && abs(eta) <1 && pt > 4.5) * 0.778 +\
347# (abs(eta)>=1 && abs(eta) <1.5 && pt > 4.5) * 0.724 +\
348# (abs(eta)>=1.5 && abs(eta) <2 && pt > 4.5) * 0.75 +\
349# (abs(eta)>=2 && abs(eta)<=2.5 && pt > 4.5) * 0.6985
350# }
351}
352
353# }
354
355
356#################
357# Muon efficiency
358#################
359
360module Efficiency MuonEfficiency {
361 set InputArray MuonMomentumSmearing/muons
362 set OutputArray muons
363
364 # set EfficiencyFormula {efficiency as a function of eta and pt}
365 #default efficiency for muons # efficiency formula for muons
366 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
367 (abs(eta) <= 1.5) * (pt > 5.0) * (.95) + \
368 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 5.0) * (.85) + \
369 (abs(eta) > 2.7)* (0.00)}
370# set EfficiencyFormula {
371# (pt <= 6.0) * (0.00) + \
372# (abs(eta) <= 0.1) * (pt > 6.0) * (0.573) + \
373# (abs(eta) > 0.1 && abs(eta) <=1.05) * (pt > 6.0) * (0.947) + \
374# (abs(eta) > 1.05 && abs(eta) <=1.5) * (pt > 6.0) * (0.943) + \
375# (abs(eta) > 1.5 && abs(eta) <=2.) * (pt > 6.0) * (0.949) + \
376# (abs(eta) > 2. && abs(eta) <=2.5) * (pt > 6.0) * (0.949) + \
377# (abs(eta) > 2.5) * (0.00)}
378}
379
380# }
381
382
383###################
384# Missing ET merger
385###################
386
387module Merger MissingET {
388# add InputArray InputArray
389 add InputArray EFlowMerger/eflow
390 set MomentumOutputArray momentum
391}
392
393
394#####################
395# MC truth jet finder
396#####################
397
398module FastJetFinder GenJetFinder {
399 set InputArray Delphes/stableParticles
400
401 set OutputArray jets
402
403 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
404 set JetAlgorithm 6
405 set ParameterR 0.4
406
407 set JetPTMin 20.0
408}
409
410############
411# Jet finder
412############
413
414module FastJetFinder FastJetFinder {
415 set InputArray Calorimeter/towers
416
417 set OutputArray jets
418
419 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
420 set JetAlgorithm 6
421 set ParameterR 0.4
422
423 set JetPTMin 20.0
424}
425
426##################
427# Jet Energy Scale
428##################
429
430module EnergyScale JetEnergyScale {
431 set InputArray FastJetFinder/jets
432 set OutputArray jets
433
434 # scale formula for jets
435 set ScaleFormula {1.08}
436}
437
438########################
439# Jet Flavor Association
440########################
441
442module JetFlavorAssociation JetFlavorAssociation {
443
444 set PartonInputArray Delphes/partons
445 set ParticleInputArray Delphes/allParticles
446 set ParticleLHEFInputArray Delphes/allParticlesLHEF
447 set JetInputArray JetEnergyScale/jets
448
449 set DeltaR 0.5
450 set PartonPTMin 1.0
451 set PartonEtaMax 6.0
452
453}
454
455
456
457###########
458# b-tagging
459###########
460
461module BTagging BTagging {
462 set PartonInputArray Delphes/partons
463 set JetInputArray JetEnergyScale/jets
464
465 set BitNumber 0
466
467 set DeltaR 0.5
468
469 set PartonPTMin 1.0
470
471 set PartonEtaMax 2.5
472
473 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
474 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
475 # gluon's PDG code has the lowest priority
476
477 # default efficiency formula (misidentification rate)
478 add EfficiencyFormula {0} { (pt <= 20.0) * (0.00) + \
479 (abs(eta) <= 2.5) * (pt >= 20.0) * (-2.657e-06*pt + 9.11e-09*pt^2 - 6.19e-3*log(pt)/pt+2.02e-3) + \
480 (abs(eta) >= 2.5) * (0.00)}
481
482 # efficiency formula for c-jets (misidentification rate)
483 add EfficiencyFormula {4} { (pt <= 20.0) * (0.00) + \
484 (abs(eta) <= 2.5) * (pt >= 20.0) * (-4.533e-04*pt + 2.595e-07*pt^2 - 1.243*log(pt)/pt+0.249) + \
485 (abs(eta) >= 2.5) * (0.00)}
486
487 # efficiency formula for b-jets
488 add EfficiencyFormula {5} { (pt <= 20.0) * (0.00) + \
489 (abs(eta) <= 2.5) * (pt >= 20) * (-1.28247e-03*pt + 1.3885e-06*pt^2 - 3.56999*log(pt)/pt+0.817471) + \
490 (abs(eta) >= 2.5) * (0.00)}
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 2.5
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.01}
508 # efficiency formula for tau-jets
509 add EfficiencyFormula {15} {0.6}
510}
511
512
513#####################################################
514# Find uniquely identified photons/electrons/tau/jets
515#####################################################
516
517module UniqueObjectFinder UniqueObjectFinder {
518# earlier arrays take precedence over later ones
519# add InputArray InputArray OutputArray
520# add InputArray PhotonIsolation/photons photons
521# add InputArray ElectronIsolation/electrons electrons
522# add InputArray MuonIsolation/muons muons
523 add InputArray JetEnergyScale/jets jets
524}
525
526##################
527# ROOT tree writer
528##################
529
530# tracks, towers and eflow objects are not stored by default in the output.
531# if needed (for jet constituent or other studies), uncomment the relevant
532# "add Branch ..." lines.
533
534# MA5 BEGIN (TO TUNE)
535
536module TreeWriter TreeWriter {
537
538 #add Branch Delphes/allParticles Particle GenParticle
539 #add Branch GenJetFinder/jets GenJet Jet
540 add Branch TrackMerger/tracks Track Track
541 add Branch Calorimeter/towers Tower Tower
542 add Branch Calorimeter/eflowTracks EFlowTrack Track
543 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
544 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
545 add Branch JetEnergyScale/jets JetMA5 Jet
546 add Branch PhotonEfficiency/photons PhotonMA5 Photon
547 add Branch ElectronEfficiency/electrons ElectronMA5 Electron
548 add Branch MuonEfficiency/muons MuonMA5 Muon
549 add Branch MissingET/momentum MissingET MissingET
550}
551
552# MA5 END (TO TUNE)