]>
Commit | Line | Data |
---|---|---|
06b2d19f | 1 | #ifndef AliRICHChamber_h |
2 | #define AliRICHChamber_h | |
2e5f0f7b | 3 | |
4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
5 | * See cxx source for full Copyright notice */ | |
6 | ||
06b2d19f | 7 | #include <TVector3.h> |
8 | #include <TMath.h> | |
853634d3 | 9 | #include <TRotation.h> |
c60862bf | 10 | #include <TLorentzVector.h> |
c64f7552 | 11 | #include "AliRICHParam.h" |
a2f7eaf6 | 12 | #include "AliSegmentation.h" |
53fd478b | 13 | class AliRICHGeometry; |
14 | class AliRICHResponse; | |
15 | class TRotMatrix; | |
2e5f0f7b | 16 | |
0ffe413c | 17 | typedef enum {kMip, kPhoton} ResponseType; |
853634d3 | 18 | class AliRICHParam; |
2e5f0f7b | 19 | |
853634d3 | 20 | class AliRICHChamber : public TNamed |
2e5f0f7b | 21 | { |
06b2d19f | 22 | public: |
3ea9cb08 | 23 | AliRICHChamber(); |
853634d3 | 24 | AliRICHChamber(Int_t iModuleN,AliRICHParam *pParam); |
3ea9cb08 | 25 | AliRICHChamber(const AliRICHChamber &chamber):TNamed(chamber) {;} |
26 | virtual ~AliRICHChamber() {;} | |
0ffe413c | 27 | AliRICHChamber& operator=(const AliRICHChamber&) {return *this;} |
28 | ||
853634d3 | 29 | TRotMatrix* RotMatrix() const{return fpRotMatrix;} |
b89735dc | 30 | TString RotMatrixName() const{return "rot"+fName;} |
53fd478b | 31 | TRotation Rot() const{return fRot;} |
853634d3 | 32 | Double_t Rho() const{return fCenterV3.Mag();} |
c60862bf | 33 | Double_t ThetaD() const{return fCenterV3.Theta()*TMath::RadToDeg();} |
34 | Double_t PhiD() const{return fCenterV3.Phi()*TMath::RadToDeg();} | |
35 | Double_t ThetaXd() const{return fRot.ThetaX()*TMath::RadToDeg();} | |
36 | Double_t PhiXd() const{return fRot.PhiX()*TMath::RadToDeg();} | |
37 | Double_t ThetaYd() const{return fRot.ThetaY()*TMath::RadToDeg();} | |
38 | Double_t PhiYd() const{return fRot.PhiY()*TMath::RadToDeg();} | |
39 | Double_t ThetaZd() const{return fRot.ThetaZ()*TMath::RadToDeg();} | |
53fd478b | 40 | Double_t PhiZd() const{return fRot.PhiZ()*TMath::RadToDeg();} |
c64f7552 | 41 | void RotateX(Double_t a) {fRot.RotateX(a);fCenterV3.RotateX(a);fPcX3.RotateX(a);} |
42 | void RotateY(Double_t a) {fRot.RotateY(a);fCenterV3.RotateY(a);fPcX3.RotateY(a);} | |
43 | void RotateZ(Double_t a) {fRot.RotateZ(a);fCenterV3.RotateZ(a);fPcX3.RotateZ(a);} | |
853634d3 | 44 | Double_t X() const{return fCenterV3.X();} |
45 | Double_t Y() const{return fCenterV3.Y();} | |
46 | Double_t Z() const{return fCenterV3.Z();} | |
c60862bf | 47 | TVector3 L2G(TVector3 x3) const{x3.Transform(fRot);x3+=fCenterV3;return x3;} |
48 | TVector3 G2L(TVector3 x3) const{x3-=fCenterV3;x3.Transform(fRot.Inverse()); return x3;} | |
3582c1f9 | 49 | inline TVector2 Glob2Loc(TVector3 x3, Bool_t isVector=kFALSE) const; |
50 | TVector2 Glob2Loc(TLorentzVector x4,Bool_t isVector=kFALSE) const{return Glob2Loc(x4.Vect(),isVector);} | |
c60862bf | 51 | TVector3 L2G(Double_t x,Double_t y,Double_t z) const{return L2G(TVector3(x,y,z));} |
52 | TVector3 G2L(TLorentzVector x4) const{return G2L(x4.Vect());} | |
53 | Float_t G2Ly(TLorentzVector x4) const{TVector3 x3=G2L(x4.Vect()); return x3.Z();} | |
54 | TVector3 G2L(Double_t x,Double_t y,Double_t z) const{return G2L(TVector3(x,y,z));} | |
55 | Float_t G2Lx(Double_t x,Double_t y,Double_t z) const{TVector3 x3=G2L(x,y,z); return x3.X();} | |
56 | Float_t G2Ly(Double_t x,Double_t y,Double_t z) const{TVector3 x3=G2L(x,y,z); return x3.Z();} | |
0ffe413c | 57 | void Print(Option_t *sOption)const;//virtual |
58 | ||
b2d211fb | 59 | void LocaltoGlobal(Float_t pos[3],Float_t localpos[3]) { |
60 | //Transformation from local to global coordinates, chamber-dependant | |
61 | TVector3 buf = L2G(localpos[0],localpos[1],localpos[2]); | |
62 | pos[0]=buf.X();pos[1]=buf.Y();pos[2]=buf.Z(); | |
63 | } | |
64 | void GlobaltoLocal(Float_t pos[3],Float_t localpos[3]) { | |
65 | //Transformation from Global to local coordinates, chamber-dependant | |
66 | TVector3 buf = G2L(pos[0],pos[1],pos[2]); | |
67 | localpos[0]=buf.X();localpos[1]=buf.Y();localpos[2]=buf.Z(); | |
68 | } | |
c64f7552 | 69 | inline void SetToZenith(); |
70 | TRotMatrix *GetRotMatrix() const{return fpRotMatrix;} | |
853634d3 | 71 | protected: |
853634d3 | 72 | TVector3 fCenterV3; //chamber center position in MRS (cm) |
c64f7552 | 73 | TVector3 fPcX3; //PC center position in MRS (cm) |
853634d3 | 74 | TRotation fRot; //chamber rotation in MRS |
75 | TRotMatrix *fpRotMatrix; //rotation matrix of the chamber with respect to MRS | |
76 | AliRICHParam *fpParam; //main RICH parameters description | |
84365c70 | 77 | ClassDef(AliRICHChamber,5) //single RICH chamber description |
853634d3 | 78 | };//class AliRICHChamber |
c64f7552 | 79 | //__________________________________________________________________________________________________ |
80 | void AliRICHChamber::SetToZenith() | |
06b2d19f | 81 | { |
84365c70 | 82 | fCenterV3.SetXYZ(0,AliRICHParam::Offset()-AliRICHParam::GapThickness()/2,0); |
c64f7552 | 83 | fPcX3.SetXYZ(0,AliRICHParam::Offset()-AliRICHParam::GapThickness()/2+5.276+0.25,0); |
853634d3 | 84 | } |
c64f7552 | 85 | //__________________________________________________________________________________________________ |
3582c1f9 | 86 | TVector2 AliRICHChamber::Glob2Loc(TVector3 x3,Bool_t isVector)const |
9e23b054 | 87 | { |
c64f7552 | 88 | if(!isVector) x3-=fPcX3; |
9e23b054 | 89 | x3.Transform(fRot.Inverse()); |
3582c1f9 | 90 | return TVector2(x3.X(),x3.Z());//attention Y and Z are misplaced! |
9e23b054 | 91 | } |
c64f7552 | 92 | //__________________________________________________________________________________________________ |
06b2d19f | 93 | #endif //AliRICHChamber_h |