diff -u -r1.10 -r1.11
--- TrackingReconstructionPlots.java	4 Jun 2012 20:17:45 -0000	1.10
+++ TrackingReconstructionPlots.java	29 Jun 2012 18:58:02 -0000	1.11
@@ -80,9 +80,11 @@

     IPlotter top1;
     IPlotter top2;
     IPlotter top3;

+    IPlotter top4;

     IPlotter bot1;
     IPlotter bot2;
     IPlotter bot3;

+    IPlotter bot4;

     double zEcal = 1500;
     double zAtDownStrPairSpec = 914.0; //mm
     double zAtColl = -1500;

@@ -396,7 +398,7 @@

 
         IHistogram1D topdistX = aida.histogram1D("Top deltaX", 50, -400, 400);
         IHistogram1D topdistY = aida.histogram1D("Top deltaY", 50, -40, 40);

-

+        

 
         IHistogram2D topxEcalVsTrk = aida.histogram2D("Top X ECal Vs Track", 100, -400, 400, 100, -400, 400);
         IHistogram2D topyEcalVsTrk = aida.histogram2D("Top Y ECal Vs Track", 100, 0, 100, 100, 0, 100);

@@ -425,8 +427,8 @@

 
         IHistogram1D BottomdistX = aida.histogram1D("Bottom deltaX", 50, -400, 400);
         IHistogram1D BottomdistY = aida.histogram1D("Bottom deltaY", 50, -40, 40);

-
-

+        
+        

         IHistogram2D BottomxEcalVsTrk = aida.histogram2D("Bottom X ECal Vs Track", 100, -400, 400, 100, -400, 400);
         IHistogram2D BottomyEcalVsTrk = aida.histogram2D("Bottom Y ECal Vs Track", 100, -100, 0, 100, -100, 0);
 

@@ -438,7 +440,45 @@

         bot2.region(2).plot(BottomxEcalVsTrk);
         bot2.region(5).plot(BottomyEcalVsTrk);
 

+        
+//   ******************************************************************

+        top3 = fac.createPlotterFactory().create("Top ECal Plots");
+        top3.setTitle("Top ECal More Correlations");
+        IPlotterStyle stop3 = top3.style();
+        stop3.dataStyle().fillStyle().setColor("green");
+        stop3.dataStyle().errorBarStyle().setVisible(false);
+        stop3.setParameter("hist2DStyle", "colorMap");
+        stop3.dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow");
+        top3.createRegions(1, 2);
+        topFrame.addPlotter(top3);
+        
+        IHistogram2D topdistXvsX = aida.histogram2D("Top deltaX vs X", 51,-400,400, 25, -400, 400);
+        IHistogram2D topdistYvsY = aida.histogram2D("Top deltaY vs Y", 51,0,100, 25, -40, 40);
+        
+        top3.region(0).plot(topdistXvsX);
+        top3.region(1).plot(topdistYvsY);
+        
+        
+        bot3 = fac.createPlotterFactory().create("Bottom ECal Plots");
+        bot3.setTitle("Bottom ECal More Correlations");
+        IPlotterStyle sbot3 = bot3.style();
+        sbot3.dataStyle().fillStyle().setColor("green");
+        sbot3.dataStyle().errorBarStyle().setVisible(false);
+        sbot3.setParameter("hist2DStyle", "colorMap");
+        sbot3.dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow");
+        bot3.createRegions(1, 2);
+        bottomFrame.addPlotter(bot3);
+        
+        
+        IHistogram2D botdistXvsX = aida.histogram2D("Bottom deltaX vs X", 51,-400,400, 25, -400, 400);
+        IHistogram2D botdistYvsY = aida.histogram2D("Bottom deltaY vs Y", 51,-100,0, 25, -40, 40);
+        
+        
+        bot3.region(0).plot(botdistXvsX);
+        bot3.region(1).plot(botdistYvsY);
+        
+        

 //   ******************************************************************
 
 

@@ -831,6 +871,9 @@

 
             //           if (clust != null) {
             if (clust != null) {

+                
+                posAtEcal = extend.extrapolateTrack(clust.getPosition()[2]);//.positionAtEcal();
+

                 aida.histogram2D("Energy Vs Momentum").fill(clust.getEnergy(), trk.getPX() * 1000.0);
                 aida.histogram1D("Energy Over Momentum").fill(clust.getEnergy() / (trk.getPX() * 1000.0));
                 aida.histogram1D("deltaX").fill(clust.getPosition()[0] - posAtEcal.y());

@@ -847,6 +890,8 @@

                     aida.histogram1D("Top Energy Over Momentum").fill(clust.getEnergy() / (trk.getPX() * 1000.0));
                     aida.histogram1D("Top deltaX").fill(clust.getPosition()[0] - posAtEcal.y());
                     aida.histogram1D("Top deltaY").fill(clust.getPosition()[1] - posAtEcal.z());

+                    aida.histogram2D("Top deltaX vs X").fill(clust.getPosition()[0],clust.getPosition()[0] - posAtEcal.y());
+                    aida.histogram2D("Top deltaY vs Y").fill(clust.getPosition()[1],clust.getPosition()[1] - posAtEcal.z());

                     aida.histogram2D("Top X ECal Vs Track").fill(clust.getPosition()[0], posAtEcal.y());
                     aida.histogram2D("Top Y ECal Vs Track").fill(clust.getPosition()[1], posAtEcal.z());
                 } else {

@@ -854,6 +899,8 @@

                     aida.histogram1D("Bottom Energy Over Momentum").fill(clust.getEnergy() / (trk.getPX() * 1000.0));
                     aida.histogram1D("Bottom deltaX").fill(clust.getPosition()[0] - posAtEcal.y());
                     aida.histogram1D("Bottom deltaY").fill(clust.getPosition()[1] - posAtEcal.z());

+                    aida.histogram2D("Bottom deltaX vs X").fill(clust.getPosition()[0],clust.getPosition()[0] - posAtEcal.y());
+                    aida.histogram2D("Bottom deltaY vs Y").fill(clust.getPosition()[1],clust.getPosition()[1] - posAtEcal.z());

                     aida.histogram2D("Bottom X ECal Vs Track").fill(clust.getPosition()[0], posAtEcal.y());
                     aida.histogram2D("Bottom Y ECal Vs Track").fill(clust.getPosition()[1], posAtEcal.z());
                 }