#include "AliMUONConstants.h"
#include "AliLog.h"
#include "TMath.h"
#include "TClass.h"
#include "AliMpConstants.h"
ClassImp(AliMUONConstants)
Int_t AliMUONConstants::fgNTrackingSt = 5;
Int_t AliMUONConstants::fgNTriggerCircuit = 234;
Int_t AliMUONConstants::fgNDetElem = 228;
Int_t AliMUONConstants::fgNGeomModules = 20;
Float_t AliMUONConstants::fgkTriggerTofLimit = 75E-9;
Float_t AliMUONConstants::fgDefaultChamberZ[14] =
{-526.16, -545.24, -676.4, -695.4,
-967.5, -998.5, -1276.5, -1307.5, -1406.6, -1437.6,
-1603.5, -1620.5, -1703.5, -1720.5};
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};
Float_t AliMUONConstants::fgSt345inclination = 0.794;
Float_t AliMUONConstants::fgDzCh = 15.5/2.;
Float_t AliMUONConstants::fgDzSlat = 8.5/2.;
Float_t AliMUONConstants::fgSqrtKx3 = 0.7131;
Float_t AliMUONConstants::fgSqrtKy3 = 0.7642;
Float_t AliMUONConstants::fgSqrtKx3St1 = 0.7000;
Float_t AliMUONConstants::fgSqrtKy3St1 = 0.7550;
Float_t AliMUONConstants::fgChargeCorrel = 0.11;
Float_t AliMUONConstants::fgChargeCorrelSt1 = 1.0;
Float_t AliMUONConstants::fgPitch = 0.25;
Float_t AliMUONConstants::fgPitchSt1 = 0.21;
Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 63.0, 79.0, 79.0, 98.8, 100.0};
Float_t AliMUONConstants::fgDmax[7] = { 176.6, 229.0, 308.84, 418.2, 522.0, 850.0, 900.0};
Int_t AliMUONConstants::fgMaxZoom = 20;
Double_t AliMUONConstants::fgCoilZ = -994.05;
Double_t AliMUONConstants::fgCoilL = 502.1;
Double_t AliMUONConstants::fgYokeZ = -986.6;
Double_t AliMUONConstants::fgYokeL = 309.4;
const Double_t AliMUONConstants::fgkAbsZBeg = -90.;
const Double_t AliMUONConstants::fgkAbsZEnd = -505.;
const Double_t AliMUONConstants::fgkTriggerNonBendingReso = 0.2;
const Double_t AliMUONConstants::fgkTriggerBendingReso = 0.2;
const Double_t AliMUONConstants::fgkMuonFilterZBeg = -1471.;
const Double_t AliMUONConstants::fgkMuonFilterZEnd = -1471.-120.;
const Double_t AliMUONConstants::fgkMuonFilterX0 = 1.76;
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};
const Float_t AliMUONConstants::fgkDefaultA0 = 1.25;
const Float_t AliMUONConstants::fgkDefaultADC2MV = 0.61;
const Float_t AliMUONConstants::fgkDefaultCapa = 0.2;
Int_t AliMUONConstants::NCh()
{
return AliMpConstants::NofChambers();
}
Int_t AliMUONConstants::NTrackingCh()
{
return AliMpConstants::NofTrackingChambers();
}
Int_t AliMUONConstants::ChamberNumber(Float_t z, bool warn)
{
if (fgDefaultChamberMinZ[0] == 0)
{
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);
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);
}
}
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
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)
{
Float_t q = qtot*1.19*(1.24-timeDif*1E6)*TMath::Exp(-(0.97-timeDif*1E6)*(0.97-timeDif*1E6)/2.42);
return q;
}
Float_t AliMUONConstants::FC2ADC() {
Float_t fc2adc = 1./(DefaultA0()*DefaultCapa()*DefaultADC2MV());
return fc2adc;
}