Fix histo name
[u/mrichter/AliRoot.git] / MUON / AliMUONConstants.cxx
index af46b1c..c3997ab 100644 (file)
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/*
-$Log$
-Revision 1.8  2001/03/16 15:32:06  morsch
-Corrections of overlap with beam shield and dipole (A. de Falco)
+/* $Id$ */
 
-Revision 1.7  2001/03/14 17:22:15  pcrochet
-Geometry of the trigger chambers : a vertical gap of has been introduced around x=0 according fig.3.27 of the TDR (P.Dupieux)
+#include "AliMUONConstants.h"
 
-Revision 1.6  2001/01/30 12:19:39  morsch
-Update chamber positions (AdTDR version update 4/12/2000).
+#include "AliLog.h"
 
-Revision 1.5  2000/10/18 13:26:10  morsch
-New z-positions of chambers after Erice
+#include "TMath.h"
+#include "TClass.h"
+#include "AliMpConstants.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).
+//-----------------------------------------------------------------------------
+/// \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.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
+/// \cond CLASSIMP
+ClassImp(AliMUONConstants)
+/// \endcond
 
-Revision 1.2  2000/06/27 09:46:57  morsch
-kMAXZOOM global constant now in AliMUONConstants
+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;
 
-Revision 1.1  2000/06/26 14:02:38  morsch
-Add class AliMUONConstants with MUON specific constants using static memeber data and access methods.
+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
 
-*/
 
-#include "AliMUONConstants.h"
+// These are used by AliMUONConstants::ChamberNumber and must be calculated once
+// by that method from fgDzCh, fgDzSlat, fgDefaultChamberZ and fgSt345inclination,
+// so for now we set everything to zero.
+Float_t AliMUONConstants::fgDefaultChamberMinZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+Float_t AliMUONConstants::fgDefaultChamberMaxZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
 
+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,  98.8,  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.;
+    
+// Default trigger chamber resolution (cm)
+const Double_t AliMUONConstants::fgkTriggerNonBendingReso = 1.;
+const Double_t AliMUONConstants::fgkTriggerBendingReso = 1.;
+
+// Defaults parameters for muon filter (19/11/07)
+const Double_t AliMUONConstants::fgkMuonFilterZBeg = -1471.;
+const Double_t AliMUONConstants::fgkMuonFilterZEnd = -1471.-120.;
+const Double_t AliMUONConstants::fgkMuonFilterX0 = 1.76;
+
+// Defaults parameters for track reconstruction
+Double_t  AliMUONConstants::fgChamberThicknessInX0 = 0.03;
+
+const Double_t AliMUONConstants::fgkMostProbBendingMomentum = 2.0;
+
+//______________________________________________________________________________
+Int_t AliMUONConstants::NCh()
+{
+  return AliMpConstants::NofChambers();
+}
+
+//______________________________________________________________________________
+Int_t AliMUONConstants::NTrackingCh()
+{
+  return AliMpConstants::NofTrackingChambers();
+}
+
+//______________________________________________________________________________
+Int_t AliMUONConstants::ChamberNumber(Float_t z, bool warn)
+{
+  // return chamber number according z position of hit. Should be taken from geometry ?
+
+  if (fgDefaultChamberMinZ[0] == 0) // Are the min/max Z arrays initialised?
+  {
+    // The min and max Z arrays need to be calculated.
+    for (Int_t i = 0; i < NCh(); i++)
+    {
+      Float_t a = 0, b = 0;
+      if (4 <= i and i < 10)
+      {
+        Float_t dzAngle = TMath::Tan(TMath::Pi()*St345Inclination()/180.) * Rmax(i/2);
+        // We add 2.5mm since Rmax is an under-estimate.
+        a = DefaultChamberZ(i) + DzSlat() + DzCh() + dzAngle + 0.25;
+        b = DefaultChamberZ(i) - DzSlat() - DzCh() - dzAngle - 0.25;
+      }
+      else
+      {
+        a = DefaultChamberZ(i) + DzSlat();
+        b = DefaultChamberZ(i) - DzSlat();
+      }
+      fgDefaultChamberMinZ[i] = TMath::Min(a, b);
+      fgDefaultChamberMaxZ[i] = TMath::Max(a, b);
+    }
+  }
 
-Float_t  AliMUONConstants::fgDmin[7] = {  36.4,  46.2,  66.0,   80.,   80., 100., 100.};    
-Float_t  AliMUONConstants::fgDmax[7]  = {183., 245., 395.,  560.,  563., 850., 900.};  
-Int_t   AliMUONConstants::fgMaxZoom = 20;
+  // We can apply a binary search for the chamber since the fgDefaultChamberMinZ and
+  // fgDefaultChamberMaxZ arrays are ordered.
+  Int_t mini = 0, maxi = NCh()-1;
+  while (mini <= maxi)
+  {
+    Int_t iChamber = (maxi + mini) / 2;
+    if (z < fgDefaultChamberMinZ[iChamber])
+      mini = iChamber+1;
+    else if (z > fgDefaultChamberMaxZ[iChamber])
+      maxi = iChamber-1;
+    else
+      // We are between min and max Z of chamber number iChamber so we found our chamber.
+      return iChamber;
+  }
 
+  if (warn) AliWarningClass(Form("No chamber number found for z = %f",z));
+  return -1;
+}