Fork me on GitHub

source: git/cards/delphes_card_IDEA.tcl@ 1c1c9c2

Last change on this file since 1c1c9c2 was 1c1c9c2, checked in by Michele Selvaggi <michele.selvaggi@…>, 3 years ago

updated IDEA card: genkT jets, PFcands in tree
  • Property mode set to 100644
File size: 28.5 KB
RevLine 
[3b02069]1#################################################################### l
[a0db751]2# FCC-ee IDEA detector model
3#
[b8a6aa3]4# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
[7c9f3da]5# email: efontane@bo.infn.it,
[a0db751]6# lorenzo.pezzotti01@universitadipavia.it,
7# m.antonello@uninsubria.it,
[b8a6aa3]8# michele.selvaggi@cern.ch
[a0db751]9#####################################################################
[b8a6aa3]10
11set B 2.0
12
[d489660]13## Drift chamber coordinates
14set DCHZMIN -2.125
15set DCHZMAX 2.125
16set DCHRMIN 0.345
17set DCHRMAX 2.02
18
19
[ef8a06d]20#######################################
21# Order of execution of various modules
22#######################################
23
24set ExecutionPath {
[058667b]25
26 TruthVertexFinder
[ef8a06d]27 ParticlePropagator
28
29 ChargedHadronTrackingEfficiency
30 ElectronTrackingEfficiency
31 MuonTrackingEfficiency
32
[a0db751]33 TrackMergerPre
34 TrackSmearing
[a95da74]35 ClusterCounting
[058667b]36 TimeSmearing
37 TimeOfFlight
[ef8a06d]38
[a0db751]39 TrackMerger
[ef8a06d]40 Calorimeter
41 EFlowMerger
42
43 PhotonEfficiency
44 PhotonIsolation
45
[b8a6aa3]46 MuonFilter
47
[ef8a06d]48 ElectronFilter
49 ElectronEfficiency
50 ElectronIsolation
51
52 MuonEfficiency
53 MuonIsolation
54
55 MissingET
56
57 NeutrinoFilter
58 GenJetFinder
59 GenMissingET
[a0db751]60
[ef8a06d]61 FastJetFinder
62
63 JetEnergyScale
64
65 JetFlavorAssociation
66
67 BTagging
68 TauTagging
69
70 UniqueObjectFinder
71
72 ScalarHT
73 TreeWriter
74}
75
[058667b]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
[ef8a06d]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
[7c9f3da]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
[ef8a06d]106
[7c9f3da]107 # magnetic field, in T
[b8a6aa3]108 set Bz $B
[ef8a06d]109}
110
111####################################
112# Charged hadron tracking efficiency
113####################################
114
115module Efficiency ChargedHadronTrackingEfficiency {
[7c9f3da]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 }
[ef8a06d]126}
127
[7c9f3da]128
129
[ef8a06d]130##############################
131# Electron tracking efficiency
132##############################
133
134module Efficiency ElectronTrackingEfficiency {
[7c9f3da]135 set InputArray ParticlePropagator/electrons
136 set OutputArray electrons
[ef8a06d]137
138
[7c9f3da]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)
[a0db751]145 }
[ef8a06d]146}
147
[7c9f3da]148
[ef8a06d]149##########################
150# Muon tracking efficiency
151##########################
152
153module Efficiency MuonTrackingEfficiency {
[7c9f3da]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 }
[ef8a06d]164}
165
[a0db751]166##############
167# Track merger
168##############
[7c9f3da]169
[a0db751]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}
[ef8a06d]177
178
[a0db751]179########################################
180# Smearing for charged tracks
181########################################
[ef8a06d]182
[a0db751]183module TrackCovariance TrackSmearing {
[b8a6aa3]184
185 set InputArray TrackMergerPre/tracks
186 set OutputArray tracks
187
[a0db751]188 ## minimum number of hits to accept a track
189 set NMinHits 6
190
[a95da74]191 ## magnetic field
192 set Bz $B
193
[a0db751]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
[b8a6aa3]212 # barrel name zmin zmax r w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
[a0db751]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
[b8a6aa3]216 1 PIPE -100 100 0.015 0.001655 0.2805 0 0 0 0 0 0
[a0db751]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
[d489660]234 1 DCHCANI $DCHZMIN $DCHZMAX $DCHRMIN 0.0002 0.237223 0 0 0 0 0 0
[a0db751]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
[d489660]347 1 DCHCANO $DCHZMIN $DCHZMAX $DCHRMAX $DCHRMAX 0.02 1.667 0 0 0 0 0 0
[a0db751]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
[d489660]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
[a0db751]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
[7c9f3da]362 }
[ef8a06d]363
364}
365
[a95da74]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
[058667b]377 ## check that these are consistent with DCHCANI/DCHNANO parameters in TrackCovariance module
[3fdfd04]378 set Rmin $DCHRMIN
379 set Rmax $DCHRMAX
380 set Zmin $DCHZMIN
381 set Zmax $DCHZMAX
[a95da74]382
[058667b]383 # gas mix option:
[a95da74]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
[058667b]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
[a95da74]425
[ef8a06d]426##############
427# Track merger
428##############
429
430module Merger TrackMerger {
431# add InputArray InputArray
[058667b]432 add InputArray TimeOfFlight/tracks
[ef8a06d]433 set OutputArray tracks
434}
435
436
[a0db751]437#############
438# Calorimeter
439#############
[ef8a06d]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
[7c9f3da]452 set HCalEnergyMin 0.5
[a1b19ea]453 set EnergyMin 0.5
[ef8a06d]454 set ECalEnergySignificanceMin 1.0
455 set HCalEnergySignificanceMin 1.0
456 set EnergySignificanceMin 1.0
457
458 set SmearTowerCenter true
[7c9f3da]459 set pi [expr {acos(-1)}]
460
[a0db751]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
[7c9f3da]467
[a0db751]468 #barrel:
[7c9f3da]469 set PhiBins {}
470 for {set i -120} {$i <= 120} {incr i} {
471 add PhiBins [expr {$i * $pi/120}]
472 }
[a0db751]473 #deta=0.02 units for |eta| <= 0.88
[7c9f3da]474 for {set i -44} {$i < 45} {incr i} {
475 set eta [expr {$i * 0.02}]
476 add EtaPhiBins $eta $PhiBins
477 }
478
[a0db751]479 #endcaps:
[7c9f3da]480 set PhiBins {}
481 for {set i -120} {$i <= 120} {incr i} {
482 add PhiBins [expr {$i* $pi/120}]
483 }
[a0db751]484 #deta=0.02 units for 0.88 < |eta| <= 3.0
485 #first, from -3.0 to -0.88
[7c9f3da]486 for {set i 1} {$i <=106} {incr i} {
487 set eta [expr {-3.00 + $i * 0.02}]
488 add EtaPhiBins $eta $PhiBins
489 }
[a0db751]490 #same for 0.88 to 3.0
[7c9f3da]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
[a0db751]496 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
[7c9f3da]497 add EnergyFraction {0} {0.0 1.0}
[a0db751]498 # energy fractions for e, gamma and pi0
[7c9f3da]499 add EnergyFraction {11} {1.0 0.0}
500 add EnergyFraction {22} {1.0 0.0}
501 add EnergyFraction {111} {1.0 0.0}
[a0db751]502 # energy fractions for muon, neutrinos and neutralinos
[7c9f3da]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}
[a0db751]512 # energy fractions for K0short and Lambda
[7c9f3da]513 add EnergyFraction {310} {0.3 0.7}
514 add EnergyFraction {3122} {0.3 0.7}
515
516
[a0db751]517 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
[7c9f3da]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
[a0db751]523 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
[7c9f3da]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 }
[ef8a06d]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
[7c9f3da]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 }
[ef8a06d]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
[7c9f3da]574 set PTRatioMax 999.
[ef8a06d]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
[a0db751]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
[ef8a06d]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
[a0db751]613 set EfficiencyFormula {
[7c9f3da]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 }
[ef8a06d]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 {
[a0db751]643 set InputArray MuonFilter/muons
[ef8a06d]644 set OutputArray muons
645
646 # set EfficiencyFormula {efficiency as a function of eta and pt}
647
648 # efficiency formula for muons
[a0db751]649 set EfficiencyFormula {
[7c9f3da]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 }
[ef8a06d]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 set OutputArray jets
724
[1c1c9c2]725 set JetAlgorithm 10
726 set ParameterR 1.5
727 set ParameterP -1.0
[7c9f3da]728 set JetPTMin 1.0
[1c1c9c2]729
[ef8a06d]730}
731
[7c9f3da]732
[ef8a06d]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
[1c1c9c2]754 set JetAlgorithm 10
755 set ParameterR 1.5
756 set ParameterP -1.0
[7c9f3da]757 set JetPTMin 1.0
[1c1c9c2]758
759
[ef8a06d]760}
761
762##################
763# Jet Energy Scale
764##################
765
766module EnergyScale JetEnergyScale {
767 set InputArray FastJetFinder/jets
768 set OutputArray jets
769
770 # scale formula for jets
[1c1c9c2]771 set ScaleFormula {1.00}
[ef8a06d]772}
773
774########################
775# Jet Flavor Association
776########################
777
778module JetFlavorAssociation JetFlavorAssociation {
779
780 set PartonInputArray Delphes/partons
781 set ParticleInputArray Delphes/allParticles
782 set ParticleLHEFInputArray Delphes/allParticlesLHEF
783 set JetInputArray JetEnergyScale/jets
784
785 set DeltaR 0.5
786 set PartonPTMin 1.0
[7c9f3da]787 set PartonEtaMax 3.0
[ef8a06d]788}
789
790###########
791# b-tagging
792###########
793
794module BTagging BTagging {
795 set JetInputArray JetEnergyScale/jets
796
797 set BitNumber 0
798
799 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
[a0db751]800
[ef8a06d]801 # default efficiency formula (misidentification rate)
[7c9f3da]802 add EfficiencyFormula {0} {0.01}
[ef8a06d]803
804 # efficiency formula for c-jets (misidentification rate)
[7c9f3da]805 add EfficiencyFormula {4} {0.10}
[ef8a06d]806
807 # efficiency formula for b-jets
[7c9f3da]808 add EfficiencyFormula {5} {0.80}
[ef8a06d]809}
810
811#############
812# tau-tagging
813#############
814
815module TauTagging TauTagging {
816 set ParticleInputArray Delphes/allParticles
817 set PartonInputArray Delphes/partons
818 set JetInputArray JetEnergyScale/jets
819
820 set DeltaR 0.5
821 set TauPTMin 1.0
[7c9f3da]822 set TauEtaMax 3.0
[ef8a06d]823
824 # default efficiency formula (misidentification rate)
[7c9f3da]825 add EfficiencyFormula {0} {0.001}
[ef8a06d]826 # efficiency formula for tau-jets
827 add EfficiencyFormula {15} {0.6}
828}
829
[7c9f3da]830
[ef8a06d]831#####################################################
832# Find uniquely identified photons/electrons/tau/jets
833#####################################################
834
835module UniqueObjectFinder UniqueObjectFinder {
[3b02069]836
837 set UseUniqueID true
[ef8a06d]838# earlier arrays take precedence over later ones
839# add InputArray InputArray OutputArray
840 add InputArray PhotonIsolation/photons photons
841 add InputArray ElectronIsolation/electrons electrons
842 add InputArray MuonIsolation/muons muons
843 add InputArray JetEnergyScale/jets jets
844}
845
[7c9f3da]846
[ef8a06d]847##################
848# ROOT tree writer
849##################
850
[7c9f3da]851# Tracks, towers and eflow objects are not stored by default in the output.
852# If needed (for jet constituent or other studies), uncomment the relevant
[ef8a06d]853# "add Branch ..." lines.
854
855module TreeWriter TreeWriter {
[7c9f3da]856 # add Branch InputArray BranchName BranchClass
[a0db751]857
[7c9f3da]858 add Branch Delphes/allParticles Particle GenParticle
859
860 add Branch TrackMerger/tracks Track Track
861 add Branch Calorimeter/towers Tower Tower
[a0db751]862
[7c9f3da]863 add Branch Calorimeter/eflowTracks EFlowTrack Track
864 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
865 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
866
[1c1c9c2]867 add Branch EFlowMerger/eflow ParticleFlowCandidate ParticleFlowCandidate
868
[7c9f3da]869 add Branch Calorimeter/photons CaloPhoton Photon
870 add Branch PhotonEfficiency/photons PhotonEff Photon
871 add Branch PhotonIsolation/photons PhotonIso Photon
[a0db751]872
[7c9f3da]873 add Branch GenJetFinder/jets GenJet Jet
874 add Branch GenMissingET/momentum GenMissingET MissingET
[a0db751]875
[7c9f3da]876 add Branch UniqueObjectFinder/jets Jet Jet
877 add Branch UniqueObjectFinder/electrons Electron Electron
878 add Branch UniqueObjectFinder/photons Photon Photon
879 add Branch UniqueObjectFinder/muons Muon Muon
[a0db751]880
[7c9f3da]881 add Branch MissingET/momentum MissingET MissingET
882 add Branch ScalarHT/energy ScalarHT ScalarHT
[b8a6aa3]883
884 # add Info InfoName InfoValue
885 add Info Bz $B
[ef8a06d]886}
Note: See TracBrowser for help on using the repository browser.