Belle II Software development
Public Member Functions | Private Member Functions | Private Attributes | List of all members
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:

Public Member Functions

 Belle2PhysicsList (const G4String &physicsListName, const G4int hadronicVerbosityLevel=0)
 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.
 
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
 
G4ProductionCuts * m_regionCuts
 Global pointer for the region cuts to avoid memory leak and side effects when deleting the pointer.
 

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 27 of file Belle2PhysicsList.h.

Constructor & Destructor Documentation

◆ Belle2PhysicsList()

Belle2PhysicsList ( const G4String &  physicsListName,
const G4int  hadronicVerbosityLevel = 0 
)
explicit

Constructor.

Definition at line 49 of file Belle2PhysicsList.cc.

50 : G4VModularPhysicsList(), m_globalCutValue(0.07)
51{
52 B2INFO("Using " << physicsListName << " physics list");
53
54 m_pxdCutValue = m_globalCutValue;
55 m_svdCutValue = m_globalCutValue;
56 m_cdcCutValue = m_globalCutValue;
57 m_arichtopCutValue = m_globalCutValue;
58 m_eclCutValue = m_globalCutValue;
59 m_klmCutValue = m_globalCutValue;
60
61 // Decay
62 RegisterPhysics(new G4DecayPhysics());
63
64 // Hadronic physics
65 RegisterPhysics(new ProtonPhysics());
66 RegisterPhysics(new NeutronPhysics());
67 RegisterPhysics(new PionPhysics());
68 RegisterPhysics(new KaonPhysics());
69 RegisterPhysics(new HyperonPhysics());
70 RegisterPhysics(new AntiBaryonPhysics());
71 RegisterPhysics(new IonPhysics());
72 RegisterPhysics(new GammaLeptoNuclearPhysics(hadronicVerbosityLevel));
73 RegisterPhysics(new Geant4ePhysics());
74
75 m_regionCuts = new G4ProductionCuts;
76}
Belle2::Simulation::Belle2PhysicsList::m_regionCuts
G4ProductionCuts * m_regionCuts
Global pointer for the region cuts to avoid memory leak and side effects when deleting the pointer.
Definition: Belle2PhysicsList.h:98
Belle2::Simulation::Belle2PhysicsList::m_klmCutValue
G4double m_klmCutValue
threshold for BKLM and EKLM
Definition: Belle2PhysicsList.h:96
Belle2::Simulation::Belle2PhysicsList::m_globalCutValue
G4double m_globalCutValue
Secondary production thresholds.
Definition: Belle2PhysicsList.h:89
Belle2::Simulation::Belle2PhysicsList::m_pxdCutValue
G4double m_pxdCutValue
threshold for PXD
Definition: Belle2PhysicsList.h:91
Belle2::Simulation::Belle2PhysicsList::m_svdCutValue
G4double m_svdCutValue
threshold for SVD
Definition: Belle2PhysicsList.h:92
Belle2::Simulation::Belle2PhysicsList::m_arichtopCutValue
G4double m_arichtopCutValue
threshold for ARICH and TOP
Definition: Belle2PhysicsList.h:94
Belle2::Simulation::Belle2PhysicsList::m_eclCutValue
G4double m_eclCutValue
threshold for ECL
Definition: Belle2PhysicsList.h:95
Belle2::Simulation::Belle2PhysicsList::m_cdcCutValue
G4double m_cdcCutValue
threshold for CDC
Definition: Belle2PhysicsList.h:93

◆ ~Belle2PhysicsList()

~Belle2PhysicsList ( )

destructor

Definition at line 79 of file Belle2PhysicsList.cc.

80{
81 delete m_regionCuts;
82}

Member Function Documentation

◆ ConstructG4eParticles()

void ConstructG4eParticles ( )
private

Construct parallel particle types needed for reco.

Definition at line 109 of file Belle2PhysicsList.cc.

