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 #include "AliMpConstants.h"
  25
  26 //-----------------------------------------------------------------------------
  27 /// \class AliMUONConstants
  28 /// This class holds various constants to be used in many places,
  29 /// such as the number of tracking and trigger chambers,
  30 /// some geometrical constants (to build the initial geometry for simulation)
  31 /// and mathieson distribution default values.
  32 /// Those constants should as much as possible replace hard-coded values
  33 /// which are to be considered strictly illegal in the MUON code (or any code,
  34 /// by the way).
  35 //-----------------------------------------------------------------------------
  36
  37 /// \cond CLASSIMP
  38 ClassImp(AliMUONConstants)
  39 /// \endcond
  40
  41 Int_t   AliMUONConstants::fgNTrackingSt = 5;
  42 Int_t   AliMUONConstants::fgNTriggerCircuit = 234;
  43 Int_t   AliMUONConstants::fgNDetElem = 228;
  44 Int_t   AliMUONConstants::fgNGeomModules = 20;
  45 Float_t AliMUONConstants::fgkTriggerTofLimit = 75E-9;
  46
  47 Float_t AliMUONConstants::fgDefaultChamberZ[14] =
  48   {-526.16, -545.24, -676.4, -695.4, // St12
  49    -967.5, -998.5, -1276.5, -1307.5, -1406.6, -1437.6,// updated 08/05, EDMS id 335328 (A. Tournaire)
  50    -1603.5, -1620.5, -1703.5, -1720.5}; // M1 & M2
  51
  52
  53 // These are used by AliMUONConstants::ChamberNumber and must be calculated once
  54 // by that method from fgDzCh, fgDzSlat, fgDefaultChamberZ and fgSt345inclination,
  55 // so for now we set everything to zero.
  56 Float_t AliMUONConstants::fgDefaultChamberMinZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
  57 Float_t AliMUONConstants::fgDefaultChamberMaxZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
  58
  59 Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] =
  60 {1., 1.01060, 1.06236, 1.07296};
  61
  62
  63 Float_t  AliMUONConstants::fgSt345inclination = 0.794; // in degrees, rotation axis is X axis
  64
  65 Float_t  AliMUONConstants::fgDzCh   = 15.5/2.;
  66 Float_t  AliMUONConstants::fgDzSlat = 8.5/2.;
  67
  68 Float_t  AliMUONConstants::fgSqrtKx3 = 0.7131;
  69 Float_t  AliMUONConstants::fgSqrtKy3 = 0.7642;
  70
  71 Float_t  AliMUONConstants::fgSqrtKx3St1 = 0.7000;
  72 Float_t  AliMUONConstants::fgSqrtKy3St1 = 0.7550;
  73
  74 Float_t  AliMUONConstants::fgChargeCorrel    = 0.11;
  75 Float_t  AliMUONConstants::fgChargeCorrelSt1 = 1.0; //???
  76 Float_t  AliMUONConstants::fgPitch     = 0.25;
  77 Float_t  AliMUONConstants::fgPitchSt1  = 0.21;
  78
  79 // From Alain TOURNAIRE
  80 // ALICE / ALICE Engineering baseline / Dimuonspectrometer (DIS) v7-1
  81 // EDMS Id 335328 for "search in EDMS
  82 // These are the diameter (Dmin == innner and Dmax - outner) values of the active surface
  83 // In the case of Dmax, the value corresponds to the maximum diameter of the active surface with 2pi coverture in phi
  84 Float_t  AliMUONConstants::fgDmin[7]  = {   36.4,  46.2,  63.0,   79.0,   79.0,  98.8,  100.0};  // cm
  85 Float_t  AliMUONConstants::fgDmax[7]  = {  176.6, 229.0, 308.84, 418.2,  522.0, 850.0, 900.0};   // cm
  86
  87 Int_t    AliMUONConstants::fgMaxZoom = 20;
  88
  89 // Defaults parameters for dipole magnet
  90 // From ALICE Dimuon - parameters / geometry table,
  91 // V7-3 (version 7 created 24/03/2004 updated 25/10/2005)
  92 Double_t AliMUONConstants::fgCoilZ = -994.05;
  93 Double_t AliMUONConstants::fgCoilL = 502.1;
  94 Double_t AliMUONConstants::fgYokeZ = -986.6;
  95 Double_t AliMUONConstants::fgYokeL = 309.4;
  96
  97 // Defaults parameters for absorber (27/06/07)
  98 const Double_t AliMUONConstants::fgkAbsZBeg = -90.;
  99 const Double_t AliMUONConstants::fgkAbsZEnd = -505.;
 100
 101 // Default trigger chamber resolution (cm)
 102 const Double_t AliMUONConstants::fgkTriggerNonBendingReso = 1.;
 103 const Double_t AliMUONConstants::fgkTriggerBendingReso = 1.;
 104
 105 // Defaults parameters for muon filter (19/11/07)
 106 const Double_t AliMUONConstants::fgkMuonFilterZBeg = -1471.;
 107 const Double_t AliMUONConstants::fgkMuonFilterZEnd = -1471.-120.;
 108 const Double_t AliMUONConstants::fgkMuonFilterX0 = 1.76;
 109
 110 // Defaults parameters for track reconstruction
 111 Double_t  AliMUONConstants::fgChamberThicknessInX0 = 0.03;
 112
 113 const Double_t AliMUONConstants::fgkMostProbBendingMomentum = 2.0;
 114
 115 //______________________________________________________________________________
 116 Int_t AliMUONConstants::NCh()
 117 {
 118   return AliMpConstants::NofChambers();
 119 }
 120
 121 //______________________________________________________________________________
 122 Int_t AliMUONConstants::NTrackingCh()
 123 {
 124   return AliMpConstants::NofTrackingChambers();
 125 }
 126
 127 //______________________________________________________________________________
 128 Int_t AliMUONConstants::ChamberNumber(Float_t z, bool warn)
 129 {
 130   // return chamber number according z position of hit. Should be taken from geometry ?
 131
 132   if (fgDefaultChamberMinZ[0] == 0) // Are the min/max Z arrays initialised?
 133   {
 134     // The min and max Z arrays need to be calculated.
 135     for (Int_t i = 0; i < NCh(); i++)
 136     {
 137       Float_t a = 0, b = 0;
 138       if (4 <= i and i < 10)
 139       {
 140         Float_t dzAngle = TMath::Tan(TMath::Pi()*St345Inclination()/180.) * Rmax(i/2);
 141         // We add 2.5mm since Rmax is an under-estimate.
 142         a = DefaultChamberZ(i) + DzSlat() + DzCh() + dzAngle + 0.25;
 143         b = DefaultChamberZ(i) - DzSlat() - DzCh() - dzAngle - 0.25;
 144       }
 145       else
 146       {
 147         a = DefaultChamberZ(i) + DzSlat();
 148         b = DefaultChamberZ(i) - DzSlat();
 149       }
 150       fgDefaultChamberMinZ[i] = TMath::Min(a, b);
 151       fgDefaultChamberMaxZ[i] = TMath::Max(a, b);
 152     }
 153   }
 154
 155   // We can apply a binary search for the chamber since the fgDefaultChamberMinZ and
 156   // fgDefaultChamberMaxZ arrays are ordered.
 157   Int_t mini = 0, maxi = NCh()-1;
 158   while (mini <= maxi)
 159   {
 160     Int_t iChamber = (maxi + mini) / 2;
 161     if (z < fgDefaultChamberMinZ[iChamber])
 162       mini = iChamber+1;
 163     else if (z > fgDefaultChamberMaxZ[iChamber])
 164       maxi = iChamber-1;
 165     else
 166       // We are between min and max Z of chamber number iChamber so we found our chamber.
 167       return iChamber;
 168   }
 169
 170   if (warn) AliWarningClass(Form("No chamber number found for z = %f",z));
 171   return -1;
 172 }