1 #ifndef ALIITSGEOMMATRIX_H
2 #define ALIITSGEOMMATRIX_H
3 /* Copyright(c) 2000, ALICE Experiment at CERN, All rights reserved. *
4 * see cxx source for full Copyright notice. */
6 ////////////////////////////////////////////////////////////////////////
7 // ITS geometry manipulation routines on the module level. This class is
8 // to replace the structure ITS_geom in the class AliITSgeom.
9 // Created May 30 2000.
12 ////////////////////////////////////////////////////////////////////////
16 class AliITSgeomMatrix : public TObject {
18 AliITSgeomMatrix(); // Default constructor
19 // Standard constructor #1
20 AliITSgeomMatrix(Int_t idt,const Int_t id[3],
21 const Double_t rot[3],const Double_t tran[3]);
22 // Standard constructor #2
23 AliITSgeomMatrix(Int_t idt,const Int_t id[3],
24 Double_t matrix[3][3],const Double_t tran[3]);
25 // Standard constructor #3
26 AliITSgeomMatrix(const Double_t rotd[6]/*degrees GEANT angles*/,
27 Int_t idt,const Int_t id[3],
28 const Double_t tran[3]);
30 AliITSgeomMatrix(const AliITSgeomMatrix &source);
32 void operator=(const AliITSgeomMatrix &sourse); // copy
33 virtual ~AliITSgeomMatrix(){}; // default constructor.
34 // Prints a line describing the output format of the function Print.
35 void PrintComment(ostream *os) const;
36 // Prints out the content of this class in ASCII format.
37 void Print(ostream *os);
38 // Prints out the content of this class in ASCII format but includes
39 // formating and strings that make it more humanly readable.
40 void PrintTitles(ostream *os) const;
41 // Reads in the content of this class in the format of Print
42 void Read(istream *is);
43 virtual void Print(Option_t *option="") const {TObject::Print(option);}
44 virtual Int_t Read(const char *name) {return TObject::Read(name);}
46 // Given the rotation angles [radians] it fills frot and computes
47 // the rotation matrix fm.
48 void SetAngles(const Double_t rot[3]){// [radians]
49 for(Int_t i=0;i<3;i++)frot[i] = rot[i];this->MatrixFromAngle();}
50 // Sets the translation vector and computes fCylR and fCylPhi.
51 void SetTranslation(const Double_t tran[3]);
52 // sets the rotation matrix and computes the rotation angles [radians]
53 void SetMatrix(Double_t matrix[3][3]){ for(Int_t i=0;i<3;i++)
54 for(Int_t j=0;j<3;j++) fm[i][j]=matrix[i][j];this->AngleFromMatrix();}
55 // Sets the detector index value
56 void SetDetectorIndex(Int_t idt) {fDetectorIndex = idt;}
57 // Sets the detector layer, ladder, detector (id) values.
58 void SetIndex(const Int_t id[3]){
59 for(Int_t i=0;i<3;i++) fid[i] = id[i];}
60 // Returns the rotation angles [radians]
61 void GetAngles(Double_t rot[3]) const {// [radians]
62 for(Int_t i=0;i<3;i++) rot[i] = frot[i];}
63 // Returns the translation vector [cm]
64 void GetTranslation(Double_t tran[3]) const {
65 for(Int_t i=0;i<3;i++) tran[i] = ftran[i];}
66 // Returns the translation vector in cylindrical
67 // coordinates [cm,radians]
68 void GetTranslationCylinderical (Double_t tran[3]) const {
72 // Returns the values of the rotation matrix
73 void GetMatrix(Double_t matrix[3][3]) const {for(Int_t i=0;i<3;i++)
74 for(Int_t j=0;j<3;j++) matrix[i][j] = fm[i][j];}
75 // Returns the detector index value.
76 Int_t GetDetectorIndex() const {return fDetectorIndex;}
77 // returns the modules index layer, ladder, detector
78 void GetIndex(Int_t id[3]) const {for(Int_t i=0;i<3;i++) id[i] = fid[i];}
79 // Sets the rotation matrix based on the 6 GEANT rotation
81 void MatrixFromSixAngles(const Double_t *ang);
82 // Returns the 6 GEANT rotation angles [radians] from the
83 // existing rotation matrix.
84 void SixAnglesFromMatrix(Double_t *ang);
86 // Given a position in Cartesian ALICE global coordinates [cm]
87 // returns the position in Cartesian detector/module local
89 void GtoLPosition(const Double_t g[3],Double_t l[3]) const;
90 // Given a position in Cartesian detector/module local coordinates [cm]
91 // returns the position in Cartesian ALICE global
93 void LtoGPosition(const Double_t l[3],Double_t g[3]) const;
94 // Given a momentum in Cartesian ALICE global coordinates
95 // returns the momentum in Cartesian detector/module local
97 void GtoLMomentum(const Double_t g[3],Double_t l[3]) const;
98 // Given a momentum in Cartesian detector/module local coordinates
99 // returns the momentum in Cartesian ALICE global coordinates
100 void LtoGMomentum(const Double_t l[3],Double_t g[3]) const;
101 // given a position error matrix in ALICE Cartesian global
102 // coordinates [cm] returns a position error matrix in detector/
103 // module local Cartesian local coordinates [cm]
104 void GtoLPositionError(Double_t g[3][3],Double_t l[3][3]) const;
105 // given a position error matrix in detector/module Cartesian local
106 // coordinates [cm] returns a position error matrix in ALICE
107 // Cartesian global coordinates [cm]
108 void LtoGPositionError(Double_t l[3][3],Double_t g[3][3]) const;
109 // Tracking Related Routines
110 void GtoLPositionTracking(const Double_t g[3],Double_t l[3]) const;
111 // Given a position in Cartesian Tracking global coordinates [cm]
112 // returns the position in Cartesian detector/module local
114 void LtoGPositionTracking(const Double_t l[3],Double_t g[3]) const;
115 // Given a position in Cartesian detector/module local coordinates [cm]
116 // returns the position in Cartesian Tracking global
118 void GtoLMomentumTracking(const Double_t g[3],Double_t l[3]) const;
119 // Given a momentum in Cartesian detector/module local coordinates
120 // returns the momentum in Cartesian Tracking global coordinates
121 void LtoGMomentumTracking(const Double_t l[3],Double_t g[3]) const;
122 // given a position error matrix in Tracking Cartesian global
123 // coordinates [cm] returns a position error matrix in detector/
124 // module local Cartesian local coordinates [cm]
125 void GtoLPositionErrorTracking(Double_t g[3][3],
126 Double_t l[3][3]) const;
127 // given a position error matrix in detector/module Cartesian local
128 // coordinates [cm] returns a position error matrix in Tracking
129 // Cartesian global coordinates [cm]
130 void LtoGPositionErrorTracking(Double_t l[3][3],
131 Double_t g[3][3]) const;
132 // Computes the distance squared [cm^2] between a point t[3] and
133 // this module/detector
134 Double_t Distance2(const Double_t t[3]) const {Double_t d=0.0,q;
135 for(Int_t i=0;i<3;i++){q = t[i]-ftran[i]; d += q*q;}
137 private: // private functions
138 // Given the rotation matrix fm it fills the rotation angles frot
139 void MatrixFromAngle();
140 // Given the rotation angles frot it fills the rotation matrix fm
141 void AngleFromMatrix();
142 private: // Data members.
143 Int_t fDetectorIndex; // Detector type index (like fShapeIndex was)
144 Int_t fid[3]; // layer, ladder, detector numbers.
145 Double_t frot[3]; //! vector of rotations about x,y,z [radians].
146 Double_t ftran[3]; // Translation vector of module x,y,z.
147 Double_t fCylR,fCylPhi; //! Translation vector in Cylindrical coord.
148 Double_t fm[3][3]; // Rotation matrix based on frot.
150 // Note, fCylR and fCylPhi are added as data members because it costs
151 // about a factor of 10 to compute them over looking them up. Since
152 // they are used in some tracking algorithms this can be a large cost
153 // in computing time. They are not written out but computed.
155 ClassDef(AliITSgeomMatrix,1) // Matrix class used by AliITSgeom.
157 // Input and output function for standard C++ input/output.
158 ostream &operator<<(ostream &os,AliITSgeomMatrix &source);
159 istream &operator>>(istream &os,AliITSgeomMatrix &source);