| 5 |
#include <algorithm> // EMILIANO |
#include <algorithm> // EMILIANO |
| 6 |
#include <OrbitalInfoStruct.h> |
#include <OrbitalInfoStruct.h> |
| 7 |
#include <TClonesArray.h> |
#include <TClonesArray.h> |
| 8 |
|
#include <TMatrixD.h> |
| 9 |
|
|
| 10 |
class OrbitalInfoTrkVar : public TObject { |
class OrbitalInfoTrkVar : public TObject { |
| 11 |
private: |
private: |
| 14 |
// |
// |
| 15 |
Int_t trkseqno; // tof sequ. number: -1=ToF standalone, 0=first Tracker track, ... |
Int_t trkseqno; // tof sequ. number: -1=ToF standalone, 0=first Tracker track, ... |
| 16 |
// |
// |
| 17 |
Float_t pitch; |
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(); |
OrbitalInfoTrkVar(); |
| 25 |
OrbitalInfoTrkVar* GetOrbitalInfoTrkVar(){return this;}; |
OrbitalInfoTrkVar* GetOrbitalInfoTrkVar(){return this;}; |
| 26 |
// |
// |
| 27 |
void Clear(Option_t *t=""); |
void Clear(Option_t *t=""); |
| 28 |
|
void Delete(Option_t *t=""); //ELENA |
| 29 |
// |
// |
| 30 |
ClassDef(OrbitalInfoTrkVar, 1); |
ClassDef(OrbitalInfoTrkVar, 2); |
| 31 |
// |
// |
| 32 |
}; |
}; |
| 33 |
|
|
| 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[17]; |
// Float_t cutoff[17]; |
| 81 |
|
Float_t cutoffsvl; |
| 82 |
|
|
| 83 |
// Quaternions |
// Quaternions |
| 84 |
Float_t q0; ///< Quaternion 0 |
Float_t q0; ///< Quaternion 0 |
| 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 |
// Pitch angles |
| 95 |
Float_t pamzenitangle; |
// Float_t pamzenitangle; |
| 96 |
Float_t pamBangle; |
// 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 |
| 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(); |
OrbitalInfo(); |
| 133 |
~OrbitalInfo(){Delete();}; //ELENA |
~OrbitalInfo(){Delete();}; //ELENA |
| 140 |
|
|
| 141 |
void Clear(Option_t *t=""); // emiliano |
void Clear(Option_t *t=""); // emiliano |
| 142 |
// |
// |
| 143 |
ClassDef(OrbitalInfo, 5); |
ClassDef(OrbitalInfo, 7); |
| 144 |
}; |
}; |
| 145 |
#endif |
#endif |