simcore::KaonPhysics Class Reference

Allows the configuration of properties of kaons produced in the simulation, in particular setting their lifetime and branching ratios. More...

#include <KaonPhysics.h>

Public Member Functions
	KaonPhysics (const G4String &name, const framework::config::Parameters &parameters)
void	setDecayProperties (G4ParticleDefinition *kaon, const std::vector< double > &branching_ratios, double lifetime_factor) const
	Set the lifetime and branching ratios for one of the kaon species.
void	ConstructParticle () override
	Construct/Update particles.
void	dumpDecayDetails (const G4ParticleDefinition *kaon) const
void	ConstructProcess () override
	Construct processes.

Private Types
enum	ChargedKaonDecayChannel {   mu_nu = 0 , pi_pi0 = 1 , pi_pi_pi = 2 , pi0_e_nu = 3 ,   pi0_mu_nu = 4 , pi_pi0_pi0 = 5 }
	Represents the 6 possible decay channels for charged kaons in Geant4 (See G4KaonMinus.cc and G4KaonPlus.cc in the Geant4 sources) in the order they appear in the decay table. More...
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 }
	Corresponding entries for the K^0_L. More...
enum	KaonZeroShortDecayChannel { pip_pim = 0 , pi0_pi0 = 1 }

Private Attributes
double	kplus_lifetime_factor_ {1}
double	kminus_lifetime_factor_ {1}
double	k0l_lifetime_factor_ {1}
double	k0s_lifetime_factor_ {1}
std::vector< double >	kplus_branching_ratios_
std::vector< double >	kminus_branching_ratios_
std::vector< double >	k0l_branching_ratios_
std::vector< double >	k0s_branching_ratios_
framework::logging::logger	the_log_

Detailed Description

Allows the configuration of properties of kaons produced in the simulation, in particular setting their lifetime and branching ratios.

Definition at line 23 of file KaonPhysics.h.

Member Enumeration Documentation

ChargedKaonDecayChannel

enum simcore::KaonPhysics::ChargedKaonDecayChannel

private

Represents the 6 possible decay channels for charged kaons in Geant4 (See G4KaonMinus.cc and G4KaonPlus.cc in the Geant4 sources) in the order they appear in the decay table.

The processes are K^+ -> \mu^+ + \nu_\mu K^+ -> \pi^+ + \pi^0 K^+ -> \pi^+ + \pi^- + \pi^+ K^+ -> \pi^0 + e^+ + \nu_e K^+ -> \pi^0 + \mu^+ + \nu_\mu K^+ -> \pi^+ + \pi^0 + \pi^0

And vice versa for K^-. The indices here correspond to the position of the branching ratio for that process in the corresponding parameter as well as the position in the decay table.

Note

: The order in the the decay table is sorted by the branching ratios of the default physics settings since G4DecayTable::insert sorts its entries by branching ratio!

Definition at line 50 of file KaonPhysics.h.

                               {
    mu_nu = 0,
    pi_pi0 = 1,
    pi_pi_pi = 2,
    pi0_e_nu = 3,
    pi0_mu_nu = 4,
    pi_pi0_pi0 = 5,
  };

KaonZeroLongDecayChannel

enum simcore::KaonPhysics::KaonZeroLongDecayChannel

private

Corresponding entries for the K^0_L.

Note that K^0_L and K^0_S decays are not symmetric like the charged ones so they need to be handled manually.

The processes are

K^0_L -> \pi^- + e^+ + \nu_e K^0_L -> \pi^+ + e^- + \nu_e K^0_L -> \pi^0 + \pi^0 + \pi^0 K^0_L -> \pi^- + \mu^+ + \nu_\mu K^0_L -> \pi^+ + \mu^- + \nu_\mu K^0_L -> \pi^0 + \pi^+ + \pi^-

and

K^0_S -> \pi^+ + \pi^- K^0_S -> \pi^0 + \pi^0

Note

: The order in the the decay table is sorted by the branching ratios of the default physics settings since G4DecayTable::insert sorts its entries by branching ratio!

Definition at line 82 of file KaonPhysics.h.

                                {
    pim_e_nu = 0,
    pip_e_nu = 1,
    pi0_pi0_pi0 = 2,
    pim_mu_nu = 3,
    pip_mu_nu = 4,
    pi0_pip_pim = 5,
  };

KaonZeroShortDecayChannel

enum simcore::KaonPhysics::KaonZeroShortDecayChannel

private

Definition at line 90 of file KaonPhysics.h.

                                 {
    pip_pim = 0,
    pi0_pi0 = 1,
  };

Constructor & Destructor Documentation

KaonPhysics()

simcore::KaonPhysics::KaonPhysics	(	const G4String &	name,
		const framework::config::Parameters &	parameters )

