350 likes | 362 Views
This workshop explores the early days of Feynman-Field phenomenology in studying high energy hadron collisions, including the analysis of "min-bias" collisions and the tuning of QCD Monte-Carlo models. It also discusses the predictions of the Quark-Quark Black-Box Model and the development of the Feynman-Field Jet Model.
E N D
1st Workshop on Energy Scalingin Hadron-Hadron Collisions Rick’s View of Hadron Collisions Fermilab 2009 Rick Field University of Florida Outline of Talk • The early days of Feynman-Field Phenomenology. • Studying “min-bias” collisions and the “underlying event” in Run 1 at CDF. • Tuning the QCD Monte-Carlo model generators. • Studying the “associated” charged particle densities in “min-bias” collisions. CDF Run 2 CMS at the LHC Rick Field – Florida/CDF/CMS
Feynman-Field Phenomenology Toward and Understanding of Hadron-Hadron Collisions 1st hat! Feynman and Field • From 7 GeV/c p0’s to 600 GeV/c Jets. The early days of trying to understand and simulate hadron-hadron collisions. Rick Field – Florida/CDF/CMS
Hadron-Hadron Collisions FF1 1977 • What happens when two hadrons collide at high energy? Feynman quote from FF1 “The model we shall choose is not a popular one, so that we will not duplicate too much of the work of others who are similarly analyzing various models (e.g. constituent interchange model, multiperipheral models, etc.). We shall assume that the high PT particles arise from direct hard collisions between constituent quarks in the incoming particles, which fragment or cascade down into several hadrons.” • Most of the time the hadrons ooze through each other and fall apart (i.e.no hard scattering). The outgoing particles continue in roughly the same direction as initial proton and antiproton. • Occasionally there will be a large transverse momentum meson. Question: Where did it come from? • We assumed it came from quark-quark elastic scattering, but we did not know how to calculate it! “Black-Box Model” Rick Field – Florida/CDF/CMS
Quark-Quark Black-Box Model No gluons! FF1 1977 Quark Distribution Functions determined from deep-inelastic lepton-hadron collisions Feynman quote from FF1 “Because of the incomplete knowledge of our functions some things can be predicted with more certainty than others. Those experimental results that are not well predicted can be “used up” to determine these functions in greater detail to permit better predictions of further experiments. Our papers will be a bit long because we wish to discuss this interplay in detail.” Quark Fragmentation Functions determined from e+e- annihilations Quark-Quark Cross-Section Unknown! Deteremined from hadron-hadron collisions. Rick Field – Florida/CDF/CMS
Quark-Quark Black-Box Model FF1 1977 Predict increase with increasing CM energy W Predict particle ratios “Beam-Beam Remnants” Predict overall event topology (FFF1 paper 1977) 7 GeV/c p0’s! Rick Field – Florida/CDF/CMS
Feynman Talk at Coral Gables(December 1976) 1st transparency Last transparency “Feynman-Field Jet Model” Rick Field – Florida/CDF/CMS
QCD Approach: Quarks & Gluons Quark & Gluon Fragmentation Functions Q2 dependence predicted from QCD FFF2 1978 Feynman quote from FFF2 “We investigate whether the present experimental behavior of mesons with large transverse momentum in hadron-hadron collisions is consistent with the theory of quantum-chromodynamics (QCD) with asymptotic freedom, at least as the theory is now partially understood.” Parton Distribution Functions Q2 dependence predicted from QCD Quark & Gluon Cross-Sections Calculated from QCD Rick Field – Florida/CDF/CMS
(bk) (ka) (ba) (cb) cc pair bb pair A Parameterization of the Properties of Jets • Assumed that jets could be analyzed on a “recursive” principle. Field-Feynman 1978 Secondary Mesons (after decay) • Let f(h)dh be the probability that the rank 1 meson leaves fractional momentum h to the remaining cascade, leaving quark “b” with momentum P1 = h1P0. Rank 2 Rank 1 • Assume that the mesons originating from quark “b” are distributed in presisely the same way as the mesons which came from quark a (i.e. same function f(h)), leaving quark “c” with momentum P2 = h2P1 = h2h1P0. Primary Mesons continue • Add in flavor dependence by letting bu = probabliity of producing u-ubar pair, bd = probability of producing d-dbar pair, etc. Calculate F(z) from f(h) and bi! • Let F(z)dz be the probability of finding a meson (independent of rank) with fractional mementum z of the original quark “a” within the jet. Original quark with flavor “a” and momentum P0 Rick Field – Florida/CDF/CMS
Feynman-Field Jet Model R. P. Feynman ISMD, Kaysersberg, France, June 12, 1977 Feynman quote from FF2 “The predictions of the model are reasonable enough physically that we expect it may be close enough to reality to be useful in designing future experiments and to serve as a reasonable approximation to compare to data. We do not think of the model as a sound physical theory, ....” Rick Field – Florida/CDF/CMS
High PT Jets CDF (2006) Feynman, Field, & Fox (1978) Predict large “jet” cross-section 30 GeV/c! Feynman quote from FFF “At the time of this writing, there is still no sharp quantitative test of QCD. An important test will come in connection with the phenomena of high PT discussed here.” 600 GeV/c Jets! Rick Field – Florida/CDF/CMS
CDF DiJet Event: M(jj) ≈ 1.4 TeV ETjet1 = 666 GeV ETjet2 = 633 GeV Esum = 1,299 GeV M(jj) = 1,364 GeV M(jj)/Ecm≈ 70%!! Rick Field – Florida/CDF/CMS
Monte-Carlo Simulationof Hadron-Hadron Collisions FF1-FFF1 (1977) “Black-Box” Model FF2 (1978) Monte-Carlo simulation of “jets” F1-FFF2 (1978) QCD Approach FFFW “FieldJet” (1980) QCD “leading-log order” simulation of hadron-hadron collisions “FF” or “FW” Fragmentation my early days yesterday ISAJET (“FF” Fragmentation) HERWIG (“FW” Fragmentation) PYTHIA (“String” Fragmentation) today HERWIG++ PYTHIA 6.4 SHERPA Rick Field – Florida/CDF/CMS
The Fermilab Tevatron CDF “SciCo” Shift December 12-19, 2008 • I joined CDF in January 1998. My wife Jimmie on shift with me! Acquired 4728 nb-1 during 8 hour “owl” shift! Rick Field – Florida/CDF/CMS
Proton-AntiProton Collisionsat the Tevatron The CDF “Min-Bias” trigger picks up most of the “hard core” cross-section plus a small amount of single & double diffraction. stot = sEL + sIN stot = sEL + sSD+sDD+sHC 1.8 TeV: 78mb = 18mb + 9mb + (4-7)mb + (47-44)mb CDF “Min-Bias” trigger 1 charged particle in forward BBC AND 1 charged particle in backward BBC The “hard core” component contains both “hard” and “soft” collisions. “Inelastic Non-Diffractive Component” Beam-Beam Counters 3.2 < |h| < 5.9 Rick Field – Florida/CDF/CMS
“Hard Scattering” Component QCD Monte-Carlo Models:High Transverse Momentum Jets • Start with the perturbative 2-to-2 (or sometimes 2-to-3) parton-parton scattering and add initial and final-state gluon radiation (in the leading log approximation or modified leading log approximation). “Underlying Event” • The “underlying event” consists of the “beam-beam remnants” and from particles arising from soft or semi-soft multiple parton interactions (MPI). The “underlying event” is an unavoidable background to most collider observables and having good understand of it leads to more precise collider measurements! • Of course the outgoing colored partons fragment into hadron “jet” and inevitably “underlying event” observables receive contributions from initial and final-state radiation. Rick Field – Florida/CDF/CMS
3 charged particles dNchg/dhdf = 3/4p = 0.24 1 charged particle Divide by 4p 1 GeV/c PTsum dNchg/dhdf = 1/4p = 0.08 3 GeV/c PTsum dPTsum/dhdf = 1/4p GeV/c = 0.08 GeV/c dPTsum/dhdf = 3/4p GeV/c = 0.24 GeV/c Particle Densities • Study the charged particles (pT > 0.5 GeV/c, |h| < 1) and form the charged particle density, dNchg/dhdf, and the charged scalar pT sum density, dPTsum/dhdf. Charged Particles pT > 0.5 GeV/c |h| < 1 CDF Run 2 “Min-Bias” DhDf = 4p = 12.6 Rick Field – Florida/CDF/CMS
CDF Run 1: Evolution of Charged Jets“Underlying Event” • Look at charged particle correlations in the azimuthal angle Df relative to the leading charged particle jet. • Define |Df| < 60o as “Toward”, 60o < |Df| < 120o as “Transverse”, and |Df| > 120o as “Away”. • All three regions have the same size in h-f space, DhxDf = 2x120o = 4p/3. Charged Particle Df Correlations PT > 0.5 GeV/c |h| < 1 Look at the charged particle density in the “transverse” region! “Transverse” region very sensitive to the “underlying event”! CDF Run 1 Analysis Rick Field – Florida/CDF/CMS
Factor of 2! Run 1 Charged Particle Density“Transverse” pT Distribution • Compares the average “transverse” charge particle density with the average “Min-Bias” charge particle density (|h|<1, pT>0.5 GeV). Shows how the “transverse” charge particle density and the Min-Bias charge particle density is distributed in pT. PT(charged jet#1) > 30 GeV/c “Transverse” <dNchg/dhdf> = 0.56 “Min-Bias” CDF Run 1 Min-Bias data <dNchg/dhdf> = 0.25 Rick Field – Florida/CDF/CMS
MPI: Multiple PartonInteractions • PYTHIA models the “soft” component of the underlying event with color string fragmentation, but in addition includes a contribution arising from multiple parton interactions (MPI) in which one interaction is hard and the other is “semi-hard”. • The probability that a hard scattering events also contains a semi-hard multiple parton interaction can be varied but adjusting the cut-off for the MPI. • One can also adjust whether the probability of a MPI depends on the PT of the hard scattering, PT(hard) (constant cross section or varying with impact parameter). • One can adjust the color connections and flavor of the MPI (singlet or nearest neighbor, q-qbar or glue-glue). • Also, one can adjust how the probability of a MPI depends on PT(hard) (single or double Gaussian matter distribution). Rick Field – Florida/CDF/CMS
Tuning PYTHIA:Multiple Parton Interaction Parameters Hard Core Determine by comparing with 630 GeV data! Affects the amount of initial-state radiation! Take E0 = 1.8 TeV Reference point at 1.8 TeV Rick Field – Florida/CDF/CMS
PYTHIA 6.206 Defaults MPI constant probability scattering • Plot shows the “Transverse” charged particle density versus PT(chgjet#1) compared to the QCD hard scattering predictions of PYTHIA 6.206 (PT(hard) > 0) using the default parameters for multiple parton interactions and CTEQ3L, CTEQ4L, and CTEQ5L. PYTHIA default parameters Default parameters give very poor description of the “underlying event”! Note Change PARP(67) = 4.0 (< 6.138) PARP(67) = 1.0 (> 6.138) Rick Field – Florida/CDF/CMS
Run 1 PYTHIA Tune A CDF Default! • Plot shows the “transverse” charged particle density versus PT(chgjet#1) compared to the QCD hard scattering predictions of two tuned versions of PYTHIA 6.206 (CTEQ5L, Set B (PARP(67)=1)andSet A(PARP(67)=4)). PYTHIA 6.206 CTEQ5L Run 1 Analysis Old PYTHIA default (more initial-state radiation) Old PYTHIA default (more initial-state radiation) New PYTHIA default (less initial-state radiation) New PYTHIA default (less initial-state radiation) Rick Field – Florida/CDF/CMS
PYTHIA Tune A Min-Bias“Soft” + ”Hard” Tuned to fit the CDF Run 1 “underlying event”! PYTHIA Tune A CDF Run 2 Default 12% of “Min-Bias” events have PT(hard) > 5 GeV/c! 1% of “Min-Bias” events have PT(hard) > 10 GeV/c! • PYTHIA regulates the perturbative 2-to-2 parton-parton cross sections with cut-off parameters which allows one to run with PT(hard) > 0. One can simulate both “hard” and “soft” collisions in one program. Lots of “hard” scattering in “Min-Bias” at the Tevatron! • The relative amount of “hard” versus “soft” depends on the cut-off and can be tuned. • This PYTHIA fit predicts that 12% of all “Min-Bias” events are a result of a hard 2-to-2 parton-parton scattering with PT(hard) > 5 GeV/c (1% with PT(hard) > 10 GeV/c)! Rick Field – Florida/CDF/CMS
PYTHIA Tune ALHC Min-Bias Predictions 12% of “Min-Bias” events have PT(hard) > 10 GeV/c! LHC? • Shows the center-of-mass energy dependence of the charged particle density, dNchg/dhdfdPT, for “Min-Bias” collisions compared with PYTHIA Tune A with PT(hard) > 0. 1% of “Min-Bias” events have PT(hard) > 10 GeV/c! • PYTHIA Tune A predicts that 1% of all “Min-Bias” events at 1.8 TeV are a result of a hard 2-to-2 parton-parton scattering with PT(hard) > 10 GeV/c which increases to 12% at 14 TeV! Rick Field – Florida/CDF/CMS
CDF Run 1 PT(Z) Tune used by the CDF-EWK group! PYTHIA 6.2 CTEQ5L • Shows the Run 1 Z-boson pT distribution (<pT(Z)> ≈ 11.5 GeV/c) compared with PYTHIA Tune A (<pT(Z)> = 9.7 GeV/c), and PYTHIA Tune AW (<pT(Z)> = 11.7 GeV/c). UE Parameters ISR Parameters Effective Q cut-off, below which space-like showers are not evolved. Intrensic KT The Q2 = kT2 in as for space-like showers is scaled by PARP(64)! Rick Field – Florida/CDF/CMS
Df Jet#1-Jet#2 Jet#1-Jet#2 Df Distribution Jet-Jet Correlations (DØ) • MidPoint Cone Algorithm (R = 0.7, fmerge = 0.5) • L= 150 pb-1 (Phys. Rev. Lett. 94 221801 (2005)) • Data/NLO agreement good. Data/HERWIG agreement good. • Data/PYTHIA agreement good provided PARP(67) = 1.0→4.0 (i.e. like Tune A, best fit 2.5). Rick Field – Florida/CDF/CMS
CDF Run 1 PT(Z) PYTHIA 6.2 CTEQ5L • Shows the Run 1 Z-boson pT distribution (<pT(Z)> ≈ 11.5 GeV/c) compared with PYTHIA Tune DW, and HERWIG. UE Parameters ISR Parameters Tune DW uses D0’s perfered value of PARP(67)! Intrensic KT Tune DW has a lower value of PARP(67) and slightly more MPI! Rick Field – Florida/CDF/CMS
PYTHIA 6.2 Tunes All use LO as with L = 192 MeV! UE Parameters Uses CTEQ6L Tune A energy dependence! ISR Parameter Intrinsic KT Rick Field – Florida/CDF/CMS
PYTHIA 6.2 Tunes These are “old” PYTHIA 6.2 tunes! There are new 6.420 tunes by Peter Skands (Tune S320, update of S0) Peter Skands (Tune N324, N0CR) Hendrik Hoeth (Tune P329, “Professor”) All use LO as with L = 192 MeV! UE Parameters Tune B Tune AW Tune BW Tune A ATLAS energy dependence! ISR Parameter Tune DW Tune D6 Tune D Tune D6T Intrinsic KT Rick Field – Florida/CDF/CMS
JIMMY at CDF JIMMY was tuned to fit the energy density in the “transverse” region for “leading jet” events! JIMMY Runs with HERWIG and adds multiple parton interactions! PT(JIM)= 2.5 GeV/c. The Drell-Yan JIMMY Tune PTJIM = 3.6 GeV/c, JMRAD(73) = 1.8 JMRAD(91) = 1.8 The Energy in the “Underlying Event” in High PT Jet Production JIMMY: MPI J. M. Butterworth J. R. Forshaw M. H. Seymour PT(JIM)= 3.25 GeV/c. “Transverse” <Densities> vs PT(jet#1) Rick Field – Florida/CDF/CMS
Min-Bias “Associated”Charged Particle Density “Associated” densities do not include PTmax! Highest pT charged particle! • Use the maximum pT charged particle in the event, PTmax, to define a direction and look at the the “associated” density, dNchg/dhdf, in “min-bias” collisions (pT > 0.5 GeV/c, |h| < 1). It is more probable to find a particle accompanying PTmax than it is to find a particle in the central region! • Shows the data on the Df dependence of the “associated” charged particle density, dNchg/dhdf, for charged particles (pT > 0.5 GeV/c, |h| < 1, not including PTmax) relative to PTmax (rotated to 180o) for “min-bias” events. Also shown is the average charged particle density, dNchg/dhdf, for “min-bias” events. Rick Field – Florida/CDF/CMS
Min-Bias “Associated”Charged Particle Density Rapid rise in the particle density in the “transverse” region as PTmax increases! PTmax > 2.0 GeV/c Transverse Region Transverse Region Ave Min-Bias 0.25 per unit h-f PTmax > 0.5 GeV/c • Shows the data on the Df dependence of the “associated” charged particle density, dNchg/dhdf, for charged particles (pT > 0.5 GeV/c, |h| < 1, not including PTmax) relative to PTmax (rotated to 180o) for “min-bias” eventswith PTmax > 0.5, 1.0, and 2.0 GeV/c. • Shows “jet structure” in “min-bias” collisions (i.e.the “birth” of the leading two jets!). Rick Field – Florida/CDF/CMS
Min-Bias “Associated”Charged Particle Density PY Tune A PTmax > 2.0 GeV/c Transverse Region Transverse Region PTmax > 0.5 GeV/c • Shows the data on the Df dependence of the “associated” charged particle density, dNchg/dhdf, for charged particles (pT > 0.5 GeV/c, |h| < 1, not including PTmax) relative to PTmax (rotated to 180o) for “min-bias” eventswith PTmax > 0.5 GeV/c and PTmax > 2.0 GeV/c compared with PYTHIA Tune A (after CDFSIM). • PYTHIA Tune A predicts a larger correlation than is seen in the “min-bias” data (i.e.Tune A “min-bias” is a bit too “jetty”). Rick Field – Florida/CDF/CMS
Min-Bias “Associated”Charged Particle Density • Shows the “associated” charged particle density in the “toward”, “away” and “transverse” regions as a function of PTmax for charged particles (pT > 0.5 GeV/c, |h| < 1, not including PTmax) for “min-bias” events at 1.96 TeVfrom PYTHIA Tune A (generator level). “Toward” Region ~ factor of 2! “Transverse” “Transverse” • Shows the Df dependence of the “associated” charged particle density, dNchg/dhdf, for charged particles (pT > 0.5 GeV/c, |h| < 1, not including PTmax) relative to PTmax (rotated to 180o) for “min-bias” events at 1.96 TeVwith PTmax > 0.5, 1.0, 2.0, 5.0, and 10.0 GeV/c from PYTHIA Tune A (generator level). Rick Field – Florida/CDF/CMS
1st Workshop on Energy Scalingin Hadron-Hadron Collisions • Rick Field Talk 2 Tomorrow at 1:30pm From Min-Bias to the Underlying Event Comparing with the 630 GeV data CDF Run 2 Underlying Event Studies • Rick Field Talk 3 Wednesday at 9:00am From CDF to CMS Tune S320 and P329 compared with Tune A, DW, and DWT Extrapolating to the LHC Rick Field – Florida/CDF/CMS