code: def _get_freq(freqfunc, m1, m2, s1z, s2z):
      
      # Convert to SI units for lalsimulation
      m1kg = float(m1) * lal.MSUN_SI
      m2kg = float(m2) * lal.MSUN_SI
      return lalsimulation.SimInspiralGetFrequency(
          m1kg, m2kg, 0, 0, float(s1z), 0, 0, float(s2z), int(freqfunc))
docs: Wrapper of the LALSimulation function returning the frequency
      for a given frequency function and template parameters.
  
      Parameters
      ----------
      freqfunc : lalsimulation FrequencyFunction wrapped object e.g.
          lalsimulation.fEOBNRv2RD
      m1 : float-ish, i.e. castable to float
          First component mass in solar masses
      m2 : float-ish
          Second component mass in solar masses
      s1z : float-ish
          First component dimensionless spin S_1/m_1^2 projected onto L
      s2z : float-ish
          Second component dimensionless spin S_2/m_2^2 projected onto L
  
      Returns
      -------
      f : float
          Frequency in Hz