MUON/AliMUONConstants.cxx

   1  /**************************************************************************
   2  * Copyright(c) 1998-1999, 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
  16 /* $Id$ */
  17
  18 #include "AliMUONConstants.h"
  19
  20 #include "AliLog.h"
  21
  22 #include "TMath.h"
  23 #include "TClass.h"
  24
  25 ///
  26 /// This class holds various constants to be used in many places,
  27 /// such as the number of tracking and trigger chambers,
  28 /// some geometrical constants (to build the initial geometry for simulation)
  29 /// and mathieson distribution default values.
  30 /// Those constants should as much as possible replace hard-coded values
  31 /// which are to be considered strictly illegal in the MUON code (or any code,
  32 /// by the way).
  33 ///
  34
  35 /// \cond CLASSIMP
  36 ClassImp(AliMUONConstants)
  37 /// \endcond
  38
  39 Int_t   AliMUONConstants::fgNCh = 14;
  40 Int_t   AliMUONConstants::fgNTrackingCh = 10;
  41 Int_t   AliMUONConstants::fgNTrackingSt = 5;
  42 Int_t   AliMUONConstants::fgNTriggerCh = 4;
  43 Int_t   AliMUONConstants::fgNTriggerCircuit = 234;
  44 Int_t   AliMUONConstants::fgNDetElem = 228;
  45 Int_t   AliMUONConstants::fgNGeomModules = 20;
  46 Float_t AliMUONConstants::fgkTriggerTofLimit = 75E-9;
  47
  48 Float_t AliMUONConstants::fgDefaultChamberZ[14] =
  49   {-526.16, -545.24, -676.4, -695.4, // St12
  50    -967.5, -998.5, -1276.5, -1307.5, -1406.6, -1437.6,// updated 08/05, EDMS id 335328 (A. Tournaire)
  51    -1603.5, -1620.5, -1703.5, -1720.5}; // M1 & M2
  52
  53 Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] =
  54 {1., 1.01060, 1.06236, 1.07296};
  55
  56
  57 Float_t  AliMUONConstants::fgSt345inclination = 0.794; // in degrees, rotation axis is X axis
  58
  59 Float_t  AliMUONConstants::fgDzCh   = 15.5/2.;
  60 Float_t  AliMUONConstants::fgDzSlat = 8.5/2.;
  61
  62 Float_t  AliMUONConstants::fgSqrtKx3 = 0.7131;
  63 Float_t  AliMUONConstants::fgSqrtKy3 = 0.7642;
  64
  65 Float_t  AliMUONConstants::fgSqrtKx3St1 = 0.7000;
  66 Float_t  AliMUONConstants::fgSqrtKy3St1 = 0.7550;
  67
  68 Float_t  AliMUONConstants::fgChargeCorrel    = 0.11;
  69 Float_t  AliMUONConstants::fgChargeCorrelSt1 = 1.0; //???
  70 Float_t  AliMUONConstants::fgPitch     = 0.25;
  71 Float_t  AliMUONConstants::fgPitchSt1  = 0.21;
  72
  73 // From Alain TOURNAIRE
  74 // ALICE / ALICE Engineering baseline / Dimuonspectrometer (DIS) v7-1
  75 // EDMS Id 335328 for "search in EDMS
  76 // These are the diameter (Dmin == innner and Dmax - outner) values of the active surface
  77 // In the case of Dmax, the value corresponds to the maximum diameter of the active surface with 2pi coverture in phi
  78 Float_t  AliMUONConstants::fgDmin[7]  = {   36.4,  46.2,  63.0,   79.0,   79.0,  99.0,  100.0};  // cm
  79 Float_t  AliMUONConstants::fgDmax[7]  = {  176.6, 229.0, 308.84, 418.2,  522.0, 850.0, 900.0};   // cm
  80
  81 Int_t    AliMUONConstants::fgMaxZoom = 20;
  82
  83 // Defaults parameters for track reconstruction
  84 Double_t  AliMUONConstants::fgDefaultChamberThicknessInX0 = 0.03;
  85
  86 //______________________________________________________________________________
  87 Int_t AliMUONConstants::ChamberNumber(Float_t z)
  88 {
  89   // return chamber number according z position of hit. Should be taken from geometry ?
  90
  91   Float_t dMaxChamber = DzSlat() + DzCh() + 0.25; // cm st 3 &4 & 5
  92   dMaxChamber += 3.00;                            // factor for inclination of chamber
  93   // dMaxChamber += Rmax(4) * TMath::Sin(fgSt345inclination*TMath::Pi()/360.);
  94                                                   // factor for inclination of chamber
  95   if ( z >  (DefaultChamberZ(4)+50.)) dMaxChamber = 7.; // cm stations 1 & 2
  96   Int_t iChamber;
  97
  98   for (iChamber = 0; iChamber < 10; iChamber++) {
  99     if (TMath::Abs(z-DefaultChamberZ(iChamber)) < dMaxChamber) {
 100       return iChamber;
 101     }
 102   }
 103
 104   if ( z > DefaultChamberZ(NTrackingCh()-1) ) {
 105     AliWarningClass(Form("No chamber number found for z = %f",z));
 106     // for (iChamber = 0; iChamber < 10; iChamber++) {
 107     //   cout << iChamber << " zpos: " << DefaultChamberZ(iChamber)
 108     //        << "  from " << DefaultChamberZ(iChamber) + dMaxChamber
 109     //        << "  to " << DefaultChamberZ(iChamber) - dMaxChamber
 110     //        << "  delta " << dMaxChamber << endl;
 111     //}
 112   }
 113   return -1;
 114 }