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 electromagnetic/TestEm8/src/PhysicsL << 27 /// \brief Implementation of the PhysicsListMe << 28 // 26 // >> 27 // $Id: PhysicsListMessenger.cc,v 1.2 2006/06/29 17:00:47 gunter Exp $ >> 28 // GEANT4 tag $Name: geant4-09-03-patch-02 $ 29 // 29 // 30 //-------------------------------------------- << 30 // 31 // << 32 // ClassName: PhysicsListMessenger << 33 // << 34 // Description: EM physics with a possibility << 35 // << 36 // Author: V.Ivanchenko 01.09.2010 << 37 // << 38 //-------------------------------------------- << 39 // << 40 31 41 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 42 //....oooOO0OOooo........oooOO0OOooo........oo 33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 43 34 44 #include "PhysicsListMessenger.hh" 35 #include "PhysicsListMessenger.hh" 45 36 46 #include "PhysicsList.hh" 37 #include "PhysicsList.hh" 47 #include "TestParameters.hh" << 48 << 49 #include "G4UIcmdWithADouble.hh" << 50 #include "G4UIcmdWithADoubleAndUnit.hh" 38 #include "G4UIcmdWithADoubleAndUnit.hh" 51 #include "G4UIcmdWithAString.hh" 39 #include "G4UIcmdWithAString.hh" 52 #include "G4UIcmdWithAnInteger.hh" << 53 #include "G4UIdirectory.hh" << 54 40 55 //....oooOO0OOooo........oooOO0OOooo........oo 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 56 42 57 PhysicsListMessenger::PhysicsListMessenger(Phy << 43 PhysicsListMessenger::PhysicsListMessenger(PhysicsList* pPhys) 58 { << 44 :pPhysicsList(pPhys) 59 fPhysDir = new G4UIdirectory("/testem/phys/" << 45 { 60 fPhysDir->SetGuidance("physics list commands << 46 gammaCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setGCut",this); 61 << 47 gammaCutCmd->SetGuidance("Set gamma cut."); 62 fECmd = new G4UIcmdWithADoubleAndUnit("/test << 48 gammaCutCmd->SetParameterName("Gcut",false); 63 fECmd->SetGuidance("Set max energy deposit") << 49 gammaCutCmd->SetUnitCategory("Length"); 64 fECmd->SetParameterName("Emax", false); << 50 gammaCutCmd->SetRange("Gcut>0.0"); 65 fECmd->SetUnitCategory("Energy"); << 51 gammaCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 66 fECmd->SetRange("Emax>0.0"); << 52 67 fECmd->AvailableForStates(G4State_PreInit, G << 53 electCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setECut",this); 68 << 54 electCutCmd->SetGuidance("Set electron cut."); 69 fEBCmd = new G4UIcmdWithAnInteger("/testem/p << 55 electCutCmd->SetParameterName("Ecut",false); 70 fEBCmd->SetGuidance("Set number of bins in e << 56 electCutCmd->SetUnitCategory("Length"); 71 fEBCmd->SetParameterName("Ebins", false); << 57 electCutCmd->SetRange("Ecut>0.0"); 72 fEBCmd->SetRange("Ebins>0"); << 58 electCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 73 fEBCmd->AvailableForStates(G4State_PreInit, << 59 74 << 60 protoCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setPCut",this); 75 fCBCmd = new G4UIcmdWithAnInteger("/testem/p << 61 protoCutCmd->SetGuidance("Set positron cut."); 76 fCBCmd->SetGuidance("Set number of bins of c << 62 protoCutCmd->SetParameterName("Pcut",false); 77 fCBCmd->SetParameterName("Cbins", false); << 63 protoCutCmd->SetUnitCategory("Length"); 78 fCBCmd->SetRange("Cbins>0"); << 64 protoCutCmd->SetRange("Pcut>0.0"); 79 fCBCmd->AvailableForStates(G4State_PreInit, << 65 protoCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 80 << 66 81 fCMCmd = new G4UIcmdWithAnInteger("/testem/p << 67 allCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setCuts",this); 82 fCMCmd->SetGuidance("Set max number of clust << 68 allCutCmd->SetGuidance("Set cut for all."); 83 fCMCmd->SetParameterName("Cmax", false); << 69 allCutCmd->SetParameterName("cut",false); 84 fCMCmd->SetRange("Cmax>0"); << 70 allCutCmd->SetUnitCategory("Length"); 85 fCMCmd->AvailableForStates(G4State_PreInit, << 71 allCutCmd->SetRange("cut>0.0"); 86 << 72 allCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 87 fListCmd = new G4UIcmdWithAString("/testem/p << 73 88 fListCmd->SetGuidance("Add modula physics li << 74 pListCmd = new G4UIcmdWithAString("/testem/phys/addPhysics",this); 89 fListCmd->SetParameterName("PList", false); << 75 pListCmd->SetGuidance("Add modula physics list."); 90 fListCmd->AvailableForStates(G4State_PreInit << 76 pListCmd->SetParameterName("PList",false); 91 << 77 pListCmd->AvailableForStates(G4State_PreInit); 92 fADCCmd = new G4UIcmdWithADoubleAndUnit("/te << 93 fADCCmd->SetGuidance("Set energy per ADC cha << 94 fADCCmd->SetParameterName("enadc", false, fa << 95 fADCCmd->SetUnitCategory("Energy"); << 96 fADCCmd->SetDefaultUnit("keV"); << 97 fADCCmd->SetRange("enadc>0."); << 98 fADCCmd->AvailableForStates(G4State_PreInit, << 99 << 100 fNorCmd = new G4UIcmdWithADouble("/testem/ph << 101 fNorCmd->SetGuidance("Set factor for histogr << 102 fNorCmd->SetParameterName("nfac", false, fal << 103 fNorCmd->SetRange("nfac>0."); << 104 fNorCmd->AvailableForStates(G4State_PreInit, << 105 << 106 fSmCmd = new G4UIcmdWithADoubleAndUnit("/tes << 107 fSmCmd->SetGuidance("Set intrinsic width of << 108 fSmCmd->SetParameterName("sm", false, false) << 109 fSmCmd->SetDefaultUnit("keV"); << 110 fSmCmd->SetRange("sm>0."); << 111 fSmCmd->AvailableForStates(G4State_PreInit, << 112 } 78 } 113 79 114 //....oooOO0OOooo........oooOO0OOooo........oo 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 115 81 116 PhysicsListMessenger::~PhysicsListMessenger() 82 PhysicsListMessenger::~PhysicsListMessenger() 117 { 83 { 118 delete fECmd; << 84 delete gammaCutCmd; 119 delete fEBCmd; << 85 delete electCutCmd; 120 delete fCBCmd; << 86 delete protoCutCmd; 121 delete fCMCmd; << 87 delete allCutCmd; 122 delete fListCmd; << 88 delete pListCmd; 123 delete fADCCmd; << 124 delete fNorCmd; << 125 delete fSmCmd; << 126 delete fPhysDir; << 127 } 89 } 128 90 129 //....oooOO0OOooo........oooOO0OOooo........oo 91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 130 92 131 void PhysicsListMessenger::SetNewValue(G4UIcom << 93 void PhysicsListMessenger::SetNewValue(G4UIcommand* command, 132 { << 94 G4String newValue) 133 TestParameters* man = TestParameters::GetPoi << 95 { >> 96 if( command == gammaCutCmd ) >> 97 { pPhysicsList->SetCutForGamma(gammaCutCmd->GetNewDoubleValue(newValue));} >> 98 >> 99 if( command == electCutCmd ) >> 100 { pPhysicsList->SetCutForElectron(electCutCmd->GetNewDoubleValue(newValue));} >> 101 >> 102 if( command == protoCutCmd ) >> 103 { pPhysicsList->SetCutForPositron(protoCutCmd->GetNewDoubleValue(newValue));} >> 104 >> 105 if( command == allCutCmd ) >> 106 { >> 107 G4double cut = allCutCmd->GetNewDoubleValue(newValue); >> 108 pPhysicsList->SetCutForGamma(cut); >> 109 pPhysicsList->SetCutForElectron(cut); >> 110 pPhysicsList->SetCutForPositron(cut); >> 111 } 134 112 135 if (command == fECmd) { << 113 if( command == pListCmd ) 136 man->SetMaxEnergy(fECmd->GetNewDoubleValue << 114 { pPhysicsList->AddPhysicsList(newValue);} 137 } << 138 if (command == fEBCmd) { << 139 man->SetNumberBins(fEBCmd->GetNewIntValue( << 140 } << 141 if (command == fCBCmd) { << 142 man->SetNumberBinsCluster(fCBCmd->GetNewIn << 143 } << 144 if (command == fCMCmd) { << 145 man->SetMaxCluster(fCMCmd->GetNewIntValue( << 146 } << 147 if (command == fListCmd) { << 148 fPhysicsList->AddPhysicsList(newValue); << 149 } << 150 if (command == fADCCmd) { << 151 man->SetEnergyPerChannel(fADCCmd->GetNewDo << 152 } << 153 if (command == fNorCmd) { << 154 man->SetNormFactor(fNorCmd->GetNewDoubleVa << 155 } << 156 if (command == fSmCmd) { << 157 man->SetEnergySmear(fSmCmd->GetNewDoubleVa << 158 } << 159 } 115 } 160 116 161 //....oooOO0OOooo........oooOO0OOooo........oo 117 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 162 118