lcsim-contrib/src/main/java/org/lcsim/contrib/SteveMagill
diff -u -r1.1 -r1.2
--- MCqqbarME.java 12 Aug 2009 18:29:28 -0000 1.1
+++ MCqqbarME.java 14 Feb 2012 19:20:35 -0000 1.2
@@ -32,6 +32,7 @@
super.process(event);
// get true qqbar E and M and get Nayeli's mass correction
+ List<Double> Ddelm = new ArrayList<Double>();
List<MCParticle> mcl = event.get(MCParticle.class,"MCParticle");
double qqbarM = 0.;
double qqbarE = 0.;
@@ -77,10 +78,12 @@
double ptot = Math.sqrt(mcpx*mcpx+mcpy*mcpy+mcpz*mcpz);
double minv = Math.sqrt(Etot*Etot-ptot*ptot);
double delm = mprime-minv;
+ Ddelm.add(delm);
aida.cloud1D("MC FS Invariant Mass").fill(minv);
aida.cloud1D("MC FS Prime Invariant Mass").fill(mprime);
aida.cloud1D("Delta M Mp M").fill(delm);
// System.out.println("MC FS Invariant Mass " + minv);
+ event.put("MCDelM", Ddelm);
aida.cloud1D("qqbar Parent Mass").fill(qqbarM);
aida.cloud1D("qqbar Parent Energy").fill(qqbarE);