hps-java/src/main/java/org/lcsim/hps/users/phansson

diff -u -r1.16 -r1.17
--- TwoTrackAnlysis.java	20 Apr 2013 04:00:25 -0000	1.16
+++ TwoTrackAnlysis.java	23 Apr 2013 00:07:37 -0000	1.17
@@ -140,9 +140,8 @@

         List<Track> tracklist = null;
         if(event.hasCollection(Track.class,trackCollectionName)) {        
             tracklist = event.get(Track.class, trackCollectionName);

-             if(_debug) {
-                System.out.println(this.getClass().getSimpleName() + ": Number of Tracks = " + tracklist.size() + " in event " + event.getEventNumber());
-             }

+            //if(_debug)
+            System.out.println(this.getClass().getSimpleName() + ": Number of Tracks = " + tracklist.size() + " in event " + event.getEventNumber());

         } else {
              if(_debug) {
                 System.out.println(this.getClass().getSimpleName() + ": No track collection in event " + event.getEventNumber());

@@ -161,9 +160,18 @@

         ArrayList<CmpTrack> tracks = new ArrayList<CmpTrack>();
         for(int i=0;i<tracklist.size();++i) {
             Track trk = tracklist.get(i);

-            if(TrackUtils.passTrackSelection(trk, tracklist, EventQuality.Quality.MEDIUM)) {
-            if(_debug) 
-                System.out.printf("%s: trk momentum (%.3f,%.3f,%.3f) chi2=%.3f\n",this.getClass().getSimpleName(),trk.getTrackStates().get(0).getMomentum()[0],trk.getTrackStates().get(0).getMomentum()[1],trk.getTrackStates().get(0).getMomentum()[2],trk.getChi2());

+            if(TrackUtils.isGoodTrack(trk, tracklist, EventQuality.Quality.MEDIUM)) {
+                if(_debug) {
+                    int cuts = TrackUtils.passTrackSelections(trk, tracklist, EventQuality.Quality.MEDIUM);
+                    System.out.printf("%s: track cuts: \n%s\n",this.getClass().getSimpleName(),EventQuality.instance().print(cuts));
+                    System.out.printf("%s: trk momentum (%.3f,%.3f,%.3f) chi2=%.3f\n",this.getClass().getSimpleName(),trk.getTrackStates().get(0).getMomentum()[0],trk.getTrackStates().get(0).getMomentum()[1],trk.getTrackStates().get(0).getMomentum()[2],trk.getChi2());
+                }
+                if(trk.getChi2()>10.) {
+                    System.out.printf("%s: trk momentum (%.3f,%.3f,%.3f) chi2=%.3f\n",this.getClass().getSimpleName(),trk.getTrackStates().get(0).getMomentum()[0],trk.getTrackStates().get(0).getMomentum()[1],trk.getTrackStates().get(0).getMomentum()[2],trk.getChi2());
+                    int cuts = TrackUtils.passTrackSelections(trk, tracklist, EventQuality.Quality.MEDIUM);
+                    System.out.printf("%s: track cuts: \n%s\n",this.getClass().getSimpleName(),EventQuality.instance().print(cuts));
+                    System.exit(0);
+                }

                 tracks.add(new CmpTrack(trk));
             } else {
                if(_debug) System.out.println(this.getClass().getSimpleName() + ": trk failed track selections " + event.getEventNumber() + "\n" + trk.toString());