KaonPhysics.cc

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#include <simulation/physicslist/KaonPhysics.h>
 
#include "G4ProcessManager.hh"
#include "G4HadronInelasticProcess.hh"
#include "G4HadronElasticProcess.hh"
#include "G4HadronicAbsorptionBertini.hh"
 
#include "G4CascadeInterface.hh"
#include "G4TheoFSGenerator.hh"
#include "G4FTFModel.hh"
#include "G4ExcitedStringDecay.hh"
#include "G4LundStringFragmentation.hh"
#include "G4GeneratorPrecompoundInterface.hh"
#include "G4HadronElastic.hh"
 
#include "G4CrossSectionElastic.hh"
#include "G4CrossSectionInelastic.hh"
#include "G4ComponentGGHadronNucleusXsc.hh"
#include "G4SystemOfUnits.hh"
 
using namespace Belle2;
using namespace Simulation;
 
 
KaonPhysics::KaonPhysics()
  : m_ftfp(nullptr), m_stringModel(nullptr), m_stringDecay(nullptr),
    m_fragModel(nullptr), m_preCompoundModel(nullptr)
{}
 
 
KaonPhysics::~KaonPhysics()
{
  delete m_stringDecay;
  delete m_stringModel;
  delete m_fragModel;
  delete m_preCompoundModel;
}
 
 
void KaonPhysics::ConstructParticle()
{}
 
 
void KaonPhysics::ConstructProcess()
{
  G4ProcessManager* procMan;
 
  // One elastic model for all kaon energies
  G4HadronElastic* elModel = new G4HadronElastic();
 
  // Use Bertini cascade for low energies
  G4CascadeInterface* loInelModel = new G4CascadeInterface;
  loInelModel->SetMinEnergy(0.0);
  loInelModel->SetMaxEnergy(12.0 * GeV);
 
  // Use FTFP for high energies   ==>>   eventually replace this with new class FTFPInterface
  m_ftfp = new G4TheoFSGenerator("FTFP");
  m_stringModel = new G4FTFModel;
  m_stringDecay =
    new G4ExcitedStringDecay(m_fragModel = new G4LundStringFragmentation);
  m_stringModel->SetFragmentationModel(m_stringDecay);
  m_preCompoundModel = new G4GeneratorPrecompoundInterface();
 
  m_ftfp->SetHighEnergyGenerator(m_stringModel);
  m_ftfp->SetTransport(m_preCompoundModel);
  m_ftfp->SetMinEnergy(10 * GeV);
  m_ftfp->SetMaxEnergy(100 * TeV);
 
  // Inelastic cross section sets
  G4VCrossSectionDataSet* kinelCS =
    new G4CrossSectionInelastic(new G4ComponentGGHadronNucleusXsc);
 
  // Elastic cross section
  G4VCrossSectionDataSet* kelCS =
    new G4CrossSectionElastic(new G4ComponentGGHadronNucleusXsc);

  //   K+                                                                     //
 
  procMan = G4KaonPlus::KaonPlus()->GetProcessManager();
 
  // elastic
  G4HadronElasticProcess* kpProcEl = new G4HadronElasticProcess;
  kpProcEl->RegisterMe(elModel);
  kpProcEl->AddDataSet(kelCS);
  procMan->AddDiscreteProcess(kpProcEl);
 
  // inelastic
  G4HadronInelasticProcess* kpProcInel = new G4HadronInelasticProcess("kaon+Inelastic", G4KaonPlus::Definition());
  kpProcInel->RegisterMe(loInelModel);
  kpProcInel->RegisterMe(m_ftfp);
  kpProcInel->AddDataSet(kinelCS);
  procMan->AddDiscreteProcess(kpProcInel);

  //   K-                                                                     //
 
  procMan = G4KaonMinus::KaonMinus()->GetProcessManager();
 
  // elastic
  G4HadronElasticProcess* kmProcEl = new G4HadronElasticProcess;
  kmProcEl->RegisterMe(elModel);
  kmProcEl->AddDataSet(kelCS);
  procMan->AddDiscreteProcess(kmProcEl);
 
  // inelastic
  G4HadronInelasticProcess* kmProcInel = new G4HadronInelasticProcess("kaon-Inelastic", G4KaonMinus::Definition());
  kmProcInel->RegisterMe(loInelModel);
  kmProcInel->RegisterMe(m_ftfp);
  kmProcInel->AddDataSet(kinelCS);
  procMan->AddDiscreteProcess(kmProcInel);
 
  // stopping
  G4HadronicAbsorptionBertini* bertAbsorb = new G4HadronicAbsorptionBertini;
  procMan->AddRestProcess(bertAbsorb);

  //   K0L                                                                    //
 
  procMan = G4KaonZeroLong::KaonZeroLong()->GetProcessManager();
 
  // elastic
  G4HadronElasticProcess* k0LProcEl = new G4HadronElasticProcess;
  k0LProcEl->RegisterMe(elModel);
  k0LProcEl->AddDataSet(kelCS);
  procMan->AddDiscreteProcess(k0LProcEl);
 
  // inelastic
  G4HadronInelasticProcess* k0LProcInel = new G4HadronInelasticProcess("kaon0LInelastic", G4KaonZeroLong::Definition());
  k0LProcInel->RegisterMe(loInelModel);
  k0LProcInel->RegisterMe(m_ftfp);
  k0LProcInel->AddDataSet(kinelCS);
  procMan->AddDiscreteProcess(k0LProcInel);

  //   K0S                                                                    //
 
  procMan = G4KaonZeroShort::KaonZeroShort()->GetProcessManager();
 
  // elastic
  G4HadronElasticProcess* k0SProcEl = new G4HadronElasticProcess;
  k0SProcEl->RegisterMe(elModel);
  k0SProcEl->AddDataSet(kelCS);
  procMan->AddDiscreteProcess(k0SProcEl);
 
  // inelastic
  G4HadronInelasticProcess* k0SProcInel = new G4HadronInelasticProcess("kaon0SInelastic", G4KaonZeroShort::Definition());
  k0SProcInel->RegisterMe(loInelModel);
  k0SProcInel->RegisterMe(m_ftfp);
  k0SProcInel->AddDataSet(kinelCS);
  procMan->AddDiscreteProcess(k0SProcInel);
}