Fork me on GitHub

Ignore:
Timestamp:
Oct 26, 2013, 10:26:57 PM (11 years ago)
Author:
Pavel Demin
Message:

synchronize delphes_card_CMS_NoVFC.tcl with other CMS cards

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/examples/delphes_card_CMS_NoVFC.tcl

    r1099 r1302  
    3232  FastJetFinder
    3333
    34   ConstituentFilter
     34  JetEnergyScale
    3535
    3636  BTagging
     
    155155  # set ResolutionFormula {resolution formula as a function of eta and energy}
    156156
    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)}
     157  # resolution formula for electrons
     158  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) + \
     159                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
     160                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)}
    160161
    161162}
     
    172173
    173174  # resolution formula for muons
    174   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    175                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    176                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    177                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    178                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
    179                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.03) + \
    180                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
    181                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
     175  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
     176                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
     177                                           (abs(eta) <= 0.5) * (pt > 1.0e2)                * (0.03) + \
     178                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
     179                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
     180                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.04) + \
     181                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
     182                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
     183                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.05)}
    182184}
    183185
     
    190192  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    191193  add InputArray ElectronEnergySmearing/electrons
     194  add InputArray MuonMomentumSmearing/muons
    192195  set OutputArray tracks
    193196}
     
    252255
    253256  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    254   set ECalResolutionFormula {sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2)}
     257  set ECalResolutionFormula {sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)}
    255258
    256259  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
     
    266269  add InputArray Calorimeter/eflowTracks
    267270  add InputArray Calorimeter/eflowTowers
    268   add InputArray MuonMomentumSmearing/muons
    269271  set OutputArray eflow
    270272}
     
    378380module Merger MissingET {
    379381# add InputArray InputArray
    380   add InputArray EFlowMerger/eflow
     382  add InputArray Calorimeter/eflowTracks
     383  add InputArray Calorimeter/eflowTowers
    381384  set MomentumOutputArray momentum
    382385}
     
    391394  add InputArray UniqueObjectFinder/electrons
    392395  add InputArray UniqueObjectFinder/photons
    393   add InputArray MuonIsolation/muons
     396  add InputArray UniqueObjectFinder/muons
    394397  set EnergyOutputArray energy
    395398}
     
    406409  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    407410  set JetAlgorithm 6
    408   set ParameterR 0.7
    409 
    410   set ConeRadius 0.5
    411   set SeedThreshold 1.0
    412   set ConeAreaFraction 1.0
    413   set AdjacencyCut 2.0
    414   set OverlapThreshold 0.75
    415 
    416   set MaxIterations 100
    417   set MaxPairSize 2
    418   set Iratch 1
     411  set ParameterR 0.5
    419412
    420413  set JetPTMin 20.0
     
    433426  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    434427  set JetAlgorithm 6
    435   set ParameterR 0.7
    436 
    437   set ConeRadius 0.5
    438   set SeedThreshold 1.0
    439   set ConeAreaFraction 1.0
    440   set AdjacencyCut 2.0
    441   set OverlapThreshold 0.75
    442 
    443   set MaxIterations 100
    444   set MaxPairSize 2
    445   set Iratch 1
     428  set ParameterR 0.5
    446429
    447430  set JetPTMin 20.0
    448431}
    449432
    450 ####################
    451 # Constituent filter
    452 ####################
    453 
    454 module ConstituentFilter ConstituentFilter {
    455 
    456 # add JetInputArray InputArray
    457   add JetInputArray GenJetFinder/jets
    458   add JetInputArray FastJetFinder/jets
    459 
    460 # add ConstituentInputArray InputArray OutputArray
    461   add ConstituentInputArray Delphes/stableParticles stableParticles
    462   add ConstituentInputArray Calorimeter/eflowTracks eflowTracks
    463   add ConstituentInputArray Calorimeter/eflowTowers eflowTowers
    464   add ConstituentInputArray MuonMomentumSmearing/muons muons
     433##################
     434# Jet Energy Scale
     435##################
     436
     437module EnergyScale JetEnergyScale {
     438  set InputArray FastJetFinder/jets
     439  set OutputArray jets
     440
     441 # scale formula for jets
     442  set ScaleFormula {1.08}
    465443}
    466444
     
    471449module BTagging BTagging {
    472450  set PartonInputArray Delphes/partons
    473   set JetInputArray FastJetFinder/jets
     451  set JetInputArray JetEnergyScale/jets
    474452
    475453  set BitNumber 0
     
    505483  set ParticleInputArray Delphes/allParticles
    506484  set PartonInputArray Delphes/partons
    507   set JetInputArray FastJetFinder/jets
     485  set JetInputArray JetEnergyScale/jets
    508486
    509487  set DeltaR 0.5
     
    530508  add InputArray PhotonIsolation/photons photons
    531509  add InputArray ElectronIsolation/electrons electrons
    532   add InputArray FastJetFinder/jets jets
     510  add InputArray MuonIsolation/muons muons
     511  add InputArray JetEnergyScale/jets jets
    533512}
    534513
     
    542521  add Branch TrackMerger/tracks Track Track
    543522  add Branch Calorimeter/towers Tower Tower
    544   add Branch ConstituentFilter/eflowTracks EFlowTrack Track
    545   add Branch ConstituentFilter/eflowTowers EFlowTower Tower
    546   add Branch ConstituentFilter/muons EFlowMuon Muon
     523  add Branch Calorimeter/eflowTracks EFlowTrack Track
     524  add Branch Calorimeter/eflowTowers EFlowTower Tower
    547525  add Branch GenJetFinder/jets GenJet Jet
    548526  add Branch UniqueObjectFinder/jets Jet Jet
    549527  add Branch UniqueObjectFinder/electrons Electron Electron
    550528  add Branch UniqueObjectFinder/photons Photon Photon
    551   add Branch MuonIsolation/muons Muon Muon
     529  add Branch UniqueObjectFinder/muons Muon Muon
    552530  add Branch MissingET/momentum MissingET MissingET
    553531  add Branch ScalarHT/energy ScalarHT ScalarHT
Note: See TracChangeset for help on using the changeset viewer.