diff -u -r1.1 -r1.2
--- SRPdijetm.java	12 Aug 2009 18:33:59 -0000	1.1
+++ SRPdijetm.java	14 Feb 2012 19:35:58 -0000	1.2
@@ -15,6 +15,7 @@

 {
     private String _jlinkmap;
     private String _scinjets;

+    private String _dmcorr;

     private AIDA aida = AIDA.defaultInstance();
 
     public SRPdijetm()

@@ -26,6 +27,8 @@

     {
         super.process(event);
 

+        //  get mass correction from tracks
+        List<Double> Delm = event.get(Double.class, _dmcorr);

         // array of ReconstructedParticles to add to the event...
         List<ReconstructedParticle> rpList = new ArrayList<ReconstructedParticle>();
 

@@ -74,7 +77,7 @@

             aida.cloud1D("C over S ratio per jet").fill(scrat);
 //           System.out.println("CS ratio " + scrat);
 //            double sjetCE = jet.getEnergy()/(.555-.200*scrat+.643*scrat*scrat);

-            double sjetCE = jet.getEnergy()/(.480+.628*scrat-1.085*scrat*scrat+.975*scrat*scrat*scrat);

+            double sjetCE = jet.getEnergy()/(.587+.173*scrat-.291*scrat*scrat+.531*scrat*scrat*scrat);

 //            double sjetCE = jet.getEnergy()/(.68+.31*scrat);
             aida.cloud1D("Correction Factor for Jets").fill(jet.getEnergy()/sjetCE);
             Hep3Vector jmom = jet.getMomentum();

@@ -114,7 +117,12 @@

         }
         double ptot = Math.sqrt(jpx*jpx+jpy*jpy+jpz*jpz);
         double jinv = Math.sqrt(Etot*Etot-ptot*ptot+mjtot2);

-        aida.cloud1D("Ceren Corrected Scin Jet Mass").fill(jinv);

+        double mcor = 0.;
+        for (Double dlm : Delm)
+        {
+            mcor = dlm.doubleValue();
+        }
+        aida.cloud1D("Ceren Corrected Scin Jet Mass").fill(jinv-mcor);

         aida.cloud1D("Ceren Corrected Scin Jet ETot").fill(Etot);
 //        System.out.println("Total S Jet Mass " + jinv);
         }

@@ -129,5 +137,10 @@

   {
       _scinjets = scinjets;
   }

+
+    public void setDeltaMCorr(String dmcorr)
+    {
+        _dmcorr = dmcorr;
+    }

   
 }