def get_soap_structure(obj, alp, bet, rCut=5.0, nMax=5, Lmax=5, crossOver=True, all_atomtypes=None, eta=1.0):
    """Get the RBF basis SOAP output for atoms in a finite structure.

    Args:
        obj(ase.Atoms): Atomic structure for which the SOAP output is
            calculated.
        alp: Alphas
        bet: Betas
        rCut: Radial cutoff.
        nMax: Maximum nmber of radial basis functions
        Lmax: Maximum spherical harmonics degree
        crossOver:
        all_atomtypes: Can be used to specify the atomic elements for which to
            calculate the output. If given the output is calculated only for the
            given species.
        eta: The gaussian smearing width.

    Returns:
        np.ndarray: SOAP output for the given structure.
    """
    Hpos = obj.get_positions()
    arrsoap = get_soap_locals(obj, Hpos, alp, bet,  rCut, nMax, Lmax, crossOver, all_atomtypes=all_atomtypes, eta=eta)

    return arrsoap
Get the RBF basis SOAP output for atoms in a finite structure.

    Args:
        obj(ase.Atoms): Atomic structure for which the SOAP output is
            calculated.
        alp: Alphas
        bet: Betas
        rCut: Radial cutoff.
        nMax: Maximum nmber of radial basis functions
        Lmax: Maximum spherical harmonics degree
        crossOver:
        all_atomtypes: Can be used to specify the atomic elements for which to
            calculate the output. If given the output is calculated only for the
            given species.
        eta: The gaussian smearing width.

    Returns:
        np.ndarray: SOAP output for the given structure.