4 #include "G4EmParameters.hh" 5 #include "G4ParticleTable.hh" 7 #include "G4ParticleDefinition.hh" 8 #include "G4LossTableManager.hh" 10 #include "G4ComptonScattering.hh" 11 #include "G4GammaConversion.hh" 12 #include "G4PhotoElectricEffect.hh" 14 #include "G4hMultipleScattering.hh" 15 #include "G4eMultipleScattering.hh" 16 #include "G4MuMultipleScattering.hh" 17 #include "G4CoulombScattering.hh" 18 #include "G4eCoulombScatteringModel.hh" 19 #include "G4WentzelVIModel.hh" 21 #include "G4eIonisation.hh" 22 #include "G4eBremsstrahlung.hh" 23 #include "G4eplusAnnihilation.hh" 24 #include "G4UAtomicDeexcitation.hh" 26 #include "G4MuIonisation.hh" 27 #include "G4MuBremsstrahlung.hh" 28 #include "G4MuPairProduction.hh" 30 #include "G4hIonisation.hh" 31 #include "G4ionIonisation.hh" 32 #include "G4hBremsstrahlung.hh" 33 #include "G4hPairProduction.hh" 36 #include "G4Electron.hh" 37 #include "G4Positron.hh" 38 #include "G4MuonPlus.hh" 39 #include "G4MuonMinus.hh" 40 #include "G4TauMinus.hh" 41 #include "G4TauPlus.hh" 42 #include "G4PionPlus.hh" 43 #include "G4PionMinus.hh" 44 #include "G4KaonPlus.hh" 45 #include "G4KaonMinus.hh" 46 #include "G4BMesonMinus.hh" 47 #include "G4BMesonPlus.hh" 48 #include "G4DMesonMinus.hh" 49 #include "G4DMesonPlus.hh" 50 #include "G4Proton.hh" 51 #include "G4AntiProton.hh" 52 #include "G4SigmaMinus.hh" 53 #include "G4AntiSigmaMinus.hh" 54 #include "G4SigmaPlus.hh" 55 #include "G4AntiSigmaPlus.hh" 56 #include "G4XiMinus.hh" 57 #include "G4AntiXiMinus.hh" 58 #include "G4OmegaMinus.hh" 59 #include "G4AntiOmegaMinus.hh" 60 #include "G4LambdacPlus.hh" 61 #include "G4AntiLambdacPlus.hh" 62 #include "G4XicPlus.hh" 63 #include "G4AntiXicPlus.hh" 64 #include "G4Deuteron.hh" 65 #include "G4Triton.hh" 68 #include "G4GenericIon.hh" 70 #include "G4PhysicsListHelper.hh" 71 #include "G4BuilderType.hh" 73 #include "G4SystemOfUnits.hh" 76 G4VPhysicsConstructor(
"CMSEmStandard_opt1"),
verbose(ver) {
77 G4EmParameters* param = G4EmParameters::Instance();
80 param->SetApplyCuts(
true);
81 param->SetMscRangeFactor(0.2);
82 param->SetMscStepLimitType(fMinimal);
83 SetPhysicsType(bElectromagnetic);
94 G4Positron::Positron();
95 G4MuonPlus::MuonPlus();
96 G4MuonMinus::MuonMinus();
97 G4TauMinus::TauMinusDefinition();
98 G4TauPlus::TauPlusDefinition();
101 G4PionPlus::PionPlusDefinition();
102 G4PionMinus::PionMinusDefinition();
103 G4KaonPlus::KaonPlusDefinition();
104 G4KaonMinus::KaonMinusDefinition();
105 G4DMesonMinus::DMesonMinusDefinition();
106 G4DMesonPlus::DMesonPlusDefinition();
107 G4BMesonMinus::BMesonMinusDefinition();
108 G4BMesonPlus::BMesonPlusDefinition();
112 G4AntiProton::AntiProton();
113 G4SigmaMinus::SigmaMinusDefinition();
114 G4AntiSigmaMinus::AntiSigmaMinusDefinition();
115 G4SigmaPlus::SigmaPlusDefinition();
116 G4AntiSigmaPlus::AntiSigmaPlusDefinition();
117 G4XiMinus::XiMinusDefinition();
118 G4AntiXiMinus::AntiXiMinusDefinition();
119 G4OmegaMinus::OmegaMinusDefinition();
120 G4AntiOmegaMinus::AntiOmegaMinusDefinition();
121 G4LambdacPlus::LambdacPlusDefinition();
122 G4AntiLambdacPlus::AntiLambdacPlusDefinition();
123 G4XicPlus::XicPlusDefinition();
124 G4AntiXicPlus::AntiXicPlusDefinition();
127 G4Deuteron::Deuteron();
131 G4GenericIon::GenericIonDefinition();
137 G4cout <<
"### " << GetPhysicsName() <<
" Construct Processes " << G4endl;
144 G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
147 G4MuBremsstrahlung* mub =
nullptr;
148 G4MuPairProduction* mup =
nullptr;
149 G4hBremsstrahlung* pib =
nullptr;
150 G4hPairProduction* pip =
nullptr;
151 G4hBremsstrahlung* kb =
nullptr;
152 G4hPairProduction*
kp =
nullptr;
153 G4hBremsstrahlung* pb =
nullptr;
154 G4hPairProduction*
pp =
nullptr;
157 G4MuMultipleScattering* mumsc=
nullptr;
158 G4hMultipleScattering* pimsc=
nullptr;
159 G4hMultipleScattering* kmsc =
nullptr;
160 G4hMultipleScattering* pmsc =
nullptr;
161 G4hMultipleScattering* hmsc =
nullptr;
164 G4double highEnergyLimit = 100*
MeV;
166 G4ParticleTable* table = G4ParticleTable::GetParticleTable();
169 G4ParticleDefinition* particle = table->FindParticle(
particleName);
173 ph->RegisterProcess(
new G4PhotoElectricEffect(), particle);
174 ph->RegisterProcess(
new G4ComptonScattering(), particle);
175 ph->RegisterProcess(
new G4GammaConversion(), particle);
179 G4eIonisation* eioni =
new G4eIonisation();
180 eioni->SetStepFunction(0.8, 1.0*mm);
182 G4eMultipleScattering*
msc =
new G4eMultipleScattering;
183 msc->SetStepLimitType(fMinimal);
185 G4WentzelVIModel* msc2 =
new G4WentzelVIModel();
186 msc1->SetHighEnergyLimit(highEnergyLimit);
187 msc2->SetLowEnergyLimit(highEnergyLimit);
188 msc->AddEmModel(0, msc1);
189 msc->AddEmModel(0, msc2);
191 G4eCoulombScatteringModel* ssm =
new G4eCoulombScatteringModel();
192 G4CoulombScattering* ss =
new G4CoulombScattering();
194 ss->SetMinKinEnergy(highEnergyLimit);
195 ssm->SetLowEnergyLimit(highEnergyLimit);
196 ssm->SetActivationLowEnergyLimit(highEnergyLimit);
198 ph->RegisterProcess(msc, particle);
199 ph->RegisterProcess(eioni, particle);
200 ph->RegisterProcess(
new G4eBremsstrahlung(), particle);
201 ph->RegisterProcess(ss, particle);
205 G4eIonisation* eioni =
new G4eIonisation();
206 eioni->SetStepFunction(0.8, 1.0*mm);
208 G4eMultipleScattering*
msc =
new G4eMultipleScattering;
209 msc->SetStepLimitType(fMinimal);
211 G4WentzelVIModel* msc2 =
new G4WentzelVIModel();
212 msc1->SetHighEnergyLimit(highEnergyLimit);
213 msc2->SetLowEnergyLimit(highEnergyLimit);
214 msc->AddEmModel(0, msc1);
215 msc->AddEmModel(0, msc2);
217 G4eCoulombScatteringModel* ssm =
new G4eCoulombScatteringModel();
218 G4CoulombScattering* ss =
new G4CoulombScattering();
220 ss->SetMinKinEnergy(highEnergyLimit);
221 ssm->SetLowEnergyLimit(highEnergyLimit);
222 ssm->SetActivationLowEnergyLimit(highEnergyLimit);
224 ph->RegisterProcess(msc, particle);
225 ph->RegisterProcess(eioni, particle);
226 ph->RegisterProcess(
new G4eBremsstrahlung(), particle);
227 ph->RegisterProcess(
new G4eplusAnnihilation(), particle);
228 ph->RegisterProcess(ss, particle);
234 mub =
new G4MuBremsstrahlung();
235 mup =
new G4MuPairProduction();
236 mumsc =
new G4MuMultipleScattering();
237 mumsc->AddEmModel(0,
new G4WentzelVIModel());
239 ph->RegisterProcess(mumsc, particle);
240 ph->RegisterProcess(
new G4MuIonisation(), particle);
241 ph->RegisterProcess(mub, particle);
242 ph->RegisterProcess(mup, particle);
243 ph->RegisterProcess(
new G4CoulombScattering(), particle);
248 ph->RegisterProcess(
new G4hMultipleScattering(), particle);
249 ph->RegisterProcess(
new G4ionIonisation(), particle);
253 if(
nullptr == hmsc) {
254 hmsc =
new G4hMultipleScattering(
"ionmsc");
256 ph->RegisterProcess(hmsc, particle);
257 ph->RegisterProcess(
new G4ionIonisation(), particle);
263 pib =
new G4hBremsstrahlung();
264 pip =
new G4hPairProduction();
265 pimsc =
new G4hMultipleScattering();
266 pimsc->AddEmModel(0,
new G4WentzelVIModel());
268 ph->RegisterProcess(pimsc, particle);
269 ph->RegisterProcess(
new G4hIonisation(), particle);
270 ph->RegisterProcess(pib, particle);
271 ph->RegisterProcess(pip, particle);
272 ph->RegisterProcess(
new G4CoulombScattering(), particle);
278 kb =
new G4hBremsstrahlung();
279 kp =
new G4hPairProduction();
280 kmsc =
new G4hMultipleScattering();
281 kmsc->AddEmModel(0,
new G4WentzelVIModel());
283 ph->RegisterProcess(kmsc, particle);
284 ph->RegisterProcess(
new G4hIonisation(), particle);
285 ph->RegisterProcess(kb, particle);
286 ph->RegisterProcess(kp, particle);
287 ph->RegisterProcess(
new G4CoulombScattering(), particle);
293 pb =
new G4hBremsstrahlung();
294 pp =
new G4hPairProduction();
295 pmsc =
new G4hMultipleScattering();
296 pmsc->AddEmModel(0,
new G4WentzelVIModel());
299 ph->RegisterProcess(pmsc, particle);
300 ph->RegisterProcess(
new G4hIonisation(), particle);
301 ph->RegisterProcess(pb, particle);
302 ph->RegisterProcess(pp, particle);
303 ph->RegisterProcess(
new G4CoulombScattering(), particle);
336 if(
nullptr == hmsc) {
337 hmsc =
new G4hMultipleScattering(
"ionmsc");
339 ph->RegisterProcess(hmsc, particle);
340 ph->RegisterProcess(
new G4hIonisation(), particle);
const std::vector< G4String > & PartNames() const
void ConstructParticle() override
void ConstructProcess() override
~CMSEmStandardPhysics() override
CMSEmStandardPhysics(G4int ver)