AliRoot Core  v5-06-15 (45dab64)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliMUONConstants.cxx
Go to the documentation of this file.
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 //-----------------------------------------------------------------------------
35 //-----------------------------------------------------------------------------
36 
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 
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 
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 
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 
116 
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 //______________________________________________________________________________
129 {
131 }
132 
133 //______________________________________________________________________________
135 {
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  if (fgDefaultChamberMinZ[0] == 0) // Are the min/max Z arrays initialised?
145  {
146  // The min and max Z arrays need to be calculated.
147  for (Int_t i = 0; i < NCh(); i++)
148  {
149  Float_t a = 0, b = 0;
150  if (4 <= i and i < 10)
151  {
152  Float_t dzAngle = TMath::Tan(TMath::Pi()*St345Inclination()/180.) * Rmax(i/2);
153  // We add 2.5mm since Rmax is an under-estimate.
154  a = DefaultChamberZ(i) + DzSlat() + DzCh() + dzAngle + 0.25;
155  b = DefaultChamberZ(i) - DzSlat() - DzCh() - dzAngle - 0.25;
156  }
157  else
158  {
159  a = DefaultChamberZ(i) + DzSlat();
160  b = DefaultChamberZ(i) - DzSlat();
161  }
162  fgDefaultChamberMinZ[i] = TMath::Min(a, b);
163  fgDefaultChamberMaxZ[i] = TMath::Max(a, b);
164  }
165  }
166 
167  // We can apply a binary search for the chamber since the fgDefaultChamberMinZ and
168  // fgDefaultChamberMaxZ arrays are ordered.
169  Int_t mini = 0, maxi = NCh()-1;
170  while (mini <= maxi)
171  {
172  Int_t iChamber = (maxi + mini) / 2;
173  if (z < fgDefaultChamberMinZ[iChamber])
174  mini = iChamber+1;
175  else if (z > fgDefaultChamberMaxZ[iChamber])
176  maxi = iChamber-1;
177  else
178  // We are between min and max Z of chamber number iChamber so we found our chamber.
179  return iChamber;
180  }
181 
182  if (warn) AliWarningClass(Form("No chamber number found for z = %f",z));
183  return -1;
184 }
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  Float_t q = qtot*1.19*(1.24-timeDif*1E6)*TMath::Exp(-(0.97-timeDif*1E6)*(0.97-timeDif*1E6)/2.42);
191  return q;
192 }
193 
194 //______________________________________________________________________________
196  // Return conversionfactor fc to adc
197  Float_t fc2adc = 1./(DefaultA0()*DefaultCapa()*DefaultADC2MV());
198 
199  return fc2adc;
200 }
static Int_t NofChambers()
Return number of chambers.
static Float_t fgDefaultChamberMinZ[14]
! Minimum z-positions of chambers
static Float_t fgPitch
Wire pitch for St2 & Slats.
static Float_t fgDefaultChamberMaxZ[14]
! Maximum z-positions of chambers
static Float_t fgPitchSt1
Wire pitch for Station 1.
static const Double_t fgkMuonFilterZEnd
z-position of the end of the muon filter
static Float_t fgDefaultRatioTriggerChamber[4]
Ratio between trigger chambers.
static Int_t NTrackingCh()
Return number of tracking chambers.
static Float_t St345Inclination()
Return Inclination with respect the vertical axis of stations 345.
static const Double_t fgkAbsZBeg
z-position of the begining of the absorber
static Int_t fgMaxZoom
Maximum Zoom for event display.
static Float_t fgSqrtKx3
SqrtKx3 for St2 & Slat.
static Int_t ChamberNumber(Float_t z, bool warn=false)
Return chamber number according z position of hit.
static Double_t fgChamberThicknessInX0[10]
default chamber thickness in X0 for reconstruction
static Double_t fgYokeL
Yoke lenght.
static const Double_t fgkMuonFilterZBeg
z-position of the begining of the muon filter
static Float_t DefaultA0()
static Float_t ReducedQTot(Float_t qtot, Float_t timeDif)
Return a reduced total charge.
static Float_t fgSqrtKx3St1
SqrtKx3 for Station 1.
ClassImp(TPCGenInfo)
Definition: AliTPCCmpNG.C:254
static Float_t fgSqrtKy3
SqrtKy3 for St2 & Slat.
static Float_t Rmax(Int_t i)
Return chamber i outer radius.
static Float_t DzCh()
Return half-distance between two half-chambers.
static const Float_t fgkDefaultADC2MV
Default ADC to mv conversion factor.
static Float_t fgAverageChamberT[14]
average arrival time to chamber
static const Double_t fgkMostProbBendingMomentum
most probable value (GeV/c) of muon momentum in bending plane (used when B = 0) needed to get some "r...
static const Double_t fgkAbsZEnd
z-position of the end of the absorber
MUON global constants.
static Float_t DefaultCapa()
static const Float_t fgkDefaultCapa
Default Capa.
static Float_t DefaultADC2MV()
static Double_t fgCoilL
Coil lenght.
static Double_t fgYokeZ
Yoke z-position.
static Int_t NCh()
Return number of chambers.
static const Float_t fgkDefaultA0
Default gain A0.
static Double_t fgCoilZ
Coil z-position.
static const Double_t fgkTriggerBendingReso
default trigger chamber resolution in the bending direction, DUE TO ALIGNMENT ONLY ...
static Float_t fgChargeCorrel
Charge correlation for St2 & Slats.
static Float_t fgDmin[7]
! Inner diameter
static Float_t * DefaultChamberZ()
Return pointer to array of positions.
static Float_t fgSqrtKy3St1
SqrtKy3 for Station 1.
static Float_t fgDzCh
Half-distance between two half-chambers.
static Float_t FC2ADC()
static const Double_t fgkTriggerNonBendingReso
default trigger chamber resolution in the non-bending direction, DUE TO ALIGNMENT ONLY ...
static Float_t fgDmax[7]
! Outer diameter
static Float_t fgDzSlat
Half-distance between two slat on the same chamber.
static Float_t DzSlat()
Return half-distance between two slats.
static Float_t fgChargeCorrelSt1
Charge correlation for Station 1.
static Float_t fgSt345inclination
! Inclination with respect the vertical axis of stations 345
static Int_t NofTrackingChambers()
Return number of tracking chambers.
static const Double_t fgkMuonFilterX0
x0 of the muon filter