Commit in lcsim-contrib/src/main/java/org/lcsim/contrib/SteveMagill on MAIN | |||
MCRPdijetm.java | +11 | -2 | 1.1 -> 1.2 |
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)
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