include fgkTriggerTofLimit (Philippe)

[u/mrichter/AliRoot.git] / MUON / AliMUONConstants.cxx

diff --git a/MUON/AliMUONConstants.cxx b/MUON/AliMUONConstants.cxx
index e66a0edbb6f58e0bb6ed0da3299c10131430e4f9..84fcf96df5c308d8f337ec707bdab3199f5408dc 100644 (file)
--- a/MUON/AliMUONConstants.cxx
+++ b/MUON/AliMUONConstants.cxx
@@ -17,43 +17,94 @@
 
 #include "AliMUONConstants.h"
 
+#include "AliLog.h"
 
+#include "TMath.h"
+#include "TClass.h"
+
+///
+/// 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).
+///
+
+/// \cond CLASSIMP
 ClassImp(AliMUONConstants)
+/// \endcond
 
 Int_t   AliMUONConstants::fgNCh = 14;
 Int_t   AliMUONConstants::fgNTrackingCh = 10;
 Int_t   AliMUONConstants::fgNTriggerCh = 4;
 Int_t   AliMUONConstants::fgNTriggerCircuit = 234;
-Int_t   AliMUONConstants::fgNofDetElements[14] =
-{ 4, 4, 4, 4, 18, 18, 26, 26, 26, 26, 18, 18, 18, 18 };
+Int_t   AliMUONConstants::fgNDetElem = 228;
+Int_t   AliMUONConstants::fgNGeomModules = 20;
+Float_t AliMUONConstants::fgkTriggerTofLimit = 75E-9;
 
 Float_t AliMUONConstants::fgDefaultChamberZ[14] = 
-  {-533.5, -546.5, -678.5, -693.5, // St12
-   -966.9, -997.9, -1274.5, -1305.5, -1408.6, -1439.6, // St345  update sept04 Ch. Finck 
+  {-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};
+
+
+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::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;
+Float_t  AliMUONConstants::fgSqrtKx3St1 = 0.7000;
+Float_t  AliMUONConstants::fgSqrtKy3St1 = 0.7550;
 
-//______________________________________________________________________________
-Int_t AliMUONConstants::GetChamberId(Int_t detElemId)
-{ 
-// Get chamber Id from detection element Id
-// ---
+Float_t  AliMUONConstants::fgChargeCorrel    = 0.11;
+Float_t  AliMUONConstants::fgChargeCorrelSt1 = 1.0; //??? 
+Float_t  AliMUONConstants::fgPitch     = 0.25;
+Float_t  AliMUONConstants::fgPitchSt1  = 0.21; 
 
-  return detElemId/100 - 1;
-}  
+// 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;
 
 //______________________________________________________________________________
-Int_t AliMUONConstants::GetFirstDetElemId(Int_t chamberId)
+Int_t AliMUONConstants::ChamberNumber(Float_t z) 
 {
-// Get first detection element Id for chamber specified by chamber Id
-// ---
+  // 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;
 
-  return (chamberId+1)*100;
-}  
+  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;
+}