1 / 41

Jets, Kinematics, and Other Variables

Jets, Kinematics, and Other Variables. A Tutorial Drew Baden University of Maryland World Scientific Int.J.Mod.Phys.A13:1817-1845,1998. Detector. x. f. q. r. Proton beam direction. Proton or anti-proton beam direction. z. Protons. Anti-Protons. Transverse E  E T. y. E.

oliver
Download Presentation

Jets, Kinematics, and Other Variables

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Jets, Kinematics, and Other Variables A Tutorial Drew Baden University of Maryland World Scientific Int.J.Mod.Phys.A13:1817-1845,1998 UCLA

  2. Detector x f q r Proton beam direction Proton or anti-proton beam direction z Protons Anti-Protons Transverse E  ET y E Coordinates UCLA

  3. Nucleon-nucleon Scattering • Forward-forward scattering, no disassociation b >> 2 rp UCLA

  4. Single-diffractive scattering • One of the 2 nucleons disassociates b ~ 2 rp UCLA

  5. Double-diffractive scattering • Both nucleons disassociates b < rp UCLA

  6. Proton-(anti)Proton Collisions • At “high” energies we are probing the nucleon structure • “High” means Ebeam >> hc/rproton ~ 1 GeV (Ebeam=1TeV@FNAL, 7TeV@LHC) • We are really doing parton–parton scattering (parton = quark, gluon) • Look for scatterings with large momentum transfer, ends up in detector “central region” (large angles wrt beam direction) • Each parton has a momentum distribution – CM of hard scattering is not fixed as in e+e- • CM of parton–parton system will be moving along z-axis with some boost • This motivates studying boosts along z • What’s “left over” from the other partons is called the “underlying event” • If no hard scattering happens, can still have disassociation • Underlying event with no hard scattering is called “minimum bias” UCLA

  7. “Total Cross-section” • By far most of the processes in nucleon-nucleon scattering are described by: • s(Total) ~ s(scattering) + s(single diffractive) + s(double diffractive) • This can be naively estimated…. • s ~ 4prp2 ~ 100mb • Total cross-section stuff is NOT the reason we do these experiments! • Examples of “interesting” physics @ Tevatron (2 TeV) • W production and decay via lepton • sBr(W en) ~ 2nb • 1 in 5x107 collisions • Z production and decay to lepton pairs • About 1/10 that of W to leptons • Top quark production • s(total) ~ 5pb • 1 in 2x1010 collisions UCLA

  8. or Phase Space • Relativistic invariant phase-space element: • Define pp/pp collision axis along z-axis: • Coordinates pm = (E,px,py,px) – Invariance with respect to boosts along z? • 2 longitudinal components: E & pz (and dpz/E) NOT invariant • 2 transverse components: px py, (and dpx, dpy) ARE invariant • Boosts along z-axis • For convenience: define pm where only 1 component is not Lorentz invariant • Choose pT, m, f as the “transverse” (invariant) coordinates • pT psin(q) and f is the azimuthal angle • For 4th coordinate define “rapidity” (y) • …How does it transform? UCLA

  9. Boosts Along beam-axis • Form a boost of velocity b along z axis • pz g(pz + bE) • E g(E+ bpz) • Transform rapidity: • Boosts along the beam axis with v=b will change y by a constant yb • (pT,y,f,m)  (pT,y+yb,f,m) with y  y+ yb , yb  ln g(1+b) simple additive to rapidity • Relationship between y, b, and q can be seen using pz = pcos(q) and p = bE or where b is the CM boost UCLA

  10. Phase Space (cont) • Transform phase space element dt from (E,px,py,pz) to (pt, y, f, m) • Gives: • Basic quantum mechanics: ds = |M |2dt • If |M |2 varies slowly with respect to rapidity, ds/dy will be ~constant in y • Origin of the “rapidity plateau” for the min bias and underlying event structure • Apply to jet fragmentation - particles should be uniform in rapidity wrt jet axis: • We expect jet fragmentation to be function of momentum perpendicular to jet axis • This is tested in detectors that have a magnetic field used to measure tracks & using UCLA

  11. Transverse Energy and Momentum Definitions • Transverse Momentum: momentum perpendicular to beam direction: • Transverse Energy defined as the energy if pz was identically 0: ETE(pz=0) • How does E and pz change with the boost along beam direction? • Using and gives • (remember boosts cause y  y + yb) • Note that the sometimes used formula is not (strictly) correct! • But it’s close – more later…. or then or which also means UCLA

  12. with and Invariant Mass M1,2 of 2 particles p1, p2 • Well defined: • Switch to pm=(pT,y,f,m) (and do some algebra…) • This gives • With bT pT/ET • Note: • For Dy  0 and Df 0, high momentum limit: M  0: angles “generate” mass • For b1 (m/p  0) This is a useful formula when analyzing data… UCLA

  13. Invariant Mass, multi particles • Extend to more than 2 particles: • In the high energy limit as m/p  0 for each particle: • Multi-particle invariant masses where each mass is negligible – no need to id • Example: t Wb and W jet+jet • Find M(jet,jet,b) by just adding the 3 2-body invariant masses • Doesn’t matter which one you call the b-jet and which the “other” jets as long as you are in the high energy limit UCLA

  14. Pseudo-rapidity UCLA

  15. CMS HCAL CMS ECAL “Pseudo”rapidity and “Real” rapidity • Definition of y: tanh(y) = b cos(q) • Can almost (but not quite) associate position in the detector (q) with rapidity (y) • But…at Tevatron and LHC, most particles in the detector (>90%) are p’s with b1 • Define “pseudo-rapidity” defined as h  y(q,b=1), or tanh(h) = cos(q) or (h=5, q=0.77°) UCLA

  16. h vs y • From tanh(h) = cos(q) = tanh(y)/b • We see that h  y • Processes “flat” in rapidity y will not be “flat” in pseudo-rapidity h 1.4 GeV p UCLA

  17. |h| – |y| and pT – Calorimater Cells • At colliders, cm can be moving with respect to detector frame • Lots of longitudinal momentum can escape down beam pipe • But transverse momentum pT is conserved in the detector • Plot h-y for constant mp, pT b(q) • For all h in DØ/CDF, can use h position to give y: • Pions: |h|-|y| < 0.1 for pT > 0.1GeV • Protons: |h|-|y| < 0.1 for pT > 2.0GeV • As b →1, y→ h (so much for “pseudo”) pT=0.1GeV DØ calorimeter cell width Dh=0.1 pT=0.2GeV pT=0.3GeV CMS HCAL cell width 0.08 CMS ECAL cell width 0.005 UCLA

  18. Rapidity “plateau” …some useful formulae… • Constant pt, rapidity plateau means ds/dy ~ k • How does that translate into ds/dh ? • Calculate dy/dh keeping m, and pt constant • After much algebra… dy/dh = b(h) • “pseudo-rapidity” plateau…only for b1 UCLA

  19. Transverse Mass UCLA

  20. Measured momentum conservation • Momentum conservation: and • What we measure using the calorimeter: and • For processes with high energy neutrinos in the final state: • We “measure” pn by “missing pT” method: • e.g. W  en or mn • Longitudinal momentum of neutrino cannot be reliably estimated • “Missing” measured longitudinal momentum also due to CM energy going down beam pipe due to the other (underlying) particles in the event • This gets a lot worse at LHC where there are multiple pp interactions per crossing • Most of the interactions don’t involve hard scattering so it looks like a busier underlying event UCLA

  21. Transverse Mass • Since we don’t measure pz of neutrino, cannot construct invariant mass of W • What measurements/constraints do we have? • Electron 4-vector • Neutrino 2-d momentum (pT) and m=0 • So construct “transverse mass” MT by: • Form “transverse” 4-momentum by ignoring pz (or set pz=0) • Form “transverse mass” from these 4-vectors: • This is equivalent to setting h1=h2=0 • For e/m and n, set me= mm= mn= 0 to get: UCLA

  22. Df=p Transverse Mass Kinematics for Ws • Transverse mass distribution? • Start with • Constrain to MW=80GeV and pT(W)=0 • cosDf = -1 • ETe = ETn • This gives you ETeETn versus Dh • Now construct transverse mass • Cleary MT=MW when Dh=0 UCLA

  23. Neutrino Rapidity • Can you constrain M(e,n) to determine the pseudo-rapidity of the n? • Would be nice, then you could veto on qn in “crack” regions • Use M(e,n) = 80GeV and • Since we know he, we know that hn=he ± Dh • Two solutions. Neutrino can be either higher or lower in rapidity than electron • Why? Because invariant mass involves the opening angle between particles. • Clean up sample of W’s by requiring both solutions are away from gaps? to get and solve for Dh: UCLA

  24. Jets UCLA

  25. and points to tower i with energy Jet Definition • How to define a “jet” using calorimeter towers so that we can use it for invariant mass calculations • And for inclusive QCD measurements (e.g. ds/dET) • QCD motivated: • Leading parton radiates gluons uniformly distributed azimuthally around jet axis • Assume zero-mass particles using calorimeter towers • 1 particle per tower • Each “particle” will have an energy perpendicular to the jet axis: • From energy conservation we expect total energy perpendicular to the jet axis to be zero on average: • Find jet axis that minimizes kT relative to that axis • Use this to define jet 4-vector from calorimeter towers • Since calorimeter towers measure total energy, make a basic assumption: • Energy of tower is from a single particle with that energy • Assume zero mass particle (assume it’s a pion and you will be right >90%!) • Momentum of the particle is then given by • Note: mi=0 does NOT mean Mjet=0 • Mass of jet is determined by opening angle between all contributors • Can see this in case of 2 “massless” particles, or energy in only 2 towers: • Mass is “generated” by opening angles. • A rule of thumb: Zero mass parents of decay have q12=0 always UCLA

  26. Quasi-analytical approach • Transform each calorimeter tower to frame of jet and minimize kT • 2-d Euler rotation (in picture, f=fjet, q=qjet, set c=0) • Tower in jet momentum frame: and apply • Check: for 1 tower, ftower=fjet, should get Exi = Exi = 0 and Ezi = Ejet • It does, after some algebra… UCLA

  27. Minimize kT to Find Jet Axis • The equation is equivalent to so… • Momentum of the jet is such that: UCLA

  28. Jet 4-momentum summary • Jet Energy: • Jet Momentum: • Jet Mass: • Jet 4-vector: • Jet is an object now! So how do we define ET? UCLA

  29. ET of a Jet • For any object, ET is well defined: • There are 2 more ways you could imagine using to define ET of a jet but neither are technically correct: • How do they compare? • Is there any ET or h dependence? correct Alternative 1 Alternative 2 or UCLA

  30. True ET vs Alternative 1 • True: • Alternative 1: • Define which is always >0 • Expand in powers of : • For small h, tanhh h so either way is fine • Alternative 1 is the equivalent to true def central jets • Agree at few% level for |h|<0.5 • For h~0.5 or greater....cone dependent • Or “mass” dependent....same thing Leading jet, |h|>0.5 UCLA

  31. True ET vs Alternative 2 Alternative 2: harder to see analytically…imagine a jet w/2 towers • TRUE: • Alternative 2: • Take difference: • So this method also underestimates “true” ET • But not as much as Alternative 1 Leading jet, |h|>0.5 Always > 0! UCLA

  32. Jet Shape • Jets are defined by but the “shape” is determined by • From Euler: • Now form for those towers close to the jet axis: dq0 and df 0 • From we get which means So… and… UCLA

  33. Jet Shape – ET Weighted • Define and • This gives: and equivalently, • Momentum of each “cell” perpendicular to jet momentum is from • Eti of particle in the detector, and • Distance from jet in hf plane • This also suggests jet shape should be roughly circular in hf plane • Providing above approximations are indicative overall…. • Shape defined: • Use energy weighting to calculate true 2nd moment in hf plane with UCLA

  34. Jet Shape – ET Weighted (cont) • Use sample of “unmerged” jets • Plot • Shape depends on cone parameter • Mean and widths scale linearly with cone parameter <sR> vs Jet Clustering Parameter (Cone Size) • “Small angle” approximation pretty good • For Cone=0.7, distribution in sR has: • Mean ± Width =.25 ± .05 • 99% of jets have sR <0.4 <sR> UCLA

  35. Jet Mass UCLA

  36. Jet Samples • Run 1 • All pathologies eliminated (Main Ring, Hot Cells, etc.) • |Zvtx|<60cm • No t, e, or g candidates in event • Checked hf coords of teg vs. jet list • Cut on cone size for jets • .025, .040, .060 for jets from cone cuttoff 0.3, 0.5, 0.7 respectively • “UNMERGED” Sample: • RECO events had 2 and only 2 jets for cones .3, .5, and .7 • Bias against merged jets but they can still be there • e.g. if merging for all cones • “MERGED” Sample: • Jet algorithm reports merging UCLA

  37. Jet Mass • Jet is a physics object, so mass is calculated using: • Either one… • Note: there is no such thing as “transverse mass” for a jet • Transverse mass is only defined for pairs (or more) of 4-vectors… • For large ET,jet we can see what happens by writing • And take limit as jet narrows and and expand ET and pT • This gives so….  using Jet mass is related to jet shape!!! (in the thin jet, high energy limit) UCLA

  38. Jet Mass (cont) • Jet Mass for unmerged sample How good is “thin jet” approximation? Low-side tail is due to lower ET jets for smaller cones (this sample has 2 and only 2 jets for all cones) UCLA

  39. Jet Merging • Does jet merging matter for physics? • For some inclusive QCD studies, it doesn’t matter • For invariant mass calculations from e.g. topWb, it will smear out mass distribution if merging two “tree-level” jets that happen to be close • Study sR…see clear correlation between sR and whether jet is merged or not • Can this be used to construct some kind of likelihood? “Unmerged”, Jet Algorithm reports merging, all cone sizes “Unmerged” v. “Merged” sample UCLA

  40. Merging Likelihood • Crude attempt at a likelihood • Can see that for this (biased) sample, can use this to pick out “unmerged” jets based on shape • Might be useful in Higgs search for H bb jet invariant mass? UCLA

  41. Merged Shape • Width in hf “assumes” circular • Large deviations due to merging? • Define should be independent of cone size • Clear broadening seen – “cigar”-shaped jets, maybe study… “Unmerged” Sample “Merged” Sample UCLA

More Related