FullScanTool.cpp

#include "FullScanTool.h"
#include <algorithm>

namespace SctAnalysis{

FullScanTool::FullScanTool(RawScanResult& raw) : RawScanResult(raw.getHeader(), 
                                   raw.getConfiguration(),
                                   raw.getPoints(),
                                   raw.getScanData(0),
                                   raw.getScanData(1)
                                   ) {}
  
  
FullScanTool::~FullScanTool() throw() {
}

void FullScanTool::add(boost::shared_ptr<const SctData::RawScanResult> raw, double offset){
  if (!raw) throw Sct::IllegalStateError("Null-pointer to raw scan result", __FILE__,__LINE__);
  results.push_back(std::make_pair(raw,offset));
}

// first sort the scan-points by value
struct full_scan_point_data{
  full_scan_point_data(int scan, int point_in_scan, double value) : scan(scan), 
                                    point_in_scan(point_in_scan), 
                                    value(value)
  {
  }
  int scan;
  int point_in_scan;
  double value;
  bool operator< (const full_scan_point_data& other) const {
    return value<other.value;
  }
};



void FullScanTool::combine(){
  using namespace SctData;

  // add all points from all scans to a big list
  std::vector<full_scan_point_data> all_points;
  for (unsigned iscan=0; iscan<results.size(); ++iscan){
    const ScanPoints& pts = results[iscan].first->getPoints();
    for (unsigned ipoint=0; ipoint < pts.getNPoints(); ++ipoint){
      // point has value which differs because of offset
      all_points.push_back( full_scan_point_data(iscan, 
                         ipoint, 
                         pts.getPoint(ipoint)+results[iscan].second )
                );
    }
  }
  
  // sort that list by value
  std::sort(all_points.begin(), all_points.end());
  
  // make a new set of scanpoints
  float points[all_points.size()];
  unsigned events[all_points.size()];
  unsigned errors[all_points.size()];
  
  for (unsigned ipoint=0;ipoint<all_points.size(); ++ipoint){
    const ScanPoints& old_points = results[all_points[ipoint].scan].first->getPoints();
    points[ipoint] = old_points.getPoint(all_points[ipoint].point_in_scan);
    events[ipoint] = old_points.getNEvents(all_points[ipoint].point_in_scan);
    errors[ipoint] = old_points.getNErrorEvents(all_points[ipoint].point_in_scan);
  }

  ScanPoints new_points(all_points.size(), points, events, errors);
  
  getPoints() = new_points;
  
  // RawScanResultIOHelper helper;
  // helper.set(*this);
  
}
}

---