KaonPhysics.h

#ifndef SIMCORE_KAON_PHYSICS_H
#define SIMCORE_KAON_PHYSICS_H
 
#include <G4DecayTable.hh>
#include <G4KaonMinus.hh>
#include <G4KaonPlus.hh>
#include <G4KaonZeroLong.hh>
#include <G4KaonZeroShort.hh>
#include <G4ParticleDefinition.hh>
#include <G4VDecayChannel.hh>
#include <G4VPhysicsConstructor.hh>
#include <iomanip>
#include <iostream>
#include <numeric>
#include <vector>
 
#include "Framework/Configure/Parameters.h"
#include "Framework/Exception/Exception.h"
#include "Framework/Logger.h"
 
namespace simcore {
 
class KaonPhysics : public G4VPhysicsConstructor {
 private:
  enum ChargedKaonDecayChannel {
    mu_nu = 0,
    pi_pi0 = 1,
    pi_pi_pi = 2,
    pi0_e_nu = 3,
    pi0_mu_nu = 4,
    pi_pi0_pi0 = 5,
  };
  enum KaonZeroLongDecayChannel {
    pim_e_nu = 0,
    pip_e_nu = 1,
    pi0_pi0_pi0 = 2,
    pim_mu_nu = 3,
    pip_mu_nu = 4,
    pi0_pip_pim = 5,
  };
  enum KaonZeroShortDecayChannel {
    pip_pim = 0,
    pi0_pi0 = 1,
  };
 
  // Factor to scale the K^+/K^-/K^0_L/K^0_S lifetimes by. To reduce the
  // lifetime of all charged kaons by a factor 50, set these to 1/50.
  double kplus_lifetime_factor{1};
  double kminus_lifetime_factor{1};
  double k0l_lifetime_factor{1};
  double k0s_lifetime_factor{1};
 
  // Branching ratios for each of the decay processes listed in the
  // KaonDecayChannel enumerators
  std::vector<double> kplus_branching_ratios;
  std::vector<double> kminus_branching_ratios;
  std::vector<double> k0l_branching_ratios;
  std::vector<double> k0s_branching_ratios;
 
  // If > 0, dump details about what was changed
  // If > 1, dump details about the initial branching ratios
  int verbosity;
 
 public:
  KaonPhysics(const G4String& name,
              const framework::config::Parameters& parameters);
  virtual ~KaonPhysics() = default;
 
  void setDecayProperties(G4ParticleDefinition* kaon,
                          const std::vector<double>& branching_ratios,
                          double lifetime_factor) const;
 
  void ConstructParticle() override;
 
  void DumpDecayDetails(const G4ParticleDefinition* kaon) const;
 
  void ConstructProcess() override{};
 
 private:
  mutable framework::logging::logger theLog_{
      framework::logging::makeLogger("KaonPhysics")};
};
}  // namespace simcore
 
#endif /* SIMCORE_KAON_PHYSICS_H */