Changes in cards/FCC/FCChh_PileUpVtx.tcl [542eb90:2c81caa] in git

File:

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

cards/FCC/FCChh_PileUpVtx.tcl

@@ -12 +12 @@
 #######################################
 
-set MaxEvents 100
-set RandomSeed 123
-
 set ExecutionPath {
 
-  BeamSpotFilter
-  PileUpMerger
   ParticlePropagator
 
@@ -30 +25 @@
 
   TrackMerger
-
-  TrackSmearing
-  TimeSmearing  
-
-  VertexFinderDA4D  
-
-  TrackTimingPileUpSubtractor  
 
   ECal
@@ -80 +68 @@
 }
 
-#######################
-# GenBeamSpotFilter
-# Saves a particle intended to represent the beamspot
-#######################
-
-module BeamSpotFilter BeamSpotFilter {
-    set InputArray Delphes/stableParticles
-    set OutputArray beamSpotParticle
-
-}
-
-###############
-# PileUp Merger
-###############
-
-module PileUpMerger PileUpMerger {
+
+#################################
+# Propagate particles in cylinder
+#################################
+
+module ParticlePropagator ParticlePropagator {
   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
-#####################################
-
-module ParticlePropagator ParticlePropagator {
-  set InputArray PileUpMerger/stableParticles
   set OutputArray stableParticles
   set ChargedHadronOutputArray chargedHadrons
@@ -265 +204 @@
 
 
-########################################
-#   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}
-}
+#### ADD TRACKSMEARING
+#### ADD TIMESMEARING
+
 
 
@@ -349 +217 @@
 module SimpleCalorimeter ECal {
   set ParticleInputArray ParticlePropagator/stableParticles
-  set TrackInputArray TimeSmearing/tracks
+  set TrackInputArray TrackMerger/tracks
 
   set TowerOutputArray ecalTowers
@@ -506 +374 @@
 }
 
-
 #################
 # Electron filter
@@ -627 +494 @@
   set MomentumOutputArray momentum
 }
+
+
 
 
@@ -1004 +873 @@
   add Branch GenMissingET/momentum GenMissingET MissingET
 
-  add Branch TimeSmearing/tracks Track Track
+  add Branch TrackMerger/tracks Track Track
   add Branch Calorimeter/towers Tower Tower
 
@@ -1023 +892 @@
   add Branch MissingET/momentum MissingET MissingET
   add Branch ScalarHT/energy ScalarHT ScalarHT
-  add Branch VertexFinderDA4D/vertices Vertex4D Vertex
-}
-
+}
+