120 theParticleIterator->reset();
121 while( (*theParticleIterator)() ){
122 G4ParticleDefinition* particle = theParticleIterator->value();
123 G4ProcessManager* pmanager = particle->GetProcessManager();
124 G4String particleName = particle->GetParticleName();
126 G4cout <<
"### " << GetPhysicsName() <<
" instantiates for "
127 << particleName << G4endl;
129 if (particleName ==
"gamma") {
131 pmanager->AddDiscreteProcess(
new G4PhotoElectricEffect);
132 pmanager->AddDiscreteProcess(
new G4ComptonScattering);
133 pmanager->AddDiscreteProcess(
new G4GammaConversion);
135 }
else if (particleName ==
"e-") {
137 G4eIonisation* eioni =
new G4eIonisation();
138 eioni->SetStepFunction(0.8, 1.0*mm);
139 G4eMultipleScattering* msc =
new G4eMultipleScattering;
140 msc->SetStepLimitType(fMinimal);
141 pmanager->AddProcess(msc, -1, 1, 1);
142 pmanager->AddProcess(eioni, -1, 2, 2);
143 pmanager->AddProcess(
new G4eBremsstrahlung, -1,-3, 3);
145 }
else if (particleName ==
"e+") {
147 G4eIonisation* eioni =
new G4eIonisation();
148 eioni->SetStepFunction(0.8, 1.0*mm);
149 G4eMultipleScattering* msc =
new G4eMultipleScattering;
150 msc->SetStepLimitType(fMinimal);
151 pmanager->AddProcess(msc, -1, 1, 1);
152 pmanager->AddProcess(eioni, -1, 2, 2);
153 pmanager->AddProcess(
new G4eBremsstrahlung, -1,-3, 3);
154 pmanager->AddProcess(
new G4eplusAnnihilation, 0,-1, 4);
156 }
else if (particleName ==
"mu+" ||
157 particleName ==
"mu-" ) {
159 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
160 pmanager->AddProcess(
new G4MuIonisation, -1, 2, 2);
161 pmanager->AddProcess(
new G4MuBremsstrahlung, -1,-3, 3);
162 pmanager->AddProcess(
new G4MuPairProduction, -1,-4, 4);
164 }
else if (particleName ==
"alpha" ||
165 particleName ==
"He3" ||
166 particleName ==
"GenericIon") {
168 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
169 pmanager->AddProcess(
new G4ionIonisation, -1, 2, 2);
171 }
else if (particleName ==
"B+" ||
172 particleName ==
"B-" ||
173 particleName ==
"D+" ||
174 particleName ==
"D-" ||
175 particleName ==
"Ds+" ||
176 particleName ==
"Ds-" ||
177 particleName ==
"anti_lambda_c+" ||
178 particleName ==
"anti_omega-" ||
179 particleName ==
"anti_proton" ||
180 particleName ==
"anti_sigma_c+" ||
181 particleName ==
"anti_sigma_c++" ||
182 particleName ==
"anti_sigma+" ||
183 particleName ==
"anti_sigma-" ||
184 particleName ==
"anti_xi_c+" ||
185 particleName ==
"anti_xi-" ||
186 particleName ==
"deuteron" ||
187 particleName ==
"kaon+" ||
188 particleName ==
"kaon-" ||
189 particleName ==
"lambda_c+" ||
190 particleName ==
"omega-" ||
191 particleName ==
"pi+" ||
192 particleName ==
"pi-" ||
193 particleName ==
"proton" ||
194 particleName ==
"sigma_c+" ||
195 particleName ==
"sigma_c++" ||
196 particleName ==
"sigma+" ||
197 particleName ==
"sigma-" ||
198 particleName ==
"tau+" ||
199 particleName ==
"tau-" ||
200 particleName ==
"triton" ||
201 particleName ==
"xi_c+" ||
202 particleName ==
"xi-" ) {
204 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
205 pmanager->AddProcess(
new G4hIonisation, -1, 2, 2);
207 if (particleName ==
"pi+" ||
208 particleName ==
"kaon+" ||
209 particleName ==
"kaon-" ||
210 particleName ==
"proton" ||
211 particleName ==
"pi-" ) {
212 pmanager->AddProcess(
new G4hBremsstrahlung(), -1,-3, 3);
213 pmanager->AddProcess(
new G4hPairProduction(), -1,-4, 4);
219 G4EmProcessOptions opt;