119 theParticleIterator->reset();
120 while( (*theParticleIterator)() ){
121 G4ParticleDefinition* particle = theParticleIterator->value();
122 G4ProcessManager* pmanager = particle->GetProcessManager();
123 G4String particleName = particle->GetParticleName();
125 G4cout <<
"### " << GetPhysicsName() <<
" instantiates for "
126 << particleName << G4endl;
128 if (particleName ==
"gamma") {
130 pmanager->AddDiscreteProcess(
new G4PhotoElectricEffect);
131 pmanager->AddDiscreteProcess(
new G4ComptonScattering);
132 pmanager->AddDiscreteProcess(
new G4GammaConversion);
134 }
else if (particleName ==
"e-") {
136 pmanager->AddProcess(
new G4eMultipleScattering,-1, 1, 1);
137 pmanager->AddProcess(
new G4eIonisation, -1, 2, 2);
138 pmanager->AddProcess(
new G4eBremsstrahlung(), -1,-3, 3);
140 }
else if (particleName ==
"e+") {
142 pmanager->AddProcess(
new G4eMultipleScattering,-1, 1, 1);
143 pmanager->AddProcess(
new G4eIonisation, -1, 2, 2);
144 pmanager->AddProcess(
new G4eBremsstrahlung, -1,-3, 3);
145 pmanager->AddProcess(
new G4eplusAnnihilation, 0,-1, 4);
147 }
else if (particleName ==
"mu+" ||
148 particleName ==
"mu-" ) {
150 pmanager->AddProcess(
new G4hMultipleScattering,-1, 1, 1);
151 pmanager->AddProcess(
new G4MuIonisation, -1, 2, 2);
152 pmanager->AddProcess(
new G4MuBremsstrahlung, -1,-3, 3);
153 pmanager->AddProcess(
new G4MuPairProduction, -1,-4, 4);
155 }
else if (particleName ==
"alpha" ||
156 particleName ==
"He3" ||
157 particleName ==
"GenericIon") {
159 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
160 pmanager->AddProcess(
new G4ionIonisation, -1, 2, 2);
163 }
else if (particleName ==
"B+" ||
164 particleName ==
"B-" ||
165 particleName ==
"D+" ||
166 particleName ==
"D-" ||
167 particleName ==
"Ds+" ||
168 particleName ==
"Ds-" ||
169 particleName ==
"anti_lambda_c+" ||
170 particleName ==
"anti_omega-" ||
171 particleName ==
"anti_proton" ||
172 particleName ==
"anti_sigma_c+" ||
173 particleName ==
"anti_sigma_c++" ||
174 particleName ==
"anti_sigma+" ||
175 particleName ==
"anti_sigma-" ||
176 particleName ==
"anti_xi_c+" ||
177 particleName ==
"anti_xi-" ||
178 particleName ==
"deuteron" ||
179 particleName ==
"kaon+" ||
180 particleName ==
"kaon-" ||
181 particleName ==
"lambda_c+" ||
182 particleName ==
"omega-" ||
183 particleName ==
"pi+" ||
184 particleName ==
"pi-" ||
185 particleName ==
"proton" ||
186 particleName ==
"sigma_c+" ||
187 particleName ==
"sigma_c++" ||
188 particleName ==
"sigma+" ||
189 particleName ==
"sigma-" ||
190 particleName ==
"tau+" ||
191 particleName ==
"tau-" ||
192 particleName ==
"triton" ||
193 particleName ==
"xi_c+" ||
194 particleName ==
"xi-" ) {
196 pmanager->AddProcess(
new G4hMultipleScattering, -1, 1, 1);
197 pmanager->AddProcess(
new G4hIonisation, -1, 2, 2);
199 if (particleName ==
"pi+" ||
200 particleName ==
"kaon+" ||
201 particleName ==
"kaon-" ||
202 particleName ==
"proton" ||
203 particleName ==
"pi-" ) {
204 pmanager->AddProcess(
new G4hBremsstrahlung(), -1,-3, 3);
205 pmanager->AddProcess(
new G4hPairProduction(), -1,-4, 4);
211 G4EmProcessOptions opt;