lcsim-contrib/src/main/java/org/lcsim/contrib/Cassell/recon/analysis
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.));