merge Jia's pt unintegrated implementations and the implementations for outputting TaTb and rho_binary profiles
add reaction plane average option
speed up a little bit
improved a little bit in the performance
optimize the speed for operation = 3