]>
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 | ||
2e5f0f7b | 7 | #include <TRotMatrix.h> |
06b2d19f | 8 | #include <TVector3.h> |
9 | #include <TMath.h> | |
853634d3 | 10 | #include <TRotation.h> |
c60862bf | 11 | #include <TLorentzVector.h> |
3ea9cb08 | 12 | #include "AliRICHConst.h" |
c64f7552 | 13 | #include "AliRICHParam.h" |
4faf338d | 14 | #include "AliRICHTresholdMap.h" |
a2f7eaf6 | 15 | #include "AliSegmentation.h" |
237c933d | 16 | #include "AliRICHGeometry.h" |
17 | #include "AliRICHResponse.h" | |
2e5f0f7b | 18 | |
0ffe413c | 19 | typedef enum {kMip, kPhoton} ResponseType; |
853634d3 | 20 | class AliRICHParam; |
2e5f0f7b | 21 | |
853634d3 | 22 | class AliRICHChamber : public TNamed |
2e5f0f7b | 23 | { |
06b2d19f | 24 | public: |
2e5f0f7b | 25 | |
06b2d19f | 26 | Int_t fIndexMap[50]; //indeces of tresholds |
27 | AliRICHTresholdMap* fTresh; //map of tresholds | |
28 | ||
29 | public: | |
3ea9cb08 | 30 | AliRICHChamber(); |
853634d3 | 31 | AliRICHChamber(Int_t iModuleN,AliRICHParam *pParam); |
3ea9cb08 | 32 | AliRICHChamber(const AliRICHChamber &chamber):TNamed(chamber) {;} |
33 | virtual ~AliRICHChamber() {;} | |
0ffe413c | 34 | AliRICHChamber& operator=(const AliRICHChamber&) {return *this;} |
35 | ||
853634d3 | 36 | TRotMatrix* RotMatrix() const{return fpRotMatrix;} |
37 | const char* RotMatrixName() const{return "rot"+fName;} | |
3ea9cb08 | 38 | TRotation Rot() {return fRot;} |
853634d3 | 39 | Double_t Rho() const{return fCenterV3.Mag();} |
c60862bf | 40 | Double_t ThetaD() const{return fCenterV3.Theta()*TMath::RadToDeg();} |
41 | Double_t PhiD() const{return fCenterV3.Phi()*TMath::RadToDeg();} | |
42 | Double_t ThetaXd() const{return fRot.ThetaX()*TMath::RadToDeg();} | |
43 | Double_t PhiXd() const{return fRot.PhiX()*TMath::RadToDeg();} | |
44 | Double_t ThetaYd() const{return fRot.ThetaY()*TMath::RadToDeg();} | |
45 | Double_t PhiYd() const{return fRot.PhiY()*TMath::RadToDeg();} | |
46 | Double_t ThetaZd() const{return fRot.ThetaZ()*TMath::RadToDeg();} | |
a277aaca | 47 | Double_t PhiZd() const{return fRot.PhiZ()*kR2d;} |
c64f7552 | 48 | void RotateX(Double_t a) {fRot.RotateX(a);fCenterV3.RotateX(a);fPcX3.RotateX(a);} |
49 | void RotateY(Double_t a) {fRot.RotateY(a);fCenterV3.RotateY(a);fPcX3.RotateY(a);} | |
50 | void RotateZ(Double_t a) {fRot.RotateZ(a);fCenterV3.RotateZ(a);fPcX3.RotateZ(a);} | |
853634d3 | 51 | Double_t X() const{return fCenterV3.X();} |
52 | Double_t Y() const{return fCenterV3.Y();} | |
53 | Double_t Z() const{return fCenterV3.Z();} | |
c60862bf | 54 | TVector3 L2G(TVector3 x3) const{x3.Transform(fRot);x3+=fCenterV3;return x3;} |
55 | TVector3 G2L(TVector3 x3) const{x3-=fCenterV3;x3.Transform(fRot.Inverse()); return x3;} | |
9e23b054 | 56 | inline TVector3 Global2Local(TVector3 x3, Bool_t isVector=kFALSE) const; |
57 | TVector3 Global2Local(TLorentzVector x4,Bool_t isVector=kFALSE) const{return Global2Local(x4.Vect(),isVector);} | |
c60862bf | 58 | TVector3 L2G(Double_t x,Double_t y,Double_t z) const{return L2G(TVector3(x,y,z));} |
59 | TVector3 G2L(TLorentzVector x4) const{return G2L(x4.Vect());} | |
60 | Float_t G2Ly(TLorentzVector x4) const{TVector3 x3=G2L(x4.Vect()); return x3.Z();} | |
61 | TVector3 G2L(Double_t x,Double_t y,Double_t z) const{return G2L(TVector3(x,y,z));} | |
62 | Float_t G2Lx(Double_t x,Double_t y,Double_t z) const{TVector3 x3=G2L(x,y,z); return x3.X();} | |
63 | Float_t G2Ly(Double_t x,Double_t y,Double_t z) const{TVector3 x3=G2L(x,y,z); return x3.Z();} | |
0ffe413c | 64 | void Print(Option_t *sOption)const;//virtual |
65 | ||
66 | void LocaltoGlobal(Float_t pos[3],Float_t Localpos[3]);//Transformation from local to global coordinates, chamber-dependant | |
67 | void GlobaltoLocal(Float_t pos[3],Float_t localpos[3]);//Transformation from Global to local coordinates, chamber-dependant | |
68 | void GenerateTresholds(); //Generate pad dependent tresholds | |
69 | void DisIntegration(Float_t eloss, Float_t xhit, Float_t yhit, Int_t&x, Float_t newclust[6][500], ResponseType res);// Cluster formation method | |
70 | void Init(Int_t id) {fSegmentation->Init(id);} // Recalculates all the values after some of them have been changed | |
71 | void SetGeometryModel(AliRICHGeometry* pRICHGeometry) {fGeometry=pRICHGeometry;} | |
72 | AliRICHGeometry* GetGeometryModel() const{return fGeometry;} | |
73 | void SetResponseModel(AliRICHResponse* pRICHResponse) {fResponse=pRICHResponse;} | |
74 | AliRICHResponse* GetResponseModel() const{return fResponse;} | |
75 | void SetSegmentationModel(AliSegmentation* pRICHSegmentation) {fSegmentation=pRICHSegmentation;} | |
76 | AliSegmentation* GetSegmentationModel() const{return fSegmentation;} | |
0ffe413c | 77 | void SigGenInit(Float_t x, Float_t y, Float_t z) {fSegmentation->SigGenInit(x, y, z) ;} |
78 | Int_t SigGenCond(Float_t x, Float_t y, Float_t z) {return fSegmentation->SigGenCond(x, y, z);} | |
79 | Int_t Sector(Float_t x, Float_t y) {return fSegmentation->Sector((Int_t)x, (Int_t)y);} // Returns number of sector containing (x,y) position | |
80 | void SetPadSize(Float_t p1, Float_t p2) {fSegmentation->SetPadSize(p1,p2);} | |
853634d3 | 81 | Double_t GetX() const{return fX;} |
82 | Double_t GetY() const{return fY;} | |
83 | Double_t GetZ() const{return fZ;} | |
c64f7552 | 84 | inline void SetToZenith(); |
85 | TRotMatrix *GetRotMatrix() const{return fpRotMatrix;} | |
853634d3 | 86 | protected: |
853634d3 | 87 | Float_t fX,fY,fZ; // Position of the center of the chamber in MRS (cm) |
237c933d | 88 | |
0ffe413c | 89 | AliSegmentation *fSegmentation; //???Segmentation model for each chamber |
90 | AliRICHResponse *fResponse; //???Response model for each chamber | |
91 | AliRICHGeometry *fGeometry; //???Geometry model for each chamber | |
06b2d19f | 92 | |
853634d3 | 93 | TVector3 fCenterV3; //chamber center position in MRS (cm) |
c64f7552 | 94 | TVector3 fPcX3; //PC center position in MRS (cm) |
853634d3 | 95 | TRotation fRot; //chamber rotation in MRS |
96 | TRotMatrix *fpRotMatrix; //rotation matrix of the chamber with respect to MRS | |
97 | AliRICHParam *fpParam; //main RICH parameters description | |
c64f7552 | 98 | ClassDef(AliRICHChamber,3) //single RICH chamber description |
853634d3 | 99 | };//class AliRICHChamber |
c64f7552 | 100 | //__________________________________________________________________________________________________ |
101 | void AliRICHChamber::SetToZenith() | |
06b2d19f | 102 | { |
c64f7552 | 103 | fCenterV3.SetXYZ(fX=0,fY=AliRICHParam::Offset()-AliRICHParam::GapThickness()/2,fZ=0); |
104 | fPcX3.SetXYZ(0,AliRICHParam::Offset()-AliRICHParam::GapThickness()/2+5.276+0.25,0); | |
853634d3 | 105 | } |
c64f7552 | 106 | //__________________________________________________________________________________________________ |
9e23b054 | 107 | TVector3 AliRICHChamber::Global2Local(TVector3 x3,Bool_t isVector)const |
108 | { | |
c64f7552 | 109 | if(!isVector) x3-=fPcX3; |
9e23b054 | 110 | x3.Transform(fRot.Inverse()); |
c64f7552 | 111 | Double_t tmp=x3.Y(); x3.SetY(x3.Z()); x3.SetZ(tmp); |
9e23b054 | 112 | return x3; |
113 | } | |
c64f7552 | 114 | //__________________________________________________________________________________________________ |
06b2d19f | 115 | #endif //AliRICHChamber_h |