Target.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
 
// Standard library includes
#include <ostream>
#include <map>
#include <vector>
 
// MARLEY includes
#include "TargetAtom.hh"
 
namespace marley {

  class Target {
 
    public:

      explicit Target( int pdg );

      explicit Target( int Z, int A );

      Target( const std::vector<TargetAtom>& nuclides,
        const std::vector<double>& atom_fracs );

      double atom_fraction( const marley::TargetAtom& atom ) const;

      bool contains( const marley::TargetAtom& atom ) const;

      inline bool has_single_nuclide() const
        { return (atom_fractions_.size() == 1u); }

      void print( std::ostream& out ) const;

      inline const std::map< marley::TargetAtom, double >& atom_fraction_map()
        const { return atom_fractions_; }
 
    protected:

      void initialize_single_nuclide( int pdg );

      void initialize( const std::vector<TargetAtom>& nuclides,
        const std::vector<double>& atom_fracs );

      std::map< marley::TargetAtom, double > atom_fractions_;
  };
 
}
 
inline std::ostream& operator<<(std::ostream& out, const marley::Target& t) {
  t.print( out );
  return out;
}