* 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
// Latest changes by Christian Holm Christensen
//
-#include <AliLog.h> // ALILOG_H
+// #include <AliLog.h> // ALILOG_H
#include "AliFMDRing.h" // ALIFMDRING_H
-#include <TMath.h> // ROOT_TMath
+// #include <TMath.h> // ROOT_TMath
#include <TVector2.h> // ROOT_TVector2
//====================================================================
fId(id),
fVerticies(0)
{
+ // CTOR
SetBondingWidth();
SetWaferRadius();
SetSiThickness();
SetPrintboardThickness();
SetCopperThickness();
SetChipThickness();
+ SetSpacing();
if (fId == 'I' || fId == 'i') {
SetLowR(4.3);
void
AliFMDRing::Init()
{
+ // Initialize
Double_t tanTheta = TMath::Tan(fTheta * TMath::Pi() / 180.);
Double_t tanTheta2 = TMath::Power(tanTheta,2);
Double_t r2 = TMath::Power(fWaferRadius,2);
TVector2*
AliFMDRing::GetVertex(Int_t i) const
{
+ // Get the i'th vertex of polygon shape
return static_cast<TVector2*>(fVerticies.At(i));
}
+//____________________________________________________________________
+Double_t
+AliFMDRing::GetStripRadius(UShort_t strip) const
+{
+ // Return the nominal strip radius
+ Double_t rmax = GetMaxR();
+ Double_t stripoff = GetMinR();
+ Double_t dstrip = (rmax - stripoff) / GetNStrips();
+ return (strip + .5) * dstrip + stripoff; // fLowR
+}
+
//____________________________________________________________________
void
AliFMDRing::Detector2XYZ(UShort_t sector,
Double_t& y,
Double_t& z) const
{
+ // Translate detector coordinates (this,sector,strip) to global
+ // coordinates (x,y,z)
if (sector >= GetNSectors()) {
Error("Detector2XYZ", "Invalid sector number %d (>=%d) in ring %c",
sector, GetNSectors(), fId);
UShort_t& sector,
UShort_t& strip) const
{
+ // Translate global coordinates (x,y,z) to detector coordinates
+ // (this,sector,strip)
sector = strip = 0;
Double_t r = TMath::Sqrt(x * x + y * y);
Int_t str = Int_t((r - fMinR) / GetPitch());