]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - FMD/AliFMD3.cxx
Update by Raphaelle
[u/mrichter/AliRoot.git] / FMD / AliFMD3.cxx
index 3fc331a74fb0deffd63a669086263f8d5d8e73de..e6a8d44d61299648f3898a39f89e61845f39d6d8 100644 (file)
  * about the suitability of this software for any purpose. It is          *
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
-
 /* $Id$ */
-
+/** @file    AliFMD3.cxx
+    @author  Christian Holm Christensen <cholm@nbi.dk>
+    @date    Sun Mar 26 18:26:12 2006
+    @brief   Concrete implementation of AliFMDDetector for FMD3
+*/
 //____________________________________________________________________
 //                                                                          
-// Concrete implementation of AliFMDSubDetector 
+// Concrete implementation of AliFMDDetector 
 //
-// This implements the geometry for FMD3
+// This implements the geometry for FMD3.
+// This has 2 rings.
+// The support of the FMD3 is a carbon-fibre cone, attached to the ITS
+// support via flanges.  The cone also supports the beam-pipe.
+// The support is a special cone of carbon-fibre made by a Danish
+// Yacht company.
 //
-#ifndef ROOT_TVirtualMC
-# include <TVirtualMC.h>
-#endif
-#ifndef ROOT_TCONS
-# include <TCONS.h>
-#endif
-#ifndef ROOT_TNode
-# include <TNode.h>
-#endif
-#ifndef ROOT_TList
-# include <TList.h>
-#endif
-#ifndef ALIFMD3_H
-# include "AliFMD3.h"
-#endif 
-#ifndef ALILOG_H
-# include "AliLog.h"
-#endif
-#ifndef ALIFMDRING_H
-# include "AliFMDRing.h"
-#endif 
-#include <iostream>
 
-//____________________________________________________________________
-ClassImp(AliFMD3);
-
-//____________________________________________________________________
-AliFMD3::AliFMD3() 
-  : AliFMDSubDetector(3), 
-    fVolumeId(0),
-    fDz(0)
-{
-  // Default constructor for the FMD3 sub-detector 
-}
+#include <TMath.h>             // ROOT_TMath
 
+#include "AliFMD3.h"           // ALIFMD3_H 
+#include "AliLog.h"            // ALILOG_H
+#include "AliFMDRing.h"                // ALIFMDRING_H 
 
-//____________________________________________________________________
-AliFMD3::~AliFMD3() 
-{
-  // Destructor - does nothing 
-}
-
+//====================================================================
+ClassImp(AliFMD3)
+#if 0
+  ; // This is here to keep Emacs for indenting the next line
+#endif
 
 //____________________________________________________________________
-void 
-AliFMD3::SetupGeometry(Int_t airId, Int_t kaptionId) 
+AliFMD3::AliFMD3(AliFMDRing* inner, AliFMDRing* outer) 
+  : AliFMDDetector(3, inner, outer),
+    fNoseZ(0),
+    fNoseLowR(0),
+    fNoseHighR(0),
+    fNoseLength(0),
+    fBackLowR(0),
+    fBackHighR(0),
+    fBackLength(0),
+    fBeamThickness(0),
+    fBeamWidth(0),
+    fConeLength(0),
+    fFlangeR(0),
+    fZ(0),
+    fAlpha(0), 
+    fNBeam(0), 
+    fNFlange(0)
 {
-  // Setup the FMD3 sub-detector geometry 
-  // 
-  // Parameters:
-  // 
-  //     airId         Id # of the Air medium 
-  //     kaptionId     Id # of the Aluminium medium 
-  // 
-  fInnerHoneyLowR  = fInner->GetLowR() + 1;
-  fInnerHoneyHighR = fInner->GetHighR() + 1;
-  fOuterHoneyLowR  = fOuter->GetLowR() + 1;
-  fOuterHoneyHighR = fOuter->GetHighR() + 1;
-
-  CalculateDz();
-  Double_t par[3];
-  par[0] = fInner->GetLowR();
-  par[1] = fOuterHoneyHighR;
-  par[2] = fDz;
-  fVolumeId = gMC->Gsvolu("FMD3", "TUBE", airId, par, 3);
-
-  gMC->Matrix(fRotationId, 90, 0, 90, 90, 0, 0); 
-  //0, 180, 90, 90, 180, 0);
-
-  AliFMDSubDetector::SetupGeometry(airId, kaptionId);
+  // Constructor. 
+  SetInnerZ(-62.8);
+  SetOuterZ(-75.2);
+  SetNoseZ();
+  SetNoseLowR();
+  SetNoseHighR();
+  SetNoseLength();
+  SetBackLowR();
+  SetBackHighR();
+  SetBackLength();
+  SetBeamThickness();
+  SetBeamWidth();
+  SetConeLength();
+  SetFlangeR();
+  SetNBeam();
+  SetNFlange();
 }
 
 //____________________________________________________________________
