code: def get_periodic_soap_locals(obj, Hpos, alp, bet, rCut=5.0, nMax=5, Lmax=5, crossOver=True, all_atomtypes=None, eta=1.0):
      
      suce = _get_supercell(obj, rCut)
      arrsoap = get_soap_locals(suce, Hpos, alp, bet, rCut, nMax=nMax, Lmax=Lmax, crossOver=crossOver, all_atomtypes=all_atomtypes, eta=eta)
  
      return arrsoap
docs: Get the RBF basis SOAP output for the given position in a periodic system.
  
      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 position.