FMD/AliFMD3.cxx

   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
  16 /* $Id$ */
  17
  18 //____________________________________________________________________
  19 //
  20 // Concrete implementation of AliFMDSubDetector
  21 //
  22 // This implements the geometry for FMD3
  23 //
  24 #include "TVirtualMC.h"         // ROOT_TVirtualMC
  25 #include "TCONS.h"              // ROOT_TCONS
  26 #include "TNode.h"              // ROOT_TNode
  27 #include "TList.h"              // ROOT_TList
  28 #include "AliFMD3.h"            // ALIFMD3_H
  29 #include "AliLog.h"             // ALILOG_H
  30 #include "AliFMDRing.h"         // ALIFMDRING_H
  31 #include <Riostream.h>          // ROOT_Riostream
  32
  33 //____________________________________________________________________
  34 ClassImp(AliFMD3)
  35
  36 //____________________________________________________________________
  37 AliFMD3::AliFMD3()
  38   : AliFMDSubDetector(3),
  39     fVolumeId(0)
  40 {
  41   // Default constructor for the FMD3 sub-detector
  42   AliDebug(10, "\t\tDefault CTOR");
  43 }
  44
  45
  46 //____________________________________________________________________
  47 AliFMD3::~AliFMD3()
  48 {
  49   // Destructor - does nothing
  50   AliDebug(10, "\t\tDTOR");
  51 }
  52
  53
  54 //____________________________________________________________________
  55 void
  56 AliFMD3::SetupGeometry(Int_t airId, Int_t alId, Int_t carbonId)
  57 {
  58   // Setup the FMD3 sub-detector geometry
  59   //
  60   // Parameters:
  61   //
  62   //     airId         Id # of the Air medium
  63   //     kaptionId     Id # of the Aluminium medium
  64   //
  65   AliDebug(10, "\tSetting up the geometry for FMD3");
  66   Double_t innerZl = fInnerZ;
  67   Double_t innerZh = (fInnerZ
  68                       - fInner->GetModuleSpacing()
  69                       - fInner->GetLegLength()
  70                       - fInner->GetSiThickness()
  71                       - fInner->GetPrintboardThickness()
  72                       - fHoneycombThickness);
  73   Double_t innerRl = fInner->GetLowR();
  74   Double_t outerZl = fOuterZ;
  75   Double_t outerZh = (fOuterZ
  76                       - fOuter->GetModuleSpacing()
  77                       - fOuter->GetLegLength()
  78                       - fOuter->GetSiThickness()
  79                       - fOuter->GetPrintboardThickness()
  80                       - fHoneycombThickness);
  81   Double_t outerRl = fOuter->GetLowR();
  82
  83   fSupport.SetupGeometry(airId, carbonId,
  84                          innerZl, innerZh, innerRl,
  85                          outerZl, outerZh, outerRl);
  86
  87   fInnerHoneyLowR  = fInner->GetLowR() + 1;
  88   fInnerHoneyHighR = fSupport.ConeR(innerZh + fHoneycombThickness, "I");
  89   fOuterHoneyLowR  = fOuter->GetLowR() + 1;
  90   fOuterHoneyHighR = fSupport.GetBackLowR();
  91
  92   // Identity matrix
  93   gMC->Matrix(fRotationId, 90, 0, 90, 90, 0, 0);
  94   //0, 180, 90, 90, 180, 0);
  95
  96
  97   AliFMDSubDetector::SetupGeometry(airId, alId, carbonId);
  98 }
  99
 100 //____________________________________________________________________
 101 void
 102 AliFMD3::Geometry(const char* mother, Int_t pbRotId,
 103                   Int_t idRotId, Double_t z)
 104 {
 105   // Position the FMD3 sub-detector volume
 106   //
 107   // Parameters
 108   //
 109   //     mother     name of the mother volume
 110   //     pbRotId    Printboard roation matrix ID
 111   //     idRotId    Identity rotation matrix ID
 112   //     z          Z position (not really used here, but passed down)
 113   //
 114   z = fSupport.GetZ();
 115   fSupport.Geometry(mother, fRotationId, z);
 116   AliDebug(10, Form("\t\tPassing z=%lf to ring volumes", z));
 117   AliFMDSubDetector::Geometry("FMD3", pbRotId, idRotId, z);
 118 }
 119
 120
 121 //____________________________________________________________________
 122 void
 123 AliFMD3::SimpleGeometry(TList* nodes,
 124                         TNode* mother,
 125                         Int_t colour,
 126                         Double_t zMother)
 127 {
 128   // We need to get the equation for the line that connects the
 129   // outer circumfrences of the two rings, as  well as for the line
 130   // that connects the inner curcumfrences, so that we can project to
 131   // where the honey-comb actually ends.
 132   //
 133   // we have
 134   //
 135   //   y = a * x + b
 136   //   b = y - a * x;
 137   //
 138   // For the outer line, we have the two equations
 139   //
 140   //    fOuterHoneyHighR = a * x1 + b;
 141   //    fInnerHoneyHighR = a * x2 + b;
 142   //
 143   // where
 144   //
 145   //    x1 = (fOuterZ + fOuter->fSiThickness + fOuter->fPrintboardThickness
 146   //          + fOuter->fLegLength + fModuleSpacing)
 147   //       = fInner - fDz + fHoneycombThickness
 148   //    x2 = (fInnerZ + fInner->fSiThickness + fInner->fPrintboardThickness
 149   //          + fInner->fLegLength + fModuleSpacing)
 150   //
 151   // and
 152   //
 153   //    a  = (fOuterHoneyHighR - fInnerHoneyHighR) / (x1 - x2)
 154   //
 155   //
 156   AliDebug(10, "\tCreating simplified geometry for FMD3");
 157   Double_t dz = (TMath::Abs(fInnerZ - fOuterZ)
 158                  + fOuter->GetSiThickness()
 159                  + fOuter->GetPrintboardThickness()
 160                  + fOuter->GetLegLength()
 161                  + fOuter->GetModuleSpacing()
 162                  + fHoneycombThickness) / 2;
 163 #if 1
 164   Double_t x1  = (fOuterZ - (fOuter->GetSiThickness()
 165                              + fOuter->GetPrintboardThickness()
 166                              + fOuter->GetLegLength()
 167                              + fOuter->GetModuleSpacing()));
 168   Double_t x2  = (fInnerZ - (fInner->GetSiThickness()
 169                              + fInner->GetPrintboardThickness()
 170                              + fInner->GetLegLength()
 171                              + fInner->GetModuleSpacing()));
 172   Double_t ao   = 0;
 173   Double_t ao1  = (fOuterHoneyHighR - fInnerHoneyHighR) / (x1 - x2);
 174   Double_t ao2  = ((fOuter->GetHighR() - fInner->GetHighR())
 175                    / (fOuterZ - fInnerZ));
 176   Double_t bo   = 0;
 177   if (ao2 > ao1) {
 178     // std::cout << "Wafer determinds the size" << std::endl;
 179     ao  = ao2;
 180     bo  = fInner->GetHighR() - ao * fInnerZ;
 181   }
 182   else {
 183     ao = ao1;
 184     bo = fOuterHoneyHighR - ao * x1;
 185   }
 186
 187   Double_t y1o = ao * (fInnerZ - 2 * dz) + bo;
 188   Double_t y2o = ao * fInnerZ + bo;
 189 #endif
 190   // We probably need to make a PCON here.
 191   TShape* shape = new TCONS("FMD3", "FMD3", "", dz,
 192                             fOuter->GetLowR(),  y1o, /* fOuterHoneyHighR, */
 193                             fInner->GetLowR(),  y2o, /* fInnerHoneyHighR, */
 194                             0, 360);
 195   mother->cd();
 196   zMother = fInnerZ - dz;
 197   TNode* node = new TNode("FMD3", "FMD3", shape, 0, 0, zMother, 0);
 198   node->SetVisibility(0);
 199   nodes->Add(node);
 200   AliFMDSubDetector::SimpleGeometry(nodes, node, colour, zMother);
 201 }
 202
 203 //____________________________________________________________________
 204 void
 205 AliFMD3::Gsatt() const
 206 {
 207   // Set draw attributes for the FMD3
 208   AliDebug(10, "Setting drawing attributes for FMD3");
 209   AliFMDSubDetector::Gsatt();
 210   fSupport.Gsatt();
 211 }
 212
 213 //____________________________________________________________________
 214 //
 215 // EOF
 216 //