CellTally.h

Go to the documentation of this file.

/********************************************************************/
/*                  SOFTWARE COPYRIGHT NOTIFICATION                 */
/*                             Cardinal                             */
/*                                                                  */
/*                  (c) 2021 UChicago Argonne, LLC                  */
/*                        ALL RIGHTS RESERVED                       */
/*                                                                  */
/*                 Prepared by UChicago Argonne, LLC                */
/*               Under Contract No. DE-AC02-06CH11357               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*             Prepared by Battelle Energy Alliance, LLC            */
/*               Under Contract No. DE-AC07-05ID14517               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*                 See LICENSE for full restrictions                */
/********************************************************************/

#pragma once

#include "TallyBase.h"
#include "OpenMCCellAverageProblem.h"

#include "openmc/tallies/filter_cell.h"

class CellTally : public TallyBase
{
public:
  static InputParameters validParams();

  CellTally(const InputParameters & parameters);

  virtual std::pair<unsigned int, openmc::Filter *> spatialFilter() override;

  const std::unordered_set<SubdomainID> & getBlocks() const { return _tally_blocks; }

protected:
  virtual Real storeResultsInner(const std::vector<unsigned int> & var_numbers,
                                 unsigned int local_score,
                                 unsigned int global_score,
                                 std::vector<xt::xtensor<double, 1>> tally_vals) override;

  void checkCellMappedSubdomains();

  std::vector<OpenMCCellAverageProblem::cellInfo> getTallyCells() const;

  std::unordered_set<SubdomainID> _tally_blocks;

  std::map<OpenMCCellAverageProblem::cellInfo, bool> _cell_has_tally;

  openmc::CellInstanceFilter * _cell_filter;

  const bool & _check_equal_mapped_tally_volumes;

  const Real & _equal_tally_volume_abs_tol;
};