A three-dimensional kinetic-fluid numerical code to study the equilibrium structure of the magnetotail: The role of electrons in the formation of the bifurcated current sheet