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