Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • cards/delphes_card_IDEA.tcl

    r3051358a rb8a6aa3  
    11####################################################################                                l
    2 # FCC-ee IDEA detector model                                                                                     
    3 #                                                                                                   
    4 # Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello                                 
     2# FCC-ee IDEA detector model
     3#
     4# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
    55# email: efontane@bo.infn.it,
    6 #        lorenzo.pezzotti01@universitadipavia.it,                                                   
    7 #        m.antonello@uninsubria.it,                                                                 
    8 #####################################################################                               
    9 #       
     6#        lorenzo.pezzotti01@universitadipavia.it,
     7#        m.antonello@uninsubria.it,
     8#        michele.selvaggi@cern.ch
     9#####################################################################
     10
     11set B 2.0
     12
    1013#######################################
    1114# Order of execution of various modules
     
    1922  MuonTrackingEfficiency
    2023
    21   ChargedHadronMomentumSmearing
    22   ElectronMomentumSmearing
    23   MuonMomentumSmearing
    24 
    25   TrackMerger
     24  TrackMergerPre
     25  TrackSmearing
     26
     27  TrackMerger
    2628  Calorimeter
    2729  EFlowMerger
     
    2931  PhotonEfficiency
    3032  PhotonIsolation
     33
     34  MuonFilter
    3135
    3236  ElectronFilter
     
    4246  GenJetFinder
    4347  GenMissingET
    44  
     48
    4549  FastJetFinder
    4650
     
    5761  TreeWriter
    5862}
    59 
    6063
    6164#################################
     
    7881
    7982  # magnetic field, in T
    80   set Bz 2.0
     83  set Bz $B
    8184}
    8285
     
    98101}
    99102
    100 #       (pt <= 0.1)                                     * (0.00) + 
     103#       (pt <= 0.1)                                     * (0.00) +
    101104#       (abs(eta) <= 3.0)               * (pt > 0.1)    * (1.00) +
    102105#       (abs(eta) > 3)                                  * (0.00)
     
    119122        (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0)    * (0.65) +
    120123        (energy < 0.3) * (abs(eta) <= 3.0)                     * (0.06)
    121     } 
     124    }
    122125}
    123126
     
    140143}
    141144
    142 
    143 ########################################
    144 # Momentum resolution for charged tracks
    145 ########################################
    146 
    147 module MomentumSmearing ChargedHadronMomentumSmearing {
    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.0001145^2 + 0.0002024^2*pt + (pt*2.093e-005)^2)
    156     }
    157 }
    158 
    159 ###################################
    160 # Momentum resolution for electrons
    161 ###################################
    162 
    163 module MomentumSmearing ElectronMomentumSmearing {
    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.0001145^2 + 0.0002024^2*pt + (pt*2.093e-005)^2)
    171    } 
    172 }
    173 
    174 ###############################
    175 # Momentum resolution for muons
    176 ###############################
    177 
    178 module MomentumSmearing MuonMomentumSmearing {
    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.0001145^2 + 0.0002024^2*pt + (pt*2.093e-005)^2)
    186     }
    187 }
    188 
    189145##############
    190146# Track merger
    191147##############
    192148
     149module Merger TrackMergerPre {
     150# add InputArray InputArray
     151  add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
     152  add InputArray ElectronTrackingEfficiency/electrons
     153  add InputArray MuonTrackingEfficiency/muons
     154  set OutputArray tracks
     155}
     156
     157
     158########################################
     159# Smearing for charged tracks
     160########################################
     161
     162module TrackCovariance TrackSmearing {
     163    set InputArray TrackMergerPre/tracks
     164    set OutputArray tracks
     165
     166
     167    set InputArray TrackMergerPre/tracks
     168    set OutputArray tracks
     169
     170    set Bz 2.0
     171
     172    ## minimum number of hits to accept a track
     173    set NMinHits 6
     174
     175    ## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt
     176    set DetectorGeometry {
     177
     178
     179      # Layer type 1 = R (barrel) or 2 = z (forward/backward)
     180      # Layer label
     181      # Minimum dimension z for barrel or R for forward
     182      # Maximum dimension z for barrel or R for forward
     183      # R/z location of layer
     184      # Thickness (meters)
     185      # Radiation length (meters)
     186      # Number of measurements in layers (1D or 2D)
     187      # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Upper side
     188      # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Lower side
     189      # Resolution Upper side (meters) - 0 = no measurement
     190      # Resolution Lower side (meters) - 0 = no measurement
     191      # measurement flag = T, scattering only = F
     192
     193      # barrel  name       zmin   zmax   r        w (m)      X0        n_meas  th_up (rad) th_down (rad)    reso_up (m)   reso_down (m)  flag
     194
     195      # barrel  name       zmin   zmax   r        w (m)      X0        n_meas  th_up (rad) th_down (rad)    reso_up (m)   reso_down (m)  flag
     196
     197      1        PIPE       -100    100    0.015    0.001655  0.2805     0        0          0                0             0              0
     198      1        VTXLOW     -0.12   0.12   0.017    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     199      1        VTXLOW     -0.16   0.16   0.023    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     200      1        VTXLOW     -0.16   0.16   0.031    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     201      1        VTXHIGH    -1      1      0.32     0.00047   0.0937     2        0          1.5708           7e-006        7e-006         1
     202      1        VTXHIGH    -1.05   1.05   0.34     0.00047   0.0937     2        0          1.5708           7e-006        7e-006         1
     203
     204      # endcap  name       rmin   rmax   z        w (m)      X0        n_meas   th_up (rad)  th_down (rad)   reso_up (m)   reso_down (m) flag
     205
     206      2        VTXDSK      0.141  0.3   -0.92     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     207      2        VTXDSK      0.138  0.3   -0.9      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     208      2        VTXDSK      0.065  0.3   -0.42     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     209      2        VTXDSK      0.062  0.3   -0.4      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     210      2        VTXDSK      0.062  0.3    0.4      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     211      2        VTXDSK      0.065  0.3    0.42     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     212      2        VTXDSK      0.138  0.3    0.9      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     213      2        VTXDSK      0.141  0.3    0.92     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     214
     215      1 DCHCANI -2.125 2.125 0.345 0.0002 0.237223 0 0 0 0 0 0
     216      1 DCH -2 2 0.36 0.0147748 1400 1 0.0203738 0 0.0001 0 1
     217      1 DCH -2 2 0.374775 0.0147748 1400 1 -0.0212097 0 0.0001 0 1
     218      1 DCH -2 2 0.38955 0.0147748 1400 1 0.0220456 0 0.0001 0 1
     219      1 DCH -2 2 0.404324 0.0147748 1400 1 -0.0228814 0 0.0001 0 1
     220      1 DCH -2 2 0.419099 0.0147748 1400 1 0.0237172 0 0.0001 0 1
     221      1 DCH -2 2 0.433874 0.0147748 1400 1 -0.024553 0 0.0001 0 1
     222      1 DCH -2 2 0.448649 0.0147748 1400 1 0.0253888 0 0.0001 0 1
     223      1 DCH -2 2 0.463423 0.0147748 1400 1 -0.0262245 0 0.0001 0 1
     224      1 DCH -2 2 0.478198 0.0147748 1400 1 0.0270602 0 0.0001 0 1
     225      1 DCH -2 2 0.492973 0.0147748 1400 1 -0.0278958 0 0.0001 0 1
     226      1 DCH -2 2 0.507748 0.0147748 1400 1 0.0287314 0 0.0001 0 1
     227      1 DCH -2 2 0.522523 0.0147748 1400 1 -0.029567 0 0.0001 0 1
     228      1 DCH -2 2 0.537297 0.0147748 1400 1 0.0304025 0 0.0001 0 1
     229      1 DCH -2 2 0.552072 0.0147748 1400 1 -0.031238 0 0.0001 0 1
     230      1 DCH -2 2 0.566847 0.0147748 1400 1 0.0320734 0 0.0001 0 1
     231      1 DCH -2 2 0.581622 0.0147748 1400 1 -0.0329088 0 0.0001 0 1
     232      1 DCH -2 2 0.596396 0.0147748 1400 1 0.0337442 0 0.0001 0 1
     233      1 DCH -2 2 0.611171 0.0147748 1400 1 -0.0345795 0 0.0001 0 1
     234      1 DCH -2 2 0.625946 0.0147748 1400 1 0.0354147 0 0.0001 0 1
     235      1 DCH -2 2 0.640721 0.0147748 1400 1 -0.0362499 0 0.0001 0 1
     236      1 DCH -2 2 0.655495 0.0147748 1400 1 0.0370851 0 0.0001 0 1
     237      1 DCH -2 2 0.67027 0.0147748 1400 1 -0.0379202 0 0.0001 0 1
     238      1 DCH -2 2 0.685045 0.0147748 1400 1 0.0387552 0 0.0001 0 1
     239      1 DCH -2 2 0.69982 0.0147748 1400 1 -0.0395902 0 0.0001 0 1
     240      1 DCH -2 2 0.714595 0.0147748 1400 1 0.0404252 0 0.0001 0 1
     241      1 DCH -2 2 0.729369 0.0147748 1400 1 -0.04126 0 0.0001 0 1
     242      1 DCH -2 2 0.744144 0.0147748 1400 1 0.0420949 0 0.0001 0 1
     243      1 DCH -2 2 0.758919 0.0147748 1400 1 -0.0429296 0 0.0001 0 1
     244      1 DCH -2 2 0.773694 0.0147748 1400 1 0.0437643 0 0.0001 0 1
     245      1 DCH -2 2 0.788468 0.0147748 1400 1 -0.044599 0 0.0001 0 1
     246      1 DCH -2 2 0.803243 0.0147748 1400 1 0.0454336 0 0.0001 0 1
     247      1 DCH -2 2 0.818018 0.0147748 1400 1 -0.0462681 0 0.0001 0 1
     248      1 DCH -2 2 0.832793 0.0147748 1400 1 0.0471025 0 0.0001 0 1
     249      1 DCH -2 2 0.847568 0.0147748 1400 1 -0.0479369 0 0.0001 0 1
     250      1 DCH -2 2 0.862342 0.0147748 1400 1 0.0487713 0 0.0001 0 1
     251      1 DCH -2 2 0.877117 0.0147748 1400 1 -0.0496055 0 0.0001 0 1
     252      1 DCH -2 2 0.891892 0.0147748 1400 1 0.0504397 0 0.0001 0 1
     253      1 DCH -2 2 0.906667 0.0147748 1400 1 -0.0512738 0 0.0001 0 1
     254      1 DCH -2 2 0.921441 0.0147748 1400 1 0.0521079 0 0.0001 0 1
     255      1 DCH -2 2 0.936216 0.0147748 1400 1 -0.0529418 0 0.0001 0 1
     256      1 DCH -2 2 0.950991 0.0147748 1400 1 0.0537757 0 0.0001 0 1
     257      1 DCH -2 2 0.965766 0.0147748 1400 1 -0.0546095 0 0.0001 0 1
     258      1 DCH -2 2 0.980541 0.0147748 1400 1 0.0554433 0 0.0001 0 1
     259      1 DCH -2 2 0.995315 0.0147748 1400 1 -0.056277 0 0.0001 0 1
     260      1 DCH -2 2 1.01009 0.0147748 1400 1 0.0571106 0 0.0001 0 1
     261      1 DCH -2 2 1.02486 0.0147748 1400 1 -0.0579441 0 0.0001 0 1
     262      1 DCH -2 2 1.03964 0.0147748 1400 1 0.0587775 0 0.0001 0 1
     263      1 DCH -2 2 1.05441 0.0147748 1400 1 -0.0596108 0 0.0001 0 1
     264      1 DCH -2 2 1.06919 0.0147748 1400 1 0.0604441 0 0.0001 0 1
     265      1 DCH -2 2 1.08396 0.0147748 1400 1 -0.0612773 0 0.0001 0 1
     266      1 DCH -2 2 1.09874 0.0147748 1400 1 0.0621104 0 0.0001 0 1
     267      1 DCH -2 2 1.11351 0.0147748 1400 1 -0.0629434 0 0.0001 0 1
     268      1 DCH -2 2 1.12829 0.0147748 1400 1 0.0637763 0 0.0001 0 1
     269      1 DCH -2 2 1.14306 0.0147748 1400 1 -0.0646092 0 0.0001 0 1
     270      1 DCH -2 2 1.15784 0.0147748 1400 1 0.0654419 0 0.0001 0 1
     271      1 DCH -2 2 1.17261 0.0147748 1400 1 -0.0662746 0 0.0001 0 1
     272      1 DCH -2 2 1.18739 0.0147748 1400 1 0.0671071 0 0.0001 0 1
     273      1 DCH -2 2 1.20216 0.0147748 1400 1 -0.0679396 0 0.0001 0 1
     274      1 DCH -2 2 1.21694 0.0147748 1400 1 0.068772 0 0.0001 0 1
     275      1 DCH -2 2 1.23171 0.0147748 1400 1 -0.0696042 0 0.0001 0 1
     276      1 DCH -2 2 1.24649 0.0147748 1400 1 0.0704364 0 0.0001 0 1
     277      1 DCH -2 2 1.26126 0.0147748 1400 1 -0.0712685 0 0.0001 0 1
     278      1 DCH -2 2 1.27604 0.0147748 1400 1 0.0721005 0 0.0001 0 1
     279      1 DCH -2 2 1.29081 0.0147748 1400 1 -0.0729324 0 0.0001 0 1
     280      1 DCH -2 2 1.30559 0.0147748 1400 1 0.0737642 0 0.0001 0 1
     281      1 DCH -2 2 1.32036 0.0147748 1400 1 -0.0745958 0 0.0001 0 1
     282      1 DCH -2 2 1.33514 0.0147748 1400 1 0.0754274 0 0.0001 0 1
     283      1 DCH -2 2 1.34991 0.0147748 1400 1 -0.0762589 0 0.0001 0 1
     284      1 DCH -2 2 1.36468 0.0147748 1400 1 0.0770903 0 0.0001 0 1
     285      1 DCH -2 2 1.37946 0.0147748 1400 1 -0.0779215 0 0.0001 0 1
     286      1 DCH -2 2 1.39423 0.0147748 1400 1 0.0787527 0 0.0001 0 1
     287      1 DCH -2 2 1.40901 0.0147748 1400 1 -0.0795837 0 0.0001 0 1
     288      1 DCH -2 2 1.42378 0.0147748 1400 1 0.0804147 0 0.0001 0 1
     289      1 DCH -2 2 1.43856 0.0147748 1400 1 -0.0812455 0 0.0001 0 1
     290      1 DCH -2 2 1.45333 0.0147748 1400 1 0.0820762 0 0.0001 0 1
     291      1 DCH -2 2 1.46811 0.0147748 1400 1 -0.0829068 0 0.0001 0 1
     292      1 DCH -2 2 1.48288 0.0147748 1400 1 0.0837373 0 0.0001 0 1
     293      1 DCH -2 2 1.49766 0.0147748 1400 1 -0.0845677 0 0.0001 0 1
     294      1 DCH -2 2 1.51243 0.0147748 1400 1 0.0853979 0 0.0001 0 1
     295      1 DCH -2 2 1.52721 0.0147748 1400 1 -0.086228 0 0.0001 0 1
     296      1 DCH -2 2 1.54198 0.0147748 1400 1 0.087058 0 0.0001 0 1
     297      1 DCH -2 2 1.55676 0.0147748 1400 1 -0.0878879 0 0.0001 0 1
     298      1 DCH -2 2 1.57153 0.0147748 1400 1 0.0887177 0 0.0001 0 1
     299      1 DCH -2 2 1.58631 0.0147748 1400 1 -0.0895474 0 0.0001 0 1
     300      1 DCH -2 2 1.60108 0.0147748 1400 1 0.0903769 0 0.0001 0 1
     301      1 DCH -2 2 1.61586 0.0147748 1400 1 -0.0912063 0 0.0001 0 1
     302      1 DCH -2 2 1.63063 0.0147748 1400 1 0.0920356 0 0.0001 0 1
     303      1 DCH -2 2 1.64541 0.0147748 1400 1 -0.0928647 0 0.0001 0 1
     304      1 DCH -2 2 1.66018 0.0147748 1400 1 0.0936937 0 0.0001 0 1
     305      1 DCH -2 2 1.67495 0.0147748 1400 1 -0.0945226 0 0.0001 0 1
     306      1 DCH -2 2 1.68973 0.0147748 1400 1 0.0953514 0 0.0001 0 1
     307      1 DCH -2 2 1.7045 0.0147748 1400 1 -0.09618 0 0.0001 0 1
     308      1 DCH -2 2 1.71928 0.0147748 1400 1 0.0970085 0 0.0001 0 1
     309      1 DCH -2 2 1.73405 0.0147748 1400 1 -0.0978369 0 0.0001 0 1
     310      1 DCH -2 2 1.74883 0.0147748 1400 1 0.0986651 0 0.0001 0 1
     311      1 DCH -2 2 1.7636 0.0147748 1400 1 -0.0994932 0 0.0001 0 1
     312      1 DCH -2 2 1.77838 0.0147748 1400 1 0.100321 0 0.0001 0 1
     313      1 DCH -2 2 1.79315 0.0147748 1400 1 -0.101149 0 0.0001 0 1
     314      1 DCH -2 2 1.80793 0.0147748 1400 1 0.101977 0 0.0001 0 1
     315      1 DCH -2 2 1.8227 0.0147748 1400 1 -0.102804 0 0.0001 0 1
     316      1 DCH -2 2 1.83748 0.0147748 1400 1 0.103632 0 0.0001 0 1
     317      1 DCH -2 2 1.85225 0.0147748 1400 1 -0.104459 0 0.0001 0 1
     318      1 DCH -2 2 1.86703 0.0147748 1400 1 0.105286 0 0.0001 0 1
     319      1 DCH -2 2 1.8818 0.0147748 1400 1 -0.106113 0 0.0001 0 1
     320      1 DCH -2 2 1.89658 0.0147748 1400 1 0.10694 0 0.0001 0 1
     321      1 DCH -2 2 1.91135 0.0147748 1400 1 -0.107766 0 0.0001 0 1
     322      1 DCH -2 2 1.92613 0.0147748 1400 1 0.108593 0 0.0001 0 1
     323      1 DCH -2 2 1.9409 0.0147748 1400 1 -0.109419 0 0.0001 0 1
     324      1 DCH -2 2 1.95568 0.0147748 1400 1 0.110246 0 0.0001 0 1
     325      1 DCH -2 2 1.97045 0.0147748 1400 1 -0.111072 0 0.0001 0 1
     326      1 DCH -2 2 1.98523 0.0147748 1400 1 0.111898 0 0.0001 0 1
     327      1 DCH -2 2 2 0.0147748 1400 1 -0.112723 0 0.0001 0 1
     328      1 DCHCANO -2.125 2.125 2.02 0.02 1.667 0 0 0 0 0 0
     329      1 BSILWRP -2.35 2.35 2.04 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     330      1 BSILWRP -2.35 2.35 2.06 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     331      1 MAG -2.5 2.5 2.25 0.05 0.0658 0 0 0 0 0 0
     332      1 BPRESH -2.55 2.55 2.45 0.02 1 2 0 1.5708 7e-005 0.01 1
     333      2 DCHWALL 0.345 2.02 2.125 0.25 5.55 0 0 0 0 0 0
     334      2 DCHWALL 0.345 2.02 -2.125 0.25 5.55 0 0 0 0 0 0
     335      2 FSILWRP 0.354 2.02 -2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     336      2 FSILWRP 0.35 2.02 -2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     337      2 FSILWRP 0.35 2.02 2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     338      2 FSILWRP 0.354 2.02 2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     339      2 FRAD 0.38 2.09 2.49 0.0043 0.005612 0 0 0 0 0 0
     340      2 FRAD 0.38 2.09 -2.49 0.0043 0.005612 0 0 0 0 0 0
     341      2 FPRESH 0.39 2.43 -2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
     342      2 FPRESH 0.39 2.43 2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
     343    }
     344
     345    set Bz $B
     346}
     347
     348##############
     349# Track merger
     350##############
     351
    193352module Merger TrackMerger {
    194353# add InputArray InputArray
    195   add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    196   add InputArray ElectronMomentumSmearing/electrons
    197   add InputArray MuonMomentumSmearing/muons
     354  add InputArray TrackSmearing/tracks
    198355  set OutputArray tracks
    199356}
    200357
    201358
    202 #############                                                                                                                         
    203 # Calorimeter                                                                                                                                           
    204 #############                                                                                                                                           
     359#############
     360# Calorimeter
     361#############
    205362module DualReadoutCalorimeter Calorimeter {
    206363  set ParticleInputArray ParticlePropagator/stableParticles
     
    224381    set pi [expr {acos(-1)}]
    225382
    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:                                                                                       
     383    # Lists of the edges of each tower in eta and phi;
     384    # each list starts with the lower edge of the first tower;
     385    # the list ends with the higher edged of the last tower.
     386    # Barrel:  deta=0.02 towers up to |eta| <= 0.88 ( up to 45°)
     387    # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad)
     388    # Cell size: about 6 cm x 6 cm
     389
     390    #barrel:
    234391    set PhiBins {}
    235392    for {set i -120} {$i <= 120} {incr i} {
    236393        add PhiBins [expr {$i * $pi/120}]
    237394    }
    238     #deta=0.02 units for |eta| <= 0.88                                                             
     395    #deta=0.02 units for |eta| <= 0.88
    239396    for {set i -44} {$i < 45} {incr i} {
    240397        set eta [expr {$i * 0.02}]
     
    242399    }
    243400
    244     #endcaps:                                                                                       
     401    #endcaps:
    245402    set PhiBins {}
    246403    for {set i -120} {$i <= 120} {incr i} {
    247404        add PhiBins [expr {$i* $pi/120}]
    248405    }
    249     #deta=0.02 units for 0.88 < |eta| <= 3.0                                                       
    250     #first, from -3.0 to -0.88                                                                     
     406    #deta=0.02 units for 0.88 < |eta| <= 3.0
     407    #first, from -3.0 to -0.88
    251408    for {set i 1} {$i <=106} {incr i} {
    252409        set eta [expr {-3.00 + $i * 0.02}]
    253410        add EtaPhiBins $eta $PhiBins
    254411    }
    255     #same for 0.88 to 3.0                                                                           
     412    #same for 0.88 to 3.0
    256413    for  {set i 1} {$i <=106} {incr i} {
    257414        set eta [expr {0.88 + $i * 0.02}]
     
    259416    }
    260417
    261     # default energy fractions {abs(PDG code)} {Fecal Fhcal}                                                                                 
     418    # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    262419    add EnergyFraction {0} {0.0 1.0}
    263     # energy fractions for e, gamma and pi0                                                                                                           
     420    # energy fractions for e, gamma and pi0
    264421    add EnergyFraction {11} {1.0 0.0}
    265422    add EnergyFraction {22} {1.0 0.0}
    266423    add EnergyFraction {111} {1.0 0.0}
    267     # energy fractions for muon, neutrinos and neutralinos                                                                             
     424    # energy fractions for muon, neutrinos and neutralinos
    268425    add EnergyFraction {12} {0.0 0.0}
    269426    add EnergyFraction {13} {0.0 0.0}
     
    275432    add EnergyFraction {1000035} {0.0 0.0}
    276433    add EnergyFraction {1000045} {0.0 0.0}
    277     # energy fractions for K0short and Lambda                                                                                                           
     434    # energy fractions for K0short and Lambda
    278435    add EnergyFraction {310} {0.3 0.7}
    279436    add EnergyFraction {3122} {0.3 0.7}
    280437
    281438
    282     # set ECalResolutionFormula {resolution formula as a function of eta and energy}                               
     439    # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    283440    set ECalResolutionFormula {
    284441    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.11^2)+
     
    286443    }
    287444
    288     # set HCalResolutionFormula {resolution formula as a function of eta and energy}                                               
     445    # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    289446    set HCalResolutionFormula {
    290447    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.30^2)+
     
    352509}
    353510
     511#################
     512# Muon filter
     513#################
     514
     515module PdgCodeFilter MuonFilter {
     516  set InputArray Calorimeter/eflowTracks
     517  set OutputArray muons
     518  set Invert true
     519  add PdgCode {13}
     520  add PdgCode {-13}
     521}
     522
     523
    354524#####################
    355525# Electron efficiency
     
    363533
    364534  # efficiency formula for electrons
    365   set EfficiencyFormula {         
     535  set EfficiencyFormula {
    366536        (energy < 2.0)                                         * (0.000)+
    367537        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
     
    393563
    394564module Efficiency MuonEfficiency {
    395   set InputArray MuonMomentumSmearing/muons
     565  set InputArray MuonFilter/muons
    396566  set OutputArray muons
    397567
     
    399569
    400570  # efficiency formula for muons
    401   set EfficiencyFormula {                                   
     571  set EfficiencyFormula {
    402572        (energy < 2.0)                                         * (0.000)+
    403573        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
     
    547717
    548718  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    549  
     719
    550720  # default efficiency formula (misidentification rate)
    551721  add EfficiencyFormula {0} {0.01}
     
    603773module TreeWriter TreeWriter {
    604774    # add Branch InputArray BranchName BranchClass
    605    
     775
    606776    add Branch Delphes/allParticles Particle GenParticle
    607777
    608778    add Branch TrackMerger/tracks Track Track
    609779    add Branch Calorimeter/towers Tower Tower
    610    
     780
    611781    add Branch Calorimeter/eflowTracks EFlowTrack Track
    612782    add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
     
    616786    add Branch PhotonEfficiency/photons PhotonEff Photon
    617787    add Branch PhotonIsolation/photons PhotonIso Photon
    618    
     788
    619789    add Branch GenJetFinder/jets GenJet Jet
    620790    add Branch GenMissingET/momentum GenMissingET MissingET
    621    
     791
    622792    add Branch UniqueObjectFinder/jets Jet Jet
    623793    add Branch UniqueObjectFinder/electrons Electron Electron
    624794    add Branch UniqueObjectFinder/photons Photon Photon
    625795    add Branch UniqueObjectFinder/muons Muon Muon
    626    
    627     add Branch JetEnergyScale/jets AntiKtJet Jet 
    628    
     796
     797    add Branch JetEnergyScale/jets AntiKtJet Jet
     798
    629799    add Branch MissingET/momentum MissingET MissingET
    630800    add Branch ScalarHT/energy ScalarHT ScalarHT
    631 }
    632 
     801
     802    # add Info InfoName InfoValue
     803    add Info Bz $B
     804}
Note: See TracChangeset for help on using the changeset viewer.