spectrum_lightcurve.h

#ifndef SPECTRUM_H
#define SPECTRUM_H

#include <array>
#include <string>
#include <vector>

#include "exspec.h"
#include "packet.h"

struct Spectra {
  double nu_min = -1.;
  double nu_max = -1.;
  std::array<float, MNUBINS> lower_freq;
  std::array<float, MNUBINS> delta_freq;
  std::vector<double> fluxalltimesteps;
  std::vector<double> absorptionalltimesteps;
  std::vector<double> emissionalltimesteps;
  std::vector<double> trueemissionalltimesteps;
  bool do_emission_res = false;
};

void write_spectrum(const std::string &spec_filename, const std::string &emission_filename,
                    const std::string &trueemission_filename, const std::string &absorption_filename,
                    const Spectra &spectra, int numtimesteps);

void write_specpol(const std::string &specpol_filename, const std::string &emission_filename,
                   const std::string &absorption_filename, const Spectra *stokes_i, const Spectra *stokes_q,
                   const Spectra *stokes_u);

void add_to_spec_res(const Packet &pkt, int current_abin, Spectra &spectra, Spectra *stokes_i, Spectra *stokes_q,
                     Spectra *stokes_u);

void init_spectra(Spectra &spectra, double nu_min, double nu_max, bool do_emission_res);
void init_spectrum_trace();
void write_partial_lightcurve_spectra(int my_rank, int nts, const Packet *pkts);

void add_to_lc_res(const Packet &pkt, int current_abin, std::vector<double> &light_curve_lum,
                   std::vector<double> &light_curve_lumcmf);

void write_light_curve(const std::string &lc_filename, int current_abin, const std::vector<double> &light_curve_lum,
                       const std::vector<double> &light_curve_lumcmf, int numtimesteps);

#endif  // SPECTRUM_H