Fork me on GitHub

source: git/cards/delphes_card_IDEA.tcl@ a3784f9

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

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