code: def get_soap_structure(obj, alp, bet, rCut=5.0, nMax=5, Lmax=5, crossOver=True, all_atomtypes=None, eta=1.0):
      
      Hpos = obj.get_positions()
      arrsoap = get_soap_locals(obj, Hpos, alp, bet,  rCut, nMax, Lmax, crossOver, all_atomtypes=all_atomtypes, eta=eta)
  
      return arrsoap
docs: 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.