Fork me on GitHub

source: git/cards/delphes_card_IDEA.tcl@ b1cb322

Last change on this file since b1cb322 was a0db751, checked in by michele <michele.selvaggi@…>, 4 years ago

updated IDEA card with latest TrackCovariance

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