lcsim/src/org/lcsim/recon/tracking/seedtracker/diagnostic

diff -u -r1.5 -r1.6
--- SeedTrackerDiagnostics.java	19 Aug 2009 22:11:33 -0000	1.5
+++ SeedTrackerDiagnostics.java	1 Sep 2009 22:55:34 -0000	1.6
@@ -5,6 +5,8 @@

 
 package org.lcsim.recon.tracking.seedtracker.diagnostic;
 

+import hep.physics.vec.Hep3Vector;
+import org.lcsim.event.MCParticle;

 import org.lcsim.fit.helicaltrack.HelicalTrackFitter;
 import org.lcsim.fit.helicaltrack.HelicalTrackHit;
 import org.lcsim.recon.tracking.seedtracker.SeedCandidate;

@@ -25,12 +27,26 @@

 
     @Override
     public void fireCheckHitPairFailed(HelicalTrackHit hit1, HelicalTrackHit hit2, SeedCandidate seed) {

-        if (seed.isTrueSeed()) System.out.println("True seed failed hit check");

+        if (seed.isTrueSeed()) {
+            System.out.println("True seed failed hit pair check");
+            for (MCParticle mcp : seed.getMCParticles()) {
+                Hep3Vector p = mcp.getMomentum();
+                double px = p.x();
+                double py = p.y();
+                double pz = p.z();
+                double pt = Math.sqrt(px*px+py*py);
+                double eta = -Math.log(Math.tan(Math.acos(pz / p.magnitude()) / 2.));
+                System.out.println("pT: "+pt+" eta: "+eta);
+                System.out.println("Number of hits in seed: "+seed.getHits().size());
+                System.out.println("Hit 1 position: "+hit1.getCorrectedPosition().toString());
+                System.out.println("Hit 2 position: "+hit2.getCorrectedPosition().toString());
+            }
+        }

     }
 
     @Override
     public void fireCheckHitFailed(HelicalTrackHit hit, SeedCandidate seed) {

-        if (seed.isTrueSeed()) System.out.println("True seed failed hit check");

+        if (seed.isTrueSeed()) System.out.println("True seed failed hit and seed check");

     }
 
     @Override