[u/mrichter/AliRoot.git] / FMD / AliFMD3.cxx

/**************************************************************************
 * 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    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 AliFMDDetector 
//
// 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.
//
#include "AliFMD3.h"		// ALIFMD3_H 
#include "AliLog.h"		// ALILOG_H
#include "AliFMDRing.h"		// ALIFMDRING_H 
// #include <TMath.h>		// ROOT_TMath

//====================================================================
ClassImp(AliFMD3)
#if 0
  ; // This is here to keep Emacs for indenting the next line
#endif

//____________________________________________________________________
AliFMD3::AliFMD3(AliFMDRing* inner, AliFMDRing* outer) 
  : AliFMDDetector(3, inner, outer)
{
  // Constructor. 
  SetInnerZ(-62.8);
  SetOuterZ(-75.2);
  SetNoseZ();
  SetNoseLowR();
  SetNoseHighR();
  SetNoseLength();
  SetBackLowR();
  SetBackHighR();
  SetBackLength();
  SetBeamThickness();
  SetBeamWidth();
  SetConeLength();
  SetFlangeR();
  SetNBeam();
  SetNFlange();
}

//____________________________________________________________________
void
AliFMD3::Init() 
{
  // 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();
}

//____________________________________________________________________
Double_t
AliFMD3::ConeR(Double_t z, Option_t* opt) const
{
  // Calculate the cone radius at Z
  if (fAlpha < 0) {
    AliWarning(Form("alpha not set: %lf", fAlpha));
    return -1;
  }
  if (z > fNoseZ) {
    AliWarning(Form("z=%lf is before start of cone %lf", z, fNoseZ));
    return -1;
  }
  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;
}


//____________________________________________________________________
//
// EOF
//
Commit	Line	Data
	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	//====================================================================
	38	ClassImp(AliFMD3)
	39	#if 0
	40	; // This is here to keep Emacs for indenting the next line
	41	#endif
	42
	43	//____________________________________________________________________
	44	AliFMD3::AliFMD3(AliFMDRing* inner, AliFMDRing* outer)
	45	: AliFMDDetector(3, inner, outer)
	46	{
	47	// Constructor.
	48	SetInnerZ(-62.8);
	49	SetOuterZ(-75.2);
	50	SetNoseZ();
	51	SetNoseLowR();
	52	SetNoseHighR();
	53	SetNoseLength();
	54	SetBackLowR();
	55	SetBackHighR();
	56	SetBackLength();
	57	SetBeamThickness();
	58	SetBeamWidth();
	59	SetConeLength();
	60	SetFlangeR();
	61	SetNBeam();
	62	SetNFlange();
	63	}
	64
	65	//____________________________________________________________________
	66	void
	67	AliFMD3::Init()
	68	{
	69	// Initialize
	70	AliFMDDetector::Init();
	71	SetInnerHoneyHighR(GetOuterHoneyHighR());
	72	Double_t zdist = fConeLength - fBackLength - fNoseLength;
	73	Double_t tdist = fBackHighR - fNoseHighR;
	74	Double_t innerZh = fInnerZ - fInner->GetRingDepth() - fHoneycombThickness;
	75	Double_t outerZh = fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness;
	76	Double_t minZ = TMath::Min(fNoseZ - fConeLength, outerZh);
	77	fAlpha = tdist / zdist;
	78	fZ = fNoseZ + (minZ - fNoseZ) / 2;
	79	fInnerHoneyHighR = ConeR(innerZh + fHoneycombThickness,"O") - 1;
	80	fOuterHoneyHighR = GetBackLowR();
	81	}
	82
	83	//____________________________________________________________________
	84	Double_t
	85	AliFMD3::ConeR(Double_t z, Option_t* opt) const
	86	{
	87	// Calculate the cone radius at Z
	88	if (fAlpha < 0) {
	89	AliWarning(Form("alpha not set: %lf", fAlpha));
	90	return -1;
	91	}
	92	if (z > fNoseZ) {
	93	AliWarning(Form("z=%lf is before start of cone %lf", z, fNoseZ));
	94	return -1;
	95	}
	96	if (z < fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness) {
	97	AliWarning(Form("z=%lf is after end of cone %lf", z,
	98	fOuterZ - fOuter->GetRingDepth() - fHoneycombThickness));
	99	return -1;
	100	}
	101	Double_t e = fBeamThickness / TMath::Cos(TMath::ATan(fAlpha));
	102	if (opt[0] == 'I' \|\| opt[1] == 'i') e *= -1;
	103	if (z > fNoseZ - fNoseLength) return fNoseHighR + e;
	104	if (z < fNoseZ - fConeLength + fBackLength) return fBackHighR + e;
	105	Double_t r = fNoseHighR + fAlpha * TMath::Abs(z - fNoseZ + fNoseLength) + e;
	106	return r;
	107	}
	108
	109
	110	//____________________________________________________________________
	111	//
	112	// EOF
	113	//