Beispiel #1
0
/**
  * Overwridden addOrReplace member to attach the name to the workspace when a
 * workspace object is added to the service.
  * This will overwrite one of the same name. If the workspace is group adds or
 * replaces its members.
  * @param name The name of the object
  * @param workspace The shared pointer to the workspace to store
  */
void AnalysisDataServiceImpl::addOrReplace(
    const std::string &name,
    const boost::shared_ptr<API::Workspace> &workspace) {
  verifyName(name);

  // Attach the name to the workspace
  if (workspace)
    workspace->setName(name);
  Kernel::DataService<API::Workspace>::addOrReplace(name, workspace);

  // if a group is added add its members as well
  auto group = boost::dynamic_pointer_cast<WorkspaceGroup>(workspace);
  if (!group)
    return;
  group->observeADSNotifications(true);
  for (size_t i = 0; i < group->size(); ++i) {
    auto ws = group->getItem(i);
    std::string wsName = ws->name();
    // make up a name for an anonymous workspace
    if (wsName.empty()) {
      wsName = name + "_" + boost::lexical_cast<std::string>(i + 1);
    } else if (doesExist(wsName)) { // if ws is already there do nothing
      wsName.clear();
    }
    // add member workspace if needed
    if (!wsName.empty()) {
      addOrReplace(wsName, ws);
    }
  }
}
  /** Construct an output TableWorkspace for fitting result (profile parameters)
    */
  TableWorkspace_sptr RefinePowderInstrumentParameters2::genOutputProfileTable(map<string, Parameter> parameters,
                                                                     double startchi2, double finalchi2)
  {
    // 1. Create TableWorkspace
    TableWorkspace_sptr tablews(new TableWorkspace);

    tablews->addColumn("str", "Name");
    tablews->addColumn("double", "Value");
    tablews->addColumn("str", "FitOrTie");
    tablews->addColumn("double", "Min");
    tablews->addColumn("double", "Max");
    tablews->addColumn("double", "StepSize");
    tablews->addColumn("double", "Error");

    // 2. For chi^2
    addOrReplace(parameters, "Chi2_Init", startchi2);
    addOrReplace(parameters, "Chi2_Result", finalchi2);

    // 3. Set values
    map<string, Parameter>::iterator pariter;
    for (pariter = parameters.begin(); pariter != parameters.end(); ++pariter)
    {
      Parameter& param = pariter->second;
      TableRow newrow = tablews->appendRow();

      string fitortie;
      if (param.fit)
        fitortie = "fit";
      else
        fitortie = "tie";

      newrow << param.name << param.value << fitortie << param.minvalue << param.maxvalue
             << param.stepsize << param.error;
    }

    return tablews;
  }