Changes in cards/FCC/FCChh_PileUpVtx.tcl [542eb90:2c81caa] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/FCC/FCChh_PileUpVtx.tcl
r542eb90 r2c81caa 12 12 ####################################### 13 13 14 set MaxEvents 10015 set RandomSeed 12316 17 14 set ExecutionPath { 18 15 19 BeamSpotFilter20 PileUpMerger21 16 ParticlePropagator 22 17 … … 30 25 31 26 TrackMerger 32 33 TrackSmearing34 TimeSmearing35 36 VertexFinderDA4D37 38 TrackTimingPileUpSubtractor39 27 40 28 ECal … … 80 68 } 81 69 82 ####################### 83 # GenBeamSpotFilter 84 # Saves a particle intended to represent the beamspot 85 ####################### 86 87 module BeamSpotFilter BeamSpotFilter { 88 set InputArray Delphes/stableParticles 89 set OutputArray beamSpotParticle 90 91 } 92 93 ############### 94 # PileUp Merger 95 ############### 96 97 module PileUpMerger PileUpMerger { 70 71 ################################# 72 # Propagate particles in cylinder 73 ################################# 74 75 module ParticlePropagator ParticlePropagator { 98 76 set InputArray Delphes/stableParticles 99 100 set ParticleOutputArray stableParticles101 set VertexOutputArray vertices102 103 # pre-generated minbias input file104 set PileUpFile MinBias.pileup105 106 # average expected pile up107 set MeanPileUp 10108 109 # 0-poisson, 1-uniform, 2-delta110 set PileUpDistribution 2111 112 # maximum spread in the beam direction in m113 set ZVertexSpread 0.25114 115 # maximum spread in time in s116 set TVertexSpread 800E-12117 118 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)119 120 #set VertexDistributionFormula {exp(-(t^2/(2*(0.063/2.99792458E8*exp(-(z^2/(2*(0.063)^2))))^2)))}121 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}122 123 # taking 5.3 cm x 160 ps124 125 #set VertexDistributionFormula { (abs(t) <= 160e-12) * (abs(z) <= 0.053) * (1.00) +126 # (abs(t) > 160e-12) * (abs(z) <= 0.053) * (0.00) +127 # (abs(t) <= 160e-12) * (abs(z) > 0.053) * (0.00) +128 # (abs(t) > 160e-12) * (abs(z) > 0.053) * (0.00)}129 130 }131 132 #####################################133 # Track propagation to calorimeters134 #####################################135 136 module ParticlePropagator ParticlePropagator {137 set InputArray PileUpMerger/stableParticles138 77 set OutputArray stableParticles 139 78 set ChargedHadronOutputArray chargedHadrons … … 265 204 266 205 267 ######################################## 268 # Smear tracks 269 ######################################## 270 271 module TrackSmearing TrackSmearing { 272 set InputArray TrackMerger/tracks 273 set OutputArray tracks 274 set ApplyToPileUp true 275 276 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html 277 source trackResolutionCMS.tcl 278 # FIXME !!!! we need to add track resolution of FCC-hh baseline detector !!!!! 279 } 280 281 ######################################## 282 # Time Smearing 283 ######################################## 284 285 module TimeSmearing TimeSmearing { 286 set InputArray TrackSmearing/tracks 287 set OutputArray tracks 288 289 # assume 20 ps resolution for now 290 set TimeResolution 20E-12 291 } 292 293 ################################## 294 # Primary vertex reconstruction 295 ################################## 296 297 298 module VertexFinderDA4D VertexFinderDA4D { 299 set InputArray TimeSmearing/tracks 300 301 set OutputArray tracks 302 set VertexOutputArray vertices 303 304 set Verbose 0 305 set MinPT 1.0 306 307 # in mm 308 set VertexSpaceSize 0.5 309 310 # in s 311 set VertexTimeSize 10E-12 312 313 set UseTc 1 314 set BetaMax 0.1 315 set BetaStop 1.0 316 set CoolingFactor 0.8 317 set MaxIterations 100 318 319 # in mm 320 set DzCutOff 40 321 set D0CutOff 30 322 323 } 324 325 ########################## 326 # Track pile-up subtractor 327 ########################## 328 329 module TrackTimingPileUpSubtractor TrackTimingPileUpSubtractor { 330 # add InputArray InputArray OutputArray 331 332 add InputArray ChargedHadronMomentumSmearing/chargedHadrons 333 add InputArray ElectronMomentumSmearing/electrons 334 add InputArray MuonMomentumSmearing/muons 335 336 set VertexInputArray VertexFinderDA4D/vertices 337 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution 338 # Z vertex resolution in m 339 set ZVertexResolution {0.0001} 340 } 206 #### ADD TRACKSMEARING 207 #### ADD TIMESMEARING 208 341 209 342 210 … … 349 217 module SimpleCalorimeter ECal { 350 218 set ParticleInputArray ParticlePropagator/stableParticles 351 set TrackInputArray T imeSmearing/tracks219 set TrackInputArray TrackMerger/tracks 352 220 353 221 set TowerOutputArray ecalTowers … … 506 374 } 507 375 508 509 376 ################# 510 377 # Electron filter … … 627 494 set MomentumOutputArray momentum 628 495 } 496 497 629 498 630 499 … … 1004 873 add Branch GenMissingET/momentum GenMissingET MissingET 1005 874 1006 add Branch T imeSmearing/tracks Track Track875 add Branch TrackMerger/tracks Track Track 1007 876 add Branch Calorimeter/towers Tower Tower 1008 877 … … 1023 892 add Branch MissingET/momentum MissingET MissingET 1024 893 add Branch ScalarHT/energy ScalarHT ScalarHT 1025 add Branch VertexFinderDA4D/vertices Vertex4D Vertex 1026 } 1027 894 } 895
Note:
See TracChangeset
for help on using the changeset viewer.