fixed bug in AliSTARTTrigger

[u/mrichter/AliRoot.git] / TFluka / TFlukaCerenkov.cxx

diff --git a/TFluka/TFlukaCerenkov.cxx b/TFluka/TFlukaCerenkov.cxx
index e662a001aafcc04e170888be504024657a9bdcd1..8c17ce3ec2cdb43d908f3a0b75ebd03bea84c737 100644 (file)
--- a/TFluka/TFlukaCerenkov.cxx
+++ b/TFluka/TFlukaCerenkov.cxx
@@ -1,8 +1,37 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/* $Id$*/
+
+//
+// Stores user defined cerenkov properties of media like 
+// absorption coefficient, refraction index and quantum efficiency. 
+// The properties are stored in arrays. The array index corresponds to discrete 
+// optical photon energies defined in fEnergy.
+// Access to the properties is provided by interpolation.
+// 
+// Author:
+// A. Morsch 
+// andreas.morsch@cern.ch
+//
+
 #include "TFlukaCerenkov.h"
 
 Double_t TFlukaCerenkov::fgGlobalMaximumEfficiency = 0.;
    
-ClassImp(TFlukaCerenkov);
+ClassImp(TFlukaCerenkov)
 
 
 TFlukaCerenkov::TFlukaCerenkov()
@@ -14,6 +43,7 @@ TFlukaCerenkov::TFlukaCerenkov()
       fAbsorptionCoefficient(0),
       fQuantumEfficiency(0),
       fRefractionIndex(0),
+      fReflectivity(0),
       fMaximumEfficiency(0.)
 {
 // Default constructor
@@ -29,6 +59,7 @@ TFlukaCerenkov::TFlukaCerenkov(Int_t npckov, Float_t *ppckov, Float_t *absco, Fl
     fAbsorptionCoefficient = new Float_t[fSamples];
     fRefractionIndex       = new Float_t[fSamples];
     fQuantumEfficiency     = new Float_t[fSamples];
+    fReflectivity          = new Float_t[fSamples];    
     
     for (Int_t i = 0; i < fSamples; i++) {
        fEnergy[i]             = ppckov[i];
@@ -46,14 +77,55 @@ TFlukaCerenkov::TFlukaCerenkov(Int_t npckov, Float_t *ppckov, Float_t *absco, Fl
        //
        // Flag is sensitive if quantum efficiency 0 < eff < 1 for at least one value.
        if (effic[i] < 1. && effic[i] > 0.) fIsSensitive = 1;
+       // G3 way to define metal
+       if (rindex[0] == 0.) {
+           fIsMetal = kTRUE;
+           fReflectivity[i] = absco[i];
+       }
     }
     // Find global  maximum quantum efficiency
     if (fMaximumEfficiency > GetGlobalMaximumEfficiency()) {
        SetGlobalMaximumEfficiency(fMaximumEfficiency);
     }
-    printf("Maximum eff. %f\n",  GetGlobalMaximumEfficiency());
-       
+}
+
+TFlukaCerenkov::TFlukaCerenkov(Int_t npckov, Float_t *ppckov, Float_t *absco, Float_t *effic, Float_t *rindex, Float_t *refl)
+{
+    // Constructor including reflectivity
+    fSamples = npckov;
+    fEnergy                = new Float_t[fSamples];
+    fWaveLength            = new Float_t[fSamples];
+    fAbsorptionCoefficient = new Float_t[fSamples];
+    fRefractionIndex       = new Float_t[fSamples];
+    fQuantumEfficiency     = new Float_t[fSamples];
+    
     
+    for (Int_t i = 0; i < fSamples; i++) {
+       fEnergy[i]             = ppckov[i];
+       fWaveLength[i]         = khc / ppckov[i];
+       if (absco[i] > 0.) {
+           fAbsorptionCoefficient[i]   = 1./ absco[i];
+       } else {
+           fAbsorptionCoefficient[i]   = 1.e15;
+       }
+       fRefractionIndex[i]    = rindex[i];
+       fQuantumEfficiency[i]  = effic[i];
+       //
+       // Find local maximum quantum efficiency
+       if (effic[i] > fMaximumEfficiency) fMaximumEfficiency = effic[i];
+       //
+       // Flag is sensitive if quantum efficiency 0 < eff < 1 for at least one value.
+       if (effic[i] < 1. && effic[i] > 0.) fIsSensitive = 1;
+       //
+
+    }
+    // Find global  maximum quantum efficiency
+    if (fMaximumEfficiency > GetGlobalMaximumEfficiency()) {
+       SetGlobalMaximumEfficiency(fMaximumEfficiency);
+    }
+    fReflectivity     = new Float_t[fSamples];
+    for (Int_t i = 0; i < fSamples; i++) fReflectivity[i] = refl[i];
+    fIsMetal = kTRUE;
 }
 
 Float_t TFlukaCerenkov::GetAbsorptionCoefficient(Float_t energy)
@@ -74,6 +146,15 @@ Float_t TFlukaCerenkov::GetRefractionIndex(Float_t energy)
     
 }
 
+Float_t TFlukaCerenkov::GetReflectivity(Float_t energy)
+{
+//
+// Get RefractionIndex for given energy 
+//
+    return Interpolate(energy, fEnergy, fReflectivity);
+    
+}
+
 Float_t TFlukaCerenkov::GetQuantumEfficiency(Float_t energy)
 {
 //
@@ -103,6 +184,15 @@ Float_t TFlukaCerenkov::GetRefractionIndexByWaveLength(Float_t wavelength)
     return Interpolate(energy, fEnergy, fRefractionIndex);
 }
 
+Float_t TFlukaCerenkov::GetReflectivityByWaveLength(Float_t wavelength)
+{
+//
+// Get RefractionIndex for given wavelenth 
+//
+    Float_t energy = khc / wavelength;    
+    return Interpolate(energy, fEnergy, fReflectivity);
+}
+
 Float_t TFlukaCerenkov::GetQuantumEfficiencyByWaveLength(Float_t wavelength)
 {
 //