EMCAL/AliEMCALShishKebabTrd1Module.h

   1 #ifndef ALIEMCALSHISHKEBABTRD1MODULE_H
   2 #define ALIEMCALSHISHKEBABTRD1MODULE_H
   3
   4 /* Copyright(c) 1998-2004, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 /* $Id$ */
   8 //_________________________________________________________________________
   9 // Main class for TRD1 geometry of Shish-Kebab case.
  10 // Author: Aleksei Pavlinov(WSU).
  11 // Nov 2004; Feb 2006
  12
  13 #include "TNamed.h"
  14 #include "TVector2.h"
  15
  16 class AliEMCALGeometry;
  17
  18 class AliEMCALShishKebabTrd1Module : public TNamed {
  19  public:
  20   AliEMCALShishKebabTrd1Module(Double_t theta=0.0, AliEMCALGeometry *g=0);
  21   AliEMCALShishKebabTrd1Module(AliEMCALShishKebabTrd1Module &leftNeighbor);
  22   void Init(Double_t A, Double_t B);
  23   AliEMCALShishKebabTrd1Module(const AliEMCALShishKebabTrd1Module& mod);
  24
  25   AliEMCALShishKebabTrd1Module & operator = (const AliEMCALShishKebabTrd1Module& /*rvalue*/)  {
  26     Fatal("operator =", "not implemented") ;
  27     return *this ;
  28   }
  29
  30   virtual ~AliEMCALShishKebabTrd1Module(void) {}
  31   Bool_t GetParameters();
  32   void DefineName(Double_t theta);
  33   void DefineFirstModule();
  34
  35   Double_t GetTheta() const {return fTheta;}
  36   TVector2& GetCenterOfModule() {return fOK;}
  37
  38   Double_t  GetPosX() const {return fOK.Y();}
  39   Double_t  GetPosZ() const {return fOK.X();}
  40   Double_t  GetPosXfromR() const {return fOK.Y() - fgr;}
  41   Double_t  GetA() const {return fA;}
  42   Double_t  GetB() const {return fB;}
  43   //  Additional offline staff
  44   //  ieta=0 or 1 - Jun 02, 2006
  45   TVector2& GetCenterOfCellInLocalCoordinateofSM(Int_t ieta)
  46   { if(ieta<=0) return fOK2;
  47     else        return fOK1;}
  48   void GetCenterOfCellInLocalCoordinateofSM(Int_t ieta, Double_t &xr, Double_t &zr)
  49   {
  50     if(ieta<=0) {xr = fOK2.Y(); zr = fOK2.X();
  51     } else      {xr = fOK1.Y(); zr = fOK1.X();
  52     }
  53   }
  54   // 15-may-06
  55   TVector2& GetCenterOfModuleFace() {return fOB;}
  56   TVector2& GetCenterOfModuleFace(Int_t ieta) {
  57     if(ieta<=0) return fOB2;
  58     else        return fOB1;
  59   }
  60   //
  61   Double_t GetTanBetta() const {return fgtanBetta;}
  62   Double_t Getb()        const {return fgb;}
  63   // service methods
  64   void PrintShish(Int_t pri=1) const;  // *MENU*
  65   Double_t GetThetaInDegree() const;
  66   Double_t  GetEtaOfCenterOfModule() const;
  67
  68  protected:
  69   // geometry info
  70   static AliEMCALGeometry *fgGeometry; //!
  71   static Double_t fga;        // 2*dx1=2*dy1
  72   static Double_t fga2;       // 2*dx2
  73   static Double_t fgb;        // 2*dz1
  74   static Double_t fgangle;    // ~1 degree
  75   static Double_t fgtanBetta; // tan(fgangle/2.)
  76   static Double_t fgr;        // radius to IP
  77
  78   TVector2 fOK;     // position the module center in ALICE system; x->y; z->x;
  79   Double_t fA;      // parameters of right line : y = A*z + B
  80   Double_t fB;      // system where zero point is IP.
  81   Double_t fThetaA; // angle coresponding fA - for convinience
  82   Double_t fTheta;  // theta angle of perependicular to SK module
  83   // position of towers(cells) with differents ieta (1 or 2) in local coordinate of SM
  84   // Nov 04,2004; Feb 19,2006
  85   TVector2 fOK1; // ieta=1
  86   TVector2 fOK2; // ieta=0
  87   // May 13, 2006; local position of module (cells) face
  88   TVector2 fOB;  // module
  89   TVector2 fOB1; // ieta=1
  90   TVector2 fOB2; // ieta=0
  91
  92   // public:
  93   ClassDef(AliEMCALShishKebabTrd1Module,0) // TRD1 Shish-Kebab module
  94 };
  95
  96 #endif