MA5SandBox: delphesMA5tune_card_CMS_sus14004.tcl

File delphesMA5tune_card_CMS_sus14004.tcl, 20.0 KB (added by Benjamin Fuks, 6 years ago)
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 TrackIsolationCalculation
18 Calorimeter
19 EFlowMerger
20
21 PhotonEfficiency
22 PhotonIsolationCalculation
23 PhotonIsolation
24
25 ElectronEfficiency
26 ElectronIsolationCalculation
27 ElectronIsolation
28
29 MuonEfficiency
30 MuonIsolationCalculation
31 MuonIsolation
32
33 MissingET
34
35 GenJetFinder
36 FastJetFinder
37
38 JetEnergyScale
39
40 BTagging
41 TauTagging
42
43 UniqueObjectFinder
44
45 ScalarHT
46
47 TreeWriter
48}
49
50#MA5 BEGIN
51#################################
52# Isolation Calculation
53#################################
54module IsolationCalculation PhotonIsolationCalculation {
55 set CandidateInputArray PhotonEfficiency/photons
56 set IsolationInputArray EFlowMerger/eflow
57
58 set OutputArray DelphesMA5tunePhotons
59 set PTMin 0.5
60}
61
62module IsolationCalculation ElectronIsolationCalculation {
63 set CandidateInputArray ElectronEfficiency/electrons
64 set IsolationInputArray TrackMerger/tracks
65
66 set OutputArray DelphesMA5tuneElectrons
67 set PTMin 0.5
68}
69
70module IsolationCalculation MuonIsolationCalculation {
71 set CandidateInputArray MuonEfficiency/muons
72 set IsolationInputArray TrackMerger/tracks
73
74 set OutputArray DelphesMA5tuneMuons
75 set PTMin 0.5
76}
77module IsolationCalculation TrackIsolationCalculation {
78 set CandidateInputArray TrackMerger/tracks
79 set IsolationInputArray TrackMerger/tracks
80
81 set OutputArray DelphesMA5tuneTracks
82 set PTMin 0.5
83}
84#MA5 END
85
86
87#################################
88# Propagate particles in cylinder
89#################################
90
91module ParticlePropagator ParticlePropagator {
92 set InputArray Delphes/stableParticles
93
94 set OutputArray stableParticles
95 set ChargedHadronOutputArray chargedHadrons
96 set ElectronOutputArray electrons
97 set MuonOutputArray muons
98
99 # radius of the magnetic field coverage, in m
100 set Radius 1.29
101 # half-length of the magnetic field coverage, in m
102 set HalfLength 3.00
103
104 # magnetic field
105 set Bz 3.8
106}
107
108####################################
109# Charged hadron tracking efficiency
110####################################
111
112module Efficiency ChargedHadronTrackingEfficiency {
113 set InputArray ParticlePropagator/chargedHadrons
114 set OutputArray chargedHadrons
115
116 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
117
118 # tracking efficiency formula for charged hadrons
119 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
120 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
121 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
124 (abs(eta) > 2.5) * (0.00)}
125}
126
127##############################
128# Electron tracking efficiency
129##############################
130
131module Efficiency ElectronTrackingEfficiency {
132 set InputArray ParticlePropagator/electrons
133 set OutputArray electrons
134
135 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
136
137 # tracking efficiency formula for electrons
138 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
139 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
140 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
141 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
145 (abs(eta) > 2.5) * (0.00)}
146}
147
148##########################
149# Muon tracking efficiency
150##########################
151
152module Efficiency MuonTrackingEfficiency {
153 set InputArray ParticlePropagator/muons
154 set OutputArray muons
155
156 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
157
158 # tracking efficiency formula for muons
159 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
160 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
161 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
162 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
163 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
164 (abs(eta) > 2.5) * (0.00)}
165}
166
167########################################
168# Momentum resolution for charged tracks
169########################################
170
171module MomentumSmearing ChargedHadronMomentumSmearing {
172 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
173 set OutputArray chargedHadrons
174
175 # set ResolutionFormula {resolution formula as a function of eta and pt}
176
177 # resolution formula for charged hadrons
178 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
179 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
180 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
181 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
183 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
184 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
186}
187
188#################################
189# Energy resolution for electrons
190#################################
191
192module EnergySmearing ElectronEnergySmearing {
193 set InputArray ElectronTrackingEfficiency/electrons
194 set OutputArray electrons
195
196 # set ResolutionFormula {resolution formula as a function of eta and energy}
197
198 # resolution formula for electrons
199 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \
200 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
201 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
202 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
203
204}
205
206###############################
207# Momentum resolution for muons
208###############################
209
210module MomentumSmearing MuonMomentumSmearing {
211 set InputArray MuonTrackingEfficiency/muons
212 set OutputArray muons
213
214 # set ResolutionFormula {resolution formula as a function of eta and pt}
215
216 # resolution formula for muons
217 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
218 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
219 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
220 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
221 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
222 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
223 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
224 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
225 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
226 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
227 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
228 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
229}
230
231##############
232# Track merger
233##############
234
235module Merger TrackMerger {
236# add InputArray InputArray
237 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
238 add InputArray ElectronEnergySmearing/electrons
239 add InputArray MuonMomentumSmearing/muons
240 set OutputArray tracks
241}
242
243#############
244# Calorimeter
245#############
246
247module Calorimeter Calorimeter {
248 set ParticleInputArray ParticlePropagator/stableParticles
249 set TrackInputArray TrackMerger/tracks
250
251 set TowerOutputArray towers
252 set PhotonOutputArray photons
253
254 set EFlowTrackOutputArray eflowTracks
255 set EFlowPhotonOutputArray eflowPhotons
256 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
257
258 set pi [expr {acos(-1)}]
259
260 # lists of the edges of each tower in eta and phi
261 # each list starts with the lower edge of the first tower
262 # the list ends with the higher edged of the last tower
263
264 # 5 degrees towers
265 set PhiBins {}
266 for {set i -36} {$i <= 36} {incr i} {
267 add PhiBins [expr {$i * $pi/36.0}]
268 }
269 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
270 add EtaPhiBins $eta $PhiBins
271 }
272
273 # 10 degrees towers
274 set PhiBins {}
275 for {set i -18} {$i <= 18} {incr i} {
276 add PhiBins [expr {$i * $pi/18.0}]
277 }
278 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
279 add EtaPhiBins $eta $PhiBins
280 }
281
282 # 20 degrees towers
283 set PhiBins {}
284 for {set i -9} {$i <= 9} {incr i} {
285 add PhiBins [expr {$i * $pi/9.0}]
286 }
287 foreach eta {-5 -4.7 -4.525 4.7 5} {
288 add EtaPhiBins $eta $PhiBins
289 }
290
291 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
292 add EnergyFraction {0} {0.0 1.0}
293 # energy fractions for e, gamma and pi0
294 add EnergyFraction {11} {1.0 0.0}
295 add EnergyFraction {22} {1.0 0.0}
296 add EnergyFraction {111} {1.0 0.0}
297 # energy fractions for muon, neutrinos and neutralinos
298 add EnergyFraction {12} {0.0 0.0}
299 add EnergyFraction {13} {0.0 0.0}
300 add EnergyFraction {14} {0.0 0.0}
301 add EnergyFraction {16} {0.0 0.0}
302 add EnergyFraction {1000022} {0.0 0.0}
303 add EnergyFraction {1000023} {0.0 0.0}
304 add EnergyFraction {1000025} {0.0 0.0}
305 add EnergyFraction {1000035} {0.0 0.0}
306 add EnergyFraction {1000045} {0.0 0.0}
307 # energy fractions for K0short and Lambda
308 add EnergyFraction {310} {0.3 0.7}
309 add EnergyFraction {3122} {0.3 0.7}
310
311 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
312 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
313 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
314
315 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
316 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
317 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
318}
319
320####################
321# Energy flow merger
322####################
323
324module Merger EFlowMerger {
325# add InputArray InputArray
326 add InputArray Calorimeter/eflowTracks
327 add InputArray Calorimeter/eflowPhotons
328 add InputArray Calorimeter/eflowNeutralHadrons
329 set OutputArray eflow
330}
331
332###################
333# Photon efficiency
334###################
335
336module Efficiency PhotonEfficiency {
337 set InputArray Calorimeter/eflowPhotons
338 set OutputArray photons
339
340 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
341
342 # efficiency formula for photons
343 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
344 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
345 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
346 (abs(eta) > 2.5) * (0.00)}
347}
348
349##################
350# Photon isolation
351##################
352
353module Isolation PhotonIsolation {
354 set CandidateInputArray PhotonEfficiency/photons
355 set IsolationInputArray EFlowMerger/eflow
356
357 set OutputArray photons
358
359 set DeltaRMax 0.5
360
361 set PTMin 0.5
362
363 set PTRatioMax 0.1
364}
365
366#####################
367# Electron efficiency
368#####################
369
370module Efficiency ElectronEfficiency {
371 set InputArray ElectronEnergySmearing/electrons
372 set OutputArray electrons
373
374 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
375
376 # efficiency formula for electrons
377 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
378 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
379 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
380 (abs(eta) > 2.5) * (0.00)}
381}
382
383####################
384# Electron isolation
385####################
386
387module Isolation ElectronIsolation {
388 set CandidateInputArray ElectronEfficiency/electrons
389 set IsolationInputArray EFlowMerger/eflow
390
391 set OutputArray electrons
392
393 set DeltaRMax 0.5
394
395 set PTMin 0.5
396
397 set PTRatioMax 0.1
398}
399
400#################
401# Muon efficiency
402#################
403
404module Efficiency MuonEfficiency {
405 set InputArray MuonMomentumSmearing/muons
406 set OutputArray muons
407
408 # set EfficiencyFormula {efficiency as a function of eta and pt}
409
410 # efficiency formula for muons
411 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
412 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
413 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
414 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
415 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
416 (abs(eta) > 2.4) * (0.00)}
417}
418
419################
420# Muon isolation
421################
422
423module Isolation MuonIsolation {
424 set CandidateInputArray MuonEfficiency/muons
425 set IsolationInputArray EFlowMerger/eflow
426
427 set OutputArray muons
428
429 set DeltaRMax 0.5
430
431 set PTMin 0.5
432
433 set PTRatioMax 0.1
434}
435
436###################
437# Missing ET merger
438###################
439
440module Merger MissingET {
441# add InputArray InputArray
442 add InputArray EFlowMerger/eflow
443 set MomentumOutputArray momentum
444}
445
446##################
447# Scalar HT merger
448##################
449
450module Merger ScalarHT {
451# add InputArray InputArray
452 add InputArray UniqueObjectFinder/jets
453 add InputArray UniqueObjectFinder/electrons
454 add InputArray UniqueObjectFinder/photons
455 add InputArray UniqueObjectFinder/muons
456 set EnergyOutputArray energy
457}
458
459#####################
460# MC truth jet finder
461#####################
462
463module FastJetFinder GenJetFinder {
464 set InputArray Delphes/stableParticles
465
466 set OutputArray jets
467
468 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
469 set JetAlgorithm 6
470 set ParameterR 0.5
471
472 set JetPTMin 20.0
473}
474
475############
476# Jet finder
477############
478
479module FastJetFinder FastJetFinder {
480# set InputArray Calorimeter/towers
481 set InputArray EFlowMerger/eflow
482
483 set OutputArray jets
484
485 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
486 set JetAlgorithm 6
487 set ParameterR 0.5
488
489 set JetPTMin 20.0
490}
491
492##################
493# Jet Energy Scale
494##################
495
496module EnergyScale JetEnergyScale {
497 set InputArray FastJetFinder/jets
498 set OutputArray jets
499
500 # scale formula for jets
501 set ScaleFormula {1.0}
502}
503
504###########
505# b-tagging
506###########
507
508module BTagging BTagging {
509 set PartonInputArray Delphes/partons
510 set JetInputArray JetEnergyScale/jets
511
512 set BitNumber 0
513
514 set DeltaR 0.5
515
516 set PartonPTMin 1.0
517
518 set PartonEtaMax 2.5
519
520 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
521 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
522 # gluon's PDG code has the lowest priority
523
524 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
525 # default efficiency formula (misidentification rate)
526 #add EfficiencyFormula {0} { 0.0036*pt*exp(-78.22*pt*pt*pt*pt*pt*pt+129.81*pt*pt*pt*pt*pt-44.16*pt*pt*pt*pt-4.13*pt*pt*pt+0.49*pt*pt+1.25*pt+0.59) }
527 add EfficiencyFormula {0} { 0.014 }
528
529 # efficiency formula for c-jets (misidentification rate)
530 add EfficiencyFormula {4} { 0.40*pt*exp(-2.19*pt*pt*pt*pt+2.79*pt*pt+3.01*pt-2.70) }
531
532 # efficiency formula for b-jets
533 #add EfficiencyFormula {5} { 0.828*pow(0.048*pt,2.386)/(1+pow(0.048*pt,2.386)) }
534 add EfficiencyFormula {5} { 0.67 }
535}
536
537module TauTagging TauTagging {
538 set ParticleInputArray Delphes/allParticles
539 set PartonInputArray Delphes/partons
540 set JetInputArray JetEnergyScale/jets
541
542 set DeltaR 0.5
543
544 set TauPTMin 1.0
545
546 set TauEtaMax 2.5
547
548 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
549
550 # default efficiency formula (misidentification rate)
551 add EfficiencyFormula {0} {0.001}
552 # efficiency formula for tau-jets
553 add EfficiencyFormula {15} {0.4}
554}
555
556#####################################################
557# Find uniquely identified photons/electrons/tau/jets
558#####################################################
559
560module UniqueObjectFinder UniqueObjectFinder {
561# earlier arrays take precedence over later ones
562# add InputArray InputArray OutputArray
563 add InputArray PhotonIsolation/photons photons
564 add InputArray ElectronIsolation/electrons electrons
565 add InputArray MuonIsolation/muons muons
566 add InputArray JetEnergyScale/jets jets
567}
568
569##################
570# ROOT tree writer
571##################
572
573# tracks, towers and eflow objects are not stored by default in the output.
574# if needed (for jet constituent or other studies), uncomment the relevant
575# "add Branch ..." lines.
576
577# MA5 BEGIN (TO TUNE)
578
579module TreeWriter TreeWriter {
580# add Branch InputArray BranchName BranchClass
581 add Branch Delphes/allParticles Particle GenParticle
582
583# add Branch TrackMerger/tracks Track Track
584# add Branch Calorimeter/towers Tower Tower
585
586# add Branch Calorimeter/eflowTracks EFlowTrack Track
587# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
588# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
589
590 add Branch GenJetFinder/jets GenJet Jet
591 add Branch UniqueObjectFinder/jets Jet Jet
592# add Branch UniqueObjectFinder/electrons Electron Electron
593# add Branch UniqueObjectFinder/photons Photon Photon
594# add Branch UniqueObjectFinder/muons Muon Muon
595 add Branch MissingET/momentum MissingET MissingET
596 add Branch ScalarHT/energy ScalarHT ScalarHT
597
598 add Branch ElectronIsolationCalculation/DelphesMA5tuneElectrons DelphesMA5tuneElectron Electron
599 add Branch MuonIsolationCalculation/DelphesMA5tuneMuons DelphesMA5tuneMuon Muon
600 add Branch PhotonIsolationCalculation/DelphesMA5tunePhotons DelphesMA5tunePhoton Photon
601 add Branch TrackIsolationCalculation/DelphesMA5tuneTracks DelphesMA5tuneTrack Track
602
603}
604
605# MA5 END (TO TUNE)