Definition at line 14 of file KaonPhysics.cxx.

    : G4VPhysicsConstructor(name) {
  kplus_branching_ratios_ =
      parameters.get<std::vector<double>>("kplus_branching_ratios");
  kminus_branching_ratios_ =
      parameters.get<std::vector<double>>("kminus_branching_ratios");
  k0l_branching_ratios_ =
      parameters.get<std::vector<double>>("k0l_branching_ratios");
  k0s_branching_ratios_ =
      parameters.get<std::vector<double>>("k0s_branching_ratios");
  kplus_lifetime_factor_ = parameters.get<double>("kplus_lifetime_factor");
  kminus_lifetime_factor_ = parameters.get<double>("kminus_lifetime_factor");
  k0l_lifetime_factor_ = parameters.get<double>("k0l_lifetime_factor");
  k0s_lifetime_factor_ = parameters.get<double>("k0s_lifetime_factor");
}

Member Function Documentation

ConstructParticle()

void simcore::KaonPhysics::ConstructParticle ( )

override

Construct/Update particles.

Update the particle definitions for charged kaons

1. Scale their lifetime by the corresponding lifetime factor
2. Set their branching ratios to those in the corresponding branching raito parameter

Definition at line 79 of file KaonPhysics.cxx.

                                    {
  auto kaon_plus{G4KaonPlus::Definition()};
  auto kaon_minus{G4KaonMinus::Definition()};
  auto kaon_long{G4KaonZeroLong::Definition()};
  auto kaon_short{G4KaonZeroShort::Definition()};
 
  if (!kaon_plus || !kaon_minus || !kaon_long || !kaon_short) {
    EXCEPTION_RAISE("KaonPhysics",
                    "Unable to get the charged kaon particle definitions, "
                    "something is very wrong with the configuration.");
  }
  setDecayProperties(kaon_plus, kplus_branching_ratios_,
                     kplus_lifetime_factor_);
  setDecayProperties(kaon_minus, kminus_branching_ratios_,
                     kminus_lifetime_factor_);
  setDecayProperties(kaon_long, k0l_branching_ratios_, k0l_lifetime_factor_);
  setDecayProperties(kaon_short, k0s_branching_ratios_, k0s_lifetime_factor_);
}

ConstructProcess()

void simcore::KaonPhysics::ConstructProcess ( )

inlineoverride

Construct processes.

We don't do anything here since we are just attaching/updating the kaon particle definitions.

Definition at line 142 of file KaonPhysics.h.

{};

dumpDecayDetails()

void simcore::KaonPhysics::dumpDecayDetails

(

const G4ParticleDefinition *

kaon

)

const

Definition at line 98 of file KaonPhysics.cxx.

                                                                         {
  ldmx_log(trace) << "Decay table details for " << kaon->GetParticleName()
                  << std::scientific << std::setprecision(15)
                  << " (PDG Lifetime " << kaon->GetPDGLifeTime() << ")"
                  << std::endl;
  auto* table{kaon->GetDecayTable()};
  if (not table) {
    ldmx_log(error) << "No Kaon decay table.";
    return;
  }
  const int entries{table->entries()};
  for (auto i{0}; i < entries; ++i) {
    const auto channel{(*table)[i]};
    const auto daughters{channel->GetNumberOfDaughters()};
    std::string products{};
    // N-1 to avoid extra " + "
    for (auto j{0}; j < daughters; ++j) {
      auto* d{channel->GetDaughter(j)};
      if (d) {
        products += d->GetParticleName();
      } else {
        products += "NULL";
      }
      if (j < daughters - 1) {
        // only add '+' when there are more to come
        products += " + ";
      }
    }
    ldmx_log(trace) << "Channel " << i << " (" << kaon->GetParticleName()
                    << " -> " << products << ") Kinematics type "
                    << channel->GetKinematicsName() << " with BR "
                    << channel->GetBR() << std::endl;
  }
}

setDecayProperties()

void simcore::KaonPhysics::setDecayProperties	(	G4ParticleDefinition *	kaon,
		const std::vector< double > &	branching_ratios,
		double	lifetime_factor ) const

Set the lifetime and branching ratios for one of the kaon species.

