Fork me on GitHub

Ignore:
Timestamp:
Aug 15, 2013, 1:32:50 PM (11 years ago)
Author:
Pavel Demin
Message:

update resolution, remove ConstituentFilter and add JetEnergyScale

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/examples/delphes_card_CMS.tcl

    r1099 r1254  
    3232  FastJetFinder
    3333
    34   ConstituentFilter
     34  JetEnergyScale
    3535
    3636  BTagging
    … …  
    135135
    136136  # resolution formula for charged hadrons
    137   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    138                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    139                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    140                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    141                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    142                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    143                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
    144                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
     137  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0e1) * (0.20) + \
     138                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \
     139                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.20) + \
     140                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.20) + \
     141                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.20) + \
     142                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.20)}
    145143}
    146144
    … …  
    155153  # set ResolutionFormula {resolution formula as a function of eta and energy}
    156154
    157   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
    158                                            (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
    159                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
     155  # resolution formula for electrons
     156  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (0.0225)*energy + \
     157                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
     158                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    160159                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    161160
    … …  
    173172
    174173  # resolution formula for muons
    175   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    176                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    177                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    178                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    179                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
    180                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.03) + \
    181                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
    182                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
     174  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
     175                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
     176                                           (abs(eta) <= 0.5) * (pt > 1.0e2)                * (0.03) + \
     177                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
     178                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
     179                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.04) + \
     180                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
     181                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
     182                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.05)}
    183183}
    184184
    … …  
    262262
    263263  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    264   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
     264  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    265265                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    266266
    267267  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    268   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
    269                              (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
     268  set HCalResolutionFormula {                  (abs(eta) <= 1.5) * sqrt(energy^2*0.055^2 + energy*1.20^2 + 3.50^2) + \
     269                                               (abs(eta) > 1.5 && abs(eta) <= 3)   * sqrt(energy^2*0.065^2 + energy*2.00^2 + 3.50^2) + \
     270                                               (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*3.13^2)}
    270271}
    271272
    … …  
    403404  add InputArray UniqueObjectFinder/electrons
    404405  add InputArray UniqueObjectFinder/photons
    405   add InputArray MuonIsolation/muons
     406  add InputArray UniqueObjectFinder/muons
    406407  set EnergyOutputArray energy
    407408}
    … …  
    418419  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    419420  set JetAlgorithm 6
    420   set ParameterR 0.7
    421 
    422   set ConeRadius 0.5
    423   set SeedThreshold 1.0
    424   set ConeAreaFraction 1.0
    425   set AdjacencyCut 2.0
    426   set OverlapThreshold 0.75
    427 
    428   set MaxIterations 100
    429   set MaxPairSize 2
    430   set Iratch 1
     421  set ParameterR 0.5
    431422
    432423  set JetPTMin 20.0
    … …  
    445436  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    446437  set JetAlgorithm 6
    447   set ParameterR 0.7
    448 
    449   set ConeRadius 0.5
    450   set SeedThreshold 1.0
    451   set ConeAreaFraction 1.0
    452   set AdjacencyCut 2.0
    453   set OverlapThreshold 0.75
    454 
    455   set MaxIterations 100
    456   set MaxPairSize 2
    457   set Iratch 1
     438  set ParameterR 0.5
    458439
    459440  set JetPTMin 20.0
    460441}
    461442
    462 ####################
    463 # Constituent filter
    464 ####################
    465 
    466 module ConstituentFilter ConstituentFilter {
    467 
    468 # add JetInputArray InputArray
    469   add JetInputArray GenJetFinder/jets
    470   add JetInputArray FastJetFinder/jets
    471 
    472 # add ConstituentInputArray InputArray OutputArray
    473   add ConstituentInputArray Delphes/stableParticles stableParticles
    474   add ConstituentInputArray Calorimeter/eflowTracks eflowTracks
    475   add ConstituentInputArray Calorimeter/eflowTowers eflowTowers
    476   add ConstituentInputArray MuonMomentumSmearing/muons muons
     443##################
     444# Jet Energy Scale
     445##################
     446
     447module EnergyScale JetEnergyScale {
     448  set InputArray FastJetFinder/jets
     449  set OutputArray jets
     450
     451 # scale formula for jets
     452  set ScaleFormula {1.08}
    477453}
    478454
    … …  
    483459module BTagging BTagging {
    484460  set PartonInputArray Delphes/partons
    485   set JetInputArray FastJetFinder/jets
     461  set JetInputArray JetEnergyScale/jets
    486462
    487463  set BitNumber 0
    … …  
    542518  add InputArray PhotonIsolation/photons photons
    543519  add InputArray ElectronIsolation/electrons electrons
    544   add InputArray FastJetFinder/jets jets
     520  add InputArray MuonIsolation/muons muons
     521  add InputArray JetEnergyScale/jets jets
    545522}
    546523
    … …  
    554531  add Branch TrackMerger/tracks Track Track
    555532  add Branch Calorimeter/towers Tower Tower
    556   add Branch ConstituentFilter/eflowTracks EFlowTrack Track
    557   add Branch ConstituentFilter/eflowTowers EFlowTower Tower
    558   add Branch ConstituentFilter/muons EFlowMuon Muon
     533  add Branch Calorimeter/eflowTracks EFlowTrack Track
     534  add Branch Calorimeter/eflowTowers EFlowTower Tower
     535  add Branch MuonMomentumSmearing/muons EFlowMuon Muon
    559536  add Branch GenJetFinder/jets GenJet Jet
    560537  add Branch UniqueObjectFinder/jets Jet Jet
    561538  add Branch UniqueObjectFinder/electrons Electron Electron
    562539  add Branch UniqueObjectFinder/photons Photon Photon
    563   add Branch MuonIsolation/muons Muon Muon
     540  add Branch UniqueObjectFinder/muons Muon Muon
    564541  add Branch MissingET/momentum MissingET MissingET
    565542  add Branch ScalarHT/energy ScalarHT ScalarHT
Note: See TracChangeset for help on using the changeset viewer.