| 1 | set MaxEvents 100
|
|---|
| 2 | set RandomSeed 123
|
|---|
| 3 |
|
|---|
| 4 |
|
|---|
| 5 | #
|
|---|
| 6 | # Phase II - Pile-Up
|
|---|
| 7 | #
|
|---|
| 8 | # Main authors: Michele Selvaggi (UCL)
|
|---|
| 9 | #
|
|---|
| 10 | # Released on:
|
|---|
| 11 | #
|
|---|
| 12 | # Version: v02 beta - test TrackSmearing and vertexing
|
|---|
| 13 | #
|
|---|
| 14 | #
|
|---|
| 15 | #######################################
|
|---|
| 16 | # Order of execution of various modules
|
|---|
| 17 | #######################################
|
|---|
| 18 |
|
|---|
| 19 | set ExecutionPath {
|
|---|
| 20 |
|
|---|
| 21 | BeamSpotFilter
|
|---|
| 22 |
|
|---|
| 23 | PileUpMerger
|
|---|
| 24 | ParticlePropagator
|
|---|
| 25 |
|
|---|
| 26 | ChargedHadronTrackingEfficiency
|
|---|
| 27 | ElectronTrackingEfficiency
|
|---|
| 28 | MuonTrackingEfficiency
|
|---|
| 29 |
|
|---|
| 30 | TrackMerger
|
|---|
| 31 | TrackSmearing
|
|---|
| 32 |
|
|---|
| 33 | TimeSmearing
|
|---|
| 34 |
|
|---|
| 35 | VertexFinderDA4D
|
|---|
| 36 | HighMassVertexRecover
|
|---|
| 37 |
|
|---|
| 38 | CandidateFilter
|
|---|
| 39 |
|
|---|
| 40 | TreeWriter
|
|---|
| 41 | }
|
|---|
| 42 |
|
|---|
| 43 | #######################
|
|---|
| 44 | # GenBeamSpotFilter
|
|---|
| 45 | # Saves a particle intended to represent the beamspot
|
|---|
| 46 | #######################
|
|---|
| 47 |
|
|---|
| 48 | module BeamSpotFilter BeamSpotFilter {
|
|---|
| 49 | set InputArray Delphes/stableParticles
|
|---|
| 50 | set OutputArray beamSpotParticle
|
|---|
| 51 |
|
|---|
| 52 | }
|
|---|
| 53 |
|
|---|
| 54 | ###############
|
|---|
| 55 | # PileUp Merger
|
|---|
| 56 | ###############
|
|---|
| 57 |
|
|---|
| 58 | module PileUpMerger PileUpMerger {
|
|---|
| 59 | set InputArray Delphes/stableParticles
|
|---|
| 60 |
|
|---|
| 61 | set ParticleOutputArray stableParticles
|
|---|
| 62 | set VertexOutputArray vertices
|
|---|
| 63 |
|
|---|
| 64 | # pre-generated minbias input file
|
|---|
| 65 | set PileUpFile MinBias.pileup
|
|---|
| 66 |
|
|---|
| 67 | # average expected pile up
|
|---|
| 68 | set MeanPileUp 5
|
|---|
| 69 |
|
|---|
| 70 | # 0-poisson, 1-uniform, 2-delta
|
|---|
| 71 | set PileUpDistribution 2
|
|---|
| 72 |
|
|---|
| 73 | # maximum spread in the beam direction in m
|
|---|
| 74 | set ZVertexSpread 0.25
|
|---|
| 75 |
|
|---|
| 76 | # maximum spread in time in s
|
|---|
| 77 | set TVertexSpread 800E-12
|
|---|
| 78 |
|
|---|
| 79 | # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
|
|---|
| 80 |
|
|---|
| 81 | #set VertexDistributionFormula {exp(-(t^2/(2*(0.063/2.99792458E8*exp(-(z^2/(2*(0.063)^2))))^2)))}
|
|---|
| 82 | set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
|
|---|
| 83 |
|
|---|
| 84 | # taking 5.3 cm x 160 ps
|
|---|
| 85 |
|
|---|
| 86 | #set VertexDistributionFormula { (abs(t) <= 160e-12) * (abs(z) <= 0.053) * (1.00) +
|
|---|
| 87 | # (abs(t) > 160e-12) * (abs(z) <= 0.053) * (0.00) +
|
|---|
| 88 | # (abs(t) <= 160e-12) * (abs(z) > 0.053) * (0.00) +
|
|---|
| 89 | # (abs(t) > 160e-12) * (abs(z) > 0.053) * (0.00)}
|
|---|
| 90 |
|
|---|
| 91 | }
|
|---|
| 92 |
|
|---|
| 93 | #################################
|
|---|
| 94 | # Propagate particles in cylinder
|
|---|
| 95 | #################################
|
|---|
| 96 |
|
|---|
| 97 | module ParticlePropagator ParticlePropagator {
|
|---|
| 98 | set InputArray PileUpMerger/stableParticles
|
|---|
| 99 |
|
|---|
| 100 | set OutputArray stableParticles
|
|---|
| 101 | set ChargedHadronOutputArray chargedHadrons
|
|---|
| 102 | set ElectronOutputArray electrons
|
|---|
| 103 | set MuonOutputArray muons
|
|---|
| 104 |
|
|---|
| 105 | # radius of the magnetic field coverage, in m
|
|---|
| 106 | set Radius 1.29
|
|---|
| 107 | # half-length of the magnetic field coverage, in m
|
|---|
| 108 | set HalfLength 3.0
|
|---|
| 109 |
|
|---|
| 110 | # magnetic field
|
|---|
| 111 | set Bz 3.8
|
|---|
| 112 | }
|
|---|
| 113 |
|
|---|
| 114 | ####################################
|
|---|
| 115 | # Charged hadron tracking efficiency
|
|---|
| 116 | ####################################
|
|---|
| 117 |
|
|---|
| 118 | module Efficiency ChargedHadronTrackingEfficiency {
|
|---|
| 119 | set InputArray ParticlePropagator/chargedHadrons
|
|---|
| 120 | set OutputArray chargedHadrons
|
|---|
| 121 |
|
|---|
| 122 | # add EfficiencyFormula {efficiency formula as a function of eta and pt}
|
|---|
| 123 |
|
|---|
| 124 | # tracking efficiency formula for charged hadrons
|
|---|
| 125 | set EfficiencyFormula { (pt <= 0.1) * (0.00) +
|
|---|
| 126 | (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
|
|---|
| 127 | (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
|
|---|
| 128 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
|
|---|
| 129 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
|
|---|
| 130 | (abs(eta) > 2.5) * (0.00)}
|
|---|
| 131 | }
|
|---|
| 132 |
|
|---|
| 133 | ##############################
|
|---|
| 134 | # Electron tracking efficiency
|
|---|
| 135 | ##############################
|
|---|
| 136 |
|
|---|
| 137 | module Efficiency ElectronTrackingEfficiency {
|
|---|
| 138 | set InputArray ParticlePropagator/electrons
|
|---|
| 139 | set OutputArray electrons
|
|---|
| 140 |
|
|---|
| 141 | # set EfficiencyFormula {efficiency formula as a function of eta and pt}
|
|---|
| 142 |
|
|---|
| 143 | # tracking efficiency formula for electrons
|
|---|
| 144 | set EfficiencyFormula { (pt <= 0.1) * (0.00) +
|
|---|
| 145 | (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
|
|---|
| 146 | (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
|
|---|
| 147 | (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
|
|---|
| 148 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
|
|---|
| 149 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
|
|---|
| 150 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
|
|---|
| 151 | (abs(eta) > 2.5) * (0.00)}
|
|---|
| 152 | }
|
|---|
| 153 |
|
|---|
| 154 | ##########################
|
|---|
| 155 | # Muon tracking efficiency
|
|---|
| 156 | ##########################
|
|---|
| 157 |
|
|---|
| 158 | module Efficiency MuonTrackingEfficiency {
|
|---|
| 159 | set InputArray ParticlePropagator/muons
|
|---|
| 160 | set OutputArray muons
|
|---|
| 161 |
|
|---|
| 162 | # set EfficiencyFormula {efficiency formula as a function of eta and pt}
|
|---|
| 163 |
|
|---|
| 164 | # tracking efficiency formula for muons
|
|---|
| 165 | set EfficiencyFormula { (pt <= 0.1) * (0.00) +
|
|---|
| 166 | (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
|
|---|
| 167 | (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
|
|---|
| 168 | (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
|
|---|
| 169 |
|
|---|
| 170 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
|
|---|
| 171 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
|
|---|
| 172 | (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
|
|---|
| 173 | (abs(eta) > 2.5) * (0.00)}
|
|---|
| 174 | }
|
|---|
| 175 |
|
|---|
| 176 | ##############
|
|---|
| 177 | # Track merger
|
|---|
| 178 | ##############
|
|---|
| 179 |
|
|---|
| 180 | module Merger TrackMerger {
|
|---|
| 181 | # add InputArray InputArray
|
|---|
| 182 | add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
|
|---|
| 183 | add InputArray ElectronTrackingEfficiency/electrons
|
|---|
| 184 | add InputArray MuonTrackingEfficiency/muons
|
|---|
| 185 | set OutputArray tracks
|
|---|
| 186 | }
|
|---|
| 187 |
|
|---|
| 188 | ########################################
|
|---|
| 189 | # Smear tracks
|
|---|
| 190 | ########################################
|
|---|
| 191 |
|
|---|
| 192 | module TrackSmearing TrackSmearing {
|
|---|
| 193 | set InputArray TrackMerger/tracks
|
|---|
| 194 | set OutputArray tracks
|
|---|
| 195 | set ApplyToPileUp true
|
|---|
| 196 |
|
|---|
| 197 | # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
|
|---|
| 198 | source trackResolutionCMS.tcl
|
|---|
| 199 | }
|
|---|
| 200 |
|
|---|
| 201 | module TimeSmearing TimeSmearing {
|
|---|
| 202 | set InputArray TrackSmearing/tracks
|
|---|
| 203 | set OutputArray tracks
|
|---|
| 204 |
|
|---|
| 205 | # assume 20 ps resolution for now
|
|---|
| 206 | set TimeResolution 20E-12
|
|---|
| 207 | }
|
|---|
| 208 |
|
|---|
| 209 | ##################################
|
|---|
| 210 | # Primary vertex reconstruction
|
|---|
| 211 | ##################################
|
|---|
| 212 |
|
|---|
| 213 |
|
|---|
| 214 | module VertexFinderDA4D VertexFinderDA4D {
|
|---|
| 215 | set InputArray TimeSmearing/tracks
|
|---|
| 216 |
|
|---|
| 217 | set OutputArray tracks
|
|---|
| 218 | set VertexOutputArray vertices
|
|---|
| 219 |
|
|---|
| 220 | set Verbose 0
|
|---|
| 221 | set MinPT 1.0
|
|---|
| 222 |
|
|---|
| 223 | # in mm
|
|---|
| 224 | set VertexSpaceSize 0.5
|
|---|
| 225 |
|
|---|
| 226 | # in s
|
|---|
| 227 | set VertexTimeSize 10E-12
|
|---|
| 228 |
|
|---|
| 229 | set UseTc 1
|
|---|
| 230 | set BetaMax 0.1
|
|---|
| 231 | set BetaStop 1.0
|
|---|
| 232 | set CoolingFactor 0.8
|
|---|
| 233 | set MaxIterations 100
|
|---|
| 234 |
|
|---|
| 235 | # in mm
|
|---|
| 236 | set DzCutOff 40
|
|---|
| 237 | set D0CutOff 30
|
|---|
| 238 |
|
|---|
| 239 | }
|
|---|
| 240 |
|
|---|
| 241 | ######################################
|
|---|
| 242 | # Heavy(slow) particles vertex recover
|
|---|
| 243 | ######################################
|
|---|
| 244 |
|
|---|
| 245 | module HighMassVertexRecover HighMassVertexRecover {
|
|---|
| 246 | set TrackInputArray VertexFinderDAClusterizerZT/tracks
|
|---|
| 247 | set VertexInputArray VertexFinderDAClusterizerZT/vertices
|
|---|
| 248 |
|
|---|
| 249 | set TrackOutputArray tracks
|
|---|
| 250 | set VertexOutputArray vertices
|
|---|
| 251 |
|
|---|
| 252 | set Verbose 0
|
|---|
| 253 | }
|
|---|
| 254 |
|
|---|
| 255 |
|
|---|
| 256 | ########################################
|
|---|
| 257 | # Remove uninteresting tracks #
|
|---|
| 258 | ########################################
|
|---|
| 259 |
|
|---|
| 260 | module CandidateFilter CandidateFilter {
|
|---|
| 261 | set InputArray HighMassVertexRecover/tracks
|
|---|
| 262 | set OutputArray tracks
|
|---|
| 263 |
|
|---|
| 264 | set PtMin 10
|
|---|
| 265 | set MassMin 0.2
|
|---|
| 266 | }
|
|---|
| 267 |
|
|---|
| 268 |
|
|---|
| 269 | ##################
|
|---|
| 270 | # ROOT tree writer
|
|---|
| 271 | ##################
|
|---|
| 272 |
|
|---|
| 273 | module TreeWriter TreeWriter {
|
|---|
| 274 | # add Branch InputArray BranchName BranchClass
|
|---|
| 275 | add Branch PileUpMerger/stableParticles Particle GenParticle
|
|---|
| 276 | add Branch PileUpMerger/vertices GenVertex Vertex
|
|---|
| 277 | add Branch TimeSmearing/tracks Track Track
|
|---|
| 278 | add Branch VertexFinderDA4D/vertices Vertex4D Vertex
|
|---|
| 279 | add Branch HighMassVertexRecover/vertices Vertex4D Vertex
|
|---|
| 280 | add Branch CandidateFilter/tracks Track Track
|
|---|
| 281 | }
|
|---|