Definition at line 30 of file KaonPhysics.cxx.

                                  {
  auto table{kaon->GetDecayTable()};
  if (!table) {
    EXCEPTION_RAISE("KaonPhysics", "Unable to get the decay table from " +
                                       kaon->GetParticleName());
  }
  ldmx_log(trace) << "Decay details (" << kaon->GetParticleName()
                  << ") before setting branching ratios and lifetimes";
  dumpDecayDetails(kaon);
  kaon->SetPDGLifeTime(kaon->GetPDGLifeTime() * lifetime_factor);
  if (kaon == G4KaonZeroLong::Definition()) {
    (*table)[KaonZeroLongDecayChannel::pi0_pi0_pi0]->SetBR(
        branching_ratios[KaonZeroLongDecayChannel::pi0_pi0_pi0]);
    (*table)[KaonZeroLongDecayChannel::pi0_pip_pim]->SetBR(
        branching_ratios[KaonZeroLongDecayChannel::pi0_pip_pim]);
    (*table)[KaonZeroLongDecayChannel::pip_e_nu]->SetBR(
        branching_ratios[KaonZeroLongDecayChannel::pip_e_nu]);
    (*table)[KaonZeroLongDecayChannel::pim_e_nu]->SetBR(
        branching_ratios[KaonZeroLongDecayChannel::pim_e_nu]);
    (*table)[KaonZeroLongDecayChannel::pim_mu_nu]->SetBR(
        branching_ratios[KaonZeroLongDecayChannel::pim_mu_nu]);
    (*table)[KaonZeroLongDecayChannel::pip_mu_nu]->SetBR(
        branching_ratios[KaonZeroLongDecayChannel::pip_mu_nu]);
  } else if (kaon == G4KaonZeroShort::Definition()) {
    (*table)[KaonZeroShortDecayChannel::pip_pim]->SetBR(
        branching_ratios[KaonZeroShortDecayChannel::pip_pim]);
    (*table)[KaonZeroShortDecayChannel::pi0_pi0]->SetBR(
        branching_ratios[KaonZeroShortDecayChannel::pi0_pi0]);
  } else {
    (*table)[ChargedKaonDecayChannel::mu_nu]->SetBR(
        branching_ratios[ChargedKaonDecayChannel::mu_nu]);
    (*table)[ChargedKaonDecayChannel::pi_pi0]->SetBR(
        branching_ratios[ChargedKaonDecayChannel::pi_pi0]);
    (*table)[ChargedKaonDecayChannel::pi_pi_pi]->SetBR(
        branching_ratios[ChargedKaonDecayChannel::pi_pi_pi]);
    (*table)[ChargedKaonDecayChannel::pi_pi0_pi0]->SetBR(
        branching_ratios[ChargedKaonDecayChannel::pi_pi0_pi0]);
    (*table)[ChargedKaonDecayChannel::pi0_e_nu]->SetBR(
        branching_ratios[ChargedKaonDecayChannel::pi0_e_nu]);
    (*table)[ChargedKaonDecayChannel::pi0_mu_nu]->SetBR(
        branching_ratios[ChargedKaonDecayChannel::pi0_mu_nu]);
  }
  ldmx_log(trace) << "Decay details (" << kaon->GetParticleName()
                  << ") after setting branching ratios and lifetimes"
                  << std::endl;
  dumpDecayDetails(kaon);
}

Member Data Documentation

k0l_branching_ratios_

std::vector<double> simcore::KaonPhysics::k0l_branching_ratios_

private

Definition at line 106 of file KaonPhysics.h.

k0l_lifetime_factor_

double simcore::KaonPhysics::k0l_lifetime_factor_ {1}

private

Definition at line 99 of file KaonPhysics.h.

{1};

k0s_branching_ratios_

std::vector<double> simcore::KaonPhysics::k0s_branching_ratios_

private

Definition at line 107 of file KaonPhysics.h.

k0s_lifetime_factor_

double simcore::KaonPhysics::k0s_lifetime_factor_ {1}

private

Definition at line 100 of file KaonPhysics.h.

{1};

kminus_branching_ratios_

std::vector<double> simcore::KaonPhysics::kminus_branching_ratios_

private

Definition at line 105 of file KaonPhysics.h.

kminus_lifetime_factor_

double simcore::KaonPhysics::kminus_lifetime_factor_ {1}

private

Definition at line 98 of file KaonPhysics.h.

{1};

kplus_branching_ratios_

std::vector<double> simcore::KaonPhysics::kplus_branching_ratios_

private

Definition at line 104 of file KaonPhysics.h.

kplus_lifetime_factor_

double simcore::KaonPhysics::kplus_lifetime_factor_ {1}

private

Definition at line 97 of file KaonPhysics.h.

{1};

the_log_

framework::logging::logger simcore::KaonPhysics::the_log_

mutableprivate

Initial value:

{
      framework::logging::makeLogger("KaonPhysics")}

Definition at line 145 of file KaonPhysics.h.

                                         {
      framework::logging::makeLogger("KaonPhysics")};

---

The documentation for this class was generated from the following files:

• SimCore/include/SimCore/KaonPhysics.h
• SimCore/src/SimCore/KaonPhysics.cxx