source: svn/trunk/Utilities/Hector/include/H_BeamLine.h@ 12

#ifndef _H_BeamLine_
#define _H_BeamLine_

/*
---- Hector the simulator ----
   A fast simulator of particles through generic beamlines.
   J. de Favereau, X. Rouby ~~~ hector_devel@cp3.phys.ucl.ac.be

        http://www.fynu.ucl.ac.be/hector.html

   Centre de Physique des Particules et de Phénoménologie (CP3)
   Université Catholique de Louvain (UCL)
*/

/// \file H_BeamLine.h
/// \brief Class aiming at retrieving features of the real beam optical elements.

// local #includes
#include "H_AbstractBeamLine.h"

// C++ #includes
#include <string>
using namespace std;

// Caution : mad conventions for vertically(horizontally) focusing quadrupoles
// are opposite to Wille's. See fill() for more details.

/// Reads external files and retrieves features of the real beam optical elements.
class H_BeamLine : public H_AbstractBeamLine {

        public:
        ///     Constructors and destructor
        //@{
                H_BeamLine():H_AbstractBeamLine() {direction=1; ips=0; };
                H_BeamLine(const H_BeamLine& );
                H_BeamLine(const int, const float);
                H_BeamLine& operator=(const H_BeamLine& );
                ~H_BeamLine() {return;};
        //@
                ///     Finds the IP position (s) from the MAD table. Should be "IP5" or "IP1".
                //@{
                void findIP(const string);
                void findIP(const string, const string);
                //@}
                ///     Reader for the external MAD table
                //@{
                void fill(const string);
                void fill(const string, const int, const string );
                //@}
                ///     Returns the IP position (s)
                double getIP() {return ips;};
                ///             Returns positions and angles of beam at IP
                double* getIPProperties();

        private:
                int direction; // to or from the IP.
                double ips; // s-position of the IP [m]
                double ipx; // x-position of the IP [µm]
                double ipy; // y-position of the IP [µm]
                double iptx; // x-angle of momentum at IP [µrad]
                double ipty; // y-angle of momentum at IP [µrad]
};

#endif