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 exoticphysics/monopole/src/G4Monopol << 26 // $Id: G4MonopolePhysicsMessenger.cc,v 1.1 2009/07/15 10:20:07 vnivanch Exp $ 27 /// \brief Implementation of the G4MonopolePhy << 27 // GEANT4 tag $Name: geant4-09-03 $ 28 // << 29 // << 30 // 12.07.10 S.Burdin (changed the magnetic a << 31 // from integer to double) << 32 // 28 // 33 //....oooOO0OOooo........oooOO0OOooo........oo 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 34 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 35 31 36 #include "G4MonopolePhysicsMessenger.hh" 32 #include "G4MonopolePhysicsMessenger.hh" 37 33 38 #include "G4MonopolePhysics.hh" 34 #include "G4MonopolePhysics.hh" 39 #include "G4UIcmdWithADouble.hh" << 40 #include "G4UIcmdWithADoubleAndUnit.hh" << 41 #include "G4UIcommand.hh" << 42 #include "G4UIdirectory.hh" 35 #include "G4UIdirectory.hh" >> 36 #include "G4UIcommand.hh" >> 37 #include "G4UIcmdWithAnInteger.hh" >> 38 #include "G4UIcmdWithADoubleAndUnit.hh" 43 39 44 #include <sstream> 40 #include <sstream> 45 41 46 //....oooOO0OOooo........oooOO0OOooo........oo 42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 47 43 48 G4MonopolePhysicsMessenger::G4MonopolePhysicsM 44 G4MonopolePhysicsMessenger::G4MonopolePhysicsMessenger(G4MonopolePhysics* p) 49 : G4UImessenger(), fPhys(p) << 45 : phys(p) 50 { 46 { 51 fPhysicsDir = new G4UIdirectory("/monopole/" << 47 mPhysicsDir = new G4UIdirectory("/monopole/"); 52 fPhysicsDir->SetGuidance("histograms control << 48 mPhysicsDir->SetGuidance("histograms control"); 53 << 49 54 fPhysicsCmd = new G4UIcommand("/monopole/set << 50 mPhysicsCmd = new G4UIcommand("/monopole/setup",this); 55 fPhysicsCmd->SetGuidance("Setup monopole"); << 51 mPhysicsCmd->SetGuidance("Setup monopole"); 56 // 52 // 57 G4UIparameter* qmag = new G4UIparameter("qma << 53 G4UIparameter* qmag = new G4UIparameter("qmag",'i',false); 58 qmag->SetGuidance("Magnetic charge"); 54 qmag->SetGuidance("Magnetic charge"); 59 qmag->SetDefaultValue("1"); 55 qmag->SetDefaultValue("1"); 60 fPhysicsCmd->SetParameter(qmag); << 56 mPhysicsCmd->SetParameter(qmag); 61 57 62 G4UIparameter* q = new G4UIparameter("qelec" << 58 G4UIparameter* q = new G4UIparameter("qelec",'i',false); 63 q->SetGuidance("Electric charge charge"); 59 q->SetGuidance("Electric charge charge"); 64 q->SetDefaultValue("0"); 60 q->SetDefaultValue("0"); 65 fPhysicsCmd->SetParameter(q); << 61 mPhysicsCmd->SetParameter(q); 66 // << 62 // 67 G4UIparameter* mass = new G4UIparameter("mas << 63 G4UIparameter* mass = new G4UIparameter("mass",'d',false); 68 mass->SetGuidance("mass"); 64 mass->SetGuidance("mass"); 69 mass->SetParameterRange("mass>0."); 65 mass->SetParameterRange("mass>0."); 70 qmag->SetDefaultValue("100"); 66 qmag->SetDefaultValue("100"); 71 fPhysicsCmd->SetParameter(mass); << 67 mPhysicsCmd->SetParameter(mass); 72 // << 68 // 73 G4UIparameter* unit = new G4UIparameter("uni << 69 G4UIparameter* unit = new G4UIparameter("unit",'s',false); 74 fPhysicsCmd->SetParameter(unit); << 70 mPhysicsCmd->SetParameter(unit); 75 qmag->SetDefaultValue("GeV"); 71 qmag->SetDefaultValue("GeV"); 76 fPhysicsCmd->AvailableForStates(G4State_PreI << 72 mPhysicsCmd->AvailableForStates(G4State_PreInit); >> 73 >> 74 mCmd = new G4UIcmdWithAnInteger("/monopole/magCharge",this); >> 75 mCmd->SetGuidance("Set monopole magnetic charge number"); >> 76 mCmd->SetParameterName("Qmag",false); >> 77 mCmd->AvailableForStates(G4State_PreInit); >> 78 >> 79 zCmd = new G4UIcmdWithAnInteger("/monopole/elCharge",this); >> 80 zCmd->SetGuidance("Set monopole electric charge number"); >> 81 zCmd->SetParameterName("Qel",false); >> 82 zCmd->AvailableForStates(G4State_PreInit); >> 83 >> 84 massCmd = new G4UIcmdWithADoubleAndUnit("/monopole/mass",this); >> 85 massCmd->SetGuidance("Set monopole mass"); >> 86 massCmd->SetParameterName("Mass",false); >> 87 massCmd->SetRange("Mass>0."); >> 88 massCmd->SetUnitCategory("Energy"); >> 89 massCmd->AvailableForStates(G4State_PreInit); >> 90 77 } 91 } 78 92 79 //....oooOO0OOooo........oooOO0OOooo........oo 93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 80 94 81 G4MonopolePhysicsMessenger::~G4MonopolePhysics 95 G4MonopolePhysicsMessenger::~G4MonopolePhysicsMessenger() 82 { 96 { 83 delete fPhysicsCmd; << 97 delete mPhysicsCmd; 84 delete fPhysicsDir; << 98 delete mCmd; >> 99 delete zCmd; >> 100 delete massCmd; >> 101 delete mPhysicsDir; 85 } 102 } 86 103 87 //....oooOO0OOooo........oooOO0OOooo........oo 104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 88 105 89 void G4MonopolePhysicsMessenger::SetNewValue(G << 106 void G4MonopolePhysicsMessenger::SetNewValue(G4UIcommand* command, 90 { << 107 G4String newValue) 91 if (command == fPhysicsCmd) { << 108 { 92 G4double q, m; << 109 if (command == mPhysicsCmd) 93 G4double mass; << 110 { G4int q, m; G4double mass; 94 G4String unts; << 111 G4String unts; 95 std::istringstream is(newValue); << 112 std::istringstream is(newValue); 96 is >> m >> q >> mass >> unts; << 113 is >> m >> q >> mass >> unts; 97 G4String unit = unts; << 114 G4String unit = unts; 98 G4double vUnit = G4UIcommand::ValueOf(unit << 115 G4double vUnit = G4UIcommand::ValueOf(unit); 99 fPhys->SetMagneticCharge(m); << 116 phys->SetMagneticCharge(m); 100 fPhys->SetElectricCharge(q); << 117 phys->SetElectricCharge(q); 101 fPhys->SetMonopoleMass(mass * vUnit); << 118 phys->SetMonopoleMass(mass*vUnit); 102 } << 119 } >> 120 if (command == mCmd) {phys->SetMagneticCharge(mCmd->GetNewIntValue(newValue));} >> 121 if (command == zCmd) {phys->SetElectricCharge(zCmd->GetNewIntValue(newValue));} >> 122 if (command == massCmd) {phys->SetMonopoleMass(massCmd->GetNewDoubleValue(newValue));} 103 } 123 } 104 124 105 //....oooOO0OOooo........oooOO0OOooo........oo 125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 106 126