add shell function to collect and retrieve particle's average vn flow
separate average pT values for rapidity cut and pseudorapidity cut spectra