X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=TRD%2FAliTRDgeometry.h;h=b08ca1f786f0c421f9e780c72cfe8b7dce49c4bd;hp=fb91195c97739c10679ff829d197a84508b0b0cf;hb=3d7b6a24571685d335860346c2b4d83f03b0d329;hpb=8230f242a54d81b01f5a6fe9c900cee395cdaf57 diff --git a/TRD/AliTRDgeometry.h b/TRD/AliTRDgeometry.h index fb91195c977..b08ca1f786f 100644 --- a/TRD/AliTRDgeometry.h +++ b/TRD/AliTRDgeometry.h @@ -1,80 +1,169 @@ -#ifndef ALITRDGEOMETRY_h -#define ALITRDGEOMETRY_h +#ifndef ALITRDGEOMETRY_H +#define ALITRDGEOMETRY_H /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * See cxx source for full Copyright notice */ /* $Id$ */ -#include -#include +/////////////////////////////////////////////////////////////////////////////// +// // +// TRD geometry class // +// // +/////////////////////////////////////////////////////////////////////////////// -#include "AliRun.h" -#include "AliRecPoint.h" +#include "AliGeometry.h" -#include "AliTRDconst.h" +class AliTRDparameter; class AliTRDgeometry : public AliGeometry { public: + enum { kNplan = 6, kNcham = 5, kNsect = 18, kNdet = 540 }; + AliTRDgeometry(); virtual ~AliTRDgeometry(); - virtual void CreateGeometry(Int_t *idtmed); - virtual Int_t IsVersion() const = 0; - virtual void Init(); - virtual Bool_t Local2Global(Int_t d, Float_t *local, Float_t *global); - virtual Bool_t Local2Global(Int_t p, Int_t c, Int_t s, Float_t *local, Float_t *global); - virtual Bool_t Rotate(Int_t d, Float_t *pos, Float_t *rot); - virtual Bool_t RotateBack(Int_t d, Float_t *rot, Float_t *pos); - - virtual void SetPHOShole() = 0; - virtual void SetRICHhole() = 0; - - virtual void SetRowPadSize(Float_t size) { fRowPadSize = size; }; - virtual void SetColPadSize(Float_t size) { fColPadSize = size; }; - virtual void SetTimeBinSize(Float_t size) { fTimeBinSize = size; }; - - virtual Bool_t GetPHOShole() = 0; - virtual Bool_t GetRICHhole() = 0; - - virtual Int_t GetDetector(Int_t p, Int_t c, Int_t s); - virtual Int_t GetPlane(Int_t d); - virtual Int_t GetChamber(Int_t d); - virtual Int_t GetSector(Int_t d); - - virtual Int_t GetRowMax(Int_t p, Int_t c, Int_t s) { return fRowMax[p][c][s]; }; - virtual Int_t GetColMax(Int_t p) { return fColMax[p]; }; - virtual Int_t GetTimeMax() { return fTimeMax; }; + 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 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 void SetOldGeometry(); + + virtual Bool_t GetPHOShole() const = 0; + virtual Bool_t GetRICHhole() const = 0; + + virtual Int_t GetDetectorSec(const Int_t p, const Int_t c) const; + virtual Int_t GetDetector(const Int_t p, const Int_t c, const Int_t s) const; + virtual Int_t GetPlane(const Int_t d) const; + virtual Int_t GetChamber(const Int_t d) const; + virtual Int_t GetSector(const Int_t d) const; + + Float_t GetChamberWidth(const Int_t p) const { return fCwidth[p]; }; + Float_t GetChamberLength(const Int_t p, const Int_t c) const { return fClength[p][c]; }; + + virtual void GetGlobal(const AliRecPoint* , TVector3& , TMatrix& ) const { }; + virtual void GetGlobal(const AliRecPoint* , TVector3& ) const { }; - virtual Float_t GetRow0(Int_t p, Int_t c, Int_t s) { return fRow0[p][c][s]; }; - virtual Float_t GetCol0(Int_t p) { return fCol0[p]; }; - virtual Float_t GetTime0(Int_t p) { return fTime0[p]; }; - - virtual Float_t GetRowPadSize() { return fRowPadSize; }; - virtual Float_t GetColPadSize() { return fColPadSize; }; - virtual Float_t GetTimeBinSize() { return fTimeBinSize; }; - - virtual void GetGlobal(const AliRecPoint *p, TVector3 &pos, TMatrix &mat); - virtual void GetGlobal(const AliRecPoint *p, TVector3 &pos); + static Double_t GetAlpha() { return 2 * 3.14159265358979323846 / fgkNsect; }; protected: - - Float_t fCwidth[kNplan]; // Width of the chambers - - Int_t fRowMax[kNplan][kNcham][kNsect]; // Number of pad-rows - Int_t fColMax[kNplan]; // Number of pad-columns - Int_t fTimeMax; // Number of time buckets - - Float_t fRow0[kNplan][kNcham][kNsect]; // Row-position of pad 0 - Float_t fCol0[kNplan]; // Column-position of pad 0 - Float_t fTime0[kNplan]; // Time-position of pad 0 - - Float_t fRowPadSize; // Pad size in z-direction - Float_t fColPadSize; // Pad size in rphi-direction - Float_t fTimeBinSize; // Size of the time buckets - - ClassDef(AliTRDgeometry,1) // 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 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 };