lcsim/test/org/lcsim/recon/tracking/trfeloss

diff -N DeDxFixed_t.java
--- /dev/null	1 Jan 1970 00:00:00 -0000
+++ DeDxFixed_t.java	25 Jul 2007 17:23:39 -0000	1.1
@@ -0,0 +1,56 @@

+/*
+ * DeDxFixed_t.java
+ *
+ * Created on July 24, 2007, 10:08 PM
+ *
+ * $Id: DeDxFixed_t.java,v 1.1 2007/07/25 17:23:39 ngraf Exp $
+ */
+
+package org.lcsim.recon.tracking.trfeloss;
+
+import junit.framework.TestCase;
+import org.lcsim.recon.tracking.trfutil.Assert;
+
+/**
+ *
+ * @author Norman Graf
+ */
+public class DeDxFixed_t extends TestCase
+{
+    private boolean debug;
+    /** Creates a new instance of DeDxFixed_t */
+    public void testDeDxFixed()
+    {
+        String component = "DeDxFixed";
+        String ok_prefix = component + " (I): ";
+        String error_prefix = component + " test (E): ";
+        
+        if(debug) System.out.println( ok_prefix
+                + "---------- Testing component " + component
+                + ". ----------");
+        
+        //********************************************************************
+        
+        double pionMass = 0.13957;
+        double density = 3.0;
+        DeDx dedx = new DeDxFixed(density);
+        double thickness = 10.0;
+        double energy = 10.0;
+        double testenergy = energy;
+        
+        double de_dx = dedx.dEdX(energy);
+        testenergy -= de_dx*thickness;
+        double finalenergy = dedx.loseEnergy(energy, thickness);
+        if(debug) System.out.println("energy= "+energy+", testenergy= "+testenergy+", finalenergy= "+finalenergy);
+        Assert.assertTrue(Math.abs(finalenergy-testenergy)<0.0001);
+        
+        double gainenergy = dedx.loseEnergy(energy, -thickness);
+        if(debug) System.out.println("energy= "+energy+", gainenergy= "+gainenergy);
+        Assert.assertTrue( gainenergy > energy);
+        
+        
+        if(debug) System.out.println( ok_prefix
+                + "------------- All tests passed. -------------");
+    }
+    
+}

lcsim/test/org/lcsim/recon/tracking/trfeloss

diff -N DeDx_t.java
--- /dev/null	1 Jan 1970 00:00:00 -0000
+++ DeDx_t.java	25 Jul 2007 17:23:39 -0000	1.1
@@ -0,0 +1,45 @@

+/*
+ * DeDx_t.java
+ *
+ * Created on July 24, 2007, 10:07 PM
+ *
+ * $Id: DeDx_t.java,v 1.1 2007/07/25 17:23:39 ngraf Exp $
+ */
+
+package org.lcsim.recon.tracking.trfeloss;
+
+import junit.framework.TestCase;
+import org.lcsim.recon.tracking.trfutil.Assert;
+
+/**
+ *
+ * @author Norman Graf
+ */
+public class DeDx_t extends TestCase
+{
+    private boolean debug;
+    /** Creates a new instance of DeDx_t */
+    public void testDeDx()
+    {
+        String component = "DeDx";
+        String ok_prefix = component + " (I): ";
+        String error_prefix = component + " test (E): ";
+        
+        if(debug) System.out.println( ok_prefix
+                + "-------- Testing component " + component
+                + ". --------" );
+        
+        DeDxTest test = new DeDxTest();
+        if(debug) System.out.println(test);
+        double energy = 1.;
+        double x = 10.;
+        Assert.assertTrue(test.dEdX(energy)==0.1*energy);
+        Assert.assertTrue(test.sigmaEnergy(energy,x)==0.01*energy*x);
+        Assert.assertTrue(test.loseEnergy(energy,x)==-energy*x);
+        
+        if(debug) System.out.println( ok_prefix
+                + "------------- All tests passed. -------------" );
+        
+    }
+    
+}

lcsim/test/org/lcsim/recon/tracking/trfeloss

diff -N DeDxBethe_t.java
--- /dev/null	1 Jan 1970 00:00:00 -0000
+++ DeDxBethe_t.java	25 Jul 2007 17:23:39 -0000	1.1
@@ -0,0 +1,58 @@

+/*
+ * DeDxBethe_t.java
+ *
+ * Created on July 24, 2007, 10:10 PM
+ *
+ * $Id: DeDxBethe_t.java,v 1.1 2007/07/25 17:23:39 ngraf Exp $
+ */
+
+package org.lcsim.recon.tracking.trfeloss;
+
+import junit.framework.TestCase;
+import org.lcsim.recon.tracking.trfutil.Assert;
+
+/**
+ *
+ * @author Norman Graf
+ */
+public class DeDxBethe_t extends TestCase
+{
+    private boolean debug;
+    /** Creates a new instance of DeDxBethe_t */
+    public void testDeDxBethe()
+    {
+        String component = "DeDxBethe";
+        String ok_prefix = component + " (I): ";
+        String error_prefix = component + " test (E): ";
+        
+        if(debug) System.out.println( ok_prefix
+                + "---------- Testing component " + component
+                + ". ----------");
+        
+        //********************************************************************
+        
+        double pionMass = 0.13957;
+        double density = 3.0;
+        DeDx dedx = new DeDxBethe(density);
+        double thickness = 10.0;
+        double energy = 10.0;
+        double testenergy = energy;
+        // not much of a test here...
+        double de_dx = dedx.dEdX(energy);
+        testenergy -= de_dx*thickness;
+        double finalenergy = dedx.loseEnergy(energy, thickness);
+        if(debug) System.out.println("energy= "+energy+", testenergy= "+testenergy+", finalenergy= "+finalenergy);
+        
+        // Check that we gain energy if we go backwards
+        Assert.assertTrue(Math.abs(testenergy-finalenergy)<0.0001);
+        
+        double gainenergy = dedx.loseEnergy(energy, -thickness);
+        Assert.assertTrue( gainenergy > energy);
+        
+        if(debug) System.out.println( ok_prefix
+                + "------------- All tests passed. -------------");
+        
+        
+    }
+    
+}