Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • cards/FCC/FCChh_PileUpVtx.tcl

    r364dbe1 r9c52415  
    1212#######################################
    1313
    14 set MaxEvents 1000
     14set MaxEvents 100
    1515set RandomSeed 123
    1616
    1717set ExecutionPath {
    1818
    19   BeamSpotFilter
    2019  PileUpMerger
    2120  ParticlePropagator
     
    2928  MuonMomentumSmearing
    3029
    31   EnergyLoss
    32 
    3330  TrackMerger
    3431
    35 
    3632  TrackSmearing
    37   TimeSmearing 
    38 
    39   VertexFinderDA4D 
    40 
    41   TrackTimingPileUpSubtractor 
     33  TimeSmearing     
    4234
    4335  ECal
     
    4739  EFlowMerger
    4840  EFlowFilter
     41
     42  TimeSmearingMIP
     43  TimeSmearingPhotons
     44  TimeSmearingNH   
     45
     46  VertexFinderDA4D
     47  TrackTimingPileUpSubtractor
     48
     49  HighMassVertexRecover   
    4950
    5051  PhotonEfficiency
     
    8182
    8283  TreeWriter
    83 }
    84 
    85 #######################
    86 # GenBeamSpotFilter
    87 # Saves a particle intended to represent the beamspot
    88 #######################
    89 
    90 module BeamSpotFilter BeamSpotFilter {
    91     set InputArray Delphes/stableParticles
    92     set OutputArray beamSpotParticle
    93 
    9484}
    9585
     
    152142  set Bz 4.0
    153143}
     144
    154145
    155146
     
    254245}
    255246
    256 #################################
    257 # Charged Energy Loss (dE/dx)
    258 #################################
    259 
    260 module EnergyLoss EnergyLoss {
    261   add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    262   add InputArray ElectronMomentumSmearing/electrons
    263   add InputArray MuonMomentumSmearing/muons
    264 
    265   set ActiveFraction 0.013
    266   set ChargeCollectionEfficiency 0.75
    267  
    268   # fixme: this number should probably be charge/energy dependent, or absolute number in MeV/cm?
    269   set Resolution 0.15
    270 
    271   # active material properties (cf. http://pdg.lbl.gov/2014/AtomicNuclearProperties/properties8.dat)
    272   set Z 14.
    273   set A 28.0855
    274   set rho 2.329
    275  
    276   # material polarisation correction parameters
    277   set a 0.1492
    278   set m 3.2546
    279   set x0 0.2015
    280   set x1 2.8716
    281   set I 173.0
    282   set c0 4.4355
    283 
    284 }
    285 
    286247##############
    287248# Track merger
     
    297258
    298259
    299 
    300260########################################
    301261#   Smear tracks
     
    321281
    322282  # assume 20 ps resolution for now
    323   set TimeResolution 20E-12
    324 }
    325 
    326 ##################################
    327 # Primary vertex reconstruction
    328 ##################################
    329 
    330 
    331 module VertexFinderDA4D VertexFinderDA4D {
    332   set InputArray TimeSmearing/tracks
    333 
    334   set OutputArray tracks
    335   set VertexOutputArray vertices
    336 
    337   set Verbose 0
    338   set MinPT 1.0
    339 
    340   # in mm
    341   set VertexSpaceSize 0.5
    342 
    343   # in s
    344   set VertexTimeSize 10E-12
    345 
    346   set UseTc 1
    347   set BetaMax 0.1
    348   set BetaStop 1.0
    349   set CoolingFactor 0.8
    350   set MaxIterations 100
    351 
    352   # in mm
    353   set DzCutOff 40
    354   set D0CutOff 30
    355 
    356 }
    357 
    358 ##########################
    359 # Track pile-up subtractor
    360 ##########################
    361 
    362 module TrackTimingPileUpSubtractor TrackTimingPileUpSubtractor {
    363 # add InputArray InputArray OutputArray
    364 
    365   add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    366   add InputArray ElectronMomentumSmearing/electrons
    367   add InputArray MuonMomentumSmearing/muons
    368  
    369   set VertexInputArray VertexFinderDA4D/vertices
    370   # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
    371   # Z vertex resolution in m
    372   set ZVertexResolution {0.0001}
     283  set TimeResolution {20E-12}
    373284}
    374285
     
    539450}
    540451
    541 
    542452#################
    543453# Electron filter
     
    606516}
    607517
     518########################################
     519#   Time Smearing Neutral MIP
     520########################################
     521
     522module TimeSmearing TimeSmearingMIP {
     523  set InputArray HCal/eflowTracks
     524  set OutputArray timeSmearingMIP
     525
     526  # assume 30 ps resolution for now
     527  set TimeResolution {30E-12}
     528}
     529
     530########################################
     531#   Time Smearing Neutral Photons
     532########################################
     533
     534module TimeSmearing TimeSmearingPhotons {
     535  set InputArray ECal/eflowPhotons
     536  set OutputArray timeSmearingPhotons
     537}
     538
     539########################################
     540#   Time Smearing Neutral NeutralHadrons
     541########################################
     542#
     543module TimeSmearing TimeSmearingNH {
     544  set InputArray HCal/eflowNeutralHadrons
     545  set OutputArray timeSmearingNH
     546
     547  # assume 30 ps resolution for now
     548  set TimeResolution {30E-12}
     549}
     550
     551
     552##################################
     553# Primary vertex reconstruction
     554##################################
     555
     556
     557module VertexFinderDA4D VertexFinderDA4D {
     558  set InputArray TimeSmearing/tracks
     559
     560  set OutputArray tracks
     561  set VertexOutputArray vertices
     562
     563  set Verbose 0
     564  set MinPT 1.0
     565
     566  # in mm
     567  set VertexSpaceSize 0.5
     568
     569  # in s
     570  set VertexTimeSize 10E-12
     571
     572  set UseTc 1
     573  set BetaMax 0.1
     574  set BetaStop 1.0
     575  set CoolingFactor 0.8
     576  set MaxIterations 100
     577
     578  # in mm
     579  set DzCutOff 40
     580  set D0CutOff 30
     581
     582}
     583
     584##########################
     585# Track pile-up subtractor
     586##########################
     587
     588module TrackTimingPileUpSubtractor TrackTimingPileUpSubtractor {
     589# add InputArray InputArray OutputArray
     590
     591  add InputArray TimeSmearing/tracks
     592  add InputArray TimeSmearingPhotons/timeSmearingPhotons
     593  add InputArray TimeSmearingNH/timeSmearingNH
     594
     595  set VertexInputArray VertexFinderDA4D/vertices
     596
     597  set fChargedMinSignificance {3}
     598  set fNeutralMinSignificance {3}
     599}
     600
     601######################################
     602# Heavy(slow) particles vertex recover
     603######################################
     604
     605module HighMassVertexRecover HighMassVertexRecover {
     606  set TrackInputArray VertexFinderDA4D/tracks
     607  set VertexInputArray VertexFinderDA4D/vertices
     608
     609  set TrackOutputArray tracks
     610  set VertexOutputArray vertices
     611
     612  set Verbose 0
     613
     614}
    608615
    609616###################
     
    10571064  add Branch ScalarHT/energy ScalarHT ScalarHT
    10581065  add Branch VertexFinderDA4D/vertices Vertex4D Vertex
    1059 }
    1060 
     1066
     1067  add Branch HighMassVertexRecover/tracks Track Track
     1068}
     1069
Note: See TracChangeset for help on using the changeset viewer.