Print

Print


Commit in lcsim-contrib/src/main/java/org/lcsim/contrib/Cassell/recon/analysis on MAIN
QqbarAnalysisDriver.java+91.1 -> 1.2
Add mean90 plot

lcsim-contrib/src/main/java/org/lcsim/contrib/Cassell/recon/analysis
QqbarAnalysisDriver.java 1.1 -> 1.2
diff -u -r1.1 -r1.2
--- QqbarAnalysisDriver.java	4 Nov 2010 15:12:56 -0000	1.1
+++ QqbarAnalysisDriver.java	7 Nov 2010 12:06:54 -0000	1.2
@@ -217,6 +217,7 @@
    }
    protected void suspend()
    {
+      IDataPointSet[][] dpsM = new IDataPointSet[rplnames.size()][cutname.length];
       IDataPointSet[][] dpsE = new IDataPointSet[rplnames.size()][cutname.length];
       IDataPointSet[][] dpsE2 = new IDataPointSet[rplnames.size()][cutname.length];
       IDataPointSet[][] dpsE3 = new IDataPointSet[rplnames.size()][cutname.length];
@@ -224,6 +225,7 @@
       {
          for(int j=0;j<cutname.length;j++)
          {
+            dpsM[i][j] = dpsf.create(rplnames.get(i)+":"+cutname[j]+":mean90 vs cmE",rplnames.get(i)+":"+cutname[j]+":m90 vs cmE",2);
             dpsE[i][j] = dpsf.create(rplnames.get(i)+":"+cutname[j]+":alpha90 vs jetE",rplnames.get(i)+":"+cutname[j]+":alph vs jetE",2);
             dpsE2[i][j] = dpsf.create(rplnames.get(i)+":"+cutname[j]+":dEoE vs cmE",rplnames.get(i)+":"+cutname[j]+":res vs cmE",2);
             dpsE3[i][j] = dpsf.create(rplnames.get(i)+":"+cutname[j]+":dEoE vs jetE",rplnames.get(i)+":"+cutname[j]+":res vs jetE",2);
@@ -251,11 +253,14 @@
              double dEoE = r90/(ZE+m90);
              double err = (1.1*r90/Math.sqrt(1.8*ent))/Math.sqrt(ZE+m90);
              double err2 = (1.1*r90/Math.sqrt(1.8*ent))/(ZE+m90);
+             double errm = 1.1*r90/Math.sqrt(.9*ent);
              double jde = r90*Math.sqrt(2.)/(ZE+m90);
              System.out.println(cutname[i]+" "+ent+" "+m+" "+r+" "+m90+" "+r90+" "+alpha+" "+jde);
+            dpsM[k][i].addPoint();
             dpsE[k][i].addPoint();
             dpsE2[k][i].addPoint();
             dpsE3[k][i].addPoint();
+            IDataPoint dpm = dpsM[k][i].point(iE);
             IDataPoint dp = dpsE[k][i].point(iE);
             IDataPoint dp2 = dpsE2[k][i].point(iE);
             IDataPoint dp3 = dpsE3[k][i].point(iE);
@@ -267,6 +272,10 @@
             dp2.coordinate(1).setValue(dEoE);
             dp2.coordinate(1).setErrorPlus(err2);
             dp2.coordinate(1).setErrorMinus(err2);
+            dpm.coordinate(0).setValue(ZE);
+            dpm.coordinate(1).setValue(m90);
+            dpm.coordinate(1).setErrorPlus(errm);
+            dpm.coordinate(1).setErrorMinus(errm);
             dp3.coordinate(0).setValue(ZE/2.);
             dp3.coordinate(1).setValue(dEoE*Math.sqrt(2.));
             dp3.coordinate(1).setErrorPlus(err2*Math.sqrt(2.));
CVSspam 0.2.8