| 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 B5/include/PrimaryGeneratorAction.hh
28 /// \brief Definition of the B5::PrimaryGenera
29
30 #ifndef B5PrimaryGeneratorAction_h
31 #define B5PrimaryGeneratorAction_h 1
32
33 #include "G4VUserPrimaryGeneratorAction.hh"
34 #include "globals.hh"
35
36 #include <CLHEP/Units/SystemOfUnits.h>
37
38 class G4ParticleGun;
39 class G4GenericMessenger;
40 class G4Event;
41 class G4ParticleDefinition;
42
43 namespace B5
44 {
45
46 /// Primary generator
47 ///
48 /// A single particle is generated.
49 /// User can select
50 /// - the initial momentum and angle
51 /// - the momentum and angle spreads
52 /// - random selection of a particle type from
53
54 class PrimaryGeneratorAction : public G4VUserP
55 {
56 public:
57 PrimaryGeneratorAction();
58 ~PrimaryGeneratorAction() override;
59
60 void GeneratePrimaries(G4Event*) override;
61
62 void SetMomentum(G4double val) { fMomentum
63 G4double GetMomentum() const { return fMom
64
65 void SetSigmaMomentum(G4double val) { fSig
66 G4double GetSigmaMomentum() const { return
67
68 void SetSigmaAngle(G4double val) { fSigmaA
69 G4double GetSigmaAngle() const { return fS
70
71 void SetRandomize(G4bool val) { fRandomize
72 G4bool GetRandomize() const { return fRand
73
74 private:
75 void DefineCommands();
76
77 G4ParticleGun* fParticleGun = nullptr;
78 G4GenericMessenger* fMessenger = nullptr;
79 G4ParticleDefinition* fPositron = nullptr;
80 G4ParticleDefinition* fMuon = nullptr;
81 G4ParticleDefinition* fPion = nullptr;
82 G4ParticleDefinition* fKaon = nullptr;
83 G4ParticleDefinition* fProton = nullptr;
84 G4double fMomentum = 1000. * CLHEP::MeV;
85 G4double fSigmaMomentum = 50. * CLHEP::MeV
86 G4double fSigmaAngle = 2. * CLHEP::deg;
87 G4bool fRandomizePrimary = true;
88 };
89
90 } // namespace B5
91
92 //....oooOO0OOooo........oooOO0OOooo........oo
93
94 #endif
95