tracker.hpp

#pragma once

#include "model.hpp"

namespace IOSKJ {

struct Tracker {

    std::ofstream file;

    Tracker(std::string path){
        file.open(path);
        file
            <<"replicate\t"
            <<"procedure\t"
            <<"year\t"
            <<"quarter\t"

            <<"recruits_determ_we\t"
            <<"recruits_determ_ma\t"
            <<"recruits_determ_ea\t"

            <<"recruits_deviation\t"

            <<"recruits_we\t"
            <<"recruits_ma\t"
            <<"recruits_ea\t"

            <<"biomass_status\t"

            <<"biomass_spawners_we\t"
            <<"biomass_spawners_ma\t"
            <<"biomass_spawners_ea\t"

            <<"biomass_spawning_we\t"
            <<"biomass_spawning_ma\t"
            <<"biomass_spawning_ea\t"

            <<"biomass_vulnerable_we_ps\t"
            <<"biomass_vulnerable_ma_pl\t"
            <<"biomass_vulnerable_ea_gn\t"

            <<"catches_total\t"
            <<"catches_we_ps\t"
            <<"catches_ma_pl\t"
            <<"catches_ea_gn\t"

            <<"effort_total\t"
            <<"effort_we_ps\t"
            <<"effort_ma_pl\t"
            <<"effort_ea_gn\t"

            <<"exp_rate_we_ps\t"
            <<"exp_rate_ma_pl\t"
            <<"exp_rate_ea_gn\t"
            <<std::endl;
    }

    void get(int replicate, int procedure, int time, const Model& model){
        uint year = IOSKJ::year(time);
        uint quarter = IOSKJ::quarter(time);
        file
            <<replicate<<"\t"
            <<procedure<<"\t"
            <<year<<"\t"
            <<quarter<<"\t"

            <<model.recruits_determ(WE)<<"\t"
            <<model.recruits_determ(MA)<<"\t"
            <<model.recruits_determ(EA)<<"\t"

            <<model.recruits_deviation<<"\t"

            <<model.recruits(WE)<<"\t"
            <<model.recruits(MA)<<"\t"
            <<model.recruits(EA)<<"\t"

            <<model.biomass_status()<<"\t"

            <<model.biomass_spawners(WE)<<"\t"
            <<model.biomass_spawners(MA)<<"\t"
            <<model.biomass_spawners(EA)<<"\t"

            <<model.biomass_spawning(WE,quarter)<<"\t"
            <<model.biomass_spawning(MA,quarter)<<"\t"
            <<model.biomass_spawning(EA,quarter)<<"\t"

            <<model.biomass_vulnerable(WE,PS)<<"\t"
            <<model.biomass_vulnerable(MA,PL)<<"\t"
            <<model.biomass_vulnerable(EA,GN)<<"\t"
            
            <<model.catches_taken(sum)<<"\t"
            <<model.catches_taken(WE,PS)<<"\t"
            <<model.catches_taken(MA,PL)<<"\t"
            <<model.catches_taken(EA,GN)<<"\t"

            <<model.effort(sum)<<"\t"
            <<model.effort(WE,PS)<<"\t"
            <<model.effort(MA,PL)<<"\t"
            <<model.effort(EA,GN)<<"\t"

            <<model.exploitation_rate(WE,PS)<<"\t"
            <<model.exploitation_rate(MA,PL)<<"\t"
            <<model.exploitation_rate(EA,GN)<<"\t"

            <<std::endl;
    }
};

} // namespace IOSKJ