NumberOfPeaks.h

Go to the documentation of this file.

//
// Struct NumberOfPeaks
//   A simple expression to check the number of turns in a circular machine. It checks
//   probe files (*.peaks) and counts the number of turns.
//
// Copyright (c) 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/>.
//
#ifndef __NUMBER_OF_PEAKS_H__
#define __NUMBER_OF_PEAKS_H__
 
#include <cmath>
#include <string>
#include <tuple>
#include <variant>
 
#include "Util/Types.h"
#include "Util/PeakReader.h"
#include "Expression/Parser/function.hpp"
 
struct NumberOfPeaks {
 
    static const std::string name;
 
    Expressions::Result_t operator()(client::function::arguments_t args) {
        if (args.size() != 1) {
            throw OptPilotException("NumberOfPeaks::operator()",
                                    "numberOfPeaks expects 1 arguments, " + std::to_string(args.size()) + " given");
        }
 
        sim_filename_  = std::get<std::string>(args[0]);
 
        bool is_valid = true;
 
        const std::unique_ptr<PeakReader> sim_peaks(new PeakReader(sim_filename_));
        std::size_t nPeaks = 0;
 
        try {
            sim_peaks->parseFile();
            nPeaks = sim_peaks->getNumberOfPeaks();
        } catch (OptPilotException &ex) {
            std::cout << "Caught exception: " << ex.what() << std::endl;
            is_valid = false;
        }
 
 
        return std::make_tuple(nPeaks, is_valid);
    }
 
private:
    std::string sim_filename_;
 
    // define a mapping to arguments in argument vector
    std::tuple<std::string> argument_types;
};
 
#endif