126 theParticleIterator->reset();
127 while( (*theParticleIterator)() ){
128 G4ParticleDefinition* particle = theParticleIterator->value();
129 G4ProcessManager* pmanager = particle->GetProcessManager();
130 G4String particleName = particle->GetParticleName();
132 G4cout <<
"### " << GetPhysicsName() <<
" instantiates for "
133 << particleName <<
" at " << particle << G4endl;
135 if (particleName ==
"gamma") {
137 pmanager->AddDiscreteProcess(
new G4PhotoElectricEffect);
138 pmanager->AddDiscreteProcess(
new G4ComptonScattering);
139 pmanager->AddDiscreteProcess(
new G4GammaConversion);
141 }
else if (particleName ==
"e-") {
143 G4eIonisation* eioni =
new G4eIonisation();
144 eioni->SetStepFunction(0.8, 1.0*mm);
145 G4eMultipleScattering* msc =
new G4eMultipleScattering;
146 msc->SetStepLimitType(fMinimal);
147 msc->AddEmModel(0,
new CMSUrbanMscModel92());
149 G4eBremsstrahlung* ebrem =
new G4eBremsstrahlung();
150 ebrem->SetEmModel(
new G4SeltzerBergerModel95(), 1);
151 ebrem->SetEmModel(
new G4eBremsstrahlungRelModel95(), 2);
152 ebrem->EmModel(2)->SetLowEnergyLimit(GeV);
154 pmanager->AddProcess(msc, -1, 1, 1);
155 pmanager->AddProcess(eioni, -1, 2, 2);
156 pmanager->AddProcess(ebrem, -1,-3, 3);
158 }
else if (particleName ==
"e+") {
160 G4eIonisation* eioni =
new G4eIonisation();
161 eioni->SetStepFunction(0.8, 1.0*mm);
162 G4eMultipleScattering* msc =
new G4eMultipleScattering;
163 msc->SetStepLimitType(fMinimal);
164 msc->AddEmModel(0,
new CMSUrbanMscModel92());
166 G4eBremsstrahlung* ebrem =
new G4eBremsstrahlung();
167 ebrem->SetEmModel(
new G4SeltzerBergerModel95(), 1);
168 ebrem->SetEmModel(
new G4eBremsstrahlungRelModel95(), 2);
169 ebrem->EmModel(2)->SetLowEnergyLimit(GeV);
171 pmanager->AddProcess(msc, -1, 1, 1);
172 pmanager->AddProcess(eioni, -1, 2, 2);
173 pmanager->AddProcess(ebrem, -1,-3, 3);
174 pmanager->AddProcess(
new G4eplusAnnihilation, 0,-1, 4);
176 }
else if (particleName ==
"mu+" ||
177 particleName ==
"mu-" ) {
179 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
180 pmanager->AddProcess(
new G4MuIonisation, -1, 2, 2);
181 pmanager->AddProcess(
new G4MuBremsstrahlung, -1,-3, 3);
182 pmanager->AddProcess(
new G4MuPairProduction, -1,-4, 4);
184 }
else if (particleName ==
"alpha" ||
185 particleName ==
"He3" ||
186 particleName ==
"GenericIon") {
188 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
189 pmanager->AddProcess(
new G4ionIonisation, -1, 2, 2);
191 }
else if (particleName ==
"pi+" ||
192 particleName ==
"kaon+" ||
193 particleName ==
"kaon-" ||
194 particleName ==
"proton" ||
195 particleName ==
"pi-" ) {
197 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
198 pmanager->AddProcess(
new G4hIonisation, -1, 2, 2);
199 pmanager->AddProcess(
new G4hBremsstrahlung(), -1,-3, 3);
200 pmanager->AddProcess(
new G4hPairProduction(), -1,-4, 4);
202 }
else if (particleName ==
"B+" ||
203 particleName ==
"B-" ||
204 particleName ==
"D+" ||
205 particleName ==
"D-" ||
206 particleName ==
"Ds+" ||
207 particleName ==
"Ds-" ||
208 particleName ==
"anti_lambda_c+" ||
209 particleName ==
"anti_omega-" ||
210 particleName ==
"anti_proton" ||
211 particleName ==
"anti_sigma_c+" ||
212 particleName ==
"anti_sigma_c++" ||
213 particleName ==
"anti_sigma+" ||
214 particleName ==
"anti_sigma-" ||
215 particleName ==
"anti_xi_c+" ||
216 particleName ==
"anti_xi-" ||
217 particleName ==
"deuteron" ||
218 particleName ==
"lambda_c+" ||
219 particleName ==
"omega-" ||
220 particleName ==
"sigma_c+" ||
221 particleName ==
"sigma_c++" ||
222 particleName ==
"sigma+" ||
223 particleName ==
"sigma-" ||
224 particleName ==
"tau+" ||
225 particleName ==
"tau-" ||
226 particleName ==
"triton" ||
227 particleName ==
"xi_c+" ||
228 particleName ==
"xi-" ) {
230 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
231 pmanager->AddProcess(
new G4hIonisation, -1, 2, 2);
237 G4EmProcessOptions opt;
241 opt.SetApplyCuts(
true);
245 opt.SetMinEnergy(100.*eV);
246 opt.SetMaxEnergy(10.*TeV);
247 opt.SetDEDXBinning(77);
248 opt.SetLambdaBinning(77);