Jet Physics - Experimental Aspects - Murilo Rangel ! ! New Trends in HEP and QCD School Outline
This lecture is an attempt to cover the experimental challenges of jet physics at LHC ! How to go from detector hits to theory comparison? ! What are the new ideas for jet tools? ! Notes: 1) jet algorithms⨁theory are not covered in this lectures 2) nice lectures already given (overlap is expected) Grégory Soyez, Albert de Roeck, Bruno Lenzi, Rikkert Frederix, Marc Besançon 3) some experimental aspects are not covered (ex: efficiencies, trigger)
Murilo Rangel - IF/UFRJ 2 Reminder
Reminder
TeV4LHC QCD WG - hep-ph/0610012
Murilo Rangel - IF/UFRJ 3 Is this really what is happening?
Murilo Rangel - IF/UFRJ 4 Epistemological Realism - Personal View
Truth is a place we can not go
Murilo Rangel - IF/UFRJ 5 Epistemological Realism - Personal View
Truth is a place we can not go But, we can take pictures of it
Murilo Rangel - IF/UFRJ 5 Epistemological Realism - Personal View
Truth is a place we can not go
Murilo Rangel - IF/UFRJ 5 Epistemological Realism - Personal View
Truth is a place we can not go And, we can paint how we think it is
Murilo Rangel - IF/UFRJ 5 Epistemological Realism - Personal View
Truth is a place we can not go
Murilo Rangel - IF/UFRJ 5 Epistemological Realism - Personal View
Truth is a place we can not go
Our job is to compare photographies with paintings Both photographers and painters are doing a great job
Murilo Rangel - IF/UFRJ 5 Jet evolution
‘’stable’' particles
decays
hadrons time hadronization
partons
proton proton
Murilo Rangel - IF/UFRJ 6 Jet Clustering
+ Algorithms that combine nearest particles ! o Cambridge/Aachen algorithm: combine particles nearest each other ! o “kT” algorithm: preference for combining lower-momentum particle pairs first, then moving on to higher-momentum pairs ! o “anti-kt” algorithm collects particles around the hardest particle first. It guarantees “cone-like geometry” with well-defined borders around the highest-kT particles and it maintains the infrared safety and collinear safety of sequential recombination family ! + These algorithms correspond to p=0, p=1 and p=-1. ! 2 2 2 ∆ ! d =min(k p,k p) ij ij ti tj R2
Murilo Rangel - IF/UFRJ 7 Jet evolution
“hadron-level” / “particle-level” jets or “stable” particles* used as inputs of the jet algorithm ‘’stable’' particles
decays
hadrons time hadronization
“parton-level” jets partons used as inputs of the jet algorithm
partons
proton proton
*Neutrinos are excluded Murilo Rangel - IF/UFRJ 8 Experimental Challenge
- Only “stable’’ particles are detected (� >10-8 s)* - Prior knowledge of their interactions are needed
http://www.particleadventure.org/
0 *These include �, K, p, n, e, �, � and KL . Neutrinos are invisible
Murilo Rangel - IF/UFRJ 9 Particle Interaction
R. Cavanaugh, HCPSS 2012
Murilo Rangel - IF/UFRJ 10 Particle Interaction
Nuclear Interaction Length Mean distance over which a hadron collides with a nuclei
� ~ 35 g cm-2 A1/3
R. Cavanaugh, HCPSS 2012
Murilo Rangel - IF/UFRJ 10 Particle Interaction
Nuclear Interaction Length Mean distance over which a hadron collides with a nuclei
� ~ 35 g cm-2 A1/3
Radiation length Mean distance over which the electron energy is reduced by a factor of 1/e due to radiation losses
R. Cavanaugh, HCPSS 2012
Murilo Rangel - IF/UFRJ 10 Particle Interaction
Nuclear Interaction Length Mean distance over which a hadron collides with a nuclei
� ~ 35 g cm-2 A1/3
Typical values for lead
� ~ 17 cm
X0 ~ 5.5 mm
Radiation length Mean distance over which the electron energy is reduced by a factor of 1/e due to radiation losses
R. Cavanaugh, HCPSS 2012
Murilo Rangel - IF/UFRJ 10 Calorimeters
a is the stochastic term b is the electronic noise c includes detector effects
Hadronic Calorimeter - high density material - total absorption of energy - � = 50-100% / √ET
EM Calorimeter - high-Z material - total absorption of energy - � = 1-10% / √ET
Murilo Rangel - IF/UFRJ 11 Calorimeters vs Tracker
R. Cavanaugh, HCPSS 2012
Murilo Rangel - IF/UFRJ 12 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! Muons
http://www.particleadventure.org/
Murilo Rangel - IF/UFRJ 13 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! Electrons / Photons
http://www.particleadventure.org/
Murilo Rangel - IF/UFRJ 14 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! Hadrons
http://www.particleadventure.org/
Murilo Rangel - IF/UFRJ 15 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! Converted photons
location of converted photons
Murilo Rangel - IF/UFRJ 16 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! Converted photons
location of converted photons
Murilo Rangel - IF/UFRJ 16 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! V0s (� ~ 10-10 s)
V0s live long enough to reconstruct its vertex
Murilo Rangel - IF/UFRJ 17 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! V0s (� ~ 10-10 s)
V0s live long enough to reconstruct its vertex
Murilo Rangel - IF/UFRJ 17 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! π0→γγ
Murilo Rangel - IF/UFRJ 18 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! Neutral Hadrons
neutral hadron
charged hadron
Primary Vertex
Murilo Rangel - IF/UFRJ 19 Particle Flow
Strategy to get most of detector is to match tracks with calorimeter clusters ! Track momentum is preferred over calorimeter energy ! Steps are ordered motivated by momentum resolution and particle identification purity ! Neutral Hadrons
neutral hadron
charged hadron
Primary Vertex
Murilo Rangel - IF/UFRJ 19 Particle Flow
Murilo Rangel - IF/UFRJ 20 Jets - Experimental Picture
Reconstructed “particles” are used as inputs of the jet algorithm - “detector-level” jets
“Detector-level” jets must be corrected/calibrated to compare with theory/models Calibration of jets to “particle-level” is necessary ! Calibrated jets - little dependence with detector effects (segmentation, response and resolution) - good resolution and no angle biases - good efficiency and low fake rate (Jet Identification) - stable with beam luminosity (pile-up) - computer time efficient ! Inputs - calorimeter cells/towers/clusters - tracks - tracks+calorimeter (particle flow) ! !
Murilo Rangel - IF/UFRJ 21 Jet Identification
At “detector-level”, the jet algorithm can reconstruct fake jet candidates: - Hadronic tau decays (electrons and photons too) - Cosmic ray - Detector noise - Pile-Up contribution
Murilo Rangel - IF/UFRJ 22 Jet Identification
At “detector-level”, the jet algorithm can reconstruct fake jet candidates: - Hadronic tau decays (electrons and photons too) - Cosmic ray - Detector noise - Pile-Up contribution
Jet characteristics can be used to reduce these background rate to O(1%) o Charged Fraction (from PV) o EM Fraction
Murilo Rangel - IF/UFRJ 22 Jet Identification
At “detector-level”, the jet algorithm can reconstruct fake jet candidates: - Hadronic tau decays (electrons and photons too) - Cosmic ray - Detector noise - Pile-Up contribution
Jet characteristics can be used to reduce these background rate to O(1%) o Charged Fraction (from PV) o EM Fraction
Fake jets sample ATLAS Preliminary 5 after Looser cuts 10 Data 2011, s = 7 TeV after Loose cuts after Medium cuts 4 pjet>150 GeV, |η|<2 after Tight cuts 10 T Good jets sample
103 Number of jets / 0.05
102
10
1 0 0.2 0.4 0.6 0.8 1 1.2 1.4
fch
Murilo Rangel - IF/UFRJ 22 Jet Identification
At “detector-level”, the jet algorithm can reconstruct fake jet candidates: - Hadronic tau decays (electrons and photons too) - Cosmic ray - Detector noise - Pile-Up contribution
Jet characteristics can be used to reduce these background rate to O(1%) o Charged Fraction (from PV) o EM Fraction
Fake jets sample Fake jets sample ATLAS Preliminary 8 ATLAS Preliminary 5 after Looser cuts 10 after Looser cuts 10 Data 2011, s = 7 TeV after Loose cuts Data 2011, s = 7 TeV after Loose cuts 107 after Medium cuts jet after Medium cuts jet p >150 GeV 4 p >150 GeV, |η|<2 after Tight cuts 6 T after Tight cuts 10 T Good jets sample 10 Good jets sample 105 103
Number of jets / 0.05 Number of jets / 0.05 104 2 10 103 102 10 10 1 1 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1
fch fEM
Murilo Rangel - IF/UFRJ 22 Jet Identification
At “detector-level”, the jet algorithm can reconstruct fake jet candidates: - Hadronic tau decays (electrons and photons too)
-1 - Cosmic ray ×103CMS Preliminary, s = 8TeV L=20 fb Z→µµ Data - Detector noise 600 All |η| < 2.5 Jet p > 25 GeV Gluon T Quark - Pile-Up contribution 500 PU Real Jet Events/0.025 400 Jet characteristics can be used to reduce these 300 background rate to O(1%) 200 o Charged Fraction (from PV) 100 o EM Fraction 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 β
Fake jets sample Fake jets sample ATLAS Preliminary 8 ATLAS Preliminary 5 after Looser cuts 10 after Looser cuts 10 Data 2011, s = 7 TeV after Loose cuts Data 2011, s = 7 TeV after Loose cuts 107 after Medium cuts jet after Medium cuts jet p >150 GeV 4 p >150 GeV, |η|<2 after Tight cuts 6 T after Tight cuts 10 T Good jets sample 10 Good jets sample 105 103
Number of jets / 0.05 Number of jets / 0.05 104 2 10 103 102 10 10 1 1 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1
fch fEM
Murilo Rangel - IF/UFRJ 22 Jet Composition
LHCb Generator
Jet at “particle level”
LHCb Simulation
Jet at "detector level” (uncorrected)
Murilo Rangel - IF/UFRJ 23 Jet Composition
Tests of parton-shower⨁hadronization models are necessary
particle energy fraction
Murilo Rangel - IF/UFRJ 24 Calibration Factorization
Atlas
CMS
Murilo Rangel - IF/UFRJ 25 Calibration Factorization
Atlas
CMS
Murilo Rangel - IF/UFRJ 25 Pile Up
Murilo Rangel - IF/UFRJ 26 Pile Up
CMS 78 collisions
Murilo Rangel - IF/UFRJ 27 Pile up
Pile Up Vertex Primary Vertex
Murilo Rangel - IF/UFRJ 28 Pile up Correction
In-time pile up activity depends on the number of Primary Vertices (PVs) Out-of-time pile up activity depends on the average number of PVs
larger R, larger PU
Murilo Rangel - IF/UFRJ 29 Pile up Correction
Jet “independent” PU correction is also possible, e.g., jet area method. ➡ adding “infinite” number of very soft 4-momentum vectors to cluster jets j jet area is defined (A ) as the space occupied by the very soft particles j j ➡ distribution of pT /A is related to the PU activity
NPVs=1
Number of PVs
Murilo Rangel - IF/UFRJ 30 Jet Energy Calibration - Simple view
Using simulation, we can match jets at “particle-level” and at “detector-level” ! ! ! ! ! ! Calibration factor is taken from the ratio pT(detector-level)/pT(particle-level) + factor is applied to 4-momentum: angle biases needed to be checked !
100 Counts 80 If your detector is perfect => delta function
60
40
20
0 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 p (detector-level) / p(particle-level) T T
Murilo Rangel - IF/UFRJ 31 Jet Energy Calibration - Samples
We do not want to be simulation dependent. ! Data-driven methods are developed using production of well calibrated object with a jet photon+jet or Z(→µ+µ-/e+e-)+jet ! ! ! ! ! ! ! ! ! Two jet production are also very useful for relative jet energy calibration
Murilo Rangel - IF/UFRJ 32 Jet Energy Calibration - Dijet
Dijet sample ➡ Both objects are the subject for calibration - this sample can be used to calibrate one region of the detector relative to another one. ➡ Jet energy (pT) resolution can be measured
Murilo Rangel - IF/UFRJ 33 Jet Energy Calibration - V+jet
�/Z+jet The reference object energy resolution is much better than the jet energy resolution. ! ➡ At LO, the �/Z is balanced with the parton that originates the jet. ! ➡ Missing transverse energy projection fraction (MPF) is used to include effects like: - additional parton radiation - underlying-event (UE) contribution - out-of-cone contribution
Murilo Rangel - IF/UFRJ 34 Jet Energy Calibration - V+jet
�/Z+jet The reference object energy resolution is much better than the jet energy resolution. ! ➡ At LO, the �/Z is balanced with the parton that originates the jet. ! ➡ Missing transverse energy projection fraction (MPF) is used to include effects like: - additional parton radiation - underlying-event (UE) contribution - out-of-cone contribution
Murilo Rangel - IF/UFRJ 34 Jet Energy Calibration - V+jet
�/Z+jet The reference object energy resolution is much better than the jet energy resolution. ! ➡ At LO, the �/Z is balanced with the parton that originates the jet. ! ➡ Missing transverse energy projection fraction (MPF) is used to include effects like: - additional parton radiation - underlying-event (UE) contribution - out-of-cone contribution
Murilo Rangel - IF/UFRJ 34 Jet Energy Calibration - V+jet
�/Z+jet The reference object energy resolution is much better than the jet energy resolution. ! ➡ At LO, the �/Z is balanced with the parton that originates the jet. ! ➡ Missing transverse energy projection fraction (MPF) is used to include effects like: - additional parton radiation - underlying-event (UE) contribution - out-of-cone contribution
Murilo Rangel - IF/UFRJ 34 Jet Energy Calibration - V+jet
�/Z+jet The reference object energy resolution is much better than the jet energy resolution. ! ➡ At LO, the �/Z is balanced with the parton that originates the jet. ! ➡ Missing transverse energy projection fraction (MPF) is used to include effects like: - additional parton radiation - underlying-event (UE) contribution - out-of-cone contribution
Murilo Rangel - IF/UFRJ 34 Jet Energy Calibration - V+jet
�/Z+jet The reference object energy resolution is much better than the jet energy resolution. ! ➡ At LO, the �/Z is balanced with the parton that originates the jet. ! ➡ Missing transverse energy projection fraction (MPF) is used to include effects like: - additional parton radiation - underlying-event (UE) contribution - out-of-cone contribution
Murilo Rangel - IF/UFRJ 34 Jet Energy Calibration - ISR/FSR
Jet2 γ By vetoing additional jets in the sample (pT >α pT ), the effect of initial and final parton radiation can be studied
leading jet
2nd jet
Photon
Murilo Rangel - IF/UFRJ 35 Jet Energy Calibration - ISR/FSR
Jet2 γ By vetoing additional jets in the sample (pT >α pT ), the effect of initial and final parton radiation can be studied
Murilo Rangel - IF/UFRJ 36 Residual Correction
Calibration derived in data may need to be corrected by residual effects ➾ data-to-MC differences ➾ different MC can provide different corrections
Murilo Rangel - IF/UFRJ 37 Sanity check using W→jets
Good knowledge of the W boson mass can be used to test the jet energy calibration ➾ W from top quark decay ➾ sensitive to jets originating from quarks
Murilo Rangel - IF/UFRJ 38 Jet Calibration - Flavour dependence
How can the calibration vary by changing the initial parton flavour (gluon)?