renamed a method to load libs so that it does not clash with earlier runs
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALShishKebabTrd1Module.h
index 7b96e7f..36af215 100644 (file)
@@ -5,52 +5,90 @@
  * See cxx source for full Copyright notice                               */
 
 /* $Id$ */
-//*-- Author: Aleksei Pavlinov (WSU)
-// TO DO : create class for Super module Geometry - 4-nov-04 
+//_________________________________________________________________________
+// Main class for TRD1 geometry of Shish-Kebab case.
+// Author: Aleksei Pavlinov(WSU).
+// Nov 2004; Feb 2006
 
-#include "TNamed.h"
-#include "TMath.h"
-#include "TVector2.h"
+#include <TNamed.h>
+#include <TMath.h>
+#include <TVector2.h>
 
 class AliEMCALGeometry;
 
 class AliEMCALShishKebabTrd1Module : public TNamed {
  public:
-  AliEMCALShishKebabTrd1Module(double theta=TMath::Pi()/2., AliEMCALGeometry *g=0);
+  AliEMCALShishKebabTrd1Module(Double_t theta=0.0, AliEMCALGeometry *g=0);
   AliEMCALShishKebabTrd1Module(AliEMCALShishKebabTrd1Module &leftNeighbor);
-  void Init(double A, double B);
+  void Init(Double_t A, Double_t B);
+  void DefineAllStaff();
+  AliEMCALShishKebabTrd1Module(const AliEMCALShishKebabTrd1Module& mod);
+
+  AliEMCALShishKebabTrd1Module & operator = (const AliEMCALShishKebabTrd1Module& /*rvalue*/)  {
+    Fatal("operator =", "not implemented") ;  
+    return *this ; 
+  }
 
   virtual ~AliEMCALShishKebabTrd1Module(void) {}
   Bool_t GetParameters();
-  void DefineName(double theta);
+  void DefineName(Double_t theta);
   void DefineFirstModule();
 
-  Double_t GetTheta() const{return fTheta;}
-  Double_t GetThetaInDegree() const {return fTheta*180./TMath::Pi();}
+  Double_t GetTheta() const {return fTheta;}
   TVector2& GetCenterOfModule() {return fOK;}
-  Double_t  GetEtaOfCenterOfModule(){return -TMath::Log(TMath::Tan(fOK.Phi()/2.));}
 
-  Double_t  GetPosX() {return fOK.Y();}
-  Double_t  GetPosZ() {return fOK.X();}
-  Double_t  GetPosXfromR() {return fOK.Y() - fgr;}
-  Double_t  GetA() {return fA;}
-  Double_t  GetB() {return fB;}
+  Double_t  GetPosX() const {return fOK.Y();}
+  Double_t  GetPosZ() const {return fOK.X();}
+  Double_t  GetPosXfromR() const {return fOK.Y() - fgr;}
+  Double_t  GetA() const {return fA;}
+  Double_t  GetB() const {return fB;}
+  Double_t  GetRadius() const {return fgr;}
   //  Additional offline staff 
+  //  ieta=0 or 1 - Jun 02, 2006
   TVector2& GetCenterOfCellInLocalCoordinateofSM(Int_t ieta)
-  { if(ieta<=1) return fOK2;
+  { if(ieta<=0) return fOK2;
     else        return fOK1;}
-  void GetCenterOfCellInLocalCoordinateofSM(Int_t ieta, Double_t &xr, Double_t &zr)
+  void GetCenterOfCellInLocalCoordinateofSM(Int_t ieta, Double_t &xr, Double_t &zr) const
   { 
-    if(ieta<=1) {xr = fOK2.Y(); zr = fOK2.X();
+    if(ieta<=0) {xr = fOK2.Y(); zr = fOK2.X();
     } else      {xr = fOK1.Y(); zr = fOK1.X();
     }
   }
+  void GetCenterOfCellInLocalCoordinateofSM_3X3(Int_t ieta, Double_t &xr, Double_t &zr) const
+  { // 3X3 case - Nov 9,2006
+    if(ieta < 0) ieta = 0;
+    if(ieta > 2) ieta = 2;
+    //    ieta = ieta<0? ieta=0 : ieta; // check index
+    //ieta = ieta>2? ieta=2 : ieta;
+    xr   = fOK3X3[2-ieta].Y();
+    zr   = fOK3X3[2-ieta].X();
+  }
+  void GetCenterOfCellInLocalCoordinateofSM_1X1(Double_t &xr, Double_t &zr) const
+  { // 1X1 case - Nov 27,2006
+    // Center of cell is center of module
+    xr   = fOK.Y() - fgr;
+    zr   = fOK.X();
+  }
+  // 15-may-06
+  TVector2& GetCenterOfModuleFace() {return fOB;}  
+  TVector2& GetCenterOfModuleFace(Int_t ieta) {
+    if(ieta<=0) return fOB2;
+    else        return fOB1;
+  }  
+  // Jul 30, 2007
+  void GetPositionAtCenterCellLine(Int_t ieta, Double_t dist, TVector2 &v);
   // 
-  Double_t GetTanBetta() {return fgtanBetta;}
-  Double_t Getb()        {return fgb;}
+  Double_t GetTanBetta() const {return fgtanBetta;}
+  Double_t Getb()        const {return fgb;}
   // service methods
-  void PrintShish(int pri=1) const;  // *MENU*
+  void PrintShish(Int_t pri=1) const;  // *MENU*
+  Double_t  GetThetaInDegree() const;
+  Double_t  GetEtaOfCenterOfModule() const;
+  Double_t  GetMaxEtaOfModule(int pri=0) const;
+  static Double_t ThetaToEta(Double_t theta) 
+  {return -TMath::Log(TMath::Tan(theta/2.));}
 
+ protected:
   // geometry info
   static AliEMCALGeometry *fgGeometry; //!
   static Double_t fga;        // 2*dx1=2*dy1
@@ -58,22 +96,28 @@ class AliEMCALShishKebabTrd1Module : public TNamed {
   static Double_t fgb;        // 2*dz1
   static Double_t fgangle;    // ~1 degree
   static Double_t fgtanBetta; // tan(fgangle/2.)
-  // radius to IP
-  static Double_t fgr;
+  static Double_t fgr;        // radius to IP
 
- protected:
-  TVector2 fOK;     // position the module center x->y; z->x;
+  TVector2 fOK;     // position the module center in ALICE system; x->y; z->x;
   Double_t fA;      // parameters of right line : y = A*z + B
   Double_t fB;      // system where zero point is IP.
   Double_t fThetaA; // angle coresponding fA - for convinience
-  Double_t fTheta;  // theta angle of perependicular to SK module
+  Double_t fTheta;  // theta angle of perpendicular to SK module
   // position of towers(cells) with differents ieta (1 or 2) in local coordinate of SM
   // Nov 04,2004; Feb 19,2006 
-  TVector2 fOK1; // ieta=2
-  TVector2 fOK2; // ieta=1
-
- public:
-  ClassDef(AliEMCALShishKebabTrd1Module,0) // Turned Shish-Kebab module 
+  TVector2 fOK1; // ieta=1
+  TVector2 fOK2; // ieta=0
+  // May 13, 2006; local position of module (cells) center face  
+  TVector2 fOB;   // module
+  TVector2 fOB1;  // ieta=1
+  TVector2 fOB2;  // ieta=0
+  // Jul 30, 2007
+  Double_t fThetaOB1; // theta of cell center line (go through OB1)
+  Double_t fThetaOB2;  // theta of cell center line (go through OB2)
+  // 3X3 case - Nov 9,2006
+  TVector2 fOK3X3[3];
+  // public:
+  ClassDef(AliEMCALShishKebabTrd1Module,2) // TRD1 Shish-Kebab module 
 };
 
 #endif