Below you will find the runlist I used for all of the studies leading up to a preliminary result.  For a more detailed look at how I arrived at this runlist please see my run

I'm concerned about an apparent deficit of g-g scattering events in our CDF Tune A Pythia samples.  I first noticed this deficit when I was looking at MC asymmetries for pi-minus production, where the difference between asymmetries from g-g and q-g events is extreme.  For reference, here is Vogelsang's prediction for the subprocess mixture for inclusive pi0 production:

Compare that to the following plot from Pythia for pi-plus production after combining partonic samples starting at 3_4:

If I add 2_3 GeV the gg contribution at low pT gets a boost, but I'm still seeing ~zero pions from g-g scattering at 8 GeV, where Vogelsang predicts 20% g-g. 

Renee suggested skipping the observables and just plotting the event partonic pT.  The turnover from gg to qg occurs at 3 GeV ... seems too low to me.  On the other hand, the mixture is ~20% g-g at 16 GeV, which seems about right.

One other tidbit for any Pythia tuners out there.  When I was running standalone Pythia to get the various partonic xsections for weighting purposes I saw a number of advisory warnings like

in this case, I was running the 5_7 sample.  The initial max value for the q-g xsec was 5.3392D-01, and the final measured xsec was 1.515D-01.  Not sure if it's significant.

Update 11-19-2007

I looked into the relationship between charged pion p_{T} and event partonic p_{T} in more detail.  Here are plot of charged pion multiplicity and mean pion p_{T} versus event p_{T}:

so g-g events actually produce more charged pions per event, but these pions are all at low p_{T}.  The second conclusion is reinforced by the PDF I've linked at the bottom of the page, which shows pion p_{T} spectra split by subprocess for a range of partonic p_{T} bins.  The g-g events have a much steeper slope.  Those plots are not normalized by the # of events per subprocess, so the conclusion about the number of pions per event is not immediately evident.

Update 11-27-2007

I tried manually rescaling the g-g scattering contribution to the PYTHIA asymmetries; it looks like increasing the g-g by a factor of 5 does a decent job of reproducing the NLO theoretical predictions:

no rescaling

multiply g-g by 3.0

multiply g-g by 5.0

References

Agreement Between "Polarized" Pythia and NLO pQCD
Updated Charged Pion A_{LL} predictions

Outdated or obsolete studies are archived here

Click on the link to download my slides for the PWG meeting.

Presentation

These histograms plot the fraction of reconstructed charged pions in each pion pT bin arising from gg, qg, and qq scattering.  I use the following cuts:

• fabs(mcVertexZ()) < 60
• fabs(vertexZ()) < 60
• geantId() == 8 or 9 (charged pions)
• fabs(etaPr()) < 1
• fabs(dcaGl()) < 1
• fitPts() > 25

I analyzed Pythia samples from the P05ih production in partonic pT bins 3_4 through 55_65 (excluded minbias and 2_3).  The samples were weighted according to partonic x-sections and numbers of events seen and then combined.  StEmcTriggerMaker provided simulations of the HT1 (96201), HT2 (96211), JP1 (96221), and JP2 (96233) triggers.  Here are the results.  The solid lines are MB and are identical in each plot, while the dashed lines are the yields from events passing a particular software trigger.  Each image is linked to a full-resolution copy:

Conclusions

1. Imposing an EMC trigger suppresses gg events and enhances qq, particularly for transverse momenta < 6 GeV/c.  The effect on qg events changes with pT.  The explanation is that the ratio (pion pT / partonic pT) is lower for EMC triggered events than for minimum bias.
2. High threshold triggers change the subprocess composition more than low-threshold triggers.
3. JP1 is the least-biased trigger according to this metric.  There aren't many JP1 triggers in the real data, though, as it was typically prescaled by ~30 during the 2005 pp run.  Most of the stats in the real data are in JP2.

Inclusive Jet Cross Sections in Proton-proton collisions at √s = 200 GeV

Tai Sakuma

Jet Pion Study

Procedure

pi^{-}

pi^{+}