NuclearReaction.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
#include <functional>
#include <map>
#include <memory>
#include <regex>
#include <string>
#include <vector>
 
#include "HepMC3/Attribute.h"
#include "HepMC3/GenEvent.h"
#include "HepMC3/GenParticle.h"
#include "marley/DecayScheme.hh"
#include "marley/Level.hh"
#include "marley/MassTable.hh"
#include "marley/MatrixElement.hh"
#include "marley/Reaction.hh"
#include "marley/StructureDatabase.hh"
 
namespace marley {
 
  class Generator;

  class NuclearReaction : public Reaction {
 
    public:

      NuclearReaction( ProcessType pt, int pdg_a, int pdg_b, int pdg_c,
        int pdg_d, int q_d, const std::string& source_file );
 
      inline virtual marley::TargetAtom atomic_target() const override final
        { return marley::TargetAtom( pdg_b_ ); }

      double threshold_kinetic_energy() const override;

      double max_level_energy( double KEa ) const;

      double weak_nuclear_charge() const;
 
    protected:

      virtual void set_description();

      void set_charge_attributes( std::shared_ptr< HepMC3::GenEvent >& event )
        const;

      void set_nuclear_residue_attributes(
        const std::shared_ptr< HepMC3::GenParticle >& residue,
        double E_level, int twoJ, const marley::Parity& P ) const;

      virtual std::shared_ptr< HepMC3::GenEvent > make_nuclear_event_object(
        double KEa, double pc_cm, double cos_theta_c_cm, double phi_c_cm,
        double Ec_cm, double Ed_cm, double E_level, int twoJ,
        const marley::Parity& P ) const;

      virtual std::shared_ptr< HepMC3::GenEvent > make_nuclear_event_object(
        double KEa, const std::shared_ptr< HepMC3::GenParticle >& ejectile,
        const std::shared_ptr< HepMC3::GenParticle >& residue,
        double E_level, int twoJ, const marley::Parity& P ) const;
 
      double md_gs_; 
 
      int Zi_; 
      int Ai_; 
      int Zf_; 
      int Af_; 

      int q_d_;

      double KEa_threshold_;
  };
 
}