-void 
-AliFMD3::Geometry(const char* mother, Int_t pbRotId, 
-                 Int_t idRotId, Double_t z) 
+void
+AliFMD3::Init() 
 {
-  // Position the FMD3 sub-detector volume 
-  // 
-  // Parameters 
-  //
-  //     mother     name of the mother volume 
-  //     pbRotId    Printboard roation matrix ID 
-  //     idRotId    Identity rotation matrix ID 
-  //     z          Z position (not really used here, but passed down)
-  //
-  z = fInnerZ - fDz;
-  gMC->Gspos("FMD3", 1, mother, 0, 0, z, fRotationId);
-  
-  AliFMDSubDetector::Geometry("FMD3", pbRotId, idRotId, z);
+  // Initialize 
+  AliFMDDetector::Init();
+  SetInnerHoneyHighR(GetOuterHoneyHighR());
+  Double_t zdist   = fConeLength - fBackLength - fNoseLength;
+  Double_t tdist   = fBackHighR - fNoseHighR;
+  Double_t innerZh = fInnerZ - fInner->GetRingDepth() - fHoneycombThickness;
+  Double_t outerZh = fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness;
+  Double_t minZ    = TMath::Min(fNoseZ - fConeLength, outerZh);
+  fAlpha           = tdist / zdist;
+  fZ               = fNoseZ + (minZ - fNoseZ) / 2;
+  fInnerHoneyHighR = ConeR(innerZh + fHoneycombThickness,"O") - 1;
+  fOuterHoneyHighR = GetBackLowR();
 }
 
-  
 //____________________________________________________________________
-void 
-AliFMD3::SimpleGeometry(TList* nodes, 
-                       TNode* mother, 
-                       Int_t colour, 
-                       Double_t zMother) 
+Double_t
+AliFMD3::ConeR(Double_t z, Option_t* opt) const
 {
-  // We need to get the equation for the line that connects the 
-  // outer circumfrences of the two rings, as  well as for the line
-  // that connects the inner curcumfrences, so that we can project to
-  // where the honey-comb actually ends. 
-  // 
-  // we have 
-  //   
-  //   y = a * x + b 
-  //   b = y - a * x;
-  // 
-  // For the outer line, we have the two equations 
-  // 
-  //    fOuterHoneyHighR = a * x1 + b;
-  //    fInnerHoneyHighR = a * x2 + b; 
-  // 
-  // where 
-  // 
-  //    x1 = (fOuterZ + fOuter->fSiThickness + fOuter->fPrintboardThickness 
-  //          + fOuter->fLegLength + fModuleSpacing) 
-  //       = fInner - fDz + fHoneycombThickness
-  //    x2 = (fInnerZ + fInner->fSiThickness + fInner->fPrintboardThickness 
-  //          + fInner->fLegLength + fModuleSpacing)
-  // 
-  // and 
-  //
-  //    a  = (fOuterHoneyHighR - fInnerHoneyHighR) / (x1 - x2)
-  //    
-  // 
-  CalculateDz();
-#if 1
-  Double_t x1  = (fOuterZ - (fOuter->GetSiThickness() 
-                            + fOuter->GetPrintboardThickness() 
-                            + fOuter->GetLegLength() 
-                            + fOuter->GetModuleSpacing()));
-  Double_t x2  = (fInnerZ - (fInner->GetSiThickness() 
-                            + fInner->GetPrintboardThickness() 
-                            + fInner->GetLegLength() 
-                            + fInner->GetModuleSpacing()));
-  Double_t ao   = 0;
-  Double_t ao1  = (fOuterHoneyHighR - fInnerHoneyHighR) / (x1 - x2);
-  Double_t ao2  = ((fOuter->GetHighR() - fInner->GetHighR()) 
-                  / (fOuterZ - fInnerZ));
-  Double_t bo   = 0;
-  if (ao2 > ao1) {
-    std::cout << "Wafer determinds the size" << std::endl;
-    ao  = ao2;
-    bo  = fInner->GetHighR() - ao * fInnerZ;
+  // Calculate the cone radius at Z
+  if (fAlpha < 0) {
+    AliWarning(Form("alpha not set: %lf", fAlpha));
+    return -1;
   }
-  else {
-    ao = ao1;
-    bo = fOuterHoneyHighR - ao * x1;
+  if (z > fNoseZ) {
+    AliWarning(Form("z=%lf is before start of cone %lf", z, fNoseZ));
+    return -1;
   }
-  
-  Double_t y1o = ao * (fInnerZ - 2 * fDz) + bo;
-  Double_t y2o = ao * fInnerZ + bo;
-#endif
-  // We probably need to make a PCON here. 
-  TShape* shape = new TCONS("FMD3", "FMD3", "", fDz, 
-                           fOuter->GetLowR(),  y1o, /* fOuterHoneyHighR, */
-                           fInner->GetLowR(),  y2o, /* fInnerHoneyHighR, */
-                           0, 360);
-  mother->cd();
-  zMother = fInnerZ - fDz;  
-  TNode* node = new TNode("FMD3", "FMD3", shape, 0, 0, zMother, 0);
-  node->SetVisibility(0);
-  nodes->Add(node);
-  AliFMDSubDetector::SimpleGeometry(nodes, node, colour, zMother);
+  if (z < fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness) {
+    AliWarning(Form("z=%lf is after end of cone %lf", z, 
+                   fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness));
+    return -1;
+  }
+  Double_t e = fBeamThickness / TMath::Cos(TMath::ATan(fAlpha));
+  if (opt[0] == 'I' || opt[1] == 'i') e *= -1;
+  if (z > fNoseZ - fNoseLength) return fNoseHighR + e;
+  if (z < fNoseZ - fConeLength + fBackLength) return fBackHighR + e;
+  Double_t r = fNoseHighR + fAlpha * TMath::Abs(z - fNoseZ + fNoseLength) + e;
+  return r;
 }
 
-//____________________________________________________________________
-void 
-AliFMD3::CalculateDz() 
-{
-  if (fDz > 0) return;
-  fDz = (TMath::Abs(fInnerZ - fOuterZ) 
-        + fOuter->GetSiThickness() 
-        + fOuter->GetPrintboardThickness() 
-        + fOuter->GetLegLength() 
-        + fOuter->GetModuleSpacing() 
-        + fHoneycombThickness) / 2;
-}
 
 //____________________________________________________________________
 //