PeakReader.cpp

Go to the documentation of this file.

//
// Class PeakReader
//   Implements a parser and value extractor for peak files (*.peaks)
//
// Copyright (c) 2017 - 2018, Matthias Frey, Paul Scherrer Institut, Villigen PSI, Switzerland
// All rights reserved
//
// Implemented as part of the PhD thesis
// "Precise Simulations of Multibunches in High Intensity Cyclotrons"
// and the paper
// "Matching of turn pattern measurements for cyclotrons using multiobjective optimization"
// (https://doi.org/10.1103/PhysRevAccelBeams.22.064602)
//
// This file is part of OPAL.
//
// OPAL is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#include <fstream>
#include <iterator>

#include "Util/PeakReader.h"
#include "Util/OptPilotException.h"

PeakReader::PeakReader(std::string filename)
    : filename_m(filename)
{ }


PeakReader::~PeakReader() { }


void PeakReader::parseFile() {
    
    std::ifstream peak_file;
    
    peak_file.open(filename_m.c_str(), std::ios::in);
    
    if ( !peak_file ) {
        throw OptPilotException("PeakReader::parseFile()",
                                "Error opening file " + filename_m);
    }
    
    // skip header
    std::string header;
    std::getline(peak_file, header);
    
    if ( header.find("# Peak") == std::string::npos ) {
        throw OptPilotException("PeakReader::parseFile()",
                                "Error reading file " + filename_m);
    }
    
    
    int nPeaks = 0;
    peaks_m.clear();
    std::istream_iterator<double> it(peak_file);
    while ( it != std::istream_iterator<double>() ) {
        peaks_m[++nPeaks] = *it;
        ++it;
    }
    
    peak_file.close();
}


void PeakReader::getPeak(int nPeak, double& radius) {
    
    if ( peaks_m.count(nPeak) > 0 ) {
        radius = peaks_m[nPeak];
    } else {
        throw OptPilotException("PeakReader::getPeak",
                                "peak not found!");
    }
}


std::size_t PeakReader::getNumberOfPeaks() {
    return peaks_m.size();
}