EventExporter.cpp

Go to the documentation of this file.

/*
 *  CepGen: a central exclusive processes event generator
 *  Copyright (C) 2016-2024  Laurent Forthomme
 *
 *  This program 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
 *  any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#include <iomanip>
#include <sstream>
 
#include "CepGen/Core/ParametersList.h"
#include "CepGen/Core/RunParameters.h"
#include "CepGen/EventFilter/EventExporter.h"
#include "CepGen/EventFilter/EventModifier.h"
#include "CepGen/Utils/String.h"
#include "CepGen/Version.h"
 
namespace cepgen {
  EventExporter::EventExporter(const ParametersList& params) : EventHandler(params) {}
 
  std::string EventExporter::banner(const std::string& prep) const {
    const size_t len = 45 + version::tag.size();
    std::ostringstream os;
    os << prep << "******* Sample generated with CepGen " << version::tag << " *******\n";
    if (runParameters().hasProcess()) {
      os << prep << " Process: " << runParameters().processName() << " ("
         << runParameters().kinematics().incomingBeams().mode() << ")\n";
      if (runParameters().kinematics().incomingBeams().mode() != mode::Kinematics::ElasticElastic)
        os << prep << " Structure functions: " << runParameters().kinematics().incomingBeams().structureFunctions()
           << "\n";
    }
    if (!runParameters().eventModifiersSequence().empty()) {
      os << prep << " " << utils::s("Event modifier", runParameters().eventModifiersSequence().size()) << ": ";
      std::string sep;
      for (const auto& mod : runParameters().eventModifiersSequence())
        os << sep << mod->name(), sep = ", ";
      os << "\n";
    }
    if (runParameters().hasProcess()) {
      const auto& cuts = runParameters().kinematics().cuts();
      auto dump_cuts = [&os, &prep](const auto& obj) {
        for (const auto& lim : obj.parameters().template keysOf<Limits>()) {
          const auto& limit = obj.parameters().template get<Limits>(lim);
          if (limit.valid() && obj.description().has(lim))
            os << prep << " " << obj.description().get(lim).description() << ": " << limit << "\n";
        }
      };
      os << prep << std::left << std::setw(len) << std::setfill('*') << "*** Incoming state "
         << "\n";
      dump_cuts(cuts.initial);
      os << prep << std::setw(len) << std::setfill('*') << "*** Central system "
         << "\n";
      dump_cuts(cuts.central);
      if (runParameters().kinematics().incomingBeams().mode() != mode::Kinematics::ElasticElastic) {
        os << prep << std::setw(len) << std::setfill('*') << "*** Remnants states "
           << "\n";
        dump_cuts(cuts.remnants);
      }
    }
    os << prep << std::string(45 + version::tag.size(), '*');
    return os.str();
  }
}  // namespace cepgen