bug fixed
[u/mrichter/AliRoot.git] / MUON / AliMUONConstants.cxx
CommitLineData
f665c1ea 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
88cb7938 16/* $Id$ */
f665c1ea 17
18#include "AliMUONConstants.h"
19
fe126bce 20#include "AliLog.h"
21
9d5f6a64 22#include "TMath.h"
fe126bce 23#include "TClass.h"
6061eac6 24#include "AliMpConstants.h"
9d5f6a64 25
3d1463c8 26//-----------------------------------------------------------------------------
c4ee792d 27/// \class AliMUONConstants
9d5f6a64 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).
3d1463c8 35//-----------------------------------------------------------------------------
f665c1ea 36
5398f946 37/// \cond CLASSIMP
38ClassImp(AliMUONConstants)
39/// \endcond
40
de2cd600 41Int_t AliMUONConstants::fgNTrackingSt = 5;
f665c1ea 42Int_t AliMUONConstants::fgNTriggerCircuit = 234;
edb96c6a 43Int_t AliMUONConstants::fgNDetElem = 228;
573dc4ee 44Int_t AliMUONConstants::fgNGeomModules = 20;
b717b957 45Float_t AliMUONConstants::fgkTriggerTofLimit = 75E-9;
89b0b88d 46
e516b01d 47Float_t AliMUONConstants::fgDefaultChamberZ[14] =
89b0b88d 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)
e516b01d 50 -1603.5, -1620.5, -1703.5, -1720.5}; // M1 & M2
51
3787f517 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.
56Float_t AliMUONConstants::fgDefaultChamberMinZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
57Float_t AliMUONConstants::fgDefaultChamberMaxZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
58
43adb63d 59Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] =
60{1., 1.01060, 1.06236, 1.07296};
61
62
2c9844e7 63Float_t AliMUONConstants::fgSt345inclination = 0.794; // in degrees, rotation axis is X axis
64
e516b01d 65Float_t AliMUONConstants::fgDzCh = 15.5/2.;
66Float_t AliMUONConstants::fgDzSlat = 8.5/2.;
67
1c334adf 68Float_t AliMUONConstants::fgSqrtKx3 = 0.7131;
69Float_t AliMUONConstants::fgSqrtKy3 = 0.7642;
7e4a628d 70
1c334adf 71Float_t AliMUONConstants::fgSqrtKx3St1 = 0.7000;
72Float_t AliMUONConstants::fgSqrtKy3St1 = 0.7550;
7e4a628d 73
1c334adf 74Float_t AliMUONConstants::fgChargeCorrel = 0.11;
f4a7360f 75Float_t AliMUONConstants::fgChargeCorrelSt1 = 1.0; //???
1c334adf 76Float_t AliMUONConstants::fgPitch = 0.25;
77Float_t AliMUONConstants::fgPitchSt1 = 0.21;
c47b11f0 78
f4a7360f 79// From Alain TOURNAIRE
80// ALICE / ALICE Engineering baseline / Dimuonspectrometer (DIS) v7-1
81// EDMS Id 335328 for "search in EDMS
124dd59f 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
32e13b89 84Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 63.0, 79.0, 79.0, 98.8, 100.0}; // cm
124dd59f 85Float_t AliMUONConstants::fgDmax[7] = { 176.6, 229.0, 308.84, 418.2, 522.0, 850.0, 900.0}; // cm
f4a7360f 86
e118b27e 87Int_t AliMUONConstants::fgMaxZoom = 20;
ecfa008b 88
208f139e 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)
92Double_t AliMUONConstants::fgCoilZ = -994.05;
93Double_t AliMUONConstants::fgCoilL = 502.1;
94Double_t AliMUONConstants::fgYokeZ = -986.6;
95Double_t AliMUONConstants::fgYokeL = 309.4;
96
ff71f4a2 97// Defaults parameters for absorber (27/06/07)
98const Double_t AliMUONConstants::fgkAbsZBeg = -90.;
99const Double_t AliMUONConstants::fgkAbsZEnd = -505.;
100
fda59e58 101// Default trigger chamber resolution (cm)
bdfb6eef 102// Warning: the resolution refers only to ALIGNMENT
103// For the total resolution the strip width should be taken into account!
104const Double_t AliMUONConstants::fgkTriggerNonBendingReso = 0.2;
105const Double_t AliMUONConstants::fgkTriggerBendingReso = 0.2;
fda59e58 106
107// Defaults parameters for muon filter (19/11/07)
108const Double_t AliMUONConstants::fgkMuonFilterZBeg = -1471.;
109const Double_t AliMUONConstants::fgkMuonFilterZEnd = -1471.-120.;
110const Double_t AliMUONConstants::fgkMuonFilterX0 = 1.76;
111
de2cd600 112// Defaults parameters for track reconstruction
4663da9f 113Double_t AliMUONConstants::fgChamberThicknessInX0[10] = {0.065, 0.065, 0.075, 0.075, 0.035, 0.035, 0.035, 0.035, 0.035, 0.035};
de2cd600 114
9bdbee64 115const Double_t AliMUONConstants::fgkMostProbBendingMomentum = 2.0;
116
2a5f75ae 117Float_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
24c9eaa0 121//______________________________________________________________________________
6061eac6 122Int_t AliMUONConstants::NCh()
123{
124 return AliMpConstants::NofChambers();
125}
126
127//______________________________________________________________________________
128Int_t AliMUONConstants::NTrackingCh()
129{
130 return AliMpConstants::NofTrackingChambers();
131}
132
133//______________________________________________________________________________
3787f517 134Int_t AliMUONConstants::ChamberNumber(Float_t z, bool warn)
29fc2c86 135{
136 // return chamber number according z position of hit. Should be taken from geometry ?
3787f517 137
138 if (fgDefaultChamberMinZ[0] == 0) // Are the min/max Z arrays initialised?
139 {
140 // The min and max Z arrays need to be calculated.
141 for (Int_t i = 0; i < NCh(); i++)
142 {
143 Float_t a = 0, b = 0;
144 if (4 <= i and i < 10)
145 {
146 Float_t dzAngle = TMath::Tan(TMath::Pi()*St345Inclination()/180.) * Rmax(i/2);
147 // We add 2.5mm since Rmax is an under-estimate.
148 a = DefaultChamberZ(i) + DzSlat() + DzCh() + dzAngle + 0.25;
149 b = DefaultChamberZ(i) - DzSlat() - DzCh() - dzAngle - 0.25;
150 }
151 else
152 {
153 a = DefaultChamberZ(i) + DzSlat();
154 b = DefaultChamberZ(i) - DzSlat();
155 }
156 fgDefaultChamberMinZ[i] = TMath::Min(a, b);
157 fgDefaultChamberMaxZ[i] = TMath::Max(a, b);
29fc2c86 158 }
159 }
3787f517 160
161 // We can apply a binary search for the chamber since the fgDefaultChamberMinZ and
162 // fgDefaultChamberMaxZ arrays are ordered.
163 Int_t mini = 0, maxi = NCh()-1;
164 while (mini <= maxi)
165 {
166 Int_t iChamber = (maxi + mini) / 2;
167 if (z < fgDefaultChamberMinZ[iChamber])
168 mini = iChamber+1;
169 else if (z > fgDefaultChamberMaxZ[iChamber])
170 maxi = iChamber-1;
171 else
172 // We are between min and max Z of chamber number iChamber so we found our chamber.
173 return iChamber;
fe126bce 174 }
3787f517 175
176 if (warn) AliWarningClass(Form("No chamber number found for z = %f",z));
29fc2c86 177 return -1;
178}
2a5f75ae 179
180//______________________________________________________________________________
181Float_t AliMUONConstants::ReducedQTot(Float_t qtot, Float_t timeDif)
182{
183 // return a reduced charge if the hit belongs to a track from a pileup event
184 Float_t q = qtot*1.19*(1.24-timeDif*1E6)*TMath::Exp(-(0.97-timeDif*1E6)*(0.97-timeDif*1E6)/2.42);
185 return q;
186}