X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=TRD%2FAliTRDgeometry.h;h=13df0d0bfe6a1fca51b40810bfc5c6038cf8885f;hb=5cb60fa2d233f1f7113bba0f6c84ab9dcd307504;hp=e31d56aa3e9be7b85f1bd6739efac7912671d72e;hpb=287c5d50ef50ac0acde514fa7ea112073f65bb12;p=u%2Fmrichter%2FAliRoot.git diff --git a/TRD/AliTRDgeometry.h b/TRD/AliTRDgeometry.h index e31d56aa3e9..13df0d0bfe6 100644 --- a/TRD/AliTRDgeometry.h +++ b/TRD/AliTRDgeometry.h @@ -13,159 +13,223 @@ #include "AliGeometry.h" +#include "TObjArray.h" + +class TGeoHMatrix; + class AliRunLoader; -class AliTRDparameter; + +class AliTRDpadPlane; class AliTRDgeometry : public AliGeometry { public: - enum { kNplan = 6, kNcham = 5, kNsect = 18, kNdet = 540 }; + enum { kNplan = 6 + , kNcham = 5 + , kNsect = 18 + , kNdet = 540 + , kNdets = 30 }; AliTRDgeometry(); + AliTRDgeometry(const AliTRDgeometry &g); virtual ~AliTRDgeometry(); + AliTRDgeometry &operator=(const AliTRDgeometry &g); - virtual void CreateGeometry(Int_t *idtmed); - virtual Int_t IsVersion() const = 0; virtual void Init(); - virtual Bool_t Impact(const TParticle* ) const { return kTRUE; }; - virtual Bool_t Local2Global(Int_t d, Float_t *local, Float_t *global, AliTRDparameter *par) const; - virtual Bool_t Local2Global(Int_t p, Int_t c, Int_t s, Float_t *local, Float_t *global, AliTRDparameter *par) const; - - virtual Bool_t Global2Local(Int_t mode, Float_t *local, Float_t *global, Int_t* index, AliTRDparameter *par) const; - virtual Bool_t Global2Detector(Float_t global[3], Int_t index[3], AliTRDparameter *par); - - virtual Bool_t Rotate(Int_t d, Float_t *pos, Float_t *rot) const; - virtual Bool_t RotateBack(Int_t d, Float_t *rot, Float_t *pos) const; - - static Int_t Nsect() { return fgkNsect; }; - static Int_t Nplan() { return fgkNplan; }; - static Int_t Ncham() { return fgkNcham; }; - static Int_t Ndet() { return fgkNdet; }; - - static Float_t Rmin() { return fgkRmin; }; - static Float_t Rmax() { return fgkRmax; }; - static Float_t Zmax1() { return fgkZmax1; }; - static Float_t Zmax2() { return fgkZmax2; }; - - static Float_t Cwidcha() { return (fgkSwidth2 - fgkSwidth1) - / fgkSheight * (fgkCH + fgkVspace); }; - static Float_t Cheight() { return fgkCH; }; - static Float_t Cspace() { return fgkVspace; }; - static Float_t CraHght() { return fgkCraH; }; - static Float_t CdrHght() { return fgkCdrH; }; - static Float_t CamHght() { return fgkCamH; }; - static Float_t CroHght() { return fgkCroH; }; - static Float_t CroWid() { return fgkCroW; }; - static Float_t MyThick() { return fgkMyThick; }; - static Float_t DrThick() { return fgkDrThick; }; - static Float_t AmThick() { return fgkAmThick; }; - static Float_t DrZpos() { return fgkDrZpos; }; - static Float_t RpadW() { return fgkRpadW; }; - static Float_t CpadW() { return fgkCpadW; }; - - virtual void SetPHOShole() = 0; - virtual void SetRICHhole() = 0; - - virtual Bool_t GetPHOShole() const = 0; - virtual Bool_t GetRICHhole() const = 0; - virtual Bool_t IsHole(Int_t /*iplan*/, Int_t /*icham*/, Int_t /*isect*/) const {return kFALSE;} - virtual Int_t GetDetectorSec(Int_t p, Int_t c) const; - virtual Int_t GetDetector(Int_t p, Int_t c, Int_t s) const; - virtual Int_t GetPlane(Int_t d) const; - virtual Int_t GetChamber(Int_t d) const; - virtual Int_t GetSector(Int_t d) const; + virtual void CreateGeometry(Int_t *idtmed); + virtual Int_t IsVersion() { return 1; } + virtual Bool_t Impact(const TParticle* ) const { return kTRUE; } + virtual Bool_t IsHole(Int_t /*p*/, Int_t /*c*/, Int_t /*s*/) const { return kFALSE; } - Float_t GetChamberWidth(Int_t p) const { return fCwidth[p]; }; - Float_t GetChamberLength(Int_t p, Int_t c) const { return fClength[p][c]; }; + virtual Bool_t RotateBack(Int_t det, Double_t *loc, Double_t *glb) const; - virtual void GetGlobal(const AliRecPoint* , TVector3& , TMatrix& ) const { }; - virtual void GetGlobal(const AliRecPoint* , TVector3& ) const { }; - - static Double_t GetAlpha() { return 2 * 3.14159265358979323846 / fgkNsect; }; + Bool_t ChamberInGeometry(Int_t det); - static AliTRDgeometry* GetGeometry(AliRunLoader* runLoader = NULL); + void GroupChamber(Int_t iplan, Int_t icham, Int_t *idtmed); + void CreateFrame(Int_t *idtmed); + void CreateServices(Int_t *idtmed); + + Bool_t CreateClusterMatrixArray(); + TGeoHMatrix *GetClusterMatrix(Int_t det) { return (TGeoHMatrix *) + fClusterMatrixArray->At(det); } + + void SetSMstatus(Int_t sm, Char_t status) { fSMstatus[sm] = status; } + + static Int_t GetDetectorSec(Int_t p, Int_t c); + static Int_t GetDetector(Int_t p, Int_t c, Int_t s); + static Int_t GetPlane(Int_t d); + virtual Int_t GetChamber(Int_t d) const; + virtual Int_t GetChamber(Double_t z, Int_t plane); + virtual Int_t GetSector(Int_t d) const; + + void CreatePadPlaneArray(); + AliTRDpadPlane *CreatePadPlane(Int_t p, Int_t c); + AliTRDpadPlane *GetPadPlane(Int_t p, Int_t c); + AliTRDpadPlane *GetPadPlane(Int_t det) { return GetPadPlane(GetPlane(det) + ,GetChamber(det)); } + Int_t GetRowMax(Int_t p, Int_t c, Int_t /*s*/); + Int_t GetColMax(Int_t p); + Double_t GetRow0(Int_t p, Int_t c, Int_t /*s*/); + Double_t GetCol0(Int_t p); + + // Translation from MCM to Pad and vice versa (these functions are now in feeParam) + //virtual Int_t GetPadRowFromMCM(Int_t irob, Int_t imcm) const; + //virtual Int_t GetPadColFromADC(Int_t irob, Int_t imcm, Int_t iadc) const; + //virtual Int_t GetMCMfromPad(Int_t irow, Int_t icol) const; + //virtual Int_t GetROBfromPad(Int_t irow, Int_t icol) const; + //virtual Int_t GetRobSide(Int_t irob) const; + //virtual Int_t GetColSide(Int_t icol) const; + + static Float_t GetTime0(Int_t p) { return fgkTime0[p]; } + + Char_t GetSMstatus(Int_t sm) const { return fSMstatus[sm]; } + Float_t GetChamberWidth(Int_t p) const { return fCwidth[p]; } + Float_t GetChamberLength(Int_t p, Int_t c) const { return fClength[p][c]; } + + virtual void GetGlobal(const AliRecPoint*, TVector3&, TMatrixF& ) const { }; + virtual void GetGlobal(const AliRecPoint*, TVector3& ) const { }; + + static Double_t GetAlpha() { return 2.0 + * 3.14159265358979324 + / fgkNsect; } + + static Int_t Nsect() { return fgkNsect; } + static Int_t Nplan() { return fgkNplan; } + static Int_t Ncham() { return fgkNcham; } + static Int_t Ndet() { return fgkNdet; } + + static Float_t Cheight() { return fgkCH; } + static Float_t Cspace() { return fgkVspace; } + static Float_t CraHght() { return fgkCraH; } + static Float_t CdrHght() { return fgkCdrH; } + static Float_t CamHght() { return fgkCamH; } + static Float_t CroHght() { return fgkCroH; } + static Float_t CroWid() { return fgkCroW; } + static Float_t MyThick() { return fgkMyThick; } + static Float_t DrThick() { return fgkDrThick; } + static Float_t AmThick() { return fgkAmThick; } + static Float_t DrZpos() { return fgkDrZpos; } + static Float_t RpadW() { return fgkRpadW; } + static Float_t CpadW() { return fgkCpadW; } + + static Float_t Cwidcha() { return (fgkSwidth2 - fgkSwidth1) + / fgkSheight + * (fgkCH + fgkVspace); } + + static Int_t MCMmax() { return fgkMCMmax; } + static Int_t MCMrow() { return fgkMCMrow; } + static Int_t ROBmaxC0() { return fgkROBmaxC0; } + static Int_t ROBmaxC1() { return fgkROBmaxC1; } + static Int_t ADCmax() { return fgkADCmax; } + static Int_t TBmax() { return fgkTBmax; } + static Int_t Padmax() { return fgkPadmax; } + static Int_t Colmax() { return fgkColmax; } + static Int_t RowmaxC0() { return fgkRowmaxC0; } + static Int_t RowmaxC1() { return fgkRowmaxC1; } protected: - - static const Int_t fgkNsect; // Number of sectors in the full detector (18) - static const Int_t fgkNplan; // Number of planes of the TRD (6) - static const Int_t fgkNcham; // Number of chambers in z-direction (5) - static const Int_t fgkNdet; // Total number of detectors (18 * 6 * 5 = 540) - - static const Float_t fgkRmin; // Minimal radius of the TRD - static const Float_t fgkRmax; // Maximal radius of the TRD - - static const Float_t fgkZmax1; // Half-length of the TRD at outer radius - static const Float_t fgkZmax2; // Half-length of the TRD at inner radius - - static const Float_t fgkSheight; // Height of the TRD-volume in spaceframe (BTR1-3) - static const Float_t fgkSwidth1; // Lower width of the TRD-volume in spaceframe (BTR1-3) - static const Float_t fgkSwidth2; // Upper width of the TRD-volume in spaceframe (BTR1-3) - static const Float_t fgkSlenTR1; // Length of the TRD-volume in spaceframe (BTR1) - static const Float_t fgkSlenTR2; // Length of the TRD-volume in spaceframe (BTR2) - static const Float_t fgkSlenTR3; // Length of the TRD-volume in spaceframe (BTR3) - - static const Float_t fgkSMpltT; // Thickness of the super module side plates - static const Float_t fgkSMgapT; // Thickness of the gap between side plates and space frame - - static const Float_t fgkCraH; // Height of the radiator part of the chambers - static const Float_t fgkCdrH; // Height of the drift region of the chambers - static const Float_t fgkCamH; // Height of the amplification region of the chambers - static const Float_t fgkCroH; // Height of the readout of the chambers - static const Float_t fgkCH; // Total height of the chambers - - static const Float_t fgkVspace; // Vertical spacing of the chambers - static const Float_t fgkHspace; // Horizontal spacing of the chambers - - static const Float_t fgkCalT; // Thickness of the lower aluminum frame - static const Float_t fgkCclsT; // Thickness of the lower G10 frame sides - static const Float_t fgkCclfT; // Thickness of the lower G10 frame front - static const Float_t fgkCcuT; // Thickness of the upper G10 frame - static const Float_t fgkCauT; // Thickness of the upper aluminum frame - - static const Float_t fgkCroW; // Additional width of the readout chamber frames - - static const Float_t fgkCpadW; // Difference of outer chamber width and pad plane width - static const Float_t fgkRpadW; // Difference of outer chamber width and pad plane width - - static const Float_t fgkRaThick; // Thickness of the radiator - static const Float_t fgkMyThick; // Thickness of the mylar-layer - static const Float_t fgkXeThick; // Thickness of the gas volume - static const Float_t fgkDrThick; // Thickness of the drift region - static const Float_t fgkAmThick; // Thickness of the amplification region - static const Float_t fgkCuThick; // Thickness of the pad plane - static const Float_t fgkSuThick; // Thickness of the HEXCEL+G10 support structure - static const Float_t fgkFeThick; // Thickness of the FEE + signal lines - static const Float_t fgkCoThick; // Thickness of the PE of the cooling device - static const Float_t fgkWaThick; // Thickness of the cooling water - - static const Float_t fgkRaZpos; // Position of the radiator - static const Float_t fgkMyZpos; // Position of the mylar-layer - static const Float_t fgkDrZpos; // Position of the drift region - static const Float_t fgkAmZpos; // Position of the amplification region - static const Float_t fgkCuZpos; // Position of the pad plane - static const Float_t fgkSuZpos; // Position of the HEXCEL+G10 support structure - static const Float_t fgkFeZpos; // Position of the FEE + signal lines - static const Float_t fgkCoZpos; // Position of the PE of the cooling device - static const Float_t fgkWaZpos; // Position of the colling water - - Float_t fCwidth[kNplan]; // Outer widths of the chambers - Float_t fClength[kNplan][kNcham]; // Outer lengths of the chambers - Float_t fClengthPH[kNplan][kNcham]; // For sectors with holes for the PHOS - Float_t fClengthRH[kNplan][kNcham]; // For sectors with holes for the RICH - - Float_t fRotA11[kNsect]; // Matrix elements for the rotation - Float_t fRotA12[kNsect]; // Matrix elements for the rotation - Float_t fRotA21[kNsect]; // Matrix elements for the rotation - Float_t fRotA22[kNsect]; // Matrix elements for the rotation - - Float_t fRotB11[kNsect]; // Matrix elements for the backward rotation - Float_t fRotB12[kNsect]; // Matrix elements for the backward rotation - Float_t fRotB21[kNsect]; // Matrix elements for the backward rotation - Float_t fRotB22[kNsect]; // Matrix elements for the backward rotation - - ClassDef(AliTRDgeometry,6) // TRD geometry base class + + static const Int_t fgkNsect; // Number of sectors in the full detector (18) + static const Int_t fgkNplan; // Number of planes of the TRD (6) + static const Int_t fgkNcham; // Number of chambers in z-direction (5) + static const Int_t fgkNdet; // Total number of detectors (18 * 6 * 5 = 540) + + static const Float_t fgkTlength; // Length of the TRD-volume in spaceframe (BTRD) + + static const Float_t fgkSheight; // Height of the supermodule + static const Float_t fgkSwidth1; // Lower width of the supermodule + static const Float_t fgkSwidth2; // Upper width of the supermodule + static const Float_t fgkSlength; // Length of the supermodule + + static const Float_t fgkFlength; // Length of the service space in front of a supermodule + + static const Float_t fgkSMpltT; // Thickness of the super module side plates + + static const Float_t fgkCraH; // Height of the radiator part of the chambers + static const Float_t fgkCdrH; // Height of the drift region of the chambers + static const Float_t fgkCamH; // Height of the amplification region of the chambers + static const Float_t fgkCroH; // Height of the readout of the chambers + static const Float_t fgkCH; // Total height of the chambers + + static const Float_t fgkVspace; // Vertical spacing of the chambers + static const Float_t fgkHspace; // Horizontal spacing of the chambers + static const Float_t fgkVrocsm; // Radial distance of the first ROC to the outer SM plates + static const Float_t fgkCalT; // Thickness of the lower aluminum frame + static const Float_t fgkCalW; // Width of additional aluminum on lower frame + static const Float_t fgkCclsT; // Thickness of the lower Wacosit frame sides + static const Float_t fgkCclfT; // Thickness of the lower Wacosit frame front + static const Float_t fgkCglT; // Thichness of the glue around the radiator + static const Float_t fgkCcuT; // Thickness of the upper Wacosit frame + static const Float_t fgkCauT; // Thickness of the aluminum frame of the back panel + + static const Float_t fgkCroW; // Additional width of the readout chamber frames + + static const Float_t fgkCpadW; // Difference of outer chamber width and pad plane width + static const Float_t fgkRpadW; // Difference of outer chamber width and pad plane width + + static const Float_t fgkMyThick; // Thickness of the mylar-layer + static const Float_t fgkRaThick; // Thickness of the radiator + static const Float_t fgkXeThick; // Thickness of the gas volume + static const Float_t fgkDrThick; // Thickness of the drift region + static const Float_t fgkAmThick; // Thickness of the amplification region + static const Float_t fgkWrThick; // Thickness of the wire planes + static const Float_t fgkCuThick; // Thickness of the pad plane + static const Float_t fgkGlThick; // Thickness of the glue layer + static const Float_t fgkSuThick; // Thickness of the NOMEX support structure + static const Float_t fgkRpThick; // Thickness of the PCB readout boards + static const Float_t fgkRcThick; // Thickness of the PCB copper layers + static const Float_t fgkRoThick; // Thickness of all other ROB componentes (caps, etc.) + + static const Float_t fgkRaZpos; // Position of the radiator + static const Float_t fgkDrZpos; // Position of the drift region + static const Float_t fgkAmZpos; // Position of the amplification region + static const Float_t fgkWrZpos; // Position of the wire planes + static const Float_t fgkCuZpos; // Position of the pad plane + static const Float_t fgkGlZpos; // Position of the glue layer + static const Float_t fgkSuZpos; // Position of the HEXCEL+G10 support structure + static const Float_t fgkRpZpos; // Position of the PCB readout boards + static const Float_t fgkRcZpos; // Position of the PCB copper layers + static const Float_t fgkRoZpos; // Position of all other ROB componentes (caps, etc.) + + static const Int_t fgkMCMmax; // Maximum number of MCMs per ROB + static const Int_t fgkMCMrow; // Maximum number of MCMs per ROB Row + static const Int_t fgkROBmaxC0; // Maximum number of ROBs per C0 chamber + static const Int_t fgkROBmaxC1; // Maximum number of ROBs per C1 chamber + static const Int_t fgkADCmax; // Maximum number of ADC channels per MCM + static const Int_t fgkTBmax; // Maximum number of Time bins + static const Int_t fgkPadmax; // Maximum number of pads per MCM + static const Int_t fgkColmax; // Maximum number of pads per padplane row + static const Int_t fgkRowmaxC0; // Maximum number of Rows per C0 chamber + static const Int_t fgkRowmaxC1; // Maximum number of Rows per C1 chamber + + Char_t fSMstatus[kNsect]; // Super module status byte + + Float_t fCwidth[kNplan]; // Outer widths of the chambers + Float_t fClength[kNplan][kNcham]; // Outer lengths of the chambers + + Float_t fRotB11[kNsect]; // Matrix elements for the backward rotation + Float_t fRotB12[kNsect]; // Matrix elements for the backward rotation + Float_t fRotB21[kNsect]; // Matrix elements for the backward rotation + Float_t fRotB22[kNsect]; // Matrix elements for the backward rotation + + static const Double_t fgkTime0Base; // Base value for calculation of Time-position of pad 0 + static const Float_t fgkTime0[kNplan]; // Time-position of pad 0 + + Float_t fChamberUAorig[3*kNdets][3]; // Volumes origin in + Float_t fChamberUDorig[3*kNdets][3]; // the chamber + Float_t fChamberUForig[3*kNdets][3]; // [3] = x, y, z + Float_t fChamberUUorig[3*kNdets][3]; // + + Float_t fChamberUAboxd[3*kNdets][3]; // Volumes box + Float_t fChamberUDboxd[3*kNdets][3]; // dimensions (half) + Float_t fChamberUFboxd[3*kNdets][3]; // [3] = x, y, z + Float_t fChamberUUboxd[3*kNdets][3]; // + + TObjArray *fClusterMatrixArray; //! Transformation matrices loc. cluster to tracking cs + TObjArray *fPadPlaneArray; //! Array of pad plane objects + + ClassDef(AliTRDgeometry,16) // TRD geometry class };