Fork me on GitHub

Changeset 3c9e6d2 in git for examples


Ignore:
Timestamp:
Aug 15, 2013, 2:11:25 PM (11 years ago)
Author:
pavel <pavel@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
3bb5979
Parents:
28f1285
Message:

synchronize CMS cards

File:
1 edited

Legend:

Unmodified
Added
Removed

  • examples/delphes_card_CMS_PileUp.tcl

    r28f1285 r3c9e6d2  
    44
    55set ExecutionPath {
    6 
    7   Weighter
    86
    97  PileUpMerger
    … …  
    2927  JetPileUpSubtractor
    3028
    31   ConstituentFilter
     29  JetEnergyScale
    3230
    3331  PhotonEfficiency
    … …  
    5048
    5149  TreeWriter
    52 }
    53 
    54 ##########
    55 # Weighter
    56 ##########
    57 
    58 module Weighter Weighter {
    59   set InputArray Delphes/allParticles
    60 
    61   set OutputArray weight
    62 
    63   # add Weight {PID1 PID2} {weight}
    64   # default weight
    65   add Weight {0} {1.0}
    66 
    67   add Weight {23} {0.5}
    68   add Weight {24} {0.5}
    69   add Weight {-24} {0.5}
    70   add Weight {23 24} {0.25}
    71   add Weight {23 -24} {0.25}
    7250}
    7351
    … …  
    181159
    182160  # resolution formula for charged hadrons
    183   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    184                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    185                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    186                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    187                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    188                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    189                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
    190                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
     161  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0e1) * (0.20) + \
     162                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \
     163                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.20) + \
     164                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.20) + \
     165                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \
     166                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.20)}
    191167}
    192168
    … …  
    201177  # set ResolutionFormula {resolution formula as a function of eta and energy}
    202178
    203   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
    204                                            (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
    205                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
     179  # resolution formula for electrons
     180  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (0.0225)*energy + \
     181                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
     182                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    206183                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    207184
    … …  
    219196
    220197  # resolution formula for muons
    221   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    222                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    223                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    224                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    225                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
    226                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.03) + \
    227                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
    228                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
     198  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
     199                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
     200                                           (abs(eta) <= 0.5) * (pt > 1.0e2)                * (0.03) + \
     201                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
     202                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
     203                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.04) + \
     204                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
     205                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
     206                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.05)}
    229207}
    230208
    … …  
    308286
    309287  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    310   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
     288  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    311289                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    312290
    313291  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    314   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
    315                              (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
     292  set HCalResolutionFormula {                  (abs(eta) <= 1.5) * sqrt(energy^2*0.055^2 + energy*1.20^2 + 3.50^2) + \
     293                                               (abs(eta) > 1.5 && abs(eta) <= 3)   * sqrt(energy^2*0.065^2 + energy*2.00^2 + 3.50^2) + \
     294                                               (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*3.13^2)}
    316295}
    317296
    … …  
    402381}
    403382
    404 ####################
    405 # Constituent filter
    406 ####################
    407 
    408 module ConstituentFilter ConstituentFilter {
    409 
    410 # add JetInputArray InputArray
    411   add JetInputArray GenJetFinder/jets
    412   add JetInputArray FastJetFinder/jets
    413 
    414 # add ConstituentInputArray InputArray OutputArray
    415   add ConstituentInputArray Delphes/stableParticles stableParticles
    416   add ConstituentInputArray TrackPileUpSubtractor/eflowTracks eflowTracks
    417   add ConstituentInputArray Calorimeter/eflowTowers eflowTowers
    418   add ConstituentInputArray MuonMomentumSmearing/muons muons
    419 }
    420 
    421383###########################
    422384# Jet Pile-Up Subtraction
    … …  
    430392
    431393  set JetPTMin 20.0
     394}
     395
     396##################
     397# Jet Energy Scale
     398##################
     399
     400module EnergyScale JetEnergyScale {
     401  set InputArray JetPileUpSubtractor/jets
     402  set OutputArray jets
     403
     404 # scale formula for jets
     405  set ScaleFormula {1.0}
    432406}
    433407
    … …  
    556530  add InputArray UniqueObjectFinder/electrons
    557531  add InputArray UniqueObjectFinder/photons
    558   add InputArray MuonIsolation/muons
     532  add InputArray UniqueObjectFinder/muons
    559533  set EnergyOutputArray energy
    560534}
    … …  
    567541module BTagging BTagging {
    568542  set PartonInputArray Delphes/partons
    569   set JetInputArray JetPileUpSubtractor/jets
     543  set JetInputArray JetEnergyScale/jets
    570544
    571545  set BitNumber 0
    … …  
    601575  set ParticleInputArray Delphes/allParticles
    602576  set PartonInputArray Delphes/partons
    603   set JetInputArray JetPileUpSubtractor/jets
     577  set JetInputArray JetEnergyScale/jets
    604578
    605579  set DeltaR 0.5
    … …  
    626600  add InputArray PhotonIsolation/photons photons
    627601  add InputArray ElectronIsolation/electrons electrons
    628   add InputArray JetPileUpSubtractor/jets jets
     602  add InputArray MuonIsolation/muons muons
     603  add InputArray JetEnergyScale/jets jets
    629604}
    630605
    … …  
    645620  add Branch UniqueObjectFinder/electrons Electron Electron
    646621  add Branch UniqueObjectFinder/photons Photon Photon
    647   add Branch MuonIsolation/muons Muon Muon
     622  add Branch UniqueObjectFinder/muons Muon Muon
    648623  add Branch MissingET/momentum MissingET MissingET
    649624  add Branch ScalarHT/energy ScalarHT ScalarHT
    650625  add Branch Rho/rho Rho Rho
    651   add Branch Weighter/weight Weight Weight
    652 }
    653 
     626}
     627
Note: See TracChangeset for help on using the changeset viewer.