Fork me on GitHub

source: svn/trunk/examples/delphes_card_CMS_PileUp_Hector.tcl@ 1362

Last change on this file since 1362 was 1361, checked in by Pavel Demin, 11 years ago

add HectorHit class and corresponding ROOT tree branch
File size: 23.3 KB
Line 
1
2#######################################
3# Order of execution of various modules
4#######################################
5
6set ExecutionPath {
7
8 PileUpMerger
9 ParticlePropagator
10
11 ChargedHadronTrackingEfficiency
12 ElectronTrackingEfficiency
13 MuonTrackingEfficiency
14
15 ChargedHadronMomentumSmearing
16 ElectronEnergySmearing
17 MuonMomentumSmearing
18
19 TrackMerger
20 Calorimeter
21 TrackPileUpSubtractor
22 NeutralTowerMerger
23 EFlowMerger
24
25 GenJetFinder
26
27 Rho
28 FastJetFinder
29 PileUpJetID
30 JetPileUpSubtractor
31
32 JetEnergyScale
33
34 PhotonEfficiency
35 PhotonIsolation
36
37 ElectronEfficiency
38 ElectronIsolation
39
40 MuonEfficiency
41 MuonIsolation
42
43 MissingET
44
45 BTagging
46 TauTagging
47
48 UniqueObjectFinder
49
50 ScalarHT
51
52 Hector420f
53 Hector420b
54
55 TreeWriter
56}
57
58###############
59# PileUp Merger
60###############
61
62module PileUpMerger PileUpMerger {
63 set InputArray Delphes/stableParticles
64
65 set ParticleOutputArray stableParticles
66 set VertexOutputArray vertices
67
68 # pre-generated minbias input file
69 set PileUpFile ../../Delphes/MinBias.pileup
70
71 # average expected pile up
72 set MeanPileUp 10
73
74 # maximum spread in the beam direction in m
75 set ZVertexSpread 0.10
76
77 # maximum spread in time in s
78 set TVertexSpread 1.5E-09
79
80 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
81
82 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
83
84 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
85 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
86 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + \
87 # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
88
89
90}
91
92#################################
93# Propagate particles in cylinder
94#################################
95
96module ParticlePropagator ParticlePropagator {
97 set InputArray PileUpMerger/stableParticles
98
99 set OutputArray stableParticles
100 set ChargedHadronOutputArray chargedHadrons
101 set ElectronOutputArray electrons
102 set MuonOutputArray muons
103
104 # radius of the magnetic field coverage, in m
105 set Radius 1.29
106 # half-length of the magnetic field coverage, in m
107 set HalfLength 3.00
108
109 # magnetic field
110 set Bz 3.8
111}
112
113####################################
114# Charged hadron tracking efficiency
115####################################
116
117module Efficiency ChargedHadronTrackingEfficiency {
118 set InputArray ParticlePropagator/chargedHadrons
119 set OutputArray chargedHadrons
120
121 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
122
123 # tracking efficiency formula for charged hadrons
124 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
125 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
126 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
127 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
128 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
129 (abs(eta) > 2.5) * (0.00)}
130}
131
132##############################
133# Electron tracking efficiency
134##############################
135
136module Efficiency ElectronTrackingEfficiency {
137 set InputArray ParticlePropagator/electrons
138 set OutputArray electrons
139
140 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
141
142 # tracking efficiency formula for electrons
143 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
144 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
145 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
146 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
148 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
150 (abs(eta) > 2.5) * (0.00)}
151}
152
153##########################
154# Muon tracking efficiency
155##########################
156
157module Efficiency MuonTrackingEfficiency {
158 set InputArray ParticlePropagator/muons
159 set OutputArray muons
160
161 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
162
163 # tracking efficiency formula for muons
164 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
165 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
166 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
167 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
168 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
169 (abs(eta) > 2.5) * (0.00)}
170}
171
172########################################
173# Momentum resolution for charged tracks
174########################################
175
176module MomentumSmearing ChargedHadronMomentumSmearing {
177 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
178 set OutputArray chargedHadrons
179
180 # set ResolutionFormula {resolution formula as a function of eta and pt}
181
182 # resolution formula for charged hadrons
183 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
184 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
185 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
186 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
188 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
189 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
190 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
191}
192
193#################################
194# Energy resolution for electrons
195#################################
196
197module EnergySmearing ElectronEnergySmearing {
198 set InputArray ElectronTrackingEfficiency/electrons
199 set OutputArray electrons
200
201 # set ResolutionFormula {resolution formula as a function of eta and energy}
202
203 # resolution formula for electrons
204 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \
205 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
206 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
207 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
208
209}
210
211###############################
212# Momentum resolution for muons
213###############################
214
215module MomentumSmearing MuonMomentumSmearing {
216 set InputArray MuonTrackingEfficiency/muons
217 set OutputArray muons
218
219 # set ResolutionFormula {resolution formula as a function of eta and pt}
220
221 # resolution formula for muons
222 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
223 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
224 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
225 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
226 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
227 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
228 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
229 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
230 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
231 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
232 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
233 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
234}
235
236##############
237# Track merger
238##############
239
240module Merger TrackMerger {
241# add InputArray InputArray
242 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
243 add InputArray ElectronEnergySmearing/electrons
244 add InputArray MuonMomentumSmearing/muons
245 set OutputArray tracks
246}
247
248#############
249# Calorimeter
250#############
251
252module Calorimeter Calorimeter {
253 set ParticleInputArray ParticlePropagator/stableParticles
254 set TrackInputArray TrackMerger/tracks
255
256 set TowerOutputArray towers
257 set PhotonOutputArray photons
258
259 set EFlowTrackOutputArray eflowTracks
260 set EFlowPhotonOutputArray eflowPhotons
261 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
262
263 set pi [expr {acos(-1)}]
264
265 # lists of the edges of each tower in eta and phi
266 # each list starts with the lower edge of the first tower
267 # the list ends with the higher edged of the last tower
268
269 # 5 degrees towers
270 set PhiBins {}
271 for {set i -36} {$i <= 36} {incr i} {
272 add PhiBins [expr {$i * $pi/36.0}]
273 }
274 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} {
275 add EtaPhiBins $eta $PhiBins
276 }
277
278 # 10 degrees towers
279 set PhiBins {}
280 for {set i -18} {$i <= 18} {incr i} {
281 add PhiBins [expr {$i * $pi/18.0}]
282 }
283 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} {
284 add EtaPhiBins $eta $PhiBins
285 }
286
287 # 20 degrees towers
288 set PhiBins {}
289 for {set i -9} {$i <= 9} {incr i} {
290 add PhiBins [expr {$i * $pi/9.0}]
291 }
292 foreach eta {-5 -4.7 -4.525 4.7 5} {
293 add EtaPhiBins $eta $PhiBins
294 }
295
296 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
297 add EnergyFraction {0} {0.0 1.0}
298 # energy fractions for e, gamma and pi0
299 add EnergyFraction {11} {1.0 0.0}
300 add EnergyFraction {22} {1.0 0.0}
301 add EnergyFraction {111} {1.0 0.0}
302 # energy fractions for muon, neutrinos and neutralinos
303 add EnergyFraction {12} {0.0 0.0}
304 add EnergyFraction {13} {0.0 0.0}
305 add EnergyFraction {14} {0.0 0.0}
306 add EnergyFraction {16} {0.0 0.0}
307 add EnergyFraction {1000022} {0.0 0.0}
308 add EnergyFraction {1000023} {0.0 0.0}
309 add EnergyFraction {1000025} {0.0 0.0}
310 add EnergyFraction {1000035} {0.0 0.0}
311 add EnergyFraction {1000045} {0.0 0.0}
312 # energy fractions for K0short and Lambda
313 add EnergyFraction {310} {0.3 0.7}
314 add EnergyFraction {3122} {0.3 0.7}
315
316 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
317 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
318 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
319
320 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
321 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
322 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
323}
324
325##########################
326# Track pile-up subtractor
327##########################
328
329module TrackPileUpSubtractor TrackPileUpSubtractor {
330# add InputArray InputArray OutputArray
331 add InputArray Calorimeter/eflowTracks eflowTracks
332 add InputArray ElectronEnergySmearing/electrons electrons
333 add InputArray MuonMomentumSmearing/muons muons
334
335 set VertexInputArray PileUpMerger/vertices
336 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
337 # Z vertex resolution in m
338 set ZVertexResolution 0.0001
339}
340
341####################
342# Neutral tower merger
343####################
344
345module Merger NeutralTowerMerger {
346# add InputArray InputArray
347 add InputArray Calorimeter/eflowPhotons
348 add InputArray Calorimeter/eflowNeutralHadrons
349 set OutputArray eflowTowers
350}
351
352
353####################
354# Energy flow merger
355####################
356
357module Merger EFlowMerger {
358# add InputArray InputArray
359 add InputArray TrackPileUpSubtractor/eflowTracks
360 add InputArray Calorimeter/eflowPhotons
361 add InputArray Calorimeter/eflowNeutralHadrons
362 set OutputArray eflow
363}
364
365
366#############
367# Rho pile-up
368#############
369
370module FastJetFinder Rho {
371# set InputArray Calorimeter/towers
372 set InputArray EFlowMerger/eflow
373
374 set ComputeRho true
375 set RhoOutputArray rho
376
377 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
378 set AreaAlgorithm 5
379
380 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
381 set JetAlgorithm 4
382 set ParameterR 0.6
383 set GhostEtaMax 5.0
384
385 add RhoEtaRange 0.0 2.5
386 add RhoEtaRange 2.5 5.0
387
388 set JetPTMin 0.0
389}
390
391#####################
392# MC truth jet finder
393#####################
394
395module FastJetFinder GenJetFinder {
396 set InputArray Delphes/stableParticles
397
398 set OutputArray jets
399
400 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
401 set JetAlgorithm 6
402 set ParameterR 0.5
403
404 set JetPTMin 20.0
405}
406
407############
408# Jet finder
409############
410
411module FastJetFinder FastJetFinder {
412# set InputArray Calorimeter/towers
413 set InputArray EFlowMerger/eflow
414
415 set OutputArray jets
416
417 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
418 set AreaAlgorithm 5
419
420 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
421 set JetAlgorithm 6
422 set ParameterR 0.5
423
424 set JetPTMin 20.0
425}
426
427###########################
428# Jet Pile-Up ID
429###########################
430
431module PileUpJetID PileUpJetID {
432 set JetInputArray FastJetFinder/jets
433 set TrackInputArray Calorimeter/eflowTracks
434 set NeutralInputArray NeutralTowerMerger/eflowTowers
435
436 set VertexInputArray PileUpMerger/vertices
437 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
438 # Z vertex resolution in m
439 set ZVertexResolution 0.0001
440
441 set OutputArray jets
442
443 set UseConstituents 0
444 set ParameterR 0.5
445
446 set JetPTMin 20.0
447}
448
449###########################
450# Jet Pile-Up Subtraction
451###########################
452
453module JetPileUpSubtractor JetPileUpSubtractor {
454 set JetInputArray PileUpJetID/jets
455 set RhoInputArray Rho/rho
456
457 set OutputArray jets
458
459 set JetPTMin 20.0
460}
461
462##################
463# Jet Energy Scale
464##################
465
466module EnergyScale JetEnergyScale {
467 set InputArray JetPileUpSubtractor/jets
468 set OutputArray jets
469
470 # scale formula for jets
471 set ScaleFormula {1.0}
472}
473
474###################
475# Photon efficiency
476###################
477
478module Efficiency PhotonEfficiency {
479 set InputArray Calorimeter/eflowPhotons
480 set OutputArray photons
481
482 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
483
484 # efficiency formula for photons
485 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
486 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
487 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
488 (abs(eta) > 2.5) * (0.00)}
489}
490
491
492##################
493# Photon isolation
494##################
495
496module Isolation PhotonIsolation {
497 set CandidateInputArray PhotonEfficiency/photons
498 set IsolationInputArray EFlowMerger/eflow
499 set RhoInputArray Rho/rho
500
501 set OutputArray photons
502
503 set DeltaRMax 0.5
504
505 set PTMin 0.5
506
507 set PTRatioMax 0.1
508}
509
510#####################
511# Electron efficiency
512#####################
513
514module Efficiency ElectronEfficiency {
515 set InputArray TrackPileUpSubtractor/electrons
516 set OutputArray electrons
517
518 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
519
520 # efficiency formula for electrons
521 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
522 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
523 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
524 (abs(eta) > 2.5) * (0.00)}
525}
526
527####################
528# Electron isolation
529####################
530
531module Isolation ElectronIsolation {
532 set CandidateInputArray ElectronEfficiency/electrons
533 set IsolationInputArray EFlowMerger/eflow
534 set RhoInputArray Rho/rho
535
536 set OutputArray electrons
537
538 set DeltaRMax 0.5
539
540 set PTMin 0.5
541
542 set PTRatioMax 0.1
543}
544
545#################
546# Muon efficiency
547#################
548
549module Efficiency MuonEfficiency {
550 set InputArray TrackPileUpSubtractor/muons
551 set OutputArray muons
552
553 # set EfficiencyFormula {efficiency as a function of eta and pt}
554
555 # efficiency formula for muons
556 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
557 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
558 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
559 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
560 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
561 (abs(eta) > 2.4) * (0.00)}
562}
563
564################
565# Muon isolation
566################
567
568module Isolation MuonIsolation {
569 set CandidateInputArray MuonEfficiency/muons
570 set IsolationInputArray EFlowMerger/eflow
571 set RhoInputArray Rho/rho
572
573 set OutputArray muons
574
575 set DeltaRMax 0.5
576
577 set PTMin 0.5
578
579 set PTRatioMax 0.1
580}
581
582###################
583# Missing ET merger
584###################
585
586module Merger MissingET {
587# add InputArray InputArray
588 add InputArray EFlowMerger/eflow
589 set MomentumOutputArray momentum
590}
591
592
593
594##################
595# Scalar HT merger
596##################
597
598module Merger ScalarHT {
599# add InputArray InputArray
600 add InputArray UniqueObjectFinder/jets
601 add InputArray UniqueObjectFinder/electrons
602 add InputArray UniqueObjectFinder/photons
603 add InputArray UniqueObjectFinder/muons
604 set EnergyOutputArray energy
605}
606
607###########
608# b-tagging
609###########
610
611module BTagging BTagging {
612 set PartonInputArray Delphes/partons
613 set JetInputArray JetEnergyScale/jets
614
615 set BitNumber 0
616
617 set DeltaR 0.5
618
619 set PartonPTMin 1.0
620
621 set PartonEtaMax 2.5
622
623 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
624 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
625 # gluon's PDG code has the lowest priority
626
627 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
628 # default efficiency formula (misidentification rate)
629 add EfficiencyFormula {0} {0.001}
630
631 # efficiency formula for c-jets (misidentification rate)
632 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
633 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
634 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
635 (abs(eta) > 2.5) * (0.000)}
636
637 # efficiency formula for b-jets
638 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
639 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
640 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
641 (abs(eta) > 2.5) * (0.000)}
642}
643
644module TauTagging TauTagging {
645 set ParticleInputArray Delphes/allParticles
646 set PartonInputArray Delphes/partons
647 set JetInputArray JetEnergyScale/jets
648
649 set DeltaR 0.5
650
651 set TauPTMin 1.0
652
653 set TauEtaMax 2.5
654
655 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
656
657 # default efficiency formula (misidentification rate)
658 add EfficiencyFormula {0} {0.001}
659 # efficiency formula for tau-jets
660 add EfficiencyFormula {15} {0.4}
661}
662
663#####################################################
664# Find uniquely identified photons/electrons/tau/jets
665#####################################################
666
667module UniqueObjectFinder UniqueObjectFinder {
668# earlier arrays take precedence over later ones
669# add InputArray InputArray OutputArray
670 add InputArray PhotonIsolation/photons photons
671 add InputArray ElectronIsolation/electrons electrons
672 add InputArray MuonIsolation/muons muons
673 add InputArray JetEnergyScale/jets jets
674}
675
676########
677# Hector
678########
679
680module Hector Hector420f {
681 set InputArray ParticlePropagator/stableParticles
682
683 set OutputArray hits
684
685 set Direction 1
686
687 set BeamLineFile examples/LHCB1IR5_5TeV.tfs
688 set IPName IP5
689
690 set BeamLineLength 430.0
691 set Distance 420.0
692
693 set OffsetX -0.97
694 set OffsetS 120
695
696 set SigmaE 0.0
697 set SigmaX 0.0
698 set SigmaY 0.0
699 set SigmaT 0.0
700}
701
702module Hector Hector420b {
703 set InputArray ParticlePropagator/stableParticles
704
705 set OutputArray hits
706
707 set Direction -1
708
709 set BeamLineFile examples/LHCB2IR5_5TeV.tfs
710 set IPName IP5
711
712 set BeamLineLength 430.0
713 set Distance 420.0
714
715 set OffsetX 0.97
716 set OffsetS 120
717
718 set SigmaE 0.0
719 set SigmaX 0.0
720 set SigmaY 0.0
721 set SigmaT 0.0
722}
723
724##################
725# ROOT tree writer
726##################
727
728# tracks, towers and eflow objects are not stored by default in the output.
729# if needed (for jet constituent or other studies), uncomment the relevant
730# "add Branch ..." lines.
731
732
733module TreeWriter TreeWriter {
734# add Branch InputArray BranchName BranchClass
735
736 add Branch Delphes/allParticles Particle GenParticle
737
738# add Branch TrackMerger/tracks Track Track
739# add Branch Calorimeter/towers Tower Tower
740# add Branch Calorimeter/eflowTracks EFlowTrack Track
741# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
742# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
743
744 add Branch GenJetFinder/jets GenJet Jet
745 add Branch UniqueObjectFinder/jets Jet Jet
746 add Branch UniqueObjectFinder/electrons Electron Electron
747 add Branch UniqueObjectFinder/photons Photon Photon
748 add Branch UniqueObjectFinder/muons Muon Muon
749 add Branch MissingET/momentum MissingET MissingET
750 add Branch ScalarHT/energy ScalarHT ScalarHT
751 add Branch Rho/rho Rho Rho
752 add Branch PileUpMerger/vertices Vertex Vertex
753 add Branch Hector420f/hits Hector420f HectorHit
754 add Branch Hector420b/hits Hector420b HectorHit
755}
756
Note: See TracBrowser for help on using the repository browser.