BinSpecification.cc

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) BDSIM Collaboration, 2001 - 2026.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "BinSpecification.hh"
#include "RBDSException.hh"
 
#include <string>
#include <vector>
 
ClassImp(BinSpecification)
 
BinSpecification::BinSpecification():
  low(0),
  high(1),
  n(1),
  edges(nullptr),
  isLogSpaced(false)
{;}
 
BinSpecification::BinSpecification(double lowIn,
                                   double highIn,
                                   int    nIn):
  low(lowIn),
  high(highIn),
  n(nIn),
  edges(nullptr),
  isLogSpaced(false)
{
  if (high <= low)
    {throw RBDSException("high end of binning <= low end -> must be >");}
  if (n < 1)
    {throw RBDSException("n bins < 1 -> must be >= 1");}
}
 
BinSpecification::BinSpecification(const std::vector<double>& edgesIn):
  low(0),
  high(1),
  n(1),
  edges(new std::vector<double>(edgesIn)),
  isLogSpaced(false)
{
  if (edges->size() < 2)
    {throw RBDSException("too few bin edges -> must be at least 2 edges to define 1 bin");}
  n    = (int)edges->size() - 1;
  low  = (*edges)[0];
  high = edges->back();
}
  
BinSpecification::BinSpecification(const BinSpecification& other):
  low(other.low),
  high(other.high),
  n(other.n),
  edges(nullptr),
  isLogSpaced(other.isLogSpaced),
  edgesFileName(other.edgesFileName)
{
  edges = other.edges ? new std::vector<double>(*other.edges) : nullptr;
}
 
BinSpecification::BinSpecification(BinSpecification&& other) noexcept:
  low(other.low),
  high(other.high),
  n(other.n),
  edges(other.edges),
  isLogSpaced(other.isLogSpaced),
  edgesFileName(other.edgesFileName)
{
  other.low   = 0;
  other.high  = 1;
  other.n     = 1;
  other.edges = nullptr;
  other.isLogSpaced = false;
  other.edgesFileName.clear();
}
 
BinSpecification::~BinSpecification()
{
  delete edges;
}
 
std::string BinSpecification::GetBinString() const
{
  if (!edgesFileName.empty())
    {return edgesFileName;}
  else
    {return std::to_string(low) + ":" + std::to_string(high);}
}