Belle2PhysicsList Class Reference

Custom Geant4 physics list for Belle II with options to add optical physics, standard EM physics and high precision neutrons. More...

#include <Belle2PhysicsList.h>

Inheritance diagram for Belle2PhysicsList:

Collaboration diagram for Belle2PhysicsList:

Public Member Functions
	Belle2PhysicsList (const G4String &physicsListName)
	Constructor.
	~Belle2PhysicsList ()
	destructor
virtual void	ConstructParticle ()
	Build all particle types used in physics list.
virtual void	SetCuts ()
	Set the secondary particle production thresholds.
void	SetVerbosity (G4int verb)
	Run/event verbosity level.
void	SetProductionCutValue (G4double)
	Use parameter to set global cut value.
void	SetPXDProductionCutValue (G4double)
	Set cut value for PXD envelope.
void	SetSVDProductionCutValue (G4double)
	Set cut value for SVD envelope.
void	SetCDCProductionCutValue (G4double)
	Set cut value for CDC envelope.
void	SetARICHTOPProductionCutValue (G4double)
	Set cut value for ARICH and TOP envelopes.
void	SetECLProductionCutValue (G4double)
	Set cut value for ECL barrel, forward and backward envelopes.
void	SetKLMProductionCutValue (G4double)
	Set cut value for BKLM and EKLM envelopes.
void	UseStandardEMPhysics (G4bool)
	Use standard EM physics instead of EM option1.
void	UseOpticalPhysics (G4bool)
	Add optical photon physics.
void	UseHighPrecisionNeutrons (G4bool)
	Use high precision neutron models below 20 MeV.
void	UseLongLivedNeutralParticles ()
	Simulate neutral long-lived particles with given pdg and mass value.

Private Member Functions
void	setRegionCuts (const std::string &name, const std::vector< std::string > &regions, double cutValue)
	Set the produciton cuts to the given value for a list of regions belonging to a sub detector. More...
void	ConstructG4eParticles ()
	Construct parallel particle types needed for reco.

Private Attributes
G4double	m_globalCutValue
	Secondary production thresholds.
G4double	m_pxdCutValue
	threshold for PXD
G4double	m_svdCutValue
	threshold for SVD
G4double	m_cdcCutValue
	threshold for CDC
G4double	m_arichtopCutValue
	threshold for ARICH and TOP
G4double	m_eclCutValue
	threshold for ECL
G4double	m_klmCutValue
	threshold for BKLM and EKLM

Detailed Description

Custom Geant4 physics list for Belle II with options to add optical physics, standard EM physics and high precision neutrons.

Definition at line 37 of file Belle2PhysicsList.h.

Member Function Documentation

setRegionCuts()

void setRegionCuts ( const std::string &  name, const std::vector< std::string > &  regions, double  cutValue  )

private

Set the produciton cuts to the given value for a list of regions belonging to a sub detector.

Parameters

name	name of the sub detector to print in messages
regions	list of region names to set the cut value
cutValue	production cut in cm, if 0 the global cut value will be used

Definition at line 236 of file Belle2PhysicsList.cc.

{
  G4RegionStore* theRegionStore = G4RegionStore::GetInstance();
  if (cutValue == 0) {
    cutValue = m_globalCutValue;
  } else {
    B2INFO("Set production cut for detector region" << LogVar("detector", name) << LogVar("production_cut", cutValue));
  }
  auto* regionCuts = new G4ProductionCuts;
  regionCuts->SetProductionCut(cutValue * cm);
  bool foundOne{false};
  for (const auto& regionName : regions) {
    auto* region = theRegionStore->GetRegion(regionName, false);
    if (!region) {
      B2WARNING("Cannot find Geant4 region for sub detector. Probably detector not present?"
                << LogVar("detector", name) << LogVar("region", regionName));
      continue;
    }
    region->SetProductionCuts(regionCuts);
    foundOne = true;
  }
  if (!foundOne) delete regionCuts;
}

---

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

• simulation/physicslist/include/Belle2PhysicsList.h
• simulation/physicslist/src/Belle2PhysicsList.cc