Changes between Version 48 and Version 49 of WorkBook/PileUp
- Timestamp:
- Dec 18, 2014, 1:46:55 AM (10 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
WorkBook/PileUp
v48 v49 23 23 24 24 {{{ZVertexSpread}}}:: Pile-up and hard scattering events are randomly distributed in time and z position according to some parametrization specified by the user (in meters and second units). It can be either continuous or binned: 25 26 {{{ 25 {{{#!tcl 27 26 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))} 28 27 }}} 29 30 {{{ 28 {{{#!tcl 31 29 set VertexDistributionFormula {(abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + 32 30 (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + … … 43 41 By default the time information is now stored in each output collection object. However, in order to account for an actual timing measurement, the module {{{TimeSmearing}}} should be called. This module simply smears (by a Gaussian) the final time according to some resolution specified by the user, controlled by the parameter 44 42 45 {{{ 43 {{{#!tcl 46 44 set TimeResolution 1.0e-10 47 45 }}} … … 53 51 First, the user must specify whether to calculate the area while clustering the jets within the {{{FastJetFinder}}} module. Several methods for the area calculation can be specified (active area, passive area, Voronoi) via the parameter {{{AreaAlgorithm}}}. By default this parameter is set to 0 (no area calculation): 54 52 55 {{{ 53 {{{#!tcl 56 54 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area 57 55 set AreaAlgorithm 5 … … 60 58 Then the median density (in GeV/A) of pile-up contamination (rho) per event can the be computed within the {{{FastJetFinder}}} module: 61 59 62 {{{ 60 {{{#!tcl 63 61 set ComputeRho true 64 62 set RhoOutputArray rho … … 69 67 It is now possible to specify several eta ranges for the computation of the pile-up density, leading to a more accurate estimate. For instance: 70 68 71 {{{ 69 {{{#!tcl 72 70 add RhoEtaRange 0.0 2.5 73 71 add RhoEtaRange 2.5 5.0 … … 79 77 This can be done by calling the module !FastJetGridMedianEstimator. You can specify different regions as before, and, in addition, 80 78 the grid-spacing in eta, phi: 81 82 83 {{{ 79 {{{#!tcl 84 80 module FastJetGridMedianEstimator Rho { 85 81 … … 105 101 {{{}}}:: For |z| < ZVertexResolution the hard interaction vertex cannot be distinguished from pile-up vertices. 106 102 {{{}}}:: For |z| > ZVertexResolution the hard interaction vertex can be distinguished from pile-up vertices. 107 108 109 {{{ 103 {{{#!tcl 110 104 module TrackPileUpSubtractor TrackPileUpSubtractor { 111 105 … … 124 118 * Jet pile-up subtraction is done via the {{{JetPileUpSubtractor}}} module that takes as input the jet collection and rho: 125 119 126 {{{ 120 {{{#!tcl 127 121 set JetInputArray FastJetFinder/jets 128 122 set RhoInputArray rho … … 131 125 * Isolation subtraction is done inside the {{{Isolation}}} module itself just by adding the line in the Delphes card: 132 126 133 {{{ 127 {{{#!tcl 134 128 set RhoInputArray rho 135 129 }}} … … 140 134 The ID variables can be computed either using the jet constituents: 141 135 142 {{{ 136 {{{#!tcl 143 137 set UseConstituents 1 144 138 }}} … … 146 140 or using all the neutral and charged particle-flow objects in a cone around the jet axis 147 141 148 {{{ 142 {{{#!tcl 149 143 set UseConstituents 0 150 144 }}} … … 164 158 Convert your minimum bias sample into binary format: 165 159 166 {{{ 160 {{{#!sh 167 161 ./stdhep2pileup MinBias.pileup MinBias.hep 168 162 ./hepmc2pileup MinBias.pileup MinBias.hepmc … … 174 168 Run Delphes on your sample X with pile-up: 175 169 176 {{{ 170 {{{#!sh 177 171 ./DelphesSTDHEP cards/delphes_card_CMS_PileUp.tcl X_PileUp.root X.hep 178 172 ./DelphesSTDHEP cards/delphes_card_ATLAS_PileUp.tcl X_PileUp.root X.hep
