SensorThickness was defined twice. Set inner chip thickness to 250mum to bypass bug...

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALShishKebabTrd1Module.h

diff --git a/EMCAL/AliEMCALShishKebabTrd1Module.h b/EMCAL/AliEMCALShishKebabTrd1Module.h
index 859ae44ba772e47197e964b48a47940555136195..f78f51ca9084749df74514c57bf180733793a41a 100644 (file)
--- a/EMCAL/AliEMCALShishKebabTrd1Module.h
+++ b/EMCAL/AliEMCALShishKebabTrd1Module.h
@@ -1,24 +1,24 @@
 #ifndef ALIEMCALSHISHKEBABTRD1MODULE_H
 #define ALIEMCALSHISHKEBABTRD1MODULE_H
 
-/* Copyright(c) 1998-2004, ALICE Experiment at CERN, All rights reserved. *
+/* Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
 /* $Id$ */
 //_________________________________________________________________________
 // Main class for TRD1 geometry of Shish-Kebab case.
-// Author: Aleksei Pavlinov(WSU).
-// Nov 2004; Feb 2006
+// Author: Alexei Pavlinov(WSU).
+// Nov 2004; Feb 2006; Apr 2010; Oct 23, 2010
 
 #include <TNamed.h>
 #include <TMath.h>
 #include <TVector2.h>
 
-class AliEMCALGeometry;
+class AliEMCALEMCGeometry;
 
 class AliEMCALShishKebabTrd1Module : public TNamed {
  public:
-  AliEMCALShishKebabTrd1Module(Double_t theta=0.0, AliEMCALGeometry *g=0);
+  AliEMCALShishKebabTrd1Module(Double_t theta=0.0, AliEMCALEMCGeometry *g=0);
   AliEMCALShishKebabTrd1Module(AliEMCALShishKebabTrd1Module &leftNeighbor);
   void Init(Double_t A, Double_t B);
   void DefineAllStaff();
@@ -32,7 +32,7 @@ class AliEMCALShishKebabTrd1Module : public TNamed {
   virtual ~AliEMCALShishKebabTrd1Module(void) {}
   Bool_t GetParameters();
   void DefineName(Double_t theta);
-  void DefineFirstModule();
+  void DefineFirstModule(const Int_t key=0); // key=0-zero tilt of first module
 
   Double_t GetTheta() const {return fTheta;}
   TVector2& GetCenterOfModule() {return fOK;}
@@ -43,27 +43,30 @@ class AliEMCALShishKebabTrd1Module : public TNamed {
   Double_t  GetA() const {return fA;}
   Double_t  GetB() const {return fB;}
   Double_t  GetRadius() const {return fgr;}
+  TVector2  GetORB() const {return fORB;}
+  TVector2  GetORT() const {return fORT;}
   //  Additional offline staff 
   //  ieta=0 or 1 - Jun 02, 2006
   TVector2& GetCenterOfCellInLocalCoordinateofSM(Int_t ieta)
-  { if(ieta<=0) return fOK2;
-    else        return fOK1;}
-  void GetCenterOfCellInLocalCoordinateofSM(Int_t ieta, Double_t &xr, Double_t &zr) const
+  { 
+    if(ieta<=0) return fOK2;
+    else        return fOK1;
+  }
+  void GetCenterOfCellInLocalCoordinateofSM(Int_t ieta, Double_t &xr, Double_t &zr, Int_t pri=0) const
   { 
     if(ieta<=0) {xr = fOK2.Y(); zr = fOK2.X();
     } else      {xr = fOK1.Y(); zr = fOK1.X();
     }
+    if(pri>0) printf(" %s ieta %2.2i xr %8.4f  zr %8.4f \n", GetName(), ieta, xr, zr);  
   }
-  void GetCenterOfCellInLocalCoordinateofSM_3X3(Int_t ieta, Double_t &xr, Double_t &zr) const
+  void GetCenterOfCellInLocalCoordinateofSM3X3(Int_t ieta, Double_t &xr, Double_t &zr) const
   { // 3X3 case - Nov 9,2006
-    ieta = ieta<0? ieta=0 : ieta; // check index
-    ieta = ieta>2? ieta=2 : ieta;
-    xr   = fOK3X3[2-ieta].Y();
-    zr   = fOK3X3[2-ieta].X();
+    if(ieta < 0) ieta = 0; //ieta = ieta<0? ieta=0 : ieta; // check index
+    if(ieta > 2) ieta = 2; //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
+  void GetCenterOfCellInLocalCoordinateofSM1X1(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();
   }
@@ -88,11 +91,11 @@ class AliEMCALShishKebabTrd1Module : public TNamed {
 
  protected:
   // geometry info
-  static AliEMCALGeometry *fgGeometry; //!
+  AliEMCALEMCGeometry *fGeometry; //!
   static Double_t fga;        // 2*dx1=2*dy1
   static Double_t fga2;       // 2*dx2
   static Double_t fgb;        // 2*dz1
-  static Double_t fgangle;    // ~1 degree
+  static Double_t fgangle;    // in rad (1.5 degree)
   static Double_t fgtanBetta; // tan(fgangle/2.)
   static Double_t fgr;        // radius to IP
 
@@ -114,8 +117,11 @@ class AliEMCALShishKebabTrd1Module : public TNamed {
   Double_t fThetaOB2;  // theta of cell center line (go through OB2)
   // 3X3 case - Nov 9,2006
   TVector2 fOK3X3[3];
+  // Apr 14, 2010 - checking of geometry
+  TVector2 fORB; // position of right/bottom point of module
+  TVector2 fORT; // position of right/top    point of module
   // public:
-  ClassDef(AliEMCALShishKebabTrd1Module,2) // TRD1 Shish-Kebab module 
+  ClassDef(AliEMCALShishKebabTrd1Module, 5) // TRD1 Shish-Kebab module 
 };
 
 #endif