1 #ifndef ALIHLTMUONCALCULATIONS_H
2 #define ALIHLTMUONCALCULATIONS_H
3 /* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 ////////////////////////////////////////////////////////////////////////////////
10 // Author: Artur Szostak
11 // Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
13 ////////////////////////////////////////////////////////////////////////////////
15 #include "AliHLTMUONDataTypes.h"
17 extern "C" struct AliHLTMUONTriggerRecordStruct;
20 * Note: this class uses static global variables so thread protection must be
21 * explicit in any multi-threaded usage. Or the class should be rewritten.
23 class AliHLTMUONCalculations
27 /// Calculates the momentum estimate given two track points behind the
28 /// dipole magnet and assuming origin is the interaction point.
29 static bool ComputeMomentum(
31 AliHLTFloat32_t y1, AliHLTFloat32_t y2,
32 AliHLTFloat32_t z1, AliHLTFloat32_t z2
35 static AliHLTFloat32_t Zf() { return fgZf; }
36 static void Zf(AliHLTFloat32_t value) { fgZf = value; }
37 static AliHLTFloat32_t QBL();
38 static void QBL(AliHLTFloat32_t value);
40 static AliHLTMUONParticleSign Sign() { return fgSign; }
41 static AliHLTFloat32_t Px() { return fgPx; }
42 static AliHLTFloat32_t Py() { return fgPy; }
43 static AliHLTFloat32_t Pz() { return fgPz; }
45 /// Calculates the invariant mass for a pair of particles.
46 static AliHLTFloat32_t ComputeMass(
47 AliHLTFloat32_t massA,
51 AliHLTFloat32_t massB,
57 static bool FitLineToTriggerRecord(const AliHLTMUONTriggerRecordStruct& trigger);
59 static bool FitLineToTriggerRecord(
60 const AliHLTMUONTriggerRecordStruct& trigger,
64 static AliHLTFloat32_t IdealZ1() { return fgIdealZ1; }
65 static void IdealZ1(AliHLTFloat32_t value) { fgIdealZ1 = value; }
66 static AliHLTFloat32_t IdealZ2() { return fgIdealZ2; }
67 static void IdealZ2(AliHLTFloat32_t value) { fgIdealZ2 = value; }
69 static AliHLTFloat32_t IdealX1() { return fgIdealX1; }
70 static AliHLTFloat32_t IdealY1() { return fgIdealY1; }
71 static AliHLTFloat32_t IdealX2() { return fgIdealX2; }
72 static AliHLTFloat32_t IdealY2() { return fgIdealY2; }
74 static bool FitLineToData(
75 const AliHLTFloat32_t* x, const AliHLTFloat32_t* y,
76 const AliHLTFloat32_t* z, AliHLTUInt32_t n
79 static AliHLTFloat32_t Mzx() { return fgMzx; }
80 static AliHLTFloat32_t Mzy() { return fgMzy; }
81 static AliHLTFloat32_t Czx() { return fgCzx; }
82 static AliHLTFloat32_t Czy() { return fgCzy; }
84 static AliHLTFloat32_t ComputeChi2(
85 const AliHLTFloat32_t* x, const AliHLTFloat32_t* y,
86 const AliHLTFloat32_t* z, AliHLTUInt32_t n
89 static AliHLTFloat32_t SigmaX2() { return fgSigmaX2; }
90 static void SigmaX2(AliHLTFloat32_t value) { fgSigmaX2 = (value != 0 ? value : 1.); }
91 static AliHLTFloat32_t SigmaY2() { return fgSigmaY2; }
92 static void SigmaY2(AliHLTFloat32_t value) { fgSigmaY2 = (value != 0 ? value : 1.); }
96 // Prevent destroying or creating of this object.
97 AliHLTMUONCalculations();
98 ~AliHLTMUONCalculations();
100 static AliHLTFloat32_t fgZf; /// The Z coordinate of the middle of the dipole magnetic field.
101 static AliHLTFloat32_t fgQBLScaled; /// The integrated field strength times units of charge (T.m.*c/1e9)
103 static AliHLTMUONParticleSign fgSign; /// The calculated sign.
104 static AliHLTFloat32_t fgPx; /// The calculated X momentum (GeV/c).
105 static AliHLTFloat32_t fgPy; /// The calculated Y momentum (GeV/c).
106 static AliHLTFloat32_t fgPz; /// The calculated Z momentum (GeV/c).
108 static AliHLTFloat32_t fgSigmaX2; /// The sigma squared value for the variance / uncertainty in X coordinates.
109 static AliHLTFloat32_t fgSigmaY2; /// The sigma squared value for the variance / uncertainty in Y coordinates.
111 static AliHLTFloat32_t fgMzx; /// Calculated slope of the line fitted to the ZX plane. (x = fgMzx * z + fgCzx)
112 static AliHLTFloat32_t fgMzy; /// Calculated slope of the line fitted to the ZY plane. (y = fgMzy * z + fgCzy)
113 static AliHLTFloat32_t fgCzx; /// Calculated coefficient of the line fitted to the ZX plane. (x = fgMzx * z + fgCzx)
114 static AliHLTFloat32_t fgCzy; /// Calculated coefficient of the line fitted to the ZY plane. (y = fgMzy * z + fgCzy)
116 static AliHLTFloat32_t fgIdealX1; /// Ideal X coordinate of the point on MT1
117 static AliHLTFloat32_t fgIdealY1; /// Ideal Y coordinate of the point on MT1
118 static AliHLTFloat32_t fgIdealZ1; /// Ideal Z coordinate of the point on MT1
119 static AliHLTFloat32_t fgIdealX2; /// Ideal X coordinate of the point on MT2
120 static AliHLTFloat32_t fgIdealY2; /// Ideal Y coordinate of the point on MT2
121 static AliHLTFloat32_t fgIdealZ2; /// Ideal Z coordinate of the point on MT2
124 #endif // ALIHLTMUONCALCULATIONS_H