110{
111 static G4ParticleDefinition* g4eParticle = NULL;
112
113 if (g4eParticle == NULL) {
114 // Bohr Magneton for positron and positive muon
115 G4double muBpositron = 0.5 * CLHEP::eplus * CLHEP::hbar_Planck / (0.51099906 * CLHEP::MeV / CLHEP::c_squared);
116 G4double muBmuon = 0.5 * CLHEP::eplus * CLHEP::hbar_Planck / (0.1056584 * CLHEP::GeV / CLHEP::c_squared);
117 G4double muNucleon = CLHEP::eplus * CLHEP::hbar_Planck / 2. / (CLHEP::proton_mass_c2 / CLHEP::c_squared);
118 // Copied from G4Gamma.cc
119 g4eParticle = new G4ParticleDefinition(
120 "g4e_gamma", 0.0 * CLHEP::MeV, 0.0 * CLHEP::MeV, 0.0,
121 2, -1, -1,
122 0, 0, 0,
123 "gamma", 0, 0, g4ePDGcode,
124 true, 0.0, NULL,
125 false, "photon", g4ePDGcode
126 );
127 // Copied from G4Electron.cc
128 new G4ParticleDefinition(
129 "g4e_e-", 0.51099906 * CLHEP::MeV, 0.0 * CLHEP::MeV, -1.0 * CLHEP::eplus,
130 1, 0, 0,
131 0, 0, 0,
132 "lepton", 1, 0, g4ePDGcode,
133 true, -1.0, NULL,
134 false, "e", g4ePDGcode, -1.0011596521859 * muBpositron
135 );
136 // Copied from G4Positron.cc
137 new G4ParticleDefinition(
138 "g4e_e+", 0.51099906 * CLHEP::MeV, 0.0 * CLHEP::MeV, +1.0 * CLHEP::eplus,
139 1, 0, 0,
140 0, 0, 0,
141 "lepton", -1, 0, g4ePDGcode,
142 true, -1.0, NULL,
143 false, "e", g4ePDGcode, 1.0011596521859 * muBpositron
144 );
145 // Copied from G4MuonPlus.cc
146 new G4ParticleDefinition(
147 "g4e_mu+", 0.1056584 * CLHEP::GeV, 2.99591e-16 * CLHEP::MeV, +1.0 * CLHEP::eplus,
148 1, 0, 0,
149 0, 0, 0,
150 "lepton", -1, 0, g4ePDGcode,
151 true, -1.0, NULL,
152 false, "mu", g4ePDGcode, 1.0011659208 * muBmuon
153 );
154 // Copied from G4MuonMinus.cc
155 new G4ParticleDefinition(
156 "g4e_mu-", 0.1056584 * CLHEP::GeV, 2.99591e-16 * CLHEP::MeV, -1.0 * CLHEP::eplus,
157 1, 0, 0,
158 0, 0, 0,
159 "lepton", 1, 0, g4ePDGcode,
160 true, -1.0, NULL,
161 false, "mu", g4ePDGcode, -1.0011659208 * muBmuon
162 );
163 // Copied from G4PionPlus.cc
164 new G4ParticleDefinition(
165 "g4e_pi+", 0.1395700 * CLHEP::GeV, 2.5284e-14 * CLHEP::MeV, +1.0 * CLHEP::eplus,
166 0, -1, 0,
167 2, +2, -1,
168 "meson", 0, 0, g4ePDGcode,
169 true, -1.0, NULL,
170 false, "pi", g4ePDGcode
171 );
172 // Copied from G4PionMinus.cc
173 new G4ParticleDefinition(
174 "g4e_pi-", 0.1395700 * CLHEP::GeV, 2.5284e-14 * CLHEP::MeV, -1.0 * CLHEP::eplus,
175 0, -1, 0,
176 2, -2, -1,
177 "meson", 0, 0, g4ePDGcode,
178 true, -1.0, NULL,
179 false, "pi", g4ePDGcode
180 );
181 // Copied from G4KaonPlus.cc
182 new G4ParticleDefinition(
183 "g4e_kaon+", 0.493677 * CLHEP::GeV, 5.315e-14 * CLHEP::MeV, +1.0 * CLHEP::eplus,
184 0, -1, 0,
185 1, +1, 0,
186 "meson", 0, 0, g4ePDGcode,
187 true, -1.0, NULL,
188 false, "kaon", g4ePDGcode
189 );
190 // Copied from G4KaonMinus.cc
191 new G4ParticleDefinition(
192 "g4e_kaon-", 0.493677 * CLHEP::GeV, 5.315e-14 * CLHEP::MeV, -1.0 * CLHEP::eplus,
193 0, -1, 0,
194 1, -1, 0,
195 "meson", 0, 0, g4ePDGcode,
196 true, -1.0, NULL,
197 false, "kaon", g4ePDGcode
198 );
199 // Copied from G4Proton.cc except use G4ParticleDefinition instead of G4Ions
200 new G4ParticleDefinition(
201 "g4e_proton", 0.9382723 * CLHEP::GeV, 0.0 * CLHEP::MeV, +1.0 * CLHEP::eplus,
202 1, +1, 0,
203 1, +1, 0,
204 "baryon", 0, +1, g4ePDGcode,
205 true, -1.0, NULL,
206 false, "nucleon", g4ePDGcode, 2.792847351 * muNucleon
207 );
208 // Copied from G4AntiProton.cc except use G4ParticleDefinition instead of G4Ions
209 new G4ParticleDefinition(
210 "g4e_anti_proton", 0.9382723 * CLHEP::GeV, 0.0 * CLHEP::MeV, -1.0 * CLHEP::eplus,
211 1, +1, 0,
212 1, -1, 0,
213 "baryon", 0, -1, g4ePDGcode,
214 true, -1.0, NULL,
215 false, "nucleon", g4ePDGcode, -2.792847351 * muNucleon
216 );
217 // copied from G4Deuteron.hh except use G4ParticleDefinition instead of G4Ions
218 new G4ParticleDefinition(
219 "g4e_deuteron", 1.875613 * CLHEP::GeV, 0.0 * CLHEP::MeV, +1.0 * CLHEP::eplus,
220 2, +1, 0,
221 0, 0, 0,
222 "nucleus", 0, +2, g4ePDGcode,
223 true, -1.0, NULL,
224 false, "static", g4ePDGcode, 0.857438230 * muNucleon
225 );
226 // copied from G4AntiDeuteron.hh except use G4ParticleDefinition instead of G4Ions
227 new G4ParticleDefinition(
228 "g4e_anti_deuteron", 1.875613 * CLHEP::GeV, 0.0 * CLHEP::MeV, -1.0 * CLHEP::eplus,
229 2, +1, 0,
230 0, 0, 0,
231 "anti_nucleus", 0, -2, g4ePDGcode,
232 true, -1.0, NULL,
233 false, "static", g4ePDGcode, -0.857438230 * muNucleon
234 );
235 }
236}

