]>
Commit | Line | Data |
---|---|---|
1963b290 | 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$ */ | |
25ab8a31 | 8 | //_________________________________________________________________________ |
9 | // Main class for TRD1 geometry of Shish-Kebab case. | |
10 | // Author: Aleksei Pavlinov(WSU). | |
11 | // Nov 2004; Feb 2006 | |
1963b290 | 12 | |
9edefa04 | 13 | #include <TNamed.h> |
14 | #include <TMath.h> | |
15 | #include <TVector2.h> | |
1963b290 | 16 | |
0c5b726e | 17 | class AliEMCALEMCGeometry; |
1963b290 | 18 | |
19 | class AliEMCALShishKebabTrd1Module : public TNamed { | |
20 | public: | |
0c5b726e | 21 | AliEMCALShishKebabTrd1Module(Double_t theta=0.0, AliEMCALEMCGeometry *g=0); |
1963b290 | 22 | AliEMCALShishKebabTrd1Module(AliEMCALShishKebabTrd1Module &leftNeighbor); |
25ab8a31 | 23 | void Init(Double_t A, Double_t B); |
1d46d1f6 | 24 | void DefineAllStaff(); |
0a4cb131 | 25 | AliEMCALShishKebabTrd1Module(const AliEMCALShishKebabTrd1Module& mod); |
26 | ||
25ab8a31 | 27 | AliEMCALShishKebabTrd1Module & operator = (const AliEMCALShishKebabTrd1Module& /*rvalue*/) { |
28 | Fatal("operator =", "not implemented") ; | |
29 | return *this ; | |
30 | } | |
1963b290 | 31 | |
32 | virtual ~AliEMCALShishKebabTrd1Module(void) {} | |
33 | Bool_t GetParameters(); | |
25ab8a31 | 34 | void DefineName(Double_t theta); |
1963b290 | 35 | void DefineFirstModule(); |
36 | ||
25ab8a31 | 37 | Double_t GetTheta() const {return fTheta;} |
1963b290 | 38 | TVector2& GetCenterOfModule() {return fOK;} |
1963b290 | 39 | |
25ab8a31 | 40 | Double_t GetPosX() const {return fOK.Y();} |
41 | Double_t GetPosZ() const {return fOK.X();} | |
42 | Double_t GetPosXfromR() const {return fOK.Y() - fgr;} | |
43 | Double_t GetA() const {return fA;} | |
44 | Double_t GetB() const {return fB;} | |
1d46d1f6 | 45 | Double_t GetRadius() const {return fgr;} |
d1ccaddd | 46 | TVector2 GetORB() const {return fORB;} |
47 | TVector2 GetORT() const {return fORT;} | |
1963b290 | 48 | // Additional offline staff |
dc7da436 | 49 | // ieta=0 or 1 - Jun 02, 2006 |
e52475ed | 50 | TVector2& GetCenterOfCellInLocalCoordinateofSM(Int_t ieta) |
dc7da436 | 51 | { if(ieta<=0) return fOK2; |
e52475ed | 52 | else return fOK1;} |
1d46d1f6 | 53 | void GetCenterOfCellInLocalCoordinateofSM(Int_t ieta, Double_t &xr, Double_t &zr) const |
e52475ed | 54 | { |
dc7da436 | 55 | if(ieta<=0) {xr = fOK2.Y(); zr = fOK2.X(); |
e52475ed | 56 | } else {xr = fOK1.Y(); zr = fOK1.X(); |
57 | } | |
58 | } | |
0c5b726e | 59 | void GetCenterOfCellInLocalCoordinateofSM3X3(Int_t ieta, Double_t &xr, Double_t &zr) const |
1d46d1f6 | 60 | { // 3X3 case - Nov 9,2006 |
0c5b726e | 61 | if(ieta < 0) ieta = 0; //ieta = ieta<0? ieta=0 : ieta; // check index |
62 | if(ieta > 2) ieta = 2; //ieta = ieta>2? ieta=2 : ieta; | |
63 | xr = fOK3X3[2-ieta].Y(); zr = fOK3X3[2-ieta].X(); | |
1d46d1f6 | 64 | } |
0c5b726e | 65 | void GetCenterOfCellInLocalCoordinateofSM1X1(Double_t &xr, Double_t &zr) const |
66 | { // 1X1 case - Nov 27,2006 // Center of cell is center of module | |
d25f2c54 | 67 | xr = fOK.Y() - fgr; |
68 | zr = fOK.X(); | |
69 | } | |
b44d5aa4 | 70 | // 15-may-06 |
71 | TVector2& GetCenterOfModuleFace() {return fOB;} | |
72 | TVector2& GetCenterOfModuleFace(Int_t ieta) { | |
dc7da436 | 73 | if(ieta<=0) return fOB2; |
b44d5aa4 | 74 | else return fOB1; |
75 | } | |
1ae500a2 | 76 | // Jul 30, 2007 |
77 | void GetPositionAtCenterCellLine(Int_t ieta, Double_t dist, TVector2 &v); | |
1963b290 | 78 | // |
25ab8a31 | 79 | Double_t GetTanBetta() const {return fgtanBetta;} |
80 | Double_t Getb() const {return fgb;} | |
1963b290 | 81 | // service methods |
25ab8a31 | 82 | void PrintShish(Int_t pri=1) const; // *MENU* |
1d46d1f6 | 83 | Double_t GetThetaInDegree() const; |
25ab8a31 | 84 | Double_t GetEtaOfCenterOfModule() const; |
1d46d1f6 | 85 | Double_t GetMaxEtaOfModule(int pri=0) const; |
86 | static Double_t ThetaToEta(Double_t theta) | |
87 | {return -TMath::Log(TMath::Tan(theta/2.));} | |
1963b290 | 88 | |
25ab8a31 | 89 | protected: |
1963b290 | 90 | // geometry info |
0c5b726e | 91 | static AliEMCALEMCGeometry *fgGeometry; //! |
1963b290 | 92 | static Double_t fga; // 2*dx1=2*dy1 |
93 | static Double_t fga2; // 2*dx2 | |
94 | static Double_t fgb; // 2*dz1 | |
95 | static Double_t fgangle; // ~1 degree | |
96 | static Double_t fgtanBetta; // tan(fgangle/2.) | |
25ab8a31 | 97 | static Double_t fgr; // radius to IP |
1963b290 | 98 | |
b44d5aa4 | 99 | TVector2 fOK; // position the module center in ALICE system; x->y; z->x; |
1963b290 | 100 | Double_t fA; // parameters of right line : y = A*z + B |
101 | Double_t fB; // system where zero point is IP. | |
102 | Double_t fThetaA; // angle coresponding fA - for convinience | |
1d46d1f6 | 103 | Double_t fTheta; // theta angle of perpendicular to SK module |
e52475ed | 104 | // position of towers(cells) with differents ieta (1 or 2) in local coordinate of SM |
105 | // Nov 04,2004; Feb 19,2006 | |
dc7da436 | 106 | TVector2 fOK1; // ieta=1 |
107 | TVector2 fOK2; // ieta=0 | |
1ae500a2 | 108 | // May 13, 2006; local position of module (cells) center face |
109 | TVector2 fOB; // module | |
110 | TVector2 fOB1; // ieta=1 | |
111 | TVector2 fOB2; // ieta=0 | |
112 | // Jul 30, 2007 | |
113 | Double_t fThetaOB1; // theta of cell center line (go through OB1) | |
114 | Double_t fThetaOB2; // theta of cell center line (go through OB2) | |
1d46d1f6 | 115 | // 3X3 case - Nov 9,2006 |
116 | TVector2 fOK3X3[3]; | |
d1ccaddd | 117 | // Apr 14, 2010 - checking of geometry |
118 | TVector2 fORB; // position of right/bottom point of module | |
119 | TVector2 fORT; // position of right/top point of module | |
25ab8a31 | 120 | // public: |
d1ccaddd | 121 | ClassDef(AliEMCALShishKebabTrd1Module, 4) // TRD1 Shish-Kebab module |
1963b290 | 122 | }; |
123 | ||
124 | #endif |