Fork me on GitHub

source: git/cards/delphes_card_CMS_NoFastJet.tcl@ 20dc57b

Last change on this file since 20dc57b was a2983ec, checked in by Pavel Demin <pavel.demin@…>, 10 years ago

add parameter SmearTowerCenter to all cards
  • Property mode set to 100644
File size: 11.1 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronEnergySmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17 Calorimeter
18 EFlowMerger
19}
20
21#################################
22# Propagate particles in cylinder
23#################################
24
25module ParticlePropagator ParticlePropagator {
26 set InputArray Delphes/stableParticles
27
28 set OutputArray stableParticles
29 set ChargedHadronOutputArray chargedHadrons
30 set ElectronOutputArray electrons
31 set MuonOutputArray muons
32
33 # radius of the magnetic field coverage, in m
34 set Radius 1.29
35 # half-length of the magnetic field coverage, in m
36 set HalfLength 3.00
37
38 # magnetic field
39 set Bz 3.8
40}
41
42####################################
43# Charged hadron tracking efficiency
44####################################
45
46module Efficiency ChargedHadronTrackingEfficiency {
47 set InputArray ParticlePropagator/chargedHadrons
48 set OutputArray chargedHadrons
49
50 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
51
52 # tracking efficiency formula for charged hadrons
53 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
54 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
55 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
56 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
57 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
58 (abs(eta) > 2.5) * (0.00)}
59}
60
61##############################
62# Electron tracking efficiency
63##############################
64
65module Efficiency ElectronTrackingEfficiency {
66 set InputArray ParticlePropagator/electrons
67 set OutputArray electrons
68
69 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
70
71 # tracking efficiency formula for electrons
72 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
73 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
74 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
75 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
76 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
77 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
78 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
79 (abs(eta) > 2.5) * (0.00)}
80}
81
82##########################
83# Muon tracking efficiency
84##########################
85
86module Efficiency MuonTrackingEfficiency {
87 set InputArray ParticlePropagator/muons
88 set OutputArray muons
89
90 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
91
92 # tracking efficiency formula for muons
93 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
94 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
95 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
96 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
97 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
98 (abs(eta) > 2.5) * (0.00)}
99}
100
101########################################
102# Momentum resolution for charged tracks
103########################################
104
105module MomentumSmearing ChargedHadronMomentumSmearing {
106 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
107 set OutputArray chargedHadrons
108
109 # set ResolutionFormula {resolution formula as a function of eta and pt}
110
111 # resolution formula for charged hadrons
112 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) +
113 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) +
114 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
115 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) +
116 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) +
117 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) +
118 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
119 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
120}
121
122#################################
123# Energy resolution for electrons
124#################################
125
126module EnergySmearing ElectronEnergySmearing {
127 set InputArray ElectronTrackingEfficiency/electrons
128 set OutputArray electrons
129
130 # set ResolutionFormula {resolution formula as a function of eta and energy}
131
132 # resolution formula for electrons
133 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) +
134 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
135 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
136 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
137
138}
139
140###############################
141# Momentum resolution for muons
142###############################
143
144module MomentumSmearing MuonMomentumSmearing {
145 set InputArray MuonTrackingEfficiency/muons
146 set OutputArray muons
147
148 # set ResolutionFormula {resolution formula as a function of eta and pt}
149
150 # resolution formula for muons
151 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) +
152 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) +
153 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
154 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) +
155 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) +
156 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) +
157 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
158 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) +
159 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) +
160 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) +
161 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
162 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
163}
164
165##############
166# Track merger
167##############
168
169module Merger TrackMerger {
170# add InputArray InputArray
171 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
172 add InputArray ElectronEnergySmearing/electrons
173 add InputArray MuonMomentumSmearing/muons
174 set OutputArray tracks
175}
176
177#############
178# Calorimeter
179#############
180
181module Calorimeter Calorimeter {
182 set ParticleInputArray ParticlePropagator/stableParticles
183 set TrackInputArray TrackMerger/tracks
184
185 set TowerOutputArray towers
186 set PhotonOutputArray photons
187
188 set EFlowTrackOutputArray eflowTracks
189 set EFlowPhotonOutputArray eflowPhotons
190 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
191
192 set ECalEnergyMin 0.5
193 set HCalEnergyMin 1.0
194
195 set ECalEnergySignificanceMin 1.0
196 set HCalEnergySignificanceMin 1.0
197
198 set SmearTowerCenter true
199
200 set pi [expr {acos(-1)}]
201
202 # lists of the edges of each tower in eta and phi
203 # each list starts with the lower edge of the first tower
204 # the list ends with the higher edged of the last tower
205
206 # 5 degrees towers
207 set PhiBins {}
208 for {set i -36} {$i <= 36} {incr i} {
209 add PhiBins [expr {$i * $pi/36.0}]
210 }
211 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
212 add EtaPhiBins $eta $PhiBins
213 }
214
215 # 10 degrees towers
216 set PhiBins {}
217 for {set i -18} {$i <= 18} {incr i} {
218 add PhiBins [expr {$i * $pi/18.0}]
219 }
220 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
221 add EtaPhiBins $eta $PhiBins
222 }
223
224 # 20 degrees towers
225 set PhiBins {}
226 for {set i -9} {$i <= 9} {incr i} {
227 add PhiBins [expr {$i * $pi/9.0}]
228 }
229 foreach eta {-5 -4.7 -4.525 4.7 5} {
230 add EtaPhiBins $eta $PhiBins
231 }
232
233 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
234 add EnergyFraction {0} {0.0 1.0}
235 # energy fractions for e, gamma and pi0
236 add EnergyFraction {11} {1.0 0.0}
237 add EnergyFraction {22} {1.0 0.0}
238 add EnergyFraction {111} {1.0 0.0}
239 # energy fractions for muon, neutrinos and neutralinos
240 add EnergyFraction {12} {0.0 0.0}
241 add EnergyFraction {13} {0.0 0.0}
242 add EnergyFraction {14} {0.0 0.0}
243 add EnergyFraction {16} {0.0 0.0}
244 add EnergyFraction {1000022} {0.0 0.0}
245 add EnergyFraction {1000023} {0.0 0.0}
246 add EnergyFraction {1000025} {0.0 0.0}
247 add EnergyFraction {1000035} {0.0 0.0}
248 add EnergyFraction {1000045} {0.0 0.0}
249 # energy fractions for K0short and Lambda
250 add EnergyFraction {310} {0.3 0.7}
251 add EnergyFraction {3122} {0.3 0.7}
252
253 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
254 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
255 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
256
257 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
258 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
259 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
260}
261
262####################
263# Energy flow merger
264####################
265
266module Merger EFlowMerger {
267# add InputArray InputArray
268 add InputArray Calorimeter/eflowTracks
269 add InputArray Calorimeter/eflowPhotons
270 add InputArray Calorimeter/eflowNeutralHadrons
271 set OutputArray eflow
272}
Note: See TracBrowser for help on using the repository browser.