Fork me on GitHub

source: git/cards/delphes_card_CMS_NoFastJet.tcl@ 6cdc544

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 6cdc544 was 8624f58, checked in by Pavel Demin <pavel.demin@…>, 10 years ago

replace Ecal with ECal and Hcal with HCal
  • 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 ECalMinEnergy 0.5
189 set HCalMinEnergy 1.0
190
191 set ECalMinSignificance 1.0
192 set HCalMinSignificance 1.0
193
194 set EFlowTrackOutputArray eflowTracks
195 set EFlowPhotonOutputArray eflowPhotons
196 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
197
198 set pi [expr {acos(-1)}]
199
200 # lists of the edges of each tower in eta and phi
201 # each list starts with the lower edge of the first tower
202 # the list ends with the higher edged of the last tower
203
204 # 5 degrees towers
205 set PhiBins {}
206 for {set i -36} {$i <= 36} {incr i} {
207 add PhiBins [expr {$i * $pi/36.0}]
208 }
209 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} {
210 add EtaPhiBins $eta $PhiBins
211 }
212
213 # 10 degrees towers
214 set PhiBins {}
215 for {set i -18} {$i <= 18} {incr i} {
216 add PhiBins [expr {$i * $pi/18.0}]
217 }
218 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} {
219 add EtaPhiBins $eta $PhiBins
220 }
221
222 # 20 degrees towers
223 set PhiBins {}
224 for {set i -9} {$i <= 9} {incr i} {
225 add PhiBins [expr {$i * $pi/9.0}]
226 }
227 foreach eta {-5 -4.7 -4.525 4.7 5} {
228 add EtaPhiBins $eta $PhiBins
229 }
230
231 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
232 add EnergyFraction {0} {0.0 1.0}
233 # energy fractions for e, gamma and pi0
234 add EnergyFraction {11} {1.0 0.0}
235 add EnergyFraction {22} {1.0 0.0}
236 add EnergyFraction {111} {1.0 0.0}
237 # energy fractions for muon, neutrinos and neutralinos
238 add EnergyFraction {12} {0.0 0.0}
239 add EnergyFraction {13} {0.0 0.0}
240 add EnergyFraction {14} {0.0 0.0}
241 add EnergyFraction {16} {0.0 0.0}
242 add EnergyFraction {1000022} {0.0 0.0}
243 add EnergyFraction {1000023} {0.0 0.0}
244 add EnergyFraction {1000025} {0.0 0.0}
245 add EnergyFraction {1000035} {0.0 0.0}
246 add EnergyFraction {1000045} {0.0 0.0}
247 # energy fractions for K0short and Lambda
248 add EnergyFraction {310} {0.3 0.7}
249 add EnergyFraction {3122} {0.3 0.7}
250
251 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
252 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
253 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
254
255 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
256 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
257 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
258}
259
260####################
261# Energy flow merger
262####################
263
264module Merger EFlowMerger {
265# add InputArray InputArray
266 add InputArray Calorimeter/eflowTracks
267 add InputArray Calorimeter/eflowPhotons
268 add InputArray Calorimeter/eflowNeutralHadrons
269 set OutputArray eflow
270}
271
272
Note: See TracBrowser for help on using the repository browser.