1 |
#include <OrbitalInfo.h> |
#include <OrbitalInfo.h> |
2 |
|
|
3 |
OrbitalInfo::OrbitalInfo(){ |
using namespace std; |
|
absTime = 0; |
|
|
OBT = 0; |
|
|
pkt_num = 0; |
|
|
|
|
|
lon = 0.; |
|
|
lat = 0.; |
|
|
alt = 0.; |
|
|
|
|
|
Bnorth = 0.; |
|
|
Beast = 0.; |
|
|
Bdown = 0.; |
|
|
Babs = 0.; |
|
|
BB0 = 0.; |
|
|
L = 0.; |
|
|
|
|
|
londip = -1000.; |
|
|
latdip = -1000.; |
|
|
altdip = -1000.; |
|
|
|
|
|
loncgm = -1000.; |
|
|
latcgm = -1000.; |
|
|
altcgm = -1000.; |
|
|
|
|
|
loncbm = -1000.; |
|
|
latcbm = -1000.; |
|
|
altcbm = -1000.; |
|
4 |
|
|
5 |
std::fill_n(cutoff, 20, 0.); |
ClassImp(OrbitalInfoTrkVar); |
6 |
|
ClassImp(OrbitalInfo); |
7 |
|
|
8 |
// Quaternions |
OrbitalInfoTrkVar::OrbitalInfoTrkVar(){ |
9 |
q0 = -1000.; |
this->Clear(); |
10 |
q1 = -1000.; |
} |
|
q2 = -1000.; |
|
|
q3 = -1000.; |
|
11 |
|
|
12 |
// Euler angles (nadir reference frame) |
void OrbitalInfoTrkVar::Clear(Option_t *t){ |
13 |
theta = -1000.; |
trkseqno = 0; |
14 |
phi = -1000.; |
pitch = -1000.; |
15 |
etha = -1000.; |
} |
16 |
|
|
|
// Euler angles (local field reference frame) |
|
|
thetamag = -1000.; |
|
|
phimag = -1000.; |
|
|
ethamag = -1000.; |
|
17 |
|
|
18 |
std::fill_n(goodAttitude, 5, 0); |
OrbitalInfo::OrbitalInfo(){ |
19 |
|
OrbitalInfoTrk = 0; //ELENA |
20 |
|
this->Clear(); |
21 |
} |
} |
22 |
|
|
23 |
|
void OrbitalInfo::Delete(Option_t *t){ //ELENA |
24 |
|
// |
25 |
|
if(OrbitalInfoTrk){ |
26 |
|
OrbitalInfoTrk->Delete(); //ELENA |
27 |
|
delete OrbitalInfoTrk; //ELENA |
28 |
|
} |
29 |
|
// |
30 |
|
}; //ELENA |
31 |
|
|
32 |
|
void OrbitalInfo::Set(){//ELENA |
33 |
|
if ( !OrbitalInfoTrk ) OrbitalInfoTrk = new TClonesArray("OrbitalInfoTrkVar",2); //ELENA |
34 |
|
}//ELENA |
35 |
|
|
36 |
|
OrbitalInfoTrkVar *OrbitalInfo::GetOrbitalInfoTrkVar(Int_t itrk){ |
37 |
|
// |
38 |
|
if(itrk >= ntrk()){ |
39 |
|
printf(" OrbitalInfo ERROR: track related variables set %i does not exists! \n",itrk); |
40 |
|
printf(" stored track related variables = %i \n",ntrk()); |
41 |
|
return(NULL); |
42 |
|
} |
43 |
|
// |
44 |
|
if(!OrbitalInfoTrk)return 0; //ELENA |
45 |
|
TClonesArray &t = *(OrbitalInfoTrk); |
46 |
|
OrbitalInfoTrkVar *orbtrack = (OrbitalInfoTrkVar*)t[itrk]; |
47 |
|
return orbtrack; |
48 |
|
} |
49 |
|
|
50 |
void OrbitalInfo::Clear(){ |
void OrbitalInfo::Clear(Option_t *t){ |
51 |
|
// |
52 |
|
if ( OrbitalInfoTrk ) OrbitalInfoTrk->Delete(); |
53 |
|
// |
54 |
absTime = 0; |
absTime = 0; |
55 |
OBT = 0; |
OBT = 0; |
56 |
pkt_num = 0; |
pkt_num = 0; |
92 |
etha = -1000.; |
etha = -1000.; |
93 |
|
|
94 |
// Euler angles (local field reference frame) |
// Euler angles (local field reference frame) |
95 |
thetamag = -1000.; |
// thetamag = -1000.; |
96 |
phimag = -1000.; |
// phimag = -1000.; |
97 |
ethamag = -1000.; |
// ethamag = -1000.; |
98 |
|
|
99 |
|
pamzenitangle = -1000.; |
100 |
|
pamBangle = -1000.; |
101 |
|
|
102 |
|
|
103 |
std::fill_n(goodAttitude, 5, 0); |
mode = 0; |
104 |
|
// std::fill_n(goodAttitude, 5, 0); |
105 |
} |
} |
106 |
|
|
107 |
/** |
/** |
148 |
l2->phi = phi; |
l2->phi = phi; |
149 |
l2->etha = etha; |
l2->etha = etha; |
150 |
|
|
151 |
|
l2->mode = mode; |
152 |
// Euler angles (local field reference frame) |
// Euler angles (local field reference frame) |
153 |
l2->thetamag = thetamag; |
// l2->thetamag = thetamag; |
154 |
l2->phimag = phimag; |
// l2->phimag = phimag; |
155 |
l2->ethamag = ethamag; |
// l2->ethamag = ethamag; |
156 |
|
|
157 |
memcpy(l2->goodAttitude, goodAttitude, sizeof(goodAttitude)); |
// memcpy(l2->goodAttitude, goodAttitude, sizeof(goodAttitude)); |
158 |
} |
} |
159 |
|
|
160 |
void OrbitalInfo::SetFromLevel2Struct(cOrbitalInfo *l2){ |
void OrbitalInfo::SetFromLevel2Struct(cOrbitalInfo *l2){ |
198 |
phi = l2->phi; |
phi = l2->phi; |
199 |
etha = l2->etha; |
etha = l2->etha; |
200 |
|
|
201 |
|
mode = l2->mode; |
202 |
// Euler angles (local field reference frame) |
// Euler angles (local field reference frame) |
203 |
thetamag = l2->thetamag; |
// thetamag = l2->thetamag; |
204 |
phimag = l2->phimag; |
// phimag = l2->phimag; |
205 |
ethamag = l2->ethamag; |
// ethamag = l2->ethamag; |
206 |
|
|
207 |
memcpy(goodAttitude, l2->goodAttitude, sizeof(l2->goodAttitude)); |
// memcpy(goodAttitude, l2->goodAttitude, sizeof(l2->goodAttitude)); |
208 |
} |
} |
209 |
|
|
|
ClassImp(OrbitalInfo) |
|