19 def __init__(self, spin_polarization=[0,0,1], spin_superposition=[1 + 0j, 0 + 0j], momentum=[1,0,0], mass=electron_mass):
20
21 """
22 @brief se normaliza el vector (al necesitarse un vector unitario)
23 """
24 norm_n = np.linalg.norm(spin_polarization)
25 self.spin_polarization = np.array(spin_polarization) / norm_n if norm_n > 0 else np.array([0,0,1])
26
27 """
28 @brief se comprueba que la suma del modulo al cuadrado sea 1
29 """
30 norm_s = np.linalg.norm(spin_superposition)
31 self.spin_superposition = np.array(spin_superposition, dtype=np.complex128) / norm_s if norm_s > 0 else np.array([1,0])
32
33 self.momentum = np.array(momentum)
34 self.mass = mass
35
36 self._base()
37
38 self._pauli_spinor()
39
40 self._constant_spinor()
41