]>
Commit | Line | Data |
---|---|---|
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 | Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] = | |
53 | {1., 1.01060, 1.06236, 1.07296}; | |
54 | ||
55 | ||
56 | Float_t AliMUONConstants::fgSt345inclination = 0.794; // in degrees, rotation axis is X axis | |
57 | ||
58 | Float_t AliMUONConstants::fgDzCh = 15.5/2.; | |
59 | Float_t AliMUONConstants::fgDzSlat = 8.5/2.; | |
60 | ||
61 | Float_t AliMUONConstants::fgSqrtKx3 = 0.7131; | |
62 | Float_t AliMUONConstants::fgSqrtKy3 = 0.7642; | |
63 | ||
64 | Float_t AliMUONConstants::fgSqrtKx3St1 = 0.7000; | |
65 | Float_t AliMUONConstants::fgSqrtKy3St1 = 0.7550; | |
66 | ||
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; | |
71 | ||
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 | |
79 | ||
80 | Int_t AliMUONConstants::fgMaxZoom = 20; | |
81 | ||
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; | |
89 | ||
90 | // Defaults parameters for absorber (27/06/07) | |
91 | const Double_t AliMUONConstants::fgkAbsZBeg = -90.; | |
92 | const Double_t AliMUONConstants::fgkAbsZEnd = -505.; | |
93 | ||
94 | // Defaults parameters for track reconstruction | |
95 | Double_t AliMUONConstants::fgChamberThicknessInX0 = 0.03; | |
96 | ||
97 | //______________________________________________________________________________ | |
98 | Int_t AliMUONConstants::NCh() | |
99 | { | |
100 | return AliMpConstants::NofChambers(); | |
101 | } | |
102 | ||
103 | //______________________________________________________________________________ | |
104 | Int_t AliMUONConstants::NTrackingCh() | |
105 | { | |
106 | return AliMpConstants::NofTrackingChambers(); | |
107 | } | |
108 | ||
109 | //______________________________________________________________________________ | |
110 | Int_t AliMUONConstants::ChamberNumber(Float_t z) | |
111 | { | |
112 | // return chamber number according z position of hit. Should be taken from geometry ? | |
113 | ||
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 | |
119 | Int_t iChamber; | |
120 | ||
121 | for (iChamber = 0; iChamber < 10; iChamber++) { | |
122 | if (TMath::Abs(z-DefaultChamberZ(iChamber)) < dMaxChamber) { | |
123 | return iChamber; | |
124 | } | |
125 | } | |
126 | ||
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; | |
134 | //} | |
135 | } | |
136 | return -1; | |
137 | } |