openmc-dev / openmc / 21489819490

void PointFilter::from_xml(pugi::xml_node node)

  auto bins = get_node_array<double>(node, "bins");

  vector<std::pair<Position, double>> detectors;

  size_t n = bins.size() / 4;

  for (int i = 0; i < n; ++i) {

    Position pos {bins[4 * i], bins[4 * i + 1], bins[4 * i + 2]};

    detectors.push_back(std::make_pair(pos, bins[4 * i + 3]));

  this->set_detectors(detectors);

}

void PointFilter::set_detectors(span<std::pair<Position, double>> detectors)

  detectors_.clear();

  detectors_.reserve(detectors.size());

  for (auto d : detectors) {

    detectors_.push_back(d);

  n_bins_ = detectors_.size();

}

void PointFilter::reset_indices()

  indices_.clear();

  for (auto det : model::active_point_detectors) {

    vector<int32_t> temp;

    for (auto i = 0; i < detectors_.size(); i++) {

      auto [d, r] = detectors_[i];

      if (d == det)

        temp.push_back(i);

    indices_.push_back(temp);

  }

}

void PointFilter::get_all_bins(

  double attenuation = p.wgt() / p.wgt_last();

  for (auto i : indices_[simulation::i_det]) {

    auto [pos, r] = detectors_[i];

    if ((p.r() - pos).norm() < FP_COINCIDENT) {

      match.bins_.push_back(i);

      double distance = (p.r_last() - pos).norm();

      if (distance > r) {

        weight = attenuation / (distance * distance);

        weight = 3.0 * -exprel(-p.macro_xs().total * r) / (r * r);

      match.weights_.push_back(weight);

}

void PointFilter::to_statepoint(hid_t filter_group) const

  Filter::to_statepoint(filter_group);

  vector<double> detectors;

  for (auto [pos, r] : detectors_) {

    detectors.push_back(pos[0]);

    detectors.push_back(pos[1]);

    detectors.push_back(pos[2]);

    detectors.push_back(r);

  write_dataset(filter_group, "bins", detectors);

}

std::string PointFilter::text_label(int bin) const

  auto [pos, r] = detectors_.at(bin);

  return fmt::format("Point: {} {}", pos, r);