X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FAliMUONConstants.cxx;h=01fec147f3aec2896cf972b192304276d17393d0;hb=cf8425acaddcd91296e40f337aeced83bc0d2f70;hp=7e2d3b9606a88687dd6ce715684f16c6ef5ecf4e;hpb=c0f6f01db6537dab6015782c47870ee286a0fda3;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/AliMUONConstants.cxx b/MUON/AliMUONConstants.cxx index 7e2d3b9606a..01fec147f3a 100644 --- a/MUON/AliMUONConstants.cxx +++ b/MUON/AliMUONConstants.cxx @@ -21,8 +21,10 @@ #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) @@ -30,16 +32,13 @@ /// 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; @@ -50,6 +49,13 @@ Float_t AliMUONConstants::fgDefaultChamberZ[14] = -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}; @@ -75,7 +81,7 @@ Float_t AliMUONConstants::fgPitchSt1 = 0.21; // 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; @@ -88,35 +94,81 @@ Double_t AliMUONConstants::fgCoilL = 502.1; 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; + +//______________________________________________________________________________ +Int_t AliMUONConstants::NCh() +{ + return AliMpConstants::NofChambers(); +} + +//______________________________________________________________________________ +Int_t AliMUONConstants::NTrackingCh() +{ + return AliMpConstants::NofTrackingChambers(); +} //______________________________________________________________________________ -Int_t AliMUONConstants::ChamberNumber(Float_t z) +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; }