Fork me on GitHub

Changeset 77249aa in git for cards/FCC/FCChh_PileUpVtx.tcl


Ignore:
Timestamp:
Jan 26, 2020, 4:53:05 PM (5 years ago)
Author:
michele <michele.selvaggi@…>
Branches:
Timing
Children:
6777565
Parents:
f59b36de
Message:

first version of card using timing information

File:
1 edited

Legend:

Unmodified
Added
Removed
  • cards/FCC/FCChh_PileUpVtx.tcl

    rf59b36de r77249aa  
    3131
    3232  TrackMerger
    33 
    3433  TrackSmearing
    3534
    3635  ECal
    3736  HCal
     37
     38  TimeSmearingMIP
     39  TimeSmearingPhotons
     40  TimeSmearingNH
     41
     42  TimeSmearingEcal
     43  TimeSmearingHcal
     44
     45  VertexFinderDA4D
     46  HighMassVertexRecover
     47  PileUpSubtractor4D
     48
     49  ElectronFilter
     50  MuonFilter
     51  ChargedHadronFilter
    3852
    3953  Calorimeter
     
    4155  EFlowFilter
    4256
    43   TimeSmearingMIP
    44   TimeSmearingPhotons
    45   TimeSmearingNH   
    46 
    47   VertexFinderDA4D
    48   PileUpSubtractor4D
    49 
    50   HighMassVertexRecover   
    51 
    5257  PhotonEfficiency
    5358  PhotonIsolation
    5459
    55   ElectronFilter
    5660  ElectronIsolation
    57 
    58   ChargedHadronFilter
    5961
    6062  MuonIsolation
     
    256258  set ActiveFraction 0.013
    257259  set ChargeCollectionEfficiency 0.75
    258  
     260
    259261  # fixme: this number should probably be charge/energy dependent, or absolute number in MeV/cm?
    260262  set Resolution 0.15
     
    262264  # active material properties (cf. http://pdg.lbl.gov/2014/AtomicNuclearProperties/properties8.dat)
    263265  set Z 14.
    264   set A 28.0855 
    265   set rho 2.329 
    266  
     266  set A 28.0855
     267  set rho 2.329
     268
    267269  # material polarisation correction parameters
    268270  set a 0.1492
     
    270272  set x0 0.2015
    271273  set x1 2.8716
    272   set I 173.0 
     274  set I 173.0
    273275  set c0 4.4355
    274276
     
    298300  set ApplyToPileUp true
    299301
    300   # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
    301   source trackResolutionFCChh.tcl 
     302  # this is CMS resolution is d0 , dZ on top of the smearing in P applied before
     303  source trackResolutionFCChh.tcl
    302304  # FIXME !!!! we need to add track resolution of FCC-hh baseline detector !!!!!
    303305}
     
    312314module SimpleCalorimeter ECal {
    313315  set ParticleInputArray ParticlePropagator/stableParticles
    314   set TrackInputArray TimeSmearing/tracks
     316  set TrackInputArray TrackSmearing/tracks
    315317
    316318  set TowerOutputArray ecalTowers
     
    469471}
    470472
    471 #################
    472 # Electron filter
    473 #################
    474 
    475 module PdgCodeFilter ElectronFilter {
    476   set InputArray HCal/eflowTracks
    477   set OutputArray electrons
    478   set Invert true
    479   add PdgCode {11}
    480   add PdgCode {-11}
    481 }
    482 
    483 
    484 ######################
    485 # ChargedHadronFilter
    486 ######################
    487 
    488 module PdgCodeFilter ChargedHadronFilter {
    489   set InputArray HCal/eflowTracks
    490   set OutputArray chargedHadrons
    491 
    492   add PdgCode {11}
    493   add PdgCode {-11}
    494   add PdgCode {13}
    495   add PdgCode {-13}
    496 }
    497 
    498 
    499 ###################################################
    500 # Tower Merger (in case not using e-flow algorithm)
    501 ###################################################
    502 
    503 module Merger Calorimeter {
    504 # add InputArray InputArray
    505   add InputArray ECal/ecalTowers
    506   add InputArray HCal/hcalTowers
    507   add InputArray MuonMomentumSmearing/muons
    508   set OutputArray towers
    509 }
    510 
    511 ####################
    512 # Energy flow merger
    513 ####################
    514 
    515 module Merger EFlowMerger {
    516 # add InputArray InputArray
    517   add InputArray HCal/eflowTracks
    518   add InputArray ECal/eflowPhotons
    519   add InputArray HCal/eflowNeutralHadrons
    520   set OutputArray eflow
    521 }
    522 
    523 ######################
    524 # EFlowFilter
    525 ######################
    526 
    527 module PdgCodeFilter EFlowFilter {
    528   set InputArray EFlowMerger/eflow
    529   set OutputArray eflow
    530 
    531   add PdgCode {11}
    532   add PdgCode {-11}
    533   add PdgCode {13}
    534   add PdgCode {-13}
    535 }
    536473
    537474########################################
     
    543480  set OutputArray tracks
    544481
    545   # assume 30 ps resolution for now
    546   set TimeResolution {30E-12}
     482  # assume constant 30 ps resolution for now
     483  set TimeResolution {
     484                       (abs(eta) > 0.0 && abs(eta) <= 6.0)* 30E-12
     485                     }
    547486}
    548487
     
    554493  set InputArray ECal/eflowPhotons
    555494  set OutputArray photons
    556   set TimeResolution {sqrt(20^2 + 150^2)/energy^2}
     495  set TimeResolution {
     496                        (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2
     497                     }
    557498}
    558499
     
    565506  set OutputArray neutralhadrons
    566507
    567   # assume 30 ps resolution for now
    568   set TimeResolution {sqrt(20^2 + 150^2)/energy^2}
    569 }
     508  # assume resolution for now
     509  set TimeResolution {
     510                        (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2
     511                     }
     512}
     513
     514########################################
     515#   Time Smearing Neutral Ecal
     516########################################
     517
     518module TimeSmearing TimeSmearingEcal {
     519  set InputArray ECal/ecalTowers
     520  set OutputArray towers
     521  set TimeResolution {
     522                        (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2
     523                     }
     524}
     525
     526
     527########################################
     528#   Time Smearing Neutral Hcal
     529########################################
     530
     531module TimeSmearing TimeSmearingHcal {
     532  set InputArray HCal/hcalTowers
     533  set OutputArray towers
     534  set TimeResolution {
     535                        (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2
     536                     }
     537}
     538
     539
    570540
    571541
     
    573543# Primary vertex reconstruction
    574544##################################
    575 
    576545
    577546module VertexFinderDA4D VertexFinderDA4D {
     
    602571}
    603572
    604 ##########################
    605 # Track pile-up subtractor
    606 ##########################
    607 
    608 module PileUpSubtractor4D PileUpSubtractor4D {
    609 # add InputArray InputArray OutputArray
    610 
    611   add InputArray TimeSmearingMIP/tracks
    612   add InputArray TimeSmearingPhotons/photons
    613   add InputArray TimeSmearingNH/neutralhadrons
    614 
    615   set VertexInputArray VertexFinderDA4D/vertices
    616 
    617   set fChargedMinSignificance {3}
    618   set fNeutralMinSignificance {3}
    619 }
    620 
    621573######################################
    622574# Heavy(slow) particles vertex recover
     
    634586
    635587}
     588
     589##########################
     590# PileUpSubtractor4D
     591##########################
     592
     593module PileUpSubtractor4D PileUpSubtractor4D {
     594# add InputArray InputArray OutputArray
     595
     596  add InputArray TimeSmearingMIP/tracks tracks
     597  add InputArray TimeSmearingPhotons/photons photons
     598  add InputArray TimeSmearingNH/neutralhadrons neutralhadrons
     599  add InputArray TimeSmearingEcal/towers towers
     600  add InputArray TimeSmearingHcal/towers towers
     601
     602  set VertexInputArray VertexFinderDA4D/vertices
     603
     604  set fChargedMinSignificance 3.0
     605  set fNeutralMinSignificance 3.0
     606}
     607
     608
     609#################
     610# Electron filter
     611#################
     612
     613module PdgCodeFilter ElectronFilter {
     614  set InputArray PileUpSubtractor4D/tracks
     615  set OutputArray electrons
     616  set Invert true
     617  add PdgCode {11}
     618  add PdgCode {-11}
     619}
     620
     621#################
     622# Muon filter
     623#################
     624
     625module PdgCodeFilter MuonFilter {
     626  set InputArray PileUpSubtractor4D/tracks
     627  set OutputArray muons
     628  set Invert true
     629  add PdgCode {13}
     630  add PdgCode {-13}
     631}
     632
     633
     634######################
     635# ChargedHadronFilter
     636######################
     637
     638module PdgCodeFilter ChargedHadronFilter {
     639  set InputArray PileUpSubtractor4D/tracks
     640  set OutputArray chargedHadrons
     641
     642  add PdgCode {11}
     643  add PdgCode {-11}
     644  add PdgCode {13}
     645  add PdgCode {-13}
     646}
     647
     648
     649###################################################
     650# Tower Merger (in case not using e-flow algorithm)
     651###################################################
     652
     653module Merger Calorimeter {
     654# add InputArray InputArray
     655  add InputArray PileUpSubtractor4D/towers
     656  add InputArray PileUpSubtractor4D/towers
     657  add InputArray MuonFilter/muons
     658  set OutputArray towers
     659}
     660
     661####################
     662# Energy flow merger
     663####################
     664
     665module Merger EFlowMerger {
     666# add InputArray InputArray
     667  add InputArray PileUpSubtractor4D/tracks
     668  add InputArray PileUpSubtractor4D/photons
     669  add InputArray PileUpSubtractor4D/neutralhadrons
     670  set OutputArray eflow
     671}
     672
     673##############
     674# EFlowFilter
     675##############
     676
     677module PdgCodeFilter EFlowFilter {
     678  set InputArray EFlowMerger/eflow
     679  set OutputArray eflow
     680
     681  add PdgCode {11}
     682  add PdgCode {-11}
     683  add PdgCode {13}
     684  add PdgCode {-13}
     685}
     686
     687
    636688
    637689###################
     
    644696  set MomentumOutputArray momentum
    645697}
    646 
    647 
    648698
    649699##################
     
    10651115  add Branch GenMissingET/momentum GenMissingET MissingET
    10661116
    1067   add Branch TimeSmearing/tracks Track Track
     1117  add Branch TrackSmearing/tracks Track Track
    10681118  add Branch Calorimeter/towers Tower Tower
    10691119
    1070   add Branch HCal/eflowTracks EFlowTrack Track
    1071   add Branch ECal/eflowPhotons EFlowPhoton Tower
    1072   add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
     1120  add Branch PileUpSubtractor4D/tracks EFlowTrack Track
     1121  add Branch PileUpSubtractor4D/photons EFlowPhoton Tower
     1122  add Branch PileUpSubtractor4D/neutralhadrons EFlowNeutralHadron Tower
    10731123
    10741124  add Branch UniqueObjectFinder/photons Photon Photon
     
    10881138  add Branch HighMassVertexRecover/tracks Track Track
    10891139}
    1090 
Note: See TracChangeset for help on using the changeset viewer.