| 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 // 26 //
27 // 27 //
28 // ------------------------------------------- 28 // ------------------------------------------------------------
29 // GEANT 4 class implementation 29 // GEANT 4 class implementation
30 // 30 //
31 // ------------- GFlashShowerModelMesseng 31 // ------------- GFlashShowerModelMessenger -------------
32 // 32 //
33 // Author: Joanna Weng - 9.11.2004 33 // Author: Joanna Weng - 9.11.2004
34 // ------------------------------------------- 34 // ------------------------------------------------------------
35 35
36 #include "GFlashShowerModelMessenger.hh" << 36 #include <iomanip>
>> 37 #include <sstream>
37 38
38 #include "G4Electron.hh" << 39 #include "GFlashShowerModelMessenger.hh"
39 #include "G4Positron.hh" << 40 #include "GFlashShowerModel.hh"
>> 41 #include "GFlashParticleBounds.hh"
>> 42 #include "globals.hh"
40 #include "G4SystemOfUnits.hh" 43 #include "G4SystemOfUnits.hh"
41 #include "G4UIcmdWithADouble.hh" << 44 #include "G4UIdirectory.hh"
42 #include "G4UIcmdWithADoubleAndUnit.hh" <<
43 #include "G4UIcmdWithAString.hh" 45 #include "G4UIcmdWithAString.hh"
>> 46 #include "G4UIcmdWithADoubleAndUnit.hh"
>> 47 #include "G4UIcmdWithADouble.hh"
44 #include "G4UIcmdWithAnInteger.hh" 48 #include "G4UIcmdWithAnInteger.hh"
45 #include "G4UIdirectory.hh" << 49 #include "G4Electron.hh"
46 #include "globals.hh" << 50 #include "G4Positron.hh"
47 <<
48 #include "GFlashParticleBounds.hh" <<
49 #include "GFlashShowerModel.hh" <<
50 <<
51 #include <iomanip> <<
52 #include <sstream> <<
53 51
54 GFlashShowerModelMessenger::GFlashShowerModelM << 52 GFlashShowerModelMessenger::
55 { << 53 GFlashShowerModelMessenger(GFlashShowerModel * aModel)
>> 54 {
56 myParaDir = new G4UIdirectory("/GFlash/"); 55 myParaDir = new G4UIdirectory("/GFlash/");
57 myParaDir->SetGuidance("Parametrisation cont 56 myParaDir->SetGuidance("Parametrisation control.");
58 myModel = aModel; << 57 myModel= aModel;
59 << 58
60 FlagCmd = new G4UIcmdWithAnInteger("/GFlash/ << 59 FlagCmd = new G4UIcmdWithAnInteger("/GFlash/flag",this);
61 FlagCmd->SetGuidance("Defines if GFlash is a 60 FlagCmd->SetGuidance("Defines if GFlash is activated");
62 FlagCmd->SetParameterName("flag", false, fal << 61 FlagCmd->SetParameterName("flag",false,false);
63 << 62
64 ContCmd = new G4UIcmdWithAnInteger("/GFlash/ << 63 ContCmd = new G4UIcmdWithAnInteger("/GFlash/containment ",this);
65 ContCmd->SetGuidance("Defines if Containment 64 ContCmd->SetGuidance("Defines if Containment is checked");
66 ContCmd->SetParameterName("flag", false, fal << 65 ContCmd->SetParameterName("flag",false,false);
67 << 66
68 StepInX0Cmd = new G4UIcmdWithADouble("/GFlas << 67 StepInX0Cmd = new G4UIcmdWithADouble("/GFlash/stepXo",this);
69 StepInX0Cmd->SetGuidance("Defines step lengh 68 StepInX0Cmd->SetGuidance("Defines step lenghts");
70 StepInX0Cmd->SetParameterName("flag", false, << 69 StepInX0Cmd->SetParameterName("flag",false,false);
71 << 70
72 EminCmd = new G4UIcmdWithADoubleAndUnit("/GF << 71 EminCmd = new G4UIcmdWithADoubleAndUnit("/GFlash/Emin",this);
73 EminCmd->SetGuidance("Set minimum kinetic en 72 EminCmd->SetGuidance("Set minimum kinetic energy to trigger parametrisation");
74 EminCmd->SetParameterName("Emin", false, fal << 73 EminCmd->SetParameterName("Emin",false,false);
75 EminCmd->SetDefaultUnit("GeV"); 74 EminCmd->SetDefaultUnit("GeV");
76 EminCmd->SetUnitCategory("Energy"); 75 EminCmd->SetUnitCategory("Energy");
77 EminCmd->AvailableForStates(G4State_PreInit, << 76 EminCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
78 << 77
79 EmaxCmd = new G4UIcmdWithADoubleAndUnit("/GF << 78 EmaxCmd = new G4UIcmdWithADoubleAndUnit("/GFlash/Emax",this);
80 EmaxCmd->SetGuidance("Set maximum kinetic en 79 EmaxCmd->SetGuidance("Set maximum kinetic energy to trigger parametrisation");
81 EmaxCmd->SetParameterName("Emax", false, fal << 80 EmaxCmd->SetParameterName("Emax",false,false);
82 EmaxCmd->SetDefaultUnit("GeV"); 81 EmaxCmd->SetDefaultUnit("GeV");
83 EmaxCmd->SetUnitCategory("Energy"); 82 EmaxCmd->SetUnitCategory("Energy");
84 EmaxCmd->AvailableForStates(G4State_PreInit, << 83 EmaxCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
85 << 84
86 EkillCmd = new G4UIcmdWithADoubleAndUnit("/G << 85 EkillCmd = new G4UIcmdWithADoubleAndUnit("/GFlash/Ekill",this);
87 EkillCmd->SetGuidance("Set maximum kinetic e 86 EkillCmd->SetGuidance("Set maximum kinetic energy for electrons to be killed");
88 EkillCmd->SetParameterName("Ekill", false, f << 87 EkillCmd->SetParameterName("Ekill",false,false);
89 EkillCmd->SetDefaultUnit("GeV"); 88 EkillCmd->SetDefaultUnit("GeV");
90 EkillCmd->SetUnitCategory("Energy"); 89 EkillCmd->SetUnitCategory("Energy");
91 EkillCmd->AvailableForStates(G4State_PreInit << 90 EkillCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
92 } 91 }
93 92
>> 93
94 GFlashShowerModelMessenger::~GFlashShowerModel 94 GFlashShowerModelMessenger::~GFlashShowerModelMessenger()
95 { 95 {
96 delete ContCmd; 96 delete ContCmd;
97 delete FlagCmd; 97 delete FlagCmd;
98 delete StepInX0Cmd; << 98 delete StepInX0Cmd;
99 delete EminCmd; 99 delete EminCmd;
100 delete EmaxCmd; 100 delete EmaxCmd;
101 delete EkillCmd; 101 delete EkillCmd;
102 } 102 }
103 103
104 void GFlashShowerModelMessenger::SetNewValue(G <<
105 { <<
106 if (command == FlagCmd) { <<
107 myModel->SetFlagParamType(FlagCmd->GetNewI <<
108 this->GetCurrentValue(command); <<
109 } <<
110 if (command == ContCmd) { <<
111 myModel->SetFlagParticleContainment(ContCm <<
112 this->GetCurrentValue(command); <<
113 } <<
114 if (command == StepInX0Cmd) { <<
115 myModel->SetStepInX0(StepInX0Cmd->GetNewDo <<
116 this->GetCurrentValue(command); <<
117 } <<
118 104
119 else if (command == EminCmd) { << 105 void GFlashShowerModelMessenger::
>> 106 SetNewValue(G4UIcommand * command,G4String newValues)
>> 107 {
>> 108
>> 109 if( command == FlagCmd ) {
>> 110 myModel->SetFlagParamType(FlagCmd->GetNewIntValue(newValues));
>> 111 this->GetCurrentValue(command);
>> 112 }
>> 113 if( command == ContCmd ) {
>> 114 myModel->SetFlagParticleContainment(ContCmd->GetNewIntValue(newValues));
>> 115 this->GetCurrentValue(command);
>> 116 }
>> 117 if( command == StepInX0Cmd ) {
>> 118 myModel->SetStepInX0(StepInX0Cmd->GetNewDoubleValue(newValues));
>> 119 this->GetCurrentValue(command);
>> 120 }
>> 121
>> 122 else if( command == EminCmd ) {
120 myModel->PBound->SetMinEneToParametrise(*G 123 myModel->PBound->SetMinEneToParametrise(*G4Electron::ElectronDefinition(),
121 Em << 124 EminCmd->GetNewDoubleValue(newValues));
122 this->GetCurrentValue(command); << 125 this->GetCurrentValue(command);
123 } 126 }
124 << 127
125 else if (command == EmaxCmd) { << 128 else if( command == EmaxCmd ) {
126 myModel->PBound->SetMaxEneToParametrise(*G 129 myModel->PBound->SetMaxEneToParametrise(*G4Electron::ElectronDefinition(),
127 Em << 130 EmaxCmd->GetNewDoubleValue(newValues));
128 this->GetCurrentValue(command); << 131 this->GetCurrentValue(command);
129 } 132 }
130 << 133
131 else if (command == EkillCmd) { << 134 else if( command == EkillCmd ) {
132 myModel->PBound->SetEneToKill(*G4Electron: 135 myModel->PBound->SetEneToKill(*G4Electron::ElectronDefinition(),
133 EkillCmd->Ge << 136 EkillCmd->GetNewDoubleValue(newValues));
134 this->GetCurrentValue(command); << 137 this->GetCurrentValue(command);
135 } 138 }
>> 139
136 } 140 }
137 141
138 G4String GFlashShowerModelMessenger::GetCurren << 142
>> 143 G4String GFlashShowerModelMessenger::GetCurrentValue(G4UIcommand * command)
139 { 144 {
140 G4String returnValue(1, '\0'); << 145 G4String returnValue(1,'\0');
141 std::ostringstream os; 146 std::ostringstream os;
142 << 147
143 if (command == FlagCmd) { << 148 if( command == FlagCmd ) {
144 os << "/GFlash/flag " << myModel->GetFlagP << 149 os << "/GFlash/flag " << myModel->GetFlagParamType() << '\0';
145 returnValue = G4String(os.str()); 150 returnValue = G4String(os.str());
146 } 151 }
147 << 152
148 else if (command == EkillCmd) { << 153 else if( command == EkillCmd ) {
149 os << "/GFlash/Ekill " << myModel->PBound- << 154 os << "/GFlash/Ekill "
>> 155 << myModel->PBound->GetEneToKill(*G4Electron::ElectronDefinition())/GeV
150 << " GeV" << '\0'; 156 << " GeV" << '\0';
151 returnValue = G4String(os.str()); 157 returnValue = G4String(os.str());
152 } 158 }
153 << 159
154 else if (command == EminCmd) { << 160 else if( command == EminCmd ) {
155 os << "/GFlash/Emin " 161 os << "/GFlash/Emin "
156 << myModel->PBound->GetMinEneToParametr << 162 << myModel->PBound->GetMinEneToParametrise(*G4Electron::ElectronDefinition())/GeV
157 << '\0'; << 163 << " GeV" << '\0';
158 returnValue = G4String(os.str()); << 164 returnValue = G4String(os.str());
159 } 165 }
160 << 166
161 else if (command == EmaxCmd) { << 167 else if( command == EmaxCmd ) {
162 os << "/GFlash/Emax " 168 os << "/GFlash/Emax "
163 << myModel->PBound->GetMaxEneToParametr << 169 << myModel->PBound->GetMaxEneToParametrise(*G4Electron::ElectronDefinition())/GeV
164 << '\0'; << 170 << " GeV" << '\0';
165 returnValue = G4String(os.str()); 171 returnValue = G4String(os.str());
166 } 172 }
167 << 173
168 return returnValue; 174 return returnValue;
169 } 175 }
170 176