NeutronVector.hh

/*
$Id: NeutronVector.hh 2332 2012-02-16 01:05:03Z jisuzuki $
*/


#ifndef NEUTRONVECTOR_HH
#define NEUTRONVECTOR_HH
#include <boost/serialization/serialization.hpp>

#include "Header.hh"


template < class T , class H >
class NeutronVector
{
private:
  void erase();

  friend class boost::serialization::access;
  template <class Archive>
  void serialize(Archive &ar, const unsigned int ver);

  friend class WriteNeXusFile;
  template <class Write>
  void NXwrite(Write &W) const;

  friend class ReadNeXusFile;
  template <class Read>
  void NXread(Read &R);

protected:
  template <class Archive>
  std::vector<std::string> presave(Archive &ar, const std::string &masterfile, const unsigned int splitnum) const;

  template <class Archive>
  void preload(Archive &ar,
                  std::vector<std::string> &S, std::vector<UInt4> &csize);

  template <class Archive>
  friend class WriteSerializationFile;

  template <class Archive>
  friend class ReadSerializationFile;

protected:
  vector< T* > v;
  H *header;

public:
  NeutronVector();
  NeutronVector( H pheader );
  ~NeutronVector(); 
  NeutronVector( const NeutronVector &ob );

  NeutronVector operator=( const NeutronVector &ob );

  void Allocate( UInt4 size );
  void Set( UInt4 index, T* value );
  

  UInt4 AddPointer( T* value );
  UInt4 Add( T value );
  UInt4 Add( NeutronVector<T,H>  &value );


  T  Put( UInt4 index=0 );
  T* PutPointer( UInt4 index=0 );

  T* operator()( UInt4 index ) { return PutPointer( index ); }

  T &ref(UInt4 index) { return *v[index]; }

  // Put(), PutPointer(), operator() are to be replaced by this method
  const T &ref(UInt4 index) const { return *v[index]; }

  void InputHeader( H pheader );
  H PutHeader();
  H *PutHeaderPointer() const { return header; }
  void reset(){ erase(); }
  UInt4 PutTableSize() const { return v.size(); }
  UInt4 EraseElement( UInt4 i );
  UInt4 PutSize() const { return PutTableSize(); }
  UInt4 ClearElement( UInt4 i );
  UInt4 copy( UInt4 i );
  void AddToHeader( string Key, Int4   value ){ header->Add( Key, value ); }
  void AddToHeader( string Key, Double value ){ header->Add( Key, value ); }
  void AddToHeader( string Key, string value ){ header->Add( Key, value ); }
  void AddToHeader( string Key, vector<Int4>   value ){ header->Add( Key, value ); }
  void AddToHeader( string Key, vector<Double> value ){ header->Add( Key, value ); }
  void AddToHeader( string Key, vector<string> value ){ header->Add( Key, value ); }
  template <class Func>
  void Transform(const std::string &Key, std::vector< Func > F, const std::string &newunit="None");

  void Reverse(void) ;

  void Resize(UInt4 newsize);

  NeutronVector<T, H> & operator+=(const NeutronVector<T, H> &r );
  NeutronVector<T, H> & operator-=(const NeutronVector<T, H> &r );
  NeutronVector<T, H> & operator*=(const NeutronVector<T, H> &r );
  NeutronVector<T, H> & operator/=(const NeutronVector<T, H> &r );

  NeutronVector<T, H> operator+( NeutronVector<T, H> &r );
  NeutronVector<T, H> operator-( NeutronVector<T, H> &r );
  NeutronVector<T, H> operator*( NeutronVector<T, H> &r );
  NeutronVector<T, H> operator/( NeutronVector<T, H> &r );

  
};

#include "NeutronVector.cc"

#endif