OrbitalInfo/inc/InclinationInfo.h

/***************************************************************************
 *   Copyright (C) 2006 by pamelaprod                                      *
 *   pamelaprod@P1.pamela                                                  *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
#ifndef FIRST_H
#define FIRST_H

#include <TObject.h>
#include <iostream>
#include <TArrayC.h>
#include <TROOT.h>

class InclinationInfoI : public TObject {
 public:
  //Quaternions Lx_y
     // The data are organized in 5 groups according to the format
     // time    -   L0       -     L1     -     L2     -     L3
     // time[0] - quat[0][0] - quat[0][1] - quat[0][2] - quat[0][3]
     // time[1] - quat[1][0] - quat[1][1] - quat[1][2] - quat[1][3]
     // time[2] - quat[2][0] - quat[2][1] - quat[2][2] - quat[2][3]
     // time[3] - quat[3][0] - quat[3][1] - quat[3][2] - quat[3][3]
     // time[4] - quat[4][0] - quat[4][1] - quat[4][2] - quat[4][3]
     // time[5] - quat[5][0] - quat[5][1] - quat[5][2] - quat[5][3]
  //
  double time[6]; 
  double quat[6][4]; 
  
  InclinationInfoI(){}; 
  //  ~InclinationInfoI() {};
  
  void fill(TArrayC* data);
  //  const char* toXML(char*);
  ClassDef(InclinationInfoI,1)
};

/**
   @author pamelaprod
*/
class Quaternions: public InclinationInfoI {
 public:
  
  Quaternions();    
  ~Quaternions();
  
  //Int_t CodeErrQua;    

  ClassDef(Quaternions,1)
};

class InclinationInfo: public TObject {
//private:
    //Double_t k1;
    //Double_t k2;
    //void coefplane(Double_t x1, Double_t y1, Double_t z1, Double_t x2, Double_t y2, Double_t z2);
    //Double_t AngBetAxes(Double_t x1, Double_t y1, Double_t z1, Double_t x2, Double_t y2, Double_t z2);
    //Double_t ValueT(Double_t x1,Double_t y1, Double_t z1, Double_t n1, Double_t n2);
    //Double_t AngBetPlan(Double_t x1,Double_t y1, Double_t z1, Double_t n1, Double_t n2);
    
public:

    InclinationInfo();
    ~InclinationInfo();
    
    // Nico's variables //
    //Int_t nicos1;
    //Int_t nicos2;
    //Int_t nicos3;
    //Int_t nicos4;
    // End nico's //
    
    // Our variables //
    Float_t A11;
    Float_t A12;
    Float_t A13;
//    float_t A21;
//    float_t A22;
//    float_t A23;
// Angles RCS of OCS    
//    float_t AXrXo;
//    float_t AXrYo;
//    float_t AXrZo;
//    float_t AYrXo;
//    float_t AYrYo;
//    float_t AYrZo;
//    float_t AZrXo;
//    float_t AZrYo;
//    float_t AZrZo;
    
//    float_t AXrZXo;
//    float_t AXrXYo;
//    float_t AXrYZo;
//    float_t AYrZXo;
//    float_t AYrXYo;
//    float_t AYrYZo;
//    float_t AZrZXo;
//    float_t AZrXYo;
//    float_t AZrYZo;
    
//    float_t ApXrZXoZ;
//    float_t ApXrXYoZ;
//    float_t ApXrYZoZ;
//    float_t ApYrZXoZ;
//    float_t ApYrXYoZ;
//    float_t ApYrYZoZ;
//    float_t ApZrZXoZ;
//    float_t ApZrXYoZ;
//    float_t ApZrYZoZ;
    
//    float_t AXrZXo_2m;
//    float_t AXrXYo_2m;
//    float_t AXrYZo_2m;
//    float_t AYrZXo_2m;
//    float_t AYrXYo_2m;
//    float_t AYrYZo_2m;
//    float_t AZrZXo_2m;
//    float_t AZrXYo_2m;
//    float_t AZrYZo_2m;
    
//    float_t ApXrZXoZ_2m;
//    float_t ApXrXYoZ_2m;
//    float_t ApXrYZoZ_2m;
//    float_t ApYrZXoZ_2m;
//    float_t ApYrXYoZ_2m;
//    float_t ApYrYZoZ_2m;
//    float_t ApZrZXoZ_2m;
//    float_t ApZrXYoZ_2m;
//    float_t ApZrYZoZ_2m;
    
