|
[Rivet-svn] r4086 - branches/2012-06-aidarivet/src/Coreblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Dec 17 13:39:53 GMT 2012
Author: hoeth Date: Mon Dec 17 13:39:52 2012 New Revision: 4086 Log: whitespace cleanup Modified: branches/2012-06-aidarivet/src/Core/Analysis.cc Modified: branches/2012-06-aidarivet/src/Core/Analysis.cc ============================================================================== --- branches/2012-06-aidarivet/src/Core/Analysis.cc Mon Dec 17 13:11:37 2012 (r4085) +++ branches/2012-06-aidarivet/src/Core/Analysis.cc Mon Dec 17 13:39:52 2012 (r4086) @@ -293,10 +293,10 @@ IHistogram2D* Analysis::bookHistogram2D(const string& hname, - size_t nxbins, double xlower, double xupper, - size_t nybins, double ylower, double yupper, - const string& title, const string& xtitle, - const string& ytitle, const string& ztitle) { + size_t nxbins, double xlower, double xupper, + size_t nybins, double ylower, double yupper, + const string& title, const string& xtitle, + const string& ytitle, const string& ztitle) { _makeHistoDir(); const string path = histoPath(hname); if (path.find(" ") != string::npos) { @@ -304,7 +304,7 @@ } IHistogram2D* hist = histogramFactory().createHistogram2D(path, title, nxbins, xlower, xupper, - nybins, ylower, yupper); + nybins, ylower, yupper); MSG_TRACE("Made 2D histogram " << hname << " for " << name()); hist->setXTitle(xtitle); hist->setYTitle(ytitle); @@ -315,10 +315,10 @@ IHistogram2D* Analysis::bookHistogram2D(const string& hname, - const vector<double>& xbinedges, - const vector<double>& ybinedges, - const string& title, const string& xtitle, - const string& ytitle, const string& ztitle) { + const vector<double>& xbinedges, + const vector<double>& ybinedges, + const string& title, const string& xtitle, + const string& ytitle, const string& ztitle) { _makeHistoDir(); const string path = histoPath(hname); if (path.find(" ") != string::npos) { @@ -568,7 +568,7 @@ for (int iyBin = 0; iyBin != nyBins; ++iyBin) { // Leaving out factor of binWidth because AIDA's "height" // already includes a width factor. - oldintg += histo->binHeight(ixBin, iyBin); // * histo->axis().binWidth(iBin); + oldintg += histo->binHeight(ixBin, iyBin); // * histo->axis().binWidth(iBin); } if (oldintg == 0.0) { MSG_WARNING("Histo " << hpath << " has null integral during normalization"); @@ -593,21 +593,21 @@ vector<double> x, y, z, ex, ey, ez; for (size_t ix = 0, Nx = histo->xAxis().bins(); ix < Nx; ++ix) for (size_t iy = 0, Ny = histo->yAxis().bins(); iy < Ny; ++iy) { - x.push_back(0.5 * (histo->xAxis().binLowerEdge(ix) + - histo->xAxis().binUpperEdge(ix))); - ex.push_back(histo->xAxis().binWidth(ix)*0.5); - y.push_back(0.5 * (histo->yAxis().binLowerEdge(iy) + - histo->yAxis().binUpperEdge(iy))); - ey.push_back(histo->yAxis().binWidth(iy)*0.5); - - // "Bin height" is a misnomer in the AIDA spec: width is neglected. - // We'd like to do this: y.push_back(histo->binHeight(i) * scale); - z.push_back(histo->binHeight(ix, iy)*scale/ - (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy))); - // "Bin error" is a misnomer in the AIDA spec: width is neglected. - // We'd like to do this: ey.push_back(histo->binError(i) * scale); - ez.push_back(histo->binError(ix, iy)*scale/ - (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy))); + x.push_back(0.5 * (histo->xAxis().binLowerEdge(ix) + + histo->xAxis().binUpperEdge(ix))); + ex.push_back(histo->xAxis().binWidth(ix)*0.5); + y.push_back(0.5 * (histo->yAxis().binLowerEdge(iy) + + histo->yAxis().binUpperEdge(iy))); + ey.push_back(histo->yAxis().binWidth(iy)*0.5); + + // "Bin height" is a misnomer in the AIDA spec: width is neglected. + // We'd like to do this: y.push_back(histo->binHeight(i) * scale); + z.push_back(histo->binHeight(ix, iy)*scale/ + (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy))); + // "Bin error" is a misnomer in the AIDA spec: width is neglected. + // We'd like to do this: ey.push_back(histo->binError(i) * scale); + ez.push_back(histo->binError(ix, iy)*scale/ + (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy))); } string title = histo->title();
More information about the Rivet-svn mailing list |