]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDsimTR.cxx
Simulation of ArCO2 mixture
[u/mrichter/AliRoot.git] / TRD / AliTRDsimTR.cxx
index b940f7b8f285599fff6bcbfd1a34e9ae02c2b817..a1a4a7a899bc4dcdc0c98d2903a18f70bcd120be 100644 (file)
@@ -445,7 +445,7 @@ Int_t AliTRDsimTR::TrPhotons(Float_t p, Float_t mass
     energykeV = fSpectrum->GetBinCenter(iBin);
     energyeV  = energykeV * 1.0e3;
 
-    sigma    = Sigma(energykeV);
+    sigma     = Sigma(energykeV);
 
     csi1      = fFoilOmega / energyeV;
     csi2      = fGapOmega  / energyeV;
@@ -698,44 +698,46 @@ Double_t AliTRDsimTR::GetMuXe(Double_t energyMeV)
                     , 5.00000E+00, 6.00000E+00, 8.00000E+00
                    , 1.00000E+01, 1.50000E+01, 2.00000E+01 };
 
 return Interpolate(energyMeV,en,mu,kN);
+ return Interpolate(energyMeV,en,mu,kN);
 
 }
 
 //_____________________________________________________________________________
-Double_t AliTRDsimTR::GetMuBu(Double_t energyMeV)
+Double_t AliTRDsimTR::GetMuAr(Double_t energyMeV)
 {
   //
-  // Returns the photon absorbtion cross section for isobutane
+  // Returns the photon absorbtion cross section for argon
   //
 
-  const Int_t kN = 36;
-
-  Double_t mu[kN] = { 0.38846E+03, 0.12291E+03, 0.53225E+02
-                    , 0.16091E+02, 0.69114E+01, 0.36541E+01
-                    , 0.22282E+01, 0.11149E+01, 0.72887E+00
-                    , 0.45053E+00, 0.38167E+00, 0.33920E+00
-                    , 0.32155E+00, 0.30949E+00, 0.29960E+00
-                    , 0.28317E+00, 0.26937E+00, 0.24228E+00
-                    , 0.22190E+00, 0.19289E+00, 0.17288E+00
-                    , 0.15789E+00, 0.14602E+00, 0.12829E+00
-                    , 0.11533E+00, 0.10310E+00, 0.93790E-01
-                    , 0.80117E-01, 0.63330E-01, 0.53229E-01
-                    , 0.46390E-01, 0.41425E-01, 0.34668E-01
-                   , 0.30267E-01, 0.23910E-01, 0.20509E-01 };
+  const Int_t kN = 38;
 
-  Double_t en[kN] = { 0.10000E-02, 0.15000E-02, 0.20000E-02
-                    , 0.30000E-02, 0.40000E-02, 0.50000E-02
-                    , 0.60000E-02, 0.80000E-02, 0.10000E-01
-                    , 0.15000E-01, 0.20000E-01, 0.30000E-01
-                    , 0.40000E-01, 0.50000E-01, 0.60000E-01
-                    , 0.80000E-01, 0.10000E+00, 0.15000E+00
-                    , 0.20000E+00, 0.30000E+00, 0.40000E+00
-                    , 0.50000E+00, 0.60000E+00, 0.80000E+00
-                    , 0.10000E+01, 0.12500E+01, 0.15000E+01
-                    , 0.20000E+01, 0.30000E+01, 0.40000E+01
-                    , 0.50000E+01, 0.60000E+01, 0.80000E+01
-                   , 0.10000E+02, 0.15000E+02, 0.20000E+02 };
+  Double_t mu[kN] = { 3.184E+03, 1.105E+03, 5.120E+02
+                    , 1.703E+02, 1.424E+02, 1.275E+03
+                    , 7.572E+02, 4.225E+02, 2.593E+02
+                    , 1.180E+02, 6.316E+01, 1.983E+01
+                    , 8.629E+00, 2.697E+00, 1.228E+00
+                    , 7.012E-01, 4.664E-01, 2.760E-01
+                    , 2.043E-01, 1.427E-01, 1.205E-01
+                    , 9.953E-02, 8.776E-02, 7.958E-02
+                    , 7.335E-02, 6.419E-02, 5.762E-02
+                    , 5.150E-02, 4.695E-02, 4.074E-02
+                    , 3.384E-02, 3.019E-02, 2.802E-02
+                    , 2.667E-02, 2.517E-02, 2.451E-02
+                    , 2.418E-02, 2.453E-02 };
+
+  Double_t en[kN] = { 1.00000E-03, 1.50000E-03, 2.00000E-03  
+                    , 3.00000E-03, 3.20290E-03, 3.20290E-03  
+                    , 4.00000E-03, 5.00000E-03, 6.00000E-03  
+                    , 8.00000E-03, 1.00000E-02, 1.50000E-02  
+                    , 2.00000E-02, 3.00000E-02, 4.00000E-02  
+                    , 5.00000E-02, 6.00000E-02, 8.00000E-02  
+                    , 1.00000E-01, 1.50000E-01, 2.00000E-01  
+                    , 3.00000E-01, 4.00000E-01, 5.00000E-01  
+                    , 6.00000E-01, 8.00000E-01, 1.00000E+00  
+                    , 1.25000E+00, 1.50000E+00, 2.00000E+00  
+                    , 3.00000E+00, 4.00000E+00, 5.00000E+00  
+                    , 6.00000E+00, 8.00000E+00, 1.00000E+01  
+                   , 1.50000E+01, 2.00000E+01 };
 
   return Interpolate(energyMeV,en,mu,kN);