JTransmittances.hh

Go to the documentation of this file.

#ifndef __JCALIBRATE_JTRANSMITTANCES__
#define __JCALIBRATE_JTRANSMITTANCES__
 
#include <istream>
#include <ostream>
#include <iomanip>
#include <array>
#include <map>
 
#include "JLang/JObjectStreamIO.hh"
#include "JLang/JManip.hh"
 
#include "JFitK40.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JCALIBRATE {}
namespace JPP { using namespace JCALIBRATE; }
 
namespace JCALIBRATE {
 
  using JLANG::JObjectStreamIO;
 
  /**
   * Auxiliary data structure for I/O of transmittances of optical module.
   */
  struct transmittance_type :
    public std::array<double, NUMBER_OF_RINGS>
  {
    /**
     * Default constructor.
     */
    transmittance_type()
    {
      for (int i = 0; i != NUMBER_OF_RINGS; ++i) {
        (*this)[i] = 1.0;
      }
    }
 
 
    /**
     * Copy constructor.
     *
     * \param  transmittance     model transmittances
     */
    transmittance_type(const JTransmittance_t& transmittance)
    {
      for (int i = 0; i != NUMBER_OF_RINGS; ++i) {
        (*this)[i] = transmittance[i];
      }
    }
 
 
    /**
     * Read transmittances from input stream.
     *
     * \param  in                input stream
     * \param  transmittance     transmittances
     * \return                   input stream
     */
    friend inline std::istream& operator>>(std::istream& in, transmittance_type& transmittance)
    {
      for (int i = 0; i != NUMBER_OF_RINGS; ++i) {
        in >> transmittance[i];
      }
 
      return in;
    }
 
 
    /**
     * Write transmittances to output stream.
     *
     * \param  out               output stream
     * \param  transmittance     transmittances
     * \return                   output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const transmittance_type& transmittance)
    {
      for (int i = 0; i != NUMBER_OF_RINGS; ++i) {
        out << (i != 0 ? " " : "") << FIXED(7,3) << transmittance[i];
      }
      return out;
    }
  };
 
 
  /**
   * Auxiliary data structure for I/O of transmittances of optical modules in detector.
   */
  struct transmittances_type :
    public std::map<int, transmittance_type>,
    public JObjectStreamIO<transmittances_type>
  {
    using std::map<int, transmittance_type>::operator[];
 
    /**
     * Get transmittances of given optical module.
     *
     * In the absence of a module identifier, default values are returned.
     *
     * \param  id                module identifier
     * \return                   transmittances
     */
    transmittance_type operator[](const int id) const
    {
      const_iterator p = find(id);
 
      if (p != this->end())
        return p->second;
      else
        return transmittance_type();
    }
 
 
    /**
     * Read transmittances from input stream.
     *
     * \param  in                input stream
     * \param  transmittances    transmittances
     * \return                   input stream
     */
    friend inline std::istream& operator>>(std::istream& in, transmittances_type& transmittances)
    {
      int                id;
      transmittance_type transmittance;
 
      while (in >> id >> transmittance) {
        transmittances[id] = transmittance;
      }
 
      return in;
    }
 
 
    /**
     * Write transmittances to output stream.
     *
     * \param  out               output stream
     * \param  transmittances    transmittances
     * \return                   output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const transmittances_type& transmittances)
    {
      using namespace std;
 
      for (const auto& i : transmittances) {
        out << i.first << ' ' << i.second << endl;
      }
 
      return out;
    }
  };
}
 
#endif