GeomConverter/src/org/lcsim/material

diff -u -r1.18 -r1.19
--- Material.java	23 Jul 2005 00:44:08 -0000	1.18
+++ Material.java	29 Jul 2005 07:43:07 -0000	1.19
@@ -79,6 +79,9 @@

     /* radiation length in cm, taking into account the density */
     double _radiationLengthWithDensity;
     

+    /* Moliere radius in cm */
+    private double _moliereRadius;
+    

     /* nuclear interaction length in g/cm2 */
     double _nuclearInteractionLength;
     

@@ -258,6 +261,16 @@

         _radiationLength = (rlinv <= 0.0 ? MAX_RADIATION_LENGTH : 1.0/rlinv);
         _radiationLengthWithDensity = _radiationLength / _density;
         

+//
+// Set the critical energy
+//
+        double criticalEnergy  = 2.66 * pow(_radiationLength * _Zeff/_Aeff, 1.1);
+//
+// Set Moliere radius
+//
+        _moliereRadius   = _radiationLengthWithDensity * 21.2052 / criticalEnergy;
+//        _moliereRadius = 0.01 * rMoliere / _density;
+        

         return _radiationLength;
     }
     

@@ -276,7 +289,14 @@

         return _radiationLength;
     }
     

+
+    public double getMoliereRadius()
+    {
+        return _moliereRadius;
+    }
+

     /** @return lambda */

+

     public double getNuclearInteractionLength()
     {
         return _nuclearInteractionLength;