-#ifndef AliRICHChamber_H
-#define AliRICHChamber_H
+#ifndef AliRICHChamber_h
+#define AliRICHChamber_h
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
-/* $Id$ */
+#include <TVector3.h>
+#include <TMath.h>
+#include <TRotation.h>
+#include <TLorentzVector.h>
+#include "AliRICHParam.h"
+class TRotMatrix;
-#include <TObjArray.h>
-#include <TRotMatrix.h>
-#include "AliRICHSegRes.h"
-
-class AliRICHClusterFinder;
-class AliRICHResponse;
-class AliRICHSegmentation;
-class AliRICHGeometry;
-typedef enum {mip, cerenkov} Response_t;
-
-class AliRICHChamber : public TObject
+class AliRICHChamber : public TNamed
{
- public:
-
-//Rotation matrices for each chamber
-
- TRotMatrix *fChamberMatrix;
- Float_t fChamberTrans[3];
-
- public:
- AliRICHChamber();
- ~AliRICHChamber(){}
-//
-// Set and get GEANT id
- Int_t GetGid() {return fGid;}
- void SetGid(Int_t id) {fGid=id;}
-//
-// Initialisation and z-Position
- void Init();
- // Set inner radius of sensitive volume
- void SetRInner(Float_t rmin) {frMin=rmin;}
-// Set outer radius of sensitive volum
- void SetROuter(Float_t rmax) {frMax=rmax;}
-
-// Return inner radius of sensitive volume
- Float_t RInner() {return frMin;}
-// Return outer radius of sensitive volum
- Float_t ROuter() {return frMax;}
-
- void SetZPOS(Float_t p1) {fzPos=p1;}
- Float_t ZPosition() {return fzPos;}
-
-//
-//Transformation from Global to local coordinates, chamber-dependant
- void LocaltoGlobal(Float_t pos[3],Float_t Localpos[3]);
- void GlobaltoLocal(Float_t pos[3],Float_t localpos[3]);
-
-//Setting chamber specific rotation matrices
-
- void SetChamberTransform(Float_t Trans1,Float_t Trans2,Float_t Trans3,TRotMatrix *Matrix)
-
- {
- fChamberMatrix=Matrix;
- fChamberTrans[0]=Trans1;
- fChamberTrans[1]=Trans2;
- fChamberTrans[2]=Trans3;
- }
-
- TRotMatrix * GetRotMatrix() {return fChamberMatrix;}
-
-//Configure geometry model
- void GeometryModel(AliRICHGeometry* thisGeometry){
- fGeometry=thisGeometry;
- }
-
-
-// Configure response model
- void ResponseModel(AliRICHResponse* thisResponse);
-
- //
-// Configure segmentation model
- void SegmentationModel(AliRICHSegmentation* thisSegmentation) {
- fSegmentation = thisSegmentation;
- }
- void ReconstructionModel(AliRICHClusterFinder *thisReconstruction) {
- fReconstruction = thisReconstruction;
- }
-
-//
-// Get reference to response model
- AliRICHResponse* GetResponseModel();
-//
-// Get reference to segmentation model
- AliRICHSegmentation* GetSegmentationModel() {
- return fSegmentation;
- }
-
-// Get reference to geometry model
- AliRICHGeometry* GetGeometryModel() {
- return fGeometry;
- }
-
+public:
+ AliRICHChamber():TNamed(),fpRotMatrix(0) {;}
+ AliRICHChamber(Int_t iChamberN);
+ AliRICHChamber(const AliRICHChamber &chamber):TNamed(chamber) {;}
+ virtual ~AliRICHChamber() {;}
+ AliRICHChamber& operator=(const AliRICHChamber&) {return *this;}
+
+ TRotMatrix* RotMatrix() const{return fpRotMatrix;}
+ TString RotMatrixName() const{return "rot"+fName;}
+ TRotation Rot() const{return fRot;}
+ Double_t Rho() const{return fCenterV3.Mag();} //gives distance to chamber center in MRS
+ Double_t ThetaD() const{return fCenterV3.Theta()*TMath::RadToDeg();} //gives polar angle of chamber center in MRS
+ Double_t PhiD() const{return fCenterV3.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 RotateX(Double_t a) {fRot.RotateX(a);fCenterV3.RotateX(a);fPcX3.RotateX(a);} //rotate chamber around X by "a" degrees
+ void RotateY(Double_t a) {fRot.RotateY(a);fCenterV3.RotateY(a);fPcX3.RotateY(a);} //rotate chamber around Y by "a" degrees
+ void RotateZ(Double_t a) {fRot.RotateZ(a);fCenterV3.RotateZ(a);fPcX3.RotateZ(a);} //rotate chamber around Z by "a" degrees
+ Double_t X() const{return fCenterV3.X();}
+ Double_t Y() const{return fCenterV3.Y();}
+ Double_t Z() const{return fCenterV3.Z();}
+ TVector2 Glob2Loc(TVector3 x3)const{x3-=fPcX3;x3.Transform(fRot.Inverse());return TVector2(x3.Z()+0.5*AliRICHParam::PcSizeX(),-x3.X()+0.5*AliRICHParam::PcSizeY());}//Y and Z are misplaced?????
+ TVector3 Loc2Glob(TVector2 x2)const{TVector3 x3(-x2.Y()+0.5*AliRICHParam::PcSizeY(),0,x2.X()-0.5*AliRICHParam::PcSizeX());x3.Transform(fRot); x3+=fPcX3;return x3;}
+
+ TVector2 Glob2Loc(TLorentzVector x4) const{return Glob2Loc(x4.Vect());}
+
+ void Print(Option_t *sOption)const;//virtual
+
- AliRICHSegmentation* GetSegmentationModel(Int_t i) {
- return fSegmentation;
- }
-
- //
- AliRICHClusterFinder* &GetReconstructionModel() {return fReconstruction;}
-
- Int_t Nsec() {return fnsec;}
- void SetNsec(Int_t nsec) {fnsec=nsec;}
-//
-// Member function forwarding to the segmentation and response models
-//
-// Calculate pulse height from energy loss
- Float_t IntPH(Float_t eloss) {return fResponse->IntPH(eloss);}
- Float_t IntPH() {return fResponse->IntPH();}
-//
-// Ask segmentation if signal should be generated
- Int_t SigGenCond(Float_t x, Float_t y, Float_t z)
- {
- return fSegmentation->SigGenCond(x, y, z);
- }
-
-// Ask segmentation sector
- Int_t Sector(Float_t x, Float_t y)
- {
- return fSegmentation->Sector(x, y);
- }
-
-//
-// Initialisation of segmentation for hit
- void SigGenInit(Float_t x, Float_t y, Float_t z)
- {
- fSegmentation->SigGenInit(x, y, z) ;
- }
-// Configuration forwarding
-//
- void SetSigmaIntegration(Float_t p)
- {
- fResponse->SetSigmaIntegration(p);
- }
- void SetChargeSlope(Float_t p)
- {
- fResponse->SetChargeSlope(p);
- }
- void SetChargeSpread(Float_t p1, Float_t p2)
- {
- fResponse->SetChargeSpread(p1,p2);
- }
- void SetMaxAdc(Float_t p)
- {
- fResponse->SetMaxAdc(p);
- }
- void SetSqrtKx3(Float_t p)
- {
- fResponse->SetSqrtKx3(p);
- }
- void SetKx2(Float_t p)
- {
- fResponse->SetKx2(p);
- }
- void SetKx4(Float_t p)
- {
- fResponse->SetKx4(p);
- }
- void SetSqrtKy3(Float_t p)
- {
- fResponse->SetSqrtKy3(p);
- }
- void SetKy2(Float_t p)
- {
- fResponse->SetKy2(p);
- }
- void SetKy4(Float_t p)
- {
- fResponse->SetKy4(p);
- }
-
- void SetPitch(Float_t p)
- {
- fResponse->SetPitch(p);
- }
-
- void SetPadSize(Float_t p1, Float_t p2)
- {
- fSegmentation->SetPadSize(p1,p2);
- }
- void SetGapThickness(Float_t thickness)
- {
- fGeometry->SetGapThickness(thickness);
- }
- void SetProximityGapThickness(Float_t thickness)
- {
- fGeometry->SetProximityGapThickness(thickness);
- }
- void SetQuartzLength(Float_t length)
- {
- fGeometry->SetQuartzLength(length);
- }
- void SetQuartzWidth(Float_t width)
- {
- fGeometry->SetQuartzWidth(width);
- }
- void SetQuartzThickness(Float_t thickness)
- {
- fGeometry->SetQuartzThickness(thickness);
- }
- void SetOuterFreonLength(Float_t length)
- {
- fGeometry->SetOuterFreonLength(length);
- }
- void SetOuterFreonWidth(Float_t width)
- {
- fGeometry->SetOuterFreonWidth(width);
- }
- void SetInnerFreonLength(Float_t length)
- {
- fGeometry->SetInnerFreonLength(length);
- }
- void SetInnerFreonWidth(Float_t width)
- {
- fGeometry->SetInnerFreonWidth(width);
- }
- void SetFreonThickness(Float_t thickness)
- {
- fGeometry->SetFreonThickness(thickness);
- }
-
-
-//
-// Cluster formation method
- void DisIntegration(Float_t, Float_t, Float_t, Int_t&x, Float_t newclust[6][500], Response_t res);
- private:
-// GEANT volume if for sensitive volume of this
-// Maximum and Minimum Chamber size
- Float_t frMin;
- Float_t frMax;
- Int_t fGid;
-// z-position of this chamber
- Float_t fzPos;
-// The segmentation models for the cathode planes
-// fnsec=1: one plane segmented, fnsec=2: both planes are segmented.
- Int_t fnsec;
-
- AliRICHSegmentation *fSegmentation;
- AliRICHResponse *fResponse;
- AliRICHGeometry *fGeometry;
- AliRICHClusterFinder *fReconstruction;
- ClassDef(AliRICHChamber,1)
-};
-#endif
+ inline void SetToZenith();
+ TRotMatrix *GetRotMatrix() const{return fpRotMatrix;}
+protected:
+ TVector3 fCenterV3; //chamber center position in MRS (cm)
+ TVector3 fPcX3; //PC center position in MRS (cm)
+ TRotation fRot; //chamber rotation in MRS
+ TRotMatrix *fpRotMatrix; //rotation matrix of the chamber with respect to MRS
+ ClassDef(AliRICHChamber,6) //single RICH chamber description
+};//class AliRICHChamber
+//__________________________________________________________________________________________________
+void AliRICHChamber::SetToZenith()
+{
+//Put the chamber to zenith. Position of PC is shifted in X-Z plane since the origin of chamber local system is in
+//left hand down coner.
+ fCenterV3.SetXYZ(0,AliRICHParam::Offset()-AliRICHParam::GapThickness()/2 ,0);
+ fPcX3.SetXYZ(0,AliRICHParam::Offset()-AliRICHParam::GapThickness()/2+5.276+0.25,0);
+}
+//__________________________________________________________________________________________________
+#endif //AliRICHChamber_h
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff