This function is to compute the centralized treatment process given the covariates Z, L at each sorted unique event time
trt_center(D_status, Z, weights)