3 #include "G4ParticleDefinition.hh"
4 #include "G4ProcessManager.hh"
5 #include "G4LossTableManager.hh"
6 #include "G4EmProcessOptions.hh"
7 #include "G4RegionStore.hh"
9 #include "G4ComptonScattering.hh"
10 #include "G4GammaConversion.hh"
11 #include "G4PhotoElectricEffect.hh"
12 #include "G4PairProductionRelModel.hh"
14 #include "G4hMultipleScattering.hh"
15 #include "G4eMultipleScattering.hh"
16 #include "G4MuMultipleScattering.hh"
17 #include "G4MscStepLimitType.hh"
18 #include "G4UrbanMscModel93.hh"
19 #include "G4WentzelVIModel.hh"
20 #include "G4CoulombScattering.hh"
22 #include "G4eIonisation.hh"
23 #include "G4eBremsstrahlung.hh"
24 #include "G4eplusAnnihilation.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"
71 G4VPhysicsConstructor(name),
verbose(ver) {
72 G4LossTableManager::Instance();
82 G4Electron::Electron();
83 G4Positron::Positron();
84 G4MuonPlus::MuonPlus();
85 G4MuonMinus::MuonMinus();
86 G4TauMinus::TauMinusDefinition();
87 G4TauPlus::TauPlusDefinition();
90 G4PionPlus::PionPlusDefinition();
91 G4PionMinus::PionMinusDefinition();
92 G4KaonPlus::KaonPlusDefinition();
93 G4KaonMinus::KaonMinusDefinition();
94 G4DMesonMinus::DMesonMinusDefinition();
95 G4DMesonPlus::DMesonPlusDefinition();
96 G4BMesonMinus::BMesonMinusDefinition();
97 G4BMesonPlus::BMesonPlusDefinition();
101 G4AntiProton::AntiProton();
102 G4SigmaMinus::SigmaMinusDefinition();
103 G4AntiSigmaMinus::AntiSigmaMinusDefinition();
104 G4SigmaPlus::SigmaPlusDefinition();
105 G4AntiSigmaPlus::AntiSigmaPlusDefinition();
106 G4XiMinus::XiMinusDefinition();
107 G4AntiXiMinus::AntiXiMinusDefinition();
108 G4OmegaMinus::OmegaMinusDefinition();
109 G4AntiOmegaMinus::AntiOmegaMinusDefinition();
110 G4LambdacPlus::LambdacPlusDefinition();
111 G4AntiLambdacPlus::AntiLambdacPlusDefinition();
112 G4XicPlus::XicPlusDefinition();
113 G4AntiXicPlus::AntiXicPlusDefinition();
116 G4Deuteron::Deuteron();
120 G4GenericIon::GenericIonDefinition();
127 theParticleIterator->reset();
128 while( (*theParticleIterator)() ){
129 G4ParticleDefinition* particle = theParticleIterator->value();
130 G4ProcessManager* pmanager = particle->GetProcessManager();
131 G4String particleName = particle->GetParticleName();
133 G4cout <<
"### " << GetPhysicsName() <<
" instantiates for "
134 << particleName <<
" at " << particle << G4endl;
136 if (particleName ==
"gamma") {
138 pmanager->AddDiscreteProcess(
new G4PhotoElectricEffect);
139 pmanager->AddDiscreteProcess(
new G4ComptonScattering);
140 G4GammaConversion*
conv =
new G4GammaConversion();
141 G4PairProductionRelModel*
mod =
new G4PairProductionRelModel();
142 mod->SetLowEnergyLimit(80*GeV);
143 conv->AddEmModel(0, mod);
144 pmanager->AddDiscreteProcess(conv);
146 }
else if (particleName ==
"e-") {
148 G4eIonisation* eioni =
new G4eIonisation();
149 eioni->SetStepFunction(0.8, 1.0*mm);
150 G4eMultipleScattering* msc =
new G4eMultipleScattering;
151 msc->SetStepLimitType(fMinimal);
152 msc->AddEmModel(0,
new G4UrbanMscModel93());
154 G4eBremsstrahlung* ebrem =
new G4eBremsstrahlung();
156 pmanager->AddProcess(msc, -1, 1, 1);
157 pmanager->AddProcess(eioni, -1, 2, 2);
158 pmanager->AddProcess(ebrem, -1,-3, 3);
160 }
else if (particleName ==
"e+") {
162 G4eIonisation* eioni =
new G4eIonisation();
163 eioni->SetStepFunction(0.8, 1.0*mm);
164 G4eMultipleScattering* msc =
new G4eMultipleScattering;
165 msc->SetStepLimitType(fMinimal);
166 msc->AddEmModel(0,
new G4UrbanMscModel93());
168 G4eBremsstrahlung* ebrem =
new G4eBremsstrahlung();
170 pmanager->AddProcess(msc, -1, 1, 1);
171 pmanager->AddProcess(eioni, -1, 2, 2);
172 pmanager->AddProcess(ebrem, -1,-3, 3);
173 pmanager->AddProcess(
new G4eplusAnnihilation, 0,-1, 4);
175 }
else if (particleName ==
"mu+" ||
176 particleName ==
"mu-" ) {
178 G4MuMultipleScattering* mumsc =
new G4MuMultipleScattering();
179 mumsc->AddEmModel(0,
new G4WentzelVIModel());
180 pmanager->AddProcess(mumsc, -1, 1, 1);
181 pmanager->AddProcess(
new G4MuIonisation, -1, 2, 2);
182 pmanager->AddProcess(
new G4MuBremsstrahlung, -1,-3, 3);
183 pmanager->AddProcess(
new G4MuPairProduction, -1,-4, 4);
184 pmanager->AddProcess(
new G4CoulombScattering, -1,-4, 5);
186 }
else if (particleName ==
"alpha" ||
187 particleName ==
"He3" ||
188 particleName ==
"GenericIon") {
190 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
191 pmanager->AddProcess(
new G4ionIonisation, -1, 2, 2);
193 }
else if (particleName ==
"pi+" ||
194 particleName ==
"kaon+" ||
195 particleName ==
"kaon-" ||
196 particleName ==
"proton" ||
197 particleName ==
"pi-" ) {
199 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
200 pmanager->AddProcess(
new G4hIonisation, -1, 2, 2);
201 pmanager->AddProcess(
new G4hBremsstrahlung(), -1,-3, 3);
202 pmanager->AddProcess(
new G4hPairProduction(), -1,-4, 4);
204 }
else if (particleName ==
"B+" ||
205 particleName ==
"B-" ||
206 particleName ==
"D+" ||
207 particleName ==
"D-" ||
208 particleName ==
"Ds+" ||
209 particleName ==
"Ds-" ||
210 particleName ==
"anti_lambda_c+" ||
211 particleName ==
"anti_omega-" ||
212 particleName ==
"anti_proton" ||
213 particleName ==
"anti_sigma_c+" ||
214 particleName ==
"anti_sigma_c++" ||
215 particleName ==
"anti_sigma+" ||
216 particleName ==
"anti_sigma-" ||
217 particleName ==
"anti_xi_c+" ||
218 particleName ==
"anti_xi-" ||
219 particleName ==
"deuteron" ||
220 particleName ==
"lambda_c+" ||
221 particleName ==
"omega-" ||
222 particleName ==
"sigma_c+" ||
223 particleName ==
"sigma_c++" ||
224 particleName ==
"sigma+" ||
225 particleName ==
"sigma-" ||
226 particleName ==
"tau+" ||
227 particleName ==
"tau-" ||
228 particleName ==
"triton" ||
229 particleName ==
"xi_c+" ||
230 particleName ==
"xi-" ) {
232 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
233 pmanager->AddProcess(
new G4hIonisation, -1, 2, 2);
239 G4EmProcessOptions opt;
247 opt.SetApplyCuts(
true);
251 opt.SetMinEnergy(100.*eV);
252 opt.SetMaxEnergy(10.*TeV);
253 opt.SetDEDXBinning(77);
254 opt.SetLambdaBinning(77);
virtual ~CMSEmStandardPhysicsLPM()
static HepMC::IO_HEPEVT conv
CMSEmStandardPhysicsLPM(const G4String &name, G4int ver)
virtual void ConstructProcess()
virtual void ConstructParticle()
T mod(const T &a, const T &b)