code: def get_freq(freqfunc, m1, m2, s1z, s2z):
      
      lalsim_ffunc = getattr(lalsimulation, freqfunc)
      return _vec_get_freq(lalsim_ffunc, m1, m2, s1z, s2z)
docs: Returns the LALSimulation function which evaluates the frequency
      for the given frequency function and template parameters.
  
      Parameters
      ----------
      freqfunc : string
          Name of the frequency function to use, e.g., 'fEOBNRv2RD'
      m1 : float or numpy.array
          First component mass in solar masses
      m2 : float or numpy.array
          Second component mass in solar masses
      s1z : float or numpy.array
          First component dimensionless spin S_1/m_1^2 projected onto L
      s2z : float or numpy.array
          Second component dimensionless spin S_2/m_2^2 projected onto L
  
      Returns
      -------
      f : float or numpy.array
          Frequency in Hz