◆ ConstructParticle()

void ConstructParticle ( )
virtual

Build all particle types used in physics list.

Definition at line 85 of file Belle2PhysicsList.cc.

86{
87 G4BosonConstructor pBosonConstructor;
88 pBosonConstructor.ConstructParticle();
89
90 G4LeptonConstructor pLeptonConstructor;
91 pLeptonConstructor.ConstructParticle();
92
93 G4MesonConstructor pMesonConstructor;
94 pMesonConstructor.ConstructParticle();
95
96 G4BaryonConstructor pBaryonConstructor;
97 pBaryonConstructor.ConstructParticle();
98
99 G4IonConstructor pIonConstructor;
100 pIonConstructor.ConstructParticle();
101
102 G4ShortLivedConstructor pShortLivedConstructor;
103 pShortLivedConstructor.ConstructParticle();
104
105 ConstructG4eParticles();
106}
Belle2::Simulation::Belle2PhysicsList::ConstructG4eParticles
void ConstructG4eParticles()
Construct parallel particle types needed for reco.
Definition: Belle2PhysicsList.cc:109

◆ SetARICHTOPProductionCutValue()

void SetARICHTOPProductionCutValue ( G4double  value)

Set cut value for ARICH and TOP envelopes.

Definition at line 309 of file Belle2PhysicsList.cc.

