Changeset 28 in JEM-EUSO

Timestamp:

Jan 15, 2013, 2:50:08 PM (11 years ago)

Author:

moretto

Message:

Revision 3029 of the official ccin2p3 repository: Bug in treatment of bg3.dat in OpticalAdaptor was fixed. Now optical adaptors with BG3 should work. BG3 data were extrapolated to shorter wavelength.

Location:

esaf_lal/trunk/esaf

Files:

• config/Optics/NOpticalSystem/bg3b.dat (added)
• packages/simulation/detector/optics/src/OpticalAdaptor.cc (modified) (3 diffs)

esaf_lal/trunk/esaf/packages/simulation/detector/optics/src/OpticalAdaptor.cc

@@ -1 @@
-// $Id: OpticalAdaptor.cc 2849 2010-11-11 11:59:41Z biktem $
+// $Id: OpticalAdaptor.cc 3029 2013-01-15 13:08:55Z sakaki $
 // Author: D.Demarco, M.Pallavicini
 
@@ -44 +44 @@
 
     if (!fgBG3) {
-        fgBG3 = new Interpolate("config/Optics/NOpticalSystem/bg3.dat",2);
+      // ATTENSION!! 
+      // The 3rd column in NOpticalSystem/bg3.dat is the imaginary part of
+      //   the refractive index.
+        fgBG3 = new Interpolate("config/Optics/NOpticalSystem/bg3b.dat",2);
         fgBG3->SetXUnit(nm);
     }
@@ -91 +94 @@
 
 //            value = Exp(fFilterThickness*Log(fgBG3->GetValue(p->wl))/(1.*mm));
-            value = Power(fgBG3->GetValue(p->wl),
-                    fFilterThickness/(1.*mm));
+//            value = Power(fgBG3->GetValue(p->wl),
+//                    fFilterThickness/(1.*mm));
+            value=TMath::Exp(-4.0*TMath::Pi()/p->wl*fgBG3->GetValue(p->wl,1)
+                             *fFilterThickness);
 
             return ( EsafRandom::Get()->Rndm() > value );