Fixing the correct syntax for AliEn spooler
[u/mrichter/AliRoot.git] / FMD / AliFMD3.cxx
... / ...
CommitLineData
1/**************************************************************************
2 * Copyright(c) 2004, ALICE Experiment at CERN, All rights reserved. *
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/* $Id$ */
16/** @file AliFMD3.cxx
17 @author Christian Holm Christensen <cholm@nbi.dk>
18 @date Sun Mar 26 18:26:12 2006
19 @brief Concrete implementation of AliFMDDetector for FMD3
20*/
21//____________________________________________________________________
22//
23// Concrete implementation of AliFMDDetector
24//
25// This implements the geometry for FMD3.
26// This has 2 rings.
27// The support of the FMD3 is a carbon-fibre cone, attached to the ITS
28// support via flanges. The cone also supports the beam-pipe.
29// The support is a special cone of carbon-fibre made by a Danish
30// Yacht company.
31//
32#include "AliFMD3.h" // ALIFMD3_H
33#include "AliLog.h" // ALILOG_H
34#include "AliFMDRing.h" // ALIFMDRING_H
35// #include <TMath.h> // ROOT_TMath
36
37//====================================================================
38ClassImp(AliFMD3)
39#if 0
40 ; // This is here to keep Emacs for indenting the next line
41#endif
42
43//____________________________________________________________________
44AliFMD3::AliFMD3(AliFMDRing* inner, AliFMDRing* outer)
45 : AliFMDDetector(3, inner, outer),
46 fNoseZ(0),
47 fNoseLowR(0),
48 fNoseHighR(0),
49 fNoseLength(0),
50 fBackLowR(0),
51 fBackHighR(0),
52 fBackLength(0),
53 fBeamThickness(0),
54 fBeamWidth(0),
55 fConeLength(0),
56 fFlangeR(0),
57 fZ(0),
58 fAlpha(0),
59 fNBeam(0),
60 fNFlange(0)
61{
62 // Constructor.
63 SetInnerZ(-62.8);
64 SetOuterZ(-75.2);
65 SetNoseZ();
66 SetNoseLowR();
67 SetNoseHighR();
68 SetNoseLength();
69 SetBackLowR();
70 SetBackHighR();
71 SetBackLength();
72 SetBeamThickness();
73 SetBeamWidth();
74 SetConeLength();
75 SetFlangeR();
76 SetNBeam();
77 SetNFlange();
78}
79
80//____________________________________________________________________
81void
82AliFMD3::Init()
83{
84 // Initialize
85 AliFMDDetector::Init();
86 SetInnerHoneyHighR(GetOuterHoneyHighR());
87 Double_t zdist = fConeLength - fBackLength - fNoseLength;
88 Double_t tdist = fBackHighR - fNoseHighR;
89 Double_t innerZh = fInnerZ - fInner->GetRingDepth() - fHoneycombThickness;
90 Double_t outerZh = fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness;
91 Double_t minZ = TMath::Min(fNoseZ - fConeLength, outerZh);
92 fAlpha = tdist / zdist;
93 fZ = fNoseZ + (minZ - fNoseZ) / 2;
94 fInnerHoneyHighR = ConeR(innerZh + fHoneycombThickness,"O") - 1;
95 fOuterHoneyHighR = GetBackLowR();
96}
97
98//____________________________________________________________________
99Double_t
100AliFMD3::ConeR(Double_t z, Option_t* opt) const
101{
102 // Calculate the cone radius at Z
103 if (fAlpha < 0) {
104 AliWarning(Form("alpha not set: %lf", fAlpha));
105 return -1;
106 }
107 if (z > fNoseZ) {
108 AliWarning(Form("z=%lf is before start of cone %lf", z, fNoseZ));
109 return -1;
110 }
111 if (z < fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness) {
112 AliWarning(Form("z=%lf is after end of cone %lf", z,
113 fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness));
114 return -1;
115 }
116 Double_t e = fBeamThickness / TMath::Cos(TMath::ATan(fAlpha));
117 if (opt[0] == 'I' || opt[1] == 'i') e *= -1;
118 if (z > fNoseZ - fNoseLength) return fNoseHighR + e;
119 if (z < fNoseZ - fConeLength + fBackLength) return fBackHighR + e;
120 Double_t r = fNoseHighR + fAlpha * TMath::Abs(z - fNoseZ + fNoseLength) + e;
121 return r;
122}
123
124
125//____________________________________________________________________
126//
127// EOF
128//