How to use Bessel beam and the Laguerre-Gaussian beam as the incident light

[WeiHuan2020]

Dear Professor：
I am honored to read your article entitled “MCmatlab: an open-source user-friendly MATLAB-integrated 3D Monte Carlo light transport solver with heat diffusion and tissue damage” and the matlab code in https://github.com/ankrh/MCmatlab/blob/Release/README.md. This work is very great and usefull for many people. But I still have some questions and find for your help.
I want to know how to use how to use Bessel beam and the Laguerre-Gaussian beam as the incident light. For example, in Example 2 in Mcmatlab，could you please show me how to set the light source as Bessel beam and the Laguerre-Gaussian beam?
The second question is I don't know the physical significance of the following variable.
model.MC.lightSource.focalPlaneIntensityDistribution.radialDistr = 0; % Radial focal plane intensity distribution - 0: Top-hat, 1: Gaussian, Array: Custom. Doesn't need to be normalized.
model.MC.lightSource.focalPlaneIntensityDistribution.radialWidth = .01; % [cm] Radial focal plane 1/e^2 radius if top-hat or Gaussian or half-width of the full distribution if custom
model.MC.lightSource.angularIntensityDistribution.radialDistr = 1; % Radial angular intensity distribution - 0: Top-hat, 1: Gaussian, 2: Cosine (Lambertian), Array: Custom. Doesn't need to be normalized.
model.MC.lightSource.angularIntensityDistribution.radialWidth = pi/6; % [rad] Radial angular 1/e^2 half-angle if top-hat or Gaussian or half-angle of the full distribution if custom. For a diffraction limited Gaussian beam, this should be set to model.MC.wavelength1e-9/(pimodel.MC.lightSource.focalPlaneIntensityDistribution.radialWidth*1e-2))
model.MC.lightSource.theta = 0; % [rad] Polar angle of beam center axis
Thank you for taking the time to read my question. If you could help me answer it, I would be extremely grateful.