warning about AddDigit(Int_t*, Int_t*) fixed as Federico asked
[u/mrichter/AliRoot.git] / FMD / AliFMDRing.cxx
CommitLineData
4347b38f 1/**************************************************************************
c4bcfd14 2 * Copyright(c) 2004, ALICE Experiment at CERN, All rights reserved. *
4347b38f 3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/* $Id$ */
17
c4bcfd14 18//__________________________________________________________________
4347b38f 19//
20// Utility class to help implement collection of FMD modules into
1a1fdef7 21// rings. This is used by AliFMDDetector and AliFMDGeometry.
4347b38f 22//
1a1fdef7 23// The AliFMDGeometry object owns the AliFMDRing objects, and the
24// AliFMDDetector objects reference these. That is, the AliFMDRing
25// objects are share amoung the AliFMDDetector objects.
4347b38f 26//
27// Latest changes by Christian Holm Christensen
28//
5edc364e 29
56b1929b 30#include <AliLog.h> // ALILOG_H
e802be3e 31#include "AliFMDRing.h" // ALIFMDRING_H
56b1929b 32#include <TMath.h> // ROOT_TMath
56b1929b 33#include <TVector2.h> // ROOT_TVector2
4347b38f 34
1a1fdef7 35//====================================================================
e6a0eb08 36ClassImp(AliFMDRing)
1a1fdef7 37#if 0
38 ; // This is here to keep Emacs for indenting the next line
39#endif
4347b38f 40
41//____________________________________________________________________
1a1fdef7 42AliFMDRing::AliFMDRing(Char_t id)
43 : TNamed(Form("FMD%c", id), "Forward multiplicity ring"),
44 fId(id),
45 fVerticies(0)
4347b38f 46{
1a1fdef7 47 SetBondingWidth();
48 SetWaferRadius();
49 SetSiThickness();
50 SetLegRadius();
51 SetLegLength();
52 SetLegOffset();
53 SetModuleSpacing();
54 SetPrintboardThickness();
4ac75127 55 SetCopperThickness();
56 SetChipThickness();
4347b38f 57
1a1fdef7 58 if (fId == 'I' || fId == 'i') {
59 SetLowR(4.3);
60 SetHighR(17.2);
61 SetTheta(36/2);
62 SetNStrips(512);
63 }
64 else if (fId == 'O' || fId == 'o') {
65 SetLowR(15.6);
66 SetHighR(28.0);
67 SetTheta(18/2);
68 SetNStrips(256);
4347b38f 69 }
4347b38f 70}
71
4347b38f 72//____________________________________________________________________
73void
1a1fdef7 74AliFMDRing::Init()
4347b38f 75{
1a1fdef7 76 Double_t tanTheta = TMath::Tan(fTheta * TMath::Pi() / 180.);
77 Double_t tanTheta2 = TMath::Power(tanTheta,2);
78 Double_t r2 = TMath::Power(fWaferRadius,2);
79 Double_t yA = tanTheta * fLowR;
80 Double_t lr2 = TMath::Power(fLowR, 2);
81 Double_t hr2 = TMath::Power(fHighR,2);
82 Double_t xD = fLowR + TMath::Sqrt(r2 - tanTheta2 * lr2);
83 Double_t xD2 = TMath::Power(xD,2);
84 Double_t yB = TMath::Sqrt(r2 - hr2 + 2 * fHighR * xD - xD2);
85 Double_t xC = ((xD + TMath::Sqrt(-tanTheta2 * xD2 + r2
86 + r2 * tanTheta2))
87 / (1 + tanTheta2));
88 Double_t yC = tanTheta * xC;
4347b38f 89
1a1fdef7 90 fVerticies.Expand(6);
91 fVerticies.AddAt(new TVector2(fLowR, -yA), 0);
92 fVerticies.AddAt(new TVector2(xC, -yC), 1);
93 fVerticies.AddAt(new TVector2(fHighR, -yB), 2);
94 fVerticies.AddAt(new TVector2(fHighR, yB), 3);
95 fVerticies.AddAt(new TVector2(xC, yC), 4);
96 fVerticies.AddAt(new TVector2(fLowR, yA), 5);
97
98 fRingDepth = (fSiThickness + fPrintboardThickness
4ac75127 99 + fCopperThickness + fChipThickness
100 + fLegLength + fModuleSpacing + fSpacing);
4347b38f 101}
102
103//____________________________________________________________________
1a1fdef7 104TVector2*
105AliFMDRing::GetVertex(Int_t i) const
4347b38f 106{
1a1fdef7 107 return static_cast<TVector2*>(fVerticies.At(i));
4347b38f 108}
109
110//____________________________________________________________________
1a1fdef7 111void
112AliFMDRing::Detector2XYZ(UShort_t sector,
113 UShort_t strip,
114 Double_t& x,
115 Double_t& y,
116 Double_t& z) const
4347b38f 117{
1a1fdef7 118 if (sector >= GetNSectors()) {
119 Error("Detector2XYZ", "Invalid sector number %d (>=%d) in ring %c",
120 sector, GetNSectors(), fId);
121 return;
4347b38f 122 }
1a1fdef7 123 if (strip >= GetNStrips()) {
124 Error("Detector2XYZ", "Invalid strip number %d (>=%d)",
125 strip, GetNStrips(), fId);
126 return;
4347b38f 127 }
1a1fdef7 128 Double_t phi = Float_t(sector + .5) / GetNSectors() * 2 * TMath::Pi();
129 Double_t r = Float_t(strip + .5) / GetNStrips() * (fHighR - fLowR) + fLowR;
130 x = r * TMath::Cos(phi);
131 y = r * TMath::Sin(phi);
132 if (((sector / 2) % 2) == 1)
133 z += TMath::Sign(fModuleSpacing, z);
4347b38f 134}
135
136//
137// EOF
138//