Changes in cards/FCC/FCChh_PileUpVtx.tcl [9c52415:364dbe1] in git

File:

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

cards/FCC/FCChh_PileUpVtx.tcl

@@ -12 +12 @@
 #######################################
 
-set MaxEvents 100
+set MaxEvents 1000
 set RandomSeed 123
 
 set ExecutionPath {
 
+  BeamSpotFilter
   PileUpMerger
   ParticlePropagator
@@ -28 +29 @@
   MuonMomentumSmearing
 
+  EnergyLoss
+
   TrackMerger
 
+
   TrackSmearing
-  TimeSmearing      
+  TimeSmearing  
+
+  VertexFinderDA4D  
+
+  TrackTimingPileUpSubtractor  
 
   ECal
@@ -39 +47 @@
   EFlowMerger
   EFlowFilter
-
-  TimeSmearingMIP
-  TimeSmearingPhotons
-  TimeSmearingNH    
-
-  VertexFinderDA4D
-  TrackTimingPileUpSubtractor
-
-  HighMassVertexRecover    
 
   PhotonEfficiency
@@ -82 +81 @@
 
   TreeWriter
+}
+
+#######################
+# GenBeamSpotFilter
+# Saves a particle intended to represent the beamspot
+#######################
+
+module BeamSpotFilter BeamSpotFilter {
+    set InputArray Delphes/stableParticles
+    set OutputArray beamSpotParticle
+
 }
 
@@ -142 +152 @@
   set Bz 4.0
 }
-
 
 
@@ -245 +254 @@
 }
 
+#################################
+# Charged Energy Loss (dE/dx)
+#################################
+
+module EnergyLoss EnergyLoss {
+  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
+  add InputArray ElectronMomentumSmearing/electrons
+  add InputArray MuonMomentumSmearing/muons
+
+  set ActiveFraction 0.013
+  set ChargeCollectionEfficiency 0.75
+  
+  # fixme: this number should probably be charge/energy dependent, or absolute number in MeV/cm?
+  set Resolution 0.15
+
+  # active material properties (cf. http://pdg.lbl.gov/2014/AtomicNuclearProperties/properties8.dat)
+  set Z 14.
+  set A 28.0855 
+  set rho 2.329 
+  
+  # material polarisation correction parameters
+  set a 0.1492
+  set m 3.2546
+  set x0 0.2015
+  set x1 2.8716
+  set I 173.0 
+  set c0 4.4355
+
+}
+
 ##############
 # Track merger
@@ -258 +297 @@
 
 
+
 ########################################
 #   Smear tracks
@@ -281 +321 @@
 
   # assume 20 ps resolution for now
-  set TimeResolution {20E-12}
+  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}
 }
 
@@ -450 +539 @@
 }
 
+
 #################
 # Electron filter
@@ -516 +606 @@
 }
 
-########################################
-#   Time Smearing Neutral MIP
-########################################
-
-module TimeSmearing TimeSmearingMIP {
-  set InputArray HCal/eflowTracks
-  set OutputArray timeSmearingMIP
-
-  # assume 30 ps resolution for now
-  set TimeResolution {30E-12}
-}
-
-########################################
-#   Time Smearing Neutral Photons
-########################################
-
-module TimeSmearing TimeSmearingPhotons {
-  set InputArray ECal/eflowPhotons
-  set OutputArray timeSmearingPhotons
-}
-
-########################################
-#   Time Smearing Neutral NeutralHadrons
-########################################
-#
-module TimeSmearing TimeSmearingNH {
-  set InputArray HCal/eflowNeutralHadrons
-  set OutputArray timeSmearingNH
-
-  # assume 30 ps resolution for now
-  set TimeResolution {30E-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 TimeSmearing/tracks
-  add InputArray TimeSmearingPhotons/timeSmearingPhotons
-  add InputArray TimeSmearingNH/timeSmearingNH
-
-  set VertexInputArray VertexFinderDA4D/vertices
-
-  set fChargedMinSignificance {3}
-  set fNeutralMinSignificance {3}
-}
-
-######################################
-# Heavy(slow) particles vertex recover
-######################################
-
-module HighMassVertexRecover HighMassVertexRecover {
-  set TrackInputArray VertexFinderDA4D/tracks
-  set VertexInputArray VertexFinderDA4D/vertices
-
-  set TrackOutputArray tracks
-  set VertexOutputArray vertices
-
-  set Verbose 0
-
-}
 
 ###################
@@ -1064 +1057 @@
   add Branch ScalarHT/energy ScalarHT ScalarHT
   add Branch VertexFinderDA4D/vertices Vertex4D Vertex
-
-  add Branch HighMassVertexRecover/tracks Track Track
-}
-
+}
+