diff -u -r1.1 -r1.2
--- HiggsAnalysisDriver.java	19 Mar 2013 19:58:49 -0000	1.1
+++ HiggsAnalysisDriver.java	19 Mar 2013 20:27:06 -0000	1.2
@@ -189,14 +189,18 @@

     private void makeThrustPlots(List<MCParticle> mcpl, String prefix) {
         EventShape es = getEventShape(mcpl);
         Hep3Vector thrust = es.thrust();

-        Hep3Vector majAxis = es.majorAxis();
-        Hep3Vector minAxis = es.minorAxis();

+        Hep3Vector ThrustAxis = es.thrustAxis();

         double ob = es.oblateness();
         

-        if(thrust.x()>0) aida.cloud1D(prefix+"/Thrust").fill(thrust.x());

+        if(thrust.x()>0) aida.cloud1D(prefix+"/Thrst").fill(thrust.x());

         aida.cloud1D(prefix+"/Major Axis").fill(thrust.y());
         aida.cloud1D(prefix+"/Minor Axis").fill(thrust.z());
         aida.cloud1D(prefix+"/Oblateness").fill(ob);

+        aida.cloud2D(prefix+"/Correl Thrust Major").fill(thrust.x(), thrust.y());
+        aida.cloud2D(prefix+"/Correl Thrust Oblateness").fill(thrust.x(), ob);
+        aida.cloud1D(prefix+"/Thrust Axis CosTheta").fill(VecOp.cosTheta(ThrustAxis));
+        aida.cloud2D(prefix+"/Correl Thrust vs Thrust Axis CosT"
+                ).fill(thrust.x(), VecOp.cosTheta(ThrustAxis));

     }
     
     private void makeJetPlots(EventHeader event, String prefix, String decay,String evt) {