X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=TRD%2FAliTRDgeometry.h;h=c73001dc7e758cb3b72fee9ed8f7dbc873ede414;hb=33d4a5802f751ff1cacd0542234357662333b13f;hp=f5fe2417e27f2a1e12a0090ee08fe967e2856159;hpb=b640260adf3ef0182cf212d83539288264997595;p=u%2Fmrichter%2FAliRoot.git diff --git a/TRD/AliTRDgeometry.h b/TRD/AliTRDgeometry.h index f5fe2417e27..c73001dc7e7 100644 --- a/TRD/AliTRDgeometry.h +++ b/TRD/AliTRDgeometry.h @@ -11,13 +11,11 @@ // // /////////////////////////////////////////////////////////////////////////////// -#include "AliGeometry.h" - #include "TObjArray.h" -class TGeoHMatrix; +#include "AliGeometry.h" -class AliRunLoader; +class TGeoHMatrix; class AliTRDpadPlane; @@ -25,11 +23,11 @@ class AliTRDgeometry : public AliGeometry { public: - enum { kNplan = 6 - , kNcham = 5 - , kNsect = 18 - , kNdet = 540 - , kNdets = 30 }; + enum { kNlayer = 6 + , kNstack = 5 + , kNsector = 18 + , kNdet = 540 + , kNdets = 30 }; AliTRDgeometry(); AliTRDgeometry(const AliTRDgeometry &g); @@ -38,96 +36,101 @@ class AliTRDgeometry : public AliGeometry { virtual void Init(); 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; } - + virtual Int_t IsVersion() { return 1; } + virtual Bool_t Impact(const TParticle* ) const { return kTRUE; } + virtual Bool_t IsHole(Int_t la, Int_t st, Int_t se) const; + virtual Bool_t IsOnBoundary(Int_t det, Float_t y, Float_t z, Float_t eps = .5) const; virtual Bool_t RotateBack(Int_t det, Double_t *loc, Double_t *glb) const; - void GroupChamber(Int_t iplan, Int_t icham, Int_t *idtmed); + Bool_t ChamberInGeometry(Int_t det); + + void AssembleChamber(Int_t ilayer, Int_t istack); void CreateFrame(Int_t *idtmed); void CreateServices(Int_t *idtmed); Bool_t CreateClusterMatrixArray(); - TGeoHMatrix *GetClusterMatrix(Int_t det) { return (TGeoHMatrix *) - fClusterMatrixArray->At(det); } + TGeoHMatrix *GetClusterMatrix(Int_t det) { return (TGeoHMatrix *) + fClusterMatrixArray->At(det); } - void SetSMstatus(Int_t sm, Char_t status) { fSMstatus[sm] = status; } + 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 GetSector(Int_t d) const; + static Int_t GetDetectorSec(Int_t layer, Int_t stack); + static Int_t GetDetector(Int_t layer, Int_t stack, Int_t sector); + static Int_t GetLayer(Int_t det); + static Int_t GetStack(Int_t det); + Int_t GetStack(Double_t z, Int_t layer); + static Int_t GetSector(Int_t det); void CreatePadPlaneArray(); - AliTRDpadPlane *CreatePadPlane(Int_t p, Int_t c); - AliTRDpadPlane *GetPadPlane(Int_t p, Int_t c); - 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]; } + AliTRDpadPlane *CreatePadPlane(Int_t layer, Int_t stack); + AliTRDpadPlane *GetPadPlane(Int_t layer, Int_t stack); + AliTRDpadPlane *GetPadPlane(Int_t det) { return GetPadPlane(GetLayer(det) + ,GetStack(det)); } + Int_t GetRowMax(Int_t layer, Int_t stack, Int_t /*sector*/); + Int_t GetColMax(Int_t layer); + Double_t GetRow0(Int_t layer, Int_t stack, Int_t /*sector*/); + Double_t GetCol0(Int_t layer); + + static Float_t GetTime0(Int_t layer) { return fgkTime0[layer]; } + + static Double_t GetXtrdBeg() { return fgkXtrdBeg; } + static Double_t GetXtrdEnd() { return fgkXtrdEnd; } + + Char_t GetSMstatus(Int_t sm) const { return fSMstatus[sm]; } + Float_t GetChamberWidth(Int_t layer) const { return fCwidth[layer] ; } + Float_t GetChamberLength(Int_t layer, Int_t stack) const { return fClength[layer][stack]; } 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; } + static Double_t GetAlpha() { return 2.0 + * 3.14159265358979324 + / fgkNsector; } + + static Int_t Nsector() { return fgkNsector; } + static Int_t Nlayer() { return fgkNlayer; } + static Int_t Nstack() { return fgkNstack; } + static Int_t Ndet() { return fgkNdet; } + + static Float_t Cheight() { return fgkCH; } + static Float_t CheightSV() { return fgkCHsv; } + 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 CsvHght() { return fgkCsvH; } + static Float_t CroWid() { return fgkCroW; } + + static Float_t AnodePos() { return fgkAnodePos; } + + static Float_t MyThick() { return fgkRMyThick; } + 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 fgkNsector; // Number of sectors in the full detector (18) + static const Int_t fgkNlayer; // Number of layers of the TRD (6) + static const Int_t fgkNstack; // Number of stacks 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) @@ -145,47 +148,59 @@ class AliTRDgeometry : public AliGeometry { 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 fgkCsvH; // Height of the services on top of the chambers + static const Float_t fgkCH; // Total height of the chambers (w/o services) + static const Float_t fgkCHsv; // Total height of the chambers (with services) + + static const Float_t fgkAnodePos; // Distance of anode wire plane relative to alignabl volume 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 fgkCalW; // Width of additional aluminum ledge on lower frame + static const Float_t fgkCalH; // Height of additional aluminum ledge on lower frame + static const Float_t fgkCalWmod; // Width of additional aluminum ledge on lower frame + static const Float_t fgkCalHmod; // Height of additional aluminum ledge on lower frame + static const Float_t fgkCwsW; // Width of additional wacosit ledge on lower frame + static const Float_t fgkCwsH; // Height of additional wacosit ledge 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 fgkCcuTa; // Thickness of the upper Wacosit frame around amp. region + static const Float_t fgkCcuTb; // Thickness of the upper Wacosit frame around amp. region 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 fgkPPdThick; // Thickness of copper of the pad plane + static const Float_t fgkPPpThick; // Thickness of PCB board of the pad plane + static const Float_t fgkPGlThick; // Thickness of the glue layer + static const Float_t fgkPCbThick; // Thickness of the carbon layers + static const Float_t fgkPHcThick; // Thickness of the honeycomb support structure + static const Float_t fgkPPcThick; // Thickness of the PCB readout boards + static const Float_t fgkPRbThick; // Thickness of the PCB copper layers + static const Float_t fgkPElThick; // Thickness of all other electronics components (caps, etc.) + + static const Float_t fgkRFbThick; // Thickness of the fiber layers in the radiator + static const Float_t fgkRRhThick; // Thickness of the rohacell layers in the radiator + static const Float_t fgkRGlThick; // Thickness of the glue layers in the radiator + static const Float_t fgkRCbThick; // Thickness of the carbon layers in the radiator + static const Float_t fgkRMyThick; // Thickness of the mylar layers in 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 Float_t fgkWrZposA; // Position of the wire planes + static const Float_t fgkWrZposB; // Position of the wire planes + static const Float_t fgkCalZpos; // Position of the additional aluminum ledges static const Int_t fgkMCMmax; // Maximum number of MCMs per ROB static const Int_t fgkMCMrow; // Maximum number of MCMs per ROB Row @@ -198,34 +213,26 @@ class AliTRDgeometry : public AliGeometry { 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[kNlayer]; // Outer widths of the chambers + Float_t fClength[kNlayer][kNstack]; // Outer lengths of the chambers - 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 + Float_t fRotB11[kNsector]; // Matrix elements for the backward rotation + Float_t fRotB12[kNsector]; // Matrix elements for the backward rotation + Float_t fRotB21[kNsector]; // Matrix elements for the backward rotation + Float_t fRotB22[kNsector]; // 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]; // + static const Float_t fgkTime0[kNlayer]; // Time-position of pad 0 + + static const Double_t fgkXtrdBeg; // X-coordinate in tracking system of begin of TRD mother volume + static const Double_t fgkXtrdEnd; // X-coordinate in tracking system of end of TRD mother volume TObjArray *fClusterMatrixArray; //! Transformation matrices loc. cluster to tracking cs TObjArray *fPadPlaneArray; //! Array of pad plane objects - ClassDef(AliTRDgeometry,16) // TRD geometry class + Char_t fSMstatus[kNsector]; // Super module status byte -}; + ClassDef(AliTRDgeometry,23) // TRD geometry class +}; #endif