RMGPhysics.hh

// Copyright (C) 2022 Luigi Pertoldi <https://orcid.org/0000-0002-0467-2571>
//
// This program is free software: you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License as published by the Free
// Software Foundation, either version 3 of the License, or (at your option) 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 Lesser General Public License for more
// details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with this program.  If not, see <https://www.gnu.org/licenses/>.
 
#ifndef _RMG_PHYSICS_HH_
#define _RMG_PHYSICS_HH_
 
#include <map>
#include <memory>
 
#include "G4GenericMessenger.hh"
#include "G4VModularPhysicsList.hh"
#include "globals.hh"

class RMGPhysics : public G4VModularPhysicsList {
 
  public:

    RMGPhysics();
 
    RMGPhysics(RMGPhysics const&) = delete;
    RMGPhysics& operator=(RMGPhysics const&) = delete;
    RMGPhysics(RMGPhysics&&) = delete;
    RMGPhysics& operator=(RMGPhysics&&) = delete;

    enum class LowEnergyEMOption {
      kOption1,
      kOption2,
      kOption3,
      kOption4,
      kPenelope,
      kLivermore,
      kLivermorePolarized,
      kNone
    };

    enum class HadronicPhysicsListOption {
      kQGSP_BIC_HP,
      kQGSP_BERT_HP,
      kFTFP_BERT_HP,
      kShielding,
      kNone
    };
 

    struct ProdCutStore {
        ProdCutStore() = default;

        ProdCutStore(double def_cut)
            : gamma(def_cut), electron(def_cut), positron(def_cut), proton(def_cut), alpha(def_cut),
              generic_ion(def_cut) {}
 
        double gamma;
        double electron;
        double positron;
        double proton;
        double alpha;
        double generic_ion;
    };

    void SetCuts() override;

    void SetLowEnergyRange(G4double low_energy) { fLowEnergyRange = low_energy; };

    void SetHighEnergyRange(G4double high_energy) { fHighEnergyRange = high_energy; };

    void SetDefaultProductionCut(double cut);

    void SetSensitiveProductionCut(double cut);

    void SetLowEnergyEMOptionString(std::string option);

    void SetHadronicPhysicsListOptionString(std::string option);

    void SetUseGammaAngCorr(bool);
 
    void SetGammaTwoJMAX(int max_two_j);
    void SetStoreICLevelData(bool);
 
    void DumpProcessesForParticles(std::string file_name);
 
  protected:
 
    void ConstructParticle() override;
    void ConstructProcess() override;
    virtual void ConstructOptical();
 
  private:
 
    ProdCutStore fProdCuts = {};
    ProdCutStore fProdCutsSensitive = {};
    bool fConstructOptical = false;
    bool fUseOpticalCustomWLS = true;
    bool fUseNeutronThermalScattering = false;
    bool fUseGrabmayrGammaCascades = false;
    bool fUseInnerBremsstrahlung = false;
    LowEnergyEMOption fLowEnergyEMOption = LowEnergyEMOption::kLivermore;
    HadronicPhysicsListOption fHadronicPhysicsListOption = HadronicPhysicsListOption::kNone;
    G4double fLowEnergyRange = 250 * CLHEP::eV;
    G4double fHighEnergyRange = 100. * CLHEP::GeV;
    std::unique_ptr<G4GenericMessenger> fMessenger;
    void DefineCommands();
};
 
#endif
 
// vim: tabstop=2 shiftwidth=2 expandtab