Commit in mcd-analysis/src/main/java/org/lcsim/mcd/analysis on MAIN | |||
MCHiggs.java | +9 | -4 | 1.1 -> 1.2 |
Simple driver to grab some MC info from a Higgs .stdhep file and show some plots - changed to histograms
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());
Use REPLY-ALL to reply to list
To unsubscribe from the LCD-CVS list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=LCD-CVS&A=1