population.h

#ifndef _POPULATION_H_
#define _POPULATION_H_

#include <string.h>
#include <iostream>
#include <sstream>
#include "data.h"
#include "sample_name.h"
#include "typedef.h"
#include "genotype.h"

class Population : public Indexed_data {
private:
    std::vector <Genotype>::iterator igl_;      
    std::vector <std::string> sample_names_;    
    static const Registration registered;
    static Data * create(const std::vector <std::string> & Columns){
        return new Population(Columns);
    }
public:
//Data block
    std::vector <Genotype> likelihoods;     

    gt_t major;                 
    gt_t minor;                 
    float_t m;                  
    float_t f;                  
//end Data block
    Population();                   
    Population(const Sample_name &);        
    Population(const std::vector <std::string> &);  
    Population(const Population &);         
    ~Population();                  
    size_t size() const;                    
    void add(const Genotype &likelihood);           
    void add(const float_t &lMM, const float_t &lMm, const float_t &lmm, const count_t &N); 
    void clear();                       
    std::string header(void) const;             

    static const std::string table_name;            
    static const std::string file_name;         
    static const bool binary;
    const bool get_binary() const;

    inline std::vector <std::string> get_sample_names(void) const {return sample_names_;};      
    inline void set_sample_names(const std::vector <std::string>& sample_names) {
        sample_names_=sample_names;
        likelihoods.resize(sample_names.size() );
    };      

    Population & operator= (const Population&);

    void write (std::ostream&) const;
    void read (std::istream&);

    const std::string get_file_name() const {return Population::file_name;};
        const std::string get_table_name() const {return Population::table_name;};

    void write_binary (std::ostream&) const;
    void read_binary (std::istream&);
};

#endif