docs: def calc_phase_space_density(self, mass, omega_array, temperature_array):
          """
          Calculates the partition function of your system at each point in time.
      
          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 the Hamiltonian
          temperature_array : array
              array which represents the temperature at every point in your time trace
              and should therefore have the same length as the Hamiltonian
          
          Requirements
          ------------
          self.position_data : array
              Already filtered for the degree of freedom of intrest and converted into meters. 
  
          Returns:
          -------
          Phasespace-density : array
              The Partition Function at every point in time over a given trap-frequency and temperature change.
          """
  
          return self.calc_hamiltonian(mass, omega_array)/calc_partition_function(mass, omega_array,temperature_array)