SIM: fix some wrong input
while calculate momentum,
double new_len = mean_length(h.preMom(),h.postMom())/hit_len;
direction *= new_len/hit_len;
length was not right to divide twice, fixed to
double new_len = mean_length(h.preMom(),h.postMom());
direction *= new_len/hit_len;