Crystal.hpp:
в классе BeamSplitting массив Facets перенесен из protected в public (стал доступен для изменений вне класса)

main.cpp:
1. (+) F_Mt (флаг: 1 = записывать в файлы 2Д-картинки для номеров диф. колец и оценки интенсивности, 0 - записывать только гистограммы интенсивности)
2. (+) betta_i, gamma_j - переменные в циклах по простр. ориентациям
3. (+) cos_angle = косинус "границы" окрестности направления назад
4. (+) lm  = длина волны
5. (-) mxd, back, forw, ThetaNumber  - не нужны
6. (+) vector<matrix> nc  - хранит номера колец, dif - оценку интенсивности (для F_Mt = 1)
7. (+) vector<double> en - хранит суммарную оценку интенсивности (для F_Mt = 0)

в main - функции

8. (-) Функцию ShowCurrentTime, вставила строки непосредственно в main-функцию (+ #include <QTime>)
9. (+) для Sorting>0 Facets = Face (в цикле)
10. цикл по ориентациям с равными шагами по bettа и гамма (! изменилась нормировка)
11. вывод в файл осуществляется в папку Data_(время старта расчетов): "dif*.dat", "nc*.dat" - оценка интенсивности и номера колец для траектории "*" (для случая F_Mt = 1); log.dat - в этот файл записываются данные в формате "betta_i%соответствующий угол бетта в градусах%число траектории пучков, полученное к данному моменту"; res.dat - суммарная оценка интенсивности.

в Handler - функции:

12. Добавлен поиск "наперед" заданных траекторий (для Sorting>0)
13. Выбираем из "оставшихся" траектории попавшие в конус обратного рассеяния, заполняем список найденных траекторий
14. Находим номер кольца и оценку интенсивности(с вероятностью P).
в main.cpp:

1. в функции MaskAppend удалила неиспользуемую переменную n
2. в Handler изменила условие "попадания" в окрестность направления назад
   if((bm.r*k)<cos_angle) return; // где cos_angle = cos(M_PI/2.0-params[7]*M_PI/180.0);
3. в Handler  заменила AreaOfBeam  на CrossSection

в PARAMS.DAT:

7-ая строка теперь  - радиус окрестности в направлении назад (<90)
в PARAMS.DAT:

10-ая строка  - длина волны
12-ая строка - флаг: "0" - означает, что считаем только гистограмму оценочной энергии, "1" означает, что рассчитываем также 2Д-картинки (по углам betta и gamma) для оценочной энергии и номера кольца, которым пучок вносит вклад в направление назад
Merge branch 'physical-optics-dev' of github.com:sasha-tvo/Beam-Splitting_GO into physical-optics-dev