Working on ScatteringFactors
implemented density conversion function...
added new + and - operator for spectrum
implemented Structure factor calculation, with matrix implementation for speed