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 // Geant4 Class file << 26 // $Id: G4PolarizationMessenger.cc 81375 2014-05-27 13:08:40Z gcosmo $ >> 27 // >> 28 // >> 29 // GEANT4 Class file 27 // 30 // 28 // File name: G4PolarizationManager 31 // File name: G4PolarizationManager 29 // 32 // 30 // Author: Andreas Schaelicke 33 // Author: Andreas Schaelicke 31 // 34 // >> 35 // Creation date: 01.05.2005 >> 36 // >> 37 // Modifications: >> 38 // 32 // Class Description: 39 // Class Description: 33 // Provides access to general polarization i << 40 // 34 // polarization for logical volumes through << 41 // Provides access to general polarization information and to >> 42 // polarization for logical volumes through macro files. 35 43 36 #include "G4PolarizationMessenger.hh" 44 #include "G4PolarizationMessenger.hh" 37 << 38 #include "G4ios.hh" << 39 #include "G4PolarizationManager.hh" 45 #include "G4PolarizationManager.hh" 40 #include "G4PolarizationHelper.hh" 46 #include "G4PolarizationHelper.hh" >> 47 >> 48 #include "G4UIdirectory.hh" 41 #include "G4Tokenizer.hh" 49 #include "G4Tokenizer.hh" 42 #include "G4UIcmdWithABool.hh" << 43 #include "G4UIcmdWithAnInteger.hh" << 44 #include "G4UIcmdWithAString.hh" << 45 #include "G4UIcmdWithoutParameter.hh" 50 #include "G4UIcmdWithoutParameter.hh" 46 #include "G4UIdirectory.hh" << 51 #include "G4UIcmdWithAString.hh" >> 52 #include "G4UIcmdWithAnInteger.hh" >> 53 #include "G4UIcmdWithABool.hh" >> 54 #include "G4ios.hh" >> 55 47 56 48 G4PolarizationMessenger::G4PolarizationMesseng << 57 G4PolarizationMessenger::G4PolarizationMessenger(G4PolarizationManager * polMgr) 49 : polarizationManager(polMgr) 58 : polarizationManager(polMgr) 50 { 59 { 51 polarizationDirectory = new G4UIdirectory("/ 60 polarizationDirectory = new G4UIdirectory("/polarization/"); 52 polarizationDirectory->SetGuidance("polariza 61 polarizationDirectory->SetGuidance("polarization control commands."); 53 62 54 managerDirectory = new G4UIdirectory("/polar 63 managerDirectory = new G4UIdirectory("/polarization/manager/"); 55 managerDirectory->SetGuidance("general polar 64 managerDirectory->SetGuidance("general polarization information."); 56 65 57 verboseCmd = new G4UIcmdWithAnInteger("/pola << 66 verboseCmd = new G4UIcmdWithAnInteger("/polarization/manager/verbose",this); 58 verboseCmd->SetGuidance("Set the Verbose lev 67 verboseCmd->SetGuidance("Set the Verbose level of G4PolarizationManager."); 59 verboseCmd->SetGuidance(" 0 : Silent (defaul 68 verboseCmd->SetGuidance(" 0 : Silent (default)"); 60 verboseCmd->SetGuidance(" 1 : Verbose"); 69 verboseCmd->SetGuidance(" 1 : Verbose"); 61 verboseCmd->SetParameterName("level", true); << 70 verboseCmd->SetParameterName("level",true); 62 verboseCmd->SetDefaultValue(0); 71 verboseCmd->SetDefaultValue(0); 63 verboseCmd->SetRange("level >=0 && level <=1 72 verboseCmd->SetRange("level >=0 && level <=1"); 64 73 65 optActivateCmd = new G4UIcmdWithABool("/pola << 74 optActivateCmd = new G4UIcmdWithABool("/polarization/manager/activate",this); 66 optActivateCmd->SetGuidance("activate/deacti 75 optActivateCmd->SetGuidance("activate/deactivate polarization treatment"); 67 optActivateCmd->SetParameterName("flag", tru << 76 optActivateCmd->SetParameterName("flag",true); 68 optActivateCmd->SetDefaultValue(true); 77 optActivateCmd->SetDefaultValue(true); 69 78 70 volumeDirectory = new G4UIdirectory("/polari 79 volumeDirectory = new G4UIdirectory("/polarization/volume/"); 71 volumeDirectory->SetGuidance( << 80 volumeDirectory->SetGuidance("Status control commands of registered polarized logical volumes."); 72 "Status control commands of registered pol << 73 81 74 printVolumeListCmd = << 82 printVolumeListCmd = new G4UIcmdWithoutParameter("/polarization/volume/list",this); 75 new G4UIcmdWithoutParameter("/polarization << 83 printVolumeListCmd->SetGuidance("print list of registered polarized logical volumes"); 76 printVolumeListCmd->SetGuidance( << 84 printVolumeListCmd->AvailableForStates(G4State_PreInit,G4State_Idle,G4State_GeomClosed); 77 "print list of registered polarized logica << 85 78 printVolumeListCmd->AvailableForStates(G4Sta << 86 setPolarizationCmd = new G4UIcommand("/polarization/volume/set",this); 79 G4Sta << 87 setPolarizationCmd->SetGuidance("set or change polarization of a logical volume"); 80 << 88 // setPolarizationCmd->SetParameterName("polarization",true); 81 setPolarizationCmd = new G4UIcommand("/polar << 89 // setPolarizationCmd->SetDefaultValue("worldVolume 0. 0. 0."); 82 setPolarizationCmd->SetGuidance( << 90 setPolarizationCmd->AvailableForStates(G4State_PreInit,G4State_Idle,G4State_GeomClosed); 83 "set or change polarization of a logical v << 84 setPolarizationCmd->AvailableForStates(G4Sta << 85 G4Sta << 86 91 87 G4UIparameter* param; 92 G4UIparameter* param; 88 param = new G4UIparameter("logicalVolumeName << 93 param = new G4UIparameter("logicalVolumeName",'s',false); 89 param->SetDefaultValue("worldVolume"); 94 param->SetDefaultValue("worldVolume"); 90 setPolarizationCmd->SetParameter(param); 95 setPolarizationCmd->SetParameter(param); 91 param = new G4UIparameter("px", 'd', true); << 96 param = new G4UIparameter("px",'d',true); 92 param->SetDefaultValue("0.0"); 97 param->SetDefaultValue("0.0"); 93 setPolarizationCmd->SetParameter(param); 98 setPolarizationCmd->SetParameter(param); 94 param = new G4UIparameter("py", 'd', true); << 99 param = new G4UIparameter("py",'d',true); 95 param->SetDefaultValue("0.0"); 100 param->SetDefaultValue("0.0"); 96 setPolarizationCmd->SetParameter(param); 101 setPolarizationCmd->SetParameter(param); 97 param = new G4UIparameter("pz", 'd', true); << 102 param = new G4UIparameter("pz",'d',true); 98 param->SetDefaultValue("0.0"); 103 param->SetDefaultValue("0.0"); 99 setPolarizationCmd->SetParameter(param); 104 setPolarizationCmd->SetParameter(param); 100 105 101 testDirectory = new G4UIdirectory("/polariza 106 testDirectory = new G4UIdirectory("/polarization/test/"); 102 testDirectory->SetGuidance("provides access 107 testDirectory->SetGuidance("provides access to some internal test routines."); 103 108 104 testPolarizationTransformationCmd = new G4UI << 109 testPolarizationTransformationCmd = new G4UIcmdWithoutParameter("/polarization/test/polarizationTransformation",this); 105 "/polarization/test/polarizationTransforma << 110 testPolarizationTransformationCmd->SetGuidance("checks definition of particle reference frame and corresponding translation routines"); 106 testPolarizationTransformationCmd->SetGuidan << 111 testPolarizationTransformationCmd->AvailableForStates(G4State_PreInit,G4State_Idle,G4State_GeomClosed); 107 "checks definition of particle reference f << 112 108 "translation routines"); << 113 testInteractionFrameCmd = new G4UIcmdWithoutParameter("/polarization/test/interactionFrame",this); 109 testPolarizationTransformationCmd->Available << 114 testInteractionFrameCmd->SetGuidance("checks definition of interaction frame"); 110 G4State_PreInit, G4State_Idle, G4State_Geo << 115 testInteractionFrameCmd->AvailableForStates(G4State_PreInit,G4State_Idle,G4State_GeomClosed); 111 << 112 testInteractionFrameCmd = << 113 new G4UIcmdWithoutParameter("/polarization << 114 testInteractionFrameCmd->SetGuidance( << 115 "checks definition of interaction frame"); << 116 testInteractionFrameCmd->AvailableForStates( << 117 << 118 } 116 } 119 117 120 G4PolarizationMessenger::~G4PolarizationMessen 118 G4PolarizationMessenger::~G4PolarizationMessenger() 121 { 119 { 122 delete verboseCmd; 120 delete verboseCmd; 123 delete testInteractionFrameCmd; 121 delete testInteractionFrameCmd; 124 delete testPolarizationTransformationCmd; 122 delete testPolarizationTransformationCmd; 125 delete testDirectory; 123 delete testDirectory; 126 delete setPolarizationCmd; 124 delete setPolarizationCmd; 127 delete printVolumeListCmd; 125 delete printVolumeListCmd; 128 delete volumeDirectory; 126 delete volumeDirectory; 129 delete optActivateCmd; 127 delete optActivateCmd; 130 delete managerDirectory; 128 delete managerDirectory; 131 delete polarizationDirectory; 129 delete polarizationDirectory; 132 } 130 } 133 131 134 void G4PolarizationMessenger::SetNewValue(G4UI << 132 void G4PolarizationMessenger::SetNewValue(G4UIcommand * command,G4String newValue) 135 G4St << 136 { 133 { 137 if(command == verboseCmd) << 134 if( command==verboseCmd ) { 138 { << 135 polarizationManager->SetVerbose(verboseCmd->GetNewIntValue(newValue)); 139 polarizationManager->SetVerbose(verboseCmd << 140 } << 141 else if(command == optActivateCmd) << 142 { << 143 polarizationManager->SetActivated( << 144 optActivateCmd->GetNewBoolValue(newValue << 145 } 136 } 146 else if(command == printVolumeListCmd) << 137 else if ( command==optActivateCmd ) { 147 { << 138 polarizationManager->SetActivated(optActivateCmd->GetNewBoolValue(newValue)); >> 139 } >> 140 else if ( command==printVolumeListCmd ) { 148 polarizationManager->ListVolumes(); 141 polarizationManager->ListVolumes(); 149 } 142 } 150 else if(command == setPolarizationCmd) << 143 else if ( command==setPolarizationCmd ) { 151 { << 144 G4Tokenizer next( newValue ); 152 G4Tokenizer next(newValue); << 145 G4String volumeName=next(); 153 G4String volumeName = next(); << 146 G4double px=0.,py=0.,pz=0.; 154 G4double px = 0., py = 0., pz = 0.; << 147 G4String dvalue=next(); 155 G4String dvalue = next(); << 148 if (!dvalue.isNull()) { 156 if(!dvalue.empty()) << 149 px=StoD(dvalue); 157 { << 150 dvalue=next(); 158 px = StoD(dvalue); << 151 if (!dvalue.isNull()) { 159 dvalue = next(); << 152 py=StoD(dvalue); 160 if(!dvalue.empty()) << 153 dvalue=next(); 161 { << 154 if (!dvalue.isNull()) pz=StoD(dvalue); 162 py = StoD(dvalue); << 163 dvalue = next(); << 164 if(!dvalue.empty()) << 165 pz = StoD(dvalue); << 166 } 155 } 167 } 156 } 168 G4ThreeVector pol(px, py, pz); << 157 G4ThreeVector pol(px,py,pz); 169 polarizationManager->SetVolumePolarization << 158 polarizationManager->SetVolumePolarization(volumeName,pol); 170 } 159 } 171 else if(command == testPolarizationTransform << 160 else if ( command==testPolarizationTransformationCmd ) { 172 { << 173 G4PolarizationHelper::TestPolarizationTran 161 G4PolarizationHelper::TestPolarizationTransformations(); 174 } 162 } 175 else if(command == testInteractionFrameCmd) << 163 else if (command==testInteractionFrameCmd ) { 176 { << 177 G4PolarizationHelper::TestInteractionFrame 164 G4PolarizationHelper::TestInteractionFrame(); 178 } 165 } 179 } 166 } 180 167 181 G4String G4PolarizationMessenger::GetCurrentVa << 168 G4String G4PolarizationMessenger::GetCurrentValue(G4UIcommand * command) 182 { 169 { 183 G4String cv; 170 G4String cv; 184 if(command == verboseCmd) << 171 if( command==verboseCmd ) 185 { << 172 { cv = verboseCmd->ConvertToString(polarizationManager->GetVerbose()); } 186 cv = verboseCmd->ConvertToString(polarizat << 187 } << 188 173 189 return cv; 174 return cv; 190 } 175 } 191 176