Changeset 1342 for trunk/examples/extended/electromagnetic/TestEm8
- Timestamp:
- Nov 5, 2010, 4:08:39 PM (14 years ago)
- Location:
- trunk/examples/extended/electromagnetic/TestEm8
- Files:
-
- 11 edited
-
GNUmakefile (modified) (2 diffs)
-
History (modified) (2 diffs)
-
README (modified) (4 diffs)
-
TestEm8.cc (modified) (6 diffs)
-
TestEm8.in (modified) (2 diffs)
-
TestEm8.large_N.in (modified) (2 diffs)
-
TestEm8.out (modified) (36 diffs)
-
include/PhysicsList.hh (modified) (4 diffs)
-
include/PhysicsListMessenger.hh (modified) (3 diffs)
-
src/PhysicsList.cc (modified) (8 diffs)
-
src/PhysicsListMessenger.cc (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/electromagnetic/TestEm8/GNUmakefile
r1230 r1342 1 # $Id: GNUmakefile,v 1. 8 2008/06/11 22:27:45 maireExp $1 # $Id: GNUmakefile,v 1.9 2010/09/28 19:07:32 vnivanch Exp $ 2 2 # -------------------------------------------------------------- 3 3 # GNUmakefile for examples module. Gabriele Cosmo, 06/04/98. … … 20 20 21 21 include $(G4INSTALL)/config/binmake.gmk 22 23 histclean: 24 rm ${G4WORKDIR}/tmp/${G4SYSTEM}/${G4TARGET}/Histo.o 25 -
trunk/examples/extended/electromagnetic/TestEm8/History
r1337 r1342 1 $Id: History,v 1. 37 2010/06/06 04:31:28 perlExp $1 $Id: History,v 1.40 2010/10/27 14:52:07 vnivanch Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 27-10-10 V.Ivant (testem8-V09-03-06) 19 - DetectorConstruction - added more printout 20 - DetectorMessenger - fixed state (added G4State_PreInit) 21 - PhysicsList - fixed UI commands for step limiter 22 23 28-09-10 V.Ivant (testem8-V09-03-05) 24 - GNUmakefile - added histclean option 25 26 08-09-10 V.Ivant (testem8-V09-03-04) 27 - The example is completely rewritten: physics, scoring, histograms, 28 README 17 29 18 30 04-06-10 J.Perl (testem8-V09-03-03) -
trunk/examples/extended/electromagnetic/TestEm8/README
r807 r1342 1 $Id: README,v 1. 4 2003/11/24 18:09:22vnivanch Exp $1 $Id: README,v 1.5 2010/09/08 11:22:01 vnivanch Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 9 9 ------- 10 10 11 Test for investigation of ionisation in thin absorbers, transition 12 and synchrotron radiations. 13 11 Example for investigation of ionisation in thin absorbers and gaseous 12 detectors 14 13 15 14 1- GEOMETRY DEFINITION 16 15 17 The "absorber" is a tube made of a given material. 16 The target is a cilinder made of a given material placed inside 17 cilindrical container, which is placed inside the world volume. 18 18 19 Three parameters define the absorber : 20 - the material of the absorber, 21 - the thickness of an absorber, 22 - the transverse size of the absorber (the input face is a square). 23 24 The volume "World" contains the "absorber". 25 In this test the parameters of the "World" can be changed , too. 19 Following parameters define the geometry: 20 - the material of the target, 21 - the thickness of the target, 22 - the radius of the target, 23 - the material of the container, 24 - the thickness of the container, 25 - the material of the world. 26 27 The list of materials used in gaseous detectors are built inside 28 the DetectorConstruction class, also NIST materials are availabe. 29 The default geometry is provided but all parameters can be changed via 30 UI commands defined in the DetectorMessenger class, for example, 26 31 27 In addition a transverse uniform magnetic field can be applied. 28 29 The default geometry is constructed in DetectorConstruction class, 30 but all the parameters can be changed via 31 the commands defined in the DetectorMessenger class. 32 /testem/setGasMat XeCH4C3H8 33 /testem/setWindowMat G4_MYLAR 34 /testem/setWorldMat G4_AIR 35 /testem/setGasThick 10 cm 36 /testem/setGasRad 20 cm 37 /testem/setWindowThick 50 um 38 32 39 33 40 2- AN EVENT : THE PRIMARY GENERATOR … … 35 42 The primary kinematic consists of a single particle which hits the 36 43 absorber perpendicular to the input face. The type of the particle 37 and its energy are set in the PrimaryGeneratorAction class, and can 38 be changed via the G4 build-in commands of ParticleGun class (see 39 the macros provided with this example). 40 44 and its energy can be set via the G4 build-in commands of ParticleGun 41 45 A RUN is a set of events. 42 46 43 47 3- DETECTOR RESPONSE 44 48 45 Here we test G4PAIionisation , G4IonisationByLogicalVolume and 46 transition radiation processes 47 48 A HIT is a record, event per event , of all the 49 informations needed to simulate and analyse the detector response. 50 51 In this example a CalorHit is defined as a set of 2 informations: 52 - the total energy deposit in the absorber, 53 - the total tracklength of all charged particles in the absorber, 54 55 Therefore the absorber is declared 56 'sensitive detector' (SD), which means they can contribute to the hit. 57 58 At the end of a run, from the histogram(s), one can study 59 different physics quantities such as : 60 - angle distribution, 61 - energy deposit, 62 - transmission/backscattering, 63 - ... 64 65 The test contains 10 built-in histograms, which can be activated by 66 interactive commands (see the macros runxx.mac for details). 49 The TargetSD class sending information about each step inside the target 50 to the HistoManager class scoring of energy deposition in the detector. 51 Additionally at each step of a particle inside the target the number of 52 ionisation clusters is sampled using G4ElectronIonPair helper class. The 53 parameter of transformation of energy into ionisation clusters can be 54 set via UI command: 67 55 68 The histogram files can be viewed using PAW e.g with the commands 69 70 paw> h/file 1 geant4.plot01 or g4.p11 71 paw> option stat 72 paw> h/pl 1 73 56 /testem/setPairEnergy 19 eV 74 57 75 58 76 4- PHYSICS DEMO59 4- PHYSICS 77 60 78 61 The particle's type and the physic processes which will be available 79 in this example are set in PhysicsList class. 80 81 The messenger classes introduce interactive commands . Using these 82 commands the geometry of the detector, the data of the primary 83 particle, the limits of the histograms , etc. can be changed. 84 62 in this example are set in PhysicsList class, which uses Geant4 63 EM physics constructors providing in the physics_list library. 64 65 The messenger classes introduce interactive commands. In particular, 66 PAI ionisation model can be added using G4EmConfurator helper class, 67 which is invokated by the UI command 68 69 /testem/phys/addPhysics pai 85 70 86 71 5- HOW TO START ? … … 91 76 92 77 - execute TestEm8 in 'batch' mode from macro files e.g. 93 % $(G4INSTALL)/bin/$(G4SYSTEM)/TestEm8 run11.mac78 % $(G4INSTALL)/bin/$(G4SYSTEM)/TestEm8 TestEm8.in 94 79 95 80 - execute TestEm8 in 'interactive' mode with visualization e.g. -
trunk/examples/extended/electromagnetic/TestEm8/TestEm8.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // 27 // $Id: TestEm8.cc,v 1.9 2010/05/21 18:15:04 maire Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 26 // $Id: TestEm8.cc,v 1.10 2010/09/08 09:12:10 vnivanch Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 29 28 // 30 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 35 34 #include "Randomize.hh" 36 35 37 #include " Em8DetectorConstruction.hh"36 #include "DetectorConstruction.hh" 38 37 #include "PhysicsList.hh" 39 #include "Em8PrimaryGeneratorAction.hh" 40 #include "Em8RunAction.hh" 41 #include "Em8EventAction.hh" 42 #include "Em8SteppingAction.hh" 43 #include "Em8SteppingVerbose.hh" 38 #include "PrimaryGeneratorAction.hh" 39 #include "RunAction.hh" 40 #include "EventAction.hh" 44 41 45 42 #ifdef G4VIS_USE … … 55 52 int main(int argc,char** argv) 56 53 { 57 58 54 //choose the Random engine 59 60 55 CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine); 61 56 62 //my Verbose output class63 64 G4VSteppingVerbose::SetInstance(new Em8SteppingVerbose);65 66 57 // Construct the default run manager 67 68 58 G4RunManager * runManager = new G4RunManager; 69 59 70 60 // set mandatory initialization classes 71 72 Em8DetectorConstruction* detector;73 detector = new Em8DetectorConstruction;74 runManager->SetUserInitialization(detector);75 61 runManager->SetUserInitialization(new PhysicsList); 62 PrimaryGeneratorAction* gun = new PrimaryGeneratorAction(); 63 runManager->SetUserInitialization(new DetectorConstruction(gun)); 76 64 77 65 // set user action classes 78 79 runManager->SetUserAction(new Em8PrimaryGeneratorAction(detector)); 80 81 Em8RunAction* runAction = new Em8RunAction; 82 83 runManager->SetUserAction(runAction); 84 85 Em8EventAction* eventAction = new Em8EventAction(runAction); 86 87 runManager->SetUserAction(eventAction); 88 89 Em8SteppingAction* steppingAction = new Em8SteppingAction(detector, 90 eventAction, 91 runAction); 92 runManager->SetUserAction(steppingAction); 66 runManager->SetUserAction(gun); 67 runManager->SetUserAction(new RunAction()); 68 runManager->SetUserAction(new EventAction()); 93 69 94 70 G4UImanager* UI = G4UImanager::GetUIpointer(); … … 96 72 if (argc!=1) // batch mode 97 73 { 98 G4String command = "/control/execute ";99 G4String fileName = argv[1];100 UI->ApplyCommand(command+fileName);74 G4String command = "/control/execute "; 75 G4String fileName = argv[1]; 76 UI->ApplyCommand(command+fileName); 101 77 } 102 78 … … 104 80 { 105 81 #ifdef G4VIS_USE 106 G4VisManager* visManager = new G4VisExecutive;107 visManager->Initialize();82 G4VisManager* visManager = new G4VisExecutive; 83 visManager->Initialize(); 108 84 #endif 109 85 … … 115 91 116 92 #ifdef G4VIS_USE 117 delete visManager;93 delete visManager; 118 94 #endif 119 95 } -
trunk/examples/extended/electromagnetic/TestEm8/TestEm8.in
r1337 r1342 13 13 #/testem/phys/addPhysics pai_photon 14 14 # 15 /testem/phys/setCuts 1 mm16 15 #/testem/stepMax 10 mm 17 16 # 18 # PLOT: energy deposit distribution in absorber 19 # 20 /plots/setnbinEn 20 21 /plots/setEnlow 0.0 MeV 22 /plots/setEnhigh 0.060 MeV 23 /plots/sethistName g4.p11 17 /testem/phys/setNbinsE 20 18 /testem/phys/setMaxE 0.060 MeV 24 19 # 25 20 /run/initialize … … 37 32 /gun/energy 200000.00 MeV 38 33 # 39 #/gun/energy 200.0 MeV40 #41 #/tracking/verbose 142 #43 /event/printModulo 10044 34 /run/beamOn 1000 45 35 # -
trunk/examples/extended/electromagnetic/TestEm8/TestEm8.large_N.in
r807 r1342 8 8 #/tracking/verbose 2 9 9 # 10 #11 10 # select EM PhysicsList 12 11 # 13 /testem/phys/addPhysics standard 12 /testem/phys/addPhysics pai 13 #/testem/phys/addPhysics standard 14 14 # 15 #/testem/phys/addPhysics model 16 # 17 /testem/phys/setCuts 1 mm 18 #/testem/phys/VertexCuts 0.1 mm 19 #/testem/phys/MuonCuts 10.0 mm 20 #/testem/stepMax 10 mm 21 # 22 #/process/eLoss/preciseRange true 23 #/process/eLoss/integral true 24 # 25 # hist file name 26 # 27 /plots/sethistName g4.p11 28 # 29 # PLOT: energy deposit distribution in absorber 30 # 31 /plots/setnbinEn 20 32 /plots/setEnlow 0.0 MeV 33 /plots/setEnhigh 0.060 MeV 34 # 15 /testem/phys/setNbinsE 20 16 /testem/phys/setMaxE 0.060 MeV 35 17 # 36 18 /run/initialize 19 # 20 /run/setCut 1 mm 21 /run/setCutForRegion GasDetector 0.2 mm 37 22 /run/particle/dumpCutValues 38 #39 #/gun/particle e-40 23 # 41 24 /gun/particle proton … … 49 32 #/tracking/verbose 1 50 33 # 51 /event/printModulo 100052 34 /run/beamOn 10000 53 #54 35 # 55 36 -
trunk/examples/extended/electromagnetic/TestEm8/TestEm8.out
r1337 r1342 5 5 6 6 ************************************************************* 7 Geant4 version Name: geant4-09-03-ref-0 6(25-June-2010)7 Geant4 version Name: geant4-09-03-ref-09 (25-June-2010) 8 8 Copyright : Geant4 Collaboration 9 9 Reference : NIM A 506 (2003), 250-303 … … 11 11 ************************************************************* 12 12 13 14 ***** Table : Nb of materials = 19 ***** 15 16 Material: G4_Ar density: 1.662 mg/cm3 RadL: 117.621 m Nucl.Int.Length: 719.260 m Imean: 188.000 eV temperature: 273.15 K pressure: 1.00 atm 17 ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole 18 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 % 19 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 % 20 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 21 22 Material: G4_Kr density: 3.478 mg/cm3 RadL: 32.694 m Nucl.Int.Length: 439.629 m Imean: 352.000 eV temperature: 273.15 K pressure: 1.00 atm 23 ---> Element: Kr (Kr) Z = 36.0 N = 83.9 A = 83.80 g/mole 24 ---> Isotope: Kr78 Z = 36 N = 78 A = 77.92 g/mole abundance: 0.35 % 25 ---> Isotope: Kr80 Z = 36 N = 80 A = 79.92 g/mole abundance: 2.28 % 26 ---> Isotope: Kr82 Z = 36 N = 82 A = 81.91 g/mole abundance: 11.58 % 27 ---> Isotope: Kr83 Z = 36 N = 83 A = 82.91 g/mole abundance: 11.49 % 28 ---> Isotope: Kr84 Z = 36 N = 84 A = 83.91 g/mole abundance: 57.00 % 29 ---> Isotope: Kr86 Z = 36 N = 86 A = 85.91 g/mole abundance: 17.30 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 30 31 Material: G4_Xe density: 5.485 mg/cm3 RadL: 15.462 m Nucl.Int.Length: 324.779 m Imean: 482.000 eV temperature: 273.15 K pressure: 1.00 atm 32 ---> Element: Xe (Xe) Z = 54.0 N = 131.4 A = 131.29 g/mole 33 ---> Isotope: Xe124 Z = 54 N = 124 A = 123.91 g/mole abundance: 0.09 % 34 ---> Isotope: Xe126 Z = 54 N = 126 A = 125.90 g/mole abundance: 0.09 % 35 ---> Isotope: Xe128 Z = 54 N = 128 A = 127.90 g/mole abundance: 1.92 % 36 ---> Isotope: Xe129 Z = 54 N = 129 A = 128.91 g/mole abundance: 26.44 % 37 ---> Isotope: Xe130 Z = 54 N = 130 A = 129.90 g/mole abundance: 4.08 % 38 ---> Isotope: Xe131 Z = 54 N = 131 A = 130.90 g/mole abundance: 21.18 % 39 ---> Isotope: Xe132 Z = 54 N = 132 A = 131.90 g/mole abundance: 26.89 % 40 ---> Isotope: Xe134 Z = 54 N = 134 A = 133.91 g/mole abundance: 10.44 % 41 ---> Isotope: Xe136 Z = 54 N = 136 A = 135.91 g/mole abundance: 8.87 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 42 43 Material: G4_CARBON_DIOXIDE CO_2 density: 1.842 mg/cm3 RadL: 196.486 m Nucl.Int.Length: 466.077 m Imean: 85.000 eV temperature: 273.15 K pressure: 1.00 atm 44 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 45 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 46 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 27.29 % ElmAbundance 33.33 % 47 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 48 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 49 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 50 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 72.71 % ElmAbundance 66.67 % 51 52 Material: G4_MYLAR density: 1.400 g/cm3 RadL: 28.536 cm Nucl.Int.Length: 55.985 cm Imean: 78.700 eV 53 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 54 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 55 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 4.20 % ElmAbundance 36.36 % 56 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 57 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 58 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 62.50 % ElmAbundance 45.46 % 59 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 60 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 61 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 62 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 33.30 % ElmAbundance 18.18 % 63 64 Material: G4_METHANE density: 0.667 kg/m3 RadL: 696.479 m Nucl.Int.Length: 910.708 m Imean: 41.700 eV temperature: 273.15 K pressure: 1.00 atm 65 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 66 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 67 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 25.13 % ElmAbundance 80.00 % 68 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 69 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 70 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 74.87 % ElmAbundance 20.00 % 71 72 Material: G4_PROPANE density: 1.879 mg/cm3 RadL: 241.434 m Nucl.Int.Length: 346.143 m Imean: 47.100 eV temperature: 273.15 K pressure: 1.00 atm 73 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 74 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 75 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 18.29 % ElmAbundance 72.73 % 76 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 77 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 78 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 81.71 % ElmAbundance 27.27 % 79 80 Material: G4_Galactic density: 0.000 kg/m3 RadL: 204310101.835 pc Nucl.Int.Length: 114327973.206 pc Imean: 21.800 eV temperature: 2.73 K pressure: 0.00 atm 81 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 82 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 83 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 84 85 Material: Kr7CH4 density: 3.491 mg/cm3 RadL: 32.924 m Nucl.Int.Length: 428.924 m Imean: 337.765 eV temperature: 273.15 K pressure: 1.00 atm 86 ---> Element: Kr (Kr) Z = 36.0 N = 83.9 A = 83.80 g/mole 87 ---> Isotope: Kr78 Z = 36 N = 78 A = 77.92 g/mole abundance: 0.35 % 88 ---> Isotope: Kr80 Z = 36 N = 80 A = 79.92 g/mole abundance: 2.28 % 89 ---> Isotope: Kr82 Z = 36 N = 82 A = 81.91 g/mole abundance: 11.58 % 90 ---> Isotope: Kr83 Z = 36 N = 83 A = 82.91 g/mole abundance: 11.49 % 91 ---> Isotope: Kr84 Z = 36 N = 84 A = 83.91 g/mole abundance: 57.00 % 92 ---> Isotope: Kr86 Z = 36 N = 86 A = 85.91 g/mole abundance: 17.30 % ElmMassFraction: 98.60 % ElmAbundance 72.95 % 93 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 94 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 95 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 0.35 % ElmAbundance 21.64 % 96 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 97 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 98 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 1.05 % ElmAbundance 5.41 % 99 100 Material: TRT_Xe density: 5.485 mg/cm3 RadL: 15.463 m Nucl.Int.Length: 324.801 m Imean: 482.000 eV temperature: 293.15 K pressure: 1.00 atm 101 ---> Element: Xe (Xe) Z = 54.0 N = 131.4 A = 131.29 g/mole 102 ---> Isotope: Xe124 Z = 54 N = 124 A = 123.91 g/mole abundance: 0.09 % 103 ---> Isotope: Xe126 Z = 54 N = 126 A = 125.90 g/mole abundance: 0.09 % 104 ---> Isotope: Xe128 Z = 54 N = 128 A = 127.90 g/mole abundance: 1.92 % 105 ---> Isotope: Xe129 Z = 54 N = 129 A = 128.91 g/mole abundance: 26.44 % 106 ---> Isotope: Xe130 Z = 54 N = 130 A = 129.90 g/mole abundance: 4.08 % 107 ---> Isotope: Xe131 Z = 54 N = 131 A = 130.90 g/mole abundance: 21.18 % 108 ---> Isotope: Xe132 Z = 54 N = 132 A = 131.90 g/mole abundance: 26.89 % 109 ---> Isotope: Xe134 Z = 54 N = 134 A = 133.91 g/mole abundance: 10.44 % 110 ---> Isotope: Xe136 Z = 54 N = 136 A = 135.91 g/mole abundance: 8.87 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 111 112 Material: TRT_CO2 density: 1.842 mg/cm3 RadL: 196.499 m Nucl.Int.Length: 466.107 m Imean: 90.958 eV temperature: 293.15 K pressure: 1.00 atm 113 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 114 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 115 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 27.29 % ElmAbundance 33.33 % 116 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 117 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 118 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 119 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 72.71 % ElmAbundance 66.67 % 120 121 Material: TRT_CF4 density: 3.900 mg/cm3 RadL: 87.158 m Nucl.Int.Length: 234.197 m Imean: 109.384 eV temperature: 293.15 K pressure: 1.00 atm 122 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 123 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 124 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 13.65 % ElmAbundance 20.00 % 125 ---> Element: F (F) Z = 9.0 N = 19.0 A = 19.00 g/mole 126 ---> Isotope: F19 Z = 9 N = 19 A = 19.00 g/mole abundance: 100.00 % ElmMassFraction: 86.35 % ElmAbundance 80.00 % 127 128 Material: XeCO2CF4 density: 4.760 mg/cm3 RadL: 20.851 m Nucl.Int.Length: 314.963 m Imean: 345.469 eV temperature: 293.15 K pressure: 1.00 atm 129 ---> Element: Xe (Xe) Z = 54.0 N = 131.4 A = 131.29 g/mole 130 ---> Isotope: Xe124 Z = 54 N = 124 A = 123.91 g/mole abundance: 0.09 % 131 ---> Isotope: Xe126 Z = 54 N = 126 A = 125.90 g/mole abundance: 0.09 % 132 ---> Isotope: Xe128 Z = 54 N = 128 A = 127.90 g/mole abundance: 1.92 % 133 ---> Isotope: Xe129 Z = 54 N = 129 A = 128.91 g/mole abundance: 26.44 % 134 ---> Isotope: Xe130 Z = 54 N = 130 A = 129.90 g/mole abundance: 4.08 % 135 ---> Isotope: Xe131 Z = 54 N = 131 A = 130.90 g/mole abundance: 21.18 % 136 ---> Isotope: Xe132 Z = 54 N = 132 A = 131.90 g/mole abundance: 26.89 % 137 ---> Isotope: Xe134 Z = 54 N = 134 A = 133.91 g/mole abundance: 10.44 % 138 ---> Isotope: Xe136 Z = 54 N = 136 A = 135.91 g/mole abundance: 8.87 % ElmMassFraction: 80.70 % ElmAbundance 35.01 % 139 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 140 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 141 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 3.17 % ElmAbundance 15.02 % 142 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 143 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 144 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 145 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 2.84 % ElmAbundance 10.10 % 146 ---> Element: F (F) Z = 9.0 N = 19.0 A = 19.00 g/mole 147 ---> Isotope: F19 Z = 9 N = 19 A = 19.00 g/mole abundance: 100.00 % ElmMassFraction: 13.30 % ElmAbundance 39.87 % 148 149 Material: C3H8 density: 3.758 mg/cm3 RadL: 120.742 m Nucl.Int.Length: 173.106 m Imean: 52.014 eV temperature: 273.15 K pressure: 1.00 atm 150 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 151 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 152 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 81.71 % ElmAbundance 27.27 % 153 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 154 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 155 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 18.29 % ElmAbundance 72.73 % 156 157 Material: XeCH4C3H8 density: 4.920 mg/cm3 RadL: 17.657 m Nucl.Int.Length: 344.114 m Imean: 438.263 eV temperature: 273.15 K pressure: 1.00 atm 158 ---> Element: Xe (Xe) Z = 54.0 N = 131.4 A = 131.29 g/mole 159 ---> Isotope: Xe124 Z = 54 N = 124 A = 123.91 g/mole abundance: 0.09 % 160 ---> Isotope: Xe126 Z = 54 N = 126 A = 125.90 g/mole abundance: 0.09 % 161 ---> Isotope: Xe128 Z = 54 N = 128 A = 127.90 g/mole abundance: 1.92 % 162 ---> Isotope: Xe129 Z = 54 N = 129 A = 128.91 g/mole abundance: 26.44 % 163 ---> Isotope: Xe130 Z = 54 N = 130 A = 129.90 g/mole abundance: 4.08 % 164 ---> Isotope: Xe131 Z = 54 N = 131 A = 130.90 g/mole abundance: 21.18 % 165 ---> Isotope: Xe132 Z = 54 N = 132 A = 131.90 g/mole abundance: 26.89 % 166 ---> Isotope: Xe134 Z = 54 N = 134 A = 133.91 g/mole abundance: 10.44 % 167 ---> Isotope: Xe136 Z = 54 N = 136 A = 135.91 g/mole abundance: 8.87 % ElmMassFraction: 97.10 % ElmAbundance 48.49 % 168 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 169 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 170 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 0.60 % ElmAbundance 38.95 % 171 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 172 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 173 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 2.30 % ElmAbundance 12.56 % 174 175 Material: Ar7CH4 density: 1.709 mg/cm3 RadL: 116.342 m Nucl.Int.Length: 680.394 m Imean: 177.716 eV temperature: 273.15 K pressure: 1.00 atm 176 ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole 177 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 % 178 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 % 179 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 97.10 % ElmAbundance 72.89 % 180 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 181 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 182 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 0.73 % ElmAbundance 21.68 % 183 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 184 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 185 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 2.17 % ElmAbundance 5.42 % 186 187 Material: ArCO2 density: 1.822 mg/cm3 RadL: 119.168 m Nucl.Int.Length: 604.526 m Imean: 158.492 eV temperature: 273.15 K pressure: 1.00 atm 188 ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole 189 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 % 190 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 % 191 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 78.30 % ElmAbundance 56.99 % 192 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 193 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 194 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 5.92 % ElmAbundance 14.34 % 195 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 196 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 197 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 198 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 15.78 % ElmAbundance 28.67 % 199 200 Material: Xe20CO2 density: 5.082 mg/cm3 RadL: 17.750 m Nucl.Int.Length: 323.450 m Imean: 412.596 eV temperature: 273.15 K pressure: 1.00 atm 201 ---> Element: Xe (Xe) Z = 54.0 N = 131.4 A = 131.29 g/mole 202 ---> Isotope: Xe124 Z = 54 N = 124 A = 123.91 g/mole abundance: 0.09 % 203 ---> Isotope: Xe126 Z = 54 N = 126 A = 125.90 g/mole abundance: 0.09 % 204 ---> Isotope: Xe128 Z = 54 N = 128 A = 127.90 g/mole abundance: 1.92 % 205 ---> Isotope: Xe129 Z = 54 N = 129 A = 128.91 g/mole abundance: 26.44 % 206 ---> Isotope: Xe130 Z = 54 N = 130 A = 129.90 g/mole abundance: 4.08 % 207 ---> Isotope: Xe131 Z = 54 N = 131 A = 130.90 g/mole abundance: 21.18 % 208 ---> Isotope: Xe132 Z = 54 N = 132 A = 131.90 g/mole abundance: 26.89 % 209 ---> Isotope: Xe134 Z = 54 N = 134 A = 133.91 g/mole abundance: 10.44 % 210 ---> Isotope: Xe136 Z = 54 N = 136 A = 135.91 g/mole abundance: 8.87 % ElmMassFraction: 92.20 % ElmAbundance 56.91 % 211 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 212 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 213 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 2.13 % ElmAbundance 14.36 % 214 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 215 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 216 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 217 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 5.67 % ElmAbundance 28.73 % 218 219 Material: Kr20CO2 density: 3.601 mg/cm3 RadL: 34.157 m Nucl.Int.Length: 391.053 m Imean: 296.925 eV temperature: 273.15 K pressure: 1.00 atm 220 ---> Element: Kr (Kr) Z = 36.0 N = 83.9 A = 83.80 g/mole 221 ---> Isotope: Kr78 Z = 36 N = 78 A = 77.92 g/mole abundance: 0.35 % 222 ---> Isotope: Kr80 Z = 36 N = 80 A = 79.92 g/mole abundance: 2.28 % 223 ---> Isotope: Kr82 Z = 36 N = 82 A = 81.91 g/mole abundance: 11.58 % 224 ---> Isotope: Kr83 Z = 36 N = 83 A = 82.91 g/mole abundance: 11.49 % 225 ---> Isotope: Kr84 Z = 36 N = 84 A = 83.91 g/mole abundance: 57.00 % 226 ---> Isotope: Kr86 Z = 36 N = 86 A = 85.91 g/mole abundance: 17.30 % ElmMassFraction: 89.00 % ElmAbundance 58.62 % 227 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 228 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 229 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 3.00 % ElmAbundance 13.79 % 230 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 231 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 232 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 233 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 8.00 % ElmAbundance 27.59 % 234 235 13 236 /run/verbose 2 14 237 #/tracking/verbose 2 … … 19 242 #/testem/phys/addPhysics pai_photon 20 243 # 21 /testem/phys/setCuts 1 mm22 244 #/testem/stepMax 10 mm 23 245 # 24 # PLOT: energy deposit distribution in absorber 25 # 26 /plots/setnbinEn 20 27 Nb of bins in Edep plot = 20 28 /plots/setEnlow 0.0 MeV 29 Elow in the Edep plot = 0 30 /plots/setEnhigh 0.060 MeV 31 Ehigh in the Edep plot = 0.06 32 /plots/sethistName g4.p11 33 hist file = g4.p11 246 /testem/phys/setNbinsE 20 247 /testem/phys/setMaxE 0.060 MeV 34 248 # 35 249 /run/initialize 36 250 userDetector->Construct() start. 37 251 38 The WORLD is made of 2 3.1022mm of Mylar, the transverse size (R) of the world is 100mm.39 The ABSORBER is made of 23mm of XeCH4C3H8, the transverse size (R) is 100 mm.40 Z position of the (middle of the) absorber 0 mm.252 The WORLD is made of 27.7224mm of G4_Galactic, the transverse size (R) of the world is 120.122 mm. 253 The CONTAINER is made of 0.051mm of G4_MYLAR 254 The TARGET is made of 23mm of XeCH4C3H8, the transverse size (R) is 100 mm. 41 255 42 256 World is registered to the default region. … … 59 273 /gun/energy 200000.00 MeV 60 274 # 61 #/gun/energy 200.0 MeV62 #63 #/tracking/verbose 164 #65 /event/printModulo 10066 275 /run/beamOn 1000 67 G4PAIModel::Initialise for e+68 G4PAIModel::Initialise for e-69 G4PAIModel::Initialise for mu+70 G4PAIModel::Initialise for mu-71 G4PAIModel::Initialise for pi+72 G4PAIModel::Initialise for pi-73 G4PAIModel::Initialise for proton74 276 75 277 phot: for gamma SubType= 12 76 278 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 77 PhotoElectric : Emin= 0 eV Emax=10 TeV279 PhotoElectric : Emin= 0 eV Emax= 10 TeV 78 280 79 281 compt: for gamma SubType= 13 80 282 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 81 283 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 82 Klein-Nishina : Emin= 0 eV Emax=10 TeV284 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 83 285 84 286 conv: for gamma SubType= 14 85 287 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 86 288 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 87 Bethe-Heitler : Emin= 0 eV Emax=10 TeV289 BetheHeitler : Emin= 0 eV Emax= 10 TeV 88 290 89 291 msc: for e- SubType= 10 … … 91 293 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1 92 294 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 93 UrbanMsc93 : Emin= 0 eV Emax=10 TeV295 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 94 296 95 297 eIoni: for e- SubType= 2 … … 98 300 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 99 301 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 100 MollerBhabha : Emin= 0 eV Emax=10 TeV302 MollerBhabha : Emin= 0 eV Emax= 10 TeV 101 303 ===== EM models for the G4Region GasDetector ====== 102 PAIModel : Emin= 0 eV Emax=10 TeV304 PAIModel : Emin= 0 eV Emax= 10 TeV 103 305 104 306 eBrem: for e- SubType= 3 … … 107 309 LPM flag: 1 for E > 1 GeV 108 310 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 109 eBrem : Emin= 0 eV Emax= 1 GeV110 eBrem Rel : Emin= 1 GeV Emax= 10 TeV311 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 312 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 111 313 112 314 eIoni: for e+ SubType= 2 … … 115 317 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 116 318 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 117 MollerBhabha : Emin= 0 eV Emax=10 TeV319 MollerBhabha : Emin= 0 eV Emax= 10 TeV 118 320 ===== EM models for the G4Region GasDetector ====== 119 PAIModel : Emin= 0 eV Emax=10 TeV321 PAIModel : Emin= 0 eV Emax= 10 TeV 120 322 121 323 eBrem: for e+ SubType= 3 … … 124 326 LPM flag: 1 for E > 1 GeV 125 327 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 126 eBrem : Emin= 0 eV Emax= 1 GeV127 eBrem Rel : Emin= 1 GeV Emax= 10 TeV328 eBrem : Emin= 0 eV Emax= 1 GeV AngularGenUrban 329 eBremLPM : Emin= 1 GeV Emax= 10 TeV AngularGenUrban 128 330 129 331 annihil: for e+ SubType= 5 130 332 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 131 333 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 132 eplus2gg : Emin= 0 eV Emax=10 TeV334 eplus2gg : Emin= 0 eV Emax= 10 TeV 133 335 134 336 msc: for proton SubType= 10 … … 136 338 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 137 339 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 138 UrbanMsc90 : Emin= 0 eV Emax=10 TeV340 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 139 341 140 342 hIoni: for proton SubType= 2 … … 143 345 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 144 346 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 145 Bragg : Emin= 0 eV Emax=2 MeV146 BetheBloch : Emin= 2 MeV Emax=10 TeV347 Bragg : Emin= 0 eV Emax= 2 MeV 348 BetheBloch : Emin= 2 MeV Emax= 10 TeV 147 349 ===== EM models for the G4Region GasDetector ====== 148 PAIModel : Emin= 0 eV Emax=10 TeV350 PAIModel : Emin= 0 eV Emax= 10 TeV 149 351 150 352 hBrems: for proton SubType= 3 … … 152 354 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 153 355 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 154 hBrem : Emin= 0 eV Emax=10 TeV356 hBrem : Emin= 0 eV Emax= 10 TeV 155 357 156 358 hPairProd: for proton SubType= 4 … … 158 360 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 159 361 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 160 hPairProd : Emin= 0 eV Emax=10 TeV362 hPairProd : Emin= 0 eV Emax= 10 TeV 161 363 162 364 msc: for GenericIon SubType= 10 163 365 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 164 366 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 165 UrbanMsc90 : Emin= 0 eV Emax=10 TeV367 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 166 368 167 369 ionIoni: for GenericIon SubType= 2 … … 169 371 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 170 372 finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02 171 Stopping Power data for 17 ion/material pairs , nuclearStopping: 1172 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 173 BraggIon : Emin= 0 eV Emax=2 MeV174 BetheBloch : Emin= 2 MeV Emax=10 TeV373 Stopping Power data for 17 ion/material pairs 374 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 375 BraggIon : Emin= 0 eV Emax= 2 MeV 376 BetheBloch : Emin= 2 MeV Emax= 10 TeV 175 377 176 378 hIoni: for anti_proton SubType= 2 … … 179 381 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 180 382 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 181 ICRU73QO : Emin= 0 eV Emax= 2 MeV 182 BetheBloch : Emin= 2 MeV Emax= 10 TeV 383 ICRU73QO : Emin= 0 eV Emax= 2 MeV 384 BetheBloch : Emin= 2 MeV Emax= 10 TeV 385 386 msc: for kaon+ SubType= 10 387 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 388 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 389 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 390 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 183 391 184 392 hIoni: for kaon+ SubType= 2 … … 187 395 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 188 396 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 189 Bragg : Emin= 0 eV Emax=1.05231 MeV190 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV397 Bragg : Emin= 0 eV Emax= 1.05231 MeV 398 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 191 399 192 400 hBrems: for kaon+ SubType= 3 … … 194 402 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 195 403 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 196 hBrem : Emin= 0 eV Emax=10 TeV404 hBrem : Emin= 0 eV Emax= 10 TeV 197 405 198 406 hPairProd: for kaon+ SubType= 4 … … 200 408 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 201 409 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 202 hPairProd : Emin= 0 eV Emax=10 TeV410 hPairProd : Emin= 0 eV Emax= 10 TeV 203 411 204 412 hIoni: for kaon- SubType= 2 … … 207 415 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 208 416 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 209 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV210 BetheBloch : Emin= 1.05231 MeV Emax=10 TeV417 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 418 BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV 211 419 212 420 hBrems: for kaon- SubType= 3 … … 214 422 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 215 423 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 216 hBrem : Emin= 0 eV Emax=10 TeV424 hBrem : Emin= 0 eV Emax= 10 TeV 217 425 218 426 hPairProd: for kaon- SubType= 4 … … 220 428 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 221 429 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 222 hPairProd : Emin= 0 eV Emax=10 TeV430 hPairProd : Emin= 0 eV Emax= 10 TeV 223 431 224 432 muMsc: for mu+ SubType= 10 225 433 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 226 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0227 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 228 UrbanMsc90 : Emin= 0 eV Emax=10 TeV434 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 11.4592 435 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 436 WentzelVIUni : Emin= 0 eV Emax= 10 TeV 229 437 230 438 muIoni: for mu+ SubType= 2 … … 233 441 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 234 442 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 235 Bragg : Emin= 0 eV Emax=200 keV236 BetheBloch : Emin= 200 keV Emax=1 GeV237 MuBetheBloch : Emin= 1 GeV Emax=10 TeV443 Bragg : Emin= 0 eV Emax= 200 keV 444 BetheBloch : Emin= 200 keV Emax= 1 GeV 445 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 238 446 ===== EM models for the G4Region GasDetector ====== 239 PAIModel : Emin= 0 eV Emax=10 TeV447 PAIModel : Emin= 0 eV Emax= 10 TeV 240 448 241 449 muBrems: for mu+ SubType= 3 … … 243 451 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 244 452 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 245 MuBrem : Emin= 0 eV Emax=10 TeV453 MuBrem : Emin= 0 eV Emax= 10 TeV 246 454 247 455 muPairProd: for mu+ SubType= 4 … … 249 457 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 250 458 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 251 muPairProd : Emin= 0 eV Emax= 10 TeV 459 muPairProd : Emin= 0 eV Emax= 10 TeV 460 461 CoulombScat: for mu+ SubType= 1 462 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 463 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 464 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 465 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 252 466 253 467 muIoni: for mu- SubType= 2 … … 256 470 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 257 471 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 258 ICRU73QO : Emin= 0 eV Emax=200 keV259 BetheBloch : Emin= 200 keV Emax=1 GeV260 MuBetheBloch : Emin= 1 GeV Emax=10 TeV472 ICRU73QO : Emin= 0 eV Emax= 200 keV 473 BetheBloch : Emin= 200 keV Emax= 1 GeV 474 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 261 475 ===== EM models for the G4Region GasDetector ====== 262 PAIModel : Emin= 0 eV Emax=10 TeV476 PAIModel : Emin= 0 eV Emax= 10 TeV 263 477 264 478 muBrems: for mu- SubType= 3 … … 266 480 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 267 481 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 268 MuBrem : Emin= 0 eV Emax=10 TeV482 MuBrem : Emin= 0 eV Emax= 10 TeV 269 483 270 484 muPairProd: for mu- SubType= 4 … … 272 486 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 273 487 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 274 muPairProd : Emin= 0 eV Emax= 10 TeV 488 muPairProd : Emin= 0 eV Emax= 10 TeV 489 490 CoulombScat: for mu- SubType= 1 491 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 492 11.4592 < Theta(degree) < 180, Eth(MeV)= 10; pLimit(GeV^1)= 0.139531 493 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 494 eCoulombScattering : Emin= 0 eV Emax= 10 TeV 275 495 276 496 hIoni: for pi+ SubType= 2 … … 279 499 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 280 500 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 281 Bragg : Emin= 0 eV Emax= 297.504keV282 BetheBloch : Emin= 297.504 keV Emax=10 TeV501 Bragg : Emin= 0 eV Emax= 297.505 keV 502 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 283 503 ===== EM models for the G4Region GasDetector ====== 284 PAIModel : Emin= 0 eV Emax=10 TeV504 PAIModel : Emin= 0 eV Emax= 10 TeV 285 505 286 506 hBrems: for pi+ SubType= 3 … … 288 508 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 289 509 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 290 hBrem : Emin= 0 eV Emax=10 TeV510 hBrem : Emin= 0 eV Emax= 10 TeV 291 511 292 512 hPairProd: for pi+ SubType= 4 … … 294 514 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 295 515 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 296 hPairProd : Emin= 0 eV Emax=10 TeV516 hPairProd : Emin= 0 eV Emax= 10 TeV 297 517 298 518 msc: for pi- SubType= 10 … … 300 520 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1 301 521 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 302 UrbanMsc90 : Emin= 0 eV Emax=10 TeV522 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 303 523 304 524 hIoni: for pi- SubType= 2 … … 307 527 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 308 528 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 309 ICRU73QO : Emin= 0 eV Emax= 297.504keV310 BetheBloch : Emin= 297.504 keV Emax=10 TeV529 ICRU73QO : Emin= 0 eV Emax= 297.505 keV 530 BetheBloch : Emin= 297.505 keV Emax= 10 TeV 311 531 ===== EM models for the G4Region GasDetector ====== 312 PAIModel : Emin= 0 eV Emax=10 TeV532 PAIModel : Emin= 0 eV Emax= 10 TeV 313 533 314 534 hBrems: for pi- SubType= 3 … … 316 536 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 317 537 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 318 hBrem : Emin= 0 eV Emax=10 TeV538 hBrem : Emin= 0 eV Emax= 10 TeV 319 539 320 540 hPairProd: for pi- SubType= 4 … … 322 542 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 323 543 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 324 hPairProd : Emin= 0 eV Emax=10 TeV544 hPairProd : Emin= 0 eV Emax= 10 TeV 325 545 326 546 Region <DefaultRegionForTheWorld> -- appears in <World> world volume 327 547 Root logical volume(s) : World 328 548 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] 329 Materials : Mylar549 Materials : G4_Galactic G4_MYLAR 330 550 Production cuts : gamma 500 um e- 500 um e+ 500 um proton 500 um 331 551 332 552 Region <GasDetector> -- appears in <World> world volume 333 Root logical volume(s) : Absorber553 Root logical volume(s) : Gas 334 554 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] 335 555 Materials : XeCH4C3H8 … … 339 559 340 560 Index : 0 used in the geometry : Yes recalculation needed : No 341 Material : Mylar561 Material : G4_Galactic 342 562 Range cuts : gamma 500 um e- 500 um e+ 500 um proton 500 um 343 Energy thresholds : gamma 2.27786 keV e- 261.805 keV e+ 258.31 keVproton 50 keV563 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 50 keV 344 564 Region(s) which use this couple : 345 565 DefaultRegionForTheWorld 346 566 347 567 Index : 1 used in the geometry : Yes recalculation needed : No 568 Material : G4_MYLAR 569 Range cuts : gamma 500 um e- 500 um e+ 500 um proton 500 um 570 Energy thresholds : gamma 2.27786 keV e- 263.57 keV e+ 258.31 keV proton 50 keV 571 Region(s) which use this couple : 572 DefaultRegionForTheWorld 573 574 Index : 2 used in the geometry : Yes recalculation needed : No 348 575 Material : XeCH4C3H8 349 576 Range cuts : gamma 1.8 mm e- 1.8 mm e+ 1.8 mm proton 1.8 mm … … 359 586 Total memory consumed for geometry optimisation: 0 kByte 360 587 Total CPU time elapsed for geometry optimisation: 0 seconds 361 ### Run 0 start. 362 363 --------- Ranecu engine status --------- 364 Initial seed (index) = 0 365 Current couple of seeds = 9876, 54321 366 ---------------------------------------- 588 ### Run 0 start 589 HistoManager: Histograms are booked and run has been started 367 590 Start Run processing. 368 369 ---> Begin of Event: 0 370 371 ---> Begin of Event: 100 372 373 ---> Begin of Event: 200 374 375 ---> Begin of Event: 300 376 377 ---> Begin of Event: 400 378 379 ---> Begin of Event: 500 380 381 ---> Begin of Event: 600 382 383 ---> Begin of Event: 700 384 385 ---> Begin of Event: 800 386 387 ---> Begin of Event: 900 591 ### PrimaryGeneratorAction::GeneratePrimaries ########## 592 ### proton E(MeV)= 200000 593 ######################################################## 388 594 Run terminated. 389 595 Run Summary 390 596 Number of events processed : 1000 391 User=0.11s Real=0.1s Sys=0s 597 User=0.11s Real=0.11s Sys=0s 598 RunAction: End of run actions are started 392 599 ================== run summary ===================== 393 end of Run TotNbofEvents = 1000 394 mean charged track length in absorber=24.2704 +- 0.310817 mm 395 396 mean energy deposit in absorber=0.0160189 +- 0.000378691 MeV 397 398 mean number of steps in absorber (charged) =0.898 +- 0.0848269 399 mean number of steps in absorber (neutral) =0.002 +- 0.0014128 400 401 mean number of charged secondaries = 0.257 +- 0.0163997 402 403 mean number of neutral secondaries = 0.002 +- 0.0014128 404 405 mean number of e-s =0.257 and e+s =0 406 407 (number) transmission coeff=1 reflection coeff=0 408 409 energy deposit distribution 410 #entries=1000 #underflows=0 #overflows=12 600 End of Run TotNbofEvents = 1000 601 602 Mean energy deposit in absorber = 15.134 +- 0.24587 keV; RMS/Emean = 0.51376 603 Mean number of steps in absorber= 1.253; mean number of ion-clusters = 794.89 604 605 ====== Energy deposit distribution Noverflows= 0.009 ====== 411 606 bin nb Elow entries normalized 412 607 0 0 0 0 413 1 3 4 0.004414 2 6 95 0.095415 3 9 27 0 0.27416 4 12 2 51 0.251417 5 15 14 6 0.146418 6 18 90 0.09419 7 21 49 0.049420 8 24 3 0 0.03421 9 27 16 0.016608 1 3 5 0.005 609 2 6 105 0.105 610 3 9 278 0.278 611 4 12 238 0.238 612 5 15 141 0.141 613 6 18 75 0.075 614 7 21 50 0.05 615 8 24 32 0.032 616 9 27 21 0.021 422 617 10 30 13 0.013 423 11 33 5 0.005424 12 36 5 0.005425 13 39 9 0.009426 14 42 1 0.001618 11 33 6 0.006 619 12 36 9 0.009 620 13 39 3 0.003 621 14 42 2 0.002 427 622 15 45 2 0.002 428 623 16 48 1 0.001 429 17 51 0 0 430 18 54 0 0 431 19 57 1 0.001 432 433 Emp = 0.009 width= 0.006 MeV 434 435 436 --------- Ranecu engine status --------- 437 Initial seed (index) = 0 438 Current couple of seeds = 1481926969, 852281434 439 ---------------------------------------- 624 17 51 2 0.002 625 18 54 5 0.005 626 19 57 3 0.003 627 ================== run end ========================== 440 628 # 441 629 UserDetectorConstruction deleted. -
trunk/examples/extended/electromagnetic/TestEm8/include/PhysicsList.hh
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsList.hh,v 1.5 2010/09/08 09:12:10 vnivanch Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 26 28 // 27 // $Id: PhysicsList.hh,v 1.4 2010/04/13 08:12:23 vnivanch Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 //--------------------------------------------------------------------------- 29 30 // 30 // ....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......31 // ClassName: PhysicsList 31 32 // 32 // 14.10.02 (V.Ivanchenko) provide modular list on base of old PhysicsList33 // Description: EM physics with a possibility to add PAI model 33 34 // 34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 35 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 35 // Author: V.Ivanchenko 01.09.2010 36 // 37 //---------------------------------------------------------------------------- 38 // 36 39 37 40 #ifndef PhysicsList_h … … 51 54 { 52 55 public: 56 53 57 PhysicsList(); 54 58 virtual ~PhysicsList(); … … 60 64 void SetCutForElectron(G4double); 61 65 void SetCutForPositron(G4double); 66 void SetCutForProton(G4double); 62 67 63 68 void AddPhysicsList(const G4String& name); 64 69 void ConstructProcess(); 65 70 66 void AddStepMax(); 67 //StepMax* GetStepMaxProcess() {return stepMaxProcess;}; 71 void AddStepMax(); 68 72 69 73 private: … … 78 82 G4double cutForElectron; 79 83 G4double cutForPositron; 84 G4double cutForProton; 80 85 81 86 G4VPhysicsConstructor* emPhysicsList; -
trunk/examples/extended/electromagnetic/TestEm8/include/PhysicsListMessenger.hh
r1337 r1342 25 25 // 26 26 // 27 // $Id: PhysicsListMessenger.hh,v 1. 2 2006/06/29 16:59:49 gunterExp $28 // GEANT4 tag $Name: geant4-09-04-beta-01$27 // $Id: PhysicsListMessenger.hh,v 1.3 2010/09/08 09:12:10 vnivanch Exp $ 28 // GEANT4 tag $Name: examples-V09-03-09 $ 29 29 // 30 //--------------------------------------------------------------------------- 31 // 32 // ClassName: PhysicsListMessenger 33 // 34 // Description: EM physics with a possibility to add PAI model 35 // 36 // Author: V.Ivanchenko 01.09.2010 37 // 38 //---------------------------------------------------------------------------- 30 39 // 31 40 … … 42 51 class G4UIcmdWithADoubleAndUnit; 43 52 class G4UIcmdWithAString; 53 class G4UIcmdWithAnInteger; 44 54 45 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 47 57 class PhysicsListMessenger: public G4UImessenger 48 58 { 49 59 public: 50 60 51 52 ~PhysicsListMessenger();61 PhysicsListMessenger(PhysicsList* ); 62 virtual ~PhysicsListMessenger(); 53 63 54 64 void SetNewValue(G4UIcommand*, G4String); 55 65 56 66 private: 57 67 58 68 PhysicsList* pPhysicsList; 59 69 60 G4UIcmdWithADoubleAndUnit* gammaCutCmd; 61 G4UIcmdWithADoubleAndUnit* electCutCmd; 62 G4UIcmdWithADoubleAndUnit* protoCutCmd; 63 G4UIcmdWithADoubleAndUnit* allCutCmd; 64 G4UIcmdWithAString* pListCmd; 70 G4UIcmdWithADoubleAndUnit* eCmd; 71 G4UIcmdWithAnInteger* ebCmd; 72 G4UIcmdWithAnInteger* cbCmd; 73 G4UIcmdWithAString* pListCmd; 65 74 66 75 }; -
trunk/examples/extended/electromagnetic/TestEm8/src/PhysicsList.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // 27 // $Id: PhysicsList.cc,v 1.19 2010/06/04 17:08:15 vnivanch Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 26 // $Id: PhysicsList.cc,v 1.21 2010/10/27 14:52:07 vnivanch Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 28 // 29 //--------------------------------------------------------------------------- 30 // 31 // ClassName: PhysicsList 32 // 33 // Description: EM physics with a possibility to add PAI model 34 // 35 // Author: V.Ivanchenko 01.09.2010 36 // 37 //---------------------------------------------------------------------------- 29 38 // 30 39 … … 39 48 #include "G4EmStandardPhysics_option2.hh" 40 49 #include "G4EmStandardPhysics_option3.hh" 50 #include "G4EmLivermorePhysics.hh" 51 #include "G4EmPenelopePhysics.hh" 41 52 #include "G4DecayPhysics.hh" 42 53 … … 47 58 #include "G4Electron.hh" 48 59 #include "G4Positron.hh" 60 #include "G4Proton.hh" 49 61 50 62 #include "G4UnitsTable.hh" … … 67 79 cutForElectron = defaultCutValue; 68 80 cutForPositron = defaultCutValue; 69 70 stepMaxProcess = 0; 81 cutForProton = defaultCutValue; 71 82 72 83 pMessenger = new PhysicsListMessenger(this); 84 85 stepMaxProcess = new StepMax(); 73 86 74 87 // Decay Physics is always defined … … 140 153 emPhysicsList = new G4EmStandardPhysics_option3(); 141 154 155 } else if (name == "emlivermore") { 156 157 emName = name; 158 delete emPhysicsList; 159 emPhysicsList = new G4EmLivermorePhysics(); 160 161 } else if (name == "empenelope") { 162 163 emName = name; 164 delete emPhysicsList; 165 emPhysicsList = new G4EmPenelopePhysics(); 166 142 167 } else if (name == "pai") { 143 168 … … 163 188 { 164 189 // Step limitation seen as a process 165 stepMaxProcess = new StepMax();166 190 167 191 theParticleIterator->reset(); … … 195 219 SetCutValue(cutForElectron, "e-"); 196 220 SetCutValue(cutForPositron, "e+"); 221 SetCutValue(cutForProton, "proton"); 197 222 198 223 if ( verboseLevel > 0 ) { DumpCutValuesTable(); } … … 221 246 cutForPositron = cut; 222 247 SetParticleCuts(cutForPositron, G4Positron::Positron()); 248 } 249 250 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 251 252 void PhysicsList::SetCutForProton(G4double cut) 253 { 254 cutForPositron = cut; 255 SetParticleCuts(cutForProton, G4Proton::Proton()); 223 256 } 224 257 -
trunk/examples/extended/electromagnetic/TestEm8/src/PhysicsListMessenger.cc
r1337 r1342 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsListMessenger.cc,v 1.3 2010/09/08 09:12:10 vnivanch Exp $ 27 // GEANT4 tag $Name: examples-V09-03-09 $ 26 28 // 27 // $Id: PhysicsListMessenger.cc,v 1.2 2006/06/29 17:00:47 gunter Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 //--------------------------------------------------------------------------- 29 30 // 31 // ClassName: PhysicsListMessenger 32 // 33 // Description: EM physics with a possibility to add PAI model 34 // 35 // Author: V.Ivanchenko 01.09.2010 36 // 37 //---------------------------------------------------------------------------- 30 38 // 31 39 … … 38 46 #include "G4UIcmdWithADoubleAndUnit.hh" 39 47 #include "G4UIcmdWithAString.hh" 48 #include "G4UIcmdWithAnInteger.hh" 49 #include "HistoManager.hh" 40 50 41 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 44 54 :pPhysicsList(pPhys) 45 55 { 46 gammaCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setGCut",this);47 gammaCutCmd->SetGuidance("Set gamma cut.");48 gammaCutCmd->SetParameterName("Gcut",false);49 gammaCutCmd->SetUnitCategory("Length");50 gammaCutCmd->SetRange("Gcut>0.0");51 gammaCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle);56 eCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setMaxE",this); 57 eCmd->SetGuidance("Set max energy deposit"); 58 eCmd->SetParameterName("Emax",false); 59 eCmd->SetUnitCategory("Energy"); 60 eCmd->SetRange("Emax>0.0"); 61 eCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 52 62 53 electCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setECut",this); 54 electCutCmd->SetGuidance("Set electron cut."); 55 electCutCmd->SetParameterName("Ecut",false); 56 electCutCmd->SetUnitCategory("Length"); 57 electCutCmd->SetRange("Ecut>0.0"); 58 electCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 63 ebCmd = new G4UIcmdWithAnInteger("/testem/phys/setNbinsE",this); 64 ebCmd->SetGuidance("Set number of bins in energy."); 65 ebCmd->SetParameterName("Ebins",false); 66 ebCmd->SetRange("Ebins>0"); 67 ebCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 59 68 60 protoCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setPCut",this); 61 protoCutCmd->SetGuidance("Set positron cut."); 62 protoCutCmd->SetParameterName("Pcut",false); 63 protoCutCmd->SetUnitCategory("Length"); 64 protoCutCmd->SetRange("Pcut>0.0"); 65 protoCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 66 67 allCutCmd = new G4UIcmdWithADoubleAndUnit("/testem/phys/setCuts",this); 68 allCutCmd->SetGuidance("Set cut for all."); 69 allCutCmd->SetParameterName("cut",false); 70 allCutCmd->SetUnitCategory("Length"); 71 allCutCmd->SetRange("cut>0.0"); 72 allCutCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 69 cbCmd = new G4UIcmdWithAnInteger("/testem/phys/setNbinsCl",this); 70 cbCmd->SetGuidance("Set max number of clusters."); 71 cbCmd->SetParameterName("Cbins",false); 72 cbCmd->SetRange("Cbins>0"); 73 cbCmd->AvailableForStates(G4State_PreInit,G4State_Idle); 73 74 74 75 pListCmd = new G4UIcmdWithAString("/testem/phys/addPhysics",this); … … 82 83 PhysicsListMessenger::~PhysicsListMessenger() 83 84 { 84 delete gammaCutCmd; 85 delete electCutCmd; 86 delete protoCutCmd; 87 delete allCutCmd; 85 delete eCmd; 86 delete ebCmd; 87 delete cbCmd; 88 88 delete pListCmd; 89 89 } … … 94 94 G4String newValue) 95 95 { 96 if( command == gammaCutCmd ) 97 { pPhysicsList->SetCutForGamma(gammaCutCmd->GetNewDoubleValue(newValue));} 96 HistoManager* man = HistoManager::GetPointer(); 97 if( command == eCmd ) 98 { man->SetMaxEnergy(eCmd->GetNewDoubleValue(newValue));} 98 99 99 if( command == e lectCutCmd )100 { pPhysicsList->SetCutForElectron(electCutCmd->GetNewDoubleValue(newValue));}100 if( command == ebCmd ) 101 { man->SetNumberBins(ebCmd->GetNewIntValue(newValue));} 101 102 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 } 103 if( command == cbCmd ) 104 { man->SetNumberBinsCluster(cbCmd->GetNewIntValue(newValue));} 112 105 113 106 if( command == pListCmd )
Note: See TracChangeset
for help on using the changeset viewer.
