Fork me on GitHub

Changeset d612dec in git for cards/delphes_card_IDEA.tcl


Ignore:
Timestamp:
Dec 9, 2021, 7:52:15 AM (3 years ago)
Author:
christinaw97 <christina.wang@…>
Children:
29b722a
Parents:
a5af1df (diff), 0c0c9af (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.
Message:

Merge branch 'master' of github.com:Christinaw97/delphes into HEAD

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_IDEA.tcl

    ra5af1df rd612dec  
     1set RandomSeed 123
     2
    13####################################################################                                l
    2 # FCC-ee IDEA detector model                                                                                     
    3 #                                                                                                   
    4 # Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello                                 
     4# FCC-ee IDEA detector model
     5#
     6# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
    57# email: efontane@bo.infn.it,
    6 #        lorenzo.pezzotti01@universitadipavia.it,                                                   
    7 #        m.antonello@uninsubria.it,                                                                 
    8 #####################################################################                               
    9 #       
     8#        lorenzo.pezzotti01@universitadipavia.it,
     9#        m.antonello@uninsubria.it,
     10#        michele.selvaggi@cern.ch
     11#####################################################################
     12
     13## MOD2: set vtx mode timing to MC truth
     14
     15set B 2.0
     16
     17## Drift chamber coordinates
     18set DCHZMIN -2.125
     19set DCHZMAX 2.125
     20set DCHRMIN 0.345
     21set DCHRMAX 2.02
     22
     23
    1024#######################################
    1125# Order of execution of various modules
    … …  
    1327
    1428set ExecutionPath {
     29
     30  TruthVertexFinder
    1531  ParticlePropagator
    1632
    … …  
    1935  MuonTrackingEfficiency
    2036
    21   ChargedHadronMomentumSmearing
    22   ElectronMomentumSmearing
    23   MuonMomentumSmearing
    24 
    25   TrackMerger
     37  TrackMergerPre
     38  TrackSmearing
     39  ClusterCounting
     40  TimeSmearing
     41  TimeOfFlight
     42
     43  TrackMerger
     44  ForwardLooperTracks
    2645  Calorimeter
     46
     47  TimeSmearingNeutrals
     48  TimeOfFlightNeutralHadron
     49
     50  EFlowTrackMerger
    2751  EFlowMerger
    2852
    2953  PhotonEfficiency
    3054  PhotonIsolation
     55
     56  MuonFilter
    3157
    3258  ElectronFilter
    … …  
    4268  GenJetFinder
    4369  GenMissingET
    44  
     70
    4571  FastJetFinder
    4672
    … …  
    5278  TauTagging
    5379
    54   UniqueObjectFinder
    55 
    56   ScalarHT
    5780  TreeWriter
    5881}
    5982
     83#################################
     84# Truth Vertex Finder
     85#################################
     86
     87module TruthVertexFinder TruthVertexFinder {
     88
     89  ## below this distance two vertices are assumed to be merged
     90  set Resolution 1E-06
     91
     92  set InputArray Delphes/stableParticles
     93  set VertexOutputArray vertices
     94}
    6095
    6196#################################
    … …  
    78113
    79114  # magnetic field, in T
    80   set Bz 2.0
     115  set Bz $B
    81116}
    82117
    … …  
    88123    set InputArray ParticlePropagator/chargedHadrons
    89124    set OutputArray chargedHadrons
    90     # We use only one efficiency, we set only 0 effincency out of eta bounds:
     125    set UseMomentumVector true
    91126
    92127    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)
     128        (abs(eta) > 2.56)                                  * (0.000) +
     129        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000) +
     130        (pt >= 0.1) * (abs(eta) <= 2.56)                   * (1.000)
     131    }
     132}
    103133
    104134
    … …  
    111141    set InputArray ParticlePropagator/electrons
    112142    set OutputArray electrons
    113 
    114 
    115     # Current full simulation with CLICdet provides for electrons:
     143    set UseMomentumVector true
     144
    116145    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     }
     146        (abs(eta) > 2.56)                                  * (0.000) +
     147        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000) +
     148        (pt >= 0.1) * (abs(eta) <= 2.56)                   * (1.000)
     149    }
    122150}
    123151
    … …  
    130158    set InputArray ParticlePropagator/muons
    131159    set OutputArray muons
    132 
    133     # Current full simulation with CLICdet provides for muons:
     160    set UseMomentumVector true
     161
    134162    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 
    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.0006^2 + (pt*7.e-5)^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.0006^2 + (pt*7.e-5)^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.0006^2 + (pt*7.e-5)^2)
     163        (abs(eta) > 2.56)                                  * (0.000) +
     164        (pt < 0.1) * (abs(eta) <= 2.56)                    * (0.000) +
     165        (pt >= 0.1) * (abs(eta) <= 2.56)                   * (1.000)
    186166    }
    187167}
    … …  
    191171##############
    192172
     173module Merger TrackMergerPre {
     174# add InputArray InputArray
     175  add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
     176  add InputArray ElectronTrackingEfficiency/electrons
     177  add InputArray MuonTrackingEfficiency/muons
     178  set OutputArray tracks
     179}
     180
     181
     182
     183########################################
     184# Smearing for charged tracks
     185########################################
     186
     187module TrackCovariance TrackSmearing {
     188
     189    set InputArray TrackMergerPre/tracks
     190    set OutputArray tracks
     191
     192    ## minimum number of hits to accept a track
     193    set NMinHits 6
     194
     195    ## magnetic field
     196    set Bz $B
     197
     198    ## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt
     199    set DetectorGeometry {
     200
     201
     202      # Layer type 1 = R (barrel) or 2 = z (forward/backward)
     203      # Layer label
     204      # Minimum dimension z for barrel or R for forward
     205      # Maximum dimension z for barrel or R for forward
     206      # R/z location of layer
     207      # Thickness (meters)
     208      # Radiation length (meters)
     209      # Number of measurements in layers (1D or 2D)
     210      # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Upper side
     211      # Stereo angle (rad) - 0(pi/2) = axial(z) layer - Lower side
     212      # Resolution Upper side (meters) - 0 = no measurement
     213      # Resolution Lower side (meters) - 0 = no measurement
     214      # measurement flag = T, scattering only = F
     215
     216      # barrel  name       zmin   zmax   r        w (m)      X0        n_meas  th_up (rad) th_down (rad)    reso_up (m)   reso_down (m)  flag
     217
     218      # barrel  name       zmin   zmax   r        w (m)      X0        n_meas  th_up (rad) th_down (rad)    reso_up (m)   reso_down (m)  flag
     219
     220      1        PIPE       -100    100    0.015    0.001655  0.2805     0        0          0                0             0              0
     221      1        VTXLOW     -0.12   0.12   0.017    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     222      1        VTXLOW     -0.16   0.16   0.023    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     223      1        VTXLOW     -0.16   0.16   0.031    0.00028   0.0937     2        0          1.5708           3e-006        3e-006         1
     224      1        VTXHIGH    -1      1      0.32     0.00047   0.0937     2        0          1.5708           7e-006        7e-006         1
     225      1        VTXHIGH    -1.05   1.05   0.34     0.00047   0.0937     2        0          1.5708           7e-006        7e-006         1
     226
     227      # endcap  name       rmin   rmax   z        w (m)      X0        n_meas   th_up (rad)  th_down (rad)   reso_up (m)   reso_down (m) flag
     228
     229      2        VTXDSK      0.141  0.3   -0.92     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     230      2        VTXDSK      0.138  0.3   -0.9      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     231      2        VTXDSK      0.065  0.3   -0.42     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     232      2        VTXDSK      0.062  0.3   -0.4      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     233      2        VTXDSK      0.062  0.3    0.4      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     234      2        VTXDSK      0.065  0.3    0.42     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     235      2        VTXDSK      0.138  0.3    0.9      0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     236      2        VTXDSK      0.141  0.3    0.92     0.00028   0.0937     2        0          1.5708           7e-006        7e-006         1
     237
     238      1 DCHCANI $DCHZMIN $DCHZMAX $DCHRMIN 0.0002 0.237223 0 0 0 0 0 0
     239      1 DCH -2 2 0.36 0.0147748 1400 1 0.0203738 0 0.0001 0 1
     240      1 DCH -2 2 0.374775 0.0147748 1400 1 -0.0212097 0 0.0001 0 1
     241      1 DCH -2 2 0.38955 0.0147748 1400 1 0.0220456 0 0.0001 0 1
     242      1 DCH -2 2 0.404324 0.0147748 1400 1 -0.0228814 0 0.0001 0 1
     243      1 DCH -2 2 0.419099 0.0147748 1400 1 0.0237172 0 0.0001 0 1
     244      1 DCH -2 2 0.433874 0.0147748 1400 1 -0.024553 0 0.0001 0 1
     245      1 DCH -2 2 0.448649 0.0147748 1400 1 0.0253888 0 0.0001 0 1
     246      1 DCH -2 2 0.463423 0.0147748 1400 1 -0.0262245 0 0.0001 0 1
     247      1 DCH -2 2 0.478198 0.0147748 1400 1 0.0270602 0 0.0001 0 1
     248      1 DCH -2 2 0.492973 0.0147748 1400 1 -0.0278958 0 0.0001 0 1
     249      1 DCH -2 2 0.507748 0.0147748 1400 1 0.0287314 0 0.0001 0 1
     250      1 DCH -2 2 0.522523 0.0147748 1400 1 -0.029567 0 0.0001 0 1
     251      1 DCH -2 2 0.537297 0.0147748 1400 1 0.0304025 0 0.0001 0 1
     252      1 DCH -2 2 0.552072 0.0147748 1400 1 -0.031238 0 0.0001 0 1
     253      1 DCH -2 2 0.566847 0.0147748 1400 1 0.0320734 0 0.0001 0 1
     254      1 DCH -2 2 0.581622 0.0147748 1400 1 -0.0329088 0 0.0001 0 1
     255      1 DCH -2 2 0.596396 0.0147748 1400 1 0.0337442 0 0.0001 0 1
     256      1 DCH -2 2 0.611171 0.0147748 1400 1 -0.0345795 0 0.0001 0 1
     257      1 DCH -2 2 0.625946 0.0147748 1400 1 0.0354147 0 0.0001 0 1
     258      1 DCH -2 2 0.640721 0.0147748 1400 1 -0.0362499 0 0.0001 0 1
     259      1 DCH -2 2 0.655495 0.0147748 1400 1 0.0370851 0 0.0001 0 1
     260      1 DCH -2 2 0.67027 0.0147748 1400 1 -0.0379202 0 0.0001 0 1
     261      1 DCH -2 2 0.685045 0.0147748 1400 1 0.0387552 0 0.0001 0 1
     262      1 DCH -2 2 0.69982 0.0147748 1400 1 -0.0395902 0 0.0001 0 1
     263      1 DCH -2 2 0.714595 0.0147748 1400 1 0.0404252 0 0.0001 0 1
     264      1 DCH -2 2 0.729369 0.0147748 1400 1 -0.04126 0 0.0001 0 1
     265      1 DCH -2 2 0.744144 0.0147748 1400 1 0.0420949 0 0.0001 0 1
     266      1 DCH -2 2 0.758919 0.0147748 1400 1 -0.0429296 0 0.0001 0 1
     267      1 DCH -2 2 0.773694 0.0147748 1400 1 0.0437643 0 0.0001 0 1
     268      1 DCH -2 2 0.788468 0.0147748 1400 1 -0.044599 0 0.0001 0 1
     269      1 DCH -2 2 0.803243 0.0147748 1400 1 0.0454336 0 0.0001 0 1
     270      1 DCH -2 2 0.818018 0.0147748 1400 1 -0.0462681 0 0.0001 0 1
     271      1 DCH -2 2 0.832793 0.0147748 1400 1 0.0471025 0 0.0001 0 1
     272      1 DCH -2 2 0.847568 0.0147748 1400 1 -0.0479369 0 0.0001 0 1
     273      1 DCH -2 2 0.862342 0.0147748 1400 1 0.0487713 0 0.0001 0 1
     274      1 DCH -2 2 0.877117 0.0147748 1400 1 -0.0496055 0 0.0001 0 1
     275      1 DCH -2 2 0.891892 0.0147748 1400 1 0.0504397 0 0.0001 0 1
     276      1 DCH -2 2 0.906667 0.0147748 1400 1 -0.0512738 0 0.0001 0 1
     277      1 DCH -2 2 0.921441 0.0147748 1400 1 0.0521079 0 0.0001 0 1
     278      1 DCH -2 2 0.936216 0.0147748 1400 1 -0.0529418 0 0.0001 0 1
     279      1 DCH -2 2 0.950991 0.0147748 1400 1 0.0537757 0 0.0001 0 1
     280      1 DCH -2 2 0.965766 0.0147748 1400 1 -0.0546095 0 0.0001 0 1
     281      1 DCH -2 2 0.980541 0.0147748 1400 1 0.0554433 0 0.0001 0 1
     282      1 DCH -2 2 0.995315 0.0147748 1400 1 -0.056277 0 0.0001 0 1
     283      1 DCH -2 2 1.01009 0.0147748 1400 1 0.0571106 0 0.0001 0 1
     284      1 DCH -2 2 1.02486 0.0147748 1400 1 -0.0579441 0 0.0001 0 1
     285      1 DCH -2 2 1.03964 0.0147748 1400 1 0.0587775 0 0.0001 0 1
     286      1 DCH -2 2 1.05441 0.0147748 1400 1 -0.0596108 0 0.0001 0 1
     287      1 DCH -2 2 1.06919 0.0147748 1400 1 0.0604441 0 0.0001 0 1
     288      1 DCH -2 2 1.08396 0.0147748 1400 1 -0.0612773 0 0.0001 0 1
     289      1 DCH -2 2 1.09874 0.0147748 1400 1 0.0621104 0 0.0001 0 1
     290      1 DCH -2 2 1.11351 0.0147748 1400 1 -0.0629434 0 0.0001 0 1
     291      1 DCH -2 2 1.12829 0.0147748 1400 1 0.0637763 0 0.0001 0 1
     292      1 DCH -2 2 1.14306 0.0147748 1400 1 -0.0646092 0 0.0001 0 1
     293      1 DCH -2 2 1.15784 0.0147748 1400 1 0.0654419 0 0.0001 0 1
     294      1 DCH -2 2 1.17261 0.0147748 1400 1 -0.0662746 0 0.0001 0 1
     295      1 DCH -2 2 1.18739 0.0147748 1400 1 0.0671071 0 0.0001 0 1
     296      1 DCH -2 2 1.20216 0.0147748 1400 1 -0.0679396 0 0.0001 0 1
     297      1 DCH -2 2 1.21694 0.0147748 1400 1 0.068772 0 0.0001 0 1
     298      1 DCH -2 2 1.23171 0.0147748 1400 1 -0.0696042 0 0.0001 0 1
     299      1 DCH -2 2 1.24649 0.0147748 1400 1 0.0704364 0 0.0001 0 1
     300      1 DCH -2 2 1.26126 0.0147748 1400 1 -0.0712685 0 0.0001 0 1
     301      1 DCH -2 2 1.27604 0.0147748 1400 1 0.0721005 0 0.0001 0 1
     302      1 DCH -2 2 1.29081 0.0147748 1400 1 -0.0729324 0 0.0001 0 1
     303      1 DCH -2 2 1.30559 0.0147748 1400 1 0.0737642 0 0.0001 0 1
     304      1 DCH -2 2 1.32036 0.0147748 1400 1 -0.0745958 0 0.0001 0 1
     305      1 DCH -2 2 1.33514 0.0147748 1400 1 0.0754274 0 0.0001 0 1
     306      1 DCH -2 2 1.34991 0.0147748 1400 1 -0.0762589 0 0.0001 0 1
     307      1 DCH -2 2 1.36468 0.0147748 1400 1 0.0770903 0 0.0001 0 1
     308      1 DCH -2 2 1.37946 0.0147748 1400 1 -0.0779215 0 0.0001 0 1
     309      1 DCH -2 2 1.39423 0.0147748 1400 1 0.0787527 0 0.0001 0 1
     310      1 DCH -2 2 1.40901 0.0147748 1400 1 -0.0795837 0 0.0001 0 1
     311      1 DCH -2 2 1.42378 0.0147748 1400 1 0.0804147 0 0.0001 0 1
     312      1 DCH -2 2 1.43856 0.0147748 1400 1 -0.0812455 0 0.0001 0 1
     313      1 DCH -2 2 1.45333 0.0147748 1400 1 0.0820762 0 0.0001 0 1
     314      1 DCH -2 2 1.46811 0.0147748 1400 1 -0.0829068 0 0.0001 0 1
     315      1 DCH -2 2 1.48288 0.0147748 1400 1 0.0837373 0 0.0001 0 1
     316      1 DCH -2 2 1.49766 0.0147748 1400 1 -0.0845677 0 0.0001 0 1
     317      1 DCH -2 2 1.51243 0.0147748 1400 1 0.0853979 0 0.0001 0 1
     318      1 DCH -2 2 1.52721 0.0147748 1400 1 -0.086228 0 0.0001 0 1
     319      1 DCH -2 2 1.54198 0.0147748 1400 1 0.087058 0 0.0001 0 1
     320      1 DCH -2 2 1.55676 0.0147748 1400 1 -0.0878879 0 0.0001 0 1
     321      1 DCH -2 2 1.57153 0.0147748 1400 1 0.0887177 0 0.0001 0 1
     322      1 DCH -2 2 1.58631 0.0147748 1400 1 -0.0895474 0 0.0001 0 1
     323      1 DCH -2 2 1.60108 0.0147748 1400 1 0.0903769 0 0.0001 0 1
     324      1 DCH -2 2 1.61586 0.0147748 1400 1 -0.0912063 0 0.0001 0 1
     325      1 DCH -2 2 1.63063 0.0147748 1400 1 0.0920356 0 0.0001 0 1
     326      1 DCH -2 2 1.64541 0.0147748 1400 1 -0.0928647 0 0.0001 0 1
     327      1 DCH -2 2 1.66018 0.0147748 1400 1 0.0936937 0 0.0001 0 1
     328      1 DCH -2 2 1.67495 0.0147748 1400 1 -0.0945226 0 0.0001 0 1
     329      1 DCH -2 2 1.68973 0.0147748 1400 1 0.0953514 0 0.0001 0 1
     330      1 DCH -2 2 1.7045 0.0147748 1400 1 -0.09618 0 0.0001 0 1
     331      1 DCH -2 2 1.71928 0.0147748 1400 1 0.0970085 0 0.0001 0 1
     332      1 DCH -2 2 1.73405 0.0147748 1400 1 -0.0978369 0 0.0001 0 1
     333      1 DCH -2 2 1.74883 0.0147748 1400 1 0.0986651 0 0.0001 0 1
     334      1 DCH -2 2 1.7636 0.0147748 1400 1 -0.0994932 0 0.0001 0 1
     335      1 DCH -2 2 1.77838 0.0147748 1400 1 0.100321 0 0.0001 0 1
     336      1 DCH -2 2 1.79315 0.0147748 1400 1 -0.101149 0 0.0001 0 1
     337      1 DCH -2 2 1.80793 0.0147748 1400 1 0.101977 0 0.0001 0 1
     338      1 DCH -2 2 1.8227 0.0147748 1400 1 -0.102804 0 0.0001 0 1
     339      1 DCH -2 2 1.83748 0.0147748 1400 1 0.103632 0 0.0001 0 1
     340      1 DCH -2 2 1.85225 0.0147748 1400 1 -0.104459 0 0.0001 0 1
     341      1 DCH -2 2 1.86703 0.0147748 1400 1 0.105286 0 0.0001 0 1
     342      1 DCH -2 2 1.8818 0.0147748 1400 1 -0.106113 0 0.0001 0 1
     343      1 DCH -2 2 1.89658 0.0147748 1400 1 0.10694 0 0.0001 0 1
     344      1 DCH -2 2 1.91135 0.0147748 1400 1 -0.107766 0 0.0001 0 1
     345      1 DCH -2 2 1.92613 0.0147748 1400 1 0.108593 0 0.0001 0 1
     346      1 DCH -2 2 1.9409 0.0147748 1400 1 -0.109419 0 0.0001 0 1
     347      1 DCH -2 2 1.95568 0.0147748 1400 1 0.110246 0 0.0001 0 1
     348      1 DCH -2 2 1.97045 0.0147748 1400 1 -0.111072 0 0.0001 0 1
     349      1 DCH -2 2 1.98523 0.0147748 1400 1 0.111898 0 0.0001 0 1
     350      1 DCH -2 2 2 0.0147748 1400 1 -0.112723 0 0.0001 0 1
     351      1 DCHCANO $DCHZMIN $DCHZMAX $DCHRMAX $DCHRMAX 0.02 1.667 0 0 0 0 0 0
     352      1 BSILWRP -2.35 2.35 2.04 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     353      1 BSILWRP -2.35 2.35 2.06 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     354      1 MAG -2.5 2.5 2.25 0.05 0.0658 0 0 0 0 0 0
     355      1 BPRESH -2.55 2.55 2.45 0.02 1 2 0 1.5708 7e-005 0.01 1
     356      2 DCHWALL $DCHRMIN $DCHRMAX $DCHZMAX 0.25 5.55 0 0 0 0 0 0
     357      2 DCHWALL $DCHRMIN $DCHRMAX $DCHZMIN 0.25 5.55 0 0 0 0 0 0
     358      2 FSILWRP 0.354 2.02 -2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     359      2 FSILWRP 0.35 2.02 -2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     360      2 FSILWRP 0.35 2.02 2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     361      2 FSILWRP 0.354 2.02 2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
     362      2 FRAD 0.38 2.09 2.49 0.0043 0.005612 0 0 0 0 0 0
     363      2 FRAD 0.38 2.09 -2.49 0.0043 0.005612 0 0 0 0 0 0
     364      2 FPRESH 0.39 2.43 -2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
     365      2 FPRESH 0.39 2.43 2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
     366    }
     367
     368}
     369
     370###################
     371# Cluster Counting
     372###################
     373
     374module ClusterCounting ClusterCounting {
     375
     376  add InputArray TrackSmearing/tracks
     377  set OutputArray tracks
     378
     379  set Bz $B
     380
     381  ## check that these are consistent with DCHCANI/DCHNANO parameters in TrackCovariance module
     382  set Rmin $DCHRMIN
     383  set Rmax $DCHRMAX
     384  set Zmin $DCHZMIN
     385  set Zmax $DCHZMAX
     386
     387  # gas mix option:
     388  # 0:  Helium 90% - Isobutane 10%
     389  # 1:  Helium 100%
     390  # 2:  Argon 50% - Ethane 50%
     391  # 3:  Argon 100%
     392
     393  set GasOption 0
     394
     395}
     396
     397
     398########################################
     399#   Time Smearing MIP
     400########################################
     401
     402module TimeSmearing TimeSmearing {
     403  set InputArray ClusterCounting/tracks
     404  set OutputArray tracks
     405
     406  # assume constant 30 ps resolution for now
     407  set TimeResolution {
     408                       (abs(eta) > 0.0 && abs(eta) <= 3.0)* 30E-12
     409                     }
     410}
     411
     412
     413########################################
     414#   Time Of Flight Measurement
     415########################################
     416
     417module TimeOfFlight TimeOfFlight {
     418  set InputArray TimeSmearing/tracks
     419  set VertexInputArray TruthVertexFinder/vertices
     420
     421  set OutputArray tracks
     422
     423  # 0: assume vertex time tV from MC Truth (ideal case)
     424  # 1: assume vertex time tV = 0
     425  # 2: calculate vertex time as vertex TOF, assuming tPV=0
     426
     427  set VertexTimeMode 0
     428}
     429
     430##############
     431# Track merger
     432##############
     433
    193434module Merger TrackMerger {
    194435# add InputArray InputArray
    195   add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    196   add InputArray ElectronMomentumSmearing/electrons
    197   add InputArray MuonMomentumSmearing/muons
     436  add InputArray TimeOfFlight/tracks
    198437  set OutputArray tracks
    199438}
    200439
    201440
    202 #############                                                                                                                         
    203 # Calorimeter                                                                                                                                           
    204 #############                                                                                                                                           
     441######################
     442# Looper Selection
     443######################
     444
     445module Efficiency ForwardLooperTracks  {
     446  set InputArray TrackMerger/tracks
     447  set OutputArray tracks
     448  set UseMomentumVector False
     449
     450  ## select looping tracks that end up in position |eta| > 3.142 (lost by calo)
     451  set EfficiencyFormula {
     452    (abs(eta) > 3.0 )                                 * (1.000) +
     453    (abs(eta) <= 3.0 )                                * (0.000)
     454  }
     455
     456}
     457
     458
     459#############
     460# Calorimeter
     461#############
    205462module DualReadoutCalorimeter Calorimeter {
    206463  set ParticleInputArray ParticlePropagator/stableParticles
    … …  
    216473  set ECalEnergyMin 0.5
    217474  set HCalEnergyMin 0.5
     475  set ECalEnergySignificanceMin 3.0
     476  set HCalEnergySignificanceMin 3.0
     477
    218478  set EnergyMin 0.5
    219   set ECalEnergySignificanceMin 1.0
    220   set HCalEnergySignificanceMin 1.0
    221   set EnergySignificanceMin 1.0
    222 
    223   set SmearTowerCenter true
     479  set EnergySignificanceMin 3.0
     480
     481  #set SmearTowerCenter true
     482  set SmearTowerCenter false
    224483    set pi [expr {acos(-1)}]
    225484
    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:                                                                                       
     485    # Lists of the edges of each tower in eta and phi;
     486    # each list starts with the lower edge of the first tower;
     487    # the list ends with the higher edged of the last tower.
     488    # Barrel:  deta=0.02 towers up to |eta| <= 0.88 ( up to 45°)
     489    # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad)
     490    # Cell size: about 6 cm x 6 cm
     491
     492    #barrel:
    234493    set PhiBins {}
    235494    for {set i -120} {$i <= 120} {incr i} {
    236495        add PhiBins [expr {$i * $pi/120}]
    237496    }
    238     #deta=0.02 units for |eta| <= 0.88                                                             
     497    #deta=0.02 units for |eta| <= 0.88
    239498    for {set i -44} {$i < 45} {incr i} {
    240499        set eta [expr {$i * 0.02}]
    … …  
    242501    }
    243502
    244     #endcaps:                                                                                       
     503    #endcaps:
    245504    set PhiBins {}
    246505    for {set i -120} {$i <= 120} {incr i} {
    247506        add PhiBins [expr {$i* $pi/120}]
    248507    }
    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} {
     508    #deta=0.02 units for 0.88 < |eta| <= 3.0
     509    #first, from -3.0 to -0.88
     510    for {set i 0} {$i <=106} {incr i} {
    252511        set eta [expr {-3.00 + $i * 0.02}]
    253512        add EtaPhiBins $eta $PhiBins
    254513    }
    255     #same for 0.88 to 3.0                                                                           
     514    #same for 0.88 to 3.0
    256515    for  {set i 1} {$i <=106} {incr i} {
    257516        set eta [expr {0.88 + $i * 0.02}]
    … …  
    259518    }
    260519
    261     # default energy fractions {abs(PDG code)} {Fecal Fhcal}                                                                                 
     520    # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    262521    add EnergyFraction {0} {0.0 1.0}
    263     # energy fractions for e, gamma and pi0                                                                                                           
     522    # energy fractions for e, gamma and pi0
    264523    add EnergyFraction {11} {1.0 0.0}
    265524    add EnergyFraction {22} {1.0 0.0}
    266525    add EnergyFraction {111} {1.0 0.0}
    267     # energy fractions for muon, neutrinos and neutralinos                                                                             
     526    # energy fractions for muon, neutrinos and neutralinos
    268527    add EnergyFraction {12} {0.0 0.0}
    269528    add EnergyFraction {13} {0.0 0.0}
    … …  
    275534    add EnergyFraction {1000035} {0.0 0.0}
    276535    add EnergyFraction {1000045} {0.0 0.0}
    277     # energy fractions for K0short and Lambda                                                                                                           
     536    # energy fractions for K0short and Lambda
    278537    add EnergyFraction {310} {0.3 0.7}
     538    add EnergyFraction {130} {0.3 0.7}
    279539    add EnergyFraction {3122} {0.3 0.7}
    280540
    281541
    282     # set ECalResolutionFormula {resolution formula as a function of eta and energy}                               
     542    # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    283543    set ECalResolutionFormula {
    284544    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.11^2)+
    … …  
    286546    }
    287547
    288     # set HCalResolutionFormula {resolution formula as a function of eta and energy}                                               
     548    # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    289549    set HCalResolutionFormula {
    290550    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.30^2)+
    … …  
    293553}
    294554
     555########################################
     556#   Time Smearing Neutrals
     557########################################
     558
     559module TimeSmearing TimeSmearingNeutrals {
     560  set InputArray Calorimeter/eflowNeutralHadrons
     561  set OutputArray eflowNeutralHadrons
     562
     563  # assume constant 30 ps resolution for now
     564  set TimeResolution {
     565                       (abs(eta) > 0.0 && abs(eta) <= 3.0)* 30E-12
     566                     }
     567}
     568
     569########################################
     570#   Time Of Flight Measurement
     571########################################
     572
     573module TimeOfFlight TimeOfFlightNeutralHadron {
     574  set InputArray TimeSmearingNeutrals/eflowNeutralHadrons
     575  set VertexInputArray TruthVertexFinder/vertices
     576
     577  set OutputArray eflowNeutralHadrons
     578
     579  # 0: assume vertex time tV from MC Truth (ideal case)
     580  # 1: assume vertex time tV = 0
     581  # 2: calculate vertex time as vertex TOF, assuming tPV=0
     582
     583  ## TBF (add option to take hard vertex time)
     584  set VertexTimeMode 1
     585}
     586
     587
     588############################
     589# Energy flow track merger
     590############################
     591
     592module Merger EFlowTrackMerger {
     593# add InputArray InputArray
     594  add InputArray Calorimeter/eflowTracks
     595  add InputArray ForwardLooperTracks/tracks
     596  set OutputArray eflowTracks
     597}
     598
     599
     600
    295601####################
    296602# Energy flow merger
    … …  
    299605module Merger EFlowMerger {
    300606# add InputArray InputArray
    301   add InputArray Calorimeter/eflowTracks
     607  add InputArray EFlowTrackMerger/eflowTracks
    302608  add InputArray Calorimeter/eflowPhotons
    303   add InputArray Calorimeter/eflowNeutralHadrons
     609  add InputArray TimeOfFlightNeutralHadron/eflowNeutralHadrons
    304610  set OutputArray eflow
    305611}
     612
    306613
    307614###################
    … …  
    337644  set PTMin 0.5
    338645
    339   set PTRatioMax 999.
     646  set PTRatioMax 9999.
    340647}
    341648
    … …  
    345652
    346653module PdgCodeFilter ElectronFilter {
    347   set InputArray Calorimeter/eflowTracks
     654  set InputArray EFlowTrackMerger/eflowTracks
    348655  set OutputArray electrons
    349656  set Invert true
    … …  
    352659}
    353660
     661#################
     662# Muon filter
     663#################
     664
     665module PdgCodeFilter MuonFilter {
     666  set InputArray EFlowTrackMerger/eflowTracks
     667  set OutputArray muons
     668  set Invert true
     669  add PdgCode {13}
     670  add PdgCode {-13}
     671}
     672
     673
    354674#####################
    355675# Electron efficiency
    … …  
    363683
    364684  # efficiency formula for electrons
    365   set EfficiencyFormula {         
     685  set EfficiencyFormula {
    366686        (energy < 2.0)                                         * (0.000)+
    367687        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
    … …  
    385705  set PTMin 0.5
    386706
    387   set PTRatioMax 0.12
     707  set PTRatioMax 9999
    388708}
    389709
    … …  
    393713
    394714module Efficiency MuonEfficiency {
    395   set InputArray MuonMomentumSmearing/muons
     715  set InputArray MuonFilter/muons
    396716  set OutputArray muons
    397717
    … …  
    399719
    400720  # efficiency formula for muons
    401   set EfficiencyFormula {                                   
     721  set EfficiencyFormula {
    402722        (energy < 2.0)                                         * (0.000)+
    403723        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
    … …  
    421741  set PTMin 0.5
    422742
    423   set PTRatioMax 0.25
     743  set PTRatioMax 9999.
    424744}
    425745
    … …  
    473793module FastJetFinder GenJetFinder {
    474794  set InputArray NeutrinoFilter/filteredParticles
    475 
    476795  set OutputArray jets
    477796
    478   # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    479   set JetAlgorithm 6
    480   set ParameterR 0.4
     797  set JetAlgorithm 10
     798  set ParameterR 1.5
     799  set ParameterP -1.0
    481800  set JetPTMin 1.0
     801
    482802}
    483803
    … …  
    504824
    505825  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    506   set JetAlgorithm 6
    507   set ParameterR 0.4
     826  set JetAlgorithm 10
     827  set ParameterR 1.5
     828  set ParameterP -1.0
    508829  set JetPTMin 1.0
     830
     831
    509832}
    510833
    … …  
    518841
    519842  # scale formula for jets
    520   set ScaleFormula {1.08}
     843  set ScaleFormula {1.00}
    521844}
    522845
    … …  
    547870
    548871  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    549  
     872
    550873  # default efficiency formula (misidentification rate)
    551874  add EfficiencyFormula {0} {0.01}
    … …  
    578901
    579902
    580 #####################################################
    581 # Find uniquely identified photons/electrons/tau/jets
    582 #####################################################
    583 
    584 module UniqueObjectFinder UniqueObjectFinder {
    585 # earlier arrays take precedence over later ones
    586 # add InputArray InputArray OutputArray
    587   add InputArray PhotonIsolation/photons photons
    588   add InputArray ElectronIsolation/electrons electrons
    589   add InputArray MuonIsolation/muons muons
    590   add InputArray JetEnergyScale/jets jets
    591 }
    592 
    593 
    594903
    595904##################
    … …  
    603912module TreeWriter TreeWriter {
    604913    # add Branch InputArray BranchName BranchClass
    605    
     914
    606915    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
     916    add Branch TruthVertexFinder/vertices GenVertex Vertex
     917
     918    add Branch EFlowTrackMerger/eflowTracks EFlowTrack Track
     919    add Branch TrackSmearing/tracks Track Track
    612920    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    
     921    add Branch TimeOfFlightNeutralHadron/eflowNeutralHadrons EFlowNeutralHadron Tower
     922
     923    add Branch EFlowMerger/eflow ParticleFlowCandidate ParticleFlowCandidate
     924
     925    add Branch ElectronEfficiency/electrons Electron Electron
     926    add Branch MuonEfficiency/muons Muon Muon
     927    add Branch PhotonEfficiency/photons Photon Photon
     928
     929    add Branch JetEnergyScale/jets Jet Jet
     930    add Branch MissingET/momentum MissingET MissingET
     931
    619932    add Branch GenJetFinder/jets GenJet Jet
    620933    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 
     934
     935    # add Info InfoName InfoValue
     936    add Info Bz $B
     937}
Note: See TracChangeset for help on using the changeset viewer.