docs: def calc_hamiltonian(self, mass, omega_array):
          """
          Calculates the standard (pot+kin) Hamiltonian of your system.
          
          Parameters
          ----------
          mass : float
              The mass of the particle in kg
          omega_array : array
              array which represents omega at every point in your time trace
              and should therefore have the same length as self.position_data
          
          Requirements
          ------------
          self.position_data : array
              Already filtered for the degree of freedom of intrest and converted into meters. 
  
          Returns
          -------
          Hamiltonian : array
              The calculated Hamiltonian
          """
          Kappa_t= mass*omega_array**2
          self.E_pot = 0.5*Kappa_t*self.position_data**2
          self.E_kin = 0.5*mass*(_np.insert(_np.diff(self.position_data), 0, (self.position_data[1]-self.position_data[0]))*self.SampleFreq)**2
          self.Hamiltonian = self.E_pot + self.E_kin
          return self.Hamiltonian