def _get_singlesite_ops(h_loc, params, bs, b_daggers):
    """
    List of the operators acting on every single site after chain mapping

    Args:
        h_loc (numpy.ndarray): Local Hamiltonian
        params (list): Parameters as returned by :func:`_get_parameters`
        bs (list): The list of annihilation operators acting on each site
            of the chain
        b_daggers (list): The list of creation operators acting on each site
            of the chain

    Returns:
        list: List of operators acting on every single site
    """
    omegas, ts, c0 = params
    singlesite_ops = [omegas[i]
                      * b_daggers[i].dot(bs[i]) for i in range(len(bs))]
    singlesite_ops = [h_loc] + singlesite_ops

    return singlesite_ops