Fork me on GitHub

Changeset a0db751 in git


Ignore:
Timestamp:
Dec 14, 2020, 10:54:07 AM (4 years ago)
Author:
michele <michele.selvaggi@…>
Branches:
master
Children:
f17e10d
Parents:
170a11d
Message:

updated IDEA card with latest TrackCovariance

File:
1 edited

Legend:

Unmodified
Added
Removed
  • cards/delphes_card_IDEA.tcl

    r170a11d ra0db751  
    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
    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#####################################################################
     9#
    1010#######################################
    1111# Order of execution of various modules
     
    1919  MuonTrackingEfficiency
    2020
    21   ChargedHadronMomentumSmearing
    22   ElectronMomentumSmearing
    23   MuonMomentumSmearing
    24 
    25   TrackMerger
     21  TrackMergerPre
     22  TrackSmearing
     23
     24  TrackMerger
    2625  Calorimeter
    2726  EFlowMerger
     
    3433  ElectronIsolation
    3534
     35  MuonFilter
    3636  MuonEfficiency
    3737  MuonIsolation
     
    4242  GenJetFinder
    4343  GenMissingET
    44  
     44
    4545  FastJetFinder
    4646
     
    9898}
    9999
    100 #       (pt <= 0.1)                                     * (0.00) + 
     100#       (pt <= 0.1)                                     * (0.00) +
    101101#       (abs(eta) <= 3.0)               * (pt > 0.1)    * (1.00) +
    102102#       (abs(eta) > 3)                                  * (0.00)
     
    119119        (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0)    * (0.65) +
    120120        (energy < 0.3) * (abs(eta) <= 3.0)                     * (0.06)
    121     } 
     121    }
    122122}
    123123
     
    140140}
    141141
    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 
    189142##############
    190143# Track merger
    191144##############
    192145
     146module Merger TrackMergerPre {
     147# add InputArray InputArray
     148  add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
     149  add InputArray ElectronTrackingEfficiency/electrons
     150  add InputArray MuonTrackingEfficiency/muons
     151  set OutputArray tracks
     152}
     153
     154
     155########################################
     156# Smearing for charged tracks
     157########################################
     158
     159module TrackCovariance TrackSmearing {
     160    set InputArray TrackMergerPre/tracks
     161    set OutputArray tracks
     162
     163    set Bz 2.0
     164
     165    ## minimum number of hits to accept a track
     166    set NMinHits 6
     167
     168    ## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt
     169    set DetectorGeometry {
     170
     171
     172      # Layer type 1 = R (barrel) or 2 = z (forward/backward)
     173      # Layer label
     174      # Minimum dimension z for barrel or R for forward
     175      # Maximum dimension z for barrel or R for forward
     176      # R/z location of layer
     177      # Thickness (meters)
     178      # Radiation length (meters)
     179      # Number of measurements in layers (1D or 2D)
     180      # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Upper side
     181      # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Lower side
     182      # Resolution Upper side (meters) - 0 = no measurement
     183      # Resolution Lower side (meters) - 0 = no measurement
     184      # measurement flag = T, scattering only = F
     185
     186
     187      # barrel  name       zmin   zmax   r        w (m)      X0        n_meas  th_up (rad) th_down (rad)    reso_up (m)   reso_down (m)  flag
     188
     189      1        PIPE       -100    100    0.015    0.0012    0.35276    0        0          0                0             0              0
     190      1        VTXLOW     -0.12   0.12   0.017    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     191      1        VTXLOW     -0.16   0.16   0.023    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     192      1        VTXLOW     -0.16   0.16   0.031    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     193      1        VTXHIGH    -1      1      0.32     0.00047   0.0937     2        0          1.5708           7e-006        7e-006         1
     194      1        VTXHIGH    -1.05   1.05   0.34     0.00047   0.0937     2        0          1.5708           7e-006        7e-006         1
     195
     196      # endcap  name       rmin   rmax   z        w (m)      X0        n_meas   th_up (rad)  th_down (rad)   reso_up (m)   reso_down (m) flag
     197
     198      2        VTXDSK      0.141  0.3   -0.92     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     199      2        VTXDSK      0.138  0.3   -0.9      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     200      2        VTXDSK      0.065  0.3   -0.42     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     201      2        VTXDSK      0.062  0.3   -0.4      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     202      2        VTXDSK      0.062  0.3    0.4      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     203      2        VTXDSK      0.065  0.3    0.42     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     204      2        VTXDSK      0.138  0.3    0.9      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     205      2        VTXDSK      0.141  0.3    0.92     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     206
     207      1 DCHCANI -2.125 2.125 0.345 0.0002 0.237223 0 0 0 0 0 0
     208      1 DCH -2 2 0.36 0.0147748 1400 1 0.0203738 0 0.0001 0 1
     209      1 DCH -2 2 0.374775 0.0147748 1400 1 -0.0212097 0 0.0001 0 1
     210      1 DCH -2 2 0.38955 0.0147748 1400 1 0.0220456 0 0.0001 0 1
     211      1 DCH -2 2 0.404324 0.0147748 1400 1 -0.0228814 0 0.0001 0 1
     212      1 DCH -2 2 0.419099 0.0147748 1400 1 0.0237172 0 0.0001 0 1
     213      1 DCH -2 2 0.433874 0.0147748 1400 1 -0.024553 0 0.0001 0 1
     214      1 DCH -2 2 0.448649 0.0147748 1400 1 0.0253888 0 0.0001 0 1
     215      1 DCH -2 2 0.463423 0.0147748 1400 1 -0.0262245 0 0.0001 0 1
     216      1 DCH -2 2 0.478198 0.0147748 1400 1 0.0270602 0 0.0001 0 1
     217      1 DCH -2 2 0.492973 0.0147748 1400 1 -0.0278958 0 0.0001 0 1
     218      1 DCH -2 2 0.507748 0.0147748 1400 1 0.0287314 0 0.0001 0 1
     219      1 DCH -2 2 0.522523 0.0147748 1400 1 -0.029567 0 0.0001 0 1
     220      1 DCH -2 2 0.537297 0.0147748 1400 1 0.0304025 0 0.0001 0 1
     221      1 DCH -2 2 0.552072 0.0147748 1400 1 -0.031238 0 0.0001 0 1
     222      1 DCH -2 2 0.566847 0.0147748 1400 1 0.0320734 0 0.0001 0 1
     223      1 DCH -2 2 0.581622 0.0147748 1400 1 -0.0329088 0 0.0001 0 1
     224      1 DCH -2 2 0.596396 0.0147748 1400 1 0.0337442 0 0.0001 0 1
     225      1 DCH -2 2 0.611171 0.0147748 1400 1 -0.0345795 0 0.0001 0 1
     226      1 DCH -2 2 0.625946 0.0147748 1400 1 0.0354147 0 0.0001 0 1
     227      1 DCH -2 2 0.640721 0.0147748 1400 1 -0.0362499 0 0.0001 0 1
     228      1 DCH -2 2 0.655495 0.0147748 1400 1 0.0370851 0 0.0001 0 1
     229      1 DCH -2 2 0.67027 0.0147748 1400 1 -0.0379202 0 0.0001 0 1
     230      1 DCH -2 2 0.685045 0.0147748 1400 1 0.0387552 0 0.0001 0 1
     231      1 DCH -2 2 0.69982 0.0147748 1400 1 -0.0395902 0 0.0001 0 1
     232      1 DCH -2 2 0.714595 0.0147748 1400 1 0.0404252 0 0.0001 0 1
     233      1 DCH -2 2 0.729369 0.0147748 1400 1 -0.04126 0 0.0001 0 1
     234      1 DCH -2 2 0.744144 0.0147748 1400 1 0.0420949 0 0.0001 0 1
     235      1 DCH -2 2 0.758919 0.0147748 1400 1 -0.0429296 0 0.0001 0 1
     236      1 DCH -2 2 0.773694 0.0147748 1400 1 0.0437643 0 0.0001 0 1
     237      1 DCH -2 2 0.788468 0.0147748 1400 1 -0.044599 0 0.0001 0 1
     238      1 DCH -2 2 0.803243 0.0147748 1400 1 0.0454336 0 0.0001 0 1
     239      1 DCH -2 2 0.818018 0.0147748 1400 1 -0.0462681 0 0.0001 0 1
     240      1 DCH -2 2 0.832793 0.0147748 1400 1 0.0471025 0 0.0001 0 1
     241      1 DCH -2 2 0.847568 0.0147748 1400 1 -0.0479369 0 0.0001 0 1
     242      1 DCH -2 2 0.862342 0.0147748 1400 1 0.0487713 0 0.0001 0 1
     243      1 DCH -2 2 0.877117 0.0147748 1400 1 -0.0496055 0 0.0001 0 1
     244      1 DCH -2 2 0.891892 0.0147748 1400 1 0.0504397 0 0.0001 0 1
     245      1 DCH -2 2 0.906667 0.0147748 1400 1 -0.0512738 0 0.0001 0 1
     246      1 DCH -2 2 0.921441 0.0147748 1400 1 0.0521079 0 0.0001 0 1
     247      1 DCH -2 2 0.936216 0.0147748 1400 1 -0.0529418 0 0.0001 0 1
     248      1 DCH -2 2 0.950991 0.0147748 1400 1 0.0537757 0 0.0001 0 1
     249      1 DCH -2 2 0.965766 0.0147748 1400 1 -0.0546095 0 0.0001 0 1
     250      1 DCH -2 2 0.980541 0.0147748 1400 1 0.0554433 0 0.0001 0 1
     251      1 DCH -2 2 0.995315 0.0147748 1400 1 -0.056277 0 0.0001 0 1
     252      1 DCH -2 2 1.01009 0.0147748 1400 1 0.0571106 0 0.0001 0 1
     253      1 DCH -2 2 1.02486 0.0147748 1400 1 -0.0579441 0 0.0001 0 1
     254      1 DCH -2 2 1.03964 0.0147748 1400 1 0.0587775 0 0.0001 0 1
     255      1 DCH -2 2 1.05441 0.0147748 1400 1 -0.0596108 0 0.0001 0 1
     256      1 DCH -2 2 1.06919 0.0147748 1400 1 0.0604441 0 0.0001 0 1
     257      1 DCH -2 2 1.08396 0.0147748 1400 1 -0.0612773 0 0.0001 0 1
     258      1 DCH -2 2 1.09874 0.0147748 1400 1 0.0621104 0 0.0001 0 1
     259      1 DCH -2 2 1.11351 0.0147748 1400 1 -0.0629434 0 0.0001 0 1
     260      1 DCH -2 2 1.12829 0.0147748 1400 1 0.0637763 0 0.0001 0 1
     261      1 DCH -2 2 1.14306 0.0147748 1400 1 -0.0646092 0 0.0001 0 1
     262      1 DCH -2 2 1.15784 0.0147748 1400 1 0.0654419 0 0.0001 0 1
     263      1 DCH -2 2 1.17261 0.0147748 1400 1 -0.0662746 0 0.0001 0 1
     264      1 DCH -2 2 1.18739 0.0147748 1400 1 0.0671071 0 0.0001 0 1
     265      1 DCH -2 2 1.20216 0.0147748 1400 1 -0.0679396 0 0.0001 0 1
     266      1 DCH -2 2 1.21694 0.0147748 1400 1 0.068772 0 0.0001 0 1
     267      1 DCH -2 2 1.23171 0.0147748 1400 1 -0.0696042 0 0.0001 0 1
     268      1 DCH -2 2 1.24649 0.0147748 1400 1 0.0704364 0 0.0001 0 1
     269      1 DCH -2 2 1.26126 0.0147748 1400 1 -0.0712685 0 0.0001 0 1
     270      1 DCH -2 2 1.27604 0.0147748 1400 1 0.0721005 0 0.0001 0 1
     271      1 DCH -2 2 1.29081 0.0147748 1400 1 -0.0729324 0 0.0001 0 1
     272      1 DCH -2 2 1.30559 0.0147748 1400 1 0.0737642 0 0.0001 0 1
     273      1 DCH -2 2 1.32036 0.0147748 1400 1 -0.0745958 0 0.0001 0 1
     274      1 DCH -2 2 1.33514 0.0147748 1400 1 0.0754274 0 0.0001 0 1
     275      1 DCH -2 2 1.34991 0.0147748 1400 1 -0.0762589 0 0.0001 0 1
     276      1 DCH -2 2 1.36468 0.0147748 1400 1 0.0770903 0 0.0001 0 1
     277      1 DCH -2 2 1.37946 0.0147748 1400 1 -0.0779215 0 0.0001 0 1
     278      1 DCH -2 2 1.39423 0.0147748 1400 1 0.0787527 0 0.0001 0 1
     279      1 DCH -2 2 1.40901 0.0147748 1400 1 -0.0795837 0 0.0001 0 1
     280      1 DCH -2 2 1.42378 0.0147748 1400 1 0.0804147 0 0.0001 0 1
     281      1 DCH -2 2 1.43856 0.0147748 1400 1 -0.0812455 0 0.0001 0 1
     282      1 DCH -2 2 1.45333 0.0147748 1400 1 0.0820762 0 0.0001 0 1
     283      1 DCH -2 2 1.46811 0.0147748 1400 1 -0.0829068 0 0.0001 0 1
     284      1 DCH -2 2 1.48288 0.0147748 1400 1 0.0837373 0 0.0001 0 1
     285      1 DCH -2 2 1.49766 0.0147748 1400 1 -0.0845677 0 0.0001 0 1
     286      1 DCH -2 2 1.51243 0.0147748 1400 1 0.0853979 0 0.0001 0 1
     287      1 DCH -2 2 1.52721 0.0147748 1400 1 -0.086228 0 0.0001 0 1
     288      1 DCH -2 2 1.54198 0.0147748 1400 1 0.087058 0 0.0001 0 1
     289      1 DCH -2 2 1.55676 0.0147748 1400 1 -0.0878879 0 0.0001 0 1
     290      1 DCH -2 2 1.57153 0.0147748 1400 1 0.0887177 0 0.0001 0 1
     291      1 DCH -2 2 1.58631 0.0147748 1400 1 -0.0895474 0 0.0001 0 1
     292      1 DCH -2 2 1.60108 0.0147748 1400 1 0.0903769 0 0.0001 0 1
     293      1 DCH -2 2 1.61586 0.0147748 1400 1 -0.0912063 0 0.0001 0 1
     294      1 DCH -2 2 1.63063 0.0147748 1400 1 0.0920356 0 0.0001 0 1
     295      1 DCH -2 2 1.64541 0.0147748 1400 1 -0.0928647 0 0.0001 0 1
     296      1 DCH -2 2 1.66018 0.0147748 1400 1 0.0936937 0 0.0001 0 1
     297      1 DCH -2 2 1.67495 0.0147748 1400 1 -0.0945226 0 0.0001 0 1
     298      1 DCH -2 2 1.68973 0.0147748 1400 1 0.0953514 0 0.0001 0 1
     299      1 DCH -2 2 1.7045 0.0147748 1400 1 -0.09618 0 0.0001 0 1
     300      1 DCH -2 2 1.71928 0.0147748 1400 1 0.0970085 0 0.0001 0 1
     301      1 DCH -2 2 1.73405 0.0147748 1400 1 -0.0978369 0 0.0001 0 1
     302      1 DCH -2 2 1.74883 0.0147748 1400 1 0.0986651 0 0.0001 0 1
     303      1 DCH -2 2 1.7636 0.0147748 1400 1 -0.0994932 0 0.0001 0 1
     304      1 DCH -2 2 1.77838 0.0147748 1400 1 0.100321 0 0.0001 0 1
     305      1 DCH -2 2 1.79315 0.0147748 1400 1 -0.101149 0 0.0001 0 1
     306      1 DCH -2 2 1.80793 0.0147748 1400 1 0.101977 0 0.0001 0 1
     307      1 DCH -2 2 1.8227 0.0147748 1400 1 -0.102804 0 0.0001 0 1
     308      1 DCH -2 2 1.83748 0.0147748 1400 1 0.103632 0 0.0001 0 1
     309      1 DCH -2 2 1.85225 0.0147748 1400 1 -0.104459 0 0.0001 0 1
     310      1 DCH -2 2 1.86703 0.0147748 1400 1 0.105286 0 0.0001 0 1
     311      1 DCH -2 2 1.8818 0.0147748 1400 1 -0.106113 0 0.0001 0 1
     312      1 DCH -2 2 1.89658 0.0147748 1400 1 0.10694 0 0.0001 0 1
     313      1 DCH -2 2 1.91135 0.0147748 1400 1 -0.107766 0 0.0001 0 1
     314      1 DCH -2 2 1.92613 0.0147748 1400 1 0.108593 0 0.0001 0 1
     315      1 DCH -2 2 1.9409 0.0147748 1400 1 -0.109419 0 0.0001 0 1
     316      1 DCH -2 2 1.95568 0.0147748 1400 1 0.110246 0 0.0001 0 1
     317      1 DCH -2 2 1.97045 0.0147748 1400 1 -0.111072 0 0.0001 0 1
     318      1 DCH -2 2 1.98523 0.0147748 1400 1 0.111898 0 0.0001 0 1
     319      1 DCH -2 2 2 0.0147748 1400 1 -0.112723 0 0.0001 0 1
     320      1 DCHCANO -2.125 2.125 2.02 0.02 1.667 0 0 0 0 0 0
     321      1 BSILWRP -2.35 2.35 2.04 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     322      1 BSILWRP -2.35 2.35 2.06 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     323      1 MAG -2.5 2.5 2.25 0.05 0.0658 0 0 0 0 0 0
     324      1 BPRESH -2.55 2.55 2.45 0.02 1 2 0 1.5708 7e-005 0.01 1
     325
     326
     327      2 DCHWALL 0.345 2.02 2.125 0.25 5.55 0 0 0 0 0 0
     328      2 DCHWALL 0.345 2.02 -2.125 0.25 5.55 0 0 0 0 0 0
     329      2 FSILWRP 0.354 2.02 -2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     330      2 FSILWRP 0.35 2.02 -2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     331      2 FSILWRP 0.35 2.02 2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     332      2 FSILWRP 0.354 2.02 2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     333      2 FRAD 0.38 2.09 2.49 0.0043 0.005612 0 0 0 0 0 0
     334      2 FRAD 0.38 2.09 -2.49 0.0043 0.005612 0 0 0 0 0 0
     335      2 FPRESH 0.39 2.43 -2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
     336      2 FPRESH 0.39 2.43 2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
     337    }
     338
     339}
     340
     341
     342##############
     343# Track merger
     344##############
     345
    193346module Merger TrackMerger {
    194347# add InputArray InputArray
    195   add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    196   add InputArray ElectronMomentumSmearing/electrons
    197   add InputArray MuonMomentumSmearing/muons
     348  add InputArray TrackSmearing/tracks
    198349  set OutputArray tracks
    199350}
    200351
    201352
    202 #############                                                                                                                         
    203 # Calorimeter                                                                                                                                           
    204 #############                                                                                                                                           
     353#############
     354# Calorimeter
     355#############
    205356module DualReadoutCalorimeter Calorimeter {
    206357  set ParticleInputArray ParticlePropagator/stableParticles
     
    224375    set pi [expr {acos(-1)}]
    225376
    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:                                                                                       
     377    # Lists of the edges of each tower in eta and phi;
     378    # each list starts with the lower edge of the first tower;
     379    # the list ends with the higher edged of the last tower.
     380    # Barrel:  deta=0.02 towers up to |eta| <= 0.88 ( up to 45°)
     381    # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad)
     382    # Cell size: about 6 cm x 6 cm
     383
     384    #barrel:
    234385    set PhiBins {}
    235386    for {set i -120} {$i <= 120} {incr i} {
    236387        add PhiBins [expr {$i * $pi/120}]
    237388    }
    238     #deta=0.02 units for |eta| <= 0.88                                                             
     389    #deta=0.02 units for |eta| <= 0.88
    239390    for {set i -44} {$i < 45} {incr i} {
    240391        set eta [expr {$i * 0.02}]
     
    242393    }
    243394
    244     #endcaps:                                                                                       
     395    #endcaps:
    245396    set PhiBins {}
    246397    for {set i -120} {$i <= 120} {incr i} {
    247398        add PhiBins [expr {$i* $pi/120}]
    248399    }
    249     #deta=0.02 units for 0.88 < |eta| <= 3.0                                                       
    250     #first, from -3.0 to -0.88                                                                     
     400    #deta=0.02 units for 0.88 < |eta| <= 3.0
     401    #first, from -3.0 to -0.88
    251402    for {set i 1} {$i <=106} {incr i} {
    252403        set eta [expr {-3.00 + $i * 0.02}]
    253404        add EtaPhiBins $eta $PhiBins
    254405    }
    255     #same for 0.88 to 3.0                                                                           
     406    #same for 0.88 to 3.0
    256407    for  {set i 1} {$i <=106} {incr i} {
    257408        set eta [expr {0.88 + $i * 0.02}]
     
    259410    }
    260411
    261     # default energy fractions {abs(PDG code)} {Fecal Fhcal}                                                                                 
     412    # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    262413    add EnergyFraction {0} {0.0 1.0}
    263     # energy fractions for e, gamma and pi0                                                                                                           
     414    # energy fractions for e, gamma and pi0
    264415    add EnergyFraction {11} {1.0 0.0}
    265416    add EnergyFraction {22} {1.0 0.0}
    266417    add EnergyFraction {111} {1.0 0.0}
    267     # energy fractions for muon, neutrinos and neutralinos                                                                             
     418    # energy fractions for muon, neutrinos and neutralinos
    268419    add EnergyFraction {12} {0.0 0.0}
    269420    add EnergyFraction {13} {0.0 0.0}
     
    275426    add EnergyFraction {1000035} {0.0 0.0}
    276427    add EnergyFraction {1000045} {0.0 0.0}
    277     # energy fractions for K0short and Lambda                                                                                                           
     428    # energy fractions for K0short and Lambda
    278429    add EnergyFraction {310} {0.3 0.7}
    279430    add EnergyFraction {3122} {0.3 0.7}
    280431
    281432
    282     # set ECalResolutionFormula {resolution formula as a function of eta and energy}                               
     433    # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    283434    set ECalResolutionFormula {
    284435    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.11^2)+
     
    286437    }
    287438
    288     # set HCalResolutionFormula {resolution formula as a function of eta and energy}                                               
     439    # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    289440    set HCalResolutionFormula {
    290441    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.30^2)+
     
    352503}
    353504
     505
     506#################
     507# Muon filter
     508#################
     509
     510module PdgCodeFilter MuonFilter {
     511  set InputArray Calorimeter/eflowTracks
     512  set OutputArray muons
     513  set Invert true
     514  add PdgCode {13}
     515  add PdgCode {-13}
     516}
     517
     518
    354519#####################
    355520# Electron efficiency
     
    363528
    364529  # efficiency formula for electrons
    365   set EfficiencyFormula {         
     530  set EfficiencyFormula {
    366531        (energy < 2.0)                                         * (0.000)+
    367532        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
     
    388553}
    389554
     555
    390556#################
    391557# Muon efficiency
     
    393559
    394560module Efficiency MuonEfficiency {
    395   set InputArray MuonMomentumSmearing/muons
     561  set InputArray MuonFilter/muons
    396562  set OutputArray muons
    397563
     
    399565
    400566  # efficiency formula for muons
    401   set EfficiencyFormula {                                   
     567  set EfficiencyFormula {
    402568        (energy < 2.0)                                         * (0.000)+
    403569        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
     
    547713
    548714  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    549  
     715
    550716  # default efficiency formula (misidentification rate)
    551717  add EfficiencyFormula {0} {0.01}
     
    603769module TreeWriter TreeWriter {
    604770    # add Branch InputArray BranchName BranchClass
    605    
     771
    606772    add Branch Delphes/allParticles Particle GenParticle
    607773
    608774    add Branch TrackMerger/tracks Track Track
    609775    add Branch Calorimeter/towers Tower Tower
    610    
     776
    611777    add Branch Calorimeter/eflowTracks EFlowTrack Track
    612778    add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
     
    616782    add Branch PhotonEfficiency/photons PhotonEff Photon
    617783    add Branch PhotonIsolation/photons PhotonIso Photon
    618    
     784
    619785    add Branch GenJetFinder/jets GenJet Jet
    620786    add Branch GenMissingET/momentum GenMissingET MissingET
    621    
     787
    622788    add Branch UniqueObjectFinder/jets Jet Jet
    623789    add Branch UniqueObjectFinder/electrons Electron Electron
    624790    add Branch UniqueObjectFinder/photons Photon Photon
    625791    add Branch UniqueObjectFinder/muons Muon Muon
    626    
    627     add Branch JetEnergyScale/jets AntiKtJet Jet 
    628    
     792
     793    add Branch JetEnergyScale/jets AntiKtJet Jet
     794
    629795    add Branch MissingET/momentum MissingET MissingET
    630796    add Branch ScalarHT/energy ScalarHT ScalarHT
    631797}
    632 
Note: See TracChangeset for help on using the changeset viewer.