lcsim/src/org/lcsim/contrib/Partridge/TrackingTest

diff -u -r1.1 -r1.2
--- ResolutionAnalysis.java	8 Nov 2008 02:33:40 -0000	1.1
+++ ResolutionAnalysis.java	10 Nov 2008 19:54:42 -0000	1.2
@@ -49,6 +49,7 @@

                 double ptrk = (new BasicHep3Vector(track.getMomentum())).magnitude();
                 double pmc = mcp.getMomentum().magnitude();
                 double pdif = ptrk - pmc;

+                double theta = Math.acos(mcp.getMomentum().z()/pmc);

                 if (mcp.getCharge() > 0.) {
                     aida.cloud1D("Momentum error for positive tracks").fill(pdif);
                 } else {

@@ -58,6 +59,7 @@

                 aida.cloud1D("Track momentum for all tracks").fill(ptrk);
                 aida.cloud1D("MC momentum for all tracks").fill(pmc);
                 aida.cloud1D("Momentum resolution for all tracks").fill(pdif/pmc);

+                aida.cloud1D("Polar angle").fill(theta);

             }
             
         }