| 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 /// \file GB02/src/GB02BOptrMultiParticleForce 26 /// \file GB02/src/GB02BOptrMultiParticleForceCollision.cc
27 /// \brief Implementation of the GB02BOptrMult 27 /// \brief Implementation of the GB02BOptrMultiParticleForceCollision class
28 // 28 //
29 #include "GB02BOptrMultiParticleForceCollision 29 #include "GB02BOptrMultiParticleForceCollision.hh"
>> 30 #include "G4BiasingProcessInterface.hh"
30 31
31 #include "G4BOptrForceCollision.hh" 32 #include "G4BOptrForceCollision.hh"
32 #include "G4BiasingProcessInterface.hh" <<
33 #include "G4ParticleDefinition.hh" 33 #include "G4ParticleDefinition.hh"
34 #include "G4ParticleTable.hh" 34 #include "G4ParticleTable.hh"
35 35
36 //....oooOO0OOooo........oooOO0OOooo........oo 36 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
37 37
38 GB02BOptrMultiParticleForceCollision::GB02BOpt 38 GB02BOptrMultiParticleForceCollision::GB02BOptrMultiParticleForceCollision()
39 : G4VBiasingOperator("TestManyForceCollision << 39 : G4VBiasingOperator("TestManyForceCollision")
40 {} << 40 {
>> 41 }
41 42
42 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 44
44 void GB02BOptrMultiParticleForceCollision::Add 45 void GB02BOptrMultiParticleForceCollision::AddParticle(G4String particleName)
45 { 46 {
46 const G4ParticleDefinition* particle = 47 const G4ParticleDefinition* particle =
47 G4ParticleTable::GetParticleTable()->FindP << 48 G4ParticleTable::GetParticleTable()->FindParticle( particleName );
48 49
49 if (particle == 0) { << 50 if ( particle == 0 )
50 G4ExceptionDescription ed; << 51 {
51 ed << "Particle `" << particleName << "' n << 52 G4ExceptionDescription ed;
52 G4Exception("GB02BOptrMultiParticleForceCo << 53 ed << "Particle `" << particleName << "' not found !" << G4endl;
53 ed); << 54 G4Exception("GB02BOptrMultiParticleForceCollision::AddParticle(...)",
54 return; << 55 "exGB02.01",
55 } << 56 JustWarning,
>> 57 ed);
>> 58 return;
>> 59 }
>> 60
>> 61 G4BOptrForceCollision* optr = new G4BOptrForceCollision(particleName,
>> 62 "ForceCollisionFor"+particleName);
>> 63 fParticlesToBias.push_back( particle );
>> 64 fBOptrForParticle[ particle ] = optr;
56 65
57 G4BOptrForceCollision* optr = <<
58 new G4BOptrForceCollision(particleName, "F <<
59 fParticlesToBias.push_back(particle); <<
60 fBOptrForParticle[particle] = optr; <<
61 } 66 }
62 67
63 //....oooOO0OOooo........oooOO0OOooo........oo 68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
64 69
65 G4VBiasingOperation* GB02BOptrMultiParticleFor << 70 G4VBiasingOperation*
66 const G4Track* track, const G4BiasingProcess << 71 GB02BOptrMultiParticleForceCollision::
>> 72 ProposeOccurenceBiasingOperation(const G4Track* track,
>> 73 const G4BiasingProcessInterface* callingProcess)
67 { 74 {
68 if (fCurrentOperator) << 75 if ( fCurrentOperator ) return fCurrentOperator->
69 return fCurrentOperator->GetProposedOccure << 76 GetProposedOccurenceBiasingOperation(track, callingProcess);
70 else << 77 else return 0;
71 return 0; <<
72 } 78 }
73 79
74 //....oooOO0OOooo........oooOO0OOooo........oo 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
75 81
76 G4VBiasingOperation* GB02BOptrMultiParticleFor << 82 G4VBiasingOperation*
77 const G4Track* track, const G4BiasingProcess << 83 GB02BOptrMultiParticleForceCollision::
>> 84 ProposeNonPhysicsBiasingOperation(const G4Track* track,
>> 85 const G4BiasingProcessInterface* callingProcess)
78 { 86 {
79 if (fCurrentOperator) << 87 if ( fCurrentOperator ) return fCurrentOperator->
80 return fCurrentOperator->GetProposedNonPhy << 88 GetProposedNonPhysicsBiasingOperation(track, callingProcess);
81 else << 89 else return 0;
82 return 0; <<
83 } 90 }
84 91
85 //....oooOO0OOooo........oooOO0OOooo........oo 92 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
86 93
87 G4VBiasingOperation* GB02BOptrMultiParticleFor << 94 G4VBiasingOperation*
88 const G4Track* track, const G4BiasingProcess << 95 GB02BOptrMultiParticleForceCollision::
>> 96 ProposeFinalStateBiasingOperation(const G4Track* track,
>> 97 const G4BiasingProcessInterface* callingProcess)
89 { 98 {
90 if (fCurrentOperator) << 99 if ( fCurrentOperator ) return fCurrentOperator->
91 return fCurrentOperator->GetProposedFinalS << 100 GetProposedFinalStateBiasingOperation(track, callingProcess);
92 else << 101 else return 0;
93 return 0; <<
94 } 102 }
95 103
96 //....oooOO0OOooo........oooOO0OOooo........oo 104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
97 105
98 void GB02BOptrMultiParticleForceCollision::Sta << 106 void GB02BOptrMultiParticleForceCollision::StartTracking( const G4Track* track )
99 { 107 {
100 const G4ParticleDefinition* definition = tra 108 const G4ParticleDefinition* definition = track->GetParticleDefinition();
101 std::map<const G4ParticleDefinition*, G4BOpt << 109 std::map < const G4ParticleDefinition*, G4BOptrForceCollision* > :: iterator
102 fBOptrForParticle.find(definition); << 110 it = fBOptrForParticle.find( definition );
103 fCurrentOperator = 0; 111 fCurrentOperator = 0;
104 if (it != fBOptrForParticle.end()) fCurrentO << 112 if ( it != fBOptrForParticle.end() ) fCurrentOperator = (*it).second;
105 } 113 }
106 114
107 //....oooOO0OOooo........oooOO0OOooo........oo 115 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
108 116
109 void GB02BOptrMultiParticleForceCollision::Ope << 117 void
110 const G4BiasingProcessInterface* callingProc << 118 GB02BOptrMultiParticleForceCollision::
111 G4VBiasingOperation* operationApplied, const << 119 OperationApplied( const G4BiasingProcessInterface* callingProcess,
>> 120 G4BiasingAppliedCase biasingCase,
>> 121 G4VBiasingOperation* operationApplied,
>> 122 const G4VParticleChange* particleChangeProduced )
112 { 123 {
113 if (fCurrentOperator) << 124 if ( fCurrentOperator ) fCurrentOperator->ReportOperationApplied( callingProcess,
114 fCurrentOperator->ReportOperationApplied(c << 125 biasingCase,
115 p << 126 operationApplied,
>> 127 particleChangeProduced );
>> 128
116 } 129 }
117 130
118 //....oooOO0OOooo........oooOO0OOooo........oo 131 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
119 132
120 void GB02BOptrMultiParticleForceCollision::Ope << 133 void
121 const G4BiasingProcessInterface* callingProc << 134 GB02BOptrMultiParticleForceCollision::
122 G4VBiasingOperation* occurenceOperationAppli << 135 OperationApplied( const G4BiasingProcessInterface* callingProcess,
123 G4VBiasingOperation* finalStateOperationAppl << 136 G4BiasingAppliedCase biasingCase,
>> 137 G4VBiasingOperation* occurenceOperationApplied,
>> 138 G4double weightForOccurenceInteraction,
>> 139 G4VBiasingOperation* finalStateOperationApplied,
>> 140 const G4VParticleChange* particleChangeProduced )
124 { 141 {
125 if (fCurrentOperator) << 142 if ( fCurrentOperator ) fCurrentOperator->ReportOperationApplied( callingProcess,
126 fCurrentOperator->ReportOperationApplied(c << 143 biasingCase,
127 w << 144 occurenceOperationApplied,
128 f << 145 weightForOccurenceInteraction,
>> 146 finalStateOperationApplied,
>> 147 particleChangeProduced );
129 } 148 }
130 149
131 //....oooOO0OOooo........oooOO0OOooo........oo 150 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
132 151
133 void GB02BOptrMultiParticleForceCollision::Exi << 152 void
134 const G4Track* track, const G4BiasingProcess << 153 GB02BOptrMultiParticleForceCollision::
>> 154 ExitBiasing( const G4Track* track,
>> 155 const G4BiasingProcessInterface* callingProcess )
135 { 156 {
136 if (fCurrentOperator) fCurrentOperator->Exit << 157 if ( fCurrentOperator ) fCurrentOperator->ExitingBiasing( track, callingProcess );
137 } 158 }
138 159
139 //....oooOO0OOooo........oooOO0OOooo........oo 160 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
140 161