JMechanics_t.hh

Go to the documentation of this file.

#ifndef __JACOUSTICS__JMECHANICS_T__
#define __JACOUSTICS__JMECHANICS_T__
 
#include <ostream>
#include <iomanip>
#include <map>
#include <cmath>
 
#include <TROOT.h>
#include <TObject.h>
 
#include "JLang/JException.hh"
#include "JLang/JManip.hh"
 
 
/**
 * \file
 *
 * Mechanical modelling of string.
 * \author mdejong
 */
namespace JACOUSTICS {}
namespace JPP { using namespace JACOUSTICS; }
 
namespace JACOUSTICS {
 
  /**
   * Auxiliary data structure for parameters of mechanical model.\n
   * This data structure provides for the implementation of the effective height conform the mechanical model of string.
   */
  struct JMechanics {
    /**
     * Default constructor.
     */
    JMechanics() :
      a(0.0),
      b(0.0)
    {}
 
 
    /**
     * Constructor.
     *
     * \param  a                        logarithmic term
     * \param  b                        linear      term
     */
    JMechanics(const double a,
               const double b) :
      a(a),
      b(b)
    {}
 
 
    /**
     * Get effective height for given actual height.
     *
     * \param  height                   height
     * \return                          height
     */
    double getHeight(const double height) const
    {
      return height  +  this->b * log(1.0 - this->a * height);
    }
 
 
    /**
     * Read parameters from input stream.
     *
     * \param  in                       input stream
     * \param  parameters               parameters
     * \return                          input stream
     */
    friend inline std::istream& operator>>(std::istream& in, JMechanics& parameters)
    {
      return in >> parameters.a >> parameters.b;
    }
 
 
    /**
     * Write parameters to output stream.
     *
     * \param  out                      output stream
     * \param  parameters               parameters
     * \return                          output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const JMechanics& parameters)
    {
      return out << FIXED(7,5) << parameters.a << ' '
                 << FIXED(7,3) << parameters.b;
    }
 
    ClassDefNV(JMechanics, 1);
 
    double a;   ///< <tt>0 <= a < (maximal height)⁻1;</tt>    [m^-1]
    double b;   ///< <tt>0 <= b;</tt>                         [m]
  };
 
 
  /**
   * Auxiliary data structure for mechanical model parameters of strings in a given detector.
   *
   * Note that the JDetectorMechanics::WILDCARD acts as default value for the string number.
   */
  struct JDetectorMechanics_t :
    public std::map<int, JMechanics>,
    public TObject
  {
    enum {
      WILDCARD = -1   //!< wild card for string number.
    };
 
 
    /**
     * Get mechanical parameters for given string.
     *
     * \param  string                  string number
     * \return                         mechanical parameters
     */
    const JMechanics& operator()(const int string) const
    {
      static const JMechanics mechanics;
 
      const_iterator p;
 
      if      ((p = this->find(string))   != this->end())
        return p->second;
      else if ((p = this->find(WILDCARD)) != this->end())
        return p->second;
      else
        return mechanics;
    }
 
 
    /**
     * Write detector mechanics to output.
     *
     * \param  out          output stream
     * \param  object       detector mechanics
     * \return              output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const JDetectorMechanics_t& object)
    {
      using namespace std;
      
      for (JDetectorMechanics_t::const_iterator i = object.begin(); i != object.end(); ++i) {
        out << setw(4) << i->first << ' ' << i->second << endl;
      }
 
      return out;
    }
 
    ClassDef(JDetectorMechanics_t, 1);
  };
}
 
#endif