| 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 // $Id: PrimaryGeneratorAction.cc 67994 2013-03-13 11:05:39Z gcosmo $
26 // 27 //
27 /// \file medical/GammaTherapy/src/PrimaryGene 28 /// \file medical/GammaTherapy/src/PrimaryGeneratorAction.cc
28 /// \brief Implementation of the PrimaryGenera 29 /// \brief Implementation of the PrimaryGeneratorAction class
29 // 30 //
30 31
31 //-------------------------------------------- 32 //---------------------------------------------------------------------------
32 // 33 //
33 // ClassName: PrimaryGeneratorAction 34 // ClassName: PrimaryGeneratorAction
34 // 35 //
35 // Description: Generate primary beam 36 // Description: Generate primary beam
36 // 37 //
37 // Authors: V.Grichine, V.Ivanchenko 38 // Authors: V.Grichine, V.Ivanchenko
38 // 39 //
39 // Modified: 40 // Modified:
40 // 41 //
41 //-------------------------------------------- 42 //----------------------------------------------------------------------------
42 // 43 //
43 44
44 #include "PrimaryGeneratorAction.hh" 45 #include "PrimaryGeneratorAction.hh"
45 <<
46 #include "DetectorConstruction.hh" 46 #include "DetectorConstruction.hh"
47 #include "PrimaryGeneratorMessenger.hh" 47 #include "PrimaryGeneratorMessenger.hh"
48 << 48 #include "Randomize.hh"
49 #include "G4ParticleDefinition.hh" <<
50 #include "G4ParticleGun.hh" 49 #include "G4ParticleGun.hh"
51 #include "G4ParticleTable.hh" 50 #include "G4ParticleTable.hh"
>> 51 #include "G4ParticleDefinition.hh"
52 #include "G4PhysicalConstants.hh" 52 #include "G4PhysicalConstants.hh"
53 #include "G4SystemOfUnits.hh" 53 #include "G4SystemOfUnits.hh"
54 #include "Randomize.hh" << 54 #include "Histo.hh"
55 55
56 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
57 //....oooOO0OOooo........oooOO0OOooo........oo 57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
58 58
59 PrimaryGeneratorAction::PrimaryGeneratorAction << 59 PrimaryGeneratorAction::PrimaryGeneratorAction(DetectorConstruction* pDet):
>> 60 fDetector(pDet)
60 { 61 {
61 InitializeMe(); 62 InitializeMe();
62 } 63 }
63 64
64 //....oooOO0OOooo........oooOO0OOooo........oo 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
65 66
66 void PrimaryGeneratorAction::InitializeMe() 67 void PrimaryGeneratorAction::InitializeMe()
67 { 68 {
68 fVerbose = fDetector->GetVerbose(); <<
69 fMessenger = new PrimaryGeneratorMessenger(t 69 fMessenger = new PrimaryGeneratorMessenger(this);
70 fParticleGun = new G4ParticleGun(); 70 fParticleGun = new G4ParticleGun();
71 fCounter = 0; 71 fCounter = 0;
>> 72 fVerbose = 0;
72 fX0 = 0.0; 73 fX0 = 0.0;
73 fY0 = 0.0; 74 fY0 = 0.0;
74 fZ0 = 0.0; 75 fZ0 = 0.0;
75 fSigmaX = 1.5 * mm; << 76 fSigmaX = 1.5*mm;
76 fSigmaY = 1.5 * mm; << 77 fSigmaY = 1.5*mm;
77 fSigmaZ = 0.0; 78 fSigmaZ = 0.0;
78 fSigmaE = 0.0; 79 fSigmaE = 0.0;
79 fRMax2 = 2.5 * 2.5 * mm * mm; << 80 fRMax2 = 2.5*2.5*mm*mm;
80 fSigmaTheta = 0.0; 81 fSigmaTheta = 0.0;
81 // fSigmaTheta = 0.17*degree; 82 // fSigmaTheta = 0.17*degree;
82 fMinCosTheta = 2.0; 83 fMinCosTheta = 2.0;
83 SetBeamEnergy(50.0 * MeV); << 84 SetBeamEnergy(50.0*MeV);
84 fPosition = G4ThreeVector(fX0, fY0, fZ0); << 85 fPosition = G4ThreeVector(fX0,fY0,fZ0);
85 fDirection = G4ThreeVector(0.0, 0.0, 1.0); << 86 fDirection = G4ThreeVector(0.0,0.0,1.0);
86 fGauss = true; 87 fGauss = true;
>> 88 if(fEnergy < (Histo::GetPointer())->GetMaxEnergy()) {
>> 89 (Histo::GetPointer())->SetMaxEnergy(fEnergy);
>> 90 }
87 } 91 }
88 92
89 //....oooOO0OOooo........oooOO0OOooo........oo 93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
90 94
91 PrimaryGeneratorAction::~PrimaryGeneratorActio 95 PrimaryGeneratorAction::~PrimaryGeneratorAction()
92 { 96 {
93 delete fParticleGun; 97 delete fParticleGun;
94 delete fMessenger; 98 delete fMessenger;
95 } 99 }
96 100
97 //....oooOO0OOooo........oooOO0OOooo........oo 101 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
98 102
99 void PrimaryGeneratorAction::GeneratePrimaries 103 void PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
100 { 104 {
101 fCounter++; << 105 fCounter++ ;
>> 106 fVerbose = (Histo::GetPointer())->GetVerbose();
102 107
103 // Simulation of beam position 108 // Simulation of beam position
104 G4double x = fX0; 109 G4double x = fX0;
105 G4double y = fY0; 110 G4double y = fY0;
106 G4double z = fDetector->GetGeneratorPosZ(); 111 G4double z = fDetector->GetGeneratorPosZ();
107 do { 112 do {
108 if (0.0 < fSigmaX) { << 113 if(0.0 < fSigmaX) { x = G4RandGauss::shoot(fX0,fSigmaX); }
109 x = G4RandGauss::shoot(fX0, fSigmaX); << 114 if(0.0 < fSigmaY) { y = G4RandGauss::shoot(fY0,fSigmaY); }
110 } << 115 } while (x*x + y*y > fRMax2);
111 if (0.0 < fSigmaY) { <<
112 y = G4RandGauss::shoot(fY0, fSigmaY); <<
113 } <<
114 } while (x * x + y * y > fRMax2); <<
115 116
116 fPosition = G4ThreeVector(x, y, z); << 117 fPosition = G4ThreeVector(x,y,z);
117 fParticleGun->SetParticlePosition(fPosition) 118 fParticleGun->SetParticlePosition(fPosition);
118 119
119 // Simulation of beam direction 120 // Simulation of beam direction
120 G4double ux = fDirection.x(); 121 G4double ux = fDirection.x();
121 G4double uy = fDirection.y(); 122 G4double uy = fDirection.y();
122 G4double uz = fDirection.z(); 123 G4double uz = fDirection.z();
123 124
124 // Beam particles are uniformly distributed 125 // Beam particles are uniformly distributed over phi, cosTheta
125 if (1.0 > fMinCosTheta) { << 126 if(1.0 > fMinCosTheta) {
126 uz = fMinCosTheta + (1.0 - fMinCosTheta) * << 127 uz = fMinCosTheta + (1.0 - fMinCosTheta)*G4UniformRand() ;
127 ux = std::sqrt((1.0 - uz) * (1.0 + uz)); << 128 ux = std::sqrt((1.0 - uz)*(1.0 + uz)) ;
128 } << 129 } else if (fSigmaTheta > 0.0) {
129 else if (fSigmaTheta > 0.0) { << 130 ux = G4RandGauss::shoot(0.0,fSigmaTheta);
130 ux = G4RandGauss::shoot(0.0, fSigmaTheta); << 131 uz = std::sqrt((1.0 - ux)*(1.0 + ux));
131 uz = std::sqrt((1.0 - ux) * (1.0 + ux)); <<
132 } 132 }
133 133
134 G4double phi = twopi * G4UniformRand(); << 134 G4double phi = twopi*G4UniformRand() ;
135 uy = ux; 135 uy = ux;
136 ux *= std::cos(phi); << 136 ux *= std::cos(phi) ;
137 uy *= std::sin(phi); << 137 uy *= std::sin(phi) ;
138 fDirection = G4ThreeVector(ux, uy, uz); << 138 fDirection = G4ThreeVector(ux,uy,uz) ;
139 139
140 fParticleGun->SetParticleMomentumDirection(f 140 fParticleGun->SetParticleMomentumDirection(fDirection);
141 141
142 // Simulation of beam kinetic energy 142 // Simulation of beam kinetic energy
143 G4double kinEnergy = fEnergy; 143 G4double kinEnergy = fEnergy;
144 144
145 if (fGauss == "flatE") { << 145 if(fGauss == "flatE") {
146 kinEnergy = fEnergy - fSigmaE + 2. * fSigm << 146 kinEnergy = fEnergy - fSigmaE + 2.*fSigmaE*G4UniformRand();
147 } << 147 } else if(0.0 < fSigmaE) {
148 else if (0.0 < fSigmaE) { << 148 kinEnergy = fEnergy + G4RandGauss::shoot(0.0,fSigmaE);
149 kinEnergy = fEnergy + G4RandGauss::shoot(0 <<
150 } 149 }
151 fParticleGun->SetParticleEnergy(kinEnergy); 150 fParticleGun->SetParticleEnergy(kinEnergy);
152 151
153 if (fVerbose > 0) { << 152 if(fVerbose > 0) {
154 G4ParticleDefinition* particle = fParticle 153 G4ParticleDefinition* particle = fParticleGun->GetParticleDefinition();
155 G4String particleName = particle->GetParti << 154 G4String particleName = particle->GetParticleName() ;
156 G4cout << "Event# " << fCounter << " Bea << 155 G4cout << "Event# " << fCounter
>> 156 << " Beam particle is generated by PrimaryGeneratorAction "
157 << G4endl; 157 << G4endl;
158 G4cout << "ParticleName= " << particleName << 158 G4cout << "ParticleName= " << particleName
159 << std::setprecision(5) << " KinE << 159 << " PDGcode= " << particle->GetPDGEncoding()
160 << " x(mm)= " << x / mm << " y(mm << 160 << std::setprecision(5)
161 << " ux= " << ux << " uy= " << uy << 161 << " KinEnergy(GeV)= "
162 } << 162 << fEnergy/GeV
>> 163 << " x(mm)= "
>> 164 << x/mm
>> 165 << " y(mm)= "
>> 166 << y/mm
>> 167 << " z(mm)= "
>> 168 << z/mm
>> 169 << " ux= "
>> 170 << ux
>> 171 << " uy= "
>> 172 << uy
>> 173 << " uz= "
>> 174 << uz
>> 175 << G4endl;
>> 176 }
163 177
164 fParticleGun->GeneratePrimaryVertex(anEvent) 178 fParticleGun->GeneratePrimaryVertex(anEvent);
165 } 179 }
166 180
167 //....oooOO0OOooo........oooOO0OOooo........oo 181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
168 182
169 void PrimaryGeneratorAction::SetBeamEnergy(G4d 183 void PrimaryGeneratorAction::SetBeamEnergy(G4double val)
170 { 184 {
171 fEnergy = val; 185 fEnergy = val;
172 if (fEnergy < fDetector->GetMaxEnergy()) fDe << 186 if(fEnergy < (Histo::GetPointer())->GetMaxEnergy()) {
>> 187 (Histo::GetPointer())->SetMaxEnergy(fEnergy);
>> 188 }
173 } 189 }
174 190
175 //....oooOO0OOooo........oooOO0OOooo........oo 191 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 192
>> 193
>> 194
176 195