Fork me on GitHub

source: git/cards/delphes_card_IDEA.tcl@ d489660

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

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