Available standard analysis which has been embedded in rivet:
https://rivet.hepforge.org/analyses/STAR_2006_S6500200.html
pT range: 1<pT<10 GeV |eta|<0.5
arXiv: nucl-ex/0601033
March 04 update: modified standard analysis, remove the BBC trigger, 
BBC trigger:have tracks in both 3.3< η <5 and -5< η <-3.3 

Available standard analysis on low pT: https://rivet.hepforge.org/analyses/STAR_2008_S7869363.html
But in this analysis, the final spectra is manually scaled with a scale factor.

My analysis on high pT, reference paper: https://arxiv.org/pdf/1110.0579.pdf
pT range: 3<pT<15 GeV |eta|<0.5
Pythia set up: PYTHIA8.3 default tune, HardQCD:all=on
To enhance the statistics at high pT, the pythia output files are split in mPtHat bins:(1,2,3,5,10,20,50,100), and run rivet analysis separately for each bin. Finally use rivet-merge to add the output Yoda files together.
https://drupal.star.bnl.gov/STAR/system/files/pythiatune0223.pdf