    Double_t XXRCS;
    Double_t YXRCS;
    Double_t ZXRCS;
    Double_t XYRCS;
    Double_t YYRCS;
    Double_t ZYRCS;
    Double_t XZRCS;
    Double_t YZRCS;
    Double_t ZZRCS;
    
    Double_t Tangazh;
    Double_t Ryskanie;
    Double_t Kren;
    
    bool mode;
    Float_t dircutoff;
    
    short int Sign_1(Double_t a, short int b){
    if (a>0){b=1;};
    if (a<0){b=-1;}
    else {b=0;};
    return b;};
    void QuaternionstoAngle(Quaternions);
    void TransAngle(Double_t x,Double_t y,Double_t z,Double_t Vx0,Double_t Vy0, Double_t Vz0, Double_t gamar, Double_t ksir, Double_t tetar, Double_t q0, Double_t q1, Double_t q2, Double_t q3);
    
//    float_t procRegidity(double_t L, double_t B, double_t R);
    
//    UInt_t absTime;
//    UInt_t OBT;
//    UInt_t pkt_num; 
    
    void Clear();
    ClassDef(InclinationInfo,1)
};

#endif
1	pamelaprod	1.1	/***************************************************************************
2			* Copyright (C) 2006 by pamelaprod *
3			* pamelaprod@P1.pamela *
4			* *
5			* This program is free software; you can redistribute it and/or modify *
6			* it under the terms of the GNU General Public License as published by *
7			* the Free Software Foundation; either version 2 of the License, or *
8			* (at your option) any later version. *
9			* *
10			* This program is distributed in the hope that it will be useful, *
11			* but WITHOUT ANY WARRANTY; without even the implied warranty of *
12			* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
13			* GNU General Public License for more details. *
14			* *
15			* You should have received a copy of the GNU General Public License *
16			* along with this program; if not, write to the *
17			* Free Software Foundation, Inc., *
18			* 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. *
19			***************************************************************************/
20			#ifndef FIRST_H
21			#define FIRST_H
22
23			#include <TObject.h>
24			#include <iostream>
25			#include <TArrayC.h>
26			#include <TROOT.h>
27
28			class InclinationInfoI : public TObject {
29			public:
30			//Quaternions Lx_y
31			// The data are organized in 5 groups according to the format
32			// time - L0 - L1 - L2 - L3
33			// time[0] - quat[0][0] - quat[0][1] - quat[0][2] - quat[0][3]
34			// time[1] - quat[1][0] - quat[1][1] - quat[1][2] - quat[1][3]
35			// time[2] - quat[2][0] - quat[2][1] - quat[2][2] - quat[2][3]
36			// time[3] - quat[3][0] - quat[3][1] - quat[3][2] - quat[3][3]
37			// time[4] - quat[4][0] - quat[4][1] - quat[4][2] - quat[4][3]
38			// time[5] - quat[5][0] - quat[5][1] - quat[5][2] - quat[5][3]
39			//
40			double time[6];
41			double quat[6][4];
42
43			InclinationInfoI(){};
44			// ~InclinationInfoI() {};
45
46			void fill(TArrayC* data);
47			// const char* toXML(char*);
48			ClassDef(InclinationInfoI,1)
49			};
50
51			/**
52			@author pamelaprod
53			*/
54			class Quaternions: public InclinationInfoI {
55			public:
56
57			Quaternions();
58			~Quaternions();
59
60			//Int_t CodeErrQua;
61
62			ClassDef(Quaternions,1)
63			};
64
65			class InclinationInfo: public TObject {
66			//private:
67			//Double_t k1;
68			//Double_t k2;
69			//void coefplane(Double_t x1, Double_t y1, Double_t z1, Double_t x2, Double_t y2, Double_t z2);
70			//Double_t AngBetAxes(Double_t x1, Double_t y1, Double_t z1, Double_t x2, Double_t y2, Double_t z2);
71			//Double_t ValueT(Double_t x1,Double_t y1, Double_t z1, Double_t n1, Double_t n2);
72			//Double_t AngBetPlan(Double_t x1,Double_t y1, Double_t z1, Double_t n1, Double_t n2);
73
74			public:
75
76			InclinationInfo();
77			~InclinationInfo();
78
79			// Nico's variables //
80			//Int_t nicos1;
81			//Int_t nicos2;
82			//Int_t nicos3;
83			//Int_t nicos4;
84			// End nico's //
85
86			// Our variables //
87			Float_t A11;
88			Float_t A12;
89			Float_t A13;
90			// float_t A21;
91			// float_t A22;
92			// float_t A23;
93			// Angles RCS of OCS
94			// float_t AXrXo;
95			// float_t AXrYo;
96			// float_t AXrZo;
97			// float_t AYrXo;
98			// float_t AYrYo;
99			// float_t AYrZo;
100			// float_t AZrXo;
101			// float_t AZrYo;
102			// float_t AZrZo;
103
104			// float_t AXrZXo;
105			// float_t AXrXYo;
106			// float_t AXrYZo;
107			// float_t AYrZXo;
108			// float_t AYrXYo;
109			// float_t AYrYZo;
110			// float_t AZrZXo;
111			// float_t AZrXYo;
112			// float_t AZrYZo;
113
114			// float_t ApXrZXoZ;
115			// float_t ApXrXYoZ;
116			// float_t ApXrYZoZ;
117			// float_t ApYrZXoZ;
118			// float_t ApYrXYoZ;
119			// float_t ApYrYZoZ;
120			// float_t ApZrZXoZ;
121			// float_t ApZrXYoZ;
122			// float_t ApZrYZoZ;
123
124			// float_t AXrZXo_2m;
125			// float_t AXrXYo_2m;
126			// float_t AXrYZo_2m;
127			// float_t AYrZXo_2m;
128			// float_t AYrXYo_2m;
129			// float_t AYrYZo_2m;
130			// float_t AZrZXo_2m;
131			// float_t AZrXYo_2m;
132			// float_t AZrYZo_2m;
133
134			// float_t ApXrZXoZ_2m;
135			// float_t ApXrXYoZ_2m;
136			// float_t ApXrYZoZ_2m;
137			// float_t ApYrZXoZ_2m;
138			// float_t ApYrXYoZ_2m;
139			// float_t ApYrYZoZ_2m;
140			// float_t ApZrZXoZ_2m;
141			// float_t ApZrXYoZ_2m;
142			// float_t ApZrYZoZ_2m;
143
144			Double_t XXRCS;
145			Double_t YXRCS;
146			Double_t ZXRCS;
147			Double_t XYRCS;
148			Double_t YYRCS;
149			Double_t ZYRCS;
150			Double_t XZRCS;
151			Double_t YZRCS;
152			Double_t ZZRCS;
153
154			Double_t Tangazh;
155			Double_t Ryskanie;
156			Double_t Kren;
157
158			bool mode;
159			Float_t dircutoff;
160
161			short int Sign_1(Double_t a, short int b){
162			if (a>0){b=1;};
163			if (a<0){b=-1;}
164			else {b=0;};
165			return b;};
166			void QuaternionstoAngle(Quaternions);
167			void TransAngle(Double_t x,Double_t y,Double_t z,Double_t Vx0,Double_t Vy0, Double_t Vz0, Double_t gamar, Double_t ksir, Double_t tetar, Double_t q0, Double_t q1, Double_t q2, Double_t q3);
168
169			// float_t procRegidity(double_t L, double_t B, double_t R);
170
171			// UInt_t absTime;
172			// UInt_t OBT;
173			// UInt_t pkt_num;
174
175			void Clear();
176			ClassDef(InclinationInfo,1)
177			};
178
179			#endif