We present a new algorithm for the identification and physical characterization of current sheets and reconnection sites as well as the update of post-reconnection particles spectra in 2D and 3D large scale relativistic magnetohydrodynamic simulations. Lagrangian particles, which follow the fluid, are used to sample plasma parameters before entering the reconnection sites that form during the evolution of the different configurations considered. With the sampled parameters and a subgrid model based on results of Particle-in-Cell simulations we introduced in the PLUTO code an algorithm able to describe the post-reconnection spectra associated to the non-thermal component.