| 1 | set MaxEvents 1
|
|---|
| 2 |
|
|---|
| 3 | #######################################
|
|---|
| 4 | # Order of execution of various modules
|
|---|
| 5 | #######################################
|
|---|
| 6 |
|
|---|
| 7 | set ExecutionPath {
|
|---|
| 8 | ParticlePropagator
|
|---|
| 9 |
|
|---|
| 10 | ElectronTrackingEfficiency
|
|---|
| 11 | ElectronMomentumSmearing
|
|---|
| 12 | TreeWriter
|
|---|
| 13 | }
|
|---|
| 14 |
|
|---|
| 15 | #################################
|
|---|
| 16 | # Propagate particles in cylinder
|
|---|
| 17 | #################################
|
|---|
| 18 |
|
|---|
| 19 | module ParticlePropagator ParticlePropagator {
|
|---|
| 20 | set InputArray Delphes/stableParticles
|
|---|
| 21 |
|
|---|
| 22 | set OutputArray stableParticles
|
|---|
| 23 | set ChargedHadronOutputArray chargedHadrons
|
|---|
| 24 | set ElectronOutputArray electrons
|
|---|
| 25 | set MuonOutputArray muons
|
|---|
| 26 |
|
|---|
| 27 | # radius of the magnetic field coverage, in m
|
|---|
| 28 | set Radius 1.29
|
|---|
| 29 | # half-length of the magnetic field coverage, in m
|
|---|
| 30 | set HalfLength 3.00
|
|---|
| 31 |
|
|---|
| 32 | # magnetic field
|
|---|
| 33 | set Bz 3.8
|
|---|
| 34 | }
|
|---|
| 35 |
|
|---|
| 36 |
|
|---|
| 37 | ##############################
|
|---|
| 38 | # Electron tracking efficiency
|
|---|
| 39 | ##############################
|
|---|
| 40 |
|
|---|
| 41 | module TrackEfficiency ElectronTrackingEfficiency {
|
|---|
| 42 | set InputArray ParticlePropagator/electrons
|
|---|
| 43 | set OutputArray electrons
|
|---|
| 44 |
|
|---|
| 45 |
|
|---|
| 46 | # tracking efficiency formula (efficiency formula as a function of pt, eta, phi, energy, d0 (mm), dz (mm), ctgTheta)
|
|---|
| 47 | set EfficiencyFormula {
|
|---|
| 48 | (d0 < 1e-3) * (pt <= 0.1) * (0.00) +
|
|---|
| 49 | (d0 < 1e-3) * (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
|
|---|
| 50 | (d0 < 1e-3) * (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
|
|---|
| 51 | (d0 < 1e-3) * (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
|
|---|
| 52 | (d0 < 1e-3) * (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
|
|---|
| 53 | (d0 < 1e-3) * (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
|
|---|
| 54 | (d0 < 1e-3) * (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
|
|---|
| 55 | (d0 < 1e-3) * (abs(eta) > 2.5) * (0.00) +
|
|---|
| 56 | (d0 > 1e-3) * * (0.00)
|
|---|
| 57 | }
|
|---|
| 58 | }
|
|---|
| 59 |
|
|---|
| 60 | ###################################
|
|---|
| 61 | # Momentum resolution for electrons
|
|---|
| 62 | ###################################
|
|---|
| 63 |
|
|---|
| 64 | module MomentumSmearing ElectronMomentumSmearing {
|
|---|
| 65 | set InputArray ElectronTrackingEfficiency/electrons
|
|---|
| 66 | set OutputArray electrons
|
|---|
| 67 |
|
|---|
| 68 | # set ResolutionFormula {resolution formula as a function of eta and energy}
|
|---|
| 69 |
|
|---|
| 70 | # resolution formula for electrons
|
|---|
| 71 | # based on arXiv:1502.02701
|
|---|
| 72 | set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
|
|---|
| 73 | (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
|
|---|
| 74 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
|
|---|
| 75 | }
|
|---|
| 76 |
|
|---|
| 77 | ##################
|
|---|
| 78 | # ROOT tree writer
|
|---|
| 79 | ##################
|
|---|
| 80 |
|
|---|
| 81 | # tracks, towers and eflow objects are not stored by default in the output.
|
|---|
| 82 | # if needed (for jet constituent or other studies), uncomment the relevant
|
|---|
| 83 | # "add Branch ..." lines.
|
|---|
| 84 |
|
|---|
| 85 | module TreeWriter TreeWriter {
|
|---|
| 86 | # add Branch InputArray BranchName BranchClass
|
|---|
| 87 | add Branch Delphes/allParticles Particle GenParticle
|
|---|
| 88 | add Branch ElectronMomentumSmearing/electrons Electron Electron
|
|---|
| 89 | }
|
|---|
