| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 #include "G4ChannelingOptrMultiParticleChangeC
27 #include "G4ChannelingOptrChangeCrossSection.h
28 #include "G4BiasingProcessInterface.hh"
29
30 #include "G4ProcessManager.hh"
31 #include "G4ParticleDefinition.hh"
32 #include "G4ParticleTable.hh"
33
34 #include "G4SystemOfUnits.hh"
35
36 //....oooOO0OOooo........oooOO0OOooo........oo
37
38 G4ChannelingOptrMultiParticleChangeCrossSectio
39 G4VBiasingOperator("ChannelingChangeXS-Many"),
40 fCurrentOperator(0),
41 fnInteractions(0){
42 AddChargedParticles();
43 }
44
45 //....oooOO0OOooo........oooOO0OOooo........oo
46
47 void G4ChannelingOptrMultiParticleChangeCrossS
48 const G4ParticleDefinition* particle =
49 G4ParticleTable::GetParticleTable()->FindP
50
51 if ( particle == 0 )
52 {
53 G4ExceptionDescription ed;
54 ed << "Particle `" << particleName << "'
55 G4Exception("G4ChannelingOptrMultiPartic
56 "G4Channeling",
57 JustWarning,
58 ed);
59 return;
60 }
61
62 G4ChannelingOptrChangeCrossSection* optr =
63 new G4ChannelingOptrChangeCrossSection(par
64 fParticlesToBias.push_back( particle );
65 fBOptrForParticle[ particle ] = optr;
66 }
67
68 //....oooOO0OOooo........oooOO0OOooo........oo
69
70 void G4ChannelingOptrMultiParticleChangeCrossS
71 G4ParticleTable::G4PTblDicIterator* aParti
72 (G4ParticleTable::GetParticleTable())-
73
74 aParticleIterator->reset();
75
76 while( (*aParticleIterator)() ){
77 G4ParticleDefinition* particle = aPart
78 if (particle->GetPDGCharge() !=0) {
79 AddParticle(particle->GetParticleN
80 }
81 }
82 }
83
84 //....oooOO0OOooo........oooOO0OOooo........oo
85
86 G4VBiasingOperation*
87 G4ChannelingOptrMultiParticleChangeCrossSectio
88 ProposeOccurenceBiasingOperation(const G4Track
89 const G4Biasi
90
91 if ( fCurrentOperator ) return fCurrentOpera
92 GetProposedOccuren
93 else return 0;
94 }
95
96 //....oooOO0OOooo........oooOO0OOooo........oo
97
98 void G4ChannelingOptrMultiParticleChangeCrossS
99 const G4ParticleDefinition* definition = tra
100 std::map < const G4ParticleDefinition*, G4Ch
101 it = fBOptrForParticle.find( definition );
102 fCurrentOperator = 0;
103 if ( it != fBOptrForParticle.end() ) fCurren
104 fnInteractions = 0;
105 }
106
107 //....oooOO0OOooo........oooOO0OOooo........oo
108
109 void
110 G4ChannelingOptrMultiParticleChangeCrossSectio
111 OperationApplied( const G4BiasingProcessInterf
112 G4BiasingAppliedCase
113 G4VBiasingOperation*
114 G4double
115 G4VBiasingOperation*
116 const G4VParticleChange*
117 fnInteractions++;
118 if ( fCurrentOperator ) fCurrentOperator->Re
119
120
121
122
123
124 }
125
126 //....oooOO0OOooo........oooOO0OOooo........oo
127