310{
311 m_arichtopCutValue = value;
312}

◆ SetCDCProductionCutValue()

void SetCDCProductionCutValue ( G4double  value)

Set cut value for CDC envelope.

Definition at line 303 of file Belle2PhysicsList.cc.

304{
305 m_cdcCutValue = value;
306}

◆ SetCuts()

void SetCuts ( )
virtual

Set the secondary particle production thresholds.

Definition at line 259 of file Belle2PhysicsList.cc.

260{
261 // Belle2 assumes input units are cm
262 B2INFO("Setting global Geant4 production cut" << LogVar("cutValue", m_globalCutValue));
263
264 SetCutValue(m_globalCutValue * cm, "proton");
265 SetCutValue(m_globalCutValue * cm, "e-");
266 SetCutValue(m_globalCutValue * cm, "e+");
267 SetCutValue(m_globalCutValue * cm, "gamma");
268
269 setRegionCuts("PXD", {"PXDEnvelope"}, m_pxdCutValue);
270 setRegionCuts("SVD", {"SVDEnvelope"}, m_svdCutValue);
271 setRegionCuts("CDC", {"CDCEnvelope"}, m_cdcCutValue);
272 setRegionCuts("ARICH", {"ARICHEnvelope"}, m_arichtopCutValue);
273 setRegionCuts("TOP", {"TOPEnvelope"}, m_arichtopCutValue);
274 setRegionCuts("ECL", {"ECLForwardEnvelope", "ECLBarrelSector", "ECLBackwardEnvelope"}, m_eclCutValue);
275 setRegionCuts("KLM", {"BKLMEnvelope", "EKLMEnvelope"}, m_klmCutValue);
276}
Belle2::Simulation::Belle2PhysicsList::setRegionCuts
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.
Definition: Belle2PhysicsList.cc:238
LogVar
Class to store variables with their name which were sent to the logging service.
Definition: LogVariableStream.h:24

◆ SetECLProductionCutValue()

void SetECLProductionCutValue ( G4double  value)

Set cut value for ECL barrel, forward and backward envelopes.

Definition at line 315 of file Belle2PhysicsList.cc.

316{
317 m_eclCutValue = value;
318}

◆ SetKLMProductionCutValue()

void SetKLMProductionCutValue ( G4double  value)

Set cut value for BKLM and EKLM envelopes.

Definition at line 321 of file Belle2PhysicsList.cc.

322{
323 m_klmCutValue = value;
324}

◆ SetProductionCutValue()

void SetProductionCutValue ( G4double  value)

Use parameter to set global cut value.

Definition at line 285 of file Belle2PhysicsList.cc.

286{
287 m_globalCutValue = value;
288}

◆ SetPXDProductionCutValue()

void SetPXDProductionCutValue ( G4double  value)

Set cut value for PXD envelope.

Definition at line 291 of file Belle2PhysicsList.cc.

292{
293 m_pxdCutValue = value;
294}

◆ 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
namename of the sub detector to print in messages
regionslist of region names to set the cut value
cutValueproduction cut in cm, if 0 the global cut value will be used

Definition at line 238 of file Belle2PhysicsList.cc.

239{
240 G4RegionStore* theRegionStore = G4RegionStore::GetInstance();
241 if (cutValue == 0) {
242 cutValue = m_globalCutValue;
243 } else {
244 B2INFO("Set production cut for detector region" << LogVar("detector", name) << LogVar("production_cut", cutValue));
245 }
246 m_regionCuts->SetProductionCut(cutValue * cm);
247 for (const auto& regionName : regions) {
248 auto* region = theRegionStore->GetRegion(regionName, false);
249 if (!region) {
250 B2WARNING("Cannot find Geant4 region for sub detector. Probably detector not present?"
251 << LogVar("detector", name) << LogVar("region", regionName));
252 continue;
253 }
254 region->SetProductionCuts(m_regionCuts);
255 }
256}

