| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 << 26 // $Id: G4VAdjointReverseReaction.cc,v 1.3 2010-11-11 11:51:56 ldesorgh Exp $
>> 27 // GEANT4 tag $Name: not supported by cvs2svn $
>> 28 //
27 #include "G4VAdjointReverseReaction.hh" 29 #include "G4VAdjointReverseReaction.hh"
28 <<
29 #include "G4AdjointCSManager.hh" 30 #include "G4AdjointCSManager.hh"
30 #include "G4ParticleChange.hh" << 31 #include "G4AdjointCSMatrix.hh"
>> 32 #include "G4AdjointInterpolator.hh"
>> 33 #include "G4AdjointCSMatrix.hh"
31 #include "G4VEmAdjointModel.hh" 34 #include "G4VEmAdjointModel.hh"
>> 35 #include "G4ElementTable.hh"
>> 36 #include "G4Element.hh"
>> 37 #include "G4Material.hh"
>> 38 #include "G4MaterialCutsCouple.hh"
>> 39 #include "G4AdjointCSManager.hh"
>> 40 #include "G4ParticleChange.hh"
>> 41 #include "G4AdjointElectron.hh"
32 42
33 G4VAdjointReverseReaction::G4VAdjointReverseRe <<
34 <<
35 : G4VDiscreteProcess(process_name) <<
36 { <<
37 fCSManager = G4AdjointCSManager::GetA <<
38 fIsScatProjToProj = whichScatCase; <<
39 fParticleChange = new G4ParticleChange(); <<
40 } <<
41 43
>> 44 G4VAdjointReverseReaction::
>> 45 G4VAdjointReverseReaction(G4String process_name, G4bool whichScatCase):
>> 46 G4VDiscreteProcess(process_name)
>> 47 {theAdjointCSManager = G4AdjointCSManager::GetAdjointCSManager();
>> 48 IsScatProjToProjCase=whichScatCase;
>> 49 fParticleChange=new G4ParticleChange();
>> 50 IsFwdCSUsed=false;
>> 51 IsIntegralModeUsed=false;
>> 52 currentMaterialIndex=100000;
>> 53 currentTcut=0.;
>> 54 lastCS=0.;
>> 55 }
42 ////////////////////////////////////////////// 56 //////////////////////////////////////////////////////////////////////////////
43 G4VAdjointReverseReaction::~G4VAdjointReverseR << 57 //
44 { << 58 G4VAdjointReverseReaction::
45 if(fParticleChange) << 59 ~G4VAdjointReverseReaction()
46 delete fParticleChange; << 60 { if (fParticleChange) delete fParticleChange;
>> 61 }
>> 62 //////////////////////////////////////////////////////////////////////////////
>> 63 //
>> 64 void G4VAdjointReverseReaction::PreparePhysicsTable(const G4ParticleDefinition&)
>> 65 {;
47 } 66 }
48 <<
49 ////////////////////////////////////////////// 67 //////////////////////////////////////////////////////////////////////////////
>> 68 //
50 void G4VAdjointReverseReaction::BuildPhysicsTa 69 void G4VAdjointReverseReaction::BuildPhysicsTable(const G4ParticleDefinition&)
51 { 70 {
52 fCSManager->BuildCrossSectionMatrices(); // <<
53 fCSManager->BuildTotalSigmaTables(); <<
54 } <<
55 71
>> 72 theAdjointCSManager->BuildCrossSectionMatrices(); //do not worry it will be done just once
>> 73 theAdjointCSManager->BuildTotalSigmaTables();
>> 74
>> 75 }
56 ////////////////////////////////////////////// 76 //////////////////////////////////////////////////////////////////////////////
57 G4VParticleChange* G4VAdjointReverseReaction:: << 77 //
58 << 78 G4VParticleChange* G4VAdjointReverseReaction::PostStepDoIt(const G4Track& track, const G4Step& )
59 { << 79 {
>> 80
>> 81
>> 82
60 fParticleChange->Initialize(track); 83 fParticleChange->Initialize(track);
61 fAdjointModel->SampleSecondaries(track, fIsS << 84
62 << 85 /* if (IsFwdCSUsed && IsIntegralModeUsed){ //INtegral mode still unstable
>> 86 G4double Tkin = step.GetPostStepPoint()->GetKineticEnergy();
>> 87 G4double fwdCS = theAdjointCSManager->GetTotalForwardCS(track.GetDefinition(), Tkin, track.GetMaterialCutsCouple());
>> 88 //G4cout<<"lastCS "<<lastCS<<G4endl;
>> 89 if (fwdCS<lastCS*G4UniformRand()) { // the reaction does not take place, same integral method as the one used for forward ionisation in G4
>> 90 ClearNumberOfInteractionLengthLeft();
>> 91 return fParticleChange;
>> 92 }
>> 93
>> 94 }
>> 95 */
>> 96
>> 97 theAdjointEMModel->SampleSecondaries(track,
>> 98 IsScatProjToProjCase,
>> 99 fParticleChange);
>> 100
63 ClearNumberOfInteractionLengthLeft(); 101 ClearNumberOfInteractionLengthLeft();
64 return fParticleChange; 102 return fParticleChange;
>> 103
>> 104
>> 105
65 } 106 }
66 <<
67 ////////////////////////////////////////////// 107 //////////////////////////////////////////////////////////////////////////////
>> 108 //
68 G4double G4VAdjointReverseReaction::GetMeanFre 109 G4double G4VAdjointReverseReaction::GetMeanFreePath(const G4Track& track,
69 << 110 G4double ,
70 << 111 G4ForceCondition* condition)
71 { << 112 { *condition = NotForced;
72 *condition = NotForced; <<
73 G4double preStepKinEnergy = track.GetKinetic 113 G4double preStepKinEnergy = track.GetKineticEnergy();
74 << 114
75 if(track.GetTrackID() != fTrackId) << 115 /*G4double Sigma =
76 { << 116 theAdjointEMModel->AdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase);*/
77 fTrackId = track.GetTrackID(); << 117
78 } << 118 G4double Sigma =
79 G4double sigma = fAdjointModel->AdjointCross << 119 theAdjointEMModel->GetAdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase);
80 track.GetMaterialCutsCouple(), preStepKinE << 120 G4double fwd_TotCS;
81 << 121 Sigma *= theAdjointCSManager->GetCrossSectionCorrection(track.GetDefinition(),preStepKinEnergy,track.GetMaterialCutsCouple(),IsFwdCSUsed, fwd_TotCS);
82 G4double corr = fCSManager->GetCrossSectionC << 122 //G4cout<<fwd_TotCS<<G4endl;
83 track.GetDefinition(), preStepKinEnergy, t << 123 /*if (IsFwdCSUsed && IsIntegralModeUsed){ //take the maximum cross section only for charged particle
84 fIsFwdCSUsed); << 124 G4double e_sigma_max, sigma_max;
85 << 125 theAdjointCSManager->GetMaxFwdTotalCS(track.GetDefinition(),
86 if(std::fabs(corr) > 100.) << 126 track.GetMaterialCutsCouple(), e_sigma_max, sigma_max);
87 { << 127 if (e_sigma_max > preStepKinEnergy){
88 sigma = 0.0; << 128 Sigma*=sigma_max/fwd_TotCS;
89 } << 129 }
90 else <<
91 { <<
92 sigma *= corr; <<
93 } 130 }
>> 131 */
94 132
95 G4double mean_free_path = 1.e60; << 133 G4double mean_free_path = 1.e60 *mm;
96 if(sigma > 0.) << 134 if (Sigma>0) mean_free_path = 1./Sigma;
97 mean_free_path = 1. / sigma; << 135 lastCS=Sigma;
>> 136
>> 137 /*G4cout<<"Sigma "<<Sigma<<G4endl;
>> 138 G4cout<<"mean_free_path [mm] "<<mean_free_path/mm<<G4endl;
>> 139 */
>> 140
98 141
99 return mean_free_path; 142 return mean_free_path;
100 } << 143 }
101 144