]> git.uio.no Git - u/mrichter/AliRoot.git/blame - RICH/AliRICHChamber.h
Using the same optimization in the compiler and in the linker
[u/mrichter/AliRoot.git] / RICH / AliRICHChamber.h
CommitLineData
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"
2e5f0f7b 13
4faf338d 14#include "AliRICHTresholdMap.h"
a2f7eaf6 15#include "AliSegmentation.h"
237c933d 16#include "AliRICHGeometry.h"
17#include "AliRICHResponse.h"
2e5f0f7b 18
0ffe413c 19typedef enum {kMip, kPhoton} ResponseType;
853634d3 20class AliRICHParam;
2e5f0f7b 21
853634d3 22class AliRICHChamber : public TNamed
2e5f0f7b 23{
06b2d19f 24public:
2e5f0f7b 25
06b2d19f 26 Int_t fIndexMap[50]; //indeces of tresholds
27 AliRICHTresholdMap* fTresh; //map of tresholds
28
29public:
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;}
853634d3 48 void RotateX(Double_t a) {fRot.RotateX(a);fCenterV3.RotateX(a);}
49 void RotateY(Double_t a) {fRot.RotateY(a);fCenterV3.RotateY(a);}
50 void RotateZ(Double_t a) {fRot.RotateZ(a);fCenterV3.RotateZ(a);}
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;}
853634d3 84 inline void SetCenter(Double_t x,Double_t y,Double_t z);
85 TRotMatrix *GetRotMatrix() const{return fpRotMatrix;}
c60862bf 86 void SetChamberTransform(Float_t x,Float_t y,Float_t z,TRotMatrix *pRotMatrix);
853634d3 87
853634d3 88protected:
853634d3 89 Float_t fX,fY,fZ; // Position of the center of the chamber in MRS (cm)
237c933d 90
0ffe413c 91 AliSegmentation *fSegmentation; //???Segmentation model for each chamber
92 AliRICHResponse *fResponse; //???Response model for each chamber
93 AliRICHGeometry *fGeometry; //???Geometry model for each chamber
06b2d19f 94
853634d3 95 TVector3 fCenterV3; //chamber center position in MRS (cm)
96 TRotation fRot; //chamber rotation in MRS
97 TRotMatrix *fpRotMatrix; //rotation matrix of the chamber with respect to MRS
98 AliRICHParam *fpParam; //main RICH parameters description
ac4ea1ec 99 ClassDef(AliRICHChamber,2) //single RICH chamber description
853634d3 100};//class AliRICHChamber
06b2d19f 101
853634d3 102void AliRICHChamber::SetCenter(Double_t x,Double_t y,Double_t z)
06b2d19f 103{
853634d3 104 fCenterV3.SetXYZ(x,y,z);
105 fX=x;fY=y;fZ=z;
106}
9e23b054 107
108TVector3 AliRICHChamber::Global2Local(TVector3 x3,Bool_t isVector)const
109{
110 if(!isVector) x3-=fCenterV3;
111 x3.Transform(fRot.Inverse());
112 Double_t tmp=x3.Y(); x3.SetY(-x3.Z()); x3.SetZ(tmp);
113 return x3;
114}
853634d3 115
06b2d19f 116#endif //AliRICHChamber_h