Geant4  10.02.p03
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 102356 2017-01-23 16:22:42Z 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_eionisation(true),fUse_pionisation(true),
54  fUse_brem(true),fUse_compton(true),fUse_ms(true),
55  fUse_egain_fluctuation(true),fUse_peeffect(true),
56  fEmin_adj_models(1.*keV), fEmax_adj_models(1.*MeV),
57  fCS_biasing_factor_compton(1.),fCS_biasing_factor_brem(1.),
58  fCS_biasing_factor_ionisation(1.),fCS_biasing_factor_PEeffect(1.)
59 {
60  defaultCutValue = 1.0*mm;
61  SetVerboseLevel(1);
62  fPhysicsMessenger = new G4AdjointPhysicsMessenger(this);
63 }
64 
65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66 
67 G4AdjointPhysicsList::~G4AdjointPhysicsList()
68 {
69 }
70 void G4AdjointPhysicsList::ConstructParticle()
71 {
72  // In this method, static member functions should be called
73  // for all particles which you want to use.
74  // This ensures that objects of these particle types will be
75  // created in the program.
76  ConstructBosons();
77  ConstructLeptons();
78  ConstructMesons();
79  ConstructBaryons();
80  ConstructAdjointParticles();
81 }
82 
83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84 
85 void G4AdjointPhysicsList::SetLossFluctuationFlag(bool aBool)
86 {
87  if (fEminusIonisation) fEminusIonisation->SetLossFluctuations(aBool);
88 }
89 
90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
91 
92 void G4AdjointPhysicsList::ConstructBosons()
93 {
94  // pseudo-particles
95  G4Geantino::GeantinoDefinition();
96  G4ChargedGeantino::ChargedGeantinoDefinition();
97 
98  // gamma
99  G4Gamma::GammaDefinition();
100 
101  // optical photon
102  G4OpticalPhoton::OpticalPhotonDefinition();
103 }
104 
105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
106 
107 void G4AdjointPhysicsList::ConstructLeptons()
108 {
109  // leptons
110  G4Electron::ElectronDefinition();
111  G4Positron::PositronDefinition();
112  G4MuonPlus::MuonPlusDefinition();
113  G4MuonMinus::MuonMinusDefinition();
114 
115  G4NeutrinoE::NeutrinoEDefinition();
116  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
117  G4NeutrinoMu::NeutrinoMuDefinition();
118  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
119 }
120 
121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
122 
123 void G4AdjointPhysicsList::ConstructMesons()
124 {
125 // mesons
126  G4PionPlus::PionPlusDefinition();
127  G4PionMinus::PionMinusDefinition();
128  G4PionZero::PionZeroDefinition();
129  G4Eta::EtaDefinition();
130  G4EtaPrime::EtaPrimeDefinition();
131  G4KaonPlus::KaonPlusDefinition();
132  G4KaonMinus::KaonMinusDefinition();
133  G4KaonZero::KaonZeroDefinition();
134  G4AntiKaonZero::AntiKaonZeroDefinition();
135  G4KaonZeroLong::KaonZeroLongDefinition();
136  G4KaonZeroShort::KaonZeroShortDefinition();
137 }
138 
139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
140 
141 void G4AdjointPhysicsList::ConstructBaryons()
142 {
143 // barions
144  G4Proton::ProtonDefinition();
145  G4AntiProton::AntiProtonDefinition();
146  G4Neutron::NeutronDefinition();
147  G4AntiNeutron::AntiNeutronDefinition();
148 }
149 
150 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
151 
152 #include"G4AdjointGamma.hh"
153 #include"G4AdjointElectron.hh"
154 #include"G4AdjointProton.hh"
155 void G4AdjointPhysicsList::ConstructAdjointParticles()
156 {
157 // adjoint_gammma
158  G4AdjointGamma::AdjointGammaDefinition();
159 
160 // adjoint_electron
161  G4AdjointElectron::AdjointElectronDefinition();
162 
163 // adjoint_proton
164  G4AdjointProton::AdjointProtonDefinition();
165 }
166 
167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
168 
169 void G4AdjointPhysicsList::ConstructProcess()
170 {
171  AddTransportation();
172  ConstructEM();
173  ConstructGeneral();
174 }
175 
176 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
177 
178 //#include "G4PEEffectFluoModel.hh"
179 #include "G4ComptonScattering.hh"
180 #include "G4GammaConversion.hh"
181 #include "G4PhotoElectricEffect.hh"
182 #include "G4eMultipleScattering.hh"
183 #include "G4hMultipleScattering.hh"
184 #include "G4eIonisation.hh"
185 #include "G4eBremsstrahlung.hh"
186 #include "G4eplusAnnihilation.hh"
187 #include "G4hIonisation.hh"
188 #include "G4ionIonisation.hh"
189 //#include "G4IonParametrisedLossModel.hh"
190 
191 #include "G4ContinuousGainOfEnergy.hh"
192 #include "G4eInverseIonisation.hh"
193 #include "G4AdjointeIonisationModel.hh"
194 #include "G4AdjointCSManager.hh"
195 #include "G4AdjointBremsstrahlungModel.hh"
196 #include "G4eInverseBremsstrahlung.hh"
197 #include "G4AdjointComptonModel.hh"
198 #include "G4eInverseCompton.hh"
199 #include "G4InversePEEffect.hh"
200 #include "G4AdjointPhotoElectricModel.hh"
201 #include "G4AdjointAlongStepWeightCorrection.hh"
202 #include "G4hInverseIonisation.hh"
203 #include "G4AdjointhIonisationModel.hh"
204 #include "G4AdjointhMultipleScattering.hh"
205 #include "G4IonInverseIonisation.hh"
206 #include "G4AdjointIonIonisationModel.hh"
207 
208 #include "G4AdjointSimManager.hh"
209 #include "G4AdjointProcessEquivalentToDirectProcess.hh"
210 
211 void G4AdjointPhysicsList::ConstructEM()
212 { G4AdjointCSManager* theCSManager =
213  G4AdjointCSManager::GetAdjointCSManager();
214 
215  G4AdjointSimManager* theAdjointSimManager =
216  G4AdjointSimManager::GetInstance();
217 
218  theCSManager->RegisterAdjointParticle(
219  G4AdjointElectron::AdjointElectron());
220 
221  if (fUse_brem || fUse_peeffect ||fUse_compton)
222  theCSManager->RegisterAdjointParticle(
223  G4AdjointGamma::AdjointGamma());
224 
225  if (fUse_eionisation) {
226  if (!fEminusIonisation) fEminusIonisation =
227  new G4eIonisation();
228  fEminusIonisation->SetLossFluctuations(fUse_egain_fluctuation);
229  }
230 
231  if (fUse_pionisation) {
232  if (!fPIonisation) fPIonisation = new G4hIonisation();
233  fPIonisation->SetLossFluctuations(fUse_egain_fluctuation);
234  theCSManager->RegisterAdjointParticle(
235  G4AdjointProton::AdjointProton());
236  }
237 
238  G4eBremsstrahlung* theeminusBremsstrahlung = 0;
239  if (fUse_brem && fUse_eionisation)
240  theeminusBremsstrahlung = new G4eBremsstrahlung();
241 
242  G4ComptonScattering* theComptonScattering =0;
243  if (fUse_compton) theComptonScattering = new G4ComptonScattering();
244 
245  G4PhotoElectricEffect* thePEEffect =0;
246  if (fUse_peeffect) thePEEffect = new G4PhotoElectricEffect();
247 
248  G4eMultipleScattering* theeminusMS = 0;
249  G4hMultipleScattering* thepMS= 0;
250  if (fUse_ms) {
251  theeminusMS = new G4eMultipleScattering();
252  thepMS = new G4hMultipleScattering();
253  }
254 
255  G4VProcess* theGammaConversion =0;
256  if (fUse_gamma_conversion) theGammaConversion = new G4GammaConversion();
257 
258  //Define adjoint e- ionisation
259  //-------------------
260  G4AdjointeIonisationModel* theeInverseIonisationModel = 0;
261  G4eInverseIonisation* theeInverseIonisationProjToProjCase = 0 ;
262  G4eInverseIonisation* theeInverseIonisationProdToProjCase = 0;
263  if (fUse_eionisation) {
264  theeInverseIonisationModel = new G4AdjointeIonisationModel();
265  theeInverseIonisationModel->SetHighEnergyLimit(
266  fEmax_adj_models);
267  theeInverseIonisationModel->SetLowEnergyLimit(
268  fEmin_adj_models);
269  theeInverseIonisationModel->SetCSBiasingFactor(
270  fCS_biasing_factor_ionisation);
271  theeInverseIonisationProjToProjCase =
272  new G4eInverseIonisation(true,"Inv_eIon",
273  theeInverseIonisationModel);
274  theeInverseIonisationProdToProjCase =
275  new G4eInverseIonisation(false,"Inv_eIon1",
276  theeInverseIonisationModel);
277  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
278  }
279 
280  //Define adjoint Bremsstrahlung
281  //-------------------------------
282  G4AdjointBremsstrahlungModel* theeInverseBremsstrahlungModel = 0;
283  G4eInverseBremsstrahlung* theeInverseBremsstrahlungProjToProjCase = 0;
284  G4eInverseBremsstrahlung* theeInverseBremsstrahlungProdToProjCase = 0;
285 
286  if (fUse_brem && fUse_eionisation) {
287  theeInverseBremsstrahlungModel = new G4AdjointBremsstrahlungModel();
288  theeInverseBremsstrahlungModel->SetHighEnergyLimit(fEmax_adj_models);
289  theeInverseBremsstrahlungModel->SetLowEnergyLimit(fEmin_adj_models);
290  theeInverseBremsstrahlungModel->SetCSBiasingFactor(
291  fCS_biasing_factor_brem);
292  theeInverseBremsstrahlungProjToProjCase = new G4eInverseBremsstrahlung(
293  true,"Inv_eBrem",theeInverseBremsstrahlungModel);
294  theeInverseBremsstrahlungProdToProjCase = new G4eInverseBremsstrahlung(
295  false,"Inv_eBrem1",theeInverseBremsstrahlungModel);
296  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
297  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("gamma"));
298  }
299 
300 
301  //Define adjoint Compton
302  //---------------------
303 
304  G4AdjointComptonModel* theeInverseComptonModel = 0;
305  G4eInverseCompton* theeInverseComptonProjToProjCase = 0;
306  G4eInverseCompton* theeInverseComptonProdToProjCase = 0;
307 
308  if (fUse_compton) {
309  theeInverseComptonModel = new G4AdjointComptonModel();
310  theeInverseComptonModel->SetHighEnergyLimit(fEmax_adj_models);
311  theeInverseComptonModel->SetLowEnergyLimit(fEmin_adj_models);
312  theeInverseComptonModel->SetDirectProcess(theComptonScattering);
313  theeInverseComptonModel->SetUseMatrix(false);
314  theeInverseComptonModel->SetCSBiasingFactor( fCS_biasing_factor_compton);
315  theeInverseComptonProjToProjCase = new G4eInverseCompton(true,"Inv_Compt",
316  theeInverseComptonModel);
317  theeInverseComptonProdToProjCase = new G4eInverseCompton(false,"Inv_Compt1",
318  theeInverseComptonModel);
319  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
320  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("gamma"));
321  }
322 
323  //Define adjoint PEEffect
324  //---------------------
325  G4AdjointPhotoElectricModel* theInversePhotoElectricModel = 0;
326  G4InversePEEffect* theInversePhotoElectricProcess = 0;
327 
328  if (fUse_peeffect) {
329  theInversePhotoElectricModel = new G4AdjointPhotoElectricModel();
330  theInversePhotoElectricModel->SetHighEnergyLimit(fEmax_adj_models);
331  theInversePhotoElectricModel->SetLowEnergyLimit(fEmin_adj_models);
332  theInversePhotoElectricModel->SetCSBiasingFactor(
333  fCS_biasing_factor_PEeffect);
334  theInversePhotoElectricProcess = new G4InversePEEffect("Inv_PEEffect",
335  theInversePhotoElectricModel);
336  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
337  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("gamma"));
338  }
339 
340 
341  //Define adjoint ionisation for protons
342  //---------------------
343  G4AdjointhIonisationModel* thepInverseIonisationModel = 0;
344  G4hInverseIonisation* thepInverseIonisationProjToProjCase = 0 ;
345  G4hInverseIonisation* thepInverseIonisationProdToProjCase = 0;
346  if (fUse_pionisation) {
347  thepInverseIonisationModel = new G4AdjointhIonisationModel(
348  G4Proton::Proton());
349  thepInverseIonisationModel->SetHighEnergyLimit(fEmax_adj_models);
350  thepInverseIonisationModel->SetLowEnergyLimit(fEmin_adj_models);
351  thepInverseIonisationModel->SetUseMatrix(false);
352  thepInverseIonisationProjToProjCase = new G4hInverseIonisation(true,
353  "Inv_pIon",thepInverseIonisationModel);
354  thepInverseIonisationProdToProjCase = new G4hInverseIonisation(false,
355  "Inv_pIon1",thepInverseIonisationModel);
356  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("e-"));
357  theAdjointSimManager->ConsiderParticleAsPrimary(G4String("proton"));
358  }
359 
360 
361 
362  //Declare the processes active for the different particles
363  //--------------------------------------------------------
364  auto particleIterator=GetParticleIterator();
365  particleIterator->reset();
366  while( (*particleIterator)() ){
367  G4ParticleDefinition* particle = particleIterator->value();
368  G4ProcessManager* pmanager = particle->GetProcessManager();
369  if (!pmanager) {
370  pmanager = new G4ProcessManager(particle);
371  particle->SetProcessManager(pmanager);
372  }
373 
374  G4String particleName = particle->GetParticleName();
375  if (particleName == "e-") {
376  if (fUse_ms && fUse_eionisation) pmanager->AddProcess(theeminusMS);
377  if (fUse_eionisation){
378  pmanager->AddProcess(fEminusIonisation);
379  G4AdjointCSManager::GetAdjointCSManager()->
380  RegisterEnergyLossProcess(fEminusIonisation,particle);
381  }
382  if (fUse_brem && fUse_eionisation) {
383  pmanager->AddProcess(theeminusBremsstrahlung);
384  G4AdjointCSManager::GetAdjointCSManager()->
385  RegisterEnergyLossProcess(theeminusBremsstrahlung,particle);
386  }
387  G4int n_order=0;
388  if (fUse_ms && fUse_eionisation) {
389  n_order++;
390  pmanager->SetProcessOrdering(theeminusMS, idxAlongStep,n_order);
391  }
392  if (fUse_eionisation) {
393  n_order++;
394  pmanager->SetProcessOrdering(fEminusIonisation,idxAlongStep,n_order);
395  }
396  if (fUse_brem && fUse_eionisation) {
397  n_order++;
398  pmanager->SetProcessOrdering(theeminusBremsstrahlung,
399  idxAlongStep,n_order);
400  }
401  n_order=0;
402  if (fUse_ms && fUse_eionisation) {
403  n_order++;
404  pmanager->SetProcessOrdering(theeminusMS,idxPostStep,n_order);
405  }
406  if (fUse_eionisation) {
407  n_order++;
408  pmanager->SetProcessOrdering(fEminusIonisation,idxPostStep,n_order);
409  }
410  if (fUse_brem && fUse_eionisation) {
411  n_order++;
412  pmanager->SetProcessOrdering(theeminusBremsstrahlung,idxPostStep,
413  n_order);
414  }
415  }
416 
417  if (particleName == "adj_e-") {
418  G4ContinuousGainOfEnergy* theContinuousGainOfEnergy =0;
419  if (fUse_eionisation ) {
420  theContinuousGainOfEnergy= new G4ContinuousGainOfEnergy();
421  theContinuousGainOfEnergy->SetLossFluctuations(
422  fUse_egain_fluctuation);
423  theContinuousGainOfEnergy->SetDirectEnergyLossProcess(
424  fEminusIonisation);
425  theContinuousGainOfEnergy->SetDirectParticle(G4Electron::Electron());
426  pmanager->AddProcess(theContinuousGainOfEnergy);
427  }
428  G4int n_order=0;
429  if (fUse_ms) {
430  n_order++;
431  pmanager->AddProcess(theeminusMS);
432  pmanager->SetProcessOrdering(theeminusMS, idxAlongStep,n_order);
433  }
434 
435  n_order++;
436  pmanager->SetProcessOrdering(theContinuousGainOfEnergy,idxAlongStep,
437  n_order);
438 
439 
440  n_order++;
441  G4AdjointAlongStepWeightCorrection* theAlongStepWeightCorrection =
442  new G4AdjointAlongStepWeightCorrection();
443  pmanager->AddProcess(theAlongStepWeightCorrection);
444  pmanager->SetProcessOrdering(theAlongStepWeightCorrection,
445  idxAlongStep,
446  n_order);
447  n_order=0;
448  if (fUse_eionisation) {
449  pmanager->AddProcess(theeInverseIonisationProjToProjCase);
450  pmanager->AddProcess(theeInverseIonisationProdToProjCase);
451  n_order++;
452  pmanager->SetProcessOrdering(theeInverseIonisationProjToProjCase,
453  idxPostStep,n_order);
454  n_order++;
455  pmanager->SetProcessOrdering(theeInverseIonisationProdToProjCase,
456  idxPostStep,n_order);
457  }
458 
459  if (fUse_brem && fUse_eionisation) {
460  pmanager->AddProcess(theeInverseBremsstrahlungProjToProjCase);
461  n_order++;
462  pmanager->SetProcessOrdering(
463  theeInverseBremsstrahlungProjToProjCase,
464  idxPostStep,n_order);
465  }
466 
467  if (fUse_compton) {
468  pmanager->AddProcess(theeInverseComptonProdToProjCase);
469  n_order++;
470  pmanager->SetProcessOrdering(theeInverseComptonProdToProjCase,
471  idxPostStep,n_order);
472  }
473  if (fUse_peeffect) {
474  pmanager->AddDiscreteProcess(theInversePhotoElectricProcess);
475  n_order++;
476  pmanager->SetProcessOrdering(theInversePhotoElectricProcess,
477  idxPostStep,n_order);
478  }
479  if (fUse_pionisation) {
480  pmanager->AddProcess(thepInverseIonisationProdToProjCase);
481  n_order++;
482  pmanager->SetProcessOrdering(thepInverseIonisationProdToProjCase,
483  idxPostStep,n_order);
484  }
485  if (fUse_ms && fUse_eionisation) {
486  n_order++;
487  pmanager->SetProcessOrdering(theeminusMS,idxPostStep,n_order);
488  }
489  }
490 
491 
492  if(particleName == "adj_gamma") {
493  G4int n_order=0;
494  G4AdjointAlongStepWeightCorrection* theAlongStepWeightCorrection =
495  new G4AdjointAlongStepWeightCorrection();
496  pmanager->AddProcess(theAlongStepWeightCorrection);
497  pmanager->SetProcessOrdering(theAlongStepWeightCorrection,
498  idxAlongStep,1);
499 
500  if (fUse_brem && fUse_eionisation) {
501  pmanager->AddProcess(theeInverseBremsstrahlungProdToProjCase);
502  n_order++;
503  pmanager->SetProcessOrdering(
504  theeInverseBremsstrahlungProdToProjCase,
505  idxPostStep,n_order);
506  }
507  if (fUse_compton) {
508  pmanager->AddDiscreteProcess(theeInverseComptonProjToProjCase);
509  n_order++;
510  pmanager->SetProcessOrdering(theeInverseComptonProjToProjCase,
511  idxPostStep,n_order);
512  }
513  }
514 
515  if (particleName == "gamma") {
516  if (fUse_compton) {
517  pmanager->AddDiscreteProcess(theComptonScattering);
518  G4AdjointCSManager::GetAdjointCSManager()->
519  RegisterEmProcess(theComptonScattering,particle);
520  }
521  if (fUse_peeffect) {
522  pmanager->AddDiscreteProcess(thePEEffect);
523  G4AdjointCSManager::GetAdjointCSManager()->
524  RegisterEmProcess(thePEEffect,particle);
525  }
526  if (fUse_gamma_conversion) {
527  pmanager->AddDiscreteProcess(theGammaConversion);
528  }
529  }
530 
531  if (particleName == "e+" && fUse_gamma_conversion) {//positron
532  G4VProcess* theeplusMultipleScattering = new G4eMultipleScattering();
533  G4VProcess* theeplusIonisation = new G4eIonisation();
534  G4VProcess* theeplusBremsstrahlung = new G4eBremsstrahlung();
535  G4VProcess* theeplusAnnihilation = new G4eplusAnnihilation();
536 
537  // add processes
538  pmanager->AddProcess(theeplusMultipleScattering);
539  pmanager->AddProcess(theeplusIonisation);
540  pmanager->AddProcess(theeplusBremsstrahlung);
541  pmanager->AddProcess(theeplusAnnihilation);
542 
543  // set ordering for AtRestDoIt
544  pmanager->SetProcessOrderingToFirst(theeplusAnnihilation, idxAtRest);
545 
546  // set ordering for AlongStepDoIt
547  pmanager->SetProcessOrdering(theeplusMultipleScattering,
548  idxAlongStep,1);
549  pmanager->SetProcessOrdering(theeplusIonisation, idxAlongStep,2);
550  pmanager->SetProcessOrdering(theeplusBremsstrahlung,idxAlongStep,3);
551 
552  // set ordering for PostStepDoIt
553  pmanager->SetProcessOrdering(theeplusMultipleScattering,
554  idxPostStep,1);
555  pmanager->SetProcessOrdering(theeplusIonisation,idxPostStep,2);
556  pmanager->SetProcessOrdering(theeplusBremsstrahlung,idxPostStep,3);
557  pmanager->SetProcessOrdering(theeplusAnnihilation,idxPostStep,4);
558  }
559  if (particleName == "proton" && fUse_pionisation) {
560  if (fUse_ms && fUse_pionisation) pmanager->AddProcess(thepMS);
561 
562  if (fUse_pionisation){
563  pmanager->AddProcess(fPIonisation);
564  G4AdjointCSManager::GetAdjointCSManager()->
565  RegisterEnergyLossProcess(fPIonisation,particle);
566  }
567 
568  G4int n_order=0;
569  if (fUse_ms && fUse_pionisation) {
570  n_order++;
571  pmanager->SetProcessOrdering(thepMS, idxAlongStep,n_order);
572  }
573 
574  if (fUse_pionisation) {
575  n_order++;
576  pmanager->SetProcessOrdering(fPIonisation,idxAlongStep,n_order);
577  }
578 
579  n_order=0;
580  if (fUse_ms && fUse_pionisation) {
581  n_order++;
582  pmanager->SetProcessOrdering(thepMS, idxPostStep,n_order);
583  }
584 
585  if (fUse_pionisation) {
586  n_order++;
587  pmanager->SetProcessOrdering(fPIonisation,idxPostStep,n_order);
588  }
589 
590  }
591 
592  if (particleName == "adj_proton" && fUse_pionisation) {
593  G4ContinuousGainOfEnergy* theContinuousGainOfEnergy =0;
594  if (fUse_pionisation ) {
595  theContinuousGainOfEnergy= new G4ContinuousGainOfEnergy();
596  theContinuousGainOfEnergy->SetLossFluctuations(
597  fUse_egain_fluctuation);
598  theContinuousGainOfEnergy->SetDirectEnergyLossProcess(fPIonisation);
599  theContinuousGainOfEnergy->SetDirectParticle(G4Proton::Proton());
600  pmanager->AddProcess(theContinuousGainOfEnergy);
601  }
602 
603  G4int n_order=0;
604  if (fUse_ms) {
605  n_order++;
606  pmanager->AddProcess(thepMS);
607  pmanager->SetProcessOrdering(thepMS, idxAlongStep,n_order);
608  }
609 
610  n_order++;
611  pmanager->SetProcessOrdering(theContinuousGainOfEnergy,idxAlongStep,
612  n_order);
613 
614  n_order++;
615  G4AdjointAlongStepWeightCorrection* theAlongStepWeightCorrection =
616  new G4AdjointAlongStepWeightCorrection();
617  pmanager->AddProcess(theAlongStepWeightCorrection);
618  pmanager->SetProcessOrdering(theAlongStepWeightCorrection,
619  idxAlongStep,
620  n_order);
621  n_order=0;
622  if (fUse_pionisation) {
623  pmanager->AddProcess(thepInverseIonisationProjToProjCase);
624  n_order++;
625  pmanager->SetProcessOrdering(
626  thepInverseIonisationProjToProjCase,
627  idxPostStep,n_order);
628  }
629 
630  if (fUse_ms && fUse_pionisation) {
631  n_order++;
632  pmanager->SetProcessOrdering(thepMS,idxPostStep,n_order);
633  }
634  }
635  }
636 }
637 
638 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
639 
640 #include "G4Decay.hh"
641 void G4AdjointPhysicsList::ConstructGeneral()
642 {
643  // Add Decay Process
644  G4Decay* theDecayProcess = new G4Decay();
645  auto particleIterator=GetParticleIterator();
646  particleIterator->reset();
647  while( (*particleIterator)() ){
648  G4ParticleDefinition* particle = particleIterator->value();
649  G4ProcessManager* pmanager = particle->GetProcessManager();
650  if (theDecayProcess->IsApplicable(*particle)) {
651  pmanager ->AddProcess(theDecayProcess);
652  // set ordering for PostStepDoIt and AtRestDoIt
653  pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
654  pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
655  }
656  }
657 }
658 
659 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
660 
661 void G4AdjointPhysicsList::SetCuts()
662 {
663  if (verboseLevel >0){
664  G4cout << "G4AdjointPhysicsList::SetCuts:";
665  G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
666  }
667 
668  // set cut values for gamma at first and for e- second and next for e+,
669  // because some processes for e+/e- need cut values for gamma
670  //
671  SetCutValue(defaultCutValue, "gamma");
672  SetCutValue(defaultCutValue, "e-");
673  SetCutValue(defaultCutValue, "e+");
674 
675  if (verboseLevel>0) DumpCutValuesTable();
676 }
677 
678 //....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
MeV
static const double MeV
Definition: G4SIunits.hh:211
G4AdjointGamma::AdjointGammaDefinition
static G4AdjointGamma * AdjointGammaDefinition()
Definition: G4AdjointGamma.cc:79
G4MuonPlus::MuonPlusDefinition
static G4MuonPlus * MuonPlusDefinition()
Definition: G4MuonPlus.cc:94
G4IonInverseIonisation.hh
G4AdjointPhysicsList::fPhysicsMessenger
G4AdjointPhysicsMessenger * fPhysicsMessenger
Definition: G4AdjointPhysicsList.hh:97
G4ParticleDefinition::SetProcessManager
void SetProcessManager(G4ProcessManager *aProcessManager)
Definition: G4ParticleDefinition.cc:450
G4VUserPhysicsList::SetCutValue
void SetCutValue(G4double aCut, const G4String &pname)
Definition: G4VUserPhysicsList.cc:415
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:70
G4hIonisation.hh
G4KaonMinus::KaonMinusDefinition
static G4KaonMinus * KaonMinusDefinition()
Definition: G4KaonMinus.cc:108
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
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
Definition: G4ParticleDefinition.cc:247
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:721
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
G4AdjointProton::AdjointProtonDefinition
static G4AdjointProton * AdjointProtonDefinition()
Definition: G4AdjointProton.cc:79
G4ionIonisation.hh
G4VEmAdjointModel::SetHighEnergyLimit
void SetHighEnergyLimit(G4double aVal)
Definition: G4VEmAdjointModel.cc:715
G4AdjointCSManager::RegisterAdjointParticle
void RegisterAdjointParticle(G4ParticleDefinition *aPartDef)
Definition: G4AdjointCSManager.cc:161
G4InversePEEffect.hh
G4AdjointPhotoElectricModel.hh
G4ParticleTypes.hh
G4AdjointPhysicsList::SetLossFluctuationFlag
void SetLossFluctuationFlag(bool aBool)
Definition: G4AdjointPhysicsList.cc:85
G4ContinuousGainOfEnergy::SetDirectEnergyLossProcess
void SetDirectEnergyLossProcess(G4VEnergyLossProcess *aProcess)
Definition: G4ContinuousGainOfEnergy.hh:112
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4VUserPhysicsList::DumpCutValuesTable
void DumpCutValuesTable(G4int flag=1)
Definition: G4VUserPhysicsList.cc:810
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:940
G4KaonZeroLong::KaonZeroLongDefinition
static G4KaonZeroLong * KaonZeroLongDefinition()
Definition: G4KaonZeroLong.cc:110
G4VUserPhysicsList::SetVerboseLevel
void SetVerboseLevel(G4int value)
Definition: G4VUserPhysicsList.cc:977
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
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
G4AdjointPhysicsList.hh
Definition of the G4AdjointPhysicsList class.
G4VUserPhysicsList::GetParticleIterator
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
Definition: G4VUserPhysicsList.cc:971
G4eInverseCompton.hh
G4AdjointSimManager.hh
G4AdjointElectron::AdjointElectronDefinition
static G4AdjointElectron * AdjointElectronDefinition()
Definition: G4AdjointElectron.cc:95
G4AdjointSimManager::ConsiderParticleAsPrimary
void ConsiderParticleAsPrimary(const G4String &particle_name)
Definition: G4AdjointSimManager.cc:550
G4VEmAdjointModel::SetCSBiasingFactor
virtual void SetCSBiasingFactor(G4double aVal)
Definition: G4VEmAdjointModel.hh:217
G4MuonMinus::MuonMinusDefinition
static G4MuonMinus * MuonMinusDefinition()
Definition: G4MuonMinus.cc:95
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
G4eBremsstrahlung.hh
G4endl
#define G4endl
Definition: G4ios.hh:61
keV
static const double keV
Definition: G4SIunits.hh:213
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
mm
static const double mm
Definition: G4SIunits.hh:114
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