2 |
#define OrbitalInfo_h |
#define OrbitalInfo_h |
3 |
|
|
4 |
#include <TObject.h> |
#include <TObject.h> |
5 |
|
#include <algorithm> // EMILIANO |
6 |
|
#include <OrbitalInfoStruct.h> |
7 |
|
#include <TClonesArray.h> |
8 |
|
#include <TMatrixD.h> |
9 |
|
#include <TVector3.h> |
10 |
|
|
11 |
|
class OrbitalInfoTrkVar : public TObject { |
12 |
|
private: |
13 |
|
|
14 |
|
public: |
15 |
|
// |
16 |
|
Int_t trkseqno; // tof sequ. number: -1=ToF standalone, 0=first Tracker track, ... |
17 |
|
// |
18 |
|
Float_t pitch; ///< Pitch angle |
19 |
|
Float_t sunangle; |
20 |
|
Float_t sunmagangle; |
21 |
|
// |
22 |
|
TMatrixD Eij; ///< vector of incoming particle respect to cartesian geographic coordinates |
23 |
|
TMatrixD Sij; ///< vector of incoming particle respect to flight coordinates |
24 |
|
// |
25 |
|
Float_t cutoff; ///< Calculated cutoff for the incoming particle taking into account particle direction |
26 |
|
// |
27 |
|
OrbitalInfoTrkVar(); |
28 |
|
OrbitalInfoTrkVar(const OrbitalInfoTrkVar&); |
29 |
|
OrbitalInfoTrkVar* GetOrbitalInfoTrkVar(){return this;}; |
30 |
|
// |
31 |
|
void Clear(Option_t *t=""); |
32 |
|
void Delete(Option_t *t=""); //ELENA |
33 |
|
// |
34 |
|
ClassDef(OrbitalInfoTrkVar, 3); |
35 |
|
// |
36 |
|
}; |
37 |
|
|
38 |
|
|
39 |
|
|
40 |
|
/** |
41 |
|
* Class that stores position, time, inclination, magnetic field and |
42 |
|
* cutoff informations. |
43 |
|
*/ |
44 |
class OrbitalInfo : public TObject { |
class OrbitalInfo : public TObject { |
45 |
|
private: |
46 |
|
|
47 |
public: |
public: |
48 |
|
TClonesArray *OrbitalInfoTrk; |
49 |
|
|
50 |
|
UInt_t absTime; ///< Absolute Time (seconds) |
51 |
|
UInt_t OBT; ///< On Board Time (ms) |
52 |
|
UInt_t pkt_num; ///< CPU packet number |
53 |
|
|
54 |
|
Float_t lon; ///< degrees from -180 to 180 |
55 |
|
Float_t lat; ///< degrees from -90 to 90 |
56 |
|
Float_t alt; ///< meters asl |
57 |
|
|
58 |
|
TVector3 V; /// velocity |
59 |
|
|
60 |
|
// B components. |
61 |
|
Float_t Bnorth; ///< gauss |
62 |
|
Float_t Beast; ///< gauss |
63 |
|
Float_t Bdown; ///< gauss |
64 |
|
|
65 |
|
Float_t Babs; ///< abs value (guass) |
66 |
|
|
67 |
|
Float_t M; ///< M |
68 |
|
|
69 |
|
Float_t BB0; ///< B abs over the B minimum on this field line |
70 |
|
|
71 |
|
Float_t L; ///< McIlwain's L shell (in earth radii) |
72 |
|
|
73 |
|
/* // Dipolar magnetic coordinates (not used). */ |
74 |
|
Float_t londip; ///< degrees from -180 to 180 |
75 |
|
Float_t latdip; ///< degrees from -90 to 90 |
76 |
|
|
77 |
|
// Float_t cutoff[17]; |
78 |
|
Float_t cutoffsvl; |
79 |
|
Float_t igrf_icode; |
80 |
|
|
81 |
|
// linear Quaternions |
82 |
|
Float_t q0; ///< Quaternion 0 |
83 |
|
Float_t q1; ///< Quaternion 1 |
84 |
|
Float_t q2; ///< Quaternion 2 |
85 |
|
Float_t q3; ///< Quaternion 3 |
86 |
|
|
87 |
|
// Euler angles (Resurs velocity reference frame) |
88 |
|
Float_t theta; ///< Euler angle theta in the velocity reference frame (pitch) |
89 |
|
Float_t phi; ///< Euler angle phi in the velocity reference frame (yaw) |
90 |
|
Float_t etha; ///< Euler angle etha in the velocity reference frame (roll) |
91 |
|
|
92 |
|
// Pitch angles |
93 |
|
|
94 |
|
// |
95 |
|
TMatrixD Iij; ///< Angle between PAMELA Z direction and cartesian geographic coordinates |
96 |
|
|
97 |
|
/** |
98 |
|
* The variable mode means a character time distant between two quaternions, inside which stay every events |
99 |
|
*/ |
100 |
|
Int_t mode; // 0 - means that time different pair of quaternions exuals to 0.25 seconds in R10 mode |
101 |
|
// (it mean that all quaternions in array is correct) |
102 |
|
// 1 - means that we have R10 mode and use just first value of quaternions array |
103 |
|
// 2 - means that we have non R10 mode and use every quaternions from array. |
104 |
|
// 3 - means normal transition from R10 to non R10 or from non R10 to R10. |
105 |
|
// 4 - means that we have eliminable hole between R10 and non R10 or between non R10 and R10 |
106 |
|
// 5 - means that we have uneliminable hole between R10 and non R10 or between non R10 and R10 |
107 |
|
// 6 - means that we have eliminable hole inside R10 (in such keys eliminable depends from other factors also) |
108 |
|
// 7 - means that we have uneliminable hole inside R10 |
109 |
|
// 8 - means that we have eliminable hole inside non R10 |
110 |
|
// 9 - means that we have uneliminable hole inside non R10 |
111 |
|
// 10 - means other unknown problems |
112 |
|
// -10 - means we use recovered quaternions |
113 |
|
|
114 |
|
Int_t qkind; // How matrix Qij was got. |
115 |
|
// 0 means that it was calculated from flight quaternion |
116 |
|
// 1 means that it was calculated from Euler angles from Rotation Table |
117 |
|
|
118 |
|
Float_t TimeGap; //Time gap between two points where interpolation have done. |
119 |
|
|
120 |
|
Int_t errq; // flag, if errq == 1 then real flight quaternion (not interpolated) incorrect |
121 |
|
Int_t azim; // 0 - means everything is ok |
122 |
|
// 1 - means azimutal rotations were performed in this moment and in case of absense of flight quaternions orientation calculated here incorrect |
123 |
|
// >1 No flight quaternions, no azimuthal rotations, no adequate data from RotationTable, unaccuracy equals to Bank angle of rotetion in this moment |
124 |
|
|
125 |
|
Int_t ntrk(){return OrbitalInfoTrk->GetEntries();}; |
126 |
|
/** |
127 |
|
\return Returns the B minimum along the field line. |
128 |
|
*/ |
129 |
|
Float_t GetB0() { return Babs/BB0; }; |
130 |
|
|
131 |
|
/** |
132 |
|
\return Returns the Stormer vertical cutoff using L shell: |
133 |
|
14.9/L^2 (GV/c). |
134 |
|
*/ |
135 |
|
Float_t GetCutoffSVL() { return cutoffsvl; }; |
136 |
|
|
137 |
|
void SetFromLevel2Struct(cOrbitalInfo *l2); |
138 |
|
void GetLevel2Struct(cOrbitalInfo *l2) const; |
139 |
|
OrbitalInfoTrkVar *GetOrbitalInfoTrkVar(Int_t notrack); |
140 |
|
// |
141 |
OrbitalInfo(); |
OrbitalInfo(); |
142 |
|
~OrbitalInfo(){Delete();}; //ELENA |
143 |
|
// |
144 |
|
OrbitalInfo* GetOrbitalInfo(){return this;}; // Elena |
145 |
|
void Delete(Option_t *t=""); //ELENA |
146 |
|
void Set();//ELENA |
147 |
|
// |
148 |
|
// |
149 |
|
|
150 |
ULong64_t absTime; //< Absolute Time |
void Clear(Option_t *t=""); // emiliano |
|
// EM: added On Board Time and CPU Packet Number |
|
|
ULong64_t OBT; //< On Board Time |
|
|
ULong64_t pkt_num; //< CPU packet number |
|
151 |
// |
// |
152 |
ClassDef(OrbitalInfo, 1) |
ClassDef(OrbitalInfo, 10); |
153 |
}; |
}; |
154 |
#endif |
#endif |