Line data Source code
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 18 : 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 : // Warning: the resolution refers only to ALIGNMENT
103 : // For the total resolution the strip width should be taken into account!
104 : const Double_t AliMUONConstants::fgkTriggerNonBendingReso = 0.2;
105 : const Double_t AliMUONConstants::fgkTriggerBendingReso = 0.2;
106 :
107 : // Defaults parameters for muon filter (19/11/07)
108 : const Double_t AliMUONConstants::fgkMuonFilterZBeg = -1471.;
109 : const Double_t AliMUONConstants::fgkMuonFilterZEnd = -1471.-120.;
110 : const Double_t AliMUONConstants::fgkMuonFilterX0 = 1.76;
111 :
112 : // Defaults parameters for track reconstruction
113 : Double_t AliMUONConstants::fgChamberThicknessInX0[10] = {0.065, 0.065, 0.075, 0.075, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035};
114 :
115 : const Double_t AliMUONConstants::fgkMostProbBendingMomentum = 2.0;
116 :
117 : Float_t AliMUONConstants::fgAverageChamberT[14]=
118 : {17.64*1E-9, 18.28*1E-9, 22.68*1E-9, 23.33*1E-9, 32.42*1E-9, 33.48*1E-9, 42.76*1E-9,
119 : 43.81*1E-9, 47.13*1E-9, 48.17*1E-9, 53.75*1E-9, 54.32*1E-9, 57.12*1E-9, 57.67*1E-9};
120 :
121 : // is equivalent to gain = 4 mV/fC
122 : const Float_t AliMUONConstants::fgkDefaultA0 = 1.25;
123 : // 1 ADC channel = 0.61 mV
124 : const Float_t AliMUONConstants::fgkDefaultADC2MV = 0.61;
125 : const Float_t AliMUONConstants::fgkDefaultCapa = 0.2;
126 :
127 : //______________________________________________________________________________
128 : Int_t AliMUONConstants::NCh()
129 : {
130 161494 : return AliMpConstants::NofChambers();
131 : }
132 :
133 : //______________________________________________________________________________
134 : Int_t AliMUONConstants::NTrackingCh()
135 : {
136 162212 : return AliMpConstants::NofTrackingChambers();
137 : }
138 :
139 : //______________________________________________________________________________
140 : Int_t AliMUONConstants::ChamberNumber(Float_t z, bool warn)
141 : {
142 : // return chamber number according z position of hit. Should be taken from geometry ?
143 :
144 0 : if (fgDefaultChamberMinZ[0] == 0) // Are the min/max Z arrays initialised?
145 : {
146 : // The min and max Z arrays need to be calculated.
147 0 : for (Int_t i = 0; i < NCh(); i++)
148 : {
149 : Float_t a = 0, b = 0;
150 0 : if (4 <= i and i < 10)
151 : {
152 0 : Float_t dzAngle = TMath::Tan(TMath::Pi()*St345Inclination()/180.) * Rmax(i/2);
153 : // We add 2.5mm since Rmax is an under-estimate.
154 0 : a = DefaultChamberZ(i) + DzSlat() + DzCh() + dzAngle + 0.25;
155 0 : b = DefaultChamberZ(i) - DzSlat() - DzCh() - dzAngle - 0.25;
156 0 : }
157 : else
158 : {
159 0 : a = DefaultChamberZ(i) + DzSlat();
160 0 : b = DefaultChamberZ(i) - DzSlat();
161 : }
162 0 : fgDefaultChamberMinZ[i] = TMath::Min(a, b);
163 0 : fgDefaultChamberMaxZ[i] = TMath::Max(a, b);
164 : }
165 0 : }
166 :
167 : // We can apply a binary search for the chamber since the fgDefaultChamberMinZ and
168 : // fgDefaultChamberMaxZ arrays are ordered.
169 0 : Int_t mini = 0, maxi = NCh()-1;
170 0 : while (mini <= maxi)
171 : {
172 0 : Int_t iChamber = (maxi + mini) / 2;
173 0 : if (z < fgDefaultChamberMinZ[iChamber])
174 0 : mini = iChamber+1;
175 0 : else if (z > fgDefaultChamberMaxZ[iChamber])
176 0 : maxi = iChamber-1;
177 : else
178 : // We are between min and max Z of chamber number iChamber so we found our chamber.
179 0 : return iChamber;
180 0 : }
181 :
182 0 : if (warn) AliWarningClass(Form("No chamber number found for z = %f",z));
183 0 : return -1;
184 0 : }
185 :
186 : //______________________________________________________________________________
187 : Float_t AliMUONConstants::ReducedQTot(Float_t qtot, Float_t timeDif)
188 : {
189 : // return a reduced charge if the hit belongs to a track from a pileup event
190 0 : Float_t q = qtot*1.19*(1.24-timeDif*1E6)*TMath::Exp(-(0.97-timeDif*1E6)*(0.97-timeDif*1E6)/2.42);
191 0 : return q;
192 : }
193 :
194 : //______________________________________________________________________________
195 : Float_t AliMUONConstants::FC2ADC() {
196 : // Return conversionfactor fc to adc
197 0 : Float_t fc2adc = 1./(DefaultA0()*DefaultCapa()*DefaultADC2MV());
198 :
199 0 : return fc2adc;
200 : }
|
