| 2 |
#define OrbitalInfo_h |
#define OrbitalInfo_h |
| 3 |
|
|
| 4 |
#include <TObject.h> |
#include <TObject.h> |
| 5 |
|
#include <algorithm> // EMILIANO |
| 6 |
#include <OrbitalInfoStruct.h> |
#include <OrbitalInfoStruct.h> |
| 7 |
|
|
| 8 |
|
/** |
| 9 |
|
* Class that stores position, time, inclination, magnetic field and |
| 10 |
|
* cutoff informations. |
| 11 |
|
*/ |
| 12 |
class OrbitalInfo : public TObject { |
class OrbitalInfo : public TObject { |
| 13 |
public: |
public: |
| 14 |
OrbitalInfo(); |
OrbitalInfo(); |
| 15 |
OrbitalInfo* GetOrbitalInfo(){return this;}; // Elena |
OrbitalInfo* GetOrbitalInfo(){return this;}; // Elena |
| 16 |
|
|
| 17 |
UInt_t absTime; //< Absolute Time |
UInt_t absTime; ///< Absolute Time (seconds) |
| 18 |
// EM: added On Board Time and CPU Packet Number |
UInt_t OBT; ///< On Board Time (ms) |
| 19 |
UInt_t OBT; //< On Board Time |
UInt_t pkt_num; ///< CPU packet number |
| 20 |
UInt_t pkt_num; //< CPU packet number |
|
| 21 |
// |
Float_t lon; ///< degrees from -180 to 180 |
| 22 |
|
Float_t lat; ///< degrees from -90 to 90 |
| 23 |
|
Float_t alt; ///< meters asl |
| 24 |
|
|
| 25 |
|
// B components. |
| 26 |
|
Float_t Bnorth; ///< gauss |
| 27 |
|
Float_t Beast; ///< gauss |
| 28 |
|
Float_t Bdown; ///< gauss |
| 29 |
|
|
| 30 |
|
Float_t Babs; ///< abs value (guass) |
| 31 |
|
|
| 32 |
|
Float_t BB0; ///< B abs over the B minimum on this field line |
| 33 |
|
|
| 34 |
|
Float_t L; ///< McIlwain's L shell (in earth radii) |
| 35 |
|
|
| 36 |
|
// Dipolar magnetic coordinates (not used). |
| 37 |
|
Float_t londip; ///< degrees from -180 to 180 |
| 38 |
|
Float_t latdip; ///< degrees from -90 to 90 |
| 39 |
|
Float_t altdip; ///< meters |
| 40 |
|
|
| 41 |
|
// Corrected magnetic coordinates (not used). |
| 42 |
|
Float_t loncgm; ///< degrees from -180 to 180 |
| 43 |
|
Float_t latcgm; ///< degrees from -90 to 90 |
| 44 |
|
Float_t altcgm; ///< meters |
| 45 |
|
|
| 46 |
|
// Corrected B min magnetic coordinates (not used). |
| 47 |
|
Float_t loncbm; ///< degrees from -180 to 180 |
| 48 |
|
Float_t latcbm; ///< degrees from -90 to 90 |
| 49 |
|
Float_t altcbm; ///< meters |
| 50 |
|
|
| 51 |
|
Float_t cutoff[20]; |
| 52 |
|
|
| 53 |
|
// Quaternions |
| 54 |
|
Float_t q0; ///< Quaternion 0 |
| 55 |
|
Float_t q1; ///< Quaternion 1 |
| 56 |
|
Float_t q2; ///< Quaternion 2 |
| 57 |
|
Float_t q3; ///< Quaternion 3 |
| 58 |
|
|
| 59 |
|
// Euler angles (Resurs velocity reference frame) |
| 60 |
|
Float_t theta; ///< Euler angle theta in the velocity reference frame (pitch) |
| 61 |
|
Float_t phi; ///< Euler angle phi in the velocity reference frame (yaw) |
| 62 |
|
Float_t etha; ///< Euler angle etha in the velocity reference frame (roll) |
| 63 |
|
|
| 64 |
|
/** |
| 65 |
|
* The variable mode means a character time distant between two quaternions, inside which stay every events |
| 66 |
|
*/ |
| 67 |
|
Int_t mode; // 0 - means that time different pair of quaternions exuals to 0.25 seconds in R10 mode |
| 68 |
|
// (it mean that all quaternions in array is correct) |
| 69 |
|
// 1 - means that we have R10 mode and use just first value of quaternions array |
| 70 |
|
// 2 - means that we have non R10 mode and use every quaternions from array. |
| 71 |
|
// 3 - means normal transition from R10 to non R10 or from non R10 to R10. |
| 72 |
|
// 4 - means that we have eliminable hole between R10 and non R10 or between non R10 and R10 |
| 73 |
|
// 5 - means that we have uneliminable hole between R10 and non R10 or between non R10 and R10 |
| 74 |
|
// 6 - means that we have eliminable hole inside R10 (in such keys eliminable depends from other factors also) |
| 75 |
|
// 7 - means that we have uneliminable hole inside R10 |
| 76 |
|
// 8 - means that we have eliminable hole inside non R10 |
| 77 |
|
// 9 - means that we have uneliminable hole inside non R10 |
| 78 |
|
|
| 79 |
|
/** |
| 80 |
|
\return Returns the B minimum along the field line. |
| 81 |
|
*/ |
| 82 |
|
Float_t GetB0() { return Babs/BB0; }; |
| 83 |
|
|
| 84 |
|
/** |
| 85 |
|
\return Returns the Stormer vertical cutoff using L shell: |
| 86 |
|
14.9/L^2 (GV/c). |
| 87 |
|
*/ |
| 88 |
|
Float_t GetCutoffSVL() { return cutoff[0]; }; |
| 89 |
|
|
| 90 |
void SetFromLevel2Struct(cOrbitalInfo *l2); |
void SetFromLevel2Struct(cOrbitalInfo *l2); |
| 91 |
void GetLevel2Struct(cOrbitalInfo *l2) const; |
void GetLevel2Struct(cOrbitalInfo *l2) const; |
| 92 |
|
|
| 93 |
void Clear(); // emiliano |
void Clear(Option_t *t=""); // emiliano |
| 94 |
// |
// |
| 95 |
ClassDef(OrbitalInfo, 1) |
ClassDef(OrbitalInfo, 3) |
| 96 |
}; |
}; |
| 97 |
#endif |
#endif |
Parent Directory
| 
Revision Log
| 
Patch