Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • examples/delphes_card_CMS_PileUp.tcl

    r6fb1a5d r282f591  
     1
    12#######################################
    23# Order of execution of various modules
     
    2021  TrackPileUpSubtractor
    2122  NeutralTowerMerger
    22   EFlowMergerAllTracks
    2323  EFlowMerger
    24  
    25   NeutrinoFilter
     24
    2625  GenJetFinder
    2726
     
    6867
    6968  # average expected pile up
    70   set MeanPileUp 50
     69  set MeanPileUp 10
    7170
    7271  # maximum spread in the beam direction in m
     
    254253  set TowerOutputArray towers
    255254  set PhotonOutputArray photons
    256 
    257   set EcalTowerMinEnergy 0.5
    258   set HcalTowerMinEnergy 1.0
    259 
    260   set EcalTowerMinSignificance 1.0
    261   set HcalTowerMinSignificance 1.0
    262255
    263256  set EFlowTrackOutputArray eflowTracks
     
    354347}
    355348
    356 ##################################
    357 # Energy flow merger (all tracks)
    358 ##################################
    359 
    360 module Merger EFlowMergerAllTracks {
    361 # add InputArray InputArray
    362   add InputArray TrackMerger/tracks
    363   add InputArray Calorimeter/eflowPhotons
    364   add InputArray Calorimeter/eflowNeutralHadrons
    365   set OutputArray eflow
    366 }
    367 
    368349
    369350####################
     
    380361
    381362
    382 
    383363#############
    384364# Rho pile-up
    385365#############
    386366
    387 module FastJetGridMedianEstimator Rho {
    388  
     367module FastJetFinder Rho {
     368#  set InputArray Calorimeter/towers
    389369  set InputArray EFlowMerger/eflow
     370
     371  set ComputeRho true
    390372  set RhoOutputArray rho
    391373
    392   # etamin etamax gridsize_eta gridsize_phi
    393  
    394   add GridRange 0.0 2.5 1.0 1.0
    395   add GridRange 2.5 5.0 1.0 1.0
    396 
    397 }
    398 
    399 #####################
    400 # Neutrino Filter
    401 #####################
    402 
    403 module PdgCodeFilter NeutrinoFilter {
    404  
    405   set InputArray Delphes/stableParticles
    406   set OutputArray filteredParticles
    407 
    408   set PTMin 0.0
    409  
    410   add PdgCode {12}
    411   add PdgCode {14}
    412   add PdgCode {16}
    413   add PdgCode {-12}
    414   add PdgCode {-14}
    415   add PdgCode {-16}
    416 
    417 }
    418 
    419 
     374  # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
     375  set AreaAlgorithm 5
     376
     377  # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     378  set JetAlgorithm 4
     379  set ParameterR 0.6
     380  set GhostEtaMax 5.0
     381
     382  add RhoEtaRange 0.0 2.5
     383  add RhoEtaRange 2.5 5.0
     384
     385  set JetPTMin 0.0
     386}
    420387
    421388#####################
     
    424391
    425392module FastJetFinder GenJetFinder {
    426   set InputArray NeutrinoFilter/filteredParticles
     393  set InputArray Delphes/stableParticles
    427394
    428395  set OutputArray jets
     
    616583module Merger MissingET {
    617584# add InputArray InputArray
    618   add InputArray EFlowMergerAllTracks/eflow
     585  add InputArray EFlowMerger/eflow
    619586  set MomentumOutputArray momentum
    620587}
Note: See TracChangeset for help on using the changeset viewer.