Previously, I took a stab at unfolding our pi0 and gamma-rich data using only the embedding sample:

https://drupal.star.bnl.gov/STAR/blog/dmawxc/update-08162019-run-9-pp-new-r-02-pi0-unfoldings
https://drupal.star.bnl.gov/STAR/blog/dmawxc/update-08192019-run-9-pp-new-r-02-gamma-unfolding

I realized, though, that I wasn't summing the systematic and statistical errors when I plotted the "final" uncertainty bands against Pythia.  The plots below fix this.  We also decided to add a levy fit to the list of systematic variations.  For the gamma triggers, we also added the gamma-direct purity variation to the systematic band, and decided to stick with the 2nd subtraction scheme (subtract then unfold) from the above post.

Pi0 Trigger:

Gamma Trigger: