| 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 //
27 /// \file XrayTESdetPrimaryGeneratorAction.cc
28 /// \brief Implementation of the PrimaryGenera
29 //
30 // Authors: P.Dondero (paolo.dondero@cern.ch),
31 //
32 //....oooOO0OOooo........oooOO0OOooo........oo
33 //....oooOO0OOooo........oooOO0OOooo........oo
34
35
36 #include "XrayTESdetPrimaryGeneratorAction.hh"
37 #include "XrayTESdetDetectorConstruction.hh"
38
39 #include "G4Event.hh"
40 #include "G4ParticleGun.hh"
41 #include "G4ThreeVector.hh"
42 #include "G4ParticleTable.hh"
43 #include "G4ParticleDefinition.hh"
44 #include "Randomize.hh"
45 #include "globals.hh"
46 #include "G4NistManager.hh"
47 #include "G4SystemOfUnits.hh"
48 #include "G4PhysicalConstants.hh"
49
50 //....oooOO0OOooo........oooOO0OOooo........oo
51
52 XrayTESdetPrimaryGeneratorAction::XrayTESdetPr
53 {
54 G4int n_particle = 1;
55 fParticleGun = new G4ParticleGun(n_particle)
56 G4ParticleTable* particleTable = G4ParticleT
57 G4String particleName;
58
59 // setting particle type
60 fParticleGun->SetParticleDefinition(particl
61 }
62
63 //....oooOO0OOooo........oooOO0OOooo........oo
64
65 XrayTESdetPrimaryGeneratorAction::~XrayTESdetP
66 {
67 delete fParticleGun;
68 }
69
70 //....oooOO0OOooo........oooOO0OOooo........oo
71
72 void XrayTESdetPrimaryGeneratorAction::Generat
73 {
74 G4double ax;
75 G4double ay;
76 G4double az;
77
78 G4double bx;
79 G4double by;
80 G4double bz;
81
82 G4double vx;
83 G4double vy;
84 G4double vz;
85
86 G4double energy = 0;
87 G4double u;
88 G4double u1;
89
90 G4double theta0;
91 G4double phi0;
92 G4double theta1;
93 G4double phi1;
94
95 G4double ene;
96 G4double flux;
97 G4double r;
98
99 // Radius to shot particles
100 r = 27*cm;
101
102 // Proton flux @L2 (Pamela2009), in p/cm2/s/
103 G4double phi = 0.379;
104 G4double tr = 0.938;
105 G4double T;
106
107 while (true)
108 {
109 u = G4UniformRand();
110 u1 = 0.268*G4UniformRand(); //max of
111 ene = 10.*MeV + u*100000.*MeV;
112 T = ene*0.001;
113 flux = 1.9*std::pow(((T + phi)*(T + phi +
114 if (u1 <= flux)
115 {
116 energy = ene;
117 break;
118 }
119 }
120
121 theta0 = std::acos(1. - 2.*G4UniformRand());
122 theta1 = std::acos(1. - 2.*G4UniformRand());
123 phi0 = twopi*G4UniformRand();
124 phi1 = twopi*G4UniformRand();
125 ax = r*std::sin(theta0)*std::cos(phi0);
126 ay = r*std::sin(theta0)*std::sin(phi0);
127 az = r*std::cos(theta0) + 80*mm;
128
129 bx = r*std::sin(theta1)*std::cos(phi1);
130 by = r*std::sin(theta1)*std::sin(phi1);
131 bz = r*std::cos(theta1) + 80*mm;
132 vx = bx - ax;
133 vy = by - ay;
134 vz = bz - az;
135
136 fParticleGun->SetParticleEnergy(energy);
137 fParticleGun->SetParticlePosition(G4ThreeVec
138 fParticleGun->SetParticleMomentumDirection(G
139 fParticleGun->GeneratePrimaryVertex(anEvent)
140 }
141
142 //....oooOO0OOooo........oooOO0OOooo........oo
143