Changeset d4e5b6a in git for cards/FCC/FCChh_PileUpVtx.tcl

Timestamp:

Dec 12, 2019, 3:16:58 PM (5 years ago)

Author:

GitHub <noreply@…>

Branches:

Timing

Children:

03b9c0f, 83ee320, b4a3c55

Parents:

2b5ff2c (diff), 584e2e9 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

git-author:

Michele Selvaggi <michele.selvaggi@…> (12/12/19 15:16:58)

git-committer:

GitHub <noreply@…> (12/12/19 15:16:58)

Message:

Merge pull request #4 from kaanyuxel/master

Addıng various timing related updates

File:

• cards/FCC/FCChh_PileUpVtx.tcl (modified) (10 diffs)

cards/FCC/FCChh_PileUpVtx.tcl

@@ -12 +12 @@
 #######################################
 
+set MaxEvents 100
+set RandomSeed 123
+
 set ExecutionPath {
 
+  BeamSpotFilter
+  PileUpMerger
   ParticlePropagator
 
@@ -25 +30 @@
 
   TrackMerger
+
+  TrackSmearing
+  TimeSmearing  
+
+  VertexFinderDA4D  
+
+  TrackTimingPileUpSubtractor  
 
   ECal
@@ -68 +80 @@
 }
 
+#######################
+# GenBeamSpotFilter
+# Saves a particle intended to represent the beamspot
+#######################
+
+module BeamSpotFilter BeamSpotFilter {
+    set InputArray Delphes/stableParticles
+    set OutputArray beamSpotParticle
+
+}
+
+###############
+# PileUp Merger
+###############
+
+module PileUpMerger PileUpMerger {
+  set InputArray Delphes/stableParticles
+
+  set ParticleOutputArray stableParticles
+  set VertexOutputArray vertices
+
+  # pre-generated minbias input file
+  set PileUpFile MinBias.pileup
+
+  # average expected pile up
+  set MeanPileUp 10
+
+  # 0-poisson, 1-uniform, 2-delta
+  set PileUpDistribution 2
+
+  # maximum spread in the beam direction in m
+  set ZVertexSpread 0.25
+
+  # maximum spread in time in s
+  set TVertexSpread 800E-12
+
+  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
+
+  #set VertexDistributionFormula {exp(-(t^2/(2*(0.063/2.99792458E8*exp(-(z^2/(2*(0.063)^2))))^2)))}
+  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
+
+  # taking 5.3 cm x 160 ps
+
+  #set VertexDistributionFormula { (abs(t) <= 160e-12) * (abs(z) <= 0.053) * (1.00) +
+  #                                (abs(t) >  160e-12) * (abs(z) <= 0.053) * (0.00) +
+  #                               (abs(t) <= 160e-12) * (abs(z) > 0.053)  * (0.00) +
+  #                           (abs(t) >  160e-12) * (abs(z) > 0.053)  * (0.00)}
+
+}
 
 #################################
@@ -74 +135 @@
 
 module ParticlePropagator ParticlePropagator {
-  set InputArray Delphes/stableParticles
+  set InputArray PileUpMerger/stableParticles
   set OutputArray stableParticles
   set ChargedHadronOutputArray chargedHadrons
@@ -204 +265 @@
 
 
-#### ADD TRACKSMEARING
-#### ADD TIMESMEARING
-
+########################################
+#   Smear tracks
+########################################
+
+module TrackSmearing TrackSmearing {
+  set InputArray TrackMerger/tracks
+  set OutputArray tracks
+  set ApplyToPileUp true
+
+  # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
+  source trackResolutionCMS.tcl 
+  # FIXME !!!! we need to add track resolution of FCC-hh baseline detector !!!!!
+}
+
+########################################
+#   Time Smearing
+########################################
+
+module TimeSmearing TimeSmearing {
+  set InputArray TrackSmearing/tracks
+  set OutputArray tracks
+
+  # assume 20 ps resolution for now
+  set TimeResolution 20E-12
+}
+
+##################################
+# Primary vertex reconstruction
+##################################
+
+
+module VertexFinderDA4D VertexFinderDA4D {
+  set InputArray TimeSmearing/tracks
+
+  set OutputArray tracks
+  set VertexOutputArray vertices
+
+  set Verbose 0
+  set MinPT 1.0
+
+  # in mm
+  set VertexSpaceSize 0.5
+
+  # in s
+  set VertexTimeSize 10E-12
+
+  set UseTc 1
+  set BetaMax 0.1
+  set BetaStop 1.0
+  set CoolingFactor 0.8
+  set MaxIterations 100
+
+  # in mm
+  set DzCutOff 40
+  set D0CutOff 30
+
+}
+
+##########################
+# Track pile-up subtractor
+##########################
+
+module TrackTimingPileUpSubtractor TrackTimingPileUpSubtractor {
+# add InputArray InputArray OutputArray
+
+  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
+  add InputArray ElectronMomentumSmearing/electrons
+  add InputArray MuonMomentumSmearing/muons
+  
+  set VertexInputArray VertexFinderDA4D/vertices
+  # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
+  # Z vertex resolution in m
+  set ZVertexResolution {0.0001}
+}
 
 
@@ -217 +349 @@
 module SimpleCalorimeter ECal {
   set ParticleInputArray ParticlePropagator/stableParticles
-  set TrackInputArray TrackMerger/tracks
+  set TrackInputArray TimeSmearing/tracks
 
   set TowerOutputArray ecalTowers
@@ -374 +506 @@
 }
 
+
 #################
 # Electron filter
@@ -494 +627 @@
   set MomentumOutputArray momentum
 }
-
-
 
 
@@ -873 +1004 @@
   add Branch GenMissingET/momentum GenMissingET MissingET
 
-  add Branch TrackMerger/tracks Track Track
+  add Branch TimeSmearing/tracks Track Track
   add Branch Calorimeter/towers Tower Tower
 
@@ -892 +1023 @@
   add Branch MissingET/momentum MissingET MissingET
   add Branch ScalarHT/energy ScalarHT ScalarHT
-}
-
+  add Branch VertexFinderDA4D/vertices Vertex4D Vertex
+}
+