update version number
[u/mrichter/AliRoot.git] / RICH / AliRICHChamber.h
1 #ifndef AliRICHChamber_h
2 #define AliRICHChamber_h
3
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 #include <TVector3.h>
8 #include <TMath.h>
9 #include <TRotation.h>
10 #include <TLorentzVector.h>
11 #include "AliRICHParam.h"
12 class TRotMatrix;
13
14
15 class AliRICHChamber : public TNamed
16 {
17 public:
18            AliRICHChamber():TNamed(),fpRotMatrix(0)                      {;}
19            AliRICHChamber(Int_t iChamberN);
20            AliRICHChamber(const AliRICHChamber &chamber):TNamed(chamber) {;}
21   virtual ~AliRICHChamber()                                              {;}
22            AliRICHChamber& operator=(const AliRICHChamber&)              {return *this;}
23
24   static Double_t AlphaFeedback(Int_t )      {return 0.02;}  //determines number of feedback photons updated to 9/11/04 by Di Mauro
25  
26   TRotMatrix* RotMatrix()          const{return fpRotMatrix;}
27   TString     RotMatrixName()      const{return "rot"+fName;}
28   TRotation   Rot()                const{return fRot;}
29   Double_t    Rho()                const{return fCenter.Mag();}                                //gives  distance to chamber center in MRS
30   Double_t    ThetaD()             const{return fCenter.Theta()*TMath::RadToDeg();}            //gives polar angle of chamber center in MRS
31   Double_t    PhiD()               const{return fCenter.Phi()  *TMath::RadToDeg();}            //gives azimuthal angle of chamber center in MRS
32   Double_t    ThetaXd()            const{return fRot.ThetaX()  *TMath::RadToDeg();}    
33   Double_t    PhiXd()              const{return fRot.PhiX()    *TMath::RadToDeg();}    
34   Double_t    ThetaYd()            const{return fRot.ThetaY()  *TMath::RadToDeg();}    
35   Double_t    PhiYd()              const{return fRot.PhiY()    *TMath::RadToDeg();}    
36   Double_t    ThetaZd()            const{return fRot.ThetaZ()  *TMath::RadToDeg();}    
37   Double_t    PhiZd()              const{return fRot.PhiZ()    *TMath::RadToDeg();}    
38   void        RotX(Double_t a)       {a*=TMath::DegToRad();fRot.RotateX(a);fRad.RotateX(a);fCenter.RotateX(a);fAnod.RotateX(a);fPc.RotateX(a);}//degrees around X
39   void        RotY(Double_t a)       {a*=TMath::DegToRad();fRot.RotateY(a);fRad.RotateY(a);fCenter.RotateY(a);fAnod.RotateY(a);fPc.RotateY(a);}//degrees around Y
40   void        RotZ(Double_t a)       {a*=TMath::DegToRad();fRot.RotateZ(a);fRad.RotateZ(a);fCenter.RotateZ(a);fAnod.RotateZ(a);fPc.RotateZ(a);}//degrees around Z
41   TVector3    Rad()               const{return fRad;}         //provides center of radiator position in MRS, cm   
42   TVector3    Anod()              const{return fAnod;}        //provides center of anod wires plane in MRS, cm   
43   TVector3    Pc()                const{return fPc;}          //provides center of photocathode position in MRS, cm
44   TVector3    Center()            const{return fCenter;}      //provides center of chamber position (exit from quartz window) in MRS, cm
45   void        Print(Option_t *sOption)const;                    
46   TVector3    PMrs2Loc(TVector3 p3)const{TVector3 ploc=Rot().Invert()*p3;ploc.SetXYZ(-ploc.Px(),ploc.Py(),ploc.Pz()); return ploc;} //momentum MARS-local 
47 //Transformations for radiator plane  
48   TVector2    Mrs2Rad(TVector3 x3)const{x3-=fRad;x3.Transform(fRot.Inverse());return TVector2(-x3.X()+0.5*AliRICHParam::PcSizeX(),x3.Y()+0.5*AliRICHParam::PcSizeY());}
49   TVector3    Rad2Mrs(TVector2 x2)const{TVector3 x3(-x2.X()+0.5*AliRICHParam::PcSizeX(),x2.Y()-0.5*AliRICHParam::PcSizeY(),0);x3.Transform(fRot); x3+=fRad;return x3;}  
50 //Transformations for anod wires plane  
51   TVector2    Mrs2Anod(TVector3 x3)const{x3-=fAnod;x3.Transform(fRot.Inverse());return TVector2(-x3.X()+0.5*AliRICHParam::PcSizeX(),x3.Y()+0.5*AliRICHParam::PcSizeY());}
52   TVector3    Anod2Mrs(TVector2 x2)const{TVector3 x3(-x2.X()+0.5*AliRICHParam::PcSizeX(),x2.Y()-0.5*AliRICHParam::PcSizeY(),0);x3.Transform(fRot); x3+=fAnod;return x3;}  
53 //Transformations for photcathode plane  
54   TVector2    Mrs2Pc(TVector3 x3)const{x3-=fPc;x3.Transform(fRot.Inverse());return TVector2(-x3.X()+0.5*AliRICHParam::PcSizeX(),x3.Y()+0.5*AliRICHParam::PcSizeY());}
55   TVector3    Pc2Mrs(TVector2 x2)const{TVector3 x3(-x2.X()+0.5*AliRICHParam::PcSizeX(),x2.Y()-0.5*AliRICHParam::PcSizeY(),0);x3.Transform(fRot); x3+=fPc;return x3;}  
56   TVector2    Mrs2Pc(TLorentzVector x4)            const{return Mrs2Pc(x4.Vect());}
57 protected:
58   TVector3      fRad;             //radiator entrance center position in MRS (cm)
59   TVector3      fCenter;          //chamber center position (quartz window exit) in MRS (cm) 
60   TVector3      fAnod;            //anod wires plane center position in MRS (cm)
61   TVector3      fPc;              //PC center position in MRS (cm)
62   TRotation     fRot;             //chamber rotation in MRS
63   TRotMatrix   *fpRotMatrix;      //rotation matrix of the chamber with respect to MRS 
64   ClassDef(AliRICHChamber,8)      //single RICH chamber description
65 };//class AliRICHChamber
66
67 #endif //AliRICHChamber_h