diff -u -r1.6 -r1.7
--- MCHiggs.java	19 Dec 2012 18:06:27 -0000	1.6
+++ MCHiggs.java	19 Dec 2012 18:13:10 -0000	1.7
@@ -79,6 +79,9 @@

                 //aida.cloud1D("H invariant mass").fill(inv_mass);
                 aida.histogram1D("H inv mass hist",100,124.990,124.995).fill(inv_mass);
                 aida.histogram1D("H inv mass hist big",200,100,130).fill(inv_mass);

+                
+                aida.cloud1D("b phi").fill(VecOp.phi(p));
+                aida.cloud1D("b cos theta").fill(VecOp.cosTheta(p));

             }
             if (particle.getPDGID() == 22) {
                 double mass = particle.getMass();

@@ -99,6 +102,11 @@

                 double eta1 = 0.5*Math.log((p_m+p.z())/(p_m-p.z()));
                 double eta2 = 0.5*Math.log((p2_m+p2.z())/(p2_m-p2.z()));
                 

+                aida.cloud1D("eta1").fill(eta1);
+                aida.cloud1D("eta2").fill(eta2);
+                aida.cloud1D("phi1").fill(phi1);
+                aida.cloud1D("phi2").fill(phi2);
+                

                 double inv_mass = Math.sqrt(
                         2*pt*pt2*(Math.cosh(eta1-eta2) - Math.cos(phi1-phi2)));