* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.6 2001/01/30 12:19:39 morsch
-Update chamber positions (AdTDR version update 4/12/2000).
+/* $Id$ */
-Revision 1.5 2000/10/18 13:26:10 morsch
-New z-positions of chambers after Erice
+#include "AliMUONConstants.h"
-Revision 1.4 2000/10/06 09:09:56 morsch
-Outer radius of chambers adjusted to accomodate slat chambers (to be checked and updated).
+#include "AliLog.h"
-Revision 1.3 2000/10/02 16:58:29 egangler
-Cleaning of the code :
--> coding conventions
--> void Streamers
--> some useless includes removed or replaced by "class" statement
+#include "TMath.h"
+#include "TClass.h"
+#include "AliMpConstants.h"
-Revision 1.2 2000/06/27 09:46:57 morsch
-kMAXZOOM global constant now in AliMUONConstants
+///
+/// \class AliMUONConstants
+/// This class holds various constants to be used in many places,
+/// such as the number of tracking and trigger chambers,
+/// some geometrical constants (to build the initial geometry for simulation)
+/// and mathieson distribution default values.
+/// Those constants should as much as possible replace hard-coded values
+/// which are to be considered strictly illegal in the MUON code (or any code,
+/// by the way).
+///
-Revision 1.1 2000/06/26 14:02:38 morsch
-Add class AliMUONConstants with MUON specific constants using static memeber data and access methods.
+/// \cond CLASSIMP
+ClassImp(AliMUONConstants)
+/// \endcond
-*/
+Int_t AliMUONConstants::fgNTrackingSt = 5;
+Int_t AliMUONConstants::fgNTriggerCircuit = 234;
+Int_t AliMUONConstants::fgNDetElem = 228;
+Int_t AliMUONConstants::fgNGeomModules = 20;
+Float_t AliMUONConstants::fgkTriggerTofLimit = 75E-9;
-#include "AliMUONConstants.h"
+Float_t AliMUONConstants::fgDefaultChamberZ[14] =
+ {-526.16, -545.24, -676.4, -695.4, // St12
+ -967.5, -998.5, -1276.5, -1307.5, -1406.6, -1437.6,// updated 08/05, EDMS id 335328 (A. Tournaire)
+ -1603.5, -1620.5, -1703.5, -1720.5}; // M1 & M2
+Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] =
+{1., 1.01060, 1.06236, 1.07296};
-ClassImp(AliMUONConstants)
-Int_t AliMUONConstants::fgNCh = 14;
-Int_t AliMUONConstants::fgNTrackingCh = 10;
-Int_t AliMUONConstants::fgNTriggerCh = 4;
-Int_t AliMUONConstants::fgNTriggerCircuit = 234;
-Float_t AliMUONConstants::fgDefaultChamberZ[14] =
-{533.5, 546.5, 678.5, 693.5, 964.0, 986.0, 1251.5, 1278.5, 1416.5, 1443.5,
- 1610, 1625., 1710., 1725.};
+Float_t AliMUONConstants::fgSt345inclination = 0.794; // in degrees, rotation axis is X axis
+
+Float_t AliMUONConstants::fgDzCh = 15.5/2.;
+Float_t AliMUONConstants::fgDzSlat = 8.5/2.;
+
+Float_t AliMUONConstants::fgSqrtKx3 = 0.7131;
+Float_t AliMUONConstants::fgSqrtKy3 = 0.7642;
+
+Float_t AliMUONConstants::fgSqrtKx3St1 = 0.7000;
+Float_t AliMUONConstants::fgSqrtKy3St1 = 0.7550;
+
+Float_t AliMUONConstants::fgChargeCorrel = 0.11;
+Float_t AliMUONConstants::fgChargeCorrelSt1 = 1.0; //???
+Float_t AliMUONConstants::fgPitch = 0.25;
+Float_t AliMUONConstants::fgPitchSt1 = 0.21;
+
+// From Alain TOURNAIRE
+// ALICE / ALICE Engineering baseline / Dimuonspectrometer (DIS) v7-1
+// EDMS Id 335328 for "search in EDMS
+// These are the diameter (Dmin == innner and Dmax - outner) values of the active surface
+// In the case of Dmax, the value corresponds to the maximum diameter of the active surface with 2pi coverture in phi
+Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 63.0, 79.0, 79.0, 99.0, 100.0}; // cm
+Float_t AliMUONConstants::fgDmax[7] = { 176.6, 229.0, 308.84, 418.2, 522.0, 850.0, 900.0}; // cm
+
+Int_t AliMUONConstants::fgMaxZoom = 20;
+
+// Defaults parameters for dipole magnet
+// From ALICE Dimuon - parameters / geometry table,
+// V7-3 (version 7 created 24/03/2004 updated 25/10/2005)
+Double_t AliMUONConstants::fgCoilZ = -994.05;
+Double_t AliMUONConstants::fgCoilL = 502.1;
+Double_t AliMUONConstants::fgYokeZ = -986.6;
+Double_t AliMUONConstants::fgYokeL = 309.4;
+
+// Defaults parameters for absorber (27/06/07)
+const Double_t AliMUONConstants::fgkAbsZBeg = -90.;
+const Double_t AliMUONConstants::fgkAbsZEnd = -505.;
+
+// Defaults parameters for track reconstruction
+Double_t AliMUONConstants::fgChamberThicknessInX0 = 0.03;
+
+//______________________________________________________________________________
+Int_t AliMUONConstants::NCh()
+{
+ return AliMpConstants::NofChambers();
+}
-Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 66.0, 80., 80., 100., 100.};
-Float_t AliMUONConstants::fgDmax[7] = {183., 245., 316.6, 560., 560., 850., 900.};
+//______________________________________________________________________________
+Int_t AliMUONConstants::NTrackingCh()
+{
+ return AliMpConstants::NofTrackingChambers();
+}
-Int_t AliMUONConstants::fgMaxZoom = 20;
+//______________________________________________________________________________
+Int_t AliMUONConstants::ChamberNumber(Float_t z)
+{
+ // return chamber number according z position of hit. Should be taken from geometry ?
+
+ Float_t dMaxChamber = DzSlat() + DzCh() + 0.25; // cm st 3 &4 & 5
+ dMaxChamber += 3.00; // factor for inclination of chamber
+ // dMaxChamber += Rmax(4) * TMath::Sin(fgSt345inclination*TMath::Pi()/360.);
+ // factor for inclination of chamber
+ if ( z > (DefaultChamberZ(4)+50.)) dMaxChamber = 7.; // cm stations 1 & 2
+ Int_t iChamber;
+ for (iChamber = 0; iChamber < 10; iChamber++) {
+ if (TMath::Abs(z-DefaultChamberZ(iChamber)) < dMaxChamber) {
+ return iChamber;
+ }
+ }
+
+ if ( z > DefaultChamberZ(NTrackingCh()-1) ) {
+ AliWarningClass(Form("No chamber number found for z = %f",z));
+ // for (iChamber = 0; iChamber < 10; iChamber++) {
+ // cout << iChamber << " zpos: " << DefaultChamberZ(iChamber)
+ // << " from " << DefaultChamberZ(iChamber) + dMaxChamber
+ // << " to " << DefaultChamberZ(iChamber) - dMaxChamber
+ // << " delta " << dMaxChamber << endl;
+ //}
+ }
+ return -1;
+}