Fork me on GitHub

source: git/cards/CMS_PhaseII/propagatorTest.tcl@ 9117aaa4

Timing
Last change on this file since 9117aaa4 was 6fc566b, checked in by Kaan Yüksel Oyulmaz <kaanyukseloyulmaz@…>, 5 years ago

new card to test 4DVertexing

  • Property mode set to 100644
File size: 7.3 KB
Line 
1set MaxEvents 100
2set RandomSeed 123
3
4
5#
6# Phase II - Pile-Up
7#
8# Main authors: Michele Selvaggi (UCL)
9#
10# Released on:
11#
12# Version: v02 beta - test TrackSmearing and vertexing
13#
14#
15#######################################
16# Order of execution of various modules
17#######################################
18
19set ExecutionPath {
20
21 BeamSpotFilter
22
23 PileUpMerger
24 ParticlePropagator
25
26 ChargedHadronTrackingEfficiency
27 ElectronTrackingEfficiency
28 MuonTrackingEfficiency
29
30 TrackMerger
31 TrackSmearing
32
33 TimeSmearing
34
35 VertexFinderDA4D
36
37 TreeWriter
38}
39
40#######################
41# GenBeamSpotFilter
42# Saves a particle intended to represent the beamspot
43#######################
44
45module BeamSpotFilter BeamSpotFilter {
46 set InputArray Delphes/stableParticles
47 set OutputArray beamSpotParticle
48
49}
50
51###############
52# PileUp Merger
53###############
54
55module PileUpMerger PileUpMerger {
56 set InputArray Delphes/stableParticles
57
58 set ParticleOutputArray stableParticles
59 set VertexOutputArray vertices
60
61 # pre-generated minbias input file
62 set PileUpFile MinBias.pileup
63
64 # average expected pile up
65 set MeanPileUp 5
66
67 # 0-poisson, 1-uniform, 2-delta
68 set PileUpDistribution 2
69
70 # maximum spread in the beam direction in m
71 set ZVertexSpread 0.25
72
73 # maximum spread in time in s
74 set TVertexSpread 800E-12
75
76 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
77
78 #set VertexDistributionFormula {exp(-(t^2/(2*(0.063/2.99792458E8*exp(-(z^2/(2*(0.063)^2))))^2)))}
79 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
80
81 # taking 5.3 cm x 160 ps
82
83 #set VertexDistributionFormula { (abs(t) <= 160e-12) * (abs(z) <= 0.053) * (1.00) +
84 # (abs(t) > 160e-12) * (abs(z) <= 0.053) * (0.00) +
85 # (abs(t) <= 160e-12) * (abs(z) > 0.053) * (0.00) +
86 # (abs(t) > 160e-12) * (abs(z) > 0.053) * (0.00)}
87
88}
89
90#################################
91# Propagate particles in cylinder
92#################################
93
94module ParticlePropagator ParticlePropagator {
95 set InputArray PileUpMerger/stableParticles
96
97 set OutputArray stableParticles
98 set ChargedHadronOutputArray chargedHadrons
99 set ElectronOutputArray electrons
100 set MuonOutputArray muons
101
102 # radius of the magnetic field coverage, in m
103 set Radius 1.29
104 # half-length of the magnetic field coverage, in m
105 set HalfLength 3.0
106
107 # magnetic field
108 set Bz 3.8
109}
110
111####################################
112# Charged hadron tracking efficiency
113####################################
114
115module Efficiency ChargedHadronTrackingEfficiency {
116 set InputArray ParticlePropagator/chargedHadrons
117 set OutputArray chargedHadrons
118
119 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
120
121 # tracking efficiency formula for charged hadrons
122 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
123 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
124 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
125 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
126 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
127 (abs(eta) > 2.5) * (0.00)}
128}
129
130##############################
131# Electron tracking efficiency
132##############################
133
134module Efficiency ElectronTrackingEfficiency {
135 set InputArray ParticlePropagator/electrons
136 set OutputArray electrons
137
138 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
139
140 # tracking efficiency formula for electrons
141 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
142 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
143 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
144 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
148 (abs(eta) > 2.5) * (0.00)}
149}
150
151##########################
152# Muon tracking efficiency
153##########################
154
155module Efficiency MuonTrackingEfficiency {
156 set InputArray ParticlePropagator/muons
157 set OutputArray muons
158
159 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
160
161 # tracking efficiency formula for muons
162 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
163 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
164 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
165 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
166
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 && pt <= 1.0e3) * (0.98) +
169 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
170 (abs(eta) > 2.5) * (0.00)}
171}
172
173##############
174# Track merger
175##############
176
177module Merger TrackMerger {
178# add InputArray InputArray
179 add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
180 add InputArray ElectronTrackingEfficiency/electrons
181 add InputArray MuonTrackingEfficiency/muons
182 set OutputArray tracks
183}
184
185########################################
186# Smear tracks
187########################################
188
189module TrackSmearing TrackSmearing {
190 set InputArray TrackMerger/tracks
191 set OutputArray tracks
192 set ApplyToPileUp true
193
194 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
195 source trackResolutionCMS.tcl
196}
197
198module TimeSmearing TimeSmearing {
199 set InputArray TrackSmearing/tracks
200 set OutputArray tracks
201
202 # assume 20 ps resolution for now
203 set TimeResolution 20E-12
204}
205
206##################################
207# Primary vertex reconstruction
208##################################
209
210
211module VertexFinderDA4D VertexFinderDA4D {
212 set InputArray TimeSmearing/tracks
213
214 set OutputArray tracks
215 set VertexOutputArray vertices
216
217 set Verbose 0
218 set MinPT 1.0
219
220 # in mm
221 set VertexSpaceSize 0.5
222
223 # in s
224 set VertexTimeSize 10E-12
225
226 set UseTc 1
227 set BetaMax 0.1
228 set BetaStop 1.0
229 set CoolingFactor 0.8
230 set MaxIterations 100
231
232 # in mm
233 set DzCutOff 40
234 set D0CutOff 30
235
236}
237
238##################
239# ROOT tree writer
240##################
241
242module TreeWriter TreeWriter {
243# add Branch InputArray BranchName BranchClass
244 add Branch PileUpMerger/stableParticles Particle GenParticle
245 add Branch TimeSmearing/tracks Track Track
246 add Branch VertexFinderDA4D/vertices Vertex4D Vertex
247 add Branch PileUpMerger/vertices GenVertex Vertex
248}
Note: See TracBrowser for help on using the repository browser.