| 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 #include "G4ChannelingOptrMultiParticleChangeC 26 #include "G4ChannelingOptrMultiParticleChangeCrossSection.hh"
27 #include "G4ChannelingOptrChangeCrossSection.h 27 #include "G4ChannelingOptrChangeCrossSection.hh"
28 #include "G4BiasingProcessInterface.hh" 28 #include "G4BiasingProcessInterface.hh"
29 29
30 #include "G4ProcessManager.hh" 30 #include "G4ProcessManager.hh"
31 #include "G4ParticleDefinition.hh" 31 #include "G4ParticleDefinition.hh"
32 #include "G4ParticleTable.hh" 32 #include "G4ParticleTable.hh"
33 33
34 #include "G4SystemOfUnits.hh" 34 #include "G4SystemOfUnits.hh"
35 35
36 //....oooOO0OOooo........oooOO0OOooo........oo 36 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
37 37
38 G4ChannelingOptrMultiParticleChangeCrossSectio 38 G4ChannelingOptrMultiParticleChangeCrossSection::G4ChannelingOptrMultiParticleChangeCrossSection():
39 G4VBiasingOperator("ChannelingChangeXS-Many"), 39 G4VBiasingOperator("ChannelingChangeXS-Many"),
40 fCurrentOperator(0), 40 fCurrentOperator(0),
41 fnInteractions(0){ 41 fnInteractions(0){
42 AddChargedParticles(); 42 AddChargedParticles();
43 } 43 }
44 44
45 //....oooOO0OOooo........oooOO0OOooo........oo 45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
46 46
47 void G4ChannelingOptrMultiParticleChangeCrossS << 47 void G4ChannelingOptrMultiParticleChangeCrossSection::AddParticle(G4String particleName){
48 const G4ParticleDefinition* particle = 48 const G4ParticleDefinition* particle =
49 G4ParticleTable::GetParticleTable()->FindP 49 G4ParticleTable::GetParticleTable()->FindParticle( particleName );
50 50
51 if ( particle == 0 ) 51 if ( particle == 0 )
52 { 52 {
53 G4ExceptionDescription ed; 53 G4ExceptionDescription ed;
54 ed << "Particle `" << particleName << "' 54 ed << "Particle `" << particleName << "' not found !" << G4endl;
55 G4Exception("G4ChannelingOptrMultiPartic 55 G4Exception("G4ChannelingOptrMultiParticleChangeCrossSection::AddParticle(...)",
56 "G4Channeling", 56 "G4Channeling",
57 JustWarning, 57 JustWarning,
58 ed); 58 ed);
59 return; 59 return;
60 } 60 }
61 61
62 G4ChannelingOptrChangeCrossSection* optr = << 62 G4ChannelingOptrChangeCrossSection* optr = new G4ChannelingOptrChangeCrossSection(particleName);
63 new G4ChannelingOptrChangeCrossSection(par <<
64 fParticlesToBias.push_back( particle ); 63 fParticlesToBias.push_back( particle );
65 fBOptrForParticle[ particle ] = optr; 64 fBOptrForParticle[ particle ] = optr;
66 } 65 }
67 66
68 //....oooOO0OOooo........oooOO0OOooo........oo 67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69 68
70 void G4ChannelingOptrMultiParticleChangeCrossS 69 void G4ChannelingOptrMultiParticleChangeCrossSection::AddChargedParticles(){
71 G4ParticleTable::G4PTblDicIterator* aParti 70 G4ParticleTable::G4PTblDicIterator* aParticleIterator =
72 (G4ParticleTable::GetParticleTable())- 71 (G4ParticleTable::GetParticleTable())->GetIterator();
73 72
74 aParticleIterator->reset(); 73 aParticleIterator->reset();
75 74
76 while( (*aParticleIterator)() ){ 75 while( (*aParticleIterator)() ){
77 G4ParticleDefinition* particle = aPart 76 G4ParticleDefinition* particle = aParticleIterator->value();
78 if (particle->GetPDGCharge() !=0) { 77 if (particle->GetPDGCharge() !=0) {
79 AddParticle(particle->GetParticleN 78 AddParticle(particle->GetParticleName());
80 } 79 }
81 } 80 }
82 } 81 }
83 82
84 //....oooOO0OOooo........oooOO0OOooo........oo 83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
85 84
86 G4VBiasingOperation* 85 G4VBiasingOperation*
87 G4ChannelingOptrMultiParticleChangeCrossSectio 86 G4ChannelingOptrMultiParticleChangeCrossSection::
88 ProposeOccurenceBiasingOperation(const G4Track 87 ProposeOccurenceBiasingOperation(const G4Track* track,
89 const G4Biasi 88 const G4BiasingProcessInterface* callingProcess){
90 89
91 if ( fCurrentOperator ) return fCurrentOpera 90 if ( fCurrentOperator ) return fCurrentOperator->
92 GetProposedOccuren 91 GetProposedOccurenceBiasingOperation(track, callingProcess);
93 else return 0; 92 else return 0;
94 } 93 }
95 94
96 //....oooOO0OOooo........oooOO0OOooo........oo 95 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
97 96
98 void G4ChannelingOptrMultiParticleChangeCrossS 97 void G4ChannelingOptrMultiParticleChangeCrossSection::StartTracking( const G4Track* track ){
99 const G4ParticleDefinition* definition = tra 98 const G4ParticleDefinition* definition = track->GetParticleDefinition();
100 std::map < const G4ParticleDefinition*, G4Ch 99 std::map < const G4ParticleDefinition*, G4ChannelingOptrChangeCrossSection* > :: iterator
101 it = fBOptrForParticle.find( definition ); 100 it = fBOptrForParticle.find( definition );
102 fCurrentOperator = 0; 101 fCurrentOperator = 0;
103 if ( it != fBOptrForParticle.end() ) fCurren 102 if ( it != fBOptrForParticle.end() ) fCurrentOperator = (*it).second;
104 fnInteractions = 0; 103 fnInteractions = 0;
105 } 104 }
106 105
107 //....oooOO0OOooo........oooOO0OOooo........oo 106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
108 107
109 void 108 void
110 G4ChannelingOptrMultiParticleChangeCrossSectio 109 G4ChannelingOptrMultiParticleChangeCrossSection::
111 OperationApplied( const G4BiasingProcessInterf 110 OperationApplied( const G4BiasingProcessInterface* callingProcess,
112 G4BiasingAppliedCase 111 G4BiasingAppliedCase biasingCase,
113 G4VBiasingOperation* 112 G4VBiasingOperation* occurenceOperationApplied,
114 G4double 113 G4double weightForOccurenceInteraction,
115 G4VBiasingOperation* 114 G4VBiasingOperation* finalStateOperationApplied,
116 const G4VParticleChange* 115 const G4VParticleChange* particleChangeProduced ){
117 fnInteractions++; 116 fnInteractions++;
118 if ( fCurrentOperator ) fCurrentOperator->Re 117 if ( fCurrentOperator ) fCurrentOperator->ReportOperationApplied( callingProcess,
119 118 biasingCase,
120 119 occurenceOperationApplied,
121 120 weightForOccurenceInteraction,
122 121 finalStateOperationApplied,
123 122 particleChangeProduced );
124 } 123 }
125 124
126 //....oooOO0OOooo........oooOO0OOooo........oo 125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
127 126