Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // G4QSSMessenger implementation 27 // 28 // Author: Leandro Gomez Vidal (Univ. Buenos Aires) - October 2021 29 // -------------------------------------------------------------------- 30 31 #include "G4QSSMessenger.hh" 32 33 G4QSSMessenger::G4QSSMessenger() 34 { 35 commandsShouldBeInMaster = true; 36 37 qssCmdDir = new G4UIdirectory("/QSS/", false); 38 qssCmdDir->SetGuidance("G4QSStepper configuration."); 39 40 dQMinCmd = new G4UIcmdWithADoubleAndUnit("/QSS/dQMin",this); 41 dQMinCmd->SetDefaultUnit("mm"); 42 dQMinCmd->SetParameterName("dQMinCmd",false); 43 dQMinCmd->SetUnitCategory("Length"); 44 45 dQRelCmd = new G4UIcmdWithADouble("/QSS/dQRel",this); 46 dQRelCmd->SetGuidance("Default is 1e-5"); 47 dQRelCmd->SetParameterName("dQRelCmd",false); 48 49 trialProposedStepModifierCmd = new G4UIcmdWithADouble("/QSS/trialProposedStepModifier",this); 50 trialProposedStepModifierCmd->SetGuidance("Default is 1"); 51 trialProposedStepModifierCmd->SetParameterName("trialProposedStepModifier", false); 52 53 stepperSelectorCmd = new G4UIcmdWithAString("/QSS/selectStepper",this); 54 stepperSelectorCmd->SetGuidance("Select stepper."); 55 stepperSelectorCmd->SetParameterName("choice", false); 56 stepperSelectorCmd->SetCandidates("TemplatedDoPri OldRK45 G4QSS2"); 57 58 } 59 60 G4QSSMessenger::~G4QSSMessenger() 61 { 62 delete qssCmdDir; 63 delete dQMinCmd; 64 delete dQRelCmd; 65 delete stepperSelectorCmd; 66 delete trialProposedStepModifierCmd; 67 //qssStats.print(); 68 } 69 70 G4QSSMessenger* G4QSSMessenger::instance() 71 { 72 static G4QSSMessenger theSingleMessengerInstance; 73 return &theSingleMessengerInstance; 74 } 75 76 void G4QSSMessenger::SetNewValue(G4UIcommand *command, G4String newValue) 77 { 78 if ( command == dQMinCmd ) { 79 dQMin = dQMinCmd->GetNewDoubleValue(newValue); 80 } 81 82 if ( command == dQRelCmd ) { 83 dQRel = dQRelCmd->GetNewDoubleValue(newValue); 84 } 85 86 if ( command == trialProposedStepModifierCmd ) { 87 trialProposedStepModifier = trialProposedStepModifierCmd->GetNewDoubleValue(newValue); 88 } 89 90 if(command == stepperSelectorCmd){ 91 this->selectStepper(newValue); 92 } 93 94 } 95 96 void G4QSSMessenger::selectStepper(const std::string &newValue) 97 { 98 const std::map<std::string, StepperSelection> stepperMapping = 99 {{"TemplatedDoPri", TemplatedDoPri}, {"OldRK45", OldRK45}, {"G4QSS2", G4QSS2}}; 100 _selectedStepper = stepperMapping.at(newValue); 101 G4cout << "G4QSSMessenger: Selecting stepper " << newValue << G4endl; 102 } 103 104 G4QSSMessenger::StepperSelection G4QSSMessenger::selectedStepper() 105 { 106 return _selectedStepper; 107 } 108