scatteractionsfinderparameters.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2022-2024
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 */
#ifndef SRC_INCLUDE_SMASH_SCATTERACTIONSFINDERPARAMETERS_H_
#define SRC_INCLUDE_SMASH_SCATTERACTIONSFINDERPARAMETERS_H_
 
#include <memory>
#include <set>
#include <utility>
 
#include "input_keys.h"
 
namespace smash {
 
struct StringTransitionParameters {
  const std::pair<double, double> sqrts_range_Npi =
      InputKeys::collTerm_stringTrans_rangeNpi.default_value();
  const std::pair<double, double> sqrts_range_NN =
      InputKeys::collTerm_stringTrans_rangeNN.default_value();
  const double sqrts_add_lower =
      InputKeys::collTerm_stringTrans_lower.default_value();
  const double sqrts_range_width =
      InputKeys::collTerm_stringTrans_range_width.default_value();
  const double pipi_offset =
      InputKeys::collTerm_stringTrans_pipiOffset.default_value();
  const double KN_offset =
      InputKeys::collTerm_stringTrans_KNOffset.default_value();
};
 
class ScatterActionsFinderParameters {
 public:
  ScatterActionsFinderParameters(Configuration& config,
                                 const ExperimentParameters& parameters);
  const double elastic_parameter;
  const double low_snn_cut;
  const double scale_xs;
  const double additional_el_xs;  // mb
  const double maximum_cross_section;
  const CollisionCriterion coll_crit;
  const NNbarTreatment nnbar_treatment;
  const ReactionsBitSet included_2to2;
  const MultiParticleReactionsBitSet included_multi;
  const int testparticles;
  const bool two_to_one;
  const bool allow_collisions_within_nucleus;
  const bool strings_switch;
  const bool use_AQM;
  const bool strings_with_probability;
  const bool only_warn_for_high_prob;
  const StringTransitionParameters transition_high_energy;
  const TotalCrossSectionStrategy total_xs_strategy;
  const PseudoResonance pseudoresonance_method;
 
  double AQM_scaling_factor(const PdgCode& pdg) const {
    return (1 - AQM_strange_suppression * pdg.frac_strange()) *
           (1 - AQM_charm_suppression * pdg.frac_charm()) *
           (1 - AQM_bottom_suppression * pdg.frac_bottom());
  }
 
 private:
  const double AQM_strange_suppression = 0.4;
  const double AQM_charm_suppression;
  const double AQM_bottom_suppression;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_SCATTERACTIONSFINDERPARAMETERS_H_