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.