1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 #include "AliMUONConstants.h"
24 #include "AliMpConstants.h"
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,
35 //-----------------------------------------------------------------------------
38 ClassImp(AliMUONConstants)
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;
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
52 Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] =
53 {1., 1.01060, 1.06236, 1.07296};
56 Float_t AliMUONConstants::fgSt345inclination = 0.794; // in degrees, rotation axis is X axis
58 Float_t AliMUONConstants::fgDzCh = 15.5/2.;
59 Float_t AliMUONConstants::fgDzSlat = 8.5/2.;
61 Float_t AliMUONConstants::fgSqrtKx3 = 0.7131;
62 Float_t AliMUONConstants::fgSqrtKy3 = 0.7642;
64 Float_t AliMUONConstants::fgSqrtKx3St1 = 0.7000;
65 Float_t AliMUONConstants::fgSqrtKy3St1 = 0.7550;
67 Float_t AliMUONConstants::fgChargeCorrel = 0.11;
68 Float_t AliMUONConstants::fgChargeCorrelSt1 = 1.0; //???
69 Float_t AliMUONConstants::fgPitch = 0.25;
70 Float_t AliMUONConstants::fgPitchSt1 = 0.21;
72 // From Alain TOURNAIRE
73 // ALICE / ALICE Engineering baseline / Dimuonspectrometer (DIS) v7-1
74 // EDMS Id 335328 for "search in EDMS
75 // These are the diameter (Dmin == innner and Dmax - outner) values of the active surface
76 // In the case of Dmax, the value corresponds to the maximum diameter of the active surface with 2pi coverture in phi
77 Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 63.0, 79.0, 79.0, 99.0, 100.0}; // cm
78 Float_t AliMUONConstants::fgDmax[7] = { 176.6, 229.0, 308.84, 418.2, 522.0, 850.0, 900.0}; // cm
80 Int_t AliMUONConstants::fgMaxZoom = 20;
82 // Defaults parameters for dipole magnet
83 // From ALICE Dimuon - parameters / geometry table,
84 // V7-3 (version 7 created 24/03/2004 updated 25/10/2005)
85 Double_t AliMUONConstants::fgCoilZ = -994.05;
86 Double_t AliMUONConstants::fgCoilL = 502.1;
87 Double_t AliMUONConstants::fgYokeZ = -986.6;
88 Double_t AliMUONConstants::fgYokeL = 309.4;
90 // Defaults parameters for absorber (27/06/07)
91 const Double_t AliMUONConstants::fgkAbsZBeg = -90.;
92 const Double_t AliMUONConstants::fgkAbsZEnd = -505.;
94 // Defaults parameters for track reconstruction
95 Double_t AliMUONConstants::fgChamberThicknessInX0 = 0.03;
97 //______________________________________________________________________________
98 Int_t AliMUONConstants::NCh()
100 return AliMpConstants::NofChambers();
103 //______________________________________________________________________________
104 Int_t AliMUONConstants::NTrackingCh()
106 return AliMpConstants::NofTrackingChambers();
109 //______________________________________________________________________________
110 Int_t AliMUONConstants::ChamberNumber(Float_t z)
112 // return chamber number according z position of hit. Should be taken from geometry ?
114 Float_t dMaxChamber = DzSlat() + DzCh() + 0.25; // cm st 3 &4 & 5
115 dMaxChamber += 3.00; // factor for inclination of chamber
116 // dMaxChamber += Rmax(4) * TMath::Sin(fgSt345inclination*TMath::Pi()/360.);
117 // factor for inclination of chamber
118 if ( z > (DefaultChamberZ(4)+50.)) dMaxChamber = 7.; // cm stations 1 & 2
121 for (iChamber = 0; iChamber < 10; iChamber++) {
122 if (TMath::Abs(z-DefaultChamberZ(iChamber)) < dMaxChamber) {
127 if ( z > DefaultChamberZ(NTrackingCh()-1) ) {
128 AliWarningClass(Form("No chamber number found for z = %f",z));
129 // for (iChamber = 0; iChamber < 10; iChamber++) {
130 // cout << iChamber << " zpos: " << DefaultChamberZ(iChamber)
131 // << " from " << DefaultChamberZ(iChamber) + dMaxChamber
132 // << " to " << DefaultChamberZ(iChamber) - dMaxChamber
133 // << " delta " << dMaxChamber << endl;