X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=RICH%2FAliRICHChamber.h;h=006bb8bd98cb383466436a5020beba99452d7603;hb=2276b46405b2518866fd676e2f144269be13ea4a;hp=a32233622a4bee59919eb0d52bb42147acf87f9b;hpb=53fd478b685fa8c993b49217c5b9f426ca9fec7c;p=u%2Fmrichter%2FAliRoot.git diff --git a/RICH/AliRICHChamber.h b/RICH/AliRICHChamber.h index a32233622a4..006bb8bd98c 100644 --- a/RICH/AliRICHChamber.h +++ b/RICH/AliRICHChamber.h @@ -9,99 +9,59 @@ #include #include #include "AliRICHParam.h" -#include "AliSegmentation.h" -class AliRICHGeometry; -class AliRICHResponse; class TRotMatrix; -typedef enum {kMip, kPhoton} ResponseType; -class AliRICHParam; class AliRICHChamber : public TNamed { public: - AliRICHChamber(); - AliRICHChamber(Int_t iModuleN,AliRICHParam *pParam); + AliRICHChamber():TNamed(),fpRotMatrix(0) {;} + AliRICHChamber(Int_t iChamberN); AliRICHChamber(const AliRICHChamber &chamber):TNamed(chamber) {;} virtual ~AliRICHChamber() {;} AliRICHChamber& operator=(const AliRICHChamber&) {return *this;} - + + static Double_t AlphaFeedback(Int_t ) {return 0.02;} //determines number of feedback photons updated to 9/11/04 by Di Mauro + TRotMatrix* RotMatrix() const{return fpRotMatrix;} - const char* RotMatrixName() const{return "rot"+fName;} + TString RotMatrixName() const{return "rot"+fName;} TRotation Rot() const{return fRot;} - Double_t Rho() const{return fCenterV3.Mag();} - Double_t ThetaD() const{return fCenterV3.Theta()*TMath::RadToDeg();} - Double_t PhiD() const{return fCenterV3.Phi()*TMath::RadToDeg();} - Double_t ThetaXd() const{return fRot.ThetaX()*TMath::RadToDeg();} - Double_t PhiXd() const{return fRot.PhiX()*TMath::RadToDeg();} - Double_t ThetaYd() const{return fRot.ThetaY()*TMath::RadToDeg();} - Double_t PhiYd() const{return fRot.PhiY()*TMath::RadToDeg();} - Double_t ThetaZd() const{return fRot.ThetaZ()*TMath::RadToDeg();} - Double_t PhiZd() const{return fRot.PhiZ()*TMath::RadToDeg();} - void RotateX(Double_t a) {fRot.RotateX(a);fCenterV3.RotateX(a);fPcX3.RotateX(a);} - void RotateY(Double_t a) {fRot.RotateY(a);fCenterV3.RotateY(a);fPcX3.RotateY(a);} - void RotateZ(Double_t a) {fRot.RotateZ(a);fCenterV3.RotateZ(a);fPcX3.RotateZ(a);} - Double_t X() const{return fCenterV3.X();} - Double_t Y() const{return fCenterV3.Y();} - Double_t Z() const{return fCenterV3.Z();} - TVector3 L2G(TVector3 x3) const{x3.Transform(fRot);x3+=fCenterV3;return x3;} - TVector3 G2L(TVector3 x3) const{x3-=fCenterV3;x3.Transform(fRot.Inverse()); return x3;} - inline TVector3 Glob2Loc(TVector3 x3, Bool_t isVector=kFALSE) const; - TVector3 Glob2Loc(TLorentzVector x4,Bool_t isVector=kFALSE) const{return Glob2Loc(x4.Vect(),isVector);} - TVector3 L2G(Double_t x,Double_t y,Double_t z) const{return L2G(TVector3(x,y,z));} - TVector3 G2L(TLorentzVector x4) const{return G2L(x4.Vect());} - Float_t G2Ly(TLorentzVector x4) const{TVector3 x3=G2L(x4.Vect()); return x3.Z();} - TVector3 G2L(Double_t x,Double_t y,Double_t z) const{return G2L(TVector3(x,y,z));} - Float_t G2Lx(Double_t x,Double_t y,Double_t z) const{TVector3 x3=G2L(x,y,z); return x3.X();} - Float_t G2Ly(Double_t x,Double_t y,Double_t z) const{TVector3 x3=G2L(x,y,z); return x3.Z();} - void Print(Option_t *sOption)const;//virtual - - void LocaltoGlobal(Float_t pos[3],Float_t Localpos[3]);//Transformation from local to global coordinates, chamber-dependant - void GlobaltoLocal(Float_t pos[3],Float_t localpos[3]);//Transformation from Global to local coordinates, chamber-dependant - void DisIntegration(Float_t eloss, Float_t xhit, Float_t yhit, Int_t&x, Float_t newclust[6][500], ResponseType res);// Cluster formation method - void Init(Int_t id) {fSegmentation->Init(id);} // Recalculates all the values after some of them have been changed - void SetGeometryModel(AliRICHGeometry* pRICHGeometry) {fGeometry=pRICHGeometry;} - AliRICHGeometry* GetGeometryModel() const{return fGeometry;} - void SetResponseModel(AliRICHResponse* pRICHResponse) {fResponse=pRICHResponse;} - AliRICHResponse* GetResponseModel() const{return fResponse;} - void SetSegmentationModel(AliSegmentation* pRICHSegmentation) {fSegmentation=pRICHSegmentation;} - AliSegmentation* GetSegmentationModel() const{return fSegmentation;} - void SigGenInit(Float_t x, Float_t y, Float_t z) {fSegmentation->SigGenInit(x, y, z) ;} - Int_t SigGenCond(Float_t x, Float_t y, Float_t z) {return fSegmentation->SigGenCond(x, y, z);} - 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 - void SetPadSize(Float_t p1, Float_t p2) {fSegmentation->SetPadSize(p1,p2);} - Double_t GetX() const{return fX;} - Double_t GetY() const{return fY;} - Double_t GetZ() const{return fZ;} - inline void SetToZenith(); - TRotMatrix *GetRotMatrix() const{return fpRotMatrix;} + Double_t Rho() const{return fCenter.Mag();} //gives distance to chamber center in MRS + Double_t ThetaD() const{return fCenter.Theta()*TMath::RadToDeg();} //gives polar angle of chamber center in MRS + Double_t PhiD() const{return fCenter.Phi() *TMath::RadToDeg();} //gives azimuthal angle of chamber center in MRS + Double_t ThetaXd() const{return fRot.ThetaX() *TMath::RadToDeg();} + Double_t PhiXd() const{return fRot.PhiX() *TMath::RadToDeg();} + Double_t ThetaYd() const{return fRot.ThetaY() *TMath::RadToDeg();} + Double_t PhiYd() const{return fRot.PhiY() *TMath::RadToDeg();} + Double_t ThetaZd() const{return fRot.ThetaZ() *TMath::RadToDeg();} + Double_t PhiZd() const{return fRot.PhiZ() *TMath::RadToDeg();} + 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 + 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 + 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 + TVector3 Rad() const{return fRad;} //provides center of radiator position in MRS, cm + TVector3 Anod() const{return fAnod;} //provides center of anod wires plane in MRS, cm + TVector3 Pc() const{return fPc;} //provides center of photocathode position in MRS, cm + TVector3 Center() const{return fCenter;} //provides center of chamber position (exit from quartz window) in MRS, cm + void Print(Option_t *sOption)const; + TVector3 PMrs2Loc(TVector3 p3)const{TVector3 ploc=Rot().Invert()*p3;ploc.SetXYZ(-ploc.Px(),ploc.Py(),ploc.Pz()); return ploc;} //momentum MARS-local +//Transformations for radiator plane + 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());} + 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;} +//Transformations for anod wires plane + 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());} + 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;} +//Transformations for photcathode plane + 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());} + 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;} + TVector2 Mrs2Pc(TLorentzVector x4) const{return Mrs2Pc(x4.Vect());} protected: - Float_t fX,fY,fZ; // Position of the center of the chamber in MRS (cm) - - AliSegmentation *fSegmentation; //???Segmentation model for each chamber - AliRICHResponse *fResponse; //???Response model for each chamber - AliRICHGeometry *fGeometry; //???Geometry model for each chamber - - TVector3 fCenterV3; //chamber center position in MRS (cm) - TVector3 fPcX3; //PC center position in MRS (cm) + TVector3 fRad; //radiator entrance center position in MRS (cm) + TVector3 fCenter; //chamber center position (quartz window exit) in MRS (cm) + TVector3 fAnod; //anod wires plane center position in MRS (cm) + TVector3 fPc; //PC center position in MRS (cm) TRotation fRot; //chamber rotation in MRS TRotMatrix *fpRotMatrix; //rotation matrix of the chamber with respect to MRS - AliRICHParam *fpParam; //main RICH parameters description - ClassDef(AliRICHChamber,3) //single RICH chamber description + ClassDef(AliRICHChamber,8) //single RICH chamber description };//class AliRICHChamber -//__________________________________________________________________________________________________ -void AliRICHChamber::SetToZenith() -{ - fCenterV3.SetXYZ(fX=0,fY=AliRICHParam::Offset()-AliRICHParam::GapThickness()/2,fZ=0); - fPcX3.SetXYZ(0,AliRICHParam::Offset()-AliRICHParam::GapThickness()/2+5.276+0.25,0); -} -//__________________________________________________________________________________________________ -TVector3 AliRICHChamber::Glob2Loc(TVector3 x3,Bool_t isVector)const -{ - if(!isVector) x3-=fPcX3; - x3.Transform(fRot.Inverse()); - Double_t tmp=x3.Y(); x3.SetY(x3.Z()); x3.SetZ(tmp); - return x3; -} -//__________________________________________________________________________________________________ + #endif //AliRICHChamber_h