TRD/AliTRDgeometry.h

   1 #ifndef ALITRDGEOMETRY_H
   2 #define ALITRDGEOMETRY_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 ///////////////////////////////////////////////////////////////////////////////
   9 //                                                                           //
  10 //  TRD geometry class                                                       //
  11 //                                                                           //
  12 ///////////////////////////////////////////////////////////////////////////////
  13
  14 #include "AliGeometry.h"
  15
  16 class AliRunLoader;
  17
  18 class AliTRDgeometry : public AliGeometry {
  19
  20  public:
  21
  22   enum { kNplan = 6, kNcham = 5, kNsect = 18, kNdet = 540, kNdets = 30 };
  23
  24   AliTRDgeometry();
  25   virtual ~AliTRDgeometry();
  26
  27   virtual void     CreateGeometry(Int_t *idtmed);
  28   virtual Int_t    IsVersion() const = 0;
  29   virtual void     Init();
  30   virtual Bool_t   Impact(const TParticle* ) const { return kTRUE; };
  31   virtual Bool_t   Local2Global(Int_t d, Double_t *local, Double_t *global) const;
  32   virtual Bool_t   Local2Global(Int_t p, Int_t c, Int_t s
  33                                 , Double_t *local, Double_t *global) const;
  34
  35   virtual Bool_t   Global2Local(Int_t mode, Double_t *local, Double_t *global
  36                                , Int_t* index) const;
  37   virtual Bool_t   Global2Detector(Double_t global[3], Int_t index[3]);
  38
  39   virtual Bool_t   Rotate(Int_t d, Double_t *pos, Double_t *rot) const;
  40   virtual Bool_t   RotateBack(Int_t d, Double_t *rot, Double_t *pos) const;
  41
  42   static  Int_t    Nsect()   { return fgkNsect; };
  43   static  Int_t    Nplan()   { return fgkNplan; };
  44   static  Int_t    Ncham()   { return fgkNcham; };
  45   static  Int_t    Ndet()    { return fgkNdet;  };
  46
  47   static  Float_t  Rmin()    { return fgkRmin;  };
  48   static  Float_t  Rmax()    { return fgkRmax;  };
  49   static  Float_t  Zmax1()   { return fgkZmax1; };
  50   static  Float_t  Zmax2()   { return fgkZmax2; };
  51
  52   static  Float_t  Cwidcha() { return (fgkSwidth2 - fgkSwidth1)
  53                              / fgkSheight * (fgkCH + fgkVspace); };
  54   static  Float_t  Cheight() { return fgkCH;      };
  55   static  Float_t  Cspace()  { return fgkVspace;  };
  56   static  Float_t  CraHght() { return fgkCraH;    };
  57   static  Float_t  CdrHght() { return fgkCdrH;    };
  58   static  Float_t  CamHght() { return fgkCamH;    };
  59   static  Float_t  CroHght() { return fgkCroH;    };
  60   static  Float_t  CroWid()  { return fgkCroW;    };
  61   static  Float_t  MyThick() { return fgkMyThick; };
  62   static  Float_t  DrThick() { return fgkDrThick; };
  63   static  Float_t  AmThick() { return fgkAmThick; };
  64   static  Float_t  DrZpos()  { return fgkDrZpos;  };
  65   static  Float_t  RpadW()   { return fgkRpadW;   };
  66   static  Float_t  CpadW()   { return fgkCpadW;   };
  67
  68   virtual void     SetPHOShole() = 0;
  69   virtual void     SetRICHhole() = 0;
  70
  71   virtual Bool_t   GetPHOShole() const = 0;
  72   virtual Bool_t   GetRICHhole() const = 0;
  73   virtual Bool_t   IsHole(Int_t /*iplan*/, Int_t /*icham*/, Int_t /*isect*/) const {return kFALSE;}
  74   static  Int_t    GetDetectorSec(Int_t p, Int_t c);
  75   static  Int_t    GetDetector(Int_t p, Int_t c, Int_t s);
  76   virtual Int_t    GetPlane(Int_t d)   const;
  77   virtual Int_t    GetChamber(Int_t d) const;
  78   virtual Int_t    GetSector(Int_t d)  const;
  79
  80           Float_t  GetChamberWidth(Int_t p) const           { return fCwidth[p];     };
  81           Float_t  GetChamberLength(Int_t p, Int_t c) const { return fClength[p][c]; };
  82
  83   virtual void     GetGlobal(const AliRecPoint* , TVector3& , TMatrixF& ) const { };
  84   virtual void     GetGlobal(const AliRecPoint* , TVector3& ) const { };
  85
  86   static  Double_t GetAlpha()  { return 2 * 3.14159265358979323846 / fgkNsect; };
  87
  88   static  AliTRDgeometry* GetGeometry(AliRunLoader* runLoader = NULL);
  89
  90   static Float_t  GetTime0(Int_t p)                        { return fgkTime0[p];          };
  91
  92  protected:
  93
  94   static const Int_t   fgkNsect;                            // Number of sectors in the full detector (18)
  95   static const Int_t   fgkNplan;                            // Number of planes of the TRD (6)
  96   static const Int_t   fgkNcham;                            // Number of chambers in z-direction (5)
  97   static const Int_t   fgkNdet;                             // Total number of detectors (18 * 6 * 5 = 540)
  98
  99   static const Float_t fgkRmin;                             // Minimal radius of the TRD
 100   static const Float_t fgkRmax;                             // Maximal radius of the TRD
 101
 102   static const Float_t fgkZmax1;                            // Half-length of the TRD at outer radius
 103   static const Float_t fgkZmax2;                            // Half-length of the TRD at inner radius
 104
 105   static const Float_t fgkSheight;                          // Height of the TRD-volume in spaceframe (BTR1-3)
 106   static const Float_t fgkSwidth1;                          // Lower width of the TRD-volume in spaceframe (BTR1-3)
 107   static const Float_t fgkSwidth2;                          // Upper width of the TRD-volume in spaceframe (BTR1-3)
 108   static const Float_t fgkSlenTR1;                          // Length of the TRD-volume in spaceframe (BTR1)
 109   static const Float_t fgkSlenTR2;                          // Length of the TRD-volume in spaceframe (BTR2)
 110   static const Float_t fgkSlenTR3;                          // Length of the TRD-volume in spaceframe (BTR3)
 111
 112   static const Float_t fgkSMpltT;                           // Thickness of the super module side plates
 113   static const Float_t fgkSMgapT;                           // Thickness of the gap between side plates and space frame
 114
 115   static const Float_t fgkCraH;                             // Height of the radiator part of the chambers
 116   static const Float_t fgkCdrH;                             // Height of the drift region of the chambers
 117   static const Float_t fgkCamH;                             // Height of the amplification region of the chambers
 118   static const Float_t fgkCroH;                             // Height of the readout of the chambers
 119   static const Float_t fgkCH;                               // Total height of the chambers
 120
 121   static const Float_t fgkVspace;                           // Vertical spacing of the chambers
 122   static const Float_t fgkHspace;                           // Horizontal spacing of the chambers
 123
 124   static const Float_t fgkCalT;                             // Thickness of the lower aluminum frame
 125   static const Float_t fgkCclsT;                            // Thickness of the lower G10 frame sides
 126   static const Float_t fgkCclfT;                            // Thickness of the lower G10 frame front
 127   static const Float_t fgkCcuT;                             // Thickness of the upper G10 frame
 128   static const Float_t fgkCauT;                             // Thickness of the upper aluminum frame
 129
 130   static const Float_t fgkCroW;                             // Additional width of the readout chamber frames
 131
 132   static const Float_t fgkCpadW;                            // Difference of outer chamber width and pad plane width
 133   static const Float_t fgkRpadW;                            // Difference of outer chamber width and pad plane width
 134
 135   static const Float_t fgkRaThick;                          // Thickness of the radiator
 136   static const Float_t fgkMyThick;                          // Thickness of the mylar-layer
 137   static const Float_t fgkXeThick;                          // Thickness of the gas volume
 138   static const Float_t fgkDrThick;                          // Thickness of the drift region
 139   static const Float_t fgkAmThick;                          // Thickness of the amplification region
 140   static const Float_t fgkCuThick;                          // Thickness of the pad plane
 141   static const Float_t fgkSuThick;                          // Thickness of the HEXCEL+G10 support structure
 142   static const Float_t fgkFeThick;                          // Thickness of the FEE + signal lines
 143   static const Float_t fgkCoThick;                          // Thickness of the PE of the cooling device
 144   static const Float_t fgkWaThick;                          // Thickness of the cooling water
 145
 146   static const Float_t fgkRaZpos;                           // Position of the radiator
 147   static const Float_t fgkMyZpos;                           // Position of the mylar-layer
 148   static const Float_t fgkDrZpos;                           // Position of the drift region
 149   static const Float_t fgkAmZpos;                           // Position of the amplification region
 150   static const Float_t fgkCuZpos;                           // Position of the pad plane
 151   static const Float_t fgkSuZpos;                           // Position of the HEXCEL+G10 support structure
 152   static const Float_t fgkFeZpos;                           // Position of the FEE + signal lines
 153   static const Float_t fgkCoZpos;                           // Position of the PE of the cooling device
 154   static const Float_t fgkWaZpos;                           // Position of the colling water
 155
 156   Float_t              fCwidth[kNplan];                     // Outer widths of the chambers
 157   Float_t              fClength[kNplan][kNcham];            // Outer lengths of the chambers
 158   Float_t              fClengthPH[kNplan][kNcham];          // For sectors with holes for the PHOS
 159   Float_t              fClengthRH[kNplan][kNcham];          // For sectors with holes for the RICH
 160
 161   Float_t              fRotA11[kNsect];                     // Matrix elements for the rotation
 162   Float_t              fRotA12[kNsect];                     // Matrix elements for the rotation
 163   Float_t              fRotA21[kNsect];                     // Matrix elements for the rotation
 164   Float_t              fRotA22[kNsect];                     // Matrix elements for the rotation
 165
 166   Float_t              fRotB11[kNsect];                     // Matrix elements for the backward rotation
 167   Float_t              fRotB12[kNsect];                     // Matrix elements for the backward rotation
 168   Float_t              fRotB21[kNsect];                     // Matrix elements for the backward rotation
 169   Float_t              fRotB22[kNsect];                     // Matrix elements for the backward rotation
 170
 171   static const Double_t fgkTime0Base;                       // Base value for calculation of Time-position of pad 0
 172   static const Float_t fgkTime0[kNplan];                    //  Time-position of pad 0
 173
 174   Float_t              fChamberUAorig[3*kNdets][3];         // Volumes origin in
 175   Float_t              fChamberUDorig[3*kNdets][3];         // the chamber
 176   Float_t              fChamberUForig[3*kNdets][3];         // [3] = x, y, z
 177   Float_t              fChamberUUorig[3*kNdets][3];
 178
 179   Float_t              fChamberUAboxd[3*kNdets][3];         // Volumes box
 180   Float_t              fChamberUDboxd[3*kNdets][3];         // dimensions (half)
 181   Float_t              fChamberUFboxd[3*kNdets][3];         // [3] = x, y, z
 182   Float_t              fChamberUUboxd[3*kNdets][3];
 183
 184   ClassDef(AliTRDgeometry,7)                                // TRD geometry base class
 185
 186 };
 187
 188 #endif