text: if __name__ == "__main__":
      # === Imports === #
      import numpy as np
      import LineshapeFuncs as lf
      import matplotlib.pyplot as plt
      import scipy.io as spio
      
      # === Declare varibales === #
      AngleEjct = [0]
      EGamma = 0
      AngleRangeBeta = [0]
      SD = 5
  
      # === statics === #
      Counter = 1
      inv = [-1,1]
  
      SKI = lf.Ski(AngleEjct)
      IntEjct=SKI[0]
      IntIon=SKI[1]
  
      SRIM = lf.Srim(IntEjct,IntIon)
      X= SRIM[0]
      S = SRIM[1]
      
      for Beta in AngleRangeBeta:
          for i in [1,-1]:
              
              Stretch = np.cos(abs(Beta))
              X2 = X*Stretch
              Xtra1 = lf.linspace( X2[0]-8, X2[0]-0.001 , 0.01) 
              Xtra2 = lf.linspace( X2[len(X2)-1]+0.001, X2[len(X2)-1]+8 , 0.01)
              X2 = np.concatenate((Xtra1, X2), axis=0)
              X2 = np.concatenate((X2, Xtra2), axis=0)
              ExtraLength =  len(X2) - len(S)
              ExtraLength=  ExtraLength /  2
              ExtraLength= int(ExtraLength) 
              zer = np.zeros(ExtraLength)
              S = np.concatenate((zer, S), axis=0)
              S = np.concatenate((S, zer), axis=0)
              if i ==-1:
                  S=np.flip(S)
              
              
              StepSize = len(X)/672; 
              X_int = lf.linspace(X2[0],X2[len(X2)-1],0.01 ) 
              S_int =lf. np.interp(X_int,X2,S)   
              gaus = lf.Gaus(X_int,SD)
  
              S_int =     np.convolve(S_int,gaus,'same')
              S_int=S_int/max(S_int);  
              X_int = 1.39*X_int;
              
              
              plt.plot(X_int,S_int)
              plt.show()
  