| 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 "PrimaryGeneratorMessenger.hh"
27 #include "PrimaryGeneratorAction.hh"
28
29 #include "G4UIcmdWithABool.hh"
30 #include "G4UIcmdWithADouble.hh"
31 #include "G4UIcmdWithADoubleAndUnit.hh"
32 #include "G4UIcmdWithAString.hh"
33 #include "G4UIcmdWithAnInteger.hh"
34 #include "G4UIdirectory.hh"
35
36 //....oooOO0OOooo........oooOO0OOooo........oo
37
38 PrimaryGeneratorMessenger::PrimaryGeneratorMes
39 PrimaryGeneratorAction *aPrimaryGeneratorA
40 : fPrimaryGenerator(aPrimaryGeneratorActio
41 fDirectory = new G4UIdirectory("/HGCalTestbe
42 fDirectory->SetGuidance("Primary generator c
43
44 #ifdef WITHROOT
45 fReadInputCmd =
46 new G4UIcmdWithABool("/HGCalTestbeam/gen
47 fReadInputCmd->SetGuidance(
48 "If instead of particle gun, particle da
49 fReadInputCmd->SetParameterName("ifInput", t
50 fReadInputCmd->SetDefaultValue(false);
51
52 fPathInputCmd =
53 new G4UIcmdWithAString("/HGCalTestbeam/g
54 fPathInputCmd->SetGuidance("Path to input fi
55 fPathInputCmd->SetParameterName("path", true
56
57 fStartFromEventCmd =
58 new G4UIcmdWithAnInteger("/HGCalTestbeam
59 fStartFromEventCmd->SetGuidance(
60 "From which event in the file HGCalTestb
61 fStartFromEventCmd->SetParameterName("startF
62 fStartFromEventCmd->SetDefaultValue(0);
63 #endif
64 fMomentumSpreadCmd =
65 new G4UIcmdWithADouble("/HGCalTestbeam/g
66 fMomentumSpreadCmd->SetGuidance("For particl
67 fMomentumSpreadCmd->SetGuidance(
68 "Gaussian momentum spread relative to gu
69 fMomentumSpreadCmd->SetGuidance("(e.g. 0.05
70 fMomentumSpreadCmd->SetParameterName("moment
71 fMomentumSpreadCmd->SetRange("momentumSpread
72 fMomentumSpreadCmd->SetDefaultValue(0);
73
74 fBeamSpreadTypeCmd =
75 new G4UIcmdWithAString("/HGCalTestbeam/g
76 fBeamSpreadTypeCmd->SetGuidance("Type of bea
77 fBeamSpreadTypeCmd->SetParameterName("beamSp
78 fBeamSpreadTypeCmd->SetCandidates("none Gaus
79 fBeamSpreadTypeCmd->SetDefaultValue("none");
80
81 fBeamSpreadXCmd = new G4UIcmdWithADoubleAndU
82 "/HGCalTestbeam/generator/beamSpreadX",
83 fBeamSpreadXCmd->SetGuidance("Defines sigma_
84 "half-side rang
85 fBeamSpreadXCmd->SetParameterName("sigmaBeam
86 fBeamSpreadXCmd->SetRange("sigmaBeamX>=0");
87 fBeamSpreadXCmd->SetDefaultUnit("cm");
88 fBeamSpreadXCmd->SetUnitCandidates("micron m
89
90 fBeamSpreadYCmd = new G4UIcmdWithADoubleAndU
91 "/HGCalTestbeam/generator/beamSpreadY",
92 fBeamSpreadYCmd->SetGuidance("Defines sigma_
93 "half-side rang
94 fBeamSpreadYCmd->SetParameterName("sigmaBeam
95 fBeamSpreadYCmd->SetRange("sigmaBeamY>=0");
96 fBeamSpreadYCmd->SetDefaultUnit("cm");
97 fBeamSpreadYCmd->SetUnitCandidates("micron m
98
99 fBeamZ0Cmd =
100 new G4UIcmdWithADoubleAndUnit("/HGCalTes
101 fBeamZ0Cmd->SetGuidance(
102 "Beam position along the beam line (defa
103 fBeamZ0Cmd->SetParameterName("z0", true, tru
104 fBeamZ0Cmd->SetDefaultUnit("cm");
105 fBeamZ0Cmd->SetUnitCandidates("micron mm cm
106 }
107
108 //....oooOO0OOooo........oooOO0OOooo........oo
109
110 PrimaryGeneratorMessenger::~PrimaryGeneratorMe
111 delete fDirectory;
112 #ifdef WITHROOT
113 delete fReadInputCmd;
114 delete fPathInputCmd;
115 delete fStartFromEventCmd;
116 #endif
117 delete fMomentumSpreadCmd;
118 delete fBeamSpreadTypeCmd;
119 delete fBeamSpreadXCmd;
120 delete fBeamSpreadYCmd;
121 delete fBeamZ0Cmd;
122 }
123
124 //....oooOO0OOooo........oooOO0OOooo........oo
125
126 void PrimaryGeneratorMessenger::SetNewValue(G4
127 G4
128 if (command == fMomentumSpreadCmd) {
129 fPrimaryGenerator->SetMomentumSpread(
130 fMomentumSpreadCmd->GetNewDoubleValue(
131 } else if (command == fBeamSpreadTypeCmd) {
132 fPrimaryGenerator->SetBeamSpreadType(newVa
133 } else if (command == fBeamSpreadXCmd) {
134 fPrimaryGenerator->SetBeamSpreadX(
135 fBeamSpreadXCmd->GetNewDoubleValue(new
136 } else if (command == fBeamSpreadYCmd) {
137 fPrimaryGenerator->SetBeamSpreadY(
138 fBeamSpreadYCmd->GetNewDoubleValue(new
139 } else if (command == fBeamZ0Cmd) {
140 fPrimaryGenerator->SetBeamZ0(fBeamZ0Cmd->G
141 }
142 #ifdef WITHROOT
143 else if (command == fReadInputCmd) {
144 fPrimaryGenerator->SetIfUseInputFiles(
145 fReadInputCmd->GetNewBoolValue(newValu
146 } else if (command == fPathInputCmd) {
147 fPrimaryGenerator->SetInputFiles(newValues
148 } else if (command == fStartFromEventCmd) {
149 fPrimaryGenerator->SetStartFromEvent(
150 fStartFromEventCmd->GetNewIntValue(new
151 }
152 #endif
153 }
154
155 //....oooOO0OOooo........oooOO0OOooo........oo
156
157 G4String PrimaryGeneratorMessenger::GetCurrent
158 G4String cv;
159 if (command == fMomentumSpreadCmd) {
160 cv = fMomentumSpreadCmd->ConvertToString(
161 fPrimaryGenerator->GetMomentumSpread()
162 } else if (command == fBeamSpreadTypeCmd) {
163 cv = fPrimaryGenerator->GetBeamSpreadType(
164 } else if (command == fBeamSpreadXCmd) {
165 cv = fBeamSpreadXCmd->ConvertToString(fPri
166 } else if (command == fBeamSpreadYCmd) {
167 cv = fBeamSpreadYCmd->ConvertToString(fPri
168 } else if (command == fBeamZ0Cmd) {
169 cv = fBeamZ0Cmd->ConvertToString(fPrimaryG
170 }
171 #ifdef WITHROOT
172 else if (command == fReadInputCmd) {
173 cv =
174 fReadInputCmd->ConvertToString(fPrimar
175 } else if (command == fPathInputCmd) {
176 cv = fPrimaryGenerator->GetInputFiles();
177 }
178 #endif
179 return cv;
180 }
181