Lib-файлы проекта взяты из Beam-Splitting_physical-optics,
main.cpp и params.dat из "Beam-Splitting_Physical-optics-estimation_with_approximation_diffraction

в main.cpp:

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

в main():

1. (-) Функцию ShowCurrentTime, вставила строки непосредственно в main-функцию (+ #include <QTime>, + #include <QDate>)
2. (+) для Sorting>0 Facets = Face (в цикле)
3. T - угол флаттера (0<T<90 град, T = 0 - случай хаотической ориентации)
4. цикл по ориентациям с равными шагами по bettа и гамма (! изменилась нормировка)
5. (+) p_Betta - массив вероятностей по нормальному закону, где сигма = Т/3 (для T>0, в противном случае p_Betta[i] = 1.0)
	выводятся в файл "probability.dat"
6. теперь P = p_Betta[betta_i]*sin(betta)
7. вывод в файл осуществляется в папку Data_(дата_время старта расчетов): "dif*.dat", "nc*.dat" - оценка интенсивности и номера колец для траектории "*" (для случая F_Mt = 1);	log.dat - в этот файл записываются данные в формате "betta_i%соответствующий угол бетта в градусах%число траектории пучков, полученное к данному моменту"; res.dat - суммарная оценка интенсивности; 	lim*.dat - границы интервалов по углу гама (при фиксированных бетта), где вклад траектории, соответствует 5 или меньшему "кольцу" дифракции

в Handler - функции:

1. Добавлен поиск "наперед" заданных траекторий (для Sorting>0)
2. Выбираем из "оставшихся" траектории попавшие в конус обратного рассеяния, заполняем список найденных траекторий if((bm.r*k)<cos_angle) return; // где cos_angle = cos(M_PI/2.0-params[7]*M_PI/180.0);
3. Находим номер кольца и точное значение интенсивности(с вероятностью P).

в PARAMS.DAT:

1. новый параметр T = params[9] - угол флаттера
2. нумерация с main-файле у параметров, начилая с 9-го сдвинулась на 1.

для теста ID-траектории пучка (уникального номера, для ускорения поиска нужной траектории):
1. в классе Beam добавлена переменная long int id (по умолчанию = 0)
2. в функции трассировок добавлены строки с вычислением ID для отраженных и преломленных пучков
3. в класс Chain - добавлена переменная ID, изменен конструктор класса