Changeset 34c28a9 in git for cards/FCC/FCChh_PileUpVtx.tcl

Timestamp:

Dec 9, 2019, 11:07:13 PM (5 years ago)

Author:

Kaan Yüksel Oyulmaz <kaanyukseloyulmaz@…>

Branches:

Timing

Children:

9117aaa4

Parents:

1363bf5

Message:

FCC-hh card with 4DVertex Algorithm

File:

• cards/FCC/FCChh_PileUpVtx.tcl (modified) (8 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  
+
+  TrackPileUpSubtractor  
 
   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)}
+
+}
+
 #####################################
 # Track propagation to calorimeters
@@ -73 +135 @@
 
 module ParticlePropagator ParticlePropagator {
-  set InputArray Delphes/stableParticles
+  set InputArray PileUpMerger/stableParticles
   set OutputArray stableParticles
   set ChargedHadronOutputArray chargedHadrons
@@ -203 +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 TrackPileUpSubtractor TrackPileUpSubtractor {
+# 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}
+}
 
 
@@ -216 +349 @@
 module SimpleCalorimeter ECal {
   set ParticleInputArray ParticlePropagator/stableParticles
-  set TrackInputArray TrackMerger/tracks
+  set TrackInputArray TrackSmearing/tracks
 
   set TowerOutputArray ecalTowers
@@ -373 +506 @@
 }
 
+
 #################
 # Electron filter
@@ -891 +1025 @@
   add Branch MissingET/momentum MissingET MissingET
   add Branch ScalarHT/energy ScalarHT ScalarHT
-}
-
+  add Branch VertexFinderDA4D/vertices Vertex4D Vertex
+}
+