#include "TMath.h"
#include "TClass.h"
+#include "AliMpConstants.h"
-///
+//-----------------------------------------------------------------------------
+/// \class AliMUONConstants
/// This class holds various constants to be used in many places,
/// such as the number of tracking and trigger chambers,
/// some geometrical constants (to build the initial geometry for simulation)
/// Those constants should as much as possible replace hard-coded values
/// which are to be considered strictly illegal in the MUON code (or any code,
/// by the way).
-///
+//-----------------------------------------------------------------------------
/// \cond CLASSIMP
ClassImp(AliMUONConstants)
/// \endcond
-Int_t AliMUONConstants::fgNCh = 14;
-Int_t AliMUONConstants::fgNTrackingCh = 10;
Int_t AliMUONConstants::fgNTrackingSt = 5;
-Int_t AliMUONConstants::fgNTriggerCh = 4;
Int_t AliMUONConstants::fgNTriggerCircuit = 234;
Int_t AliMUONConstants::fgNDetElem = 228;
Int_t AliMUONConstants::fgNGeomModules = 20;
-967.5, -998.5, -1276.5, -1307.5, -1406.6, -1437.6,// updated 08/05, EDMS id 335328 (A. Tournaire)
-1603.5, -1620.5, -1703.5, -1720.5}; // M1 & M2
+
+// These are used by AliMUONConstants::ChamberNumber and must be calculated once
+// by that method from fgDzCh, fgDzSlat, fgDefaultChamberZ and fgSt345inclination,
+// so for now we set everything to zero.
+Float_t AliMUONConstants::fgDefaultChamberMinZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+Float_t AliMUONConstants::fgDefaultChamberMaxZ[14] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+
Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] =
{1., 1.01060, 1.06236, 1.07296};
// EDMS Id 335328 for "search in EDMS
// These are the diameter (Dmin == innner and Dmax - outner) values of the active surface
// In the case of Dmax, the value corresponds to the maximum diameter of the active surface with 2pi coverture in phi
-Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 63.0, 79.0, 79.0, 99.0, 100.0}; // cm
+Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 63.0, 79.0, 79.0, 98.8, 100.0}; // cm
Float_t AliMUONConstants::fgDmax[7] = { 176.6, 229.0, 308.84, 418.2, 522.0, 850.0, 900.0}; // cm
Int_t AliMUONConstants::fgMaxZoom = 20;
Double_t AliMUONConstants::fgYokeZ = -986.6;
Double_t AliMUONConstants::fgYokeL = 309.4;
+// Defaults parameters for absorber (27/06/07)
+const Double_t AliMUONConstants::fgkAbsZBeg = -90.;
+const Double_t AliMUONConstants::fgkAbsZEnd = -505.;
+
+// Default trigger chamber resolution (cm)
+// Warning: the resolution refers only to ALIGNMENT
+// For the total resolution the strip width should be taken into account!
+const Double_t AliMUONConstants::fgkTriggerNonBendingReso = 0.2;
+const Double_t AliMUONConstants::fgkTriggerBendingReso = 0.2;
+
+// Defaults parameters for muon filter (19/11/07)
+const Double_t AliMUONConstants::fgkMuonFilterZBeg = -1471.;
+const Double_t AliMUONConstants::fgkMuonFilterZEnd = -1471.-120.;
+const Double_t AliMUONConstants::fgkMuonFilterX0 = 1.76;
+
// Defaults parameters for track reconstruction
-Double_t AliMUONConstants::fgChamberThicknessInX0 = 0.03;
+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};
+
+const Double_t AliMUONConstants::fgkMostProbBendingMomentum = 2.0;
+
+Float_t AliMUONConstants::fgAverageChamberT[14]=
+ {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,
+ 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};
+
+//______________________________________________________________________________
+Int_t AliMUONConstants::NCh()
+{
+ return AliMpConstants::NofChambers();
+}
//______________________________________________________________________________
-Int_t AliMUONConstants::ChamberNumber(Float_t z)
+Int_t AliMUONConstants::NTrackingCh()
+{
+ return AliMpConstants::NofTrackingChambers();
+}
+
+//______________________________________________________________________________
+Int_t AliMUONConstants::ChamberNumber(Float_t z, bool warn)
{
// return chamber number according z position of hit. Should be taken from geometry ?
-
- Float_t dMaxChamber = DzSlat() + DzCh() + 0.25; // cm st 3 &4 & 5
- dMaxChamber += 3.00; // factor for inclination of chamber
- // dMaxChamber += Rmax(4) * TMath::Sin(fgSt345inclination*TMath::Pi()/360.);
- // factor for inclination of chamber
- if ( z > (DefaultChamberZ(4)+50.)) dMaxChamber = 7.; // cm stations 1 & 2
- Int_t iChamber;
-
- for (iChamber = 0; iChamber < 10; iChamber++) {
- if (TMath::Abs(z-DefaultChamberZ(iChamber)) < dMaxChamber) {
- return iChamber;
+
+ if (fgDefaultChamberMinZ[0] == 0) // Are the min/max Z arrays initialised?
+ {
+ // The min and max Z arrays need to be calculated.
+ for (Int_t i = 0; i < NCh(); i++)
+ {
+ Float_t a = 0, b = 0;
+ if (4 <= i and i < 10)
+ {
+ Float_t dzAngle = TMath::Tan(TMath::Pi()*St345Inclination()/180.) * Rmax(i/2);
+ // We add 2.5mm since Rmax is an under-estimate.
+ a = DefaultChamberZ(i) + DzSlat() + DzCh() + dzAngle + 0.25;
+ b = DefaultChamberZ(i) - DzSlat() - DzCh() - dzAngle - 0.25;
+ }
+ else
+ {
+ a = DefaultChamberZ(i) + DzSlat();
+ b = DefaultChamberZ(i) - DzSlat();
+ }
+ fgDefaultChamberMinZ[i] = TMath::Min(a, b);
+ fgDefaultChamberMaxZ[i] = TMath::Max(a, b);
}
}
-
- if ( z > DefaultChamberZ(NTrackingCh()-1) ) {
- AliWarningClass(Form("No chamber number found for z = %f",z));
- // for (iChamber = 0; iChamber < 10; iChamber++) {
- // cout << iChamber << " zpos: " << DefaultChamberZ(iChamber)
- // << " from " << DefaultChamberZ(iChamber) + dMaxChamber
- // << " to " << DefaultChamberZ(iChamber) - dMaxChamber
- // << " delta " << dMaxChamber << endl;
- //}
+
+ // We can apply a binary search for the chamber since the fgDefaultChamberMinZ and
+ // fgDefaultChamberMaxZ arrays are ordered.
+ Int_t mini = 0, maxi = NCh()-1;
+ while (mini <= maxi)
+ {
+ Int_t iChamber = (maxi + mini) / 2;
+ if (z < fgDefaultChamberMinZ[iChamber])
+ mini = iChamber+1;
+ else if (z > fgDefaultChamberMaxZ[iChamber])
+ maxi = iChamber-1;
+ else
+ // We are between min and max Z of chamber number iChamber so we found our chamber.
+ return iChamber;
}
+
+ if (warn) AliWarningClass(Form("No chamber number found for z = %f",z));
return -1;
}
+
+//______________________________________________________________________________
+Float_t AliMUONConstants::ReducedQTot(Float_t qtot, Float_t timeDif)
+{
+ // return a reduced charge if the hit belongs to a track from a pileup event
+ Float_t q = qtot*1.19*(1.24-timeDif*1E6)*TMath::Exp(-(0.97-timeDif*1E6)*(0.97-timeDif*1E6)/2.42);
+ return q;
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff