ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
PrimaryGeneratorAction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PrimaryGeneratorAction.cc
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
28 //
29 //
30 //
31 //
32 
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
35 
36 #include "PrimaryGeneratorAction.hh"
37 #include "DetectorConstruction.hh"
38 
39 #include "G4Event.hh"
40 #include "G4ParticleTable.hh"
41 #include "G4IonTable.hh"
42 #include "G4ParticleDefinition.hh"
43 #include "G4Geantino.hh"
44 #include "G4PhysicalConstants.hh"
45 #include "G4SystemOfUnits.hh"
46 #include "Randomize.hh"
47 
48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49 
51 : G4VUserPrimaryGeneratorAction(), fParticleGun(0), fDetector(det)
52 {
53  G4int n_particle = 1;
54  fParticleGun = new G4ParticleGun(n_particle);
55 
59 
60 }
61 
62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
63 
65 {
66  delete fParticleGun;
67 }
68 
69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
70 
72 {
73  //particle type
75  G4int Z = 95, A = 241;
76  G4double ionCharge = 0.*eplus;
77  G4double excitEnergy = 0.*keV;
78 
80  = G4IonTable::GetIonTable()->GetIon(Z,A,excitEnergy);
82  fParticleGun->SetParticleCharge(ionCharge);
83  }
84 
85  //vertex position uniform within the absorber
86  //
88  G4double Rmax2 = Rmax*Rmax;
89  G4double R2 = G4UniformRand()*Rmax2;
90  G4double R = std::sqrt(R2);
91 
93  G4double ux = R*std::cos(phi),
94  uy = R*std::sin(phi);
95  G4double uz = (fDetector->GetAbsorLength())*(G4UniformRand() - 0.5);
96 
98 
99  //distribution uniform in solid angle
100  //
101  G4double cosTheta = 2*G4UniformRand() - 1.;
102  G4double sinTheta = std::sqrt(1. - cosTheta*cosTheta);
103  phi = twopi*G4UniformRand();
104  G4double vx = sinTheta*std::cos(phi),
105  vy = sinTheta*std::sin(phi),
106  vz = cosTheta;
107 
109 
111 }
112 
113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
114