hps-java/src/main/java/org/lcsim/hps/users/phansson
diff -u -r1.12 -r1.13
--- MPAlignmentParameters.java 12 Oct 2012 05:59:36 -0000 1.12
+++ MPAlignmentParameters.java 17 Oct 2012 23:00:38 -0000 1.13
@@ -79,7 +79,6 @@
boolean _includeMS = false;
private AIDAFrame plotterFrame;
- private AIDAFrame plotterFrameHTH;
private AIDAFrame plotterFrameSimHit;
private AIDAFrame plotterFrameSummary;
private boolean hideFrame = false;
@@ -87,16 +86,8 @@
private IAnalysisFactory af = aida.analysisFactory();
IDataPointSet dps_t;
IDataPointSet dps_b;
- IDataPointSet dps_hth_y_b;
- IDataPointSet dps_hth_y_t;
- IDataPointSet dps_hth_z_b;
- IDataPointSet dps_hth_z_t;
IDataPointSet dps_pull_t;
IDataPointSet dps_pull_b;
- IDataPointSet dps_pull_hth_y_b;
- IDataPointSet dps_pull_hth_y_t;
- IDataPointSet dps_pull_hth_z_b;
- IDataPointSet dps_pull_hth_z_t;
IPlotter plotter_resuydiff_t;
IPlotter plotter_resuydiff_b;
IDataPointSet[] dps_resuydiff_t = new IDataPointSet[10];
@@ -160,6 +151,10 @@
HelicalTrackHit htc = (HelicalTrackHit) hit;
+ if(!(htc instanceof HelicalTrackCross)) {
+ if(_DEBUG) System.out.println(this.getClass().getSimpleName() + ": this hit is not a cross");
+ continue;
+ }
HelicalTrackCross cross = (HelicalTrackCross) htc;
cross.setTrackDirection(_trk);
@@ -167,15 +162,15 @@
double msdz = msmap.get(htc).dz();
double sHit = _trk.PathMap().get(htc);
- calculateTrackHitResidual(htc,msmap);
List<HelicalTrackStrip> clusterlist = cross.getStrips();
if(_DEBUG) System.out.println(this.getClass().getSimpleName() + ": Loop over " + clusterlist.size() + " clusterlist for this hitontrack");
for (HelicalTrackStrip cl : clusterlist) {
+
if(_DEBUG) System.out.println(this.getClass().getSimpleName() + ": cluster size " + cl.rawhits().size());
- //if(cl.rawhits().size()>1) continue;
+
if(type=="GLOBAL") {
//CalculateLocalDerivativesGLOBAL(cl);
@@ -188,10 +183,7 @@
if(_doSimHitsResiduals) CalculateResidualSim(cl, msdrphi, msdz,simhits);
CalculateLocalDerivatives(cl);
- CalculateGlobalDerivatives(cl);
-
- // CalculateResidual(cl, 0,0);
- //PrintStripResiduals(cl);
+ CalculateGlobalDerivatives(cl);
PrintStripResidualsNew(cl);
}
}
@@ -1047,137 +1039,8 @@
}
- public void calculateTrackHitResidual(HelicalTrackHit hth,Map<HelicalTrackHit, MultipleScatter> msmap) {
-
- double msdrphi = msmap.get(hth).drphi();
- double msdz = msmap.get(hth).dz();
-
- //Calculate the residuals that are being used in the track fit
-
- //Start with the bendplane y
- double drphi_res = hth.drphi();
- double wrphi = Math.sqrt(drphi_res*drphi_res + msdrphi*msdrphi);
- double s = _trk.PathMap().get(hth);
- Hep3Vector posOnHelix = HelixUtils.PointOnHelix(_trk, s);
- double resy = hth.y() - posOnHelix.y();
- double resypull = this._includeMS ? resy/wrphi : resy/drphi_res;
-
- String side = "top";
- if (hth.getPosition()[2]<0) side = "bottom";
- int layer = hth.Layer();
- if(layer % 2 == 0) {
- layer = layer/2;
- } else {
- layer = (layer+1)/2;
- }
-
-
- aida.histogram1D("hthres_y_layer" + layer + "_" + side).fill(resy);
- aida.histogram1D("hthrespull_y_layer" + layer + "_" + side).fill(resypull);
- aida.histogram1D("hthreserr_y_layer" + layer + "_" + side).fill(wrphi);
- aida.histogram2D("hthrespull_slope_y_layer" + layer + "_" + side).fill(_trk.slope(),resypull);
-
-
- //Now the residual for the "measurement" direction z
- double resz = hth.z() - posOnHelix.z();
- //double dz = HitUtils.zres(hth, msmap, _trk);
-
- double slope = hth.z() / hth.r();
- // Don't use the helix slope if the magnitude of the slope is smaller than its uncertainty
- if (Math.abs(_trk.slope()) > _trk.getSlopeError()) slope = _trk.slope();
- // Take the resolution uncertainty to be dr * |slope|
- double dzres = hth.dr() * Math.abs(slope);
- // Combine resolution and multiple scattering uncertainties in quadrature
- double dz = dzres; //Math.sqrt(dzres*dzres + msdz*msdz);
-
- aida.histogram1D("hthres_z_layer" + layer + "_" + side).fill(resz);
- aida.histogram1D("hthrespull_z_layer" + layer + "_" + side).fill(resz/dz);
- aida.histogram1D("hthreserr_z_layer" + layer + "_" + side).fill(dz);
- aida.histogram1D("hthreserrms_z_layer" + layer + "_" + side).fill(msdz);
- aida.histogram1D("hthreserrres_z_layer" + layer + "_" + side).fill(dzres);
-
-
- //if(this._DEBUG) {
- if(resy<-0.15) {
- System.out.println(this.getClass().getSimpleName() + ": LARGE TRACK HIT RESIDUAL yres " + resy + " zres " + resz );
- double[] pos = hth.getPosition();
- double[] pos_cor = hth.getCorrectedPosition().v();
- System.out.println(this.getClass().getSimpleName() + ": HelicalTrackHit at [" + pos[0] + "," + pos[1] + "," + pos[2] + "] cor : [" + pos_cor[0] + "," + pos_cor[1] + "," + pos_cor[2] + "]");
- System.exit(1);
- }
- //}
-
-
- //
-// // Store the coordinates and errors for the line fit
-// for (int i = 0; i < npix; i++) {
-// HelicalTrackHit hit = pixel_hits.get(i);
-// z[i] = hit.z();
-// dz[i] = HitUtils.zres(hit, msmap, oldhelix);
-// s[i] = smap.get(hit);
-// }
-//
-//
-
-
- //From HelicalTrackFitter.java
- // // Store the hit coordinates and weights in arrays for the circle fitter
-// for (int i = 0; i < nc; i++) {
-// HelicalTrackHit hit = circle_hits.get(i);
-// // Store the hit position
-// x[i] = hit.x();
-// y[i] = hit.y();
-// // Find the weight (= 1/uncertainty^2) for this hit
-// // First get the multiple scattering uncertainty
-// double drphi_ms = 0.;
-// if (msmap.containsKey(hit)) {
-// drphi_ms = msmap.get(hit).drphi();
-// }
-// // Get the hit resolution and combine uncertainties in quadrature
-// double drphi_res = hit.drphi();
-// wrphi[i] = 1. / (drphi_res * drphi_res + drphi_ms * drphi_ms);
-// }
- }
-
-/*
- public void AddTarget(double beamdy, double beamdz) {
- double[][] dfdp = new double[3][1];
- double d0 = _trk.dca();
- double z0 = _trk.z0();
- double slope = _trk.slope();
- double phi0 = _trk.phi0();
- double R = _trk.R();
- double xint = 0; //target
- double s = HelixUtils.PathToXPlane(_trk, xint, smax, _nlc).get(0);
- Hep3Vector ptAtTarget = HelixUtils.PointOnHelix(_trk, s);
- double[] mydydq = dydq(R, d0, phi0, xint, s);
- double[] mydzdq = dzdq(R, d0, phi0, xint, slope, s);
- _resid[0] = ptAtTarget.z();
- _resid[1] = ptAtTarget.y();
- _resid[2] = ptAtTarget.x();
- _error[0] = beamdz;
- _error[1] = beamdy;
- _error[2] = 666;
- dfdp[0][0] = 1;
- dfdp[1][0] = 0;
- dfdp[2][0] = 0;
- _dfdp.add(FillMatrix(dfdp, 3, 1));
- _globalLabel[0] = 666;
- pWriter.printf("%4d\n", 666);
- pWriter.printf("%5.5e %5.5e %5.5e\n", _resid[0], _resid[1], _resid[2]);
- pWriter.printf("%5.5e %5.5e %5.5e\n", _error[0], _error[1], _error[2]);
- for (int i = 0; i < _nlc; i++) {
- pWriter.printf("%5.5e %5.5e -1.0\n", mydzdq[i], mydydq[i]);
- }
- for (int j = 0; j < _ngl; j++) {
- pWriter.printf("%5.5e %5.5e %5.5e %5d\n", _dfdp.get(_dfdp.size()-1).e(0, j), _dfdp.get(_dfdp.size()-1).e(1, j), _dfdp.get(_dfdp.size()-1).e(2, j), _globalLabel.get(_globalLabel.size()-1)[j]);
- }
-
- }
-*/
-
private void PrintStripResidualsNew(HelicalTrackStrip strip) {
if (_DEBUG) {
@@ -1398,9 +1261,6 @@
plotterFrame = new AIDAFrame();
plotterFrame.setTitle("Residuals");
- plotterFrameHTH = new AIDAFrame();
- plotterFrameHTH.setTitle("Residuals HTH");
-
plotterFrameSummary = new AIDAFrame();
plotterFrameSummary.setTitle("Summary");
@@ -1522,167 +1382,15 @@
plotterFrame.addPlotter(plotter_respull);
plotterFrame.addPlotter(plotter_respull_slope);
}
- }
-
- for (int s=0;s<2;++s) {
- IPlotter plotter_hthresy = af.createPlotterFactory().create();
- plotter_hthresy.createRegions(5,1,0);
- plotter_hthresy.setTitle("hthres_y " + side[s]);
- IPlotter plotter_hthresypull = af.createPlotterFactory().create();
- plotter_hthresypull.createRegions(5,1,0);
- plotter_hthresypull.setTitle("hthrespull_y " + side[s]);
- IPlotter plotter_hthresyerr = af.createPlotterFactory().create();
- plotter_hthresyerr.createRegions(5,1,0);
- plotter_hthresyerr.setTitle("hthreserr_y " + side[s]);
- IPlotter plotter_hthresz = af.createPlotterFactory().create();
- plotter_hthresz.createRegions(5,1,0);
- plotter_hthresz.setTitle("hthres_z " + side[s]);
- IPlotter plotter_hthreszerr = af.createPlotterFactory().create();
- plotter_hthreszerr.createRegions(5,1,0);
- plotter_hthreszerr.setTitle("hthreserr_z " + side[s]);
- IPlotter plotter_hthreszpull = af.createPlotterFactory().create();
- plotter_hthreszpull.createRegions(5,1,0);
- plotter_hthreszpull.setTitle("hthrespull_z " + side[s]);
- IPlotter plotter_hthreszerrms = af.createPlotterFactory().create();
- plotter_hthreszerrms.createRegions(5,1,0);
- plotter_hthreszerrms.setTitle("hthreserrms_z " + side[s]);
- IPlotter plotter_hthreszerrres = af.createPlotterFactory().create();
- plotter_hthreszerrres.createRegions(5,1,0);
- plotter_hthreszerrres.setTitle("hthreserrres_z " + side[s]);
-
- IPlotter plotter_hthresypullslope = af.createPlotterFactory().create();
- plotter_hthresypullslope.createRegions(5,1,0);
- plotter_hthresypullslope.setTitle("hthrespullslope_y " + side[s]);
-
- double xbins_hthy_res[][] = {
- {-0.05,0.05},
- {-0.05,0.05},
- {-0.05,0.05},
- {-0.2,0.2},
- {-0.05,0.05}
-// {-0.4,0.4},
-// {-0.5,0.5},
-// {-1.5,1.5},
-// {-2.5,2.5},
-// {-2.5,2.5}
- };
- double xbins_hthy_reserr[][] = {
- {0,0.5},
- {0,0.5},
- {0,1},
- {0,2},
- {0,2}
- };
- double xbins_hthy_respull[][] = {
- {-4,4},
- {-4,4},
- {-4,4},
- {-4,4},
- {-4,4}
- };
- double xbins_hthz_res[][] = {
- {-0.01,0.01},
- {-0.01,0.01},
- {-0.01,0.01},
- {-0.01,0.01},
- {-0.01,0.01}
-// {-0.03,0.03},
-// {-0.8,0.8},
-// {-0.8,0.8},
-// {-2,2},
-// {-2,2}
- };
- double xbins_hthz_reserr[][] = {
- {0,0.05},
- {0,0.5},
- {0,0.7},
- {0,1.5},
- {0,2}
- };
- double xbins_hthz_respull[][] = {
- {-4,4},
- {-4,4},
- {-4,4},
- {-4,4},
- {-4,4}
- };
+
-
- for (int iLayer=1;iLayer<6;++iLayer) {
-
-
- IHistogram h = aida.histogram1D("hthres_y_layer" + (iLayer) + "_" + side[s] , 50, xbins_hthy_res[iLayer-1][0], xbins_hthy_res[iLayer-1][1]);
- IHistogram h1 = aida.histogram1D("hthreserr_y_layer" + (iLayer) + "_" + side[s] , 50, xbins_hthy_reserr[iLayer-1][0], xbins_hthy_reserr[iLayer-1][1]);
- IHistogram h2 = aida.histogram1D("hthrespull_y_layer" + (iLayer) + "_" + side[s] , 50, xbins_hthy_respull[iLayer-1][0], xbins_hthy_respull[iLayer-1][1]);
-
- IHistogram hz = aida.histogram1D("hthres_z_layer" + (iLayer) + "_" + side[s] , 50, xbins_hthz_res[iLayer-1][0], xbins_hthz_res[iLayer-1][1]);
- IHistogram hz1 = aida.histogram1D("hthreserr_z_layer" + (iLayer) + "_" + side[s] , 50, xbins_hthz_reserr[iLayer-1][0], xbins_hthz_reserr[iLayer-1][1]);
- IHistogram hz2 = aida.histogram1D("hthrespull_z_layer" + (iLayer) + "_" + side[s] , 50, xbins_hthz_respull[iLayer-1][0], xbins_hthz_respull[iLayer-1][1]);
- IHistogram hz11 = aida.histogram1D("hthreserrms_z_layer" + (iLayer) + "_" + side[s] , 50, 0, 1);
- IHistogram hz111 = aida.histogram1D("hthreserrres_z_layer" + (iLayer) + "_" + side[s] , 50, 0, 1);
-
- IHistogram h2_slope;
- if("top".equals(side[s])) h2_slope = aida.histogram2D("hthrespull_slope_y_layer" + (iLayer) + "_" + side[s],50,0,0.1, 50, xbins_hthy_respull[iLayer-1][0], xbins_hthy_respull[iLayer-1][1]);
- else h2_slope = aida.histogram2D("hthrespull_slope_y_layer" + (iLayer) + "_" + side[s],50,-0.1,0, 50, xbins_hthy_respull[iLayer-1][0], xbins_hthy_respull[iLayer-1][1]);
-
-
-
-
- int region = (iLayer-1);//*2+iSide;
-
- plotter_hthresy.region(region).plot(h);
- plotter_hthresyerr.region(region).plot(h1);
- plotter_hthresypull.region(region).plot(h2);
- plotter_hthresypullslope.region(region).plot(h2_slope);
-
- plotter_hthresz.region(region).plot(hz);
- plotter_hthreszerr.region(region).plot(hz1);
- plotter_hthreszpull.region(region).plot(hz2);
-
- plotter_hthreszerrms.region(region).plot(hz11);
- plotter_hthreszerrres.region(region).plot(hz111);
-
-
- ((PlotterRegion) plotter_hthresy.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthresy.region(region)).getPlot().setAllowPopupMenus(true);
- ((PlotterRegion) plotter_hthresyerr.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthresyerr.region(region)).getPlot().setAllowPopupMenus(true);
- ((PlotterRegion) plotter_hthresypull.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthresypull.region(region)).getPlot().setAllowPopupMenus(true);
- ((PlotterRegion) plotter_hthresypullslope.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthresypullslope.region(region)).getPlot().setAllowPopupMenus(true);
- plotter_hthresypull.style().dataStyle().fillStyle().setParameter("colorMapScheme", "rainbow");
-
- ((PlotterRegion) plotter_hthresz.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthresz.region(region)).getPlot().setAllowPopupMenus(true);
- ((PlotterRegion) plotter_hthreszerr.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthreszerr.region(region)).getPlot().setAllowPopupMenus(true);
- ((PlotterRegion) plotter_hthreszpull.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthreszpull.region(region)).getPlot().setAllowPopupMenus(true);
-
- ((PlotterRegion) plotter_hthreszerrms.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthreszerrms.region(region)).getPlot().setAllowPopupMenus(true);
- ((PlotterRegion) plotter_hthreszerrres.region(region)).getPlot().setAllowUserInteraction(true);
- ((PlotterRegion) plotter_hthreszerrres.region(region)).getPlot().setAllowPopupMenus(true);
-
-
-
- }
-
- plotterFrameHTH.addPlotter(plotter_hthresy);
- plotterFrameHTH.addPlotter(plotter_hthresyerr);
- plotterFrameHTH.addPlotter(plotter_hthresypull);
- plotterFrameHTH.addPlotter(plotter_hthresypullslope);
- plotterFrameHTH.addPlotter(plotter_hthresz);
- plotterFrameHTH.addPlotter(plotter_hthreszerr);
- plotterFrameHTH.addPlotter(plotter_hthreszpull);
- plotterFrameHTH.addPlotter(plotter_hthreszerrms);
- plotterFrameHTH.addPlotter(plotter_hthreszerrres);
+
+
}
@@ -1814,7 +1522,7 @@
IPlotter plotter_prf = af.createPlotterFactory().create();
- plotter_prf.createRegions(2,3,0);
+ plotter_prf.createRegions(1,2,0);
plotter_prf.setTitle("<Residuals>");
IDataPointSetFactory dpsf = aida.analysisFactory().createDataPointSetFactory(null);
dps_t = dpsf.create("dps_t", "Mean of u residuals top",2);
@@ -1823,31 +1531,11 @@
plotter_prf.region(0).plot(dps_b);
dps_pull_t = dpsf.create("dps_pull_t", "Mean of u pulls top",2);
- plotter_prf.region(3).plot(dps_pull_t);
+ plotter_prf.region(1).plot(dps_pull_t);
dps_pull_b = dpsf.create("dps_pull_b", "Mean of u pulls bottom",2);
- plotter_prf.region(3).plot(dps_pull_b);
+ plotter_prf.region(1).plot(dps_pull_b);
- dps_hth_y_t = dpsf.create("dps_hth_y_t", "Mean of y residuals top",2);
- plotter_prf.region(1).plot(dps_hth_y_t);
- dps_hth_y_b = dpsf.create("dps_hth_y_b", "Mean of y residuals bottom",2);
- plotter_prf.region(1).plot(dps_hth_y_b);
-
- dps_pull_hth_y_t = dpsf.create("dps_pull_hth_y_t", "Mean of y pulls top",2);
- plotter_prf.region(4).plot(dps_pull_hth_y_t);
- dps_pull_hth_y_b = dpsf.create("dps_pull_hth_y_b", "Mean of y pulls bottom",2);
- plotter_prf.region(4).plot(dps_pull_hth_y_b);
-
- dps_hth_z_t = dpsf.create("dps_hth_z_t", "Mean of z residuals top",2);
- plotter_prf.region(2).plot(dps_hth_z_t);
- dps_hth_z_b = dpsf.create("dps_hth_z_b", "Mean of z residuals bottom",2);
- plotter_prf.region(2).plot(dps_hth_z_b);
-
- dps_pull_hth_z_t = dpsf.create("dps_pull_hth_z_t", "Mean of z pulls top",2);
- plotter_prf.region(5).plot(dps_pull_hth_z_t);
- dps_pull_hth_z_b = dpsf.create("dps_pull_hth_z_b", "Mean of z pulls bottom",2);
- plotter_prf.region(5).plot(dps_pull_hth_z_b);
-
- for(int region=0;region<6;++region) {
+ for(int region=0;region<2;++region) {
((PlotterRegion) plotter_prf.region(region)).getPlot().setAllowUserInteraction(true);
((PlotterRegion) plotter_prf.region(region)).getPlot().setAllowPopupMenus(true);
}
@@ -1901,8 +1589,6 @@
plotterFrame.pack();
plotterFrame.setVisible(!hideFrame);
- plotterFrameHTH.pack();
- plotterFrameHTH.setVisible(!hideFrame);
plotterFrameSummary.pack();
plotterFrameSummary.setVisible(!hideFrame);
@@ -1962,101 +1648,8 @@
dps_pull_b.point(i-1).coordinate(0).setErrorPlus(0);
}
- dps_hth_y_t.clear();
- dps_hth_y_b.clear();
- dps_hth_z_t.clear();
- dps_hth_z_b.clear();
- dps_pull_hth_y_t.clear();
- dps_pull_hth_y_b.clear();
- dps_pull_hth_z_t.clear();
- dps_pull_hth_z_b.clear();
- for(int i=1;i<6;++i) {
-
- double mean = aida.histogram1D("hthres_y_layer" + i + "_bottom").mean();
- double stddev = aida.histogram1D("hthres_y_layer" + i + "_bottom").rms();
- double N = aida.histogram1D("hthres_y_layer" + i + "_bottom").entries();
- double error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_hth_y_b.addPoint();
- dps_hth_y_b.point(i-1).coordinate(1).setValue(mean);
- dps_hth_y_b.point(i-1).coordinate(1).setErrorPlus(error);
- dps_hth_y_b.point(i-1).coordinate(0).setValue(i);
- dps_hth_y_b.point(i-1).coordinate(0).setErrorPlus(0);
-
- mean = aida.histogram1D("hthres_y_layer" + i + "_top").mean();
- stddev = aida.histogram1D("hthres_y_layer" + i + "_top").rms();
- N = aida.histogram1D("hthres_y_layer" + i + "_top").entries();
- error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_hth_y_t.addPoint();
- dps_hth_y_t.point(i-1).coordinate(1).setValue(mean);
- dps_hth_y_t.point(i-1).coordinate(1).setErrorPlus(error);
- dps_hth_y_t.point(i-1).coordinate(0).setValue(i);
- dps_hth_y_t.point(i-1).coordinate(0).setErrorPlus(0);
-
- mean = aida.histogram1D("hthres_z_layer" + i + "_top").mean();
- stddev = aida.histogram1D("hthres_z_layer" + i + "_top").rms();
- N = aida.histogram1D("hthres_z_layer" + i + "_top").entries();
- error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_hth_z_t.addPoint();
- dps_hth_z_t.point(i-1).coordinate(1).setValue(mean);
- dps_hth_z_t.point(i-1).coordinate(1).setErrorPlus(error);
- dps_hth_z_t.point(i-1).coordinate(0).setValue(i);
- dps_hth_z_t.point(i-1).coordinate(0).setErrorPlus(0);
-
- mean = aida.histogram1D("hthres_z_layer" + i + "_bottom").mean();
- stddev = aida.histogram1D("hthres_z_layer" + i + "_bottom").rms();
- N = aida.histogram1D("hthres_z_layer" + i + "_bottom").entries();
- error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_hth_z_b.addPoint();
- dps_hth_z_b.point(i-1).coordinate(1).setValue(mean);
- dps_hth_z_b.point(i-1).coordinate(1).setErrorPlus(error);
- dps_hth_z_b.point(i-1).coordinate(0).setValue(i);
- dps_hth_z_b.point(i-1).coordinate(0).setErrorPlus(0);
-
-
- mean = aida.histogram1D("hthrespull_y_layer" + i + "_bottom").mean();
- stddev = aida.histogram1D("hthrespull_y_layer" + i + "_bottom").rms();
- N = aida.histogram1D("hthrespull_y_layer" + i + "_bottom").entries();
- error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_pull_hth_y_b.addPoint();
- dps_pull_hth_y_b.point(i-1).coordinate(1).setValue(mean);
- dps_pull_hth_y_b.point(i-1).coordinate(1).setErrorPlus(error);
- dps_pull_hth_y_b.point(i-1).coordinate(0).setValue(i);
- dps_pull_hth_y_b.point(i-1).coordinate(0).setErrorPlus(0);
-
- mean = aida.histogram1D("hthrespull_y_layer" + i + "_top").mean();
- stddev = aida.histogram1D("hthrespull_y_layer" + i + "_top").rms();
- N = aida.histogram1D("hthrespull_y_layer" + i + "_top").entries();
- error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_pull_hth_y_t.addPoint();
- dps_pull_hth_y_t.point(i-1).coordinate(1).setValue(mean);
- dps_pull_hth_y_t.point(i-1).coordinate(1).setErrorPlus(error);
- dps_pull_hth_y_t.point(i-1).coordinate(0).setValue(i);
- dps_pull_hth_y_t.point(i-1).coordinate(0).setErrorPlus(0);
-
- mean = aida.histogram1D("hthrespull_z_layer" + i + "_top").mean();
- stddev = aida.histogram1D("hthrespull_z_layer" + i + "_top").rms();
- N = aida.histogram1D("hthrespull_z_layer" + i + "_top").entries();
- error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_pull_hth_z_t.addPoint();
- dps_pull_hth_z_t.point(i-1).coordinate(1).setValue(mean);
- dps_pull_hth_z_t.point(i-1).coordinate(1).setErrorPlus(error);
- dps_pull_hth_z_t.point(i-1).coordinate(0).setValue(i);
- dps_pull_hth_z_t.point(i-1).coordinate(0).setErrorPlus(0);
-
- mean = aida.histogram1D("hthrespull_z_layer" + i + "_bottom").mean();
- stddev = aida.histogram1D("hthrespull_z_layer" + i + "_bottom").rms();
- N = aida.histogram1D("hthrespull_z_layer" + i + "_bottom").entries();
- error = N >0 ? stddev/Math.sqrt(N) : 0;
- dps_pull_hth_z_b.addPoint();
- dps_pull_hth_z_b.point(i-1).coordinate(1).setValue(mean);
- dps_pull_hth_z_b.point(i-1).coordinate(1).setErrorPlus(error);
- dps_pull_hth_z_b.point(i-1).coordinate(0).setValue(i);
- dps_pull_hth_z_b.point(i-1).coordinate(0).setErrorPlus(0);
-
- }
-
/*
for(int iLayer=1;iLayer<11;++iLayer) {
dps_resuydiff_t[iLayer-1].clear();