Go from pointer to ifstream to ifstream.

[u/mrichter/AliRoot.git] / MUON / AliMUONConstants.cxx

diff --git a/MUON/AliMUONConstants.cxx b/MUON/AliMUONConstants.cxx
index c3997abd0a296f432c93af0d1c447fe014fe4c44..6ba88928659a256e71b7920a1c218e6c615b18ff 100644 (file)
--- a/MUON/AliMUONConstants.cxx
+++ b/MUON/AliMUONConstants.cxx
@@ -99,8 +99,10 @@ const Double_t AliMUONConstants::fgkAbsZBeg = -90.;
 const Double_t AliMUONConstants::fgkAbsZEnd = -505.;
     
 // Default trigger chamber resolution (cm)
-const Double_t AliMUONConstants::fgkTriggerNonBendingReso = 1.;
-const Double_t AliMUONConstants::fgkTriggerBendingReso = 1.;
+// 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.;
@@ -108,10 +110,20 @@ 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};
+
+// is equivalent to gain = 4 mV/fC
+const Float_t AliMUONConstants::fgkDefaultA0 = 1.25; 
+// 1 ADC channel = 0.61 mV
+const Float_t AliMUONConstants::fgkDefaultADC2MV = 0.61; 
+const Float_t AliMUONConstants::fgkDefaultCapa = 0.2; 
+
 //______________________________________________________________________________
 Int_t AliMUONConstants::NCh()
 {
@@ -170,3 +182,19 @@ Int_t AliMUONConstants::ChamberNumber(Float_t z, bool warn)
   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;
+}
+
+//______________________________________________________________________________
+Float_t AliMUONConstants::FC2ADC() {
+  // Return conversionfactor fc to adc
+  Float_t fc2adc = 1./(DefaultA0()*DefaultCapa()*DefaultADC2MV());
+  
+  return fc2adc;
+}