4 |
#include <TObject.h> |
#include <TObject.h> |
5 |
#include <algorithm> // EMILIANO |
#include <algorithm> // EMILIANO |
6 |
#include <OrbitalInfoStruct.h> |
#include <OrbitalInfoStruct.h> |
7 |
|
#include <TClonesArray.h> |
8 |
|
#include <TMatrixD.h> |
9 |
|
|
10 |
|
class OrbitalInfoTrkVar : public TObject { |
11 |
|
private: |
12 |
|
|
13 |
|
public: |
14 |
|
// |
15 |
|
Int_t trkseqno; // tof sequ. number: -1=ToF standalone, 0=first Tracker track, ... |
16 |
|
// |
17 |
|
Float_t pitch; ///< Pitch angle |
18 |
|
// |
19 |
|
TMatrixD Eij; ///< vector of incoming particle respect to cartesian geographic coordinates |
20 |
|
TMatrixD Sij; ///< vector of incoming particle respect to flight coordinates |
21 |
|
// |
22 |
|
Float_t cutoff; ///< Calculated cutoff for the incoming particle taking into account particle direction |
23 |
|
// |
24 |
|
OrbitalInfoTrkVar(); |
25 |
|
OrbitalInfoTrkVar* GetOrbitalInfoTrkVar(){return this;}; |
26 |
|
// |
27 |
|
void Clear(Option_t *t=""); |
28 |
|
void Delete(Option_t *t=""); //ELENA |
29 |
|
// |
30 |
|
ClassDef(OrbitalInfoTrkVar, 2); |
31 |
|
// |
32 |
|
}; |
33 |
|
|
34 |
|
|
35 |
|
|
36 |
/** |
/** |
37 |
* Class that stores position, time, inclination, magnetic field and |
* Class that stores position, time, inclination, magnetic field and |
38 |
* cutoff informations. |
* cutoff informations. |
39 |
*/ |
*/ |
40 |
class OrbitalInfo : public TObject { |
class OrbitalInfo : public TObject { |
41 |
|
private: |
42 |
|
|
43 |
public: |
public: |
44 |
OrbitalInfo(); |
TClonesArray *OrbitalInfoTrk; |
|
OrbitalInfo* GetOrbitalInfo(){return this;}; // Elena |
|
45 |
|
|
46 |
UInt_t absTime; ///< Absolute Time (seconds) |
UInt_t absTime; ///< Absolute Time (seconds) |
47 |
UInt_t OBT; ///< On Board Time (ms) |
UInt_t OBT; ///< On Board Time (ms) |
62 |
|
|
63 |
Float_t L; ///< McIlwain's L shell (in earth radii) |
Float_t L; ///< McIlwain's L shell (in earth radii) |
64 |
|
|
65 |
// Dipolar magnetic coordinates (not used). |
/* // Dipolar magnetic coordinates (not used). */ |
66 |
Float_t londip; ///< degrees from -180 to 180 |
/* Float_t londip; ///< degrees from -180 to 180 */ |
67 |
Float_t latdip; ///< degrees from -90 to 90 |
/* Float_t latdip; ///< degrees from -90 to 90 */ |
68 |
Float_t altdip; ///< meters |
/* Float_t altdip; ///< meters */ |
69 |
|
|
70 |
// Corrected magnetic coordinates (not used). |
/* // Corrected magnetic coordinates (not used). */ |
71 |
Float_t loncgm; ///< degrees from -180 to 180 |
/* Float_t loncgm; ///< degrees from -180 to 180 */ |
72 |
Float_t latcgm; ///< degrees from -90 to 90 |
/* Float_t latcgm; ///< degrees from -90 to 90 */ |
73 |
Float_t altcgm; ///< meters |
/* Float_t altcgm; ///< meters */ |
74 |
|
|
75 |
// Corrected B min magnetic coordinates (not used). |
/* // Corrected B min magnetic coordinates (not used). */ |
76 |
Float_t loncbm; ///< degrees from -180 to 180 |
/* Float_t loncbm; ///< degrees from -180 to 180 */ |
77 |
Float_t latcbm; ///< degrees from -90 to 90 |
/* Float_t latcbm; ///< degrees from -90 to 90 */ |
78 |
Float_t altcbm; ///< meters |
/* Float_t altcbm; ///< meters */ |
79 |
|
|
80 |
Float_t cutoff[20]; |
// Float_t cutoff[17]; |
81 |
|
Float_t cutoffsvl; |
82 |
|
|
83 |
// Quaternions |
// Quaternions |
84 |
Float_t q0; ///< Quaternion 0 |
Float_t q0; ///< Quaternion 0 |
91 |
Float_t phi; ///< Euler angle phi in the velocity reference frame (yaw) |
Float_t phi; ///< Euler angle phi in the velocity reference frame (yaw) |
92 |
Float_t etha; ///< Euler angle etha in the velocity reference frame (roll) |
Float_t etha; ///< Euler angle etha in the velocity reference frame (roll) |
93 |
|
|
94 |
|
// Pitch angles |
95 |
|
// Float_t pamzenitangle; |
96 |
|
// Float_t pamBangle; |
97 |
|
// |
98 |
|
TMatrixD Iij; ///< Angle between PAMELA Z direction and cartesian geographic coordinates |
99 |
|
|
100 |
/** |
/** |
101 |
* The variable mode means a character time distant between two quaternions, inside which stay every events |
* The variable mode means a character time distant between two quaternions, inside which stay every events |
102 |
*/ |
*/ |
111 |
// 7 - means that we have uneliminable hole inside R10 |
// 7 - means that we have uneliminable hole inside R10 |
112 |
// 8 - means that we have eliminable hole inside non R10 |
// 8 - means that we have eliminable hole inside non R10 |
113 |
// 9 - means that we have uneliminable hole inside non R10 |
// 9 - means that we have uneliminable hole inside non R10 |
114 |
|
// 10 - means other unknown problems |
115 |
|
|
116 |
|
Int_t ntrk(){return OrbitalInfoTrk->GetEntries();}; |
117 |
/** |
/** |
118 |
\return Returns the B minimum along the field line. |
\return Returns the B minimum along the field line. |
119 |
*/ |
*/ |
123 |
\return Returns the Stormer vertical cutoff using L shell: |
\return Returns the Stormer vertical cutoff using L shell: |
124 |
14.9/L^2 (GV/c). |
14.9/L^2 (GV/c). |
125 |
*/ |
*/ |
126 |
Float_t GetCutoffSVL() { return cutoff[0]; }; |
Float_t GetCutoffSVL() { return cutoffsvl; }; |
127 |
|
|
128 |
void SetFromLevel2Struct(cOrbitalInfo *l2); |
void SetFromLevel2Struct(cOrbitalInfo *l2); |
129 |
void GetLevel2Struct(cOrbitalInfo *l2) const; |
void GetLevel2Struct(cOrbitalInfo *l2) const; |
130 |
|
OrbitalInfoTrkVar *GetOrbitalInfoTrkVar(Int_t notrack); |
131 |
|
// |
132 |
|
OrbitalInfo(); |
133 |
|
~OrbitalInfo(){Delete();}; //ELENA |
134 |
|
// |
135 |
|
OrbitalInfo* GetOrbitalInfo(){return this;}; // Elena |
136 |
|
void Delete(Option_t *t=""); //ELENA |
137 |
|
void Set();//ELENA |
138 |
|
// |
139 |
|
// |
140 |
|
|
141 |
void Clear(Option_t *t=""); // emiliano |
void Clear(Option_t *t=""); // emiliano |
142 |
// |
// |
143 |
ClassDef(OrbitalInfo, 3) |
ClassDef(OrbitalInfo, 7); |
144 |
}; |
}; |
145 |
#endif |
#endif |