Fork me on GitHub

Changeset 7c9f3da in git


Ignore:
Timestamp:
Aug 19, 2019, 10:18:18 AM (5 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, master
Children:
a5ff49f
Parents:
0e7d64a
Message:

update IDEA card

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_IDEA.tcl

    r0e7d64a r7c9f3da  
    1 set MaxEvents 10
    2 
     1####################################################################                                l
     2# FCC-ee IDEA detector model                                                                                     
     3#                                                                                                   
     4# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello                                 
     5# email: efontane@bo.infn.it,
     6#        lorenzo.pezzotti01@universitadipavia.it,                                                   
     7#        m.antonello@uninsubria.it,                                                                 
     8#####################################################################                               
     9#       
    310#######################################
    411# Order of execution of various modules
    … …  
    1623  MuonMomentumSmearing
    1724
    18   TrackMerger
     25  TrackMerger 
    1926  Calorimeter
    2027  EFlowMerger
    … …  
    4855
    4956  ScalarHT
    50 
    5157  TreeWriter
    5258}
     59
    5360
    5461#################################
    … …  
    6471  set MuonOutputArray muons
    6572
    66   # radius of the magnetic field coverage, in m
    67   set Radius 1.29
    68   # half-length of the magnetic field coverage, in m
    69   set HalfLength 3.00
    70 
    71   # magnetic field
    72   set Bz 3.8
     73  # inner radius of the solenoid, in m
     74  set Radius 2.25
     75
     76  # half-length: z of the solenoid, in m
     77  set HalfLength 2.5
     78
     79  # magnetic field, in T
     80  set Bz 2.0
    7381}
    7482
    … …  
    7886
    7987module Efficiency ChargedHadronTrackingEfficiency {
    80   set InputArray ParticlePropagator/chargedHadrons
    81   set OutputArray chargedHadrons
    82 
    83   # add EfficiencyFormula {efficiency formula as a function of eta and pt}
    84 
    85   # tracking efficiency formula for charged hadrons
    86   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    87                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
    88                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
    89                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
    90                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    91                          (abs(eta) > 2.5)                                                  * (0.00)}
    92 }
     88    set InputArray ParticlePropagator/chargedHadrons
     89    set OutputArray chargedHadrons
     90    # We use only one efficiency, we set only 0 effincency out of eta bounds:
     91
     92    set EfficiencyFormula {
     93        (abs(eta) > 3.0)                                       * (0.000) +
     94        (energy >= 0.5) * (abs(eta) <= 3.0)                    * (0.997) +
     95        (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0)    * (0.65) +
     96        (energy < 0.3) * (abs(eta) <= 3.0)                     * (0.06)
     97    }
     98}
     99
     100#       (pt <= 0.1)                                     * (0.00) +
     101#       (abs(eta) <= 3.0)               * (pt > 0.1)    * (1.00) +
     102#       (abs(eta) > 3)                                  * (0.00)
     103
     104
    93105
    94106##############################
    … …  
    97109
    98110module Efficiency ElectronTrackingEfficiency {
    99   set InputArray ParticlePropagator/electrons
    100   set OutputArray electrons
    101 
    102   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    103 
    104   # tracking efficiency formula for electrons
    105   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    106                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
    107                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
    108                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
    109                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
    110                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
    111                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    112                          (abs(eta) > 2.5)                                                  * (0.00)}
    113 }
     111    set InputArray ParticlePropagator/electrons
     112    set OutputArray electrons
     113
     114
     115    # Current full simulation with CLICdet provides for electrons:
     116    set EfficiencyFormula {
     117        (abs(eta) > 3.0)                                       * (0.000) +
     118        (energy >= 0.5) * (abs(eta) <= 3.0)                    * (0.997) +
     119        (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0)    * (0.65) +
     120        (energy < 0.3) * (abs(eta) <= 3.0)                     * (0.06)
     121    }
     122}
     123
    114124
    115125##########################
    … …  
    118128
    119129module Efficiency MuonTrackingEfficiency {
    120   set InputArray ParticlePropagator/muons
    121   set OutputArray muons
    122 
    123   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    124 
    125   # tracking efficiency formula for muons
    126   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    127                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
    128                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
    129                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
    130                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    131                          (abs(eta) > 2.5)                                                  * (0.00)}
    132 }
     130    set InputArray ParticlePropagator/muons
     131    set OutputArray muons
     132
     133    # Current full simulation with CLICdet provides for muons:
     134    set EfficiencyFormula {
     135        (abs(eta) > 3.0)                                       * (0.000) +
     136        (energy >= 0.5) * (abs(eta) <= 3.0)                    * (0.997) +
     137        (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0)    * (0.65) +
     138        (energy < 0.3) * (abs(eta) <= 3.0)                     * (0.06)
     139    }
     140}
     141
    133142
    134143########################################
    … …  
    137146
    138147module MomentumSmearing ChargedHadronMomentumSmearing {
    139   set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
    140   set OutputArray chargedHadrons
    141 
    142   # set ResolutionFormula {resolution formula as a function of eta and pt}
    143 
    144   # resolution formula for charged hadrons
    145   # based on arXiv:1405.6569
    146   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
    147                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
    148                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
     148    set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
     149    set OutputArray chargedHadrons
     150
     151
     152    # Resolution given in dpT/pT.
     153    # IDEAdet
     154    set ResolutionFormula {
     155        (abs(eta) <= 3.0)                   * sqrt(0.0006^2 + (pt*7.e-5)^2)
     156    }
    149157}
    150158
    … …  
    154162
    155163module MomentumSmearing ElectronMomentumSmearing {
    156   set InputArray ElectronTrackingEfficiency/electrons
    157   set OutputArray electrons
    158 
    159   # set ResolutionFormula {resolution formula as a function of eta and energy}
    160 
    161   # resolution formula for electrons
    162   # based on arXiv:1405.6569
    163   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
    164                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
    165                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
     164    set InputArray ElectronTrackingEfficiency/electrons
     165    set OutputArray electrons
     166
     167    # Resolution given in dpT/pT.
     168    # IDEAdet
     169    set ResolutionFormula {
     170        (abs(eta) <= 3.0)                   * sqrt(0.0006^2 + (pt*7.e-5)^2)
     171    } 
    166172}
    167173
    … …  
    171177
    172178module MomentumSmearing MuonMomentumSmearing {
    173   set InputArray MuonTrackingEfficiency/muons
    174   set OutputArray muons
    175 
    176   # set ResolutionFormula {resolution formula as a function of eta and pt}
    177 
    178   # resolution formula for muons
    179   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
    180                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
    181                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
     179    set InputArray MuonTrackingEfficiency/muons
     180    set OutputArray muons
     181
     182    # Resolution given in dpT/pT.
     183    # IDEAdet
     184    set ResolutionFormula {
     185        (abs(eta) <= 3.0)                   * sqrt(0.0006^2 + (pt*7.e-5)^2)
     186    }
    182187}
    183188
    … …  
    194199}
    195200
    196 #############
    197 # Calorimeter
    198 #############
    199 
     201
     202#############                                                                                                                         
     203# Calorimeter                                                                                                                                           
     204#############                                                                                                                                           
    200205module DualReadoutCalorimeter Calorimeter {
    201206  set ParticleInputArray ParticlePropagator/stableParticles
    … …  
    210215
    211216  set ECalEnergyMin 0.5
    212   set HCalEnergyMin 1.0
     217  set HCalEnergyMin 0.5
    213218  set EnergyMin 0.5
    214 
    215219  set ECalEnergySignificanceMin 1.0
    216220  set HCalEnergySignificanceMin 1.0
    … …  
    218222
    219223  set SmearTowerCenter true
    220 
    221   set pi [expr {acos(-1)}]
    222 
    223   # lists of the edges of each tower in eta and phi
    224   # each list starts with the lower edge of the first tower
    225   # the list ends with the higher edged of the last tower
    226 
    227   # 5 degrees towers
    228   set PhiBins {}
    229   for {set i -36} {$i <= 36} {incr i} {
    230     add PhiBins [expr {$i * $pi/36.0}]
    231   }
    232   foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
    233     add EtaPhiBins $eta $PhiBins
    234   }
    235 
    236   # 10 degrees towers
    237   set PhiBins {}
    238   for {set i -18} {$i <= 18} {incr i} {
    239     add PhiBins [expr {$i * $pi/18.0}]
    240   }
    241   foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
    242     add EtaPhiBins $eta $PhiBins
    243   }
    244 
    245   # 20 degrees towers
    246   set PhiBins {}
    247   for {set i -9} {$i <= 9} {incr i} {
    248     add PhiBins [expr {$i * $pi/9.0}]
    249   }
    250   foreach eta {-5 -4.7 -4.525 4.7 5} {
    251     add EtaPhiBins $eta $PhiBins
    252   }
    253 
    254   # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    255   add EnergyFraction {0} {0.0 1.0}
    256   # energy fractions for e, gamma and pi0
    257   add EnergyFraction {11} {1.0 0.0}
    258   add EnergyFraction {22} {1.0 0.0}
    259   add EnergyFraction {111} {1.0 0.0}
    260   # energy fractions for muon, neutrinos and neutralinos
    261   add EnergyFraction {12} {0.0 0.0}
    262   add EnergyFraction {13} {0.0 0.0}
    263   add EnergyFraction {14} {0.0 0.0}
    264   add EnergyFraction {16} {0.0 0.0}
    265   add EnergyFraction {1000022} {0.0 0.0}
    266   add EnergyFraction {1000023} {0.0 0.0}
    267   add EnergyFraction {1000025} {0.0 0.0}
    268   add EnergyFraction {1000035} {0.0 0.0}
    269   add EnergyFraction {1000045} {0.0 0.0}
    270   # energy fractions for K0short and Lambda
    271   add EnergyFraction {310} {0.3 0.7}
    272   add EnergyFraction {3122} {0.3 0.7}
    273 
    274   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    275   # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
    276   set ECalResolutionFormula {                  (abs(eta) <= 1.5) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
    277                              (abs(eta) > 1.5 && abs(eta) <= 2.5) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
    278                              (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    279 
    280   # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    281   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*0.30^2) +
    282                              (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
     224    set pi [expr {acos(-1)}]
     225
     226    # Lists of the edges of each tower in eta and phi;                                                                                           
     227    # each list starts with the lower edge of the first tower;                                                                                                   
     228    # the list ends with the higher edged of the last tower.                                                                                       
     229    # Barrel:  deta=0.02 towers up to |eta| <= 0.88 ( up to 45°)                                                                                       
     230    # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad)                                                                                           
     231    # Cell size: about 6 cm x 6 cm                                                 
     232
     233    #barrel:                                                                                       
     234    set PhiBins {}
     235    for {set i -120} {$i <= 120} {incr i} {
     236        add PhiBins [expr {$i * $pi/120}]
     237    }
     238    #deta=0.02 units for |eta| <= 0.88                                                             
     239    for {set i -44} {$i < 45} {incr i} {
     240        set eta [expr {$i * 0.02}]
     241        add EtaPhiBins $eta $PhiBins
     242    }
     243
     244    #endcaps:                                                                                       
     245    set PhiBins {}
     246    for {set i -120} {$i <= 120} {incr i} {
     247        add PhiBins [expr {$i* $pi/120}]
     248    }
     249    #deta=0.02 units for 0.88 < |eta| <= 3.0                                                       
     250    #first, from -3.0 to -0.88                                                                     
     251    for {set i 1} {$i <=106} {incr i} {
     252        set eta [expr {-3.00 + $i * 0.02}]
     253        add EtaPhiBins $eta $PhiBins
     254    }
     255    #same for 0.88 to 3.0                                                                           
     256    for  {set i 1} {$i <=106} {incr i} {
     257        set eta [expr {0.88 + $i * 0.02}]
     258        add EtaPhiBins $eta $PhiBins
     259    }
     260
     261    # default energy fractions {abs(PDG code)} {Fecal Fhcal}                                                                                 
     262    add EnergyFraction {0} {0.0 1.0}
     263    # energy fractions for e, gamma and pi0                                                                                                           
     264    add EnergyFraction {11} {1.0 0.0}
     265    add EnergyFraction {22} {1.0 0.0}
     266    add EnergyFraction {111} {1.0 0.0}
     267    # energy fractions for muon, neutrinos and neutralinos                                                                             
     268    add EnergyFraction {12} {0.0 0.0}
     269    add EnergyFraction {13} {0.0 0.0}
     270    add EnergyFraction {14} {0.0 0.0}
     271    add EnergyFraction {16} {0.0 0.0}
     272    add EnergyFraction {1000022} {0.0 0.0}
     273    add EnergyFraction {1000023} {0.0 0.0}
     274    add EnergyFraction {1000025} {0.0 0.0}
     275    add EnergyFraction {1000035} {0.0 0.0}
     276    add EnergyFraction {1000045} {0.0 0.0}
     277    # energy fractions for K0short and Lambda                                                                                                           
     278    add EnergyFraction {310} {0.3 0.7}
     279    add EnergyFraction {3122} {0.3 0.7}
     280
     281
     282    # set ECalResolutionFormula {resolution formula as a function of eta and energy}                               
     283    set ECalResolutionFormula {
     284    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.11^2)+
     285    (abs(eta) > 0.88 && abs(eta) <= 3.0)    * sqrt(energy^2*0.01^2 + energy*0.11^2)
     286    }
     287
     288    # set HCalResolutionFormula {resolution formula as a function of eta and energy}                                               
     289    set HCalResolutionFormula {
     290    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.30^2)+
     291    (abs(eta) > 0.88 && abs(eta) <= 3.0)    * sqrt(energy^2*0.01^2 + energy*0.30^2)
     292    }
    283293}
    284294
    … …  
    304314
    305315  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    306 
    307316  # efficiency formula for photons
    308   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
    309                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
    310                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    311                          (abs(eta) > 2.5)                                   * (0.00)}
     317  set EfficiencyFormula {
     318        (energy < 2.0)                                        * (0.000)+
     319        (energy >= 2.0) * (abs(eta) <= 0.88)                  * (0.99) +
     320        (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
     321        (abs(eta) > 3.0)                                      * (0.000)
     322  }
    312323}
    313324
    … …  
    326337  set PTMin 0.5
    327338
    328   set PTRatioMax 0.12
     339  set PTRatioMax 999.
    329340}
    330341
    … …  
    352363
    353364  # efficiency formula for electrons
    354   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
    355                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
    356                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    357                          (abs(eta) > 2.5)                                   * (0.00)}
     365  set EfficiencyFormula {         
     366        (energy < 2.0)                                         * (0.000)+
     367        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
     368        (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0)  * (0.99) +
     369        (abs(eta) > 3.0)                                       * (0.000)
     370  }
    358371}
    359372
    … …  
    386399
    387400  # efficiency formula for muons
    388   set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
    389                                            (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
    390                                            (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
    391                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
    392                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
    393                          (abs(eta) > 2.4)                                                 * (0.00)}
     401  set EfficiencyFormula {                                   
     402        (energy < 2.0)                                         * (0.000)+
     403        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
     404        (energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0)  * (0.99) +
     405        (abs(eta) > 3.0)                                       * (0.000)
     406  }
    394407}
    395408
    … …  
    434447}
    435448
    436 
    437449#####################
    438450# Neutrino Filter
    … …  
    452464  add PdgCode {-14}
    453465  add PdgCode {-16}
    454 
    455466}
    456467
    … …  
    467478  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    468479  set JetAlgorithm 6
    469   set ParameterR 0.5
    470 
    471   set JetPTMin 20.0
    472 }
     480  set ParameterR 0.4
     481  set JetPTMin 1.0
     482}
     483
    473484
    474485#########################
    … …  
    482493}
    483494
    484 
    485 
    486495############
    487496# Jet finder
    … …  
    496505  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    497506  set JetAlgorithm 6
    498   set ParameterR 0.5
    499 
    500   set JetPTMin 20.0
     507  set ParameterR 0.4
     508  set JetPTMin 1.0
    501509}
    502510
    … …  
    510518
    511519  # scale formula for jets
    512   set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
     520  set ScaleFormula {1.08}
    513521}
    514522
    … …  
    526534  set DeltaR 0.5
    527535  set PartonPTMin 1.0
    528   set PartonEtaMax 2.5
    529 
     536  set PartonEtaMax 3.0
    530537}
    531538
    … …  
    540547
    541548  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    542   # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
    543   # gluon's PDG code has the lowest priority
    544 
    545   # based on arXiv:1211.4462
    546549 
    547550  # default efficiency formula (misidentification rate)
    548   add EfficiencyFormula {0} {0.01+0.000038*pt}
     551  add EfficiencyFormula {0} {0.01}
    549552
    550553  # efficiency formula for c-jets (misidentification rate)
    551   add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
     554  add EfficiencyFormula {4} {0.10}
    552555
    553556  # efficiency formula for b-jets
    554   add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
     557  add EfficiencyFormula {5} {0.80}
    555558}
    556559
    … …  
    565568
    566569  set DeltaR 0.5
    567 
    568570  set TauPTMin 1.0
    569 
    570   set TauEtaMax 2.5
    571 
    572   # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     571  set TauEtaMax 3.0
    573572
    574573  # default efficiency formula (misidentification rate)
    575   add EfficiencyFormula {0} {0.01}
     574  add EfficiencyFormula {0} {0.001}
    576575  # efficiency formula for tau-jets
    577576  add EfficiencyFormula {15} {0.6}
    578577}
     578
    579579
    580580#####################################################
    … …  
    591591}
    592592
     593
     594
    593595##################
    594596# ROOT tree writer
    595597##################
    596598
    597 # tracks, towers and eflow objects are not stored by default in the output.
    598 # if needed (for jet constituent or other studies), uncomment the relevant
     599# Tracks, towers and eflow objects are not stored by default in the output.
     600# If needed (for jet constituent or other studies), uncomment the relevant
    599601# "add Branch ..." lines.
    600602
    601603module TreeWriter TreeWriter {
    602 # add Branch InputArray BranchName BranchClass
    603   add Branch Delphes/allParticles Particle GenParticle
    604 
    605   add Branch TrackMerger/tracks Track Track
    606   add Branch Calorimeter/towers Tower Tower
    607 
    608   add Branch Calorimeter/eflowTracks EFlowTrack Track
    609   add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
    610   add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
    611 
    612   add Branch GenJetFinder/jets GenJet Jet
    613   add Branch GenMissingET/momentum GenMissingET MissingET
    614 
    615   add Branch UniqueObjectFinder/jets Jet Jet
    616   add Branch UniqueObjectFinder/electrons Electron Electron
    617   add Branch UniqueObjectFinder/photons Photon Photon
    618   add Branch UniqueObjectFinder/muons Muon Muon
    619   add Branch MissingET/momentum MissingET MissingET
    620   add Branch ScalarHT/energy ScalarHT ScalarHT
    621 }
     604    # add Branch InputArray BranchName BranchClass
     605   
     606    add Branch Delphes/allParticles Particle GenParticle
     607
     608    add Branch TrackMerger/tracks Track Track
     609    add Branch Calorimeter/towers Tower Tower
     610   
     611    add Branch Calorimeter/eflowTracks EFlowTrack Track
     612    add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
     613    add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
     614
     615    add Branch Calorimeter/photons CaloPhoton Photon
     616    add Branch PhotonEfficiency/photons PhotonEff Photon
     617    add Branch PhotonIsolation/photons PhotonIso Photon
     618   
     619    add Branch GenJetFinder/jets GenJet Jet
     620    add Branch GenMissingET/momentum GenMissingET MissingET
     621   
     622    add Branch UniqueObjectFinder/jets Jet Jet
     623    add Branch UniqueObjectFinder/electrons Electron Electron
     624    add Branch UniqueObjectFinder/photons Photon Photon
     625    add Branch UniqueObjectFinder/muons Muon Muon
     626   
     627    add Branch JetEnergyScale/jets AntiKtJet Jet 
     628   
     629    add Branch MissingET/momentum MissingET MissingET
     630    add Branch ScalarHT/energy ScalarHT ScalarHT
     631}
     632
Note: See TracChangeset for help on using the changeset viewer.