Fork me on GitHub

Changeset d4e5b6a in git for cards


Ignore:
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

Location:
cards
Files:
3 added
2 edited

Legend:

Unmodified
Added
Removed
  • cards/CMS_PhaseII/propagatorTest.tcl

    r2b5ff2c rd4e5b6a  
    3434
    3535  VertexFinderDA4D
     36  HighMassVertexRecover
     37
     38  CandidateFilter
    3639
    3740  TreeWriter 
     
    236239}
    237240
     241######################################
     242# Heavy(slow) particles vertex recover
     243######################################
     244
     245module HighMassVertexRecover HighMassVertexRecover {
     246  set TrackInputArray VertexFinderDAClusterizerZT/tracks
     247  set VertexInputArray VertexFinderDAClusterizerZT/vertices
     248
     249  set TrackOutputArray tracks
     250  set VertexOutputArray vertices
     251
     252  set Verbose 0
     253}
     254
     255
     256########################################
     257# Remove uninteresting tracks          #
     258########################################
     259
     260module CandidateFilter CandidateFilter {
     261  set InputArray HighMassVertexRecover/tracks
     262  set OutputArray tracks
     263
     264  set PtMin 10
     265  set MassMin 0.2
     266}
     267
     268
    238269##################
    239270# ROOT tree writer
     
    243274# add Branch InputArray BranchName BranchClass
    244275  add Branch PileUpMerger/stableParticles Particle GenParticle
     276  add Branch PileUpMerger/vertices GenVertex Vertex 
    245277  add Branch TimeSmearing/tracks Track Track
    246278  add Branch VertexFinderDA4D/vertices Vertex4D Vertex
    247   add Branch PileUpMerger/vertices GenVertex Vertex
    248 }
     279  add Branch HighMassVertexRecover/vertices Vertex4D Vertex
     280  add Branch CandidateFilter/tracks Track Track 
     281}
  • cards/FCC/FCChh_PileUpVtx.tcl

    r2b5ff2c rd4e5b6a  
    1212#######################################
    1313
     14set MaxEvents 100
     15set RandomSeed 123
     16
    1417set ExecutionPath {
    1518
     19  BeamSpotFilter
     20  PileUpMerger
    1621  ParticlePropagator
    1722
     
    2530
    2631  TrackMerger
     32
     33  TrackSmearing
     34  TimeSmearing 
     35
     36  VertexFinderDA4D 
     37
     38  TrackTimingPileUpSubtractor 
    2739
    2840  ECal
     
    6880}
    6981
     82#######################
     83# GenBeamSpotFilter
     84# Saves a particle intended to represent the beamspot
     85#######################
     86
     87module BeamSpotFilter BeamSpotFilter {
     88    set InputArray Delphes/stableParticles
     89    set OutputArray beamSpotParticle
     90
     91}
     92
     93###############
     94# PileUp Merger
     95###############
     96
     97module PileUpMerger PileUpMerger {
     98  set InputArray Delphes/stableParticles
     99
     100  set ParticleOutputArray stableParticles
     101  set VertexOutputArray vertices
     102
     103  # pre-generated minbias input file
     104  set PileUpFile MinBias.pileup
     105
     106  # average expected pile up
     107  set MeanPileUp 10
     108
     109  # 0-poisson, 1-uniform, 2-delta
     110  set PileUpDistribution 2
     111
     112  # maximum spread in the beam direction in m
     113  set ZVertexSpread 0.25
     114
     115  # maximum spread in time in s
     116  set TVertexSpread 800E-12
     117
     118  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
     119
     120  #set VertexDistributionFormula {exp(-(t^2/(2*(0.063/2.99792458E8*exp(-(z^2/(2*(0.063)^2))))^2)))}
     121  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
     122
     123  # taking 5.3 cm x 160 ps
     124
     125  #set VertexDistributionFormula { (abs(t) <= 160e-12) * (abs(z) <= 0.053) * (1.00) +
     126  #                                (abs(t) >  160e-12) * (abs(z) <= 0.053) * (0.00) +
     127  #                               (abs(t) <= 160e-12) * (abs(z) > 0.053)  * (0.00) +
     128  #                           (abs(t) >  160e-12) * (abs(z) > 0.053)  * (0.00)}
     129
     130}
    70131
    71132#################################
     
    74135
    75136module ParticlePropagator ParticlePropagator {
    76   set InputArray Delphes/stableParticles
     137  set InputArray PileUpMerger/stableParticles
    77138  set OutputArray stableParticles
    78139  set ChargedHadronOutputArray chargedHadrons
     
    204265
    205266
    206 #### ADD TRACKSMEARING
    207 #### ADD TIMESMEARING
    208 
     267########################################
     268#   Smear tracks
     269########################################
     270
     271module TrackSmearing TrackSmearing {
     272  set InputArray TrackMerger/tracks
     273  set OutputArray tracks
     274  set ApplyToPileUp true
     275
     276  # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
     277  source trackResolutionCMS.tcl
     278  # FIXME !!!! we need to add track resolution of FCC-hh baseline detector !!!!!
     279}
     280
     281########################################
     282#   Time Smearing
     283########################################
     284
     285module TimeSmearing TimeSmearing {
     286  set InputArray TrackSmearing/tracks
     287  set OutputArray tracks
     288
     289  # assume 20 ps resolution for now
     290  set TimeResolution 20E-12
     291}
     292
     293##################################
     294# Primary vertex reconstruction
     295##################################
     296
     297
     298module VertexFinderDA4D VertexFinderDA4D {
     299  set InputArray TimeSmearing/tracks
     300
     301  set OutputArray tracks
     302  set VertexOutputArray vertices
     303
     304  set Verbose 0
     305  set MinPT 1.0
     306
     307  # in mm
     308  set VertexSpaceSize 0.5
     309
     310  # in s
     311  set VertexTimeSize 10E-12
     312
     313  set UseTc 1
     314  set BetaMax 0.1
     315  set BetaStop 1.0
     316  set CoolingFactor 0.8
     317  set MaxIterations 100
     318
     319  # in mm
     320  set DzCutOff 40
     321  set D0CutOff 30
     322
     323}
     324
     325##########################
     326# Track pile-up subtractor
     327##########################
     328
     329module TrackTimingPileUpSubtractor TrackTimingPileUpSubtractor {
     330# add InputArray InputArray OutputArray
     331
     332  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
     333  add InputArray ElectronMomentumSmearing/electrons
     334  add InputArray MuonMomentumSmearing/muons
     335 
     336  set VertexInputArray VertexFinderDA4D/vertices
     337  # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
     338  # Z vertex resolution in m
     339  set ZVertexResolution {0.0001}
     340}
    209341
    210342
     
    217349module SimpleCalorimeter ECal {
    218350  set ParticleInputArray ParticlePropagator/stableParticles
    219   set TrackInputArray TrackMerger/tracks
     351  set TrackInputArray TimeSmearing/tracks
    220352
    221353  set TowerOutputArray ecalTowers
     
    374506}
    375507
     508
    376509#################
    377510# Electron filter
     
    494627  set MomentumOutputArray momentum
    495628}
    496 
    497 
    498629
    499630
     
    8731004  add Branch GenMissingET/momentum GenMissingET MissingET
    8741005
    875   add Branch TrackMerger/tracks Track Track
     1006  add Branch TimeSmearing/tracks Track Track
    8761007  add Branch Calorimeter/towers Tower Tower
    8771008
     
    8921023  add Branch MissingET/momentum MissingET MissingET
    8931024  add Branch ScalarHT/energy ScalarHT ScalarHT
    894 }
    895 
     1025  add Branch VertexFinderDA4D/vertices Vertex4D Vertex
     1026}
     1027
Note: See TracChangeset for help on using the changeset viewer.