-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathnorm_distr_lum_model.cpp
140 lines (111 loc) · 4.82 KB
/
norm_distr_lum_model.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
#include <cmath>
#include <iostream>
#include <cstdlib>
#include "stars.h"
using namespace std;
//-----------------------------------------------------------------------//
// Модель светимости пульсара, при которой энергия в конусе излучения
// распределена по нормальному закону.
//-----------------------------------------------------------------------//
//-----------------------------------------------------------------------//
// Декларация функций
double S_min (double, double, float, double, double, double, float, TMap *);
//-----------------------------------------------------------------------//
double LMExpon::is_pulsar_visible (double t, SpecialStar * sun, TMap * T_copy, double x, double y, double z, double i_incl, double P, double dot_P, float DM) {
double res, time_tmpl;
float l, b, sm;
double dist_to_sun, w50, lum_0, lum_min = 0.3, omega, theta;
double up_border, down_border;
double first[3], second[2];
double chance_1, chance_2, chance_3;
// sun->move_to(t);
dist_to_sun = sqrt(pow(sun->get_position_x() - x, 2) + pow(sun->get_position_y() - y, 2) + pow(sun->get_position_z() - z, 2));
// Вектор от солнца к пульсару направленный
first [0] = sun->get_position_x() - x;
first [1] = sun->get_position_y() - y;
first [2] = sun->get_position_z() - z;
// Вектор от Солнца к центру Галактики
second[0] = - sun->get_position_x();
second[1] = - sun->get_position_y();
b = first[2] / sqrt(pow(first[0], 2) + pow(first[1], 2) + pow(first[2], 2));
b = asin (b)/2./pi*360.;
w50 = 1.9*log10(3.2e19*sqrt(P*dot_P))-21.5 + ds;
w50 = pow(10, w50) / 180 * pi;
// Встроенная модель светимости
time_tmpl = t;
// if (time_tmpl > 1000) {
lum_0 = -0.646*log10(time_tmpl) + 7.524 + dlum;
lum_0 = pow(10, lum_0); // }
// else if (time_tmpl < 1000 && time_tmpl > 0) {
// lum_0 = -0.646*3. + 7.524 + param->get_dlum();
// lum_0 = 1000*param->get_dlum();
// lum_0 = pow(10.,lum_0); }
//-------------------------------//
lum_0 /= pow(dist_to_sun, 2);
// DM = get_DM (t, sun, &l, &b, &sm);
lum_min = S_min (l, b, sm, dist_to_sun, w50, P, DM, T_copy);
if (lum_0 < lum_min) {
return 0.;
} else {
// omega - угол, на который мы можем отклониться от центра консуса,
// чтобы ещё можно было увидеть пульсар.
omega = w50 * log(lum_0 / lum_min);
if (!does_axis_set) {
//-----------------------------------------------------------------
// Генерация оси вращения пульсара и расположения магнитного полюса
chance_1 = rand()/rand_high_board;
chance_1 *= 2*pi;
do {
chance_2 = rand()/rand_high_board;
chance_2 *= pi;
chance_2 -= pi/2.;
chance_3 = rand()/rand_high_board;
} while (cos(chance_2) <= chance_3);
x_axis = cos(chance_2) * cos(chance_1);
y_axis = cos(chance_2) * sin(chance_1);
z_axis = sin(chance_2);
does_axis_set = true;
}
theta = (x_axis * first[0] + y_axis * first[1] + z_axis * first[2])/(dist_to_sun);
if ((omega + i_incl)>pi/2.) {
down_border = pi/2;
} else {
down_border = omega + i_incl;
}
if ((i_incl - omega)<0) {
up_border = 0;
} else {
up_border = i_incl - omega;
}
if (((cos(down_border) <= theta && cos(up_border) >= theta) ||
(cos(pi - down_border) >= theta && cos(pi - up_border) <= theta))
&& abs(b)<15. ) {
res = abs(abs(acos (theta)) - abs(i_incl));
res /= abs(omega);
res = lum_0 * pow(e, -res);
return res;
} else {
return 0.;
}
}
return res;
}
LMExpon::LMExpon (vector <double> * values) {
ds = values->at(1);
dlum = values->at(3);
does_axis_set = false;
}
bool LMExpon::is_beam_on(double P) {
return true;
}
void LMExpon::print_description (ostream * out) {
*out<<"// Normal model of energy distribution inside a cone //"<<endl;
*out<<"//----------------------------------------------------------//"<<endl;
}
void LMExpon::print_parameters (ostream * out) {
*out<<"// Parameters of luminosity model //"<<endl;
*out<<"//----------------------------------------------------------//"<<endl;
*out<<"// ds - "<<ds<<endl;
*out<<"// dlum - "<<dlum<<endl;
*out<<"//----------------------------------------------------------//"<<endl;
}