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

diff -u -r1.11 -r1.12
--- TwoTrackAnlysis.java	14 Mar 2013 03:05:13 -0000	1.11
+++ TwoTrackAnlysis.java	15 Mar 2013 02:52:35 -0000	1.12
@@ -28,6 +28,7 @@

 import org.lcsim.hps.analysis.ecal.HPSMCParticlePlotsDriver;
 import org.lcsim.hps.evio.TriggerData;
 import org.lcsim.hps.recon.ecal.HPSEcalCluster;

+import org.lcsim.hps.recon.tracking.HPSTrack;

 import org.lcsim.hps.recon.tracking.SvtTrackExtrapolator;
 import org.lcsim.hps.recon.tracking.TrackUtils;
 import org.lcsim.hps.recon.vertexing.SimpleVertexer;

@@ -432,6 +433,8 @@

             br_line+="vtx_fr_x/F:vtx_fr_y/F:vtx_fr_z/F:";
             br_line+="trk1_conv_x/F:trk1_conv_y/F:trk1_conv_z/F:";
             br_line+="trk2_conv_x/F:trk2_conv_y/F:trk2_conv_z/F:";

+            br_line+="trk1_fr_conv_x/F:trk1_fr_conv_y/F:trk1_fr_conv_z/F:";
+            br_line+="trk2_fr_conv_x/F:trk2_fr_conv_y/F:trk2_fr_conv_z/F:";

             br_line+="cl1_E/F:cl1_ix/I:cl1_iy/I:";
             br_line+="cl2_E/F:cl2_ix/I:cl2_iy/I:";
             br_line+="cl3_E/F:cl3_ix/I:cl3_iy/I:";

@@ -597,6 +600,12 @@

         this.vertexer.extrapolator().setTrack(trk2);
         posAtConverter = this.vertexer.extrapolator().extrapolateTrack(SvtTrackExtrapolator.HARP_POSITION);        
         printWriter.format("%5.5f %5.5f %5.5f ", posAtConverter.z(),posAtConverter.x(),posAtConverter.y()); //note rotation from JLab->tracking

+        HPSTrack hpstrk1 = new HPSTrack(helix1);
+        Hep3Vector posAtConverterFringe1 = hpstrk1.getPositionAtZMap(100., SvtTrackExtrapolator.HARP_POSITION, 5.0)[0];
+        printWriter.format("%5.5f %5.5f %5.5f ", posAtConverterFringe1.z(),posAtConverterFringe1.x(),posAtConverterFringe1.y()); //note rotation from JLab->tracking
+        HPSTrack hpstrk2 = new HPSTrack(helix2);
+        Hep3Vector posAtConverterFringe2 = hpstrk2.getPositionAtZMap(100., SvtTrackExtrapolator.HARP_POSITION, 5.0)[0];
+        printWriter.format("%5.5f %5.5f %5.5f ", posAtConverterFringe2.z(),posAtConverterFringe2.x(),posAtConverterFringe2.y()); //note rotation from JLab->tracking

         int ncl_t=0; int ncl_b=0;
         for(int i=0;i<3;++i) {
             if(clusters==null) {