lcsim-contrib/src/main/java/org/lcsim/contrib/SteveMagill
diff -u -r1.1 -r1.2
--- MCRPdijetm.java 12 Aug 2009 18:33:49 -0000 1.1
+++ MCRPdijetm.java 14 Feb 2012 19:34:00 -0000 1.2
@@ -25,6 +25,8 @@
{
super.process(event);
+ // get delta M correction
+// List<Double> Delm = event.get(Double.class,"MCDelM");
// get dijet mass if 2 jets
List<ReconstructedParticle> jets = event.get(ReconstructedParticle.class,_jname);
double jpx = 0.;
@@ -41,7 +43,8 @@
{
Etot += jet.getEnergy();
Hep3Vector jmom = jet.getMomentum();
- aida.cloud1D("MCParticle Jet pT").fill(Math.sqrt(jmom.x()*jmom.x()+jmom.y()*jmom.y()));
+// aida.cloud1D("MC RP pT per Jet").fill(Math.sqrt(jmom.x()*jmom.x()+jmom.y()*jmom.y()));
+ aida.cloud1D("MC RP Energy per Jet").fill(jet.getEnergy());
jpx += jmom.x();
jpy += jmom.y();
jpz += jmom.z();
@@ -65,9 +68,15 @@
}
double ptot = Math.sqrt(jpx*jpx+jpy*jpy+jpz*jpz);
double jinv = Math.sqrt(Etot*Etot-ptot*ptot+jm2);
+// double mcor = 0.;
+// for (Double dlm : Delm)
+// {
+// mcor = dlm.doubleValue();
+// }
aida.cloud1D("Total MC Jet Mass").fill(jinv);
- aida.cloud1D("Total MC E jets").fill(Etot);
+ aida.cloud1D("Total MC Jet E").fill(Etot);
// System.out.println("Total MC Jet Mass " + jinv);
+// System.out.println("DeltaM Mass Correction " + mcor);
}
public void setInputJets(String jname)