Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm11/src/DetectorMessenger.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/electromagnetic/TestEm11/src/DetectorMessenger.cc (Version 11.3.0) and /examples/extended/electromagnetic/TestEm11/src/DetectorMessenger.cc (Version 9.6)


  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 electromagnetic/TestEm11/src/Detecto     26 /// \file electromagnetic/TestEm11/src/DetectorMessenger.cc
 27 /// \brief Implementation of the DetectorMesse     27 /// \brief Implementation of the DetectorMessenger class
 28 //                                                 28 //
                                                   >>  29 // $Id$
 29 //                                                 30 //
 30 //....oooOO0OOooo........oooOO0OOooo........oo     31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 31 //....oooOO0OOooo........oooOO0OOooo........oo     32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32                                                    33 
 33 #include "DetectorMessenger.hh"                    34 #include "DetectorMessenger.hh"
 34                                                    35 
 35 #include "DetectorConstruction.hh"             <<  36 #include <sstream>
 36                                                    37 
 37 #include "G4UIcmdWithADoubleAndUnit.hh"        <<  38 #include "DetectorConstruction.hh"
 38 #include "G4UIcmdWithAnInteger.hh"             << 
 39 #include "G4UIcmdWithoutParameter.hh"          << 
 40 #include "G4UIcommand.hh"                      << 
 41 #include "G4UIdirectory.hh"                        39 #include "G4UIdirectory.hh"
                                                   >>  40 #include "G4UIcommand.hh"
 42 #include "G4UIparameter.hh"                        41 #include "G4UIparameter.hh"
 43                                                <<  42 #include "G4UIcmdWithAnInteger.hh"
 44 #include <sstream>                             <<  43 #include "G4UIcmdWithADoubleAndUnit.hh"
                                                   >>  44 #include "G4UIcmdWithoutParameter.hh"
 45                                                    45 
 46 //....oooOO0OOooo........oooOO0OOooo........oo     46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 47                                                    47 
 48 DetectorMessenger::DetectorMessenger(DetectorC <<  48 DetectorMessenger::DetectorMessenger(DetectorConstruction * Det)
 49 {                                              <<  49 :fDetector(Det)
                                                   >>  50 { 
 50   fTestemDir = new G4UIdirectory("/testem/");      51   fTestemDir = new G4UIdirectory("/testem/");
 51   fTestemDir->SetGuidance(" detector control."     52   fTestemDir->SetGuidance(" detector control.");
 52                                                <<  53   
 53   fDetDir = new G4UIdirectory("/testem/det/");     54   fDetDir = new G4UIdirectory("/testem/det/");
 54   fDetDir->SetGuidance("detector construction      55   fDetDir->SetGuidance("detector construction commands");
 55                                                <<  56   
 56   fNbAbsorCmd = new G4UIcmdWithAnInteger("/tes <<  57   fNbAbsorCmd = new G4UIcmdWithAnInteger("/testem/det/setNbOfAbsor",this);
 57   fNbAbsorCmd->SetGuidance("Set number of Abso     58   fNbAbsorCmd->SetGuidance("Set number of Absorbers.");
 58   fNbAbsorCmd->SetParameterName("NbAbsor", fal <<  59   fNbAbsorCmd->SetParameterName("NbAbsor",false);
 59   fNbAbsorCmd->SetRange("NbAbsor>0");              60   fNbAbsorCmd->SetRange("NbAbsor>0");
 60   fNbAbsorCmd->AvailableForStates(G4State_PreI <<  61   fNbAbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
 61   fNbAbsorCmd->SetToBeBroadcasted(false);      <<  62    
 62                                                <<  63   fAbsorCmd = new G4UIcommand("/testem/det/setAbsor",this);
 63   fAbsorCmd = new G4UIcommand("/testem/det/set << 
 64   fAbsorCmd->SetGuidance("Set the absor nb, th     64   fAbsorCmd->SetGuidance("Set the absor nb, the material, the thickness.");
 65   fAbsorCmd->SetGuidance("  absor number : fro     65   fAbsorCmd->SetGuidance("  absor number : from 1 to NbOfAbsor");
 66   fAbsorCmd->SetGuidance("  material name");       66   fAbsorCmd->SetGuidance("  material name");
 67   fAbsorCmd->SetGuidance("  thickness (with un     67   fAbsorCmd->SetGuidance("  thickness (with unit) : t>0.");
 68   //                                               68   //
 69   G4UIparameter* AbsNbPrm = new G4UIparameter( <<  69   G4UIparameter* AbsNbPrm = new G4UIparameter("AbsorNb",'i',false);
 70   AbsNbPrm->SetGuidance("absor number : from 1     70   AbsNbPrm->SetGuidance("absor number : from 1 to NbOfAbsor");
 71   AbsNbPrm->SetParameterRange("AbsorNb>0");        71   AbsNbPrm->SetParameterRange("AbsorNb>0");
 72   fAbsorCmd->SetParameter(AbsNbPrm);               72   fAbsorCmd->SetParameter(AbsNbPrm);
 73   //                                               73   //
 74   G4UIparameter* MatPrm = new G4UIparameter("m <<  74   G4UIparameter* MatPrm = new G4UIparameter("material",'s',false);
 75   MatPrm->SetGuidance("material name");            75   MatPrm->SetGuidance("material name");
 76   fAbsorCmd->SetParameter(MatPrm);                 76   fAbsorCmd->SetParameter(MatPrm);
 77   //                                           <<  77   //    
 78   G4UIparameter* ThickPrm = new G4UIparameter( <<  78   G4UIparameter* ThickPrm = new G4UIparameter("thickness",'d',false);
 79   ThickPrm->SetGuidance("thickness of absorber     79   ThickPrm->SetGuidance("thickness of absorber");
 80   ThickPrm->SetParameterRange("thickness>0.");     80   ThickPrm->SetParameterRange("thickness>0.");
 81   fAbsorCmd->SetParameter(ThickPrm);               81   fAbsorCmd->SetParameter(ThickPrm);
 82   //                                               82   //
 83   G4UIparameter* unitPrm = new G4UIparameter(" <<  83   G4UIparameter* unitPrm = new G4UIparameter("unit",'s',false);
 84   unitPrm->SetGuidance("unit of thickness");       84   unitPrm->SetGuidance("unit of thickness");
 85   G4String unitList = G4UIcommand::UnitsList(G     85   G4String unitList = G4UIcommand::UnitsList(G4UIcommand::CategoryOf("mm"));
 86   unitPrm->SetParameterCandidates(unitList);       86   unitPrm->SetParameterCandidates(unitList);
 87   fAbsorCmd->SetParameter(unitPrm);                87   fAbsorCmd->SetParameter(unitPrm);
 88   //                                               88   //
 89   fAbsorCmd->AvailableForStates(G4State_PreIni <<  89   fAbsorCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
 90   fAbsorCmd->SetToBeBroadcasted(false);        <<  90   
 91                                                <<  91   fNdivCmd = new G4UIcommand("/testem/det/nDivAbsor",this);
 92   fNdivCmd = new G4UIcommand("/testem/det/nDiv << 
 93   fNdivCmd->SetGuidance("Divide the absor nb :     92   fNdivCmd->SetGuidance("Divide the absor nb : number of divisions");
 94   fNdivCmd->SetGuidance("  absor number : from     93   fNdivCmd->SetGuidance("  absor number : from 1 to NbOfAbsor");
 95   fNdivCmd->SetGuidance("  number of divisions     94   fNdivCmd->SetGuidance("  number of divisions >= 0");
 96   //                                               95   //
 97   G4UIparameter* AbsNbPar = new G4UIparameter( <<  96   G4UIparameter* AbsNbPar = new G4UIparameter("AbsorNb",'i',false);
 98   AbsNbPar->SetGuidance("absor number : from 1     97   AbsNbPar->SetGuidance("absor number : from 1 to NbOfAbsor");
 99   AbsNbPar->SetParameterRange("AbsorNb>0");        98   AbsNbPar->SetParameterRange("AbsorNb>0");
100   fNdivCmd->SetParameter(AbsNbPar);                99   fNdivCmd->SetParameter(AbsNbPar);
101   //                                              100   //
102   G4UIparameter* NdivPrm = new G4UIparameter(" << 101   G4UIparameter* NdivPrm = new G4UIparameter("NdivNb",'i',false);
103   NdivPrm->SetGuidance("nb of divisions > 0");    102   NdivPrm->SetGuidance("nb of divisions > 0");
104   NdivPrm->SetParameterRange("NdivNb>0");         103   NdivPrm->SetParameterRange("NdivNb>0");
105   fNdivCmd->SetParameter(NdivPrm);                104   fNdivCmd->SetParameter(NdivPrm);
106   //                                              105   //
107   fNdivCmd->AvailableForStates(G4State_PreInit << 106   fNdivCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
108   fNdivCmd->SetToBeBroadcasted(false);         << 107       
109                                                << 108   fSizeYZCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setSizeYZ",this);
110   fSizeYZCmd = new G4UIcmdWithADoubleAndUnit(" << 
111   fSizeYZCmd->SetGuidance("Set sizeYZ of the a    109   fSizeYZCmd->SetGuidance("Set sizeYZ of the absorber");
112   fSizeYZCmd->SetParameterName("SizeYZ", false << 110   fSizeYZCmd->SetParameterName("SizeYZ",false);
113   fSizeYZCmd->SetRange("SizeYZ>0.");              111   fSizeYZCmd->SetRange("SizeYZ>0.");
114   fSizeYZCmd->SetUnitCategory("Length");          112   fSizeYZCmd->SetUnitCategory("Length");
115   fSizeYZCmd->AvailableForStates(G4State_PreIn << 113   fSizeYZCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
116   fSizeYZCmd->SetToBeBroadcasted(false);       << 114             
                                                   >> 115   fMagFieldCmd = new G4UIcmdWithADoubleAndUnit("/testem/det/setField",this);  
                                                   >> 116   fMagFieldCmd->SetGuidance("Define magnetic field.");
                                                   >> 117   fMagFieldCmd->SetGuidance("Magnetic field will be in Z direction.");
                                                   >> 118   fMagFieldCmd->SetParameterName("Bz",false);
                                                   >> 119   fMagFieldCmd->SetUnitCategory("Magnetic flux density");
                                                   >> 120   fMagFieldCmd->AvailableForStates(G4State_PreInit,G4State_Idle);
                                                   >> 121       
                                                   >> 122   fUpdateCmd = new G4UIcmdWithoutParameter("/testem/det/update",this);
                                                   >> 123   fUpdateCmd->SetGuidance("Update calorimeter geometry.");
                                                   >> 124   fUpdateCmd->SetGuidance("This command MUST be applied before \"beamOn\" ");
                                                   >> 125   fUpdateCmd->SetGuidance("if you changed geometrical value(s).");
                                                   >> 126   fUpdateCmd->AvailableForStates(G4State_Idle);
117 }                                                 127 }
118                                                   128 
119 //....oooOO0OOooo........oooOO0OOooo........oo    129 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
120                                                   130 
121 DetectorMessenger::~DetectorMessenger()           131 DetectorMessenger::~DetectorMessenger()
122 {                                                 132 {
123   delete fNbAbsorCmd;                             133   delete fNbAbsorCmd;
124   delete fAbsorCmd;                               134   delete fAbsorCmd;
125   delete fNdivCmd;                             << 135   delete fNdivCmd;    
126   delete fSizeYZCmd;                           << 136   delete fSizeYZCmd;       
127   delete fDetDir;                              << 137   delete fMagFieldCmd;
                                                   >> 138   delete fUpdateCmd;
                                                   >> 139   delete fDetDir;  
128   delete fTestemDir;                              140   delete fTestemDir;
129 }                                                 141 }
130                                                   142 
131 //....oooOO0OOooo........oooOO0OOooo........oo    143 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
132                                                   144 
133 void DetectorMessenger::SetNewValue(G4UIcomman << 145 void DetectorMessenger::SetNewValue(G4UIcommand* command,G4String newValue)
134 {                                              << 146 {    
135   if (command == fNbAbsorCmd) {                << 147   if( command == fNbAbsorCmd )
136     fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNe << 148    { fDetector->SetNbOfAbsor(fNbAbsorCmd->GetNewIntValue(newValue));}
137   }                                            << 149    
138                                                << 150   if (command == fAbsorCmd)
139   if (command == fAbsorCmd) {                  << 151    {
140     G4int num;                                 << 152      G4int num; G4double tick;
141     G4double tick;                             << 153      G4String unt, mat;
142     G4String unt, mat;                         << 154      std::istringstream is(newValue);
143     std::istringstream is(newValue);           << 155      is >> num >> mat >> tick >> unt;
144     is >> num >> mat >> tick >> unt;           << 156      G4String material=mat;
145     G4String material = mat;                   << 157      tick *= G4UIcommand::ValueOf(unt);
146     tick *= G4UIcommand::ValueOf(unt);         << 158      fDetector->SetAbsorMaterial (num,material);
147     fDetector->SetAbsorMaterial(num, material) << 159      fDetector->SetAbsorThickness(num,tick);
148     fDetector->SetAbsorThickness(num, tick);   << 160    }
149   }                                            << 161    
150                                                << 162   if (command == fNdivCmd)
151   if (command == fNdivCmd) {                   << 163    {
152     G4int num, ndiv;                           << 164      G4int num, ndiv;
153     std::istringstream is(newValue);           << 165      std::istringstream is(newValue);
154     is >> num >> ndiv;                         << 166      is >> num >> ndiv;
155     fDetector->SetNbOfDivisions(num, ndiv);    << 167      fDetector->SetNbOfDivisions (num,ndiv);
156   }                                            << 168    }
157                                                << 169       
158   if (command == fSizeYZCmd) {                 << 170   if( command == fSizeYZCmd )
159     fDetector->SetAbsorSizeYZ(fSizeYZCmd->GetN << 171    { fDetector->SetAbsorSizeYZ(fSizeYZCmd->GetNewDoubleValue(newValue));}   
160   }                                            << 172                
                                                   >> 173   if( command == fMagFieldCmd )
                                                   >> 174    { fDetector->SetMagField(fMagFieldCmd->GetNewDoubleValue(newValue));}
                                                   >> 175               
                                                   >> 176   if( command == fUpdateCmd )
                                                   >> 177    { fDetector->UpdateGeometry();}
161 }                                                 178 }
162                                                   179 
163 //....oooOO0OOooo........oooOO0OOooo........oo    180 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
164                                                   181