Print

Print


Commit in lcsim-contrib/src/main/java/org/lcsim/contrib/SteveMagill on MAIN
MCRPdijetm.java+11-21.1 -> 1.2

lcsim-contrib/src/main/java/org/lcsim/contrib/SteveMagill
MCRPdijetm.java 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)
CVSspam 0.2.12


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