| 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 /// \file hadronic/Hadr01/src/PhysicsListMesse 26 /// \file hadronic/Hadr01/src/PhysicsListMessenger.cc
27 /// \brief Implementation of the PhysicsListMe 27 /// \brief Implementation of the PhysicsListMessenger class
28 // 28 //
29 // 29 //
>> 30 // $Id: PhysicsListMessenger.cc 70761 2013-06-05 12:30:51Z gcosmo $
>> 31 //
>> 32 //
30 ////////////////////////////////////////////// 33 /////////////////////////////////////////////////////////////////////////
31 // 34 //
32 // PhysicsListMessenger 35 // PhysicsListMessenger
33 // 36 //
34 // Created: 31.01.2006 V.Ivanchenko 37 // Created: 31.01.2006 V.Ivanchenko
35 // 38 //
36 // Modified: 39 // Modified:
37 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchen 40 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchenko)
38 // 41 //
39 ////////////////////////////////////////////// 42 ////////////////////////////////////////////////////////////////////////
40 // 43 //
41 // << 44 //
42 45
43 #include "PhysicsListMessenger.hh" 46 #include "PhysicsListMessenger.hh"
44 47
45 #include "PhysicsList.hh" 48 #include "PhysicsList.hh"
46 << 49 #include "G4UIcmdWithADoubleAndUnit.hh"
47 #include "G4UIcmdWithAString.hh" 50 #include "G4UIcmdWithAString.hh"
48 #include "G4UIcmdWithoutParameter.hh" 51 #include "G4UIcmdWithoutParameter.hh"
>> 52 #include "G4UImanager.hh"
49 53
50 //....oooOO0OOooo........oooOO0OOooo........oo 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
51 55
52 PhysicsListMessenger::PhysicsListMessenger(Phy 56 PhysicsListMessenger::PhysicsListMessenger(PhysicsList* pPhys)
53 : G4UImessenger(), fPhysicsList(pPhys) << 57 :G4UImessenger(), fPhysicsList(pPhys),
54 { << 58 fGammaCutCmd(0), fElectCutCmd(0), fPosCutCmd(0), fCutCmd(0), fAllCutCmd(0),
55 fPListCmd = new G4UIcmdWithAString("/testhad << 59 fPListCmd(0), fListCmd(0)
56 fPListCmd->SetGuidance("Add modular physics << 60 {
57 fPListCmd->SetParameterName("PList", false); << 61 fGammaCutCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/CutGamma",this);
>> 62 fGammaCutCmd->SetGuidance("Set gamma cut.");
>> 63 fGammaCutCmd->SetParameterName("Gcut",false);
>> 64 fGammaCutCmd->SetUnitCategory("Length");
>> 65 fGammaCutCmd->SetRange("Gcut>=0.0");
>> 66 fGammaCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 67
>> 68 fElectCutCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/CutEl",this);
>> 69 fElectCutCmd->SetGuidance("Set electron cut.");
>> 70 fElectCutCmd->SetParameterName("Ecut",false);
>> 71 fElectCutCmd->SetUnitCategory("Length");
>> 72 fElectCutCmd->SetRange("Ecut>=0.0");
>> 73 fElectCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 74
>> 75 fPosCutCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/CutPos",this);
>> 76 fPosCutCmd->SetGuidance("Set positron cut.");
>> 77 fPosCutCmd->SetParameterName("Pcut",false);
>> 78 fPosCutCmd->SetUnitCategory("Length");
>> 79 fPosCutCmd->SetRange("Pcut>=0.0");
>> 80 fPosCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 81
>> 82 fCutCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/CutProt",this);
>> 83 fCutCmd->SetGuidance("Set proton cut.");
>> 84 fCutCmd->SetParameterName("ProtCut",false);
>> 85 fCutCmd->SetUnitCategory("Length");
>> 86 fCutCmd->SetRange("ProtCut>=0.0");
>> 87 fCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 88
>> 89 fAllCutCmd = new G4UIcmdWithADoubleAndUnit("/testhadr/CutsAll",this);
>> 90 fAllCutCmd->SetGuidance("Set cut for all.");
>> 91 fAllCutCmd->SetParameterName("cut",false);
>> 92 fAllCutCmd->SetUnitCategory("Length");
>> 93 fAllCutCmd->SetRange("cut>=0.0");
>> 94 fAllCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 95
>> 96 fPListCmd = new G4UIcmdWithAString("/testhadr/Physics",this);
>> 97 fPListCmd->SetGuidance("Add modula physics list.");
>> 98 fPListCmd->SetParameterName("PList",false);
58 fPListCmd->AvailableForStates(G4State_PreIni 99 fPListCmd->AvailableForStates(G4State_PreInit);
59 100
60 fListCmd = new G4UIcmdWithoutParameter("/tes << 101 fListCmd = new G4UIcmdWithoutParameter("/testhadr/ListPhysics",this);
61 fListCmd->SetGuidance("Available Physics Lis 102 fListCmd->SetGuidance("Available Physics Lists");
62 fListCmd->AvailableForStates(G4State_PreInit << 103 fListCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 104
63 } 105 }
64 106
65 //....oooOO0OOooo........oooOO0OOooo........oo 107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66 108
67 PhysicsListMessenger::~PhysicsListMessenger() 109 PhysicsListMessenger::~PhysicsListMessenger()
68 { 110 {
>> 111 delete fGammaCutCmd;
>> 112 delete fElectCutCmd;
>> 113 delete fPosCutCmd;
>> 114 delete fCutCmd;
>> 115 delete fAllCutCmd;
69 delete fPListCmd; 116 delete fPListCmd;
70 delete fListCmd; 117 delete fListCmd;
71 } 118 }
72 119
73 //....oooOO0OOooo........oooOO0OOooo........oo 120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
74 121
75 void PhysicsListMessenger::SetNewValue(G4UIcom 122 void PhysicsListMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
76 { 123 {
77 if (command == fPListCmd) { << 124 G4UImanager* UI = G4UImanager::GetUIpointer();
78 if (fPhysicsList) { << 125 if( command == fGammaCutCmd ) {
>> 126 if(fPhysicsList) {
>> 127 fPhysicsList->SetCutForGamma(fGammaCutCmd->GetNewDoubleValue(newValue));
>> 128 } else {
>> 129 UI->ApplyCommand("/run/setCutForAGivenParticle gamma " + newValue);
>> 130 }
>> 131
>> 132 } else if( command == fElectCutCmd ) {
>> 133 if(fPhysicsList) {
>> 134 fPhysicsList->SetCutForElectron(
>> 135 fElectCutCmd->GetNewDoubleValue(newValue));
>> 136 } else {
>> 137 UI->ApplyCommand("/run/setCutForAGivenParticle e- " + newValue);
>> 138 }
>> 139
>> 140 } else if( command == fPosCutCmd ) {
>> 141 if(fPhysicsList) {
>> 142 fPhysicsList->SetCutForPositron(fPosCutCmd->GetNewDoubleValue(newValue));
>> 143 } else {
>> 144 UI->ApplyCommand("/run/setCutForAGivenParticle e+ " + newValue);
>> 145 }
>> 146
>> 147 } else if( command == fCutCmd ) {
>> 148 if(fPhysicsList) {
>> 149 fPhysicsList->SetCutForProton(fCutCmd->GetNewDoubleValue(newValue));
>> 150 } else {
>> 151 UI->ApplyCommand("/run/setCutForAGivenParticle proton " + newValue);
>> 152 }
>> 153
>> 154 } else if( command == fAllCutCmd ) {
>> 155
>> 156 if(fPhysicsList) {
>> 157 G4double cut = fAllCutCmd->GetNewDoubleValue(newValue);
>> 158 fPhysicsList->SetCutForGamma(cut);
>> 159 fPhysicsList->SetCutForElectron(cut);
>> 160 fPhysicsList->SetCutForPositron(cut);
>> 161 fPhysicsList->SetCutForProton(cut);
>> 162 } else {
>> 163 UI->ApplyCommand("/run/setCut " + newValue);
>> 164 }
>> 165
>> 166 } else if( command == fPListCmd ) {
>> 167 if(fPhysicsList) {
79 G4String name = newValue; 168 G4String name = newValue;
80 if (name == "PHYSLIST") { << 169 if(name == "PHYSLIST") {
81 char* path = std::getenv(name); << 170 char* path = getenv(name);
82 if (path) << 171 if (path) name = G4String(path);
83 name = G4String(path); <<
84 else { 172 else {
85 G4cout << "### PhysicsListMessenger 173 G4cout << "### PhysicsListMessenger WARNING: "
86 << " environment variable PHY << 174 << " environment variable PHYSLIST is not defined"
87 return; << 175 << G4endl;
>> 176 return;
88 } 177 }
89 } 178 }
90 fPhysicsList->AddPhysicsList(name); 179 fPhysicsList->AddPhysicsList(name);
91 } << 180 } else {
92 else { <<
93 G4cout << "### PhysicsListMessenger WARN 181 G4cout << "### PhysicsListMessenger WARNING: "
94 << " /testhadr/Physics UI command 182 << " /testhadr/Physics UI command is not available "
95 << "for reference Physics List" < 183 << "for reference Physics List" << G4endl;
96 } 184 }
97 } << 185
98 else if (command == fListCmd) { << 186 } else if( command == fListCmd ) {
99 if (fPhysicsList) { << 187 if(fPhysicsList) {
100 fPhysicsList->List(); 188 fPhysicsList->List();
101 } << 189 } else {
102 else { <<
103 G4cout << "### PhysicsListMessenger WARN 190 G4cout << "### PhysicsListMessenger WARNING: "
104 << " /testhadr/ListPhysics UI com 191 << " /testhadr/ListPhysics UI command is not available "
105 << "for reference Physics List" < 192 << "for reference Physics List" << G4endl;
106 } 193 }
107 } 194 }
108 } 195 }
109 196
110 //....oooOO0OOooo........oooOO0OOooo........oo 197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
111 198