PeakFinder.h

Go to the documentation of this file.

#ifndef PEAKFINDER_H
#define PEAKFINDER_H

#include "Utility/IpplInfo.h"
#include "Algorithms/Vektor.h"

#include <fstream>
#include <string>
#include <vector>
#include <list>

class PeakFinder {

public:
    using container_t = std::vector<double>;

public:

    PeakFinder() = delete;

    PeakFinder(std::string elem, double min, double max, double binwidth, bool singlemode);

    void addParticle(const Vector_t& R);

    void evaluate(const int& turn);

    void save();

private:

    // compute global histogram, involves some inter-node communication
    void createHistogram_m();

    /***************
     * Output file *
     ***************/
    void open_m();
    void append_m();
    void close_m();
    void saveASCII_m();

    void computeCentroid_m();

private:
    container_t radius_m;
    container_t globHist_m;

    std::string fn_m;

    std::string hist_m;

    std::ofstream os_m;

    std::ofstream hos_m;

    std::string element_m;

    // Histogram details
    unsigned int nBins_m;
    double binWidth_m;
    double min_m, max_m;

    int turn_m;
    double peakRadius_m;
    int registered_m;
    std::list<double> peaks_m;
    bool singlemode_m;
    bool first_m;
    bool finished_m;
};

#endif