X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EMCAL%2FAliEMCALShishKebabTrd1Module.h;h=f78f51ca9084749df74514c57bf180733793a41a;hb=132422328f016462e8858650098796f4d0ebfd66;hp=a97015c743e02d68a9b09be0336c099714d950d0;hpb=c63c3c5d365bcab29e17aca25ea8c04edc250e2d;p=u%2Fmrichter%2FAliRoot.git diff --git a/EMCAL/AliEMCALShishKebabTrd1Module.h b/EMCAL/AliEMCALShishKebabTrd1Module.h index a97015c743e..f78f51ca908 100644 --- a/EMCAL/AliEMCALShishKebabTrd1Module.h +++ b/EMCAL/AliEMCALShishKebabTrd1Module.h @@ -1,75 +1,127 @@ #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$ */ -//*-- 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: Alexei Pavlinov(WSU). +// Nov 2004; Feb 2006; Apr 2010; Oct 23, 2010 -#include "TNamed.h" -#include "TMath.h" -#include "TVector2.h" +#include +#include +#include -class AliEMCALGeometry; +class AliEMCALEMCGeometry; class AliEMCALShishKebabTrd1Module : public TNamed { public: - AliEMCALShishKebabTrd1Module(double theta=TMath::Pi()/2.); + AliEMCALShishKebabTrd1Module(Double_t theta=0.0, AliEMCALEMCGeometry *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 DefineFirstModule(); + void DefineName(Double_t theta); + void DefineFirstModule(const Int_t key=0); // key=0-zero tilt of first module - 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;} + TVector2 GetORB() const {return fORB;} + TVector2 GetORT() const {return fORT;} // Additional offline staff - TVector2& GetCenterOfCell(Int_t ieta) - { if(ieta<=1) return fOK1; - else return fOK2;} + // 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, 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 GetCenterOfCellInLocalCoordinateofSM3X3(Int_t ieta, Double_t &xr, Double_t &zr) const + { // 3X3 case - Nov 9,2006 + 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 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(); + } + // 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; //! + 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.) - // 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 - // position of towers with differents ieta (1 or 2) - 4-nov-04 - TVector2 fOK1; - TVector2 fOK2; - - public: - ClassDef(AliEMCALShishKebabTrd1Module,0) // Turned Shish-Kebab 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=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]; + // 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, 5) // TRD1 Shish-Kebab module }; #endif -/* To do - 1. Insert position the center of towers - 2 additional TVector2 - */