Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4AdjointPhysicsList.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
28 //
29 // $Id: G4AdjointPhysicsList.cc 101303 2016-11-14 11:21:18Z gcosmo $
30 //
32 // Class Name: G4AdjointPhysicsList
33 // Author: L. Desorgher
34 // Organisation: SpaceIT GmbH
35 // Contract: ESA contract 21435/08/NL/AT
36 // Customer: ESA/ESTEC
38 
39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
41 
42 #include "G4AdjointPhysicsList.hh"
43 #include "G4ProcessManager.hh"
44 #include "G4ParticleTypes.hh"
45 #include "G4AdjointPhysicsMessenger.hh"
46 #include "G4SystemOfUnits.hh"
47 
48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49 
50 G4AdjointPhysicsList::G4AdjointPhysicsList()
51  :G4VUserPhysicsList(),
52  fEminusIonisation(0),fPIonisation(0),
53  fUse_forced_interaction(true),
54  fUse_eionisation(true),fUse_pionisation(true),
55  fUse_brem(true),fUse_compton(true),fUse_ms(true),
56  fUse_egain_fluctuation(true),fUse_peeffect(true),
57  fEmin_adj_models(1.*keV), fEmax_adj_models(1.*MeV),
58  fCS_biasing_factor_compton(1.),fCS_biasing_factor_brem(1.),
59  fCS_biasing_factor_ionisation(1.),fCS_biasing_factor_PEeffect(1.)
60 {
61  defaultCutValue = 1.0*mm;
62  SetVerboseLevel(1);
63  fPhysicsMessenger = new G4AdjointPhysicsMessenger(this);
64 }
65 
66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
67 
68 G4AdjointPhysicsList::~G4AdjointPhysicsList()
69 {
70 }
71 void G4AdjointPhysicsList::ConstructParticle()
72 {
73  // In this method, static member functions should be called
74  // for all particles which you want to use.
75  // This ensures that objects of these particle types will be
76  // created in the program.
77  ConstructBosons();
78  ConstructLeptons();
79  ConstructMesons();
80  ConstructBaryons();
81  ConstructAdjointParticles();
82 }
83 
84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
85 
86 void G4AdjointPhysicsList::SetLossFluctuationFlag(bool aBool)
87 {
88  if (fEminusIonisation) fEminusIonisation->SetLossFluctuations(aBool);
89 }
90 
91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
92 
93 void G4AdjointPhysicsList::ConstructBosons()
94 {
95  // pseudo-particles
96  G4Geantino::GeantinoDefinition();
97  G4ChargedGeantino::ChargedGeantinoDefinition();
98 
99  // gamma
100  G4Gamma::GammaDefinition();
101 
102  // optical photon
103  G4OpticalPhoton::OpticalPhotonDefinition();
104 }
105 
106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
107 
108 void G4AdjointPhysicsList::ConstructLeptons()
109 {
110  // leptons
111  G4Electron::ElectronDefinition();
112  G4Positron::PositronDefinition();
113  G4MuonPlus::MuonPlusDefinition();
114  G4MuonMinus::MuonMinusDefinition();
115 
116  G4NeutrinoE::NeutrinoEDefinition();
117  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
118  G4NeutrinoMu::NeutrinoMuDefinition();
119  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
120 }
121 
122 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
123 
124 void G4AdjointPhysicsList::ConstructMesons()
125 {
126 // mesons
127  G4PionPlus::PionPlusDefinition();
128  G4PionMinus::PionMinusDefinition();
129  G4PionZero::PionZeroDefinition();
130  G4Eta::EtaDefinition();
131  G4EtaPrime::EtaPrimeDefinition();
132  G4KaonPlus::KaonPlusDefinition();
133  G4KaonMinus::KaonMinusDefinition();
134  G4KaonZero::KaonZeroDefinition();
135  G4AntiKaonZero::AntiKaonZeroDefinition();
136  G4KaonZeroLong::KaonZeroLongDefinition();
137  G4KaonZeroShort::KaonZeroShortDefinition();
138 }
139 
140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141 
142 void G4AdjointPhysicsList::ConstructBaryons()
143 {
144 // barions
145  G4Proton::ProtonDefinition();
146  G4AntiProton::AntiProtonDefinition();
147  G4Neutron::NeutronDefinition();
148  G4AntiNeutron::AntiNeutronDefinition();
149 }
150 
151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
152 
153 #include"G4AdjointGamma.hh"
154 #include"G4AdjointElectron.hh"
155 #include"G4AdjointProton.hh"
156 void G4AdjointPhysicsList::ConstructAdjointParticles()
157 {
158 // adjoint_gammma
159  G4AdjointGamma::AdjointGammaDefinition();
160 
161 // adjoint_electron
162  G4AdjointElectron::AdjointElectronDefinition();
163 
164 // adjoint_proton
165  G4AdjointProton::AdjointProtonDefinition();
166 }
167 
168 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
169 
170 void G4AdjointPhysicsList::ConstructProcess()
171 {
172  AddTransportation();
173  ConstructEM();
174  ConstructGeneral();
175 }
176 
177 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
178 
179 //#include "G4PEEffectFluoModel.hh"
180 #include "G4ComptonScattering.hh"
181 #include "G4GammaConversion.hh"
182 #include "G4PhotoElectricEffect.hh"
183 #include "G4eMultipleScattering.hh"
184 #include "G4eAdjointMultipleScattering.hh"
185 #include "G4hMultipleScattering.hh"
186 #include "G4eIonisation.hh"
187 #include "G4eBremsstrahlung.hh"
188 #include "G4eplusAnnihilation.hh"
189 #include "G4hIonisation.hh"
190 #include "G4ionIonisation.hh"
191 //#include "G4IonParametrisedLossModel.hh"
192 
193 #include "G4eBremsstrahlung.hh"
194 #include "G4ContinuousGainOfEnergy.hh"
195 #include "G4eInverseIonisation.hh"
196 #include "G4AdjointeIonisationModel.hh"
197 #include "G4AdjointCSManager.hh"
198 #include "G4AdjointBremsstrahlungModel.hh"
199 #include "G4eInverseBremsstrahlung.hh"
200 #include "G4AdjointComptonModel.hh"
201 #include "G4eInverseCompton.hh"
202 #include "G4InversePEEffect.hh"
203 #include "G4AdjointPhotoElectricModel.hh"
204 #include "G4AdjointAlongStepWeightCorrection.hh"
205 #include "G4hInverseIonisation.hh"
206 #include "G4AdjointhIonisationModel.hh"
207 #include "G4AdjointhMultipleScattering.hh"
208 #include "G4IonInverseIonisation.hh"
209 #include "G4AdjointIonIonisationModel.hh"
210 
211 #include "G4AdjointSimManager.hh"
212 #include "G4AdjointProcessEquivalentToDirectProcess.hh"
213 
214 #include "G4SystemOfUnits.hh"
215 #include "G4PhysicalConstants.hh"
216 #include "G4UrbanMscModel.hh"
217 #include "G4UrbanAdjointMscModel.hh"
218 #include "G4UrbanAdjointMscModel.hh"
219 #include "G4AdjointForcedInteractionForGamma.hh"
220 
221 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
222 
223 void G4AdjointPhysicsList::ConstructEM()
224 {
225  G4AdjointCSManager* theCSManager =
226  G4AdjointCSManager::GetAdjointCSManager();
227  G4AdjointSimManager* theAdjointSimManager =
228  G4AdjointSimManager::GetInstance();
229 
230  theCSManager->RegisterAdjointParticle(
231  G4AdjointElectron::AdjointElectron());
232 
233  if (fUse_brem || fUse_peeffect ||fUse_compton)
234  theCSManager->RegisterAdjointParticle(
235  G4AdjointGamma::AdjointGamma());
236 
237  if (fUse_eionisation) {
238  if (!fEminusIonisation) fEminusIonisation =
239  new G4eIonisation();
240  fEminusIonisation->SetLossFluctuations(fUse_egain_fluctuation);
241  }
242  if (fUse_pionisation) {
243  if (!fPIonisation) fPIonisation = new G4hIonisation();
244  fPIonisation->SetLossFluctuations(fUse_egain_fluctuation);
245  theCSManager->RegisterAdjointParticle(
246  G4AdjointProton::AdjointProton());
247  }
248 
249  G4eBremsstrahlung* theeminusBremsstrahlung = 0;
250  if (fUse_brem && fUse_eionisation)
251  theeminusBremsstrahlung = new G4eBremsstrahlung();
252 
253  G4ComptonScattering* theComptonScattering =0;
254  if (fUse_compton) theComptonScattering = new G4ComptonScattering();
255 
256  G4PhotoElectricEffect* thePEEffect =0;
257  if (fUse_peeffect) thePEEffect = new G4PhotoElectricEffect();
258 
259  G4eMultipleScattering* theeminusMS = 0;
260  G4hMultipleScattering* thepMS= 0;
261  G4eAdjointMultipleScattering* theeminusAdjointMS = 0;
262  if (fUse_ms) {
263  theeminusMS = new G4eMultipleScattering();
264  G4UrbanMscModel* msc1 = new G4UrbanMscModel();
265  msc1->SetNewDisplacementFlag(false);
266  theeminusMS->AddEmModel(0, msc1);
267  theeminusAdjointMS = new G4eAdjointMultipleScattering();
268  G4UrbanAdjointMscModel* msc2 = new G4UrbanAdjointMscModel();
269  msc2->SetNewDisplacementFlag(false);
270  theeminusAdjointMS->AddEmModel(0, msc2);
271  thepMS = new G4hMultipleScattering();
272  }
273 
274  G4VProcess* theGammaConversion =0;
275  if (fUse_gamma_conversion) theGammaConversion = new G4GammaConversion();
276  //Define adjoint e- ionisation
277  //-------------------
278  G4AdjointeIonisationModel* theeInverseIonisationModel = 0;
279  G4eInverseIonisation* theeInverseIonisationProjToProjCase = 0 ;
280  G4eInverseIonisation* theeInverseIonisationProdToProjCase = 0;
281  if (fUse_eionisation) {
282  theeInverseIonisationModel = new G4AdjointeIonisationModel();
283  theeInverseIonisationModel->SetHighEnergyLimit(
284  fEmax_adj_models);
285  theeInverseIonisationModel->SetLowEnergyLimit(
286  fEmin_adj_models);
287  theeInverseIonisationModel->SetCSBiasingFactor(
288  fCS_biasing_factor_ionisation);
289  theeInverseIonisationProjToProjCase =
290  new G4eInverseIonisation(true,"Inv_eIon",
291  theeInverseIonisationModel);
292  theeInverseIonisationProdToProjCase =
293  new G4eInverseIonisation(false,"Inv_eIon1",
294  theeInverseIonisationModel);
295  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
296  }
297 
298  //Define adjoint Bremsstrahlung
299  //-------------------------------
300  G4AdjointBremsstrahlungModel* theeInverseBremsstrahlungModel = 0;
301  G4eInverseBremsstrahlung* theeInverseBremsstrahlungProjToProjCase = 0;
302  G4eInverseBremsstrahlung* theeInverseBremsstrahlungProdToProjCase = 0;
303  G4AdjointForcedInteractionForGamma* theForcedInteractionForGamma = 0;
304  if (fUse_brem && fUse_eionisation) {
305  theeInverseBremsstrahlungModel = new G4AdjointBremsstrahlungModel();
306  theeInverseBremsstrahlungModel->SetHighEnergyLimit(fEmax_adj_models*1.01);
307  theeInverseBremsstrahlungModel->SetLowEnergyLimit(fEmin_adj_models);
308  theeInverseBremsstrahlungModel->SetCSBiasingFactor(
309  fCS_biasing_factor_brem);
310  theeInverseBremsstrahlungProjToProjCase = new G4eInverseBremsstrahlung(
311  true,"Inv_eBrem",theeInverseBremsstrahlungModel);
312  theeInverseBremsstrahlungProdToProjCase = new G4eInverseBremsstrahlung(
313  false,"Inv_eBrem1",theeInverseBremsstrahlungModel);
314  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
315  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("gamma"));
316 
317  if (!fUse_forced_interaction) theeInverseBremsstrahlungProdToProjCase
318  = new G4eInverseBremsstrahlung(false,G4String("Inv_eBrem1"),
319  theeInverseBremsstrahlungModel);
320  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
321  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("gamma"));
322  if (fUse_forced_interaction){
323  theForcedInteractionForGamma =
324  new G4AdjointForcedInteractionForGamma("ReverseGammaForcedInteraction");
325  theForcedInteractionForGamma->RegisterAdjointBremModel(
326  theeInverseBremsstrahlungModel);
327  }
328  }
329 
330 
331  //Define adjoint Compton
332  //---------------------
333 
334  G4AdjointComptonModel* theeInverseComptonModel = 0;
335  G4eInverseCompton* theeInverseComptonProjToProjCase = 0;
336  G4eInverseCompton* theeInverseComptonProdToProjCase = 0;
337 
338  if (fUse_compton) {
339  theeInverseComptonModel = new G4AdjointComptonModel();
340  theeInverseComptonModel->SetHighEnergyLimit(fEmax_adj_models);
341  theeInverseComptonModel->SetLowEnergyLimit(fEmin_adj_models);
342  theeInverseComptonModel->SetDirectProcess(theComptonScattering);
343  theeInverseComptonModel->SetUseMatrix(false);
344 
345  theeInverseComptonModel->SetCSBiasingFactor( fCS_biasing_factor_compton);
346  if (!fUse_forced_interaction) theeInverseComptonProjToProjCase =
347  new G4eInverseCompton(true,"Inv_Compt",theeInverseComptonModel);
348  theeInverseComptonProdToProjCase = new G4eInverseCompton(false,"Inv_Compt1",
349  theeInverseComptonModel);
350  if (fUse_forced_interaction){
351  if (!theForcedInteractionForGamma ) theForcedInteractionForGamma =
352  new G4AdjointForcedInteractionForGamma("ReverseGammaForcedInteraction");
353  theForcedInteractionForGamma->
354  RegisterAdjointComptonModel(theeInverseComptonModel);
355  }
356  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
357  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("gamma"));
358  }
359 
360  //Define adjoint PEEffect
361  //---------------------
362  G4AdjointPhotoElectricModel* theInversePhotoElectricModel = 0;
363  G4InversePEEffect* theInversePhotoElectricProcess = 0;
364 
365  if (fUse_peeffect) {
366  theInversePhotoElectricModel = new G4AdjointPhotoElectricModel();
367  theInversePhotoElectricModel->SetHighEnergyLimit(fEmax_adj_models);
368  theInversePhotoElectricModel->SetLowEnergyLimit(fEmin_adj_models);
369  theInversePhotoElectricModel->SetCSBiasingFactor(
370  fCS_biasing_factor_PEeffect);
371  theInversePhotoElectricProcess = new G4InversePEEffect("Inv_PEEffect",
372  theInversePhotoElectricModel);
373  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
374  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("gamma"));
375  }
376 
377 
378  //Define adjoint ionisation for protons
379  //---------------------
380  G4AdjointhIonisationModel* thepInverseIonisationModel = 0;
381  G4hInverseIonisation* thepInverseIonisationProjToProjCase = 0 ;
382  G4hInverseIonisation* thepInverseIonisationProdToProjCase = 0;
383  if (fUse_pionisation) {
384  thepInverseIonisationModel = new G4AdjointhIonisationModel(
385  G4Proton::Proton());
386  thepInverseIonisationModel->SetHighEnergyLimit(fEmax_adj_models);
387  thepInverseIonisationModel->SetLowEnergyLimit(fEmin_adj_models);
388  thepInverseIonisationModel->SetUseMatrix(false);
389  thepInverseIonisationProjToProjCase = new G4hInverseIonisation(true,
390  "Inv_pIon",thepInverseIonisationModel);
391  thepInverseIonisationProdToProjCase = new G4hInverseIonisation(false,
392  "Inv_pIon1",thepInverseIonisationModel);
393  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
394  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("proton"));
395  }
396 
397  //Declare the processes active for the different particles
398  //--------------------------------------------------------
399  auto particleIterator=GetParticleIterator();
400  particleIterator->reset();
401  while( (*particleIterator)() ){
402  G4ParticleDefinition* particle = particleIterator->value();
403  G4ProcessManager* pmanager = particle->GetProcessManager();
404  if (!pmanager) {
405  pmanager = new G4ProcessManager(particle);
406  particle->SetProcessManager(pmanager);
407  }
408 
409  G4String particleName = particle->GetParticleName();
410  if (particleName == "e-") {
411  if (fUse_ms && fUse_eionisation) pmanager->AddProcess(theeminusMS);
412  if (fUse_eionisation){
413  pmanager->AddProcess(fEminusIonisation);
414  G4AdjointCSManager::GetAdjointCSManager()->
415  RegisterEnergyLossProcess(fEminusIonisation,particle);
416  }
417  if (fUse_brem && fUse_eionisation) {
418  pmanager->AddProcess(theeminusBremsstrahlung);
419  G4AdjointCSManager::GetAdjointCSManager()->
420  RegisterEnergyLossProcess(theeminusBremsstrahlung,particle);
421  }
422  G4int n_order=0;
423  if (fUse_ms && fUse_eionisation) {
424  n_order++;
425  pmanager->SetProcessOrdering(theeminusMS, idxAlongStep,n_order);
426  }
427  if (fUse_eionisation) {
428  n_order++;
429  pmanager->SetProcessOrdering(fEminusIonisation,idxAlongStep,n_order);
430  }
431  if (fUse_brem && fUse_eionisation) {
432  n_order++;
433  pmanager->SetProcessOrdering(theeminusBremsstrahlung,
434  idxAlongStep,n_order);
435  }
436  n_order=0;
437  if (fUse_ms && fUse_eionisation) {
438  n_order++;
439  pmanager->SetProcessOrdering(theeminusMS,idxPostStep,n_order);
440  }
441  if (fUse_eionisation) {
442  n_order++;
443  pmanager->SetProcessOrdering(fEminusIonisation,idxPostStep,n_order);
444  }
445  if (fUse_brem && fUse_eionisation) {
446  n_order++;
447  pmanager->SetProcessOrdering(theeminusBremsstrahlung,idxPostStep,
448  n_order);
449  }
450  }
451 
452  if (particleName == "adj_e-") {
453  G4ContinuousGainOfEnergy* theContinuousGainOfEnergy =0;
454  if (fUse_eionisation ) {
455  theContinuousGainOfEnergy= new G4ContinuousGainOfEnergy();
456  theContinuousGainOfEnergy->SetLossFluctuations(
457  fUse_egain_fluctuation);
458  theContinuousGainOfEnergy->SetDirectEnergyLossProcess(
459  fEminusIonisation);
460  theContinuousGainOfEnergy->SetDirectParticle(G4Electron::Electron());
461  pmanager->AddProcess(theContinuousGainOfEnergy);
462  }
463  G4int n_order=0;
464  if (fUse_ms) {
465  n_order++;
466  pmanager->AddProcess(theeminusAdjointMS);
467  pmanager->SetProcessOrdering(theeminusAdjointMS,
468  idxAlongStep,n_order);
469  }
470  n_order++;
471  pmanager->SetProcessOrdering(theContinuousGainOfEnergy,idxAlongStep,
472  n_order);
473 
474  n_order++;
475  G4AdjointAlongStepWeightCorrection* theAlongStepWeightCorrection =
476  new G4AdjointAlongStepWeightCorrection();
477  pmanager->AddProcess(theAlongStepWeightCorrection);
478  pmanager->SetProcessOrdering(theAlongStepWeightCorrection,
479  idxAlongStep,
480  n_order);
481  n_order=0;
482  if (fUse_eionisation) {
483  pmanager->AddProcess(theeInverseIonisationProjToProjCase);
484  pmanager->AddProcess(theeInverseIonisationProdToProjCase);
485  n_order++;
486  pmanager->SetProcessOrdering(theeInverseIonisationProjToProjCase,
487  idxPostStep,n_order);
488  n_order++;
489  pmanager->SetProcessOrdering(theeInverseIonisationProdToProjCase,
490  idxPostStep,n_order);
491  }
492  if (fUse_brem && fUse_eionisation) {
493  pmanager->AddProcess(theeInverseBremsstrahlungProjToProjCase);
494  n_order++;
495  pmanager->SetProcessOrdering(
496  theeInverseBremsstrahlungProjToProjCase,
497  idxPostStep,n_order);
498  }
499 
500  if (fUse_compton) {
501  pmanager->AddProcess(theeInverseComptonProdToProjCase);
502  n_order++;
503  pmanager->SetProcessOrdering(theeInverseComptonProdToProjCase,
504  idxPostStep,n_order);
505  }
506  if (fUse_peeffect) {
507  pmanager->AddDiscreteProcess(theInversePhotoElectricProcess);
508  n_order++;
509  pmanager->SetProcessOrdering(theInversePhotoElectricProcess,
510  idxPostStep,n_order);
511  }
512  if (fUse_pionisation) {
513  pmanager->AddProcess(thepInverseIonisationProdToProjCase);
514  n_order++;
515  pmanager->SetProcessOrdering(thepInverseIonisationProdToProjCase,
516  idxPostStep,n_order);
517  }
518  if (fUse_ms && fUse_eionisation) {
519  n_order++;
520  pmanager->SetProcessOrdering(theeminusAdjointMS,
521  idxPostStep,n_order);
522  }
523  }
524 
525 
526  if(particleName == "adj_gamma") {
527  G4int n_order=0;
528  if (!fUse_forced_interaction){
529  G4AdjointAlongStepWeightCorrection* theAlongStepWeightCorrection =
530  new G4AdjointAlongStepWeightCorrection();
531  pmanager->AddProcess(theAlongStepWeightCorrection);
532  pmanager->SetProcessOrdering(theAlongStepWeightCorrection,
533  idxAlongStep,1);
534 
535  if (fUse_brem && fUse_eionisation) {
536  pmanager->AddProcess(theeInverseBremsstrahlungProdToProjCase);
537  n_order++;
538  pmanager->SetProcessOrdering(
539  theeInverseBremsstrahlungProdToProjCase,
540  idxPostStep,n_order);
541  }
542  if (fUse_compton) {
543  pmanager->AddDiscreteProcess(theeInverseComptonProjToProjCase);
544  n_order++;
545  pmanager->SetProcessOrdering(theeInverseComptonProjToProjCase,
546  idxPostStep,n_order);
547  }
548  }
549  else {
550  if (theForcedInteractionForGamma) {
551  pmanager->AddProcess(theForcedInteractionForGamma);
552  n_order++;
553  pmanager->SetProcessOrdering(theForcedInteractionForGamma,
554  idxPostStep,n_order);
555  pmanager->SetProcessOrdering(theForcedInteractionForGamma,
556  idxAlongStep,n_order);
557  }
558  }
559  }
560 
561  if (particleName == "gamma") {
562  if (fUse_compton) {
563  pmanager->AddDiscreteProcess(theComptonScattering);
564  G4AdjointCSManager::GetAdjointCSManager()->
565  RegisterEmProcess(theComptonScattering,particle);
566  }
567  if (fUse_peeffect) {
568  pmanager->AddDiscreteProcess(thePEEffect);
569  G4AdjointCSManager::GetAdjointCSManager()->
570  RegisterEmProcess(thePEEffect,particle);
571  }
572  if (fUse_gamma_conversion) {
573  pmanager->AddDiscreteProcess(theGammaConversion);
574  }
575  }
576 
577  if (particleName == "e+" && fUse_gamma_conversion) {//positron
578  G4VProcess* theeplusMultipleScattering = new G4eMultipleScattering();
579  G4VProcess* theeplusIonisation = new G4eIonisation();
580  G4VProcess* theeplusBremsstrahlung = new G4eBremsstrahlung();
581  G4VProcess* theeplusAnnihilation = new G4eplusAnnihilation();
582 
583  // add processes
584  pmanager->AddProcess(theeplusMultipleScattering);
585  pmanager->AddProcess(theeplusIonisation);
586  pmanager->AddProcess(theeplusBremsstrahlung);
587  pmanager->AddProcess(theeplusAnnihilation);
588 
589  // set ordering for AtRestDoIt
590  pmanager->SetProcessOrderingToFirst(theeplusAnnihilation, idxAtRest);
591 
592  // set ordering for AlongStepDoIt
593  pmanager->SetProcessOrdering(theeplusMultipleScattering,
594  idxAlongStep,1);
595  pmanager->SetProcessOrdering(theeplusIonisation, idxAlongStep,2);
596  pmanager->SetProcessOrdering(theeplusBremsstrahlung,idxAlongStep,3);
597 
598  // set ordering for PostStepDoIt
599  pmanager->SetProcessOrdering(theeplusMultipleScattering,
600  idxPostStep,1);
601  pmanager->SetProcessOrdering(theeplusIonisation,idxPostStep,2);
602  pmanager->SetProcessOrdering(theeplusBremsstrahlung,idxPostStep,3);
603  pmanager->SetProcessOrdering(theeplusAnnihilation,idxPostStep,4);
604  }
605  if (particleName == "proton" && fUse_pionisation) {
606  if (fUse_ms && fUse_pionisation) pmanager->AddProcess(thepMS);
607 
608  if (fUse_pionisation){
609  pmanager->AddProcess(fPIonisation);
610  G4AdjointCSManager::GetAdjointCSManager()->
611  RegisterEnergyLossProcess(fPIonisation,particle);
612  }
613 
614  G4int n_order=0;
615  if (fUse_ms && fUse_pionisation) {
616  n_order++;
617  pmanager->SetProcessOrdering(thepMS, idxAlongStep,n_order);
618  }
619 
620  if (fUse_pionisation) {
621  n_order++;
622  pmanager->SetProcessOrdering(fPIonisation,idxAlongStep,n_order);
623  }
624 
625  n_order=0;
626  if (fUse_ms && fUse_pionisation) {
627  n_order++;
628  pmanager->SetProcessOrdering(thepMS, idxPostStep,n_order);
629  }
630 
631  if (fUse_pionisation) {
632  n_order++;
633  pmanager->SetProcessOrdering(fPIonisation,idxPostStep,n_order);
634  }
635 
636  }
637 
638  if (particleName == "adj_proton" && fUse_pionisation) {
639  G4ContinuousGainOfEnergy* theContinuousGainOfEnergy =0;
640  if (fUse_pionisation ) {
641  theContinuousGainOfEnergy= new G4ContinuousGainOfEnergy();
642  theContinuousGainOfEnergy->SetLossFluctuations(
643  fUse_egain_fluctuation);
644  theContinuousGainOfEnergy->SetDirectEnergyLossProcess(fPIonisation);
645  theContinuousGainOfEnergy->SetDirectParticle(G4Proton::Proton());
646  pmanager->AddProcess(theContinuousGainOfEnergy);
647  }
648 
649  G4int n_order=0;
650  if (fUse_ms) {
651  n_order++;
652  pmanager->AddProcess(thepMS);
653  pmanager->SetProcessOrdering(thepMS, idxAlongStep,n_order);
654  }
655 
656  n_order++;
657  pmanager->SetProcessOrdering(theContinuousGainOfEnergy,idxAlongStep,
658  n_order);
659 
660  n_order++;
661  G4AdjointAlongStepWeightCorrection* theAlongStepWeightCorrection =
662  new G4AdjointAlongStepWeightCorrection();
663  pmanager->AddProcess(theAlongStepWeightCorrection);
664  pmanager->SetProcessOrdering(theAlongStepWeightCorrection,
665  idxAlongStep,
666  n_order);
667  n_order=0;
668  if (fUse_pionisation) {
669  pmanager->AddProcess(thepInverseIonisationProjToProjCase);
670  n_order++;
671  pmanager->SetProcessOrdering(
672  thepInverseIonisationProjToProjCase,
673  idxPostStep,n_order);
674  }
675 
676  if (fUse_ms && fUse_pionisation) {
677  n_order++;
678  pmanager->SetProcessOrdering(thepMS,idxPostStep,n_order);
679  }
680  }
681  }
682 }
683 
684 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
685 
686 #include "G4Decay.hh"
687 void G4AdjointPhysicsList::ConstructGeneral()
688 {
689  // Add Decay Process
690  G4Decay* theDecayProcess = new G4Decay();
691  auto particleIterator=GetParticleIterator();
692  particleIterator->reset();
693  while( (*particleIterator)() ){
694  G4ParticleDefinition* particle = particleIterator->value();
695  G4ProcessManager* pmanager = particle->GetProcessManager();
696  if (theDecayProcess->IsApplicable(*particle)) {
697  pmanager ->AddProcess(theDecayProcess);
698  // set ordering for PostStepDoIt and AtRestDoIt
699  pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
700  pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
701  }
702  }
703 }
704 
705 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
706 
707 void G4AdjointPhysicsList::SetCuts()
708 {
709  if (verboseLevel >0){
710  G4cout << "G4AdjointPhysicsList::SetCuts:";
711  G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
712  }
713 
714  // set cut values for gamma at first and for e- second and next for e+,
715  // because some processes for e+/e- need cut values for gamma
716  //
717  SetCutValue(defaultCutValue, "gamma");
718  SetCutValue(defaultCutValue, "e-");
719  SetCutValue(defaultCutValue, "e+");
720 
721  if (verboseLevel>0) DumpCutValuesTable();
722 }
723 
724 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
G4hMultipleScattering.hh
G4Geantino::GeantinoDefinition
static G4Geantino * GeantinoDefinition()
Definition: G4Geantino.cc:82
G4AdjointSimManager::GetInstance
static G4AdjointSimManager * GetInstance()
Definition: G4AdjointSimManager.cc:134
G4AdjointGamma::AdjointGamma
static G4AdjointGamma * AdjointGamma()
Definition: G4AdjointGamma.cc:84
G4Electron::ElectronDefinition
static G4Electron * ElectronDefinition()
Definition: G4Electron.cc:89
G4KaonPlus::KaonPlusDefinition
static G4KaonPlus * KaonPlusDefinition()
Definition: G4KaonPlus.cc:108
G4eIonisation.hh
G4UrbanAdjointMscModel.hh
G4AdjointGamma::AdjointGammaDefinition
static G4AdjointGamma * AdjointGammaDefinition()
Definition: G4AdjointGamma.cc:79
G4MuonPlus::MuonPlusDefinition
static G4MuonPlus * MuonPlusDefinition()
Definition: G4MuonPlus.cc:94
G4IonInverseIonisation.hh
mm
static constexpr double mm
Definition: G4SIunits.hh:115
G4ParticleDefinition::SetProcessManager
void SetProcessManager(G4ProcessManager *aProcessManager)
Definition: G4ParticleDefinition.cc:450
G4VUserPhysicsList::SetCutValue
void SetCutValue(G4double aCut, const G4String &pname)
Definition: G4VUserPhysicsList.cc:416
G4ProcessManager::SetProcessOrderingToFirst
void SetProcessOrderingToFirst(G4VProcess *aProcess, G4ProcessVectorDoItIndex idDoIt)
Definition: G4ProcessManager.cc:688
G4ComptonScattering.hh
G4Decay::IsApplicable
virtual G4bool IsApplicable(const G4ParticleDefinition &)
Definition: G4Decay.cc:89
G4AdjointComptonModel.hh
G4AdjointPhysicsList::ConstructParticle
virtual void ConstructParticle()
Definition: G4AdjointPhysicsList.cc:71
G4hIonisation.hh
G4KaonMinus::KaonMinusDefinition
static G4KaonMinus * KaonMinusDefinition()
Definition: G4KaonMinus.cc:108
G4UrbanMscModel::SetNewDisplacementFlag
void SetNewDisplacementFlag(G4bool)
Definition: G4UrbanMscModel.hh:195
G4eAdjointMultipleScattering.hh
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4ProcessManager::AddDiscreteProcess
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
G4eMultipleScattering.hh
G4AdjointhMultipleScattering.hh
G4KaonZero::KaonZeroDefinition
static G4KaonZero * KaonZeroDefinition()
Definition: G4KaonZero.cc:99
G4AntiKaonZero::AntiKaonZeroDefinition
static G4AntiKaonZero * AntiKaonZeroDefinition()
Definition: G4AntiKaonZero.cc:98
G4AdjointElectron::AdjointElectron
static G4AdjointElectron * AdjointElectron()
Definition: G4AdjointElectron.cc:100
G4KaonZeroShort::KaonZeroShortDefinition
static G4KaonZeroShort * KaonZeroShortDefinition()
Definition: G4KaonZeroShort.cc:99
G4BestUnit
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
Definition: G4SteppingVerbose.cc:53
G4AdjointhIonisationModel.hh
G4AntiProton::AntiProtonDefinition
static G4AntiProton * AntiProtonDefinition()
Definition: G4AntiProton.cc:88
G4VEmAdjointModel::SetLowEnergyLimit
void SetLowEnergyLimit(G4double aVal)
Definition: G4VEmAdjointModel.cc:722
G4int
int G4int
Definition: G4Types.hh:78
G4ContinuousGainOfEnergy.hh
G4AntiNeutron::AntiNeutronDefinition
static G4AntiNeutron * AntiNeutronDefinition()
Definition: G4AntiNeutron.cc:98
G4PionZero::PionZeroDefinition
static G4PionZero * PionZeroDefinition()
Definition: G4PionZero.cc:103
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4AdjointProton::AdjointProtonDefinition
static G4AdjointProton * AdjointProtonDefinition()
Definition: G4AdjointProton.cc:79
G4ionIonisation.hh
G4VEmAdjointModel::SetHighEnergyLimit
void SetHighEnergyLimit(G4double aVal)
Definition: G4VEmAdjointModel.cc:716
G4AdjointCSManager::RegisterAdjointParticle
void RegisterAdjointParticle(G4ParticleDefinition *aPartDef)
Definition: G4AdjointCSManager.cc:161
G4UrbanMscModel.hh
G4InversePEEffect.hh
G4AdjointPhotoElectricModel.hh
G4ParticleTypes.hh
G4AdjointPhysicsList::SetLossFluctuationFlag
void SetLossFluctuationFlag(bool aBool)
Definition: G4AdjointPhysicsList.cc:86
G4ContinuousGainOfEnergy::SetDirectEnergyLossProcess
void SetDirectEnergyLossProcess(G4VEnergyLossProcess *aProcess)
Definition: G4ContinuousGainOfEnergy.hh:112
G4VUserPhysicsList::DumpCutValuesTable
void DumpCutValuesTable(G4int flag=1)
Definition: G4VUserPhysicsList.cc:811
G4cout
G4GLOB_DLL std::ostream G4cout
G4AdjointeIonisationModel.hh
G4AdjointIonIonisationModel.hh
G4GammaConversion.hh
G4PionMinus::PionMinusDefinition
static G4PionMinus * PionMinusDefinition()
Definition: G4PionMinus.cc:93
idxPostStep
Definition: G4ProcessManager.hh:94
G4VEnergyLossProcess::SetLossFluctuations
void SetLossFluctuations(G4bool val)
Definition: G4VEnergyLossProcess.hh:942
G4KaonZeroLong::KaonZeroLongDefinition
static G4KaonZeroLong * KaonZeroLongDefinition()
Definition: G4KaonZeroLong.cc:110
G4VUserPhysicsList::SetVerboseLevel
void SetVerboseLevel(G4int value)
Definition: G4VUserPhysicsList.cc:978
G4ProcessManager::AddProcess
G4int AddProcess(G4VProcess *aProcess, G4int ordAtRestDoIt=ordInActive, G4int ordAlongSteptDoIt=ordInActive, G4int ordPostStepDoIt=ordInActive)
Definition: G4ProcessManager.cc:410
G4NeutrinoE::NeutrinoEDefinition
static G4NeutrinoE * NeutrinoEDefinition()
Definition: G4NeutrinoE.cc:80
G4Decay.hh
G4ContinuousGainOfEnergy::SetDirectParticle
void SetDirectParticle(G4ParticleDefinition *p)
Definition: G4ContinuousGainOfEnergy.cc:99
G4AntiNeutrinoMu::AntiNeutrinoMuDefinition
static G4AntiNeutrinoMu * AntiNeutrinoMuDefinition()
Definition: G4AntiNeutrinoMu.cc:80
G4PionPlus::PionPlusDefinition
static G4PionPlus * PionPlusDefinition()
Definition: G4PionPlus.cc:93
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93
G4ProcessManager.hh
G4AdjointProton.hh
G4VEmAdjointModel::SetUseMatrix
void SetUseMatrix(G4bool aBool)
Definition: G4VEmAdjointModel.hh:205
G4ProcessManager::SetProcessOrdering
void SetProcessOrdering(G4VProcess *aProcess, G4ProcessVectorDoItIndex idDoIt, G4int ordDoIt=ordDefault)
Definition: G4ProcessManager.cc:607
G4PhotoElectricEffect.hh
G4AdjointComptonModel::SetDirectProcess
void SetDirectProcess(G4VEmProcess *aProcess)
Definition: G4AdjointComptonModel.hh:95
G4PhysicalConstants.hh
G4AdjointGamma.hh
G4Positron::PositronDefinition
static G4Positron * PositronDefinition()
Definition: G4Positron.cc:89
G4EtaPrime::EtaPrimeDefinition
static G4EtaPrime * EtaPrimeDefinition()
Definition: G4EtaPrime.cc:106
G4AdjointCSManager.hh
G4eplusAnnihilation.hh
G4hInverseIonisation.hh
G4AdjointProton::AdjointProton
static G4AdjointProton * AdjointProton()
Definition: G4AdjointProton.cc:84
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
Definition: G4ParticleDefinition.cc:247
G4AdjointPhysicsList.hh
Definition of the G4AdjointPhysicsList class.
G4eInverseCompton.hh
G4AdjointSimManager.hh
G4AdjointElectron::AdjointElectronDefinition
static G4AdjointElectron * AdjointElectronDefinition()
Definition: G4AdjointElectron.cc:95
G4AdjointForcedInteractionForGamma.hh
G4AdjointSimManager::ConsiderParticleAsPrimary
void ConsiderParticleAsPrimary(const G4String &particle_name)
Definition: G4AdjointSimManager.cc:581
G4AdjointForcedInteractionForGamma::RegisterAdjointBremModel
void RegisterAdjointBremModel(G4VEmAdjointModel *aAdjointBremModel)
Definition: G4AdjointForcedInteractionForGamma.hh:92
G4VEmAdjointModel::SetCSBiasingFactor
virtual void SetCSBiasingFactor(G4double aVal)
Definition: G4VEmAdjointModel.hh:217
G4MuonMinus::MuonMinusDefinition
static G4MuonMinus * MuonMinusDefinition()
Definition: G4MuonMinus.cc:95
G4VMultipleScattering::AddEmModel
void AddEmModel(G4int order, G4VEmModel *, const G4Region *region=nullptr)
Definition: G4VMultipleScattering.cc:144
G4AdjointPhysicsList::fEminusIonisation
G4eIonisation * fEminusIonisation
Definition: G4AdjointPhysicsList.hh:93
G4AdjointElectron.hh
particleIterator
std::vector< G4InuclElementaryParticle >::iterator particleIterator
Definition: G4BigBanger.cc:65
G4ChargedGeantino::ChargedGeantinoDefinition
static G4ChargedGeantino * ChargedGeantinoDefinition()
Definition: G4ChargedGeantino.cc:82
G4Electron::Electron
static G4Electron * Electron()
Definition: G4Electron.cc:94
G4VUserPhysicsList::GetParticleIterator
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
Definition: G4VUserPhysicsList.cc:972
G4eBremsstrahlung.hh
G4endl
#define G4endl
Definition: G4ios.hh:61
MeV
static constexpr double MeV
Definition: G4SIunits.hh:214
G4AntiNeutrinoE::AntiNeutrinoEDefinition
static G4AntiNeutrinoE * AntiNeutrinoEDefinition()
Definition: G4AntiNeutrinoE.cc:80
G4OpticalPhoton::OpticalPhotonDefinition
static G4OpticalPhoton * OpticalPhotonDefinition()
Definition: G4OpticalPhoton.cc:82
G4SystemOfUnits.hh
G4AdjointBremsstrahlungModel.hh
G4AdjointProcessEquivalentToDirectProcess.hh
G4eInverseBremsstrahlung.hh
G4NeutrinoMu::NeutrinoMuDefinition
static G4NeutrinoMu * NeutrinoMuDefinition()
Definition: G4NeutrinoMu.cc:80
G4AdjointAlongStepWeightCorrection.hh
keV
static constexpr double keV
Definition: G4SIunits.hh:216
G4Neutron::NeutronDefinition
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
G4AdjointCSManager::GetAdjointCSManager
static G4AdjointCSManager * GetAdjointCSManager()
Definition: G4AdjointCSManager.cc:61
G4eInverseIonisation.hh
idxAtRest
Definition: G4ProcessManager.hh:92
G4AdjointPhysicsMessenger.hh
Definition of the G4AdjointPhysicsMessenger class.
G4AdjointPhysicsList::fPIonisation
G4hIonisation * fPIonisation
Definition: G4AdjointPhysicsList.hh:94
G4Eta::EtaDefinition
static G4Eta * EtaDefinition()
Definition: G4Eta.cc:104
G4Gamma::GammaDefinition
static G4Gamma * GammaDefinition()
Definition: G4Gamma.cc:81
idxAlongStep
Definition: G4ProcessManager.hh:93