/scratch0/jsnuveri/BDSIM/BDSIMgit/bdsim/include/BDSBeamline.hh

#ifndef BDSBEAMLINE_H
#define BDSBEAMLINE_H

#include "globals.hh" // geant4 globals / types

#include "BDSBeamlineElement.hh"

#include <iterator>
#include <ostream>
#include <utility>    //for std::pair
#include <vector>


class BDSBeamline;

class BDSLine;
class BDSTransform3D;

typedef std::vector<BDSBeamlineElement*>::const_iterator BDSBeamlineIterator;

class BDSBeamline{
public:
  BDSBeamline();

  BDSBeamline(G4ThreeVector     initialGlobalPosition,
              G4RotationMatrix* initialGlobalRotation);
  
  ~BDSBeamline();

  void AddSingleComponent(BDSAcceleratorComponent* component);

  void AddComponent(BDSAcceleratorComponent* component);

  void ApplyTransform3D(BDSTransform3D* component);

  void PrintAllComponents(std::ostream& out) const;
  
  BDSBeamlineElement* GetFirstItem();             
  BDSBeamlineElement* GetLastItem();              

  inline G4double     GetTotalChordLength() const;

  inline G4double     GetTotalArcLength() const;

  std::vector<BDSBeamlineElement*>::size_type size() const;

  G4ThreeVector GetMaximumExtentPositive() const;

  G4ThreeVector GetMaximumExtentNegative() const;

  G4ThreeVector GetMaximumExtentAbsolute() const;

  // Accessors in a similar style to std::vector
  BDSBeamlineElement* front() const;
  BDSBeamlineElement* back()  const;
  inline std::vector<BDSBeamlineElement*>::iterator begin();
  inline std::vector<BDSBeamlineElement*>::iterator end();
  inline std::vector<BDSBeamlineElement*>::reverse_iterator rbegin();
  inline std::vector<BDSBeamlineElement*>::reverse_iterator rend();
  inline std::vector<BDSBeamlineElement*>::const_iterator begin() const;
  inline std::vector<BDSBeamlineElement*>::const_iterator end()   const;
  inline std::vector<BDSBeamlineElement*>::const_reverse_iterator rbegin() const;
  inline std::vector<BDSBeamlineElement*>::const_reverse_iterator rend()   const;
  inline G4bool empty() const;

  friend std::ostream& operator<< (std::ostream &out, BDSBeamline const &bl);
  
private: 
  std::vector<BDSBeamlineElement*> beamline; 

  G4double totalChordLength;
  G4double totalArcLength;

  G4ThreeVector maximumExtentPositive; 
  G4ThreeVector maximumExtentNegative; 

  G4ThreeVector xARS, yARS, zARS, xARM, yARM, zARM, xARE, yARE, zARE;

  G4RotationMatrix* previousReferenceRotationEnd;

  G4ThreeVector previousReferencePositionEnd;

  G4double previousSPositionEnd;
  
  BDSBeamline& operator=(const BDSBeamline&);
  BDSBeamline(BDSBeamline&);
};

inline BDSBeamlineElement* BDSBeamline::GetFirstItem()
{return front();}

inline BDSBeamlineElement* BDSBeamline::GetLastItem()
{return back();}

inline G4double BDSBeamline::GetTotalChordLength() const
{return totalChordLength;}

inline G4double BDSBeamline::GetTotalArcLength() const
{return totalArcLength;}

inline std::vector<BDSBeamlineElement*>::size_type BDSBeamline::size() const
{return beamline.size();}

inline G4ThreeVector BDSBeamline::GetMaximumExtentPositive() const
{return maximumExtentPositive;}

inline G4ThreeVector BDSBeamline::GetMaximumExtentNegative() const
{return maximumExtentNegative;}

inline BDSBeamlineElement* BDSBeamline::front() const
{return beamline.front();}

inline BDSBeamlineElement* BDSBeamline::back() const
{return beamline.back();}

inline std::vector<BDSBeamlineElement*>::iterator BDSBeamline::begin()
{return beamline.begin();}

inline std::vector<BDSBeamlineElement*>::iterator BDSBeamline::end()
{return beamline.end();}

inline std::vector<BDSBeamlineElement*>::reverse_iterator BDSBeamline::rbegin()
{return beamline.rbegin();}

inline std::vector<BDSBeamlineElement*>::reverse_iterator BDSBeamline::rend()
{return beamline.rend();}

inline std::vector<BDSBeamlineElement*>::const_iterator BDSBeamline::begin() const
{return beamline.begin();}

inline std::vector<BDSBeamlineElement*>::const_iterator BDSBeamline::end() const
{return beamline.end();}

inline std::vector<BDSBeamlineElement*>::const_reverse_iterator BDSBeamline::rbegin() const
{return beamline.rbegin();}

inline std::vector<BDSBeamlineElement*>::const_reverse_iterator BDSBeamline::rend() const
{return beamline.rend();}

inline G4bool BDSBeamline::empty() const
{return beamline.empty();}

#endif

---