mcd-analysis/src/main/java/org/lcsim/mcd/analysis
diff -u -r1.1 -r1.2
--- MCHiggs.java 19 Dec 2012 17:14:06 -0000 1.1
+++ MCHiggs.java 19 Dec 2012 17:26:34 -0000 1.2
@@ -46,6 +46,7 @@
* 5 == b
* 15 == tau-
* 21 == gluon
+ * 22 == gamma
* 23 == Zo
* 24 == W
*/
@@ -66,13 +67,17 @@
double pt = VecOp.dot(p, z_unit);
double pt_bar = VecOp.dot(p_bar, z_unit);
- aida.cloud1D("b momentum").fill(p_m);
- aida.cloud1D("b transverse momentum").fill(pt);
- aida.cloud1D("b_bar transverse momentum").fill(pt_bar);
+ //aida.cloud1D("b momentum").fill(p_m);
+ aida.histogram1D("b momentum hist",50,62.306,62.313).fill(p_m);
+ //aida.cloud1D("b transverse momentum").fill(pt);
+ aida.histogram1D("b transverse momentum hist", 100,0,63).fill(pt);
+ //aida.cloud1D("b_bar transverse momentum").fill(pt_bar);
+ aida.histogram1D("b transverse momentum hist", 100,0,63).fill(pt_bar);
// b and b_bar are produced with equal energy, opposite momenta
double inv_mass = Math.sqrt(2.*mass*mass + 2.*(en*en-VecOp.dot(p, p_bar)));
- aida.cloud1D("H invariant mass").fill(inv_mass);
+ //aida.cloud1D("H invariant mass").fill(inv_mass);
+ aida.histogram1D("H inv mass hist",50,124.980,124.995).fill(inv_mass);
}
}
aida.cloud1D("num children").fill(children.size());