| 1 | #
|
|---|
| 2 | # Official Delphes card prepared by FCC-hh collaboration
|
|---|
| 3 | #
|
|---|
| 4 | # Main authors: Heather Gray (CERN)
|
|---|
| 5 | # Filip Moortgat (CERN)
|
|---|
| 6 | #
|
|---|
| 7 | # Released on: Sep 15th, 2015
|
|---|
| 8 | #
|
|---|
| 9 | # Version: v00
|
|---|
| 10 | #
|
|---|
| 11 | # Configuration: FCC-hh generic detector with dipole region, no pile-up
|
|---|
| 12 | #
|
|---|
| 13 |
|
|---|
| 14 | #######################################
|
|---|
| 15 | # Order of execution of various modules
|
|---|
| 16 | #######################################
|
|---|
| 17 |
|
|---|
| 18 | set ExecutionPath {
|
|---|
| 19 | PileUpMerger
|
|---|
| 20 | ParticlePropagator
|
|---|
| 21 |
|
|---|
| 22 | ChargedHadronTrackingEfficiency
|
|---|
| 23 |
|
|---|
| 24 | ChargedHadronMomentumSmearing
|
|---|
| 25 |
|
|---|
| 26 |
|
|---|
| 27 | ElectronEfficiency
|
|---|
| 28 | MuonEfficiency
|
|---|
| 29 |
|
|---|
| 30 | ChargedHadronMomentumSmearing
|
|---|
| 31 | ElectronEnergySmearing
|
|---|
| 32 | MuonMomentumSmearing
|
|---|
| 33 |
|
|---|
| 34 | TrackMerger
|
|---|
| 35 | ImpactParameterSmearing
|
|---|
| 36 |
|
|---|
| 37 | Ecal
|
|---|
| 38 | Hcal
|
|---|
| 39 | TrackPileUpSubtractor
|
|---|
| 40 |
|
|---|
| 41 | NeutralTowerMerger
|
|---|
| 42 | TowerMerger
|
|---|
| 43 | EFlowMergerAllTracks
|
|---|
| 44 | EFlowMerger
|
|---|
| 45 |
|
|---|
| 46 | NeutrinoFilter
|
|---|
| 47 | GenJetFinder
|
|---|
| 48 | GenMissingET
|
|---|
| 49 |
|
|---|
| 50 | Rho
|
|---|
| 51 | FastJetFinder
|
|---|
| 52 | PileUpJetID
|
|---|
| 53 | JetPileUpSubtractor
|
|---|
| 54 |
|
|---|
| 55 | JetEnergyScale
|
|---|
| 56 |
|
|---|
| 57 | MissingET
|
|---|
| 58 |
|
|---|
| 59 | PhotonEfficiency
|
|---|
| 60 | PhotonIsolation
|
|---|
| 61 |
|
|---|
| 62 | ElectronIsolation
|
|---|
| 63 |
|
|---|
| 64 | MuonIsolation
|
|---|
| 65 |
|
|---|
| 66 | JetFlavorAssociation
|
|---|
| 67 |
|
|---|
| 68 | BTagging
|
|---|
| 69 | CTagging
|
|---|
| 70 | TauTagging
|
|---|
| 71 |
|
|---|
| 72 | UniqueObjectFinder
|
|---|
| 73 |
|
|---|
| 74 | ScalarHT
|
|---|
| 75 |
|
|---|
| 76 | TreeWriter
|
|---|
| 77 | }
|
|---|
| 78 |
|
|---|
| 79 | ###############
|
|---|
| 80 | # PileUp Merger
|
|---|
| 81 | ###############
|
|---|
| 82 |
|
|---|
| 83 | module PileUpMerger PileUpMerger {
|
|---|
| 84 | set InputArray Delphes/stableParticles
|
|---|
| 85 |
|
|---|
| 86 | set ParticleOutputArray stableParticles
|
|---|
| 87 | set VertexOutputArray vertices
|
|---|
| 88 |
|
|---|
| 89 | # pre-generated minbias input file
|
|---|
| 90 | set PileUpFile /data02/testa/Delphes/Delphes-3.3.2/MinBias.pileup
|
|---|
| 91 |
|
|---|
| 92 | # average expected pile up
|
|---|
| 93 | set MeanPileUp 50
|
|---|
| 94 |
|
|---|
| 95 | # maximum spread in the beam direction in m
|
|---|
| 96 | set ZVertexSpread 0.10
|
|---|
| 97 |
|
|---|
| 98 | # maximum spread in time in s
|
|---|
| 99 | set TVertexSpread 1.5E-09
|
|---|
| 100 |
|
|---|
| 101 | # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
|
|---|
| 102 |
|
|---|
| 103 | set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
|
|---|
| 104 |
|
|---|
| 105 | #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
|
|---|
| 106 | # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) +
|
|---|
| 107 | # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) +
|
|---|
| 108 | # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
|
|---|
| 109 |
|
|---|
| 110 |
|
|---|
| 111 | }
|
|---|
| 112 |
|
|---|
| 113 |
|
|---|
| 114 | #################################
|
|---|
| 115 | # Propagate particles in cylinder
|
|---|
| 116 | #################################
|
|---|
| 117 |
|
|---|
| 118 | module ParticlePropagator ParticlePropagator {
|
|---|
| 119 | set InputArray Delphes/stableParticles
|
|---|
| 120 |
|
|---|
| 121 | set OutputArray stableParticles
|
|---|
| 122 | set ChargedHadronOutputArray chargedHadrons
|
|---|
| 123 | set ElectronOutputArray electrons
|
|---|
| 124 | set MuonOutputArray muons
|
|---|
| 125 |
|
|---|
| 126 | # radius of the magnetic field coverage, in m
|
|---|
| 127 | set Radius 6.00
|
|---|
| 128 | # half-length of the magnetic field coverage, in m
|
|---|
| 129 | set HalfLength 11.5
|
|---|
| 130 |
|
|---|
| 131 | # magnetic field
|
|---|
| 132 | set Bz 6.0
|
|---|
| 133 | }
|
|---|
| 134 |
|
|---|
| 135 | ####################################
|
|---|
| 136 | # Charged hadron tracking efficiency
|
|---|
| 137 | ####################################
|
|---|
| 138 |
|
|---|
| 139 | module Efficiency ChargedHadronTrackingEfficiency {
|
|---|
| 140 | set InputArray ParticlePropagator/chargedHadrons
|
|---|
| 141 | set OutputArray chargedHadrons
|
|---|
| 142 |
|
|---|
| 143 | # add EfficiencyFormula {efficiency formula as a function of eta and pt}
|
|---|
| 144 |
|
|---|
| 145 | # tracking efficiency formula for charged hadrons
|
|---|
| 146 | set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
|
|---|
| 147 | (abs(eta) <= 1.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
|
|---|
| 148 | (abs(eta) <= 1.5) * (pt > 1) * (0.95) + \
|
|---|
| 149 | (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
|
|---|
| 150 | (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
|
|---|
| 151 | (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
|
|---|
| 152 | (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
|
|---|
| 153 | (abs(eta) > 6.0) * (0.00)}
|
|---|
| 154 | }
|
|---|
| 155 |
|
|---|
| 156 | ##############################
|
|---|
| 157 | # Electron tracking efficiency
|
|---|
| 158 | ##############################
|
|---|
| 159 |
|
|---|
| 160 | module Efficiency ElectronEfficiency {
|
|---|
| 161 | set InputArray ParticlePropagator/electrons
|
|---|
| 162 | #set InputArray TrackPileUpSubtractor/electrons
|
|---|
| 163 | set OutputArray electrons
|
|---|
| 164 |
|
|---|
| 165 | # set EfficiencyFormula {efficiency formula as a function of eta and pt}
|
|---|
| 166 |
|
|---|
| 167 | # tracking efficiency formula for electrons
|
|---|
| 168 | # tracking efficiency formula for electrons
|
|---|
| 169 | set EfficiencyFormula { (pt <= 10) * (0.00) + \
|
|---|
| 170 | (abs(eta) <= 1.5) * (pt > 10 && pt <= 50) * (0.80) + \
|
|---|
| 171 | (abs(eta) <= 1.5) * (pt > 50) * (0.90) + \
|
|---|
| 172 | (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 10 && pt <= 50) * (0.80) + \
|
|---|
| 173 | (abs(eta) > 1.5 && abs(eta) <= 4) * (pt > 50) * (0.90) + \
|
|---|
| 174 | (abs(eta) > 4 && abs(eta) <= 6) * (pt > 10 && pt <= 50) * (0.70) + \
|
|---|
| 175 | (abs(eta) > 4 && abs(eta) <= 6) * (pt > 50) * (0.80) + \
|
|---|
| 176 | (abs(eta) > 6) * (0.00)}
|
|---|
| 177 | }
|
|---|
| 178 |
|
|---|
| 179 | ##########################
|
|---|
| 180 | # Muon tracking efficiency
|
|---|
| 181 | ##########################
|
|---|
| 182 |
|
|---|
| 183 | module Efficiency MuonEfficiency {
|
|---|
| 184 | set InputArray ParticlePropagator/muons
|
|---|
| 185 | # set InputArray TrackPileUpSubtractor/muons
|
|---|
| 186 | set OutputArray muons
|
|---|
| 187 |
|
|---|
| 188 | # set EfficiencyFormula {efficiency formula as a function of eta and pt}
|
|---|
| 189 |
|
|---|
| 190 | # tracking efficiency formula for muons
|
|---|
| 191 | set EfficiencyFormula { (pt <= 10) * (0.00) + \
|
|---|
| 192 | (abs(eta) <= 1.5) * (pt > 10 ) * (0.99) + \
|
|---|
| 193 | (abs(eta) <= 4 & abs(eta) > 1.5) * (pt > 10 ) * (0.99) + \
|
|---|
| 194 | (abs(eta) > 4) * (0.00)}
|
|---|
| 195 | }
|
|---|
| 196 |
|
|---|
| 197 | ########################################
|
|---|
| 198 | # Momentum resolution for charged tracks
|
|---|
| 199 | ########################################
|
|---|
| 200 |
|
|---|
| 201 | module MomentumSmearing ChargedHadronMomentumSmearing {
|
|---|
| 202 | set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
|
|---|
| 203 | set OutputArray chargedHadrons
|
|---|
| 204 |
|
|---|
| 205 | # set ResolutionFormula {resolution formula as a function of eta and pt}
|
|---|
| 206 |
|
|---|
| 207 | # resolution formula for charged hadrons
|
|---|
| 208 | #
|
|---|
| 209 | # Automatically generated tracker resolution formula for layout: FCChh_Option2 (by Z. Drasal, CERN)
|
|---|
| 210 | #
|
|---|
| 211 | # Central and Dipole
|
|---|
| 212 | #
|
|---|
| 213 | set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00194513) + \
|
|---|
| 214 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001945 + (pt-5.000000)* 0.000003) + \
|
|---|
| 215 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.001962 + (pt-10.000000)* 0.000006) + \
|
|---|
| 216 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.002547 + (pt-100.000000)* 0.000014) + \
|
|---|
| 217 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.014708 + (pt-1000.000000)* 0.000015) + \
|
|---|
| 218 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10000.0000) * (0.145553*pt/10000.000000) + \
|
|---|
| 219 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00247067) + \
|
|---|
| 220 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002471 + (pt-5.000000)* -0.000002) + \
|
|---|
| 221 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002462 + (pt-10.000000)* 0.000008) + \
|
|---|
| 222 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.003170 + (pt-100.000000)* 0.000015) + \
|
|---|
| 223 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.016485 + (pt-1000.000000)* 0.000016) + \
|
|---|
| 224 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10000.0000) * (0.162501*pt/10000.000000) + \
|
|---|
| 225 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00242524) + \
|
|---|
| 226 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002425 + (pt-5.000000)* 0.000013) + \
|
|---|
| 227 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002490 + (pt-10.000000)* 0.000042) + \
|
|---|
| 228 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.006273 + (pt-100.000000)* 0.000049) + \
|
|---|
| 229 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.050079 + (pt-1000.000000)* 0.000050) + \
|
|---|
| 230 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10000.0000) * (0.498773*pt/10000.000000) + \
|
|---|
| 231 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00168511) + \
|
|---|
| 232 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001685 + (pt-5.000000)* 0.000115) + \
|
|---|
| 233 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002262 + (pt-10.000000)* 0.000167) + \
|
|---|
| 234 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.017274 + (pt-100.000000)* 0.000172) + \
|
|---|
| 235 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.172095 + (pt-1000.000000)* 0.000172) + \
|
|---|
| 236 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10000.0000) * (1.720887*pt/10000.000000) + \
|
|---|
| 237 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00277124) + \
|
|---|
| 238 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002771 + (pt-5.000000)* 0.000430) + \
|
|---|
| 239 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.004921 + (pt-10.000000)* 0.000467) + \
|
|---|
| 240 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.046933 + (pt-100.000000)* 0.000469) + \
|
|---|
| 241 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.469096 + (pt-1000.000000)* 0.000469) + \
|
|---|
| 242 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10000.0000) * (4.690935*pt/10000.000000) + \
|
|---|
| 243 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00653887) + \
|
|---|
| 244 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.006539 + (pt-5.000000)* 0.001256) + \
|
|---|
| 245 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.012821 + (pt-10.000000)* 0.001272) + \
|
|---|
| 246 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.127345 + (pt-100.000000)* 0.001273) + \
|
|---|
| 247 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (1.273365 + (pt-1000.000000)* 0.001273) + \
|
|---|
| 248 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10000.0000) * (12.733637*pt/10000.000000)
|
|---|
| 249 | }
|
|---|
| 250 | }
|
|---|
| 251 | #################################
|
|---|
| 252 | # Energy resolution for electrons
|
|---|
| 253 | #################################
|
|---|
| 254 |
|
|---|
| 255 | module EnergySmearing ElectronEnergySmearing {
|
|---|
| 256 | set InputArray ElectronEfficiency/electrons
|
|---|
| 257 | set OutputArray electrons
|
|---|
| 258 |
|
|---|
| 259 | # set ResolutionFormula {resolution formula as a function of eta and energy}
|
|---|
| 260 |
|
|---|
| 261 | # check compared to compare to tracket resolution for low
|
|---|
| 262 | # resolution formula for electrons
|
|---|
| 263 | set ResolutionFormula {
|
|---|
| 264 | (abs(eta) <= 1.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
|
|---|
| 265 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
|
|---|
| 266 | (abs(eta) > 2.5 && abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
|
|---|
| 267 | (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
|
|---|
| 268 | }
|
|---|
| 269 |
|
|---|
| 270 |
|
|---|
| 271 | #################################
|
|---|
| 272 | # Energy resolution for photons
|
|---|
| 273 | #################################
|
|---|
| 274 | #not used for the moment (has been removed from the ExecutionPath)
|
|---|
| 275 |
|
|---|
| 276 | module EnergySmearing PhotonEnergySmearing {
|
|---|
| 277 | set InputArray PhotonEfficiency/photons
|
|---|
| 278 | set OutputArray photons
|
|---|
| 279 |
|
|---|
| 280 | # set ResolutionFormula {resolution formula as a function of eta and energy}
|
|---|
| 281 |
|
|---|
| 282 | # check compared to compare to track resolution for low
|
|---|
| 283 | # resolution formula for photons
|
|---|
| 284 | set ResolutionFormula {
|
|---|
| 285 | (abs(eta) <= 1.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
|
|---|
| 286 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
|
|---|
| 287 | (abs(eta) > 2.5 && abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
|
|---|
| 288 | (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
|
|---|
| 289 | }
|
|---|
| 290 |
|
|---|
| 291 |
|
|---|
| 292 | ###############################
|
|---|
| 293 | # Momentum resolution for muons
|
|---|
| 294 | ###############################
|
|---|
| 295 |
|
|---|
| 296 | module MomentumSmearing MuonMomentumSmearing {
|
|---|
| 297 | set InputArray MuonEfficiency/muons
|
|---|
| 298 | set OutputArray muons
|
|---|
| 299 |
|
|---|
| 300 | # set ResolutionFormula {resolution formula as a function of eta and pt}
|
|---|
| 301 |
|
|---|
| 302 | # resolution formula for muons
|
|---|
| 303 | # Automatically generated tracker resolution formula for layout: FCChh_Option2 (by Z. Drasal, CERN)
|
|---|
| 304 | # Central and Dipole
|
|---|
| 305 | #
|
|---|
| 306 | set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00194513) + \
|
|---|
| 307 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001945 + (pt-5.000000)* 0.000003) + \
|
|---|
| 308 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.001962 + (pt-10.000000)* 0.000006) + \
|
|---|
| 309 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.002547 + (pt-100.000000)* 0.000014) + \
|
|---|
| 310 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.014708 + (pt-1000.000000)* 0.000015) + \
|
|---|
| 311 | (abs(eta) >= 0.0000 && abs(eta) < 1.0000) * (pt >= 10000.0000) * (0.145553*pt/10000.000000) + \
|
|---|
| 312 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00247067) + \
|
|---|
| 313 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002471 + (pt-5.000000)* -0.000002) + \
|
|---|
| 314 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002462 + (pt-10.000000)* 0.000008) + \
|
|---|
| 315 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.003170 + (pt-100.000000)* 0.000015) + \
|
|---|
| 316 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.016485 + (pt-1000.000000)* 0.000016) + \
|
|---|
| 317 | (abs(eta) >= 1.0000 && abs(eta) < 2.0000) * (pt >= 10000.0000) * (0.162501*pt/10000.000000) + \
|
|---|
| 318 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00242524) + \
|
|---|
| 319 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002425 + (pt-5.000000)* 0.000013) + \
|
|---|
| 320 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002490 + (pt-10.000000)* 0.000042) + \
|
|---|
| 321 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.006273 + (pt-100.000000)* 0.000049) + \
|
|---|
| 322 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.050079 + (pt-1000.000000)* 0.000050) + \
|
|---|
| 323 | (abs(eta) >= 2.0000 && abs(eta) < 3.0000) * (pt >= 10000.0000) * (0.498773*pt/10000.000000) + \
|
|---|
| 324 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00168511) + \
|
|---|
| 325 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.001685 + (pt-5.000000)* 0.000115) + \
|
|---|
| 326 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.002262 + (pt-10.000000)* 0.000167) + \
|
|---|
| 327 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.017274 + (pt-100.000000)* 0.000172) + \
|
|---|
| 328 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.172095 + (pt-1000.000000)* 0.000172) + \
|
|---|
| 329 | (abs(eta) >= 3.0000 && abs(eta) < 4.0000) * (pt >= 10000.0000) * (1.720887*pt/10000.000000) + \
|
|---|
| 330 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00277124) + \
|
|---|
| 331 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.002771 + (pt-5.000000)* 0.000430) + \
|
|---|
| 332 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.004921 + (pt-10.000000)* 0.000467) + \
|
|---|
| 333 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.046933 + (pt-100.000000)* 0.000469) + \
|
|---|
| 334 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (0.469096 + (pt-1000.000000)* 0.000469) + \
|
|---|
| 335 | (abs(eta) >= 4.0000 && abs(eta) < 5.0000) * (pt >= 10000.0000) * (4.690935*pt/10000.000000) + \
|
|---|
| 336 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 0.0000 && pt < 5.0000) * (0.00653887) + \
|
|---|
| 337 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 5.0000 && pt < 10.0000) * (0.006539 + (pt-5.000000)* 0.001256) + \
|
|---|
| 338 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.012821 + (pt-10.000000)* 0.001272) + \
|
|---|
| 339 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 100.0000 && pt < 1000.0000) * (0.127345 + (pt-100.000000)* 0.001273) + \
|
|---|
| 340 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 1000.0000 && pt < 10000.0000) * (1.273365 + (pt-1000.000000)* 0.001273) + \
|
|---|
| 341 | (abs(eta) >= 5.0000 && abs(eta) < 6.0000) * (pt >= 10000.0000) * (12.733637*pt/10000.000000)
|
|---|
| 342 | }
|
|---|
| 343 | }
|
|---|
| 344 | ####################
|
|---|
| 345 | # Electron isolation
|
|---|
| 346 | ####################
|
|---|
| 347 |
|
|---|
| 348 | module Isolation ElectronIsolation {
|
|---|
| 349 | set CandidateInputArray ElectronEnergySmearing/electrons
|
|---|
| 350 | set IsolationInputArray EFlowMerger/eflow
|
|---|
| 351 |
|
|---|
| 352 | set OutputArray electrons
|
|---|
| 353 |
|
|---|
| 354 | set DeltaRMax 0.3
|
|---|
| 355 |
|
|---|
| 356 | set PTMin 0.5
|
|---|
| 357 |
|
|---|
| 358 | set PTRatioMax 0.1
|
|---|
| 359 | }
|
|---|
| 360 |
|
|---|
| 361 | ################
|
|---|
| 362 | # Muon isolation
|
|---|
| 363 | ################
|
|---|
| 364 |
|
|---|
| 365 | module Isolation MuonIsolation {
|
|---|
| 366 | set CandidateInputArray MuonMomentumSmearing/muons
|
|---|
| 367 | set IsolationInputArray EFlowMerger/eflow
|
|---|
| 368 |
|
|---|
| 369 | set OutputArray muons
|
|---|
| 370 |
|
|---|
| 371 | set DeltaRMax 0.3
|
|---|
| 372 |
|
|---|
| 373 | set PTMin 0.5
|
|---|
| 374 |
|
|---|
| 375 | set PTRatioMax 0.1
|
|---|
| 376 | }
|
|---|
| 377 |
|
|---|
| 378 | ###################
|
|---|
| 379 | # Photon efficiency
|
|---|
| 380 | ###################
|
|---|
| 381 |
|
|---|
| 382 | module Efficiency PhotonEfficiency {
|
|---|
| 383 | set InputArray Ecal/eflowPhotons
|
|---|
| 384 | set OutputArray photons
|
|---|
| 385 |
|
|---|
| 386 | # set EfficiencyFormula {efficiency formula as a function of eta and pt}
|
|---|
| 387 | # efficiency formula for photons
|
|---|
| 388 | set EfficiencyFormula {
|
|---|
| 389 | (pt <= 10.0) * (0.00) + \
|
|---|
| 390 | (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
|
|---|
| 391 | (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) + \
|
|---|
| 392 | (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
|
|---|
| 393 | (abs(eta) > 6.0) * (0.00)}
|
|---|
| 394 | }
|
|---|
| 395 |
|
|---|
| 396 | ##################
|
|---|
| 397 | # Photon isolation
|
|---|
| 398 | ##################
|
|---|
| 399 |
|
|---|
| 400 | module Isolation PhotonIsolation {
|
|---|
| 401 | set CandidateInputArray PhotonEfficiency/photons
|
|---|
| 402 | set IsolationInputArray EFlowMerger/eflow
|
|---|
| 403 |
|
|---|
| 404 | set OutputArray photons
|
|---|
| 405 |
|
|---|
| 406 | set DeltaRMax 0.3
|
|---|
| 407 |
|
|---|
| 408 | set PTMin 0.5
|
|---|
| 409 |
|
|---|
| 410 | set PTRatioMax 0.1
|
|---|
| 411 | }
|
|---|
| 412 |
|
|---|
| 413 | ##############
|
|---|
| 414 | # Track merger
|
|---|
| 415 | ##############
|
|---|
| 416 |
|
|---|
| 417 | module Merger TrackMerger {
|
|---|
| 418 | # add InputArray InputArray
|
|---|
| 419 | add InputArray ChargedHadronMomentumSmearing/chargedHadrons
|
|---|
| 420 | add InputArray ElectronEnergySmearing/electrons
|
|---|
| 421 | add InputArray MuonMomentumSmearing/muons
|
|---|
| 422 | set OutputArray tracks
|
|---|
| 423 | }
|
|---|
| 424 |
|
|---|
| 425 | ################################
|
|---|
| 426 | # Track impact parameter smearing
|
|---|
| 427 | ################################
|
|---|
| 428 |
|
|---|
| 429 | #not used at the moment
|
|---|
| 430 |
|
|---|
| 431 | module ImpactParameterSmearing ImpactParameterSmearing {
|
|---|
| 432 | set InputArray TrackMerger/tracks
|
|---|
| 433 | set OutputArray tracks
|
|---|
| 434 |
|
|---|
| 435 |
|
|---|
| 436 | # absolute impact parameter smearing formula (in mm) as a function of pt and eta
|
|---|
| 437 | set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) + \
|
|---|
| 438 | (pt > 5.0) * (0.005)}
|
|---|
| 439 |
|
|---|
| 440 | }
|
|---|
| 441 |
|
|---|
| 442 | #############
|
|---|
| 443 | # ECAL
|
|---|
| 444 | #############
|
|---|
| 445 |
|
|---|
| 446 | module SimpleCalorimeter Ecal {
|
|---|
| 447 | set ParticleInputArray ParticlePropagator/stableParticles
|
|---|
| 448 | set TrackInputArray ImpactParameterSmearing/tracks
|
|---|
| 449 |
|
|---|
| 450 | set TowerOutputArray ecalTowers
|
|---|
| 451 | set EFlowTowerOutputArray eflowPhotons
|
|---|
| 452 |
|
|---|
| 453 | set pi [expr {acos(-1)}]
|
|---|
| 454 |
|
|---|
| 455 | # lists of the edges of each tower in eta and phi
|
|---|
| 456 | # each list starts with the lower edge of the first tower
|
|---|
| 457 | # the list ends with the higher edged of the last tower
|
|---|
| 458 | # 0.0125 x 0.0125 for eta between 0 and 2.5
|
|---|
| 459 | set PhiBins {}
|
|---|
| 460 | for {set i -256} {$i <= 256} {incr i} {
|
|---|
| 461 | add PhiBins [expr {$i * $pi/256.0}]
|
|---|
| 462 | }
|
|---|
| 463 | foreach eta {-2.5 -2.4875 -2.475 -2.4625 -2.45 -2.4375 -2.425 -2.4125 -2.4 -2.3875 -2.375 -2.3625 -2.35 -2.3375 -2.325 -2.3125 -2.3 -2.2875 -2.275 -2.2625 -2.25 -2.2375 -2.225 -2.2125 -2.2 -2.1875 -2.175 -2.1625 -2.15 -2.1375 -2.125 -2.1125 -2.1 -2.0875 -2.075 -2.0625 -2.05 -2.0375 -2.025 -2.0125 -2 -1.9875 -1.975 -1.9625 -1.95 -1.9375 -1.925 -1.9125 -1.9 -1.8875 -1.875 -1.8625 -1.85 -1.8375 -1.825 -1.8125 -1.8 -1.7875 -1.775 -1.7625 -1.75 -1.7375 -1.725 -1.7125 -1.7 -1.6875 -1.675 -1.6625 -1.65 -1.6375 -1.625 -1.6125 -1.6 -1.5875 -1.575 -1.5625 -1.55 -1.5375 -1.525 -1.5125 -1.5 -1.4875 -1.475 -1.4625 -1.45 -1.4375 -1.425 -1.4125 -1.4 -1.3875 -1.375 -1.3625 -1.35 -1.3375 -1.325 -1.3125 -1.3 -1.2875 -1.275 -1.2625 -1.25 -1.2375 -1.225 -1.2125 -1.2 -1.1875 -1.175 -1.1625 -1.15 -1.1375 -1.125 -1.1125 -1.1 -1.0875 -1.075 -1.0625 -1.05 -1.0375 -1.025 -1.0125 -1.0 -0.9875 -0.975 -0.9625 -0.95 -0.9375 -0.925 -0.9125 -0.9 -0.8875 -0.875 -0.8625 -0.85 -0.8375 -0.825 -0.8125 -0.8 -0.7875 -0.775 -0.7625 -0.75 -0.7375 -0.725 -0.7125 -0.7 -0.6875 -0.675 -0.6625 -0.65 -0.6375 -0.625 -0.6125 -0.6 -0.5875 -0.575 -0.5625 -0.55 -0.5375 -0.525 -0.5125 -0.5 -0.4875 -0.475 -0.4625 -0.45 -0.4375 -0.425 -0.4125 -0.4 -0.3875 -0.375 -0.3625 -0.35 -0.3375 -0.325 -0.3125 -0.3 -0.2875 -0.275 -0.2625 -0.25 -0.2375 -0.225 -0.2125 -0.2 -0.1875 -0.175 -0.1625 -0.15 -0.1375 -0.125 -0.1125 -0.1 -0.0875 -0.075 -0.0625 -0.05 -0.0375 -0.025 -0.0125 0 0.0125 0.025 0.0375 0.05 0.0625 0.075 0.0875 0.1 0.1125 0.125 0.1375 0.15 0.1625 0.175 0.1875 0.2 0.2125 0.225 0.2375 0.25 0.2625 0.275 0.2875 0.3 0.3125 0.325 0.3375 0.35 0.3625 0.375 0.3875 0.4 0.4125 0.425 0.4375 0.45 0.4625 0.475 0.4875 0.5 0.5125 0.525 0.5375 0.55 0.5625 0.575 0.5875 0.6 0.6125 0.625 0.6375 0.65 0.6625 0.675 0.6875 0.7 0.7125 0.725 0.7375 0.75 0.7625 0.775 0.7875 0.8 0.8125 0.825 0.8375 0.85 0.8625 0.875 0.8875 0.9 0.9125 0.925 0.9375 0.95 0.9625 0.975 0.9875 1.0 1.0125 1.025 1.0375 1.05 1.0625 1.075 1.0875 1.1 1.1125 1.125 1.1375 1.15 1.1625 1.175 1.1875 1.2 1.2125 1.225 1.2375 1.25 1.2625 1.275 1.2875 1.3 1.3125 1.325 1.3375 1.35 1.3625 1.375 1.3875 1.4 1.4125 1.425 1.4375 1.45 1.4625 1.475 1.4875 1.5 1.5125 1.525 1.5375 1.55 1.5625 1.575 1.5875 1.6 1.6125 1.625 1.6375 1.65 1.6625 1.675 1.6875 1.7 1.7125 1.725 1.7375 1.75 1.7625 1.775 1.7875 1.8 1.8125 1.825 1.8375 1.85 1.8625 1.875 1.8875 1.9 1.9125 1.925 1.9375 1.95 1.9625 1.975 1.9875 2 2.0125 2.025 2.0375 2.05 2.0625 2.075 2.0875 2.1 2.1125 2.125 2.1375 2.15 2.1625 2.175 2.1875 2.2 2.2125 2.225 2.2375 2.25 2.2625 2.275 2.2875 2.3 2.3125 2.325 2.3375 2.35 2.3625 2.375 2.3875 2.4 2.4125 2.425 2.4375 2.45 2.4625 2.475 2.4875 2.5} {
|
|---|
| 464 | add EtaPhiBins $eta $PhiBins
|
|---|
| 465 | }
|
|---|
| 466 | # 0.025 eta x 0.025 for eta between 2.5 and 4.0
|
|---|
| 467 | set PhiBins {}
|
|---|
| 468 | for {set i -128} {$i <= 128} {incr i} {
|
|---|
| 469 | add PhiBins [expr {$i * $pi/128.0}]
|
|---|
| 470 | }
|
|---|
| 471 | foreach eta {-4.0 -3.975 -3.95 -3.925 -3.9 -3.875 -3.85 -3.825 -3.8 -3.775 -3.75 -3.725 -3.7 -3.675 -3.65 -3.625 -3.6 -3.575 -3.55 -3.525 -3.5 -3.475 -3.45 -3.425 -3.4 -3.375 -3.35 -3.325 -3.3 -3.275 -3.25 -3.225 -3.2 -3.175 -3.15 -3.125 -3.1 -3.075 -3.05 -3.025 -3.0 -2.975 -2.95 -2.925 -2.9 -2.875 -2.85 -2.825 -2.8 -2.775 -2.75 -2.725 -2.7 -2.675 -2.65 -2.625 -2.6 -2.575 -2.55 -2.525 2.525 2.55 2.575 2.6 2.625 2.65 2.675 2.7 2.725 2.75 2.775 2.8 2.825 2.85 2.875 2.9 2.925 2.95 2.975 3.0 3.025 3.05 3.075 3.1 3.125 3.15 3.175 3.2 3.225 3.25 3.275 3.3 3.325 3.35 3.375 3.4 3.425 3.45 3.475 3.5 3.525 3.55 3.575 3.6 3.625 3.65 3.675 3.7 3.725 3.75 3.775 3.8 3.825 3.85 3.875 3.9 3.925 3.95 3.975 4.0} {
|
|---|
| 472 | add EtaPhiBins $eta $PhiBins
|
|---|
| 473 | }
|
|---|
| 474 |
|
|---|
| 475 | # 0.05 x 0.05 for eta between 4.0 and 6.0
|
|---|
| 476 | set PhiBins {}
|
|---|
| 477 | for {set i -64} {$i <= 64} {incr i} {
|
|---|
| 478 | add PhiBins [expr {$i * $pi/64.0}]
|
|---|
| 479 | }
|
|---|
| 480 | foreach eta {-6.0 -5.95 -5.9 -5.85 -5.8 -5.75 -5.7 -5.65 -5.6 -5.55 -5.5 -5.45 -5.4 -5.35 -5.3 -5.25 -5.2 -5.15 -5.1 -5.05 -5.0 -4.95 -4.9 -4.85 -4.8 -4.75 -4.7 -4.65 -4.6 -4.55 -4.5 -4.45 -4.4 -4.35 -4.3 -4.25 -4.2 -4.15 -4.1 -4.05 4.05 4.1 4.15 4.2 4.25 4.3 4.35 4.4 4.45 4.5 4.55 4.6 4.65 4.7 4.75 4.8 4.85 4.9 4.95 5.0 5.05 5.1 5.15 5.2 5.25 5.3 5.35 5.4 5.45 5.5 5.55 5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6.0} {
|
|---|
| 481 | add EtaPhiBins $eta $PhiBins
|
|---|
| 482 | }
|
|---|
| 483 | # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
|
|---|
| 484 |
|
|---|
| 485 | add EnergyFraction {0} {0.4}
|
|---|
| 486 | # energy fractions for e, gamma and pi0
|
|---|
| 487 | add EnergyFraction {11} {1.0}
|
|---|
| 488 | add EnergyFraction {22} {1.0}
|
|---|
| 489 | add EnergyFraction {111} {1.0}
|
|---|
| 490 | # energy fractions for muon, neutrinos and neutralinos
|
|---|
| 491 | add EnergyFraction {12} {0.0}
|
|---|
| 492 | add EnergyFraction {13} {0.0}
|
|---|
| 493 | add EnergyFraction {14} {0.0}
|
|---|
| 494 | add EnergyFraction {16} {0.0}
|
|---|
| 495 | add EnergyFraction {1000022} {0.0}
|
|---|
| 496 | add EnergyFraction {1000023} {0.0}
|
|---|
| 497 | add EnergyFraction {1000025} {0.0}
|
|---|
| 498 | add EnergyFraction {1000035} {0.0}
|
|---|
| 499 | add EnergyFraction {1000045} {0.0}
|
|---|
| 500 | # energy fractions for K0short and Lambda
|
|---|
| 501 | # add EnergyFraction {310} {0.3}
|
|---|
| 502 | # add EnergyFraction {3122} {0.3}
|
|---|
| 503 |
|
|---|
| 504 | # set ECalResolutionFormula {resolution formula as a function of eta and energy}
|
|---|
| 505 | set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
|
|---|
| 506 | (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
|
|---|
| 507 |
|
|---|
| 508 |
|
|---|
| 509 | }
|
|---|
| 510 | ####################
|
|---|
| 511 | # Neutral tower merger
|
|---|
| 512 | ####################
|
|---|
| 513 |
|
|---|
| 514 | module Merger NeutralTowerMerger {
|
|---|
| 515 | # add InputArray InputArray
|
|---|
| 516 | add InputArray Ecal/eflowPhotons
|
|---|
| 517 | add InputArray Hcal/eflowNeutralHadrons
|
|---|
| 518 | set OutputArray eflowTowers
|
|---|
| 519 | }
|
|---|
| 520 |
|
|---|
| 521 |
|
|---|
| 522 | #############
|
|---|
| 523 | # HCAL
|
|---|
| 524 | #############
|
|---|
| 525 |
|
|---|
| 526 | module SimpleCalorimeter Hcal {
|
|---|
| 527 | set ParticleInputArray ParticlePropagator/stableParticles
|
|---|
| 528 | set TrackInputArray ImpactParameterSmearing/tracks
|
|---|
| 529 |
|
|---|
| 530 | set TowerOutputArray hcalTowers
|
|---|
| 531 | set EFlowTowerOutputArray eflowNeutralHadrons
|
|---|
| 532 |
|
|---|
| 533 | set pi [expr {acos(-1)}]
|
|---|
| 534 |
|
|---|
| 535 | # lists of the edges of each tower in eta and phi
|
|---|
| 536 | # each list starts with the lower edge of the first tower
|
|---|
| 537 | # the list ends with the higher edged of the last tower
|
|---|
| 538 |
|
|---|
| 539 | # 0.05 x 0.05 between eta 0 and 2.5
|
|---|
| 540 | set PhiBins {}
|
|---|
| 541 | for {set i -128} {$i <= 128} {incr i} {
|
|---|
| 542 | add PhiBins [expr {$i * $pi/128.0}]
|
|---|
| 543 | }
|
|---|
| 544 | foreach eta {-2.5 -2.45 -2.4 -2.35 -2.3 -2.25 -2.2 -2.15 -2.1 -2.05 -2.0 -1.95 -1.9 -1.85 -1.8 -1.75 -1.7 -1.65 -1.6 -1.55 -1.5 -1.45 -1.4 -1.35 -1.3 -1.25 -1.2 -1.15 -1.1 -1.05 -1.0 -0.95 -0.9 -0.85 -0.8 -0.75 -0.7 -0.65 -0.6 -0.55 -0.5 -0.45 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1.0 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2.0 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5} {
|
|---|
| 545 | add EtaPhiBins $eta $PhiBins
|
|---|
| 546 | }
|
|---|
| 547 |
|
|---|
| 548 | # 0.1 x 0.1 between eta 2.5 and 4.0
|
|---|
| 549 | set PhiBins {}
|
|---|
| 550 | for {set i -64} {$i <= 64} {incr i} {
|
|---|
| 551 | add PhiBins [expr {$i * $pi/64.0}]
|
|---|
| 552 | }
|
|---|
| 553 | foreach eta {-4.0 -3.9 -3.8 -3.7 -3.6 -3.5 -3.4 -3.3 -3.2 -3.1 -3.0 -2.9 -2.8 -2.7 -2.6 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0} {
|
|---|
| 554 | add EtaPhiBins $eta $PhiBins
|
|---|
| 555 | }
|
|---|
| 556 |
|
|---|
| 557 | # 0.2 x 0.2 between eta 4.0 and 6.0
|
|---|
| 558 | set PhiBins {}
|
|---|
| 559 | for {set i -32} {$i <= 32} {incr i} {
|
|---|
| 560 | add PhiBins [expr {$i * $pi/32.0}]
|
|---|
| 561 | }
|
|---|
| 562 | foreach eta {-6.0 -5.8 -5.6 -5.4 -5.2 -5.0 -4.8 -4.6 -4.4 -4.2 4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0} {
|
|---|
| 563 | add EtaPhiBins $eta $PhiBins
|
|---|
| 564 | }
|
|---|
| 565 |
|
|---|
| 566 | # default energy fractions {abs(PDG code)} {Fecal Fhcal}
|
|---|
| 567 | add EnergyFraction {0} {0.6}
|
|---|
| 568 | # energy fractions for e, gamma and pi0
|
|---|
| 569 | add EnergyFraction {11} {0.0}
|
|---|
| 570 | add EnergyFraction {22} {0.0}
|
|---|
| 571 | add EnergyFraction {111} {0.0}
|
|---|
| 572 | # energy fractions for muon, neutrinos and neutralinos
|
|---|
| 573 | add EnergyFraction {12} {0.0}
|
|---|
| 574 | add EnergyFraction {13} {0.0}
|
|---|
| 575 | add EnergyFraction {14} {0.0}
|
|---|
| 576 | add EnergyFraction {16} {0.0}
|
|---|
| 577 | add EnergyFraction {1000022} {0.0}
|
|---|
| 578 | add EnergyFraction {1000023} {0.0}
|
|---|
| 579 | add EnergyFraction {1000025} {0.0}
|
|---|
| 580 | add EnergyFraction {1000035} {0.0}
|
|---|
| 581 | add EnergyFraction {1000045} {0.0}
|
|---|
| 582 | # energy fractions for K0short and Lambda
|
|---|
| 583 | # add EnergyFraction {310} {0.7}
|
|---|
| 584 | # add EnergyFraction {3122} {0.7}
|
|---|
| 585 |
|
|---|
| 586 | # set HCalResolutionFormula {resolution formula as a function of eta and energy}
|
|---|
| 587 | set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
|
|---|
| 588 | (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
|
|---|
| 589 | }
|
|---|
| 590 |
|
|---|
| 591 |
|
|---|
| 592 | ##########################
|
|---|
| 593 | # Track pile-up subtractor
|
|---|
| 594 | ##########################
|
|---|
| 595 |
|
|---|
| 596 | module TrackPileUpSubtractor TrackPileUpSubtractor {
|
|---|
| 597 | # add InputArray InputArray OutputArray
|
|---|
| 598 | # add InputArray Calorimeter/eflowTracks eflowTracks
|
|---|
| 599 | add InputArray ImpactParameterSmearing/tracks tracks
|
|---|
| 600 | # add InputArray ElectronFilter/electrons electrons
|
|---|
| 601 | add InputArray ElectronEfficiency/electrons electrons
|
|---|
| 602 | add InputArray MuonMomentumSmearing/muons muons
|
|---|
| 603 |
|
|---|
| 604 | set VertexInputArray PileUpMerger/vertices
|
|---|
| 605 | # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
|
|---|
| 606 | # Z vertex resolution in m
|
|---|
| 607 | set ZVertexResolution 0.0001
|
|---|
| 608 | set OutputArray tracksPV
|
|---|
| 609 | }
|
|---|
| 610 |
|
|---|
| 611 |
|
|---|
| 612 | ####################
|
|---|
| 613 | # Tower Merger (in case not using e-flow algorithm)
|
|---|
| 614 | ####################
|
|---|
| 615 |
|
|---|
| 616 | module Merger TowerMerger {
|
|---|
| 617 | # add InputArray InputArray
|
|---|
| 618 | add InputArray Ecal/ecalTowers
|
|---|
| 619 | add InputArray Hcal/hcalTowers
|
|---|
| 620 | set OutputArray towers
|
|---|
| 621 | }
|
|---|
| 622 | ##################################
|
|---|
| 623 | # Energy flow merger (all tracks)
|
|---|
| 624 | ##################################
|
|---|
| 625 |
|
|---|
| 626 | module Merger EFlowMergerAllTracks {
|
|---|
| 627 | # add InputArray InputArray
|
|---|
| 628 | add InputArray ImpactParameterSmearing/tracks
|
|---|
| 629 | add InputArray Ecal/eflowPhotons
|
|---|
| 630 | add InputArray Hcal/eflowNeutralHadrons
|
|---|
| 631 | set OutputArray eflow
|
|---|
| 632 | }
|
|---|
| 633 | ####################
|
|---|
| 634 | # Energy flow merger
|
|---|
| 635 | ####################
|
|---|
| 636 |
|
|---|
| 637 | module Merger EFlowMerger {
|
|---|
| 638 | # add InputArray InputArray
|
|---|
| 639 | #add InputArray ImpactParameterSmearing/tracks
|
|---|
| 640 | add InputArray TrackPileUpSubtractor/tracks
|
|---|
| 641 | add InputArray Ecal/eflowPhotons
|
|---|
| 642 | add InputArray Hcal/eflowNeutralHadrons
|
|---|
| 643 | set OutputArray eflow
|
|---|
| 644 | }
|
|---|
| 645 | #############
|
|---|
| 646 | # Rho pile-up
|
|---|
| 647 | #############
|
|---|
| 648 |
|
|---|
| 649 | module FastJetGridMedianEstimator Rho {
|
|---|
| 650 |
|
|---|
| 651 | set InputArray EFlowMerger/eflow
|
|---|
| 652 | set RhoOutputArray rho
|
|---|
| 653 |
|
|---|
| 654 | # add GridRange rapmin rapmax drap dphi
|
|---|
| 655 | # rapmin - the minimum rapidity extent of the grid
|
|---|
| 656 | # rapmax - the maximum rapidity extent of the grid
|
|---|
| 657 | # drap - the grid spacing in rapidity
|
|---|
| 658 | # dphi - the grid spacing in azimuth
|
|---|
| 659 |
|
|---|
| 660 | add GridRange -5.0 -2.5 1.0 1.0
|
|---|
| 661 | add GridRange -2.5 2.5 1.0 1.0
|
|---|
| 662 | add GridRange 2.5 5.0 1.0 1.0
|
|---|
| 663 |
|
|---|
| 664 | }
|
|---|
| 665 |
|
|---|
| 666 | #####################
|
|---|
| 667 | # Neutrino Filter
|
|---|
| 668 | #####################
|
|---|
| 669 |
|
|---|
| 670 | module PdgCodeFilter NeutrinoFilter {
|
|---|
| 671 |
|
|---|
| 672 | set InputArray Delphes/stableParticles
|
|---|
| 673 | set OutputArray filteredParticles
|
|---|
| 674 |
|
|---|
| 675 | set PTMin 0.0
|
|---|
| 676 |
|
|---|
| 677 | add PdgCode {12}
|
|---|
| 678 | add PdgCode {14}
|
|---|
| 679 | add PdgCode {16}
|
|---|
| 680 | add PdgCode {-12}
|
|---|
| 681 | add PdgCode {-14}
|
|---|
| 682 | add PdgCode {-16}
|
|---|
| 683 |
|
|---|
| 684 | }
|
|---|
| 685 |
|
|---|
| 686 |
|
|---|
| 687 |
|
|---|
| 688 | ###################
|
|---|
| 689 | # Missing ET merger
|
|---|
| 690 | ###################
|
|---|
| 691 |
|
|---|
| 692 | module Merger MissingET {
|
|---|
| 693 | # add InputArray InputArray
|
|---|
| 694 | #default
|
|---|
| 695 | #add InputArray EFlowMerger/eflow
|
|---|
| 696 | add InputArray JetPileUpSubtractor/jets
|
|---|
| 697 |
|
|---|
| 698 | set MomentumOutputArray momentum
|
|---|
| 699 | }
|
|---|
| 700 |
|
|---|
| 701 |
|
|---|
| 702 | ##################
|
|---|
| 703 | # Scalar HT merger
|
|---|
| 704 | ##################
|
|---|
| 705 |
|
|---|
| 706 | module Merger ScalarHT {
|
|---|
| 707 | # add InputArray InputArray
|
|---|
| 708 | add InputArray EFlowMerger/eflow
|
|---|
| 709 | set EnergyOutputArray energy
|
|---|
| 710 | }
|
|---|
| 711 | ########################
|
|---|
| 712 | # Jet Flavor Association
|
|---|
| 713 | ########################
|
|---|
| 714 |
|
|---|
| 715 | module JetFlavorAssociation JetFlavorAssociation {
|
|---|
| 716 |
|
|---|
| 717 | set PartonInputArray Delphes/partons
|
|---|
| 718 | set ParticleInputArray Delphes/allParticles
|
|---|
| 719 | set ParticleLHEFInputArray Delphes/allParticlesLHEF
|
|---|
| 720 | set JetInputArray JetEnergyScale/jets
|
|---|
| 721 |
|
|---|
| 722 | set DeltaR 0.5
|
|---|
| 723 | set PartonPTMin 1.0
|
|---|
| 724 | set PartonEtaMax 6.0
|
|---|
| 725 |
|
|---|
| 726 | }
|
|---|
| 727 | #####################
|
|---|
| 728 | # MC truth jet finder
|
|---|
| 729 | #####################
|
|---|
| 730 |
|
|---|
| 731 | module FastJetFinder GenJetFinder {
|
|---|
| 732 | set InputArray Delphes/stableParticles
|
|---|
| 733 |
|
|---|
| 734 | set OutputArray jets
|
|---|
| 735 |
|
|---|
| 736 | # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 7 antikt wta algorithm
|
|---|
| 737 | set JetAlgorithm 7
|
|---|
| 738 | set ParameterR 0.4
|
|---|
| 739 |
|
|---|
| 740 | set JetPTMin 5.0
|
|---|
| 741 | }
|
|---|
| 742 |
|
|---|
| 743 | #########################
|
|---|
| 744 | # Gen Missing ET merger
|
|---|
| 745 | ########################
|
|---|
| 746 |
|
|---|
| 747 | module Merger GenMissingET {
|
|---|
| 748 | # add InputArray InputArray
|
|---|
| 749 | add InputArray NeutrinoFilter/filteredParticles
|
|---|
| 750 | set MomentumOutputArray momentum
|
|---|
| 751 | }
|
|---|
| 752 |
|
|---|
| 753 | ############
|
|---|
| 754 | # Jet finder
|
|---|
| 755 | ############
|
|---|
| 756 |
|
|---|
| 757 | module FastJetFinder FastJetFinder {
|
|---|
| 758 | # set InputArray Calorimeter/towers
|
|---|
| 759 | set InputArray EFlowMerger/eflow
|
|---|
| 760 |
|
|---|
| 761 | set OutputArray jets
|
|---|
| 762 |
|
|---|
| 763 | # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 antikt wta algorithm
|
|---|
| 764 | set JetAlgorithm 7
|
|---|
| 765 | set ParameterR 0.4
|
|---|
| 766 |
|
|---|
| 767 | set JetPTMin 30.0
|
|---|
| 768 | }
|
|---|
| 769 |
|
|---|
| 770 | ##################
|
|---|
| 771 | # Jet Energy Scale
|
|---|
| 772 | ##################
|
|---|
| 773 | # not used at the moment
|
|---|
| 774 |
|
|---|
| 775 | module EnergyScale JetEnergyScale {
|
|---|
| 776 | set InputArray FastJetFinder/jets
|
|---|
| 777 | set OutputArray jets
|
|---|
| 778 |
|
|---|
| 779 | # scale formula for jets
|
|---|
| 780 | set ScaleFormula {1.00}
|
|---|
| 781 | }
|
|---|
| 782 |
|
|---|
| 783 | ###########################
|
|---|
| 784 | # Jet Pile-Up ID
|
|---|
| 785 | ###########################
|
|---|
| 786 |
|
|---|
| 787 | module PileUpJetID PileUpJetID {
|
|---|
| 788 | set JetInputArray FastJetFinder/jets
|
|---|
| 789 | # set TrackInputArray Calorimeter/eflowTracks
|
|---|
| 790 | set TrackInputArray ImpactParameterSmearing/tracks
|
|---|
| 791 | set NeutralInputArray NeutralTowerMerger/eflowTowers
|
|---|
| 792 |
|
|---|
| 793 | set VertexInputArray PileUpMerger/vertices
|
|---|
| 794 | # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
|
|---|
| 795 | # Z vertex resolution in m
|
|---|
| 796 | set ZVertexResolution 0.0001
|
|---|
| 797 |
|
|---|
| 798 | set OutputArray jets
|
|---|
| 799 |
|
|---|
| 800 | set UseConstituents 0
|
|---|
| 801 | set ParameterR 0.5
|
|---|
| 802 |
|
|---|
| 803 | set JetPTMin 20.0
|
|---|
| 804 | }
|
|---|
| 805 |
|
|---|
| 806 | ###########################
|
|---|
| 807 | # Jet Pile-Up Subtraction
|
|---|
| 808 | ###########################
|
|---|
| 809 |
|
|---|
| 810 | module JetPileUpSubtractor JetPileUpSubtractor {
|
|---|
| 811 | set JetInputArray PileUpJetID/jets
|
|---|
| 812 | set RhoInputArray Rho/rho
|
|---|
| 813 |
|
|---|
| 814 | set OutputArray jets
|
|---|
| 815 |
|
|---|
| 816 | set JetPTMin 20.0
|
|---|
| 817 | }
|
|---|
| 818 |
|
|---|
| 819 |
|
|---|
| 820 |
|
|---|
| 821 | ##########################
|
|---|
| 822 | # tau-tagging
|
|---|
| 823 | ##########################
|
|---|
| 824 |
|
|---|
| 825 |
|
|---|
| 826 | module TauTagging TauTagging {
|
|---|
| 827 | set ParticleInputArray Delphes/allParticles
|
|---|
| 828 | set PartonInputArray Delphes/partons
|
|---|
| 829 | set JetInputArray JetEnergyScale/jets
|
|---|
| 830 | set DeltaR 0.4
|
|---|
| 831 | set TauPTMin 15.0
|
|---|
| 832 | set TauEtaMax 6
|
|---|
| 833 | # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
|
|---|
| 834 | add EfficiencyFormula {0} {0.01}
|
|---|
| 835 | add EfficiencyFormula {11} {0.005}
|
|---|
| 836 | add EfficiencyFormula {15} {0.6}
|
|---|
| 837 | }
|
|---|
| 838 |
|
|---|
| 839 |
|
|---|
| 840 | ##########################
|
|---|
| 841 | # b-tagging
|
|---|
| 842 | ##########################
|
|---|
| 843 |
|
|---|
| 844 | module BTagging BTagging { # need to check that this module exists or write it
|
|---|
| 845 | set JetInputArray JetEnergyScale/jets
|
|---|
| 846 |
|
|---|
| 847 | # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
|
|---|
| 848 |
|
|---|
| 849 | # default efficiency formula (misidentification rate)
|
|---|
| 850 | add EfficiencyFormula {0} {0.001}
|
|---|
| 851 | add EfficiencyFormula {4} {0.04}
|
|---|
| 852 | add EfficiencyFormula {5} {0.75}
|
|---|
| 853 | }
|
|---|
| 854 |
|
|---|
| 855 | ##########################
|
|---|
| 856 | # c-tagging
|
|---|
| 857 | ##########################
|
|---|
| 858 |
|
|---|
| 859 |
|
|---|
| 860 | module BTagging CTagging {
|
|---|
| 861 | set ParticleInputArray Delphes/allParticles
|
|---|
| 862 | set PartonInputArray Delphes/partons
|
|---|
| 863 | set JetInputArray JetEnergyScale/jets
|
|---|
| 864 |
|
|---|
| 865 | set DeltaR 0.4
|
|---|
| 866 | set BPTMin 30.0
|
|---|
| 867 | set BEtaMax 6.0
|
|---|
| 868 |
|
|---|
| 869 | # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
|
|---|
| 870 | add EfficiencyFormula {0} {0.01}
|
|---|
| 871 | add EfficiencyFormula {5} {0.10}
|
|---|
| 872 | add EfficiencyFormula {4} {0.25}
|
|---|
| 873 | }
|
|---|
| 874 |
|
|---|
| 875 | #####################################################
|
|---|
| 876 | # Find uniquely identified photons/electrons/tau/jets
|
|---|
| 877 | #####################################################
|
|---|
| 878 |
|
|---|
| 879 | module UniqueObjectFinder UniqueObjectFinder {
|
|---|
| 880 | # earlier arrays take precedence over later ones
|
|---|
| 881 | # add InputArray InputArray OutputArray
|
|---|
| 882 | add InputArray MuonIsolation/muons muons
|
|---|
| 883 | add InputArray ElectronIsolation/electrons electrons
|
|---|
| 884 | add InputArray PhotonIsolation/photons photons
|
|---|
| 885 | add InputArray JetEnergyScale/jets jets
|
|---|
| 886 | }
|
|---|
| 887 |
|
|---|
| 888 | ##################
|
|---|
| 889 | # ROOT tree writer
|
|---|
| 890 | ##################
|
|---|
| 891 |
|
|---|
| 892 | module TreeWriter TreeWriter {
|
|---|
| 893 | # add Branch InputArray BranchName BranchClass
|
|---|
| 894 |
|
|---|
| 895 | #my
|
|---|
| 896 | #add Branch TrackMerger/tracks Track Track
|
|---|
| 897 | add Branch TowerMerger/towers Tower Tower
|
|---|
| 898 | add Branch ImpactParameterSmearing/tracks TrackWithIPSmearing Track
|
|---|
| 899 | add Branch TrackPileUpSubtractor/tracksPV TrackPV Track
|
|---|
| 900 | add Branch EFlowMerger/eflow PFO Eflow
|
|---|
| 901 | #
|
|---|
| 902 | add Branch Delphes/allParticles Particle GenParticle
|
|---|
| 903 | add Branch GenJetFinder/jets GenJet Jet
|
|---|
| 904 |
|
|---|
| 905 | add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
|
|---|
| 906 | add Branch Hcal/eflowNeutralHadrons NeutralHadron Tower
|
|---|
| 907 | add Branch UniqueObjectFinder/photons Photon Photon
|
|---|
| 908 |
|
|---|
| 909 | #my
|
|---|
| 910 | #add Branch Calorimeter/eflowTracks EFlowTrack Track
|
|---|
| 911 | add Branch Ecal/eflowPhotons EFlowPhoton Tower
|
|---|
| 912 |
|
|---|
| 913 | #add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
|
|---|
| 914 | #
|
|---|
| 915 | add Branch GenJetFinder/jets GenJet Jet
|
|---|
| 916 | add Branch GenMissingET/momentum GenMissingET MissingET
|
|---|
| 917 |
|
|---|
| 918 |
|
|---|
| 919 | add Branch UniqueObjectFinder/electrons Electron Electron
|
|---|
| 920 | add Branch UniqueObjectFinder/muons Muon Muon
|
|---|
| 921 | add Branch UniqueObjectFinder/jets Jet Jet
|
|---|
| 922 | #add Branch JetEnergyScale/jets Jet Jet
|
|---|
| 923 | add Branch MissingET/momentum MissingET MissingET
|
|---|
| 924 | add Branch ScalarHT/energy ScalarHT ScalarHT
|
|---|
| 925 | add Branch Rho/rho Rho Rho
|
|---|
| 926 | add Branch PileUpMerger/vertices Vertex Vertex
|
|---|
| 927 | }
|
|---|
| 928 |
|
|---|
| 929 |
|
|---|
| 930 |
|
|---|
