Fork me on GitHub

Changeset 221abf3 in git for cards/delphes_card_FCC.tcl


Ignore:
Timestamp:
Oct 9, 2015, 12:05:39 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
679a87e
Parents:
233a019
Message:

update values for FCC, move back to Calorimeter

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_FCC.tcl

    r233a019 r221abf3  
    1515
    1616  TrackMerger
    17   AngularSmearing
    18   ImpactParameterSmearing
    19 
    20   ECal
    21   HCal
     17  Calorimeter
     18  EFlowMerger
     19
     20  PhotonEfficiency
     21  PhotonIsolation
    2222
    2323  ElectronFilter
    24 
    25   TowerMerger
    26   EFlowMerger
     24  ElectronEfficiency
     25  ElectronIsolation
     26
     27  MuonEfficiency
     28  MuonIsolation
    2729
    2830  MissingET
    … …  
    3436  JetEnergyScale
    3537
    36   TrackCountingBTagging
     38  JetFlavorAssociation
     39
     40  BTagging
    3741  TauTagging
     42
     43  UniqueObjectFinder
    3844
    3945  ScalarHT
    … …  
    5561
    5662  # radius of the magnetic field coverage, in m
    57   set Radius 2.60
     63  set Radius 6.00
    5864  # half-length of the magnetic field coverage, in m
    59   set HalfLength 6.00
     65  set HalfLength 12.00
    6066
    6167  # magnetic field
    … …  
    7480
    7581  # tracking efficiency formula for charged hadrons
    76   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    77                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     82  set EfficiencyFormula {                                                    (pt <= 0.5)   * (0.00) +
     83                                           (abs(eta) <= 1.5) * (pt > 0.5   && pt <= 1.0)   * (0.85) +
    7884                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
    79                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     85                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.5   && pt <= 1.0)   * (0.80) +
    8086                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.90) +
    81                          (abs(eta) > 4.0)                                                  * (0.00)}
    82 
     87                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.5   && pt <= 1.0)   * (0.75) +
     88                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0)                  * (0.80) +
     89                         (abs(eta) > 6.0)                                                  * (0.00)
     90                        }
    8391}
    8492
    … …  
    94102
    95103  # tracking efficiency formula for electrons
    96   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    97                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
    98                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
    99                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
    100                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) +
    101                          (abs(eta) > 4.0)                                                  * (0.00)}
     104 set EfficiencyFormula {                                                    (pt <= 0.5)    * (0.00) +
     105                                           (abs(eta) <= 1.5) * (pt > 0.5   && pt <= 1.0)   * (0.80) +
     106                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.90) +
     107                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.5   && pt <= 1.0)   * (0.75) +
     108                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.85) +
     109                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.5   && pt <= 1.0)   * (0.70) +
     110                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0)                  * (0.80) +
     111                         (abs(eta) > 6.0)                                                  * (0.00)
     112                       }
    102113}
    103114
    … …  
    111122
    112123  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    113 
    114   # tracking efficiency formula for muons
    115   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    116                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     124  set EfficiencyFormula {                                                   (pt <= 0.5)    * (0.00) +
     125                                           (abs(eta) <= 1.5) * (pt > 0.5   && pt <= 1.0)   * (0.90) +
    117126                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
    118                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
    119                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) +
    120                          (abs(eta) > 4.0)                                                  * (0.00)}
     127                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.5   && pt <= 1.0)   * (0.85) +
     128                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.95) +
     129                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.5   && pt <= 1.0)   * (0.80) +
     130                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0)                  * (0.90) +
     131                         (abs(eta) > 6.0)                                                  * (0.00)
     132                       }
     133 
    121134}
    122135
    … …  
    133146  # resolution formula for charged hadrons
    134147  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) +
    135                              (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)}
     148                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5) +
     149                             (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * (0.03 + pt*5.e-5)}
    136150
    137151
    … …  
    150164  # resolution formula for electrons
    151165  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) +
    152                              (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)}
     166                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5) +
     167                             (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * (0.03 + pt*5.e-5)}
     168
    153169}
    154170
    … …  
    164180
    165181  # resolution formula for muons
    166   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*5.e-6) +
    167                              (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.e-5)}
     182  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*5.0e-6) +
     183                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.0e-5) +
     184                             (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 0.1) * (0.03 + pt*1.5e-5)}
    168185
    169186
    … …  
    183200
    184201
    185 ########################
    186 # Track angular smearing
    187 ########################
    188 
    189 module AngularSmearing AngularSmearing {
    190   set InputArray TrackMerger/tracks
    191   set OutputArray tracks
    192 
    193 
    194   # angular smearing  in eta formula as a function of pt and eta
    195   set EtaResolutionFormula { 0.001 }
    196 
    197   # angular smearing  in phi formula as a function of pt and eta
    198   set PhiResolutionFormula { 0.001 }
    199 
    200 }
    201 
    202 #################################
    203 # Track impact parameter smearing
    204 #################################
    205 
    206 module ImpactParameterSmearing ImpactParameterSmearing {
    207   set InputArray AngularSmearing/tracks
    208   set OutputArray tracks
    209 
    210 
    211   # 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) +
    213                          (pt > 5.0)                 * (0.005)}
    214 
    215 }
    216 
    217 #############
    218 #   ECAL
    219 #############
    220 
    221 module SimpleCalorimeter ECal {
     202#################
     203# Calorimeter
     204#################
     205
     206module Calorimeter Calorimeter {
    222207  set ParticleInputArray ParticlePropagator/stableParticles
    223   set TrackInputArray ImpactParameterSmearing/tracks
    224 
    225   set TowerOutputArray ecalTowers
     208  set TrackInputArray TrackMerger/tracks
     209
     210  set TowerOutputArray towers
     211  set PhotonOutputArray photons
     212
    226213  set EFlowTrackOutputArray eflowTracks
    227   set EFlowTowerOutputArray eflowPhotons
    228 
    229   set IsEcal true
    230  
    231   set EnergyMin 0.5
    232   set EnergySignificanceMin 1.0
     214  set EFlowPhotonOutputArray eflowPhotons
     215  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     216
     217  set ECalEnergyMin 0.5
     218  set HCalEnergyMin 1.0
     219
     220  set ECalEnergySignificanceMin 1.0
     221  set HCalEnergySignificanceMin 1.0
    233222
    234223  set SmearTowerCenter true
    … …  
    246235  }
    247236
    248   # 0.01 unit in eta up to eta = 3
    249   for {set i -300} {$i <= 300} {incr i} {
     237  # 0.01 unit in eta up to eta = 6
     238  for {set i -600} {$i <= 600} {incr i} {
    250239    set eta [expr {$i * 0.01}]
    251240    add EtaPhiBins $eta $PhiBins
    … …  
    254243  # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
    255244
    256   add EnergyFraction {0} {0.0}
     245  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
     246  add EnergyFraction {0} {0.0 1.0}
    257247  # energy fractions for e, gamma and pi0
    258   add EnergyFraction {11} {1.0}
    259   add EnergyFraction {22} {1.0}
    260   add EnergyFraction {111} {1.0}
     248  add EnergyFraction {11} {1.0 0.0}
     249  add EnergyFraction {22} {1.0 0.0}
     250  add EnergyFraction {111} {1.0 0.0}
    261251  # energy fractions for muon, neutrinos and neutralinos
    262   add EnergyFraction {12} {0.0}
    263   add EnergyFraction {13} {0.0}
    264   add EnergyFraction {14} {0.0}
    265   add EnergyFraction {16} {0.0}
    266   add EnergyFraction {1000022} {0.0}
    267   add EnergyFraction {1000023} {0.0}
    268   add EnergyFraction {1000025} {0.0}
    269   add EnergyFraction {1000035} {0.0}
    270   add EnergyFraction {1000045} {0.0}
     252  add EnergyFraction {12} {0.0 0.0}
     253  add EnergyFraction {13} {0.0 0.0}
     254  add EnergyFraction {14} {0.0 0.0}
     255  add EnergyFraction {16} {0.0 0.0}
     256  add EnergyFraction {1000022} {0.0 0.0}
     257  add EnergyFraction {1000023} {0.0 0.0}
     258  add EnergyFraction {1000025} {0.0 0.0}
     259  add EnergyFraction {1000035} {0.0 0.0}
     260  add EnergyFraction {1000045} {0.0 0.0}
    271261  # energy fractions for K0short and Lambda
    272   add EnergyFraction {310} {0.3}
    273   add EnergyFraction {3122} {0.3}
    274 
    275   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    276 
    277   # This is the CMS ECAL resolution, extended up eta = 6.0
    278   set ResolutionFormula { (abs(eta) <= 3.0)                   * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2)  +
    279                           (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    280 
    281 
    282 }
    283 
    284 #############
    285 #   HCAL
    286 #############
    287 
    288 module SimpleCalorimeter HCal {
    289   set ParticleInputArray ParticlePropagator/stableParticles
    290   set TrackInputArray ECal/eflowTracks
    291 
    292   set TowerOutputArray hcalTowers
    293   set EFlowTrackOutputArray eflowTracks
    294   set EFlowTowerOutputArray eflowNeutralHadrons
    295 
    296   set IsEcal false
    297  
    298   set EnergyMin 1.0
    299   set EnergySignificanceMin 1.0
    300 
    301   set SmearTowerCenter true
    302 
    303   set pi [expr {acos(-1)}]
    304 
    305   # lists of the edges of each tower in eta and phi
    306   # each list starts with the lower edge of the first tower
    307   # the list ends with the higher edged of the last tower
    308 
    309 
    310   # 5 degree towers
    311   set PhiBins {}
    312   for {set i -72} {$i <= 72} {incr i} {
    313     add PhiBins [expr {$i * $pi/75.0}]
    314   }
    315 
    316   # 0.05 unit in eta up to eta = 3
    317   for {set i -60} {$i <= 60} {incr i} {
    318     set eta [expr {$i * 0.05}]
    319     add EtaPhiBins $eta $PhiBins
    320   }
    321 
    322 
    323   # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    324   add EnergyFraction {0} {1.0}
    325   # energy fractions for e, gamma and pi0
    326   add EnergyFraction {11} {0.0}
    327   add EnergyFraction {22} {0.0}
    328   add EnergyFraction {111} {0.0}
    329   # energy fractions for muon, neutrinos and neutralinos
    330   add EnergyFraction {12} {0.0}
    331   add EnergyFraction {13} {0.0}
    332   add EnergyFraction {14} {0.0}
    333   add EnergyFraction {16} {0.0}
    334   add EnergyFraction {1000022} {0.0}
    335   add EnergyFraction {1000023} {0.0}
    336   add EnergyFraction {1000025} {0.0}
    337   add EnergyFraction {1000035} {0.0}
    338   add EnergyFraction {1000045} {0.0}
    339   # energy fractions for K0short and Lambda
    340   add EnergyFraction {310} {0.7}
    341   add EnergyFraction {3122} {0.7}
    342 
    343   # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    344 
    345   # This is the ATLAS HCAL resolution, extended up eta = 6.0
    346   set ResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
    347                              (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
    348                              (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
    349 
     262  add EnergyFraction {310} {0.3 0.7}
     263  add EnergyFraction {3122} {0.3 0.7}
     264
     265
     266  set ECalResolutionFormula { (abs(eta) <= 3.0)                   * sqrt(energy^2*0.01^2 + energy*0.10^2)  +
     267                                    (abs(eta) > 3.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.10^2 + energy*2.00^2)}
     268
     269
     270  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.03^2 + energy*0.50^2) +
     271                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.05^2 + energy*0.70^2) +
     272                             (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09^2 + energy*1.00^2)}
     273
     274
     275}
     276
     277
     278####################
     279# Energy flow merger
     280####################
     281
     282module Merger EFlowMerger {
     283# add InputArray InputArray
     284  add InputArray Calorimeter/eflowTracks
     285  add InputArray Calorimeter/eflowPhotons
     286  add InputArray Calorimeter/eflowNeutralHadrons
     287  set OutputArray eflow
     288}
     289
     290###################
     291# Photon efficiency
     292###################
     293
     294module Efficiency PhotonEfficiency {
     295  set InputArray Calorimeter/eflowPhotons
     296  set OutputArray photons
     297
     298  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     299
     300  # efficiency formula for photons
     301  set EfficiencyFormula {                                      (pt <= 20.0) * (0.00) +
     302                                           (abs(eta) <= 1.5) * (pt > 20.0)  * (0.95) +
     303                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 20.0)  * (0.90) +
     304                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 20.0)  * (0.85) +
     305                         (abs(eta) > 6.0)                                   * (0.00)}
     306}
     307
     308##################
     309# Photon isolation
     310##################
     311
     312module Isolation PhotonIsolation {
     313  set CandidateInputArray PhotonEfficiency/photons
     314  set IsolationInputArray EFlowMerger/eflow
     315
     316  set OutputArray photons
     317
     318  set DeltaRMax 0.3
     319
     320  set PTMin 0.5
     321
     322  set PTRatioMax 0.1
    350323}
    351324
    … …  
    362335}
    363336
    364 ###################################################
    365 # Tower Merger (in case not using e-flow algorithm)
    366 ###################################################
    367 
    368 module Merger TowerMerger {
    369 # add InputArray InputArray
    370   add InputArray ECal/ecalTowers
    371   add InputArray HCal/hcalTowers
    372   set OutputArray towers
     337#####################
     338# Electron efficiency
     339#####################
     340
     341module Efficiency ElectronEfficiency {
     342  set InputArray ElectronFilter/electrons
     343  set OutputArray electrons
     344
     345  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     346
     347  # efficiency formula for electrons
     348  set EfficiencyFormula {                                      (pt <= 20.0) * (0.00) +
     349                                           (abs(eta) <= 1.5) * (pt > 20.0)  * (0.95) +
     350                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 20.0)  * (0.90) +
     351                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 20.0)  * (0.85) +
     352                         (abs(eta) > 6.0)                                   * (0.00)}
    373353}
    374354
    375355####################
    376 # Energy flow merger
     356# Electron isolation
    377357####################
    378358
    379 module Merger EFlowMerger {
    380 # add InputArray InputArray
    381   add InputArray HCal/eflowTracks
    382   add InputArray ECal/eflowPhotons
    383   add InputArray HCal/eflowNeutralHadrons
    384   set OutputArray eflow
    385 }
    386 
     359module Isolation ElectronIsolation {
     360  set CandidateInputArray ElectronEfficiency/electrons
     361  set IsolationInputArray EFlowMerger/eflow
     362
     363  set OutputArray electrons
     364
     365  set DeltaRMax 0.3
     366
     367  set PTMin 0.5
     368
     369  set PTRatioMax 0.1
     370}
     371
     372#################
     373# Muon efficiency
     374#################
     375
     376module Efficiency MuonEfficiency {
     377  set InputArray MuonMomentumSmearing/muons
     378  set OutputArray muons
     379
     380  # set EfficiencyFormula {efficiency as a function of eta and pt}
     381
     382  # efficiency formula for muons
     383  set EfficiencyFormula {                                      (pt <= 20.0) * (0.00) +
     384                                           (abs(eta) <= 1.5) * (pt > 20.0)  * (0.99) +
     385                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 20.0)  * (0.95) +
     386                         (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 20.0)  * (0.90) +
     387                         (abs(eta) > 6.0)                                   * (0.00)}
     388}
     389
     390################
     391# Muon isolation
     392################
     393
     394module Isolation MuonIsolation {
     395  set CandidateInputArray MuonEfficiency/muons
     396  set IsolationInputArray EFlowMerger/eflow
     397
     398  set OutputArray muons
     399
     400  set DeltaRMax 0.3
     401
     402  set PTMin 0.5
     403
     404  set PTRatioMax 0.1
     405}
    387406
    388407###################
    … …  
    396415}
    397416
    398 
    399417##################
    400418# Scalar HT merger
    … …  
    403421module Merger ScalarHT {
    404422# add InputArray InputArray
    405   add InputArray EFlowMerger/eflow
     423  add InputArray UniqueObjectFinder/jets
     424  add InputArray UniqueObjectFinder/electrons
     425  add InputArray UniqueObjectFinder/photons
     426  add InputArray UniqueObjectFinder/muons
    406427  set EnergyOutputArray energy
    407428}
    408429
    409 #################
     430
     431#####################
    410432# Neutrino Filter
    411 #################
     433#####################
    412434
    413435module PdgCodeFilter NeutrinoFilter {
    … …  
    441463  set ParameterR 0.5
    442464
    443   set JetPTMin 100.0
    444 }
     465  set JetPTMin 50.0
     466}
     467
    445468
    446469############
    … …  
    449472
    450473module FastJetFinder FastJetFinder {
    451 #  set InputArray TowerMerger/towers
     474#  set InputArray Calorimeter/towers
    452475  set InputArray EFlowMerger/eflow
    453476
    … …  
    458481  set ParameterR 0.5
    459482
    460   set JetPTMin 100.0
     483  set JetPTMin 50.0
    461484}
    462485
    … …  
    469492  set OutputArray jets
    470493
    471  # scale formula for jets
     494  # scale formula for jets
    472495  set ScaleFormula {1.00}
    473496}
    474497
    475 ##########################
    476 # Track Counting b-tagging
    477 ##########################
    478 
    479 module TrackCountingBTagging TrackCountingBTagging {
    480   set TrackInputArray ImpactParameterSmearing/tracks
     498########################
     499# Jet Flavor Association
     500########################
     501
     502module JetFlavorAssociation JetFlavorAssociation {
     503
     504  set PartonInputArray Delphes/partons
     505  set ParticleInputArray Delphes/allParticles
     506  set ParticleLHEFInputArray Delphes/allParticlesLHEF
    481507  set JetInputArray JetEnergyScale/jets
    482508
     509  set DeltaR 0.5
     510  set PartonPTMin 1.0
     511  set PartonEtaMax 6.0
     512
     513}
     514
     515###########
     516# b-tagging
     517###########
     518
     519module BTagging BTagging {
     520  set JetInputArray JetEnergyScale/jets
     521
    483522  set BitNumber 0
    484523
    485   # maximum distance between jet and track
    486   set DeltaR 0.3
    487 
    488   # minimum pt of tracks
    489   set TrackPTMin 1.0
    490 
    491   # minimum transverse impact parameter (in mm)
    492   set TrackIPMax 2.0
    493 
    494   # minimum ip significance for the track to be counted
    495   set SigMin 6.5
    496 
    497   # minimum number of tracks (high efficiency n=2, high purity n=3)
    498   set Ntracks 3
    499 }
    500 
     524  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     525  # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
     526  # gluon's PDG code has the lowest priority
     527
     528  # default efficiency formula (misidentification rate)
     529  add EfficiencyFormula {0} {0.001}
     530
     531  # efficiency formula for c-jets (misidentification rate)
     532  add EfficiencyFormula {4} {0.05}
     533
     534  # efficiency formula for b-jets
     535  add EfficiencyFormula {5} {0.75}
     536}
    501537
    502538#############
    503539# tau-tagging
    504540#############
    505 
    506541
    507542module TauTagging TauTagging {
    … …  
    514549  set TauPTMin 1.0
    515550
    516   set TauEtaMax 4.0
     551  set TauEtaMax 6.0
    517552
    518553  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    … …  
    521556  add EfficiencyFormula {0} {0.001}
    522557  # efficiency formula for tau-jets
    523   add EfficiencyFormula {15} {0.4}
     558  add EfficiencyFormula {15} {0.6}
     559}
     560
     561#####################################################
     562# Find uniquely identified photons/electrons/tau/jets
     563#####################################################
     564
     565module UniqueObjectFinder UniqueObjectFinder {
     566# earlier arrays take precedence over later ones
     567# add InputArray InputArray OutputArray
     568  add InputArray PhotonIsolation/photons photons
     569  add InputArray ElectronIsolation/electrons electrons
     570  add InputArray MuonIsolation/muons muons
     571  add InputArray JetEnergyScale/jets jets
    524572}
    525573
    … …  
    527575# ROOT tree writer
    528576##################
     577
     578# tracks, towers and eflow objects are not stored by default in the output.
     579# if needed (for jet constituent or other studies), uncomment the relevant
     580# "add Branch ..." lines.
    529581
    530582module TreeWriter TreeWriter {
    531583# add Branch InputArray BranchName BranchClass
    532584  add Branch Delphes/allParticles Particle GenParticle
     585
     586  add Branch TrackMerger/tracks Track Track
     587  add Branch Calorimeter/towers Tower Tower
     588
     589  add Branch Calorimeter/eflowTracks EFlowTrack Track
     590  add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
     591  add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
     592
    533593  add Branch GenJetFinder/jets GenJet Jet
    534 
    535   add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
    536   add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
    537   add Branch ECal/eflowPhotons Photon Photon
    538 
    539   add Branch ElectronFilter/electrons Electron Electron
    540   add Branch MuonMomentumSmearing/muons Muon Muon
    541   add Branch JetEnergyScale/jets Jet Jet
     594  add Branch UniqueObjectFinder/jets Jet Jet
     595  add Branch UniqueObjectFinder/electrons Electron Electron
     596  add Branch UniqueObjectFinder/photons Photon Photon
     597  add Branch UniqueObjectFinder/muons Muon Muon
    542598  add Branch MissingET/momentum MissingET MissingET
    543599  add Branch ScalarHT/energy ScalarHT ScalarHT
Note: See TracChangeset for help on using the changeset viewer.