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 ////////////////////////////////////////////// 26 //////////////////////////////////////////////////////////////////////////////// 27 // << 27 // // 28 // GEANT4 Class source file << 28 // File: G4RadioactivationMessenger.cc // 29 // << 29 // Author: D.H. Wright (SLAC) // 30 // G4RadioactivationMessenger << 30 // Date: 29 August 2017 // 31 // << 31 // Description: messenger class for biased version of G4RadioactiveDecay. // 32 // Author: D.H. Wright (SLAC) << 32 // Based on the code of F. Lei and P.R. Truscott. // 33 // Date: 29 August 2017 << 33 // // 34 // << 35 // Based on the code of F. Lei and P.R. Trusc << 36 // << 37 ////////////////////////////////////////////// 34 //////////////////////////////////////////////////////////////////////////////// 38 35 39 #include "G4RadioactivationMessenger.hh" 36 #include "G4RadioactivationMessenger.hh" 40 #include "G4NuclearLevelData.hh" 37 #include "G4NuclearLevelData.hh" 41 #include "G4RadioactiveDecay.hh" << 42 << 43 #include <sstream> 38 #include <sstream> 44 #include "G4HadronicException.hh" 39 #include "G4HadronicException.hh" 45 40 46 41 47 G4RadioactivationMessenger::G4RadioactivationM << 42 G4RadioactivationMessenger::G4RadioactivationMessenger(G4Radioactivation* theRadioactivationContainer1) 48 :theRadDecay(ptr) << 43 :theRadioactivationContainer(theRadioactivationContainer1) 49 { 44 { 50 rdmDirectory = new G4UIdirectory("/process/h 45 rdmDirectory = new G4UIdirectory("/process/had/rdm/"); 51 rdmDirectory->SetGuidance("Controls the bias 46 rdmDirectory->SetGuidance("Controls the biased version of radioactive decay"); 52 47 53 // Command to turn on/off variance reduction 48 // Command to turn on/off variance reduction options 54 analoguemcCmd = new G4UIcmdWithABool("/proce 49 analoguemcCmd = new G4UIcmdWithABool("/process/had/rdm/analogueMC",this); 55 analoguemcCmd->SetGuidance("false: variance 50 analoguemcCmd->SetGuidance("false: variance reduction method; true: analogue method"); 56 analoguemcCmd->SetParameterName("AnalogueMC" 51 analoguemcCmd->SetParameterName("AnalogueMC",true); 57 analoguemcCmd->SetDefaultValue(true); 52 analoguemcCmd->SetDefaultValue(true); 58 53 59 // Command to use branching ratio biasing or 54 // Command to use branching ratio biasing or not 60 brbiasCmd = new G4UIcmdWithABool("/process/h 55 brbiasCmd = new G4UIcmdWithABool("/process/had/rdm/BRbias",this); 61 brbiasCmd->SetGuidance("false: no biasing; t 56 brbiasCmd->SetGuidance("false: no biasing; true: all branches are treated as equal"); 62 brbiasCmd->SetParameterName("BRBias",true); 57 brbiasCmd->SetParameterName("BRBias",true); 63 brbiasCmd->SetDefaultValue(true); 58 brbiasCmd->SetDefaultValue(true); 64 59 65 // Command to set the half-life thresold for 60 // Command to set the half-life thresold for isomer production 66 hlthCmd = new G4UIcmdWithADoubleAndUnit("/pr 61 hlthCmd = new G4UIcmdWithADoubleAndUnit("/process/had/rdm/hlThreshold",this); 67 hlthCmd->SetGuidance("Set the h-l threshold 62 hlthCmd->SetGuidance("Set the h-l threshold for isomer production"); 68 hlthCmd->SetParameterName("hlThreshold",fals 63 hlthCmd->SetParameterName("hlThreshold",false); 69 hlthCmd->SetUnitCategory("Time"); 64 hlthCmd->SetUnitCategory("Time"); 70 65 71 // Command to define the incident particle s 66 // Command to define the incident particle source time profile 72 sourcetimeprofileCmd = new G4UIcmdWithAStrin 67 sourcetimeprofileCmd = new G4UIcmdWithAString("/process/had/rdm/sourceTimeProfile",this); 73 sourcetimeprofileCmd->SetGuidance 68 sourcetimeprofileCmd->SetGuidance 74 ("Supply the name of the ascii file contai 69 ("Supply the name of the ascii file containing the source particle time profile"); 75 sourcetimeprofileCmd->SetParameterName("STim 70 sourcetimeprofileCmd->SetParameterName("STimeProfile",true); 76 sourcetimeprofileCmd->SetDefaultValue("sourc 71 sourcetimeprofileCmd->SetDefaultValue("source.data"); 77 72 78 // Command to define the incident particle s 73 // Command to define the incident particle source time profile 79 decaybiasprofileCmd = new G4UIcmdWithAString 74 decaybiasprofileCmd = new G4UIcmdWithAString("/process/had/rdm/decayBiasProfile",this); 80 decaybiasprofileCmd->SetGuidance 75 decaybiasprofileCmd->SetGuidance 81 ("Supply the name of the ascii file contai 76 ("Supply the name of the ascii file containing the decay bias time profile"); 82 decaybiasprofileCmd->SetParameterName("DBias 77 decaybiasprofileCmd->SetParameterName("DBiasProfile",true); 83 decaybiasprofileCmd->SetDefaultValue("bias.d 78 decaybiasprofileCmd->SetDefaultValue("bias.data"); 84 79 85 // Command to set nuclei splitting parameter 80 // Command to set nuclei splitting parameter 86 splitnucleiCmd = new G4UIcmdWithAnInteger("/ 81 splitnucleiCmd = new G4UIcmdWithAnInteger("/process/had/rdm/splitNuclei",this); 87 splitnucleiCmd->SetGuidance("Set number of s 82 splitnucleiCmd->SetGuidance("Set number of splitting for the isotopes."); 88 splitnucleiCmd->SetParameterName("NSplit",tr 83 splitnucleiCmd->SetParameterName("NSplit",true); 89 splitnucleiCmd->SetDefaultValue(1); 84 splitnucleiCmd->SetDefaultValue(1); 90 splitnucleiCmd->SetRange("NSplit>=1"); 85 splitnucleiCmd->SetRange("NSplit>=1"); 91 } 86 } 92 87 93 88 94 G4RadioactivationMessenger::~G4Radioactivation 89 G4RadioactivationMessenger::~G4RadioactivationMessenger() 95 { 90 { 96 delete rdmDirectory; 91 delete rdmDirectory; 97 delete analoguemcCmd; 92 delete analoguemcCmd; 98 delete sourcetimeprofileCmd; 93 delete sourcetimeprofileCmd; 99 delete decaybiasprofileCmd; 94 delete decaybiasprofileCmd; 100 delete brbiasCmd; 95 delete brbiasCmd; 101 delete splitnucleiCmd; 96 delete splitnucleiCmd; 102 delete hlthCmd; 97 delete hlthCmd; 103 } 98 } 104 99 105 100 106 void G4RadioactivationMessenger::SetNewValue(G 101 void G4RadioactivationMessenger::SetNewValue(G4UIcommand* command, G4String newValues) 107 { 102 { 108 if ( command == analoguemcCmd ) { theRadDeca << 103 if ( command == analoguemcCmd ) { theRadioactivationContainer-> 109 SetAnalogueMonteCarlo( analoguemcCmd->GetN 104 SetAnalogueMonteCarlo( analoguemcCmd->GetNewBoolValue( newValues ) ); 110 } else if ( command == brbiasCmd ) { theRadD << 105 } else if ( command == brbiasCmd ) { theRadioactivationContainer-> 111 SetBRBias( brbiasCmd->GetNewBoolValue( new 106 SetBRBias( brbiasCmd->GetNewBoolValue( newValues ) ); 112 } else if ( command == sourcetimeprofileCmd << 107 } else if ( command == sourcetimeprofileCmd ) { theRadioactivationContainer-> 113 SetSourceTimeProfile( newValues ); 108 SetSourceTimeProfile( newValues ); 114 } else if ( command == decaybiasprofileCmd ) << 109 } else if ( command == decaybiasprofileCmd ) { theRadioactivationContainer-> 115 SetDecayBias( newValues ); 110 SetDecayBias( newValues ); 116 } else if ( command == splitnucleiCmd ) { th << 111 } else if ( command == splitnucleiCmd ) { theRadioactivationContainer-> 117 SetSplitNuclei( splitnucleiCmd->GetNewIntV 112 SetSplitNuclei( splitnucleiCmd->GetNewIntValue( newValues ) ); 118 } else if ( command == hlthCmd ) { theRadDec << 113 } else if ( command == hlthCmd ) { theRadioactivationContainer-> 119 SetHLThreshold( hlthCmd->GetNewDoubleValue 114 SetHLThreshold( hlthCmd->GetNewDoubleValue( newValues ) ); 120 } 115 } 121 } 116 } 122 117 123 118