Fork me on GitHub

Changeset a2983ec in git


Ignore:
Timestamp:
Dec 22, 2014, 11:38:33 AM (10 years ago)
Author:
Pavel Demin <pavel.demin@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
8c72f09
Parents:
4e09c3a
Message:

add parameter SmearTowerCenter to all cards

Location:
cards
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • cards/delphes_card_ATLAS.tcl

    r4e09c3a ra2983ec  
    7777
    7878  # tracking efficiency formula for charged hadrons
    79   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    80                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    81                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    82                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    83                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
     79  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     80                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     81                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     82                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     83                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    8484                         (abs(eta) > 2.5)                                                  * (0.00)}
    8585}
     
    9696
    9797  # tracking efficiency formula for electrons
    98   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    99                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
    100                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
    101                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
    102                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
    103                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
    104                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
     98  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     99                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     100                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     101                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     102                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     103                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     104                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    105105                         (abs(eta) > 2.5)                                                  * (0.00)}
    106106}
     
    117117
    118118  # tracking efficiency formula for muons
    119   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    120                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    121                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    122                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    123                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
     119  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     120                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     121                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     122                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     123                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    124124                         (abs(eta) > 2.5)                                                  * (0.00)}
    125125}
     
    136136
    137137  # resolution formula for charged hadrons
    138   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    139                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    140                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    141                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    142                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    143                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    144                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
     138  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     139                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     140                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     141                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     142                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     143                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     144                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    145145                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    146146}
     
    156156  # set ResolutionFormula {resolution formula as a function of eta and energy}
    157157
    158   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
    159                                            (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
    160                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
     158  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) +
     159                                           (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
     160                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
    161161                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    162162
     
    174174
    175175  # resolution formula for muons
    176   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    177                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) + \
    178                                            (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
    179                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) + \
    180                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
    181                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) + \
    182                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
     176  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     177                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) +
     178                                           (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) +
     179                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) +
     180                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) +
     181                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) +
     182                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) +
    183183                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.10)}
    184184}
     
    207207  set PhotonOutputArray photons
    208208
    209   set ECalEnergyMin 0.50
    210   set HCalEnergyMin 1.00
    211 
    212   set ECalEnergySignificanceMin 1.0
    213   set HCalEnergySignificanceMin 1.0
    214  
    215209  set EFlowTrackOutputArray eflowTracks
    216210  set EFlowPhotonOutputArray eflowPhotons
    217211  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     212
     213  set ECalEnergyMin 0.5
     214  set HCalEnergyMin 1.0
     215
     216  set ECalEnergySignificanceMin 1.0
     217  set HCalEnergySignificanceMin 1.0
     218
     219  set SmearTowerCenter true
    218220
    219221  set pi [expr {acos(-1)}]
     
    265267  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    266268  # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
    267   set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
     269  set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
    268270                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
    269271
     
    271273  # http://arxiv.org/pdf/hep-ex/0004009v1
    272274  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    273   set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
    274                              (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
     275  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
     276                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
    275277                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
    276278}
     
    299301
    300302  # efficiency formula for photons
    301   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    302                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    303                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     303  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     304                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     305                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    304306                         (abs(eta) > 2.5)                                   * (0.00)}
    305307}
     
    333335
    334336  # efficiency formula for electrons
    335   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    336                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    337                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     337  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     338                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     339                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    338340                         (abs(eta) > 2.5)                                   * (0.00)}
    339341}
     
    367369
    368370  # efficiency formula for muons
    369   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    370                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    371                          (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) + \
     371  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     372                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     373                         (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) +
    372374                         (abs(eta) > 2.7)                                   * (0.00)}
    373375}
     
    419421
    420422module PdgCodeFilter NeutrinoFilter {
    421  
     423
    422424  set InputArray Delphes/stableParticles
    423425  set OutputArray filteredParticles
    424426
    425427  set PTMin 0.0
    426  
     428
    427429  add PdgCode {12}
    428430  add PdgCode {14}
     
    475477  set OutputArray jets
    476478
    477   # scale formula for jets 
     479  # scale formula for jets
    478480  set ScaleFormula {  sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 )  }
    479481}
     
    503505
    504506  # efficiency formula for c-jets (misidentification rate)
    505   add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
    506                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
    507                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
     507  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) +
     508                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
     509                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
    508510                              (abs(eta) > 2.5)                                  * (0.000)}
    509511
    510512  # efficiency formula for b-jets
    511   add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
    512                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
    513                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
     513  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) +
     514                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
     515                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
    514516                              (abs(eta) > 2.5)                                  * (0.000)}
    515517}
  • cards/delphes_card_ATLAS_PileUp.tcl

    r4e09c3a ra2983ec  
    7979  set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
    8080
    81   #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
    82   #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
    83   #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) + \
     81  #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
     82  #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) +
     83  #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) +
    8484  #                               (abs(t) >  1.0e-09) * (abs(z) > 0.15)  * (0.00)}
    8585
     
    119119
    120120  # tracking efficiency formula for charged hadrons
    121   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    122                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    123                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    124                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    125                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
     121  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     122                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     123                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     124                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     125                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    126126                         (abs(eta) > 2.5)                                                  * (0.00)}
    127127}
     
    138138
    139139  # tracking efficiency formula for electrons
    140   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    141                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
    142                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
    143                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
    144                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
    145                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
    146                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
     140  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     141                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     142                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     143                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     144                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     145                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     146                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    147147                         (abs(eta) > 2.5)                                                  * (0.00)}
    148148}
     
    159159
    160160  # tracking efficiency formula for muons
    161   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    162                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    163                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    164                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    165                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
     161  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     162                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     163                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     164                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     165                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    166166                         (abs(eta) > 2.5)                                                  * (0.00)}
    167167}
     
    178178
    179179  # resolution formula for charged hadrons
    180   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    181                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    182                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    183                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    184                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    185                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    186                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
     180  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     181                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     182                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     183                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     184                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     185                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     186                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    187187                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    188188}
     
    198198  # set ResolutionFormula {resolution formula as a function of eta and energy}
    199199
    200   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
    201                                            (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
    202                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
     200  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) +
     201                                           (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
     202                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
    203203                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    204204
     
    216216
    217217  # resolution formula for muons
    218   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    219                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) + \
    220                                            (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
    221                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) + \
    222                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
    223                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) + \
    224                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
     218  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     219                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) +
     220                                           (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) +
     221                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) +
     222                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) +
     223                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) +
     224                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) +
    225225                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.10)}
    226226}
     
    249249  set PhotonOutputArray photons
    250250
    251   set ECalEnergyMin 0.50
    252   set HCalEnergyMin 1.00
    253 
    254   set ECalEnergySignificanceMin 1.0
    255   set HCalEnergySignificanceMin 1.0
    256 
    257251  set EFlowTrackOutputArray eflowTracks
    258252  set EFlowPhotonOutputArray eflowPhotons
    259253  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     254
     255  set ECalEnergyMin 0.5
     256  set HCalEnergyMin 1.0
     257
     258  set ECalEnergySignificanceMin 1.0
     259  set HCalEnergySignificanceMin 1.0
     260
     261  set SmearTowerCenter true
    260262
    261263  set pi [expr {acos(-1)}]
     
    307309  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    308310  # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
    309   set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
     311  set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
    310312                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
    311313
     
    313315  # http://arxiv.org/pdf/hep-ex/0004009v1
    314316  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    315   set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
    316                              (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
     317  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
     318                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
    317319                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
    318320}
     
    482484
    483485  # efficiency formula for photons
    484   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    485                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    486                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     486  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     487                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     488                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    487489                         (abs(eta) > 2.5)                                   * (0.00)}
    488490}
     
    517519
    518520  # efficiency formula for electrons
    519   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    520                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    521                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     521  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     522                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     523                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    522524                         (abs(eta) > 2.5)                                   * (0.00)}
    523525}
     
    552554
    553555  # efficiency formula for muons
    554   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    555                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    556                          (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) + \
     556  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     557                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     558                         (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) +
    557559                         (abs(eta) > 2.7)                                   * (0.00)}
    558560}
     
    624626
    625627  # efficiency formula for c-jets (misidentification rate)
    626   add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
    627                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
    628                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
     628  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) +
     629                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
     630                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
    629631                              (abs(eta) > 2.5)                                  * (0.000)}
    630632
    631633  # efficiency formula for b-jets
    632   add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
    633                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
    634                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
     634  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) +
     635                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
     636                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
    635637                              (abs(eta) > 2.5)                                  * (0.000)}
    636638}
     
    698700
    699701}
    700 
    701 
    702 
  • cards/delphes_card_CMS.tcl

    r4e09c3a ra2983ec  
    7777
    7878  # tracking efficiency formula for charged hadrons
    79   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    80                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    81                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    82                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    83                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
     79  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     80                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     81                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     82                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     83                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    8484                         (abs(eta) > 2.5)                                                  * (0.00)}
    8585}
     
    9696
    9797  # tracking efficiency formula for electrons
    98   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    99                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
    100                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
    101                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
    102                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
    103                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
    104                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
     98  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     99                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     100                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     101                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     102                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     103                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     104                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    105105                         (abs(eta) > 2.5)                                                  * (0.00)}
    106106}
     
    117117
    118118  # tracking efficiency formula for muons
    119   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    120                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    121                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    122                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    123                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
     119  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     120                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     121                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     122                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     123                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    124124                         (abs(eta) > 2.5)                                                  * (0.00)}
    125125}
     
    136136
    137137  # resolution formula for charged hadrons
    138   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    139                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    140                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    141                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    142                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    143                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    144                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
     138  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     139                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     140                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     141                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     142                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     143                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     144                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    145145                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    146146}
     
    157157
    158158  # resolution formula for electrons
    159   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) + \
    160                                            (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    161                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
     159  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
     160                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
     161                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    162162                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    163163
     
    175175
    176176  # resolution formula for muons
    177   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
    178                                            (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
    179                                            (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
    180                                            (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    181                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
    182                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
    183                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
    184                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    185                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
    186                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
    187                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
     177  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
     178                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
     179                                           (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
     180                                           (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     181                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
     182                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
     183                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
     184                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     185                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
     186                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
     187                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    188188                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
    189189}
     
    212212  set PhotonOutputArray photons
    213213
    214   set ECalEnergyMin 0.5
    215   set HCalEnergyMin 1.0
    216 
    217   set ECalEnergySignificanceMin 1.0
    218   set HCalEnergySignificanceMin 1.0
    219    
    220214  set EFlowTrackOutputArray eflowTracks
    221215  set EFlowPhotonOutputArray eflowPhotons
    222216  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     217
     218  set ECalEnergyMin 0.5
     219  set HCalEnergyMin 1.0
     220
     221  set ECalEnergySignificanceMin 1.0
     222  set HCalEnergySignificanceMin 1.0
     223
     224  set SmearTowerCenter true
    223225
    224226  set pi [expr {acos(-1)}]
     
    276278
    277279  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    278   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  + \
     280  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  +
    279281                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    280282
    281283  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    282   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
     284  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    283285                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    284286}
     
    307309
    308310  # efficiency formula for photons
    309   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    310                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    311                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     311  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     312                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     313                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    312314                         (abs(eta) > 2.5)                                   * (0.00)}
    313315}
     
    341343
    342344  # efficiency formula for electrons
    343   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    344                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    345                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     345  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     346                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     347                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    346348                         (abs(eta) > 2.5)                                   * (0.00)}
    347349}
     
    375377
    376378  # efficiency formula for muons
    377   set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) + \
    378                                            (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
    379                                            (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
    380                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
    381                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
     379  set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
     380                                           (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     381                                           (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
     382                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     383                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
    382384                         (abs(eta) > 2.4)                                                 * (0.00)}
    383385}
     
    429431
    430432module PdgCodeFilter NeutrinoFilter {
    431  
     433
    432434  set InputArray Delphes/stableParticles
    433435  set OutputArray filteredParticles
    434436
    435437  set PTMin 0.0
    436  
     438
    437439  add PdgCode {12}
    438440  add PdgCode {14}
     
    518520
    519521  # efficiency formula for c-jets (misidentification rate)
    520   add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
    521                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
    522                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
     522  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) +
     523                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
     524                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
    523525                              (abs(eta) > 2.5)                                  * (0.000)}
    524526
    525527  # efficiency formula for b-jets
    526   add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
    527                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
    528                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
     528  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) +
     529                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
     530                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
    529531                              (abs(eta) > 2.5)                                  * (0.000)}
    530532}
     
    589591  add Branch ScalarHT/energy ScalarHT ScalarHT
    590592}
    591 
  • cards/delphes_card_CMS_NoFastJet.tcl

    r4e09c3a ra2983ec  
    5151
    5252  # tracking efficiency formula for charged hadrons
    53   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    54                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    55                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    56                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    57                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
     53  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     54                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     55                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     56                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     57                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    5858                         (abs(eta) > 2.5)                                                  * (0.00)}
    5959}
     
    7070
    7171  # tracking efficiency formula for electrons
    72   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    73                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
    74                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
    75                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
    76                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
    77                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
    78                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
     72  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     73                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     74                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     75                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     76                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     77                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     78                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    7979                         (abs(eta) > 2.5)                                                  * (0.00)}
    8080}
     
    9191
    9292  # tracking efficiency formula for muons
    93   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    94                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    95                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    96                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    97                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
     93  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     94                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     95                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     96                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     97                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    9898                         (abs(eta) > 2.5)                                                  * (0.00)}
    9999}
     
    110110
    111111  # resolution formula for charged hadrons
    112   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    113                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    114                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    115                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    116                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    117                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    118                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
     112  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     113                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     114                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     115                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     116                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     117                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     118                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    119119                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    120120}
     
    131131
    132132  # resolution formula for electrons
    133   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) + \
    134                                            (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    135                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
     133  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
     134                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
     135                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    136136                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    137137
     
    149149
    150150  # resolution formula for muons
    151   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
    152                                            (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
    153                                            (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
    154                                            (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    155                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
    156                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
    157                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
    158                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    159                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
    160                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
    161                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
     151  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
     152                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
     153                                           (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
     154                                           (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     155                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
     156                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
     157                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
     158                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     159                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
     160                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
     161                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    162162                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
    163163}
     
    186186  set PhotonOutputArray photons
    187187
    188   set ECalEnergyMin 0.5
    189   set HCalEnergyMin 1.0
    190 
    191   set ECalEnergySignificanceMin 1.0
    192   set HCalEnergySignificanceMin 1.0
    193    
    194188  set EFlowTrackOutputArray eflowTracks
    195189  set EFlowPhotonOutputArray eflowPhotons
    196190  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     191
     192  set ECalEnergyMin 0.5
     193  set HCalEnergyMin 1.0
     194
     195  set ECalEnergySignificanceMin 1.0
     196  set HCalEnergySignificanceMin 1.0
     197
     198  set SmearTowerCenter true
    197199
    198200  set pi [expr {acos(-1)}]
     
    250252
    251253  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    252   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  + \
     254  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  +
    253255                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    254256
    255257  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    256   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
     258  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    257259                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    258260}
     
    269271  set OutputArray eflow
    270272}
    271 
    272 
  • cards/delphes_card_CMS_PileUp.tcl

    r4e09c3a ra2983ec  
    8080  set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
    8181
    82   #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
    83   #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
    84   #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) + \
     82  #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
     83  #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) +
     84  #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) +
    8585  #                               (abs(t) >  1.0e-09) * (abs(z) > 0.15)  * (0.00)}
    8686
     
    120120
    121121  # tracking efficiency formula for charged hadrons
    122   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    123                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    124                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    125                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    126                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
     122  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     123                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     124                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     125                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     126                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    127127                         (abs(eta) > 2.5)                                                  * (0.00)}
    128128}
     
    139139
    140140  # tracking efficiency formula for electrons
    141   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    142                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
    143                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
    144                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
    145                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
    146                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
    147                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
     141  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     142                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     143                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     144                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     145                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     146                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     147                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    148148                         (abs(eta) > 2.5)                                                  * (0.00)}
    149149}
     
    160160
    161161  # tracking efficiency formula for muons
    162   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    163                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    164                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    165                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    166                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
     162  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     163                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     164                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     165                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     166                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    167167                         (abs(eta) > 2.5)                                                  * (0.00)}
    168168}
     
    179179
    180180  # resolution formula for charged hadrons
    181   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    182                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    183                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    184                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    185                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    186                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    187                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
     181  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     182                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     183                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     184                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     185                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     186                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     187                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    188188                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    189189}
     
    200200
    201201  # resolution formula for electrons
    202   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) + \
    203                                            (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    204                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
     202  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
     203                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
     204                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    205205                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    206206
     
    218218
    219219  # resolution formula for muons
    220   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
    221                                            (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
    222                                            (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
    223                                            (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    224                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
    225                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
    226                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
    227                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    228                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
    229                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
    230                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
     220  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
     221                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
     222                                           (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
     223                                           (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     224                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
     225                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
     226                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
     227                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     228                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
     229                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
     230                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    231231                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
    232232}
     
    255255  set PhotonOutputArray photons
    256256
    257   set ECalEnergyMin 0.5
    258   set HCalEnergyMin 1.0
    259 
    260   set ECalEnergySignificanceMin 1.0
    261   set HCalEnergySignificanceMin 1.0
    262 
    263257  set EFlowTrackOutputArray eflowTracks
    264258  set EFlowPhotonOutputArray eflowPhotons
    265259  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     260
     261  set ECalEnergyMin 0.5
     262  set HCalEnergyMin 1.0
     263
     264  set ECalEnergySignificanceMin 1.0
     265  set HCalEnergySignificanceMin 1.0
     266
     267  set SmearTowerCenter true
    266268
    267269  set pi [expr {acos(-1)}]
     
    319321
    320322  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    321   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  + \
     323  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    322324                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    323325
    324326  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    325   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
     327  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    326328                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    327329}
     
    518520
    519521  # efficiency formula for photons
    520   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    521                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    522                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     522  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     523                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     524                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    523525                         (abs(eta) > 2.5)                                   * (0.00)}
    524526}
     
    554556
    555557  # efficiency formula for electrons
    556   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    557                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    558                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     558  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     559                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     560                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    559561                         (abs(eta) > 2.5)                                   * (0.00)}
    560562}
     
    589591
    590592  # efficiency formula for muons
    591   set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) + \
    592                                            (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
    593                                            (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
    594                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
    595                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
     593  set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
     594                                           (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     595                                           (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
     596                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     597                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
    596598                         (abs(eta) > 2.4)                                                 * (0.00)}
    597599}
     
    665667
    666668  # efficiency formula for c-jets (misidentification rate)
    667   add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
    668                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
    669                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
     669  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) +
     670                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
     671                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
    670672                              (abs(eta) > 2.5)                                  * (0.000)}
    671673
    672674  # efficiency formula for b-jets
    673   add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
    674                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
    675                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
     675  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) +
     676                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
     677                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
    676678                              (abs(eta) > 2.5)                                  * (0.000)}
    677679}
     
    738740  add Branch PileUpMerger/vertices Vertex Vertex
    739741}
    740 
  • cards/delphes_card_FCC_basic.tcl

    r4e09c3a ra2983ec  
    1717  AngularSmearing
    1818  ImpactParameterSmearing
    19  
     19
    2020  ECal
    2121  HCal
    22  
     22
    2323  TowerMerger
    2424  EFlowMerger
    2525
    2626  MissingET
    27  
     27
    2828  NeutrinoFilter
    2929  GenJetFinder
     
    3131
    3232  JetEnergyScale
    33  
     33
    3434  TrackCountingBTagging
    3535  TauTagging
     
    7272
    7373  # tracking efficiency formula for charged hadrons
    74   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    75                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    76                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    77                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    78                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.90) + \
     74  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     75                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     76                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     77                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     78                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.90) +
    7979                         (abs(eta) > 4.0)                                                  * (0.00)}
    80                          
     80
    8181}
    8282
     
    9292
    9393  # tracking efficiency formula for electrons
    94   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    95                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    96                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    97                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    98                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) + \
     94  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     95                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     96                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     97                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     98                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) +
    9999                         (abs(eta) > 4.0)                                                  * (0.00)}
    100100}
     
    111111
    112112  # tracking efficiency formula for muons
    113   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    114                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    115                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    116                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    117                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) + \
     113  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     114                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     115                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     116                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     117                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) +
    118118                         (abs(eta) > 4.0)                                                  * (0.00)}
    119119}
     
    130130
    131131  # resolution formula for charged hadrons
    132   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) + \
     132  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) +
    133133                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)}
    134134
     
    147147
    148148  # resolution formula for electrons
    149   set ResolutionFormula {                  (abs(eta) <= 4.0) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.007) + \
    150                                            (abs(eta) <= 4.0) * (energy > 2.0e1)                    * sqrt(energy^2*0.005^2 + energy*0.02^2) + \
     149  set ResolutionFormula {                  (abs(eta) <= 4.0) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.007) +
     150                                           (abs(eta) <= 4.0) * (energy > 2.0e1)                    * sqrt(energy^2*0.005^2 + energy*0.02^2) +
    151151                                           (abs(eta) > 4.0 && abs(eta) <= 6.0)                     * sqrt(energy^2*0.05^2 + energy*1.00^2)}
    152152
     
    164164
    165165  # resolution formula for muons
    166   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*5.e-6) + \
     166  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*5.e-6) +
    167167                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.e-5)}
    168168
     
    210210
    211211  # absolute impact parameter smearing formula (in mm) as a function of pt and eta
    212   set ResolutionFormula {(pt > 0.1  && pt <= 5.0)   * (0.010) + \
     212  set ResolutionFormula {(pt > 0.1  && pt <= 5.0)   * (0.010) +
    213213                         (pt > 5.0)                 * (0.005)}
    214214
     
    225225  set TowerOutputArray ecalTowers
    226226  set EFlowTowerOutputArray eflowPhotons
    227  
     227
    228228  set EnergyMin 0.5
    229229  set EnergySignificanceMin 1.0
    230  
     230
     231  set SmearTowerCenter true
     232
    231233  set pi [expr {acos(-1)}]
    232234
     
    234236  # each list starts with the lower edge of the first tower
    235237  # the list ends with the higher edged of the last tower
    236    
     238
    237239  # 0.5 degree towers
    238240  set PhiBins {}
     
    246248    add EtaPhiBins $eta $PhiBins
    247249  }
    248  
     250
    249251  # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
    250  
     252
    251253  add EnergyFraction {0} {0.0}
    252254  # energy fractions for e, gamma and pi0
     
    269271
    270272  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    271  
     273
    272274  # This is the CMS ECAL resolution, extended up eta = 6.0
    273   set ResolutionFormula { (abs(eta) <= 3.0)                   * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2)  + \
     275  set ResolutionFormula { (abs(eta) <= 3.0)                   * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2)  +
    274276                          (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    275277
    276  
     278
    277279}
    278280
     
    287289  set TowerOutputArray hcalTowers
    288290  set EFlowTowerOutputArray eflowNeutralHadrons
    289  
     291
    290292  set EnergyMin 1.0
    291293  set EnergySignificanceMin 1.0
    292  
     294
     295  set SmearTowerCenter true
     296
    293297  set pi [expr {acos(-1)}]
    294298
     
    296300  # each list starts with the lower edge of the first tower
    297301  # the list ends with the higher edged of the last tower
    298  
    299  
     302
     303
    300304  # 5 degree towers
    301305  set PhiBins {}
     
    309313    add EtaPhiBins $eta $PhiBins
    310314  }
    311  
    312    
     315
     316
    313317  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    314318  add EnergyFraction {0} {1.0}
     
    332336
    333337  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    334  
     338
    335339  # This is the ATLAS HCAL resolution, extended up eta = 6.0
    336   set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
    337                              (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
     340  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
     341                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
    338342                             (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
    339343
     
    373377  set MomentumOutputArray momentum
    374378}
    375                    
     379
    376380
    377381##################
     
    390394
    391395module PdgCodeFilter NeutrinoFilter {
    392  
     396
    393397  set InputArray Delphes/stableParticles
    394398  set OutputArray filteredParticles
    395399
    396400  set PTMin 0.0
    397  
     401
    398402  add PdgCode {12}
    399403  add PdgCode {14}
     
    463467  # maximum distance between jet and track
    464468  set DeltaR 0.3
    465  
    466   # minimum pt of tracks 
     469
     470  # minimum pt of tracks
    467471  set TrackPTMin 1.0
    468  
     472
    469473  # minimum transverse impact parameter (in mm)
    470474  set TrackIPMax 2.0
    471  
     475
    472476  # minimum ip significance for the track to be counted
    473477  set SigMin 6.5
    474  
     478
    475479  # minimum number of tracks (high efficiency n=2, high purity n=3)
    476480  set Ntracks 3
     
    510514  add Branch Delphes/allParticles Particle GenParticle
    511515  add Branch GenJetFinder/jets GenJet Jet
    512  
     516
    513517  add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
    514518  add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
     
    521525  add Branch ScalarHT/energy ScalarHT ScalarHT
    522526}
    523  
  • cards/delphes_card_LHCb.tcl

    r4e09c3a ra2983ec  
    88
    99set ExecutionPath {
    10  
     10
    1111  ParticlePropagator
    12  
     12
    1313  ChargedHadronMomentumSmearing
    1414  ElectronEnergySmearing
     
    1717  TrackMerger
    1818  ImpactParameterSmearing
    19  
     19
    2020  IdentificationMap
    21  
     21
    2222  ECal
    2323  HCal
     
    4141
    4242  # radius of the magnetic field coverage, in m
    43  
     43
    4444  set Radius 3.31
    45  
    46  
     45
     46
    4747  # half-length of the magnetic field coverage, in m
    4848  set HalfLength 12.0
     
    5353  # Need to veto anything with theta > 0.269 rad  -> eta = 2
    5454  #                            theta < 0.0135 rad -> eta = 5
    55  
    56   # tracker and calos are at approx 0.269 rad, R = 12*tan(0.269) 
    57 
    58 }
    59 
    60  
     55
     56  # tracker and calos are at approx 0.269 rad, R = 12*tan(0.269)
     57
     58}
     59
     60
    6161########################################
    6262# Momentum resolution for charged tracks
     
    8484
    8585  # resolution formula for electrons
    86   set ResolutionFormula { (eta > 2.0  && eta <= 5.0) * (energy > 0.1   && energy <= 8.0) * (energy*0.05) + \
     86  set ResolutionFormula { (eta > 2.0  && eta <= 5.0) * (energy > 0.1   && energy <= 8.0) * (energy*0.05) +
    8787                          (eta > 2.0  && eta <= 5.0) * (energy > 8.0)                    *  sqrt(energy^2*0.015^2 + energy*0.10^2)}
    8888  }
     
    139139  set InputArray ImpactParameterSmearing/tracks
    140140  set OutputArray tracks
    141  
     141
    142142  # {PID in} {PID out} {formula}
    143143  # make sure "PID in" and "PID out" have the same charge (e.g {-13} {211} or {-321} {211})
    144144  # {211} {-13} is equivalent to {-211} {13} (and needs to be written once only...)
    145  
    146  
    147  
    148  
    149  
    150  
    151   # --- pions --- 
    152  
    153   add EfficiencyFormula {211} {211} {      (eta <= 2.0)                                  * (0.00) + \
    154                                            (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) + \
    155                                            (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) + \
     145
     146
     147
     148
     149
     150
     151  # --- pions ---
     152
     153  add EfficiencyFormula {211} {211} {      (eta <= 2.0)                                  * (0.00) +
     154                                           (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) +
     155                                           (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) +
    156156                                           (eta > 5.0)                                   * (0.00)}
    157  
    158   add EfficiencyFormula {211} {-13} {      (eta <= 2.0)                                 * (0.00) + \
    159                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) + \
    160                                            (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.0663*exp(-0.13*pt*cosh(eta))) + \
     157
     158  add EfficiencyFormula {211} {-13} {      (eta <= 2.0)                                 * (0.00) +
     159                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) +
     160                                           (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.0663*exp(-0.13*pt*cosh(eta))) +
    161161                                           (eta > 5.0)                                  * (0.00)}
    162162
    163  
     163
    164164 # --- kaons ---
    165  
    166  
    167   add EfficiencyFormula {321} {321} {      (eta <= 2.0)                                  * (0.00) + \
    168                                            (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) + \
    169                                            (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) + \
     165
     166
     167  add EfficiencyFormula {321} {321} {      (eta <= 2.0)                                  * (0.00) +
     168                                           (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) +
     169                                           (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) +
    170170                                           (eta > 5.0)                                   * (0.00)}
    171                                            
    172   add EfficiencyFormula {321} {-13} {      (eta <= 2.0)                                 * (0.00) + \
    173                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) + \
    174                                            (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.086*exp(-0.11*pt*cosh(eta))) + \
     171
     172  add EfficiencyFormula {321} {-13} {      (eta <= 2.0)                                 * (0.00) +
     173                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) +
     174                                           (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.086*exp(-0.11*pt*cosh(eta))) +
    175175                                           (eta > 5.0)                                  * (0.00)}
    176176
    177177
    178178 # --- protons ---
    179  
    180  
    181   add EfficiencyFormula {2212} {2212} {    (eta <= 2.0)                                  * (0.00) + \
    182                                            (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) + \
    183                                            (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) + \
     179
     180
     181  add EfficiencyFormula {2212} {2212} {    (eta <= 2.0)                                  * (0.00) +
     182                                           (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) +
     183                                           (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) +
    184184                                           (eta > 5.0)                                   * (0.00)}
    185                                            
    186   add EfficiencyFormula {2212} {-13} {     (eta <= 2.0)                                 * (0.00) + \
    187                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) + \
    188                                            (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.002) + \
     185
     186  add EfficiencyFormula {2212} {-13} {     (eta <= 2.0)                                 * (0.00) +
     187                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) +
     188                                           (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.002) +
    189189                                           (eta > 5.0)                                  * (0.00)}
    190190
     
    195195
    196196
    197   add EfficiencyFormula {-13} {-13} {      (eta <= 2.0)                                * (0.00) + \
    198                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8)* (0.00) + \
    199                                            (eta > 2.0  && eta <= 5.0) *     (pt >= 0.8)* (0.97) + \
     197  add EfficiencyFormula {-13} {-13} {      (eta <= 2.0)                                * (0.00) +
     198                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8)* (0.00) +
     199                                           (eta > 2.0  && eta <= 5.0) *     (pt >= 0.8)* (0.97) +
    200200                                           (eta > 5.0)                                 * (0.00)}
    201201
    202  
     202
    203203 # efficiency for other charged particles (should be always {0} {0} {formula})
    204  
    205   add EfficiencyFormula {0} {0}     {      (eta <= 2.0)                                * (0.00) + \
    206                                            (eta > 2.0  && eta <= 5.0) *     (pt < 0.8) * (0.00) + \
    207                                            (eta > 2.0  && eta <= 5.0) *     (pt > 0.8) * (0.95) + \
     204
     205  add EfficiencyFormula {0} {0}     {      (eta <= 2.0)                                * (0.00) +
     206                                           (eta > 2.0  && eta <= 5.0) *     (pt < 0.8) * (0.00) +
     207                                           (eta > 2.0  && eta <= 5.0) *     (pt > 0.8) * (0.95) +
    208208                                           (eta > 5.0)                                 * (0.00)}
    209209
     
    223223  set TowerOutputArray ecalTowers
    224224  set EFlowTowerOutputArray eflowPhotons
    225  
     225
     226  set EnergyMin 0.0
     227  set EnergySignificanceMin 0.0
     228
     229  set SmearTowerCenter true
     230
    226231  set pi [expr {acos(-1)}]
    227232
     
    241246    add EtaPhiBins $eta $PhiBins
    242247  }
    243  
     248
    244249    # 1.25 degree towers
    245250  set PhiBins {}
     
    285290  add EnergyFraction {3122} {0.3}
    286291
    287    set ResolutionFormula {(eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)}
     292  set ResolutionFormula {(eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)}
    288293}
    289294
     
    298303  set TowerOutputArray hcalTowers
    299304  set EFlowTowerOutputArray eflowNeutralHadrons
    300  
     305
     306  set EnergyMin 0.0
     307  set EnergySignificanceMin 0.0
     308
     309  set SmearTowerCenter true
     310
    301311  set pi [expr {acos(-1)}]
    302312
     
    330340
    331341
    332    
     342
    333343  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    334344  add EnergyFraction {0} {1.0}
     
    352362
    353363  set ResolutionFormula { (eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.05^2 + energy*0.80^2)}
    354   }
     364}
    355365
    356366
     
    360370
    361371# tracks, towers and eflow objects are not stored by default in the output.
    362 # if needed (for jet constituent or other studies), uncomment the relevant 
    363 # "add Branch ..." lines. 
     372# if needed (for jet constituent or other studies), uncomment the relevant
     373# "add Branch ..." lines.
    364374
    365375module TreeWriter TreeWriter {
    366376# add Branch InputArray BranchName BranchClass
    367  
     377
    368378  add Branch Delphes/allParticles Particle GenParticle
    369379
     
    371381  add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
    372382  add Branch ECal/eflowPhotons Photon Photon
    373  
    374 }
    375 
     383
     384}
     385
Note: See TracChangeset for help on using the changeset viewer.