Fork me on GitHub

source: git/cards/delphes_card_IDEAtrkCov.tcl@ b438187

Last change on this file since b438187 was dc883b4, checked in by Pavel Demin <pavel.demin@…>, 4 years ago

add info to TreeWriter

  • Property mode set to 100644
File size: 24.8 KB
Line 
1#################################################################### l
2# FCC-ee IDEA detector model
3#
4# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
5# email: efontane@bo.infn.it,
6# lorenzo.pezzotti01@universitadipavia.it,
7# m.antonello@uninsubria.it,
8# michele.selvaggi@cern.ch
9#####################################################################
10
11set B 2.0
12
13#######################################
14# Order of execution of various modules
15#######################################
16
17set ExecutionPath {
18 ParticlePropagator
19
20 ChargedHadronTrackingEfficiency
21 ElectronTrackingEfficiency
22 MuonTrackingEfficiency
23
24 TrackMergerPre
25 TrackSmearing
26
27 TrackMerger
28 Calorimeter
29 EFlowMerger
30
31 PhotonEfficiency
32 PhotonIsolation
33
34 MuonFilter
35
36 ElectronFilter
37 ElectronEfficiency
38 ElectronIsolation
39
40 MuonEfficiency
41 MuonIsolation
42
43 MissingET
44
45 NeutrinoFilter
46 GenJetFinder
47 GenMissingET
48
49 FastJetFinder
50
51 JetEnergyScale
52
53 JetFlavorAssociation
54
55 BTagging
56 TauTagging
57
58 UniqueObjectFinder
59
60 ScalarHT
61 TreeWriter
62}
63
64#################################
65# Propagate particles in cylinder
66#################################
67
68module ParticlePropagator ParticlePropagator {
69 set InputArray Delphes/stableParticles
70
71 set OutputArray stableParticles
72 set ChargedHadronOutputArray chargedHadrons
73 set ElectronOutputArray electrons
74 set MuonOutputArray muons
75
76 # inner radius of the solenoid, in m
77 set Radius 2.25
78
79 # half-length: z of the solenoid, in m
80 set HalfLength 2.5
81
82 # magnetic field, in T
83 set Bz $B
84}
85
86####################################
87# Charged hadron tracking efficiency
88####################################
89
90module Efficiency ChargedHadronTrackingEfficiency {
91 set InputArray ParticlePropagator/chargedHadrons
92 set OutputArray chargedHadrons
93 # We use only one efficiency, we set only 0 effincency out of eta bounds:
94
95 set EfficiencyFormula {
96 (abs(eta) > 3.0) * (0.000) +
97 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
98 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
99 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
100 }
101}
102
103# (pt <= 0.1) * (0.00) +
104# (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
105# (abs(eta) > 3) * (0.00)
106
107
108
109##############################
110# Electron tracking efficiency
111##############################
112
113module Efficiency ElectronTrackingEfficiency {
114 set InputArray ParticlePropagator/electrons
115 set OutputArray electrons
116
117
118 # Current full simulation with CLICdet provides for electrons:
119 set EfficiencyFormula {
120 (abs(eta) > 3.0) * (0.000) +
121 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
122 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
123 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
124 }
125}
126
127
128##########################
129# Muon tracking efficiency
130##########################
131
132module Efficiency MuonTrackingEfficiency {
133 set InputArray ParticlePropagator/muons
134 set OutputArray muons
135
136 # Current full simulation with CLICdet provides for muons:
137 set EfficiencyFormula {
138 (abs(eta) > 3.0) * (0.000) +
139 (energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
140 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
141 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
142 }
143}
144
145##############
146# Track merger
147##############
148
149module Merger TrackMergerPre {
150# add InputArray InputArray
151 add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
152 add InputArray ElectronTrackingEfficiency/electrons
153 add InputArray MuonTrackingEfficiency/muons
154 set OutputArray tracks
155}
156
157
158########################################
159# Smearing for charged tracks
160########################################
161
162module TrackCovariance TrackSmearing {
163 set InputArray TrackMergerPre/tracks
164 set OutputArray tracks
165
166 ## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt
167 set DetectorGeometry {
168
169 1 PIPE -100 100 0.015 0.0012 0.35276 0 0 0 0 0 0
170 1 VTXLOW -0.12 0.12 0.017 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
171 1 VTXLOW -0.16 0.16 0.023 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
172 1 VTXLOW -0.16 0.16 0.031 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
173 1 VTXHIGH -1 1 0.32 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
174 1 VTXHIGH -1.05 1.05 0.34 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
175 1 DCHCANI -2.125 2.125 0.345 0.0002 0.237223 0 0 0 0 0 0
176 1 DCH -2 2 0.36 0.0147748 1400 1 0.0203738 0 0.0001 0 1
177 1 DCH -2 2 0.374775 0.0147748 1400 1 -0.0212097 0 0.0001 0 1
178 1 DCH -2 2 0.38955 0.0147748 1400 1 0.0220456 0 0.0001 0 1
179 1 DCH -2 2 0.404324 0.0147748 1400 1 -0.0228814 0 0.0001 0 1
180 1 DCH -2 2 0.419099 0.0147748 1400 1 0.0237172 0 0.0001 0 1
181 1 DCH -2 2 0.433874 0.0147748 1400 1 -0.024553 0 0.0001 0 1
182 1 DCH -2 2 0.448649 0.0147748 1400 1 0.0253888 0 0.0001 0 1
183 1 DCH -2 2 0.463423 0.0147748 1400 1 -0.0262245 0 0.0001 0 1
184 1 DCH -2 2 0.478198 0.0147748 1400 1 0.0270602 0 0.0001 0 1
185 1 DCH -2 2 0.492973 0.0147748 1400 1 -0.0278958 0 0.0001 0 1
186 1 DCH -2 2 0.507748 0.0147748 1400 1 0.0287314 0 0.0001 0 1
187 1 DCH -2 2 0.522523 0.0147748 1400 1 -0.029567 0 0.0001 0 1
188 1 DCH -2 2 0.537297 0.0147748 1400 1 0.0304025 0 0.0001 0 1
189 1 DCH -2 2 0.552072 0.0147748 1400 1 -0.031238 0 0.0001 0 1
190 1 DCH -2 2 0.566847 0.0147748 1400 1 0.0320734 0 0.0001 0 1
191 1 DCH -2 2 0.581622 0.0147748 1400 1 -0.0329088 0 0.0001 0 1
192 1 DCH -2 2 0.596396 0.0147748 1400 1 0.0337442 0 0.0001 0 1
193 1 DCH -2 2 0.611171 0.0147748 1400 1 -0.0345795 0 0.0001 0 1
194 1 DCH -2 2 0.625946 0.0147748 1400 1 0.0354147 0 0.0001 0 1
195 1 DCH -2 2 0.640721 0.0147748 1400 1 -0.0362499 0 0.0001 0 1
196 1 DCH -2 2 0.655495 0.0147748 1400 1 0.0370851 0 0.0001 0 1
197 1 DCH -2 2 0.67027 0.0147748 1400 1 -0.0379202 0 0.0001 0 1
198 1 DCH -2 2 0.685045 0.0147748 1400 1 0.0387552 0 0.0001 0 1
199 1 DCH -2 2 0.69982 0.0147748 1400 1 -0.0395902 0 0.0001 0 1
200 1 DCH -2 2 0.714595 0.0147748 1400 1 0.0404252 0 0.0001 0 1
201 1 DCH -2 2 0.729369 0.0147748 1400 1 -0.04126 0 0.0001 0 1
202 1 DCH -2 2 0.744144 0.0147748 1400 1 0.0420949 0 0.0001 0 1
203 1 DCH -2 2 0.758919 0.0147748 1400 1 -0.0429296 0 0.0001 0 1
204 1 DCH -2 2 0.773694 0.0147748 1400 1 0.0437643 0 0.0001 0 1
205 1 DCH -2 2 0.788468 0.0147748 1400 1 -0.044599 0 0.0001 0 1
206 1 DCH -2 2 0.803243 0.0147748 1400 1 0.0454336 0 0.0001 0 1
207 1 DCH -2 2 0.818018 0.0147748 1400 1 -0.0462681 0 0.0001 0 1
208 1 DCH -2 2 0.832793 0.0147748 1400 1 0.0471025 0 0.0001 0 1
209 1 DCH -2 2 0.847568 0.0147748 1400 1 -0.0479369 0 0.0001 0 1
210 1 DCH -2 2 0.862342 0.0147748 1400 1 0.0487713 0 0.0001 0 1
211 1 DCH -2 2 0.877117 0.0147748 1400 1 -0.0496055 0 0.0001 0 1
212 1 DCH -2 2 0.891892 0.0147748 1400 1 0.0504397 0 0.0001 0 1
213 1 DCH -2 2 0.906667 0.0147748 1400 1 -0.0512738 0 0.0001 0 1
214 1 DCH -2 2 0.921441 0.0147748 1400 1 0.0521079 0 0.0001 0 1
215 1 DCH -2 2 0.936216 0.0147748 1400 1 -0.0529418 0 0.0001 0 1
216 1 DCH -2 2 0.950991 0.0147748 1400 1 0.0537757 0 0.0001 0 1
217 1 DCH -2 2 0.965766 0.0147748 1400 1 -0.0546095 0 0.0001 0 1
218 1 DCH -2 2 0.980541 0.0147748 1400 1 0.0554433 0 0.0001 0 1
219 1 DCH -2 2 0.995315 0.0147748 1400 1 -0.056277 0 0.0001 0 1
220 1 DCH -2 2 1.01009 0.0147748 1400 1 0.0571106 0 0.0001 0 1
221 1 DCH -2 2 1.02486 0.0147748 1400 1 -0.0579441 0 0.0001 0 1
222 1 DCH -2 2 1.03964 0.0147748 1400 1 0.0587775 0 0.0001 0 1
223 1 DCH -2 2 1.05441 0.0147748 1400 1 -0.0596108 0 0.0001 0 1
224 1 DCH -2 2 1.06919 0.0147748 1400 1 0.0604441 0 0.0001 0 1
225 1 DCH -2 2 1.08396 0.0147748 1400 1 -0.0612773 0 0.0001 0 1
226 1 DCH -2 2 1.09874 0.0147748 1400 1 0.0621104 0 0.0001 0 1
227 1 DCH -2 2 1.11351 0.0147748 1400 1 -0.0629434 0 0.0001 0 1
228 1 DCH -2 2 1.12829 0.0147748 1400 1 0.0637763 0 0.0001 0 1
229 1 DCH -2 2 1.14306 0.0147748 1400 1 -0.0646092 0 0.0001 0 1
230 1 DCH -2 2 1.15784 0.0147748 1400 1 0.0654419 0 0.0001 0 1
231 1 DCH -2 2 1.17261 0.0147748 1400 1 -0.0662746 0 0.0001 0 1
232 1 DCH -2 2 1.18739 0.0147748 1400 1 0.0671071 0 0.0001 0 1
233 1 DCH -2 2 1.20216 0.0147748 1400 1 -0.0679396 0 0.0001 0 1
234 1 DCH -2 2 1.21694 0.0147748 1400 1 0.068772 0 0.0001 0 1
235 1 DCH -2 2 1.23171 0.0147748 1400 1 -0.0696042 0 0.0001 0 1
236 1 DCH -2 2 1.24649 0.0147748 1400 1 0.0704364 0 0.0001 0 1
237 1 DCH -2 2 1.26126 0.0147748 1400 1 -0.0712685 0 0.0001 0 1
238 1 DCH -2 2 1.27604 0.0147748 1400 1 0.0721005 0 0.0001 0 1
239 1 DCH -2 2 1.29081 0.0147748 1400 1 -0.0729324 0 0.0001 0 1
240 1 DCH -2 2 1.30559 0.0147748 1400 1 0.0737642 0 0.0001 0 1
241 1 DCH -2 2 1.32036 0.0147748 1400 1 -0.0745958 0 0.0001 0 1
242 1 DCH -2 2 1.33514 0.0147748 1400 1 0.0754274 0 0.0001 0 1
243 1 DCH -2 2 1.34991 0.0147748 1400 1 -0.0762589 0 0.0001 0 1
244 1 DCH -2 2 1.36468 0.0147748 1400 1 0.0770903 0 0.0001 0 1
245 1 DCH -2 2 1.37946 0.0147748 1400 1 -0.0779215 0 0.0001 0 1
246 1 DCH -2 2 1.39423 0.0147748 1400 1 0.0787527 0 0.0001 0 1
247 1 DCH -2 2 1.40901 0.0147748 1400 1 -0.0795837 0 0.0001 0 1
248 1 DCH -2 2 1.42378 0.0147748 1400 1 0.0804147 0 0.0001 0 1
249 1 DCH -2 2 1.43856 0.0147748 1400 1 -0.0812455 0 0.0001 0 1
250 1 DCH -2 2 1.45333 0.0147748 1400 1 0.0820762 0 0.0001 0 1
251 1 DCH -2 2 1.46811 0.0147748 1400 1 -0.0829068 0 0.0001 0 1
252 1 DCH -2 2 1.48288 0.0147748 1400 1 0.0837373 0 0.0001 0 1
253 1 DCH -2 2 1.49766 0.0147748 1400 1 -0.0845677 0 0.0001 0 1
254 1 DCH -2 2 1.51243 0.0147748 1400 1 0.0853979 0 0.0001 0 1
255 1 DCH -2 2 1.52721 0.0147748 1400 1 -0.086228 0 0.0001 0 1
256 1 DCH -2 2 1.54198 0.0147748 1400 1 0.087058 0 0.0001 0 1
257 1 DCH -2 2 1.55676 0.0147748 1400 1 -0.0878879 0 0.0001 0 1
258 1 DCH -2 2 1.57153 0.0147748 1400 1 0.0887177 0 0.0001 0 1
259 1 DCH -2 2 1.58631 0.0147748 1400 1 -0.0895474 0 0.0001 0 1
260 1 DCH -2 2 1.60108 0.0147748 1400 1 0.0903769 0 0.0001 0 1
261 1 DCH -2 2 1.61586 0.0147748 1400 1 -0.0912063 0 0.0001 0 1
262 1 DCH -2 2 1.63063 0.0147748 1400 1 0.0920356 0 0.0001 0 1
263 1 DCH -2 2 1.64541 0.0147748 1400 1 -0.0928647 0 0.0001 0 1
264 1 DCH -2 2 1.66018 0.0147748 1400 1 0.0936937 0 0.0001 0 1
265 1 DCH -2 2 1.67495 0.0147748 1400 1 -0.0945226 0 0.0001 0 1
266 1 DCH -2 2 1.68973 0.0147748 1400 1 0.0953514 0 0.0001 0 1
267 1 DCH -2 2 1.7045 0.0147748 1400 1 -0.09618 0 0.0001 0 1
268 1 DCH -2 2 1.71928 0.0147748 1400 1 0.0970085 0 0.0001 0 1
269 1 DCH -2 2 1.73405 0.0147748 1400 1 -0.0978369 0 0.0001 0 1
270 1 DCH -2 2 1.74883 0.0147748 1400 1 0.0986651 0 0.0001 0 1
271 1 DCH -2 2 1.7636 0.0147748 1400 1 -0.0994932 0 0.0001 0 1
272 1 DCH -2 2 1.77838 0.0147748 1400 1 0.100321 0 0.0001 0 1
273 1 DCH -2 2 1.79315 0.0147748 1400 1 -0.101149 0 0.0001 0 1
274 1 DCH -2 2 1.80793 0.0147748 1400 1 0.101977 0 0.0001 0 1
275 1 DCH -2 2 1.8227 0.0147748 1400 1 -0.102804 0 0.0001 0 1
276 1 DCH -2 2 1.83748 0.0147748 1400 1 0.103632 0 0.0001 0 1
277 1 DCH -2 2 1.85225 0.0147748 1400 1 -0.104459 0 0.0001 0 1
278 1 DCH -2 2 1.86703 0.0147748 1400 1 0.105286 0 0.0001 0 1
279 1 DCH -2 2 1.8818 0.0147748 1400 1 -0.106113 0 0.0001 0 1
280 1 DCH -2 2 1.89658 0.0147748 1400 1 0.10694 0 0.0001 0 1
281 1 DCH -2 2 1.91135 0.0147748 1400 1 -0.107766 0 0.0001 0 1
282 1 DCH -2 2 1.92613 0.0147748 1400 1 0.108593 0 0.0001 0 1
283 1 DCH -2 2 1.9409 0.0147748 1400 1 -0.109419 0 0.0001 0 1
284 1 DCH -2 2 1.95568 0.0147748 1400 1 0.110246 0 0.0001 0 1
285 1 DCH -2 2 1.97045 0.0147748 1400 1 -0.111072 0 0.0001 0 1
286 1 DCH -2 2 1.98523 0.0147748 1400 1 0.111898 0 0.0001 0 1
287 1 DCH -2 2 2 0.0147748 1400 1 -0.112723 0 0.0001 0 1
288 1 DCHCANO -2.125 2.125 2.02 0.02 1.667 0 0 0 0 0 0
289 1 BSILWRP -2.35 2.35 2.04 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
290 1 BSILWRP -2.35 2.35 2.06 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
291 1 MAG -2.5 2.5 2.25 0.05 0.0658 0 0 0 0 0 0
292 1 BPRESH -2.55 2.55 2.45 0.02 1 2 0 1.5708 7e-005 0.01 1
293 2 VTXDSK 0.141 0.3 -0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
294 2 VTXDSK 0.138 0.3 -0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
295 2 VTXDSK 0.065 0.3 -0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
296 2 VTXDSK 0.062 0.3 -0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
297 2 VTXDSK 0.062 0.3 0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
298 2 VTXDSK 0.065 0.3 0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
299 2 VTXDSK 0.138 0.3 0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
300 2 VTXDSK 0.141 0.3 0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
301 2 DCHWALL 0.345 2.02 2.125 0.25 5.55 0 0 0 0 0 0
302 2 DCHWALL 0.345 2.02 -2.125 0.25 5.55 0 0 0 0 0 0
303 2 FSILWRP 0.354 2.02 -2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
304 2 FSILWRP 0.35 2.02 -2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
305 2 FSILWRP 0.35 2.02 2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
306 2 FSILWRP 0.354 2.02 2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
307 2 FRAD 0.38 2.09 2.49 0.0043 0.005612 0 0 0 0 0 0
308 2 FRAD 0.38 2.09 -2.49 0.0043 0.005612 0 0 0 0 0 0
309 2 FPRESH 0.39 2.43 -2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
310 2 FPRESH 0.39 2.43 2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
311 }
312
313 set Bz $B
314}
315
316##############
317# Track merger
318##############
319
320module Merger TrackMerger {
321# add InputArray InputArray
322 add InputArray TrackSmearing/tracks
323 set OutputArray tracks
324}
325
326
327#############
328# Calorimeter
329#############
330module DualReadoutCalorimeter Calorimeter {
331 set ParticleInputArray ParticlePropagator/stableParticles
332 set TrackInputArray TrackMerger/tracks
333
334 set TowerOutputArray towers
335 set PhotonOutputArray photons
336
337 set EFlowTrackOutputArray eflowTracks
338 set EFlowPhotonOutputArray eflowPhotons
339 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
340
341 set ECalEnergyMin 0.5
342 set HCalEnergyMin 0.5
343 set EnergyMin 0.5
344 set ECalEnergySignificanceMin 1.0
345 set HCalEnergySignificanceMin 1.0
346 set EnergySignificanceMin 1.0
347
348 set SmearTowerCenter true
349 set pi [expr {acos(-1)}]
350
351 # Lists of the edges of each tower in eta and phi;
352 # each list starts with the lower edge of the first tower;
353 # the list ends with the higher edged of the last tower.
354 # Barrel: deta=0.02 towers up to |eta| <= 0.88 ( up to 45°)
355 # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad)
356 # Cell size: about 6 cm x 6 cm
357
358 #barrel:
359 set PhiBins {}
360 for {set i -120} {$i <= 120} {incr i} {
361 add PhiBins [expr {$i * $pi/120}]
362 }
363 #deta=0.02 units for |eta| <= 0.88
364 for {set i -44} {$i < 45} {incr i} {
365 set eta [expr {$i * 0.02}]
366 add EtaPhiBins $eta $PhiBins
367 }
368
369 #endcaps:
370 set PhiBins {}
371 for {set i -120} {$i <= 120} {incr i} {
372 add PhiBins [expr {$i* $pi/120}]
373 }
374 #deta=0.02 units for 0.88 < |eta| <= 3.0
375 #first, from -3.0 to -0.88
376 for {set i 1} {$i <=106} {incr i} {
377 set eta [expr {-3.00 + $i * 0.02}]
378 add EtaPhiBins $eta $PhiBins
379 }
380 #same for 0.88 to 3.0
381 for {set i 1} {$i <=106} {incr i} {
382 set eta [expr {0.88 + $i * 0.02}]
383 add EtaPhiBins $eta $PhiBins
384 }
385
386 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
387 add EnergyFraction {0} {0.0 1.0}
388 # energy fractions for e, gamma and pi0
389 add EnergyFraction {11} {1.0 0.0}
390 add EnergyFraction {22} {1.0 0.0}
391 add EnergyFraction {111} {1.0 0.0}
392 # energy fractions for muon, neutrinos and neutralinos
393 add EnergyFraction {12} {0.0 0.0}
394 add EnergyFraction {13} {0.0 0.0}
395 add EnergyFraction {14} {0.0 0.0}
396 add EnergyFraction {16} {0.0 0.0}
397 add EnergyFraction {1000022} {0.0 0.0}
398 add EnergyFraction {1000023} {0.0 0.0}
399 add EnergyFraction {1000025} {0.0 0.0}
400 add EnergyFraction {1000035} {0.0 0.0}
401 add EnergyFraction {1000045} {0.0 0.0}
402 # energy fractions for K0short and Lambda
403 add EnergyFraction {310} {0.3 0.7}
404 add EnergyFraction {3122} {0.3 0.7}
405
406
407 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
408 set ECalResolutionFormula {
409 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.11^2)+
410 (abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.11^2)
411 }
412
413 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
414 set HCalResolutionFormula {
415 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.30^2)+
416 (abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.30^2)
417 }
418}
419
420####################
421# Energy flow merger
422####################
423
424module Merger EFlowMerger {
425# add InputArray InputArray
426 add InputArray Calorimeter/eflowTracks
427 add InputArray Calorimeter/eflowPhotons
428 add InputArray Calorimeter/eflowNeutralHadrons
429 set OutputArray eflow
430}
431
432###################
433# Photon efficiency
434###################
435
436module Efficiency PhotonEfficiency {
437 set InputArray Calorimeter/eflowPhotons
438 set OutputArray photons
439
440 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
441 # efficiency formula for photons
442 set EfficiencyFormula {
443 (energy < 2.0) * (0.000)+
444 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
445 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
446 (abs(eta) > 3.0) * (0.000)
447 }
448}
449
450##################
451# Photon isolation
452##################
453
454module Isolation PhotonIsolation {
455 set CandidateInputArray PhotonEfficiency/photons
456 set IsolationInputArray EFlowMerger/eflow
457
458 set OutputArray photons
459
460 set DeltaRMax 0.5
461
462 set PTMin 0.5
463
464 set PTRatioMax 999.
465}
466
467#################
468# Electron filter
469#################
470
471module PdgCodeFilter ElectronFilter {
472 set InputArray Calorimeter/eflowTracks
473 set OutputArray electrons
474 set Invert true
475 add PdgCode {11}
476 add PdgCode {-11}
477}
478
479#################
480# Muon filter
481#################
482
483module PdgCodeFilter MuonFilter {
484 set InputArray Calorimeter/eflowTracks
485 set OutputArray muons
486 set Invert true
487 add PdgCode {13}
488 add PdgCode {-13}
489}
490
491
492#####################
493# Electron efficiency
494#####################
495
496module Efficiency ElectronEfficiency {
497 set InputArray ElectronFilter/electrons
498 set OutputArray electrons
499
500 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
501
502 # efficiency formula for electrons
503 set EfficiencyFormula {
504 (energy < 2.0) * (0.000)+
505 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
506 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
507 (abs(eta) > 3.0) * (0.000)
508 }
509}
510
511####################
512# Electron isolation
513####################
514
515module Isolation ElectronIsolation {
516 set CandidateInputArray ElectronEfficiency/electrons
517 set IsolationInputArray EFlowMerger/eflow
518
519 set OutputArray electrons
520
521 set DeltaRMax 0.5
522
523 set PTMin 0.5
524
525 set PTRatioMax 0.12
526}
527
528#################
529# Muon efficiency
530#################
531
532module Efficiency MuonEfficiency {
533 set InputArray MuonFilter/muons
534 set OutputArray muons
535
536 # set EfficiencyFormula {efficiency as a function of eta and pt}
537
538 # efficiency formula for muons
539 set EfficiencyFormula {
540 (energy < 2.0) * (0.000)+
541 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
542 (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
543 (abs(eta) > 3.0) * (0.000)
544 }
545}
546
547################
548# Muon isolation
549################
550
551module Isolation MuonIsolation {
552 set CandidateInputArray MuonEfficiency/muons
553 set IsolationInputArray EFlowMerger/eflow
554
555 set OutputArray muons
556
557 set DeltaRMax 0.5
558
559 set PTMin 0.5
560
561 set PTRatioMax 0.25
562}
563
564###################
565# Missing ET merger
566###################
567
568module Merger MissingET {
569# add InputArray InputArray
570 add InputArray EFlowMerger/eflow
571 set MomentumOutputArray momentum
572}
573
574##################
575# Scalar HT merger
576##################
577
578module Merger ScalarHT {
579# add InputArray InputArray
580 add InputArray UniqueObjectFinder/jets
581 add InputArray UniqueObjectFinder/electrons
582 add InputArray UniqueObjectFinder/photons
583 add InputArray UniqueObjectFinder/muons
584 set EnergyOutputArray energy
585}
586
587#####################
588# Neutrino Filter
589#####################
590
591module PdgCodeFilter NeutrinoFilter {
592
593 set InputArray Delphes/stableParticles
594 set OutputArray filteredParticles
595
596 set PTMin 0.0
597
598 add PdgCode {12}
599 add PdgCode {14}
600 add PdgCode {16}
601 add PdgCode {-12}
602 add PdgCode {-14}
603 add PdgCode {-16}
604}
605
606
607#####################
608# MC truth jet finder
609#####################
610
611module FastJetFinder GenJetFinder {
612 set InputArray NeutrinoFilter/filteredParticles
613
614 set OutputArray jets
615
616 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
617 set JetAlgorithm 6
618 set ParameterR 0.4
619 set JetPTMin 1.0
620}
621
622
623#########################
624# Gen Missing ET merger
625########################
626
627module Merger GenMissingET {
628# add InputArray InputArray
629 add InputArray NeutrinoFilter/filteredParticles
630 set MomentumOutputArray momentum
631}
632
633############
634# Jet finder
635############
636
637module FastJetFinder FastJetFinder {
638# set InputArray Calorimeter/towers
639 set InputArray EFlowMerger/eflow
640
641 set OutputArray jets
642
643 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
644 set JetAlgorithm 6
645 set ParameterR 0.4
646 set JetPTMin 1.0
647}
648
649##################
650# Jet Energy Scale
651##################
652
653module EnergyScale JetEnergyScale {
654 set InputArray FastJetFinder/jets
655 set OutputArray jets
656
657 # scale formula for jets
658 set ScaleFormula {1.08}
659}
660
661########################
662# Jet Flavor Association
663########################
664
665module JetFlavorAssociation JetFlavorAssociation {
666
667 set PartonInputArray Delphes/partons
668 set ParticleInputArray Delphes/allParticles
669 set ParticleLHEFInputArray Delphes/allParticlesLHEF
670 set JetInputArray JetEnergyScale/jets
671
672 set DeltaR 0.5
673 set PartonPTMin 1.0
674 set PartonEtaMax 3.0
675}
676
677###########
678# b-tagging
679###########
680
681module BTagging BTagging {
682 set JetInputArray JetEnergyScale/jets
683
684 set BitNumber 0
685
686 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
687
688 # default efficiency formula (misidentification rate)
689 add EfficiencyFormula {0} {0.01}
690
691 # efficiency formula for c-jets (misidentification rate)
692 add EfficiencyFormula {4} {0.10}
693
694 # efficiency formula for b-jets
695 add EfficiencyFormula {5} {0.80}
696}
697
698#############
699# tau-tagging
700#############
701
702module TauTagging TauTagging {
703 set ParticleInputArray Delphes/allParticles
704 set PartonInputArray Delphes/partons
705 set JetInputArray JetEnergyScale/jets
706
707 set DeltaR 0.5
708 set TauPTMin 1.0
709 set TauEtaMax 3.0
710
711 # default efficiency formula (misidentification rate)
712 add EfficiencyFormula {0} {0.001}
713 # efficiency formula for tau-jets
714 add EfficiencyFormula {15} {0.6}
715}
716
717
718#####################################################
719# Find uniquely identified photons/electrons/tau/jets
720#####################################################
721
722module UniqueObjectFinder UniqueObjectFinder {
723# earlier arrays take precedence over later ones
724# add InputArray InputArray OutputArray
725 add InputArray PhotonIsolation/photons photons
726 add InputArray ElectronIsolation/electrons electrons
727 add InputArray MuonIsolation/muons muons
728 add InputArray JetEnergyScale/jets jets
729}
730
731
732
733##################
734# ROOT tree writer
735##################
736
737# Tracks, towers and eflow objects are not stored by default in the output.
738# If needed (for jet constituent or other studies), uncomment the relevant
739# "add Branch ..." lines.
740
741module TreeWriter TreeWriter {
742 # add Branch InputArray BranchName BranchClass
743
744 add Branch Delphes/allParticles Particle GenParticle
745
746 add Branch TrackMerger/tracks Track Track
747 add Branch Calorimeter/towers Tower Tower
748
749 add Branch Calorimeter/eflowTracks EFlowTrack Track
750 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
751 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
752
753 add Branch Calorimeter/photons CaloPhoton Photon
754 add Branch PhotonEfficiency/photons PhotonEff Photon
755 add Branch PhotonIsolation/photons PhotonIso Photon
756
757 add Branch GenJetFinder/jets GenJet Jet
758 add Branch GenMissingET/momentum GenMissingET MissingET
759
760 add Branch UniqueObjectFinder/jets Jet Jet
761 add Branch UniqueObjectFinder/electrons Electron Electron
762 add Branch UniqueObjectFinder/photons Photon Photon
763 add Branch UniqueObjectFinder/muons Muon Muon
764
765 add Branch JetEnergyScale/jets AntiKtJet Jet
766
767 add Branch MissingET/momentum MissingET MissingET
768 add Branch ScalarHT/energy ScalarHT ScalarHT
769
770 # add Info InfoName InfoValue
771 add Info Bz $B
772}
Note: See TracBrowser for help on using the repository browser.