| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 /// \file electromagnetic/TestEm3/src/Detector << 23 // $Id: DetectorMessenger.cc,v 1.7 2004/11/23 14:05:31 maire Exp $
27 /// \brief Implementation of the DetectorMesse << 24 // GEANT4 tag $Name: geant4-07-00-patch-01 $
28 // <<
29 // 25 //
30 //....oooOO0OOooo........oooOO0OOooo........oo 26 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oo 27 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 28
33 #include "DetectorMessenger.hh" 29 #include "DetectorMessenger.hh"
34 30
35 #include "DetectorConstruction.hh" 31 #include "DetectorConstruction.hh"
36 <<
37 #include "G4UIcmdWithADoubleAndUnit.hh" <<
38 #include "G4UIcmdWithAnInteger.hh" <<
39 #include "G4UIcmdWithoutParameter.hh" <<
40 #include "G4UIcommand.hh" <<
41 #include "G4UIdirectory.hh" 32 #include "G4UIdirectory.hh"
>> 33 #include "G4UIcommand.hh"
42 #include "G4UIparameter.hh" 34 #include "G4UIparameter.hh"
43 << 35 #include "G4UIcmdWithAnInteger.hh"
44 #include <sstream> << 36 #include "G4UIcmdWithADoubleAndUnit.hh"
>> 37 #include "G4UIcmdWithoutParameter.hh"
45 38
46 //....oooOO0OOooo........oooOO0OOooo........oo 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47 40
48 DetectorMessenger::DetectorMessenger(DetectorC << 41 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det)
49 { << 42 :Detector(Det)
50 fTestemDir = new G4UIdirectory("/testem/"); << 43 {
51 fTestemDir->SetGuidance("UI commands specifi << 44 testemDir = new G4UIdirectory("/testem/");
52 << 45 testemDir->SetGuidance("UI commands specific to this example");
53 fDetDir = new G4UIdirectory("/testem/det/"); << 46
54 fDetDir->SetGuidance("detector construction << 47 detDir = new G4UIdirectory("/testem/det/");
55 << 48 detDir->SetGuidance("detector construction commands");
56 fSizeYZCmd = new G4UIcmdWithADoubleAndUnit(" << 49
57 fSizeYZCmd->SetGuidance("Set tranverse size << 50 SizeYZCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setSizeYZ",this);
58 fSizeYZCmd->SetParameterName("Size", false); << 51 SizeYZCmd->SetGuidance("Set tranverse size of the calorimeter");
59 fSizeYZCmd->SetRange("Size>0."); << 52 SizeYZCmd->SetParameterName("Size",false);
60 fSizeYZCmd->SetUnitCategory("Length"); << 53 SizeYZCmd->SetRange("Size>0.");
61 fSizeYZCmd->AvailableForStates(G4State_PreIn << 54 SizeYZCmd->SetUnitCategory("Length");
62 fSizeYZCmd->SetToBeBroadcasted(false); << 55 SizeYZCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
63 << 56
64 fNbLayersCmd = new G4UIcmdWithAnInteger("/te << 57 NbLayersCmd = new G4UIcmdWithAnInteger("/testem/det/setNbOfLayers",this);
65 fNbLayersCmd->SetGuidance("Set number of lay << 58 NbLayersCmd->SetGuidance("Set number of layers.");
66 fNbLayersCmd->SetParameterName("NbLayers", f << 59 NbLayersCmd->SetParameterName("NbLayers",false);
67 fNbLayersCmd->SetRange("NbLayers>0"); << 60 NbLayersCmd->SetRange("NbLayers>0 && NbLayers<500");
68 fNbLayersCmd->AvailableForStates(G4State_Pre << 61 NbLayersCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
69 fNbLayersCmd->SetToBeBroadcasted(false); << 62
70 << 63 NbAbsorCmd = new G4UIcmdWithAnInteger("/testem/det/setNbOfAbsor",this);
71 fNbAbsorCmd = new G4UIcmdWithAnInteger("/tes << 64 NbAbsorCmd->SetGuidance("Set number of Absorbers.");
72 fNbAbsorCmd->SetGuidance("Set number of Abso << 65 NbAbsorCmd->SetParameterName("NbAbsor",false);
73 fNbAbsorCmd->SetParameterName("NbAbsor", fal << 66 NbAbsorCmd->SetRange("NbAbsor>0");
74 fNbAbsorCmd->SetRange("NbAbsor>0"); << 67 NbAbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
75 fNbAbsorCmd->AvailableForStates(G4State_PreI << 68
76 fNbAbsorCmd->SetToBeBroadcasted(false); << 69 AbsorCmd = new G4UIcommand("/testem/det/setAbsor",this);
77 << 70 AbsorCmd->SetGuidance("Set the absor nb, the material, the thickness.");
78 fAbsorCmd = new G4UIcommand("/testem/det/set << 71 AbsorCmd->SetGuidance(" absor number : from 1 to NbOfAbsor");
79 fAbsorCmd->SetGuidance("Set the absor nb, th << 72 AbsorCmd->SetGuidance(" material name");
80 fAbsorCmd->SetGuidance(" absor number : fro << 73 AbsorCmd->SetGuidance(" thickness (with unit) : t>0.");
81 fAbsorCmd->SetGuidance(" material name"); <<
82 fAbsorCmd->SetGuidance(" thickness (with un <<
83 // 74 //
84 G4UIparameter* AbsNbPrm = new G4UIparameter( << 75 G4UIparameter* AbsNbPrm = new G4UIparameter("AbsorNb",'i',false);
85 AbsNbPrm->SetGuidance("absor number : from 1 76 AbsNbPrm->SetGuidance("absor number : from 1 to NbOfAbsor");
86 AbsNbPrm->SetParameterRange("AbsorNb>0"); 77 AbsNbPrm->SetParameterRange("AbsorNb>0");
87 fAbsorCmd->SetParameter(AbsNbPrm); << 78 AbsorCmd->SetParameter(AbsNbPrm);
88 // 79 //
89 G4UIparameter* MatPrm = new G4UIparameter("m << 80 G4UIparameter* MatPrm = new G4UIparameter("material",'s',false);
90 MatPrm->SetGuidance("material name"); 81 MatPrm->SetGuidance("material name");
91 fAbsorCmd->SetParameter(MatPrm); << 82 AbsorCmd->SetParameter(MatPrm);
92 // << 83 //
93 G4UIparameter* ThickPrm = new G4UIparameter( << 84 G4UIparameter* ThickPrm = new G4UIparameter("thickness",'d',false);
94 ThickPrm->SetGuidance("thickness of absorber 85 ThickPrm->SetGuidance("thickness of absorber");
95 ThickPrm->SetParameterRange("thickness>0."); 86 ThickPrm->SetParameterRange("thickness>0.");
96 fAbsorCmd->SetParameter(ThickPrm); << 87 AbsorCmd->SetParameter(ThickPrm);
97 // 88 //
98 G4UIparameter* unitPrm = new G4UIparameter(" << 89 G4UIparameter* unitPrm = new G4UIparameter("unit",'s',false);
99 unitPrm->SetGuidance("unit of thickness"); 90 unitPrm->SetGuidance("unit of thickness");
100 G4String unitList = G4UIcommand::UnitsList(G 91 G4String unitList = G4UIcommand::UnitsList(G4UIcommand::CategoryOf("mm"));
101 unitPrm->SetParameterCandidates(unitList); 92 unitPrm->SetParameterCandidates(unitList);
102 fAbsorCmd->SetParameter(unitPrm); << 93 AbsorCmd->SetParameter(unitPrm);
103 // 94 //
104 fAbsorCmd->AvailableForStates(G4State_PreIni << 95 AbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
105 fAbsorCmd->SetToBeBroadcasted(false); << 96
>> 97 MagFieldCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setField",this);
>> 98 MagFieldCmd->SetGuidance("Define magnetic field.");
>> 99 MagFieldCmd->SetGuidance("Magnetic field will be in Z direction.");
>> 100 MagFieldCmd->SetParameterName("Bz",false);
>> 101 MagFieldCmd->SetUnitCategory("Magnetic flux density");
>> 102 MagFieldCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
>> 103
>> 104 UpdateCmd = new G4UIcmdWithoutParameter("/testem/det/update",this);
>> 105 UpdateCmd->SetGuidance("Update calorimeter geometry.");
>> 106 UpdateCmd->SetGuidance("This command MUST be applied before \"beamOn\" ");
>> 107 UpdateCmd->SetGuidance("if you changed geometrical value(s).");
>> 108 UpdateCmd->AvailableForStates(G4State_Idle);
106 } 109 }
107 110
108 //....oooOO0OOooo........oooOO0OOooo........oo 111 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
109 112
110 DetectorMessenger::~DetectorMessenger() 113 DetectorMessenger::~DetectorMessenger()
111 { 114 {
112 delete fSizeYZCmd; << 115 delete SizeYZCmd;
113 delete fNbLayersCmd; << 116 delete NbLayersCmd;
114 delete fNbAbsorCmd; << 117 delete NbAbsorCmd;
115 delete fAbsorCmd; << 118 delete AbsorCmd;
116 delete fDetDir; << 119 delete MagFieldCmd;
117 delete fTestemDir; << 120 delete UpdateCmd;
>> 121 delete detDir;
>> 122 delete testemDir;
118 } 123 }
119 124
120 //....oooOO0OOooo........oooOO0OOooo........oo 125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
121 126
122 void DetectorMessenger::SetNewValue(G4UIcomman << 127 void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
123 { 128 {
124 if (command == fSizeYZCmd) { << 129 if( command == SizeYZCmd )
125 fDetector->SetCalorSizeYZ(fSizeYZCmd->GetN << 130 { Detector->SetCalorSizeYZ(SizeYZCmd->GetNewDoubleValue(newValue));}
126 } << 131
127 << 132 if( command == NbLayersCmd )
128 if (command == fNbLayersCmd) { << 133 { Detector->SetNbOfLayers(NbLayersCmd->GetNewIntValue(newValue));}
129 fDetector->SetNbOfLayers(fNbLayersCmd->Get << 134
130 } << 135 if( command == NbAbsorCmd )
131 << 136 { Detector->SetNbOfAbsor(NbAbsorCmd->GetNewIntValue(newValue));}
132 if (command == fNbAbsorCmd) { << 137
133 fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNe << 138 if (command == AbsorCmd)
134 } << 139 {
135 << 140 G4int num; G4double tick;
136 if (command == fAbsorCmd) { << 141 char mat[30],unts[30];
137 G4int num; << 142 const char* t = newValue;
138 G4double tick; << 143 std::istrstream is((char*)t);
139 G4String unt, mat; << 144 is >> num >> mat >> tick >> unts;
140 std::istringstream is(newValue); << 145 G4String material=mat, unt=unts;
141 is >> num >> mat >> tick >> unt; << 146 tick *= G4UIcommand::ValueOf(unt);
142 G4String material = mat; << 147 Detector->SetAbsorMaterial (num,material);
143 tick *= G4UIcommand::ValueOf(unt); << 148 Detector->SetAbsorThickness(num,tick);
144 fDetector->SetAbsorMaterial(num, material) << 149 }
145 fDetector->SetAbsorThickness(num, tick); << 150
146 } << 151 if( command == MagFieldCmd )
>> 152 { Detector->SetMagField(MagFieldCmd->GetNewDoubleValue(newValue));}
>> 153
>> 154 if( command == UpdateCmd )
>> 155 { Detector->UpdateGeometry();}
147 } 156 }
148 157
149 //....oooOO0OOooo........oooOO0OOooo........oo 158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
150 159