Earthquakes clustering based on maximum likelihood estimation of point process conditional intensity function