+++ /dev/null
-/**************************************************************************
- * Copyright(c) 2004, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-/* $Id$ */
-/** @file AliFMDRing.cxx
- @author Christian Holm Christensen <cholm@nbi.dk>
- @date Mon Mar 27 12:47:43 2006
- @brief FMD ring geometry parameters
-*/
-//__________________________________________________________________
-//
-// Utility class to help implement collection of FMD modules into
-// rings. This is used by AliFMDDetector and AliFMDGeometry.
-// The AliFMDGeometry object owns the AliFMDRing objects, and the
-// AliFMDDetector objects reference these. That is, the AliFMDRing
-// objects are share amoung the AliFMDDetector objects.
-//
-// Latest changes by Christian Holm Christensen
-//
-
-#include <TMath.h> // ROOT_TMath
-#include <TVector2.h> // ROOT_TVector2
-
-// #include <AliLog.h> // ALILOG_H
-#include "AliFMDRing.h" // ALIFMDRING_H
-
-//====================================================================
-ClassImp(AliFMDRing)
-#if 0
- ; // This is here to keep Emacs for indenting the next line
-#endif
-
-//____________________________________________________________________
-AliFMDRing::AliFMDRing(Char_t id)
- : TNamed(Form("FMD%c", id), "Forward multiplicity ring"),
- fId(id),
- fBondingWidth(0),
- fWaferRadius(0),
- fSiThickness(0),
- fLowR(0),
- fHighR(0),
- fMinR(0),
- fMaxR(0),
- fTheta(0),
- fNStrips(0),
- fRingDepth(0),
- fLegRadius(0),
- fLegLength(0),
- fLegOffset(0),
- fModuleSpacing(0),
- fPrintboardThickness(0),
- fCopperThickness(0),
- fChipThickness(0),
- fSpacing(0),
- fHoneycombThickness(0.),
- fAlThickness(0.),
- fVerticies(0),
- fSensorVerticies(0),
- fHybridVerticies(0),
- fFeetPositions(0)
-{
- // Constructor
- //
- // Parameters:
- // id Type of ring (either 'I' or 'O')
- //
- SetBondingWidth();
- SetWaferRadius();
- SetSiThickness();
- SetLegRadius();
- SetLegLength();
- SetLegOffset();
- SetModuleSpacing();
- SetPrintboardThickness();
- SetCopperThickness();
- SetChipThickness();
- SetSpacing();
- SetHoneycombThickness();
- SetAlThickness();
-
- if (fId == 'I' || fId == 'i') {
- SetLowR(4.3);
- SetHighR(17.2);
- SetTheta(36/2);
- SetNStrips(512);
- Double_t base = 0; // 4.1915;
- fFeetPositions.Add(new TVector2( 0.0551687, 8.0534-base));
- fFeetPositions.Add(new TVector2( 2.9993, 12.9457-base));
- fFeetPositions.Add(new TVector2(-2.9062, 12.9508-base));
-
- fHybridVerticies.Add(new TVector2(0.0000, 4.1700));
- fHybridVerticies.Add(new TVector2(1.0574, 4.1700));
- fHybridVerticies.Add(new TVector2(4.6614, 15.2622));
- fHybridVerticies.Add(new TVector2(0.9643, 17.4000));
- fHybridVerticies.Add(new TVector2(0.0000, 17.4000));
-
- fSensorVerticies.Add(new TVector2(0.0000, 4.1915));
- fSensorVerticies.Add(new TVector2(1.5793, 4.1915));
- fSensorVerticies.Add(new TVector2(5.2293, 15.4251));
- fSensorVerticies.Add(new TVector2(1.9807, 17.3035));
- fSensorVerticies.Add(new TVector2(0.0000, 17.3035));
-
- fVerticies.Add(new TVector2(0.0000, 4.3000));
- fVerticies.Add(new TVector2(1.3972, 4.3000));
- fVerticies.Add(new TVector2(4.9895, 15.3560));
- fVerticies.Add(new TVector2(1.8007, 17.2000));
- fVerticies.Add(new TVector2(0.0000, 17.2000));
- }
- else if (fId == 'O' || fId == 'o') {
- SetLowR(15.6);
- SetHighR(28.0);
- SetTheta(18/2);
- SetNStrips(256);
- Double_t base = 0; // 14.9104;
- fFeetPositions.Add(new TVector2(-1.72540000, 20.6267-base));
- fFeetPositions.Add(new TVector2( 1.72900000, 20.6267-base));
- fFeetPositions.Add(new TVector2( 0.00177616, 26.6007-base));
-
- fHybridVerticies.Add(new TVector2(0.0000, 14.9104));
- fHybridVerticies.Add(new TVector2(2.0783, 14.9104));
- fHybridVerticies.Add(new TVector2(3.9202, 26.5395));
- fHybridVerticies.Add(new TVector2(0.6784, 28.2500));
- fHybridVerticies.Add(new TVector2(0.0000, 28.2500));
-
- fSensorVerticies.Add(new TVector2(0.0000, 15.0104));
- fSensorVerticies.Add(new TVector2(2.5799, 15.0104));
- fSensorVerticies.Add(new TVector2(4.4439, 26.7766));
- fSensorVerticies.Add(new TVector2(1.8350, 28.1500));
- fSensorVerticies.Add(new TVector2(0.0000, 28.1500));
-
- fVerticies.Add(new TVector2(0.0000, 15.2104));
- fVerticies.Add(new TVector2(2.4091, 15.2104));
- fVerticies.Add(new TVector2(4.2231, 26.6638));
- fVerticies.Add(new TVector2(1.8357, 27.9500));
- fVerticies.Add(new TVector2(0.0000, 27.9500));
- }
-}
-
-//____________________________________________________________________
-void
-AliFMDRing::Init()
-{
- // Initialize
- //
- // All derived quantities are calculated here.
- //
-#if 0
- Double_t tanTheta = TMath::Tan(fTheta * TMath::Pi() / 180.);
- Double_t tanTheta2 = TMath::Power(tanTheta,2);
- Double_t r2 = TMath::Power(fWaferRadius,2);
- Double_t yA = tanTheta * fLowR;
- Double_t lr2 = TMath::Power(fLowR, 2);
- Double_t hr2 = TMath::Power(fHighR,2);
- Double_t xD = fLowR + TMath::Sqrt(r2 - tanTheta2 * lr2);
- Double_t xD2 = TMath::Power(xD,2);
- Double_t yB = TMath::Sqrt(r2 - hr2 + 2 * fHighR * xD - xD2);
- Double_t xC = ((xD + TMath::Sqrt(-tanTheta2 * xD2 + r2
- + r2 * tanTheta2))
- / (1 + tanTheta2));
- Double_t yC = tanTheta * xC;
-
- fVerticies.Expand(6);
- fVerticies.AddAt(new TVector2(fLowR, -yA), 0);
- fVerticies.AddAt(new TVector2(xC, -yC), 1);
- fVerticies.AddAt(new TVector2(fHighR, -yB), 2);
- fVerticies.AddAt(new TVector2(fHighR, yB), 3);
- fVerticies.AddAt(new TVector2(xC, yC), 4);
- fVerticies.AddAt(new TVector2(fLowR, yA), 5);
-#endif
-
- // A's length. Corresponds to distance from nominal beam line to the
- // cornor of the active silicon element.
- fMinR = GetVertex(1)->Mod(); // GetVertex(5)->Mod();
- // A's length. Corresponds to distance from nominal beam line to the
- // cornor of the active silicon element.
- fMaxR = fHighR;
-
- fRingDepth = (fSiThickness + fPrintboardThickness
- + fCopperThickness + fChipThickness
- + fLegLength + fModuleSpacing + fSpacing);
-}
-
-//____________________________________________________________________
-TVector2*
-AliFMDRing::GetVertex(Int_t i) const
-{
- // Get the i'th vertex of polygon shape
- //
- // the polygon shape describes the shape of the rings' sensors
- //
- // Parameters:
- // i The vertex number to get (from 0 to 5)
- return static_cast<TVector2*>(fVerticies.At(i));
-}
-
-//____________________________________________________________________
-TVector2*
-AliFMDRing::GetSensorVertex(Int_t i) const
-{
- // Get the i'th vertex of polygon shape
- //
- // the polygon shape describes the shape of the rings' sensors
- //
- // Parameters:
- // i The vertex number to get (from 0 to 5)
- return static_cast<TVector2*>(fSensorVerticies.At(i));
-}
-
-//____________________________________________________________________
-TVector2*
-AliFMDRing::GetHybridVertex(Int_t i) const
-{
- // Get the i'th vertex of polygon shape
- //
- // the polygon shape describes the shape of the rings' hybrid cards
- //
- // Parameters:
- // i The vertex number to get (from 0 to 5)
- return static_cast<TVector2*>(fHybridVerticies.At(i));
-}
-
-//____________________________________________________________________
-TVector2*
-AliFMDRing::GetFootPosition(Int_t i) const
-{
- // Get the i'th vertex of polygon shape
- //
- // The feet are attached to the hybrid cards
- //
- // Parameters:
- // i The foot number to get (from 0 to 2)
- return static_cast<TVector2*>(fFeetPositions.At(i));
-}
-
-//____________________________________________________________________
-Double_t
-AliFMDRing::GetStripRadius(UShort_t strip) const
-{
- // Return the nominal strip radius
- //
- // Parameter
- // strip Strip number (0-511 for inners, 0-255 for outers)
- Double_t rmax = GetMaxR();
- Double_t stripoff = GetMinR();
- Double_t dstrip = (rmax - stripoff) / GetNStrips();
- return (strip + .5) * dstrip + stripoff; // fLowR
-}
-
-//____________________________________________________________________
-Double_t
-AliFMDRing::GetModuleDepth() const
-{
- // Get the total depth of a module (sensor + hybrid card)
- //
- // The depth is the sum of
- //
- // The silicon thickness
- // The thickness of spacers between the silicon and hybrid
- // The thickness of the hybrid PCB
- // The thickness of the copper layer in the PCB
- // The thickness of the chip layer in the PCB
- // The height of the legs
- return (GetSiThickness()
- + GetSpacing()
- + GetPrintboardThickness()
- + GetCopperThickness()
- + GetChipThickness()
- + GetLegLength());
-
-}
-
-//____________________________________________________________________
-Double_t
-AliFMDRing::GetFullDepth() const
-{
- // Get the full depth of this ring, including the honeycomb,
- // digitizer and card.
- return (GetModuleDepth()
- + GetModuleSpacing()
- + GetHoneycombThickness()
- + GetFMDDPrintboardThickness()
- + GetFMDDCopperThickness()
- + GetFMDDChipThickness()
- + 0.5);
-}
-
-//____________________________________________________________________
-void
-AliFMDRing::Detector2XYZ(UShort_t sector,
- UShort_t strip,
- Double_t& x,
- Double_t& y,
- Double_t& z) const
-{
- // Translate detector coordinates (this,sector,strip) to global
- // coordinates (x,y,z)
- //
- // Parameters
- // sector Sector number in this ring
- // strip Strip number in this ring
- // x On return, the global X coordinate
- // y On return, the global Y coordinate
- // z On return, the z coordinate in the ring plane
- //
- // The ring plane is the plane half way between the two sensor
- // layers.
- if (sector >= GetNSectors()) {
- Error("Detector2XYZ", "Invalid sector number %d (>=%d) in ring %c",
- sector, GetNSectors(), fId);
- return;
- }
- if (strip >= GetNStrips()) {
- Error("Detector2XYZ", "Invalid strip number %d (>=%d) for ring type '%c'",
- strip, GetNStrips(), fId);
- return;
- }
- Double_t phi = Float_t(sector + .5) / GetNSectors() * 2 * TMath::Pi();
- Double_t r = Float_t(strip + .5) / GetNStrips() * (fHighR - fLowR) + fLowR;
- x = r * TMath::Cos(phi);
- y = r * TMath::Sin(phi);
- if (((sector / 2) % 2) == 1)
- z += TMath::Sign(fModuleSpacing, z);
-}
-
-//____________________________________________________________________
-Bool_t
-AliFMDRing::XYZ2Detector(Double_t x,
- Double_t y,
- Double_t z,
- UShort_t& sector,
- UShort_t& strip) const
-{
- // Translate global coordinates (x,y,z) to detector coordinates
- // (this,sector,strip)
- //
- // Parameters:
- // x Global x coordinate
- // y Global y coordinate
- // z Global y coordinate
- // sector On return, the sector number in this ring
- // strip On return, the strip number in this ring
- //
- sector = strip = 0;
- Double_t r = TMath::Sqrt(x * x + y * y);
- Int_t str = Int_t((r - fMinR) / GetPitch());
- if (str < 0 || str >= GetNStrips()) return kFALSE;
-
- Double_t phi = TMath::ATan2(y, x) * 180. / TMath::Pi();
- if (phi < 0) phi = 360. + phi;
- Int_t sec = Int_t(phi / fTheta);
- if (sec < 0 || sec >= GetNSectors()) return kFALSE;
- if ((sec / 2) % 2 == 1) {
- if (TMath::Abs(z - TMath::Sign(fModuleSpacing, z)) >= 0.01)
- return kFALSE;
- }
- else if (TMath::Abs(z) >= 0.01) return kFALSE;
-
- strip = str;
- sector = sec;
- return kTRUE;
-}
-//____________________________________________________________________
-Float_t
-AliFMDRing::GetStripLength(UShort_t strip) const
-{
- // Get the length of a strip
- //
- // Parameters:
- // strip Strip number (0-511 for inners, 0-255 for outers)
- //
- if(strip >= GetNStrips())
- Error("GetStripLength", "Invalid strip number %d (>=%d) for ring type %c",
- strip, GetNStrips(), fId);
-
- Float_t rad = GetMaxR()-GetMinR();
-
- Float_t segment = rad / GetNStrips();
-
- TVector2* corner1 = GetVertex(2);
- TVector2* corner2 = GetVertex(3);
-
- Float_t slope = ((corner1->Y() - corner2->Y()) /
- (corner1->X() - corner2->X()));
- Float_t constant = ((corner2->Y() * corner1->X() -
- (corner2->X()*corner1->Y())) /
- (corner1->X() - corner2->X()));
- Float_t radius = GetMinR() + strip*segment;
-
- Float_t d = (TMath::Power(TMath::Abs(radius*slope),2) +
- TMath::Power(radius,2) - TMath::Power(constant,2));
-
- Float_t arclength = GetBaseStripLength(strip);
- if(d>0) {
- Float_t x = ((-1 * TMath::Sqrt(d) -slope*constant) /
- (1 + TMath::Power(slope,2)));
- Float_t y = slope*x + constant;
- Float_t theta = TMath::ATan2(x,y);
-
- if(x < corner1->X() && y > corner1->Y()) {
- //One sector since theta is by definition half-hybrid
- arclength = radius*theta;
- }
- }
-
- return arclength;
-
-
-}
-//____________________________________________________________________
-Float_t
-AliFMDRing::GetBaseStripLength(UShort_t strip) const
-{
- // Get the basic strip length
- //
- // Parameters:
- // strip Strip number
- Float_t rad = GetMaxR()-GetMinR();
- Float_t segment = rad / GetNStrips();
- Float_t basearc = 2*TMath::Pi() / (0.5*GetNSectors());
- Float_t radius = GetMinR() + strip*segment;
- Float_t basearclength = 0.5*basearc * radius;
-
- return basearclength;
-}
-//
-// EOF
-//