Here's a fragmentation study looking at the ratio of reconstructed charged pion p_{T} and the event partonic p_{T} in PYTHIA. Cuts are

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

Error bars are just the errors on the mean partonic p_{T} in each reconstructed pion p_{T} bin. Next step is to look at the jet simulations to come up with a plot that is more directly comparable to a real data measurement.