◆ SetSVDProductionCutValue()

void SetSVDProductionCutValue ( G4double  value)

Set cut value for SVD envelope.

Definition at line 297 of file Belle2PhysicsList.cc.

298{
299 m_svdCutValue = value;
300}

◆ SetVerbosity()

void SetVerbosity ( G4int  verb)

Run/event verbosity level.

Definition at line 279 of file Belle2PhysicsList.cc.

280{
281 SetVerboseLevel(verb);
282}

◆ UseHighPrecisionNeutrons()

void UseHighPrecisionNeutrons ( G4bool  yesno)

Use high precision neutron models below 20 MeV.

Definition at line 343 of file Belle2PhysicsList.cc.

344{
345 if (yesno) G4cout << " High precision neutron option not yet ready " << G4endl;
346}

◆ UseLongLivedNeutralParticles()

void UseLongLivedNeutralParticles ( )

Simulate neutral long-lived particles with given pdg and mass value.

Definition at line 348 of file Belle2PhysicsList.cc.

349{
350 G4LongLivedNeutralPhysics* pLongLivedNeutral = new G4LongLivedNeutralPhysics();
351 RegisterPhysics(pLongLivedNeutral);
352 pLongLivedNeutral->ConstructParticle();
353}
Belle2::G4LongLivedNeutralPhysics
LongLivedNeutral physics Class – to be registered in the physics list.
Definition: G4LongLivedNeutralPhysics.h:27
Belle2::G4LongLivedNeutralPhysics::ConstructParticle
virtual void ConstructParticle()
Adds monopole and anti-monopole to GEANT4 with a pdg of +/-99666 and parameters taken from current cl...
Definition: G4LongLivedNeutralPhysics.cc:35

◆ UseOpticalPhysics()

void UseOpticalPhysics ( G4bool  yesno)

Add optical photon physics.

Definition at line 337 of file Belle2PhysicsList.cc.

338{
339 if (yesno) RegisterPhysics(new G4OpticalPhysics());
340}

◆ UseStandardEMPhysics()

void UseStandardEMPhysics ( G4bool  yesno)

Use standard EM physics instead of EM option1.

Definition at line 327 of file Belle2PhysicsList.cc.

328{
329 if (yesno) {
330 RegisterPhysics(new G4EmStandardPhysics());
331 } else {
332 RegisterPhysics(new G4EmStandardPhysics_option1());
333 }
334}

Member Data Documentation

◆ m_arichtopCutValue

G4double m_arichtopCutValue
private

threshold for ARICH and TOP

Definition at line 94 of file Belle2PhysicsList.h.

◆ m_cdcCutValue

G4double m_cdcCutValue
private

threshold for CDC

Definition at line 93 of file Belle2PhysicsList.h.

◆ m_eclCutValue

G4double m_eclCutValue
private

threshold for ECL

Definition at line 95 of file Belle2PhysicsList.h.

◆ m_globalCutValue

G4double m_globalCutValue
private

Secondary production thresholds.

Definition at line 89 of file Belle2PhysicsList.h.

◆ m_klmCutValue

G4double m_klmCutValue
private

threshold for BKLM and EKLM

Definition at line 96 of file Belle2PhysicsList.h.

◆ m_pxdCutValue

G4double m_pxdCutValue
private

threshold for PXD

Definition at line 91 of file Belle2PhysicsList.h.

◆ m_regionCuts

G4ProductionCuts* m_regionCuts
private

Global pointer for the region cuts to avoid memory leak and side effects when deleting the pointer.

Definition at line 98 of file Belle2PhysicsList.h.

◆ m_svdCutValue

G4double m_svdCutValue
private

threshold for SVD

Definition at line 92 of file Belle2PhysicsList.h.

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