Changeset 1337 for trunk/examples/extended/electromagnetic/TestEm11
- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- Location:
- trunk/examples/extended/electromagnetic/TestEm11
- Files:
-
- 41 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/electromagnetic/TestEm11/History
r807 r1337 1 $Id: History,v 1. 16 2007/11/07 17:22:15 maireExp $1 $Id: History,v 1.33 2010/06/07 05:40:45 perl Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 06-06-10 J.Perl (testem11-V09-03-03) 19 - Remove unused variable in EventAction 20 21 03-06-10 J.Perl (testem11-V09-03-02) 22 - Updated vis usage 23 24 25-05-10 mma (testem11-V09-03-01) 25 - add PhysListStandardGS, PhysListStandardSS, PhysListStandardWVI 26 27 20-05-10 mma (testem11-V09-03-00) 28 - TestEm11.cc : introduction of G4UIExecutive 29 30 16-11-09 mma (testem11-V09-02-01) 31 - PhysListEmStandard copied from TestEm5 32 33 22-09-09 mma (testem11-V09-02-00) 34 - PhysicsList.cc : 35 - renamed "standard" -> "local"; "emstandard" -> "emstandard_opt0" 36 - Substituted obsolete low-energy physics builders by builders from 37 the physics_list library; 38 - updated README and macro files 39 40 16-11-08 mma (testem11-V09-01-07) 41 - Em options: 20bins/decade LinLossLimit=default=0.01 42 43 30-10-08 mma (testem11-V09-01-06) 44 - update PhysListEmStandard : options and hMultipleScattering for muon 45 46 24-10-08 mma (testem11-V09-01-05) 47 - update PhysListEmStandard and README 48 49 22-09-08 mma (testem11-V09-01-04) 50 - HistoManager : root format by default 51 52 12-06-08 mma (testem11-V09-01-03) 53 - Remove AIDA from GNUmakefile 54 55 07-05-08 mma (testem11-V09-01-02) 56 - README : update Aida informations 57 58 11-04-08 mma 59 - HistoManager::Scale() : suppress warning 60 61 07-04-08 mma (testem11-V09-01-01) 62 - PhysListEmStandard : add G4eMultipleScattering; 63 - HistoManager::saveAscii() : redefinition of binEdge 64 65 10-01-08 j.jacquemier (testem11-V09-01-00) 66 - add setfileOption in HistoMessenger 67 68 06-12-07 mma 69 - new subdirectory sandia (from O. Kadri) 17 70 18 71 07-11-07 mma (testem11-V09-00-01) -
trunk/examples/extended/electromagnetic/TestEm11/README
r807 r1337 1 $Id: README,v 1. 7 2007/11/07 17:22:15maire Exp $1 $Id: README,v 1.11 2009/09/22 14:20:31 maire Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 31 31 2- PHYSICS LIST 32 32 33 The particle list is the one of novice/exampleN02. 34 The physics list contains the 'standard' electromagnetic processes, 35 and decay. 36 37 Few commands have been added to PhysicsList, in order to set the 38 production threshold for secondaries either in range for gamma, e-/e+. 33 Physics lists can be local (eg. in this example) or from G4 kernel 34 physics_lists subdirectory. 35 36 Local physics lists: 37 - "local" standard EM physics with current 'best' options setting. 38 these options are explicited in PhysListEmStandard 39 40 From geant4/source/physics_lists/builders: 41 - "emstandard_opt0" recommended standard EM physics for LHC 42 - "emstandard_opt1" best CPU performance standard physics for LHC 43 - "emstandard_opt2" 44 - "emstandard_opt3" best current advanced EM options. 45 analog to "local" above 46 - "emlivermore" low-energy EM physics using Livermore data 47 - "empenelope" low-energy EM physics implementing Penelope models 48 49 Physics lists and options can be (re)set with UI commands 50 51 Please, notice that options set through G4EmProcessOPtions are global, eg 52 for all particle types. In G4 builders, it is shown how to set options per 53 particle type. 54 39 55 40 56 3- AN EVENT : THE PRIMARY GENERATOR … … 128 144 129 145 It is possible to choose the format of the histogram file (hbook, root, XML) 130 with the command /testem/histo/setFileType ( hbookby default)146 with the command /testem/histo/setFileType (root by default) 131 147 132 148 It is also possible to print selected histograms on an ascii file: … … 141 157 142 158 To use histograms, at least one of the AIDA implementations should be 143 available (see http://aida.freehep.org). 144 145 8a - PI 146 147 A package including AIDA and extended interfaces also using Python is PI, 148 available from: http://cern.ch/pi 149 150 Once installed PI or PI-Lite in a specified local area $MYPY, it is required 151 to add the installation path to $PATH, i.e. for example, for release 1.2.1 of 152 PI: 153 setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin 154 155 CERN users can use the PATH to the LCG area on AFS. 156 Before running the example the command should be issued: 157 eval `aida-config --runtime csh` 158 159 8b - OpenScientist 160 161 OpenScientist is available at http://OpenScientist.lal.in2p3.fr. 162 163 You have to "setup" the OpenScientist AIDA implementation before compiling 164 (then with G4ANALYSIS_USE set) and running your Geant4 application. 165 166 On UNIX you setup, with a csh flavoured shell : 167 csh> source <<OpenScientist install path>/aida-setup.csh 168 or with a sh flavoured shell : 169 sh> . <<OpenScientist install path>/aida-setup.sh 170 On Windows : 171 DOS> call <<OpenScientist install path>/aida-setup.bat 172 173 You can use various file formats for writing (AIDA-XML, hbook, root). 174 These formats are readable by the Lab onx interactive program 175 or the OpenPAW application. See the web pages. 176 177 178 With OpenPAW, on a run.hbook file, one can view the histograms 179 with something like : 180 OS> opaw 181 opaw> h/file 1 run.hbook ( or opaw> h/file 1 run.aida or run.root) 182 opaw> zone 2 2 183 opaw> h/plot 1 184 opaw> h/plot 2 159 available. See InstallAida.txt -
trunk/examples/extended/electromagnetic/TestEm11/TestEm11.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: TestEm11.cc,v 1. 3 2006/06/29 16:39:18 gunterExp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: TestEm11.cc,v 1.4 2010/05/20 13:02:41 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 32 32 #include "G4RunManager.hh" 33 33 #include "G4UImanager.hh" 34 #include "G4UIterminal.hh"35 #include "G4UItcsh.hh"36 34 #include "Randomize.hh" 37 35 … … 49 47 #ifdef G4VIS_USE 50 48 #include "G4VisExecutive.hh" 49 #endif 50 51 #ifdef G4UI_USE 52 #include "G4UIExecutive.hh" 51 53 #endif 52 54 … … 71 73 runManager->SetUserInitialization(phys = new PhysicsList); 72 74 runManager->SetUserAction(kin = new PrimaryGeneratorAction(det)); 73 74 #ifdef G4VIS_USE75 //visualization manager76 G4VisManager* visManager = new G4VisExecutive;77 visManager->Initialize();78 #endif79 75 80 76 //set user action classes … … 89 85 90 86 //get the pointer to the User Interface manager 91 G4UImanager* UI = G4UImanager::GetUIpointer(); 92 93 if (argc==1) // Define UI terminal for interactive mode 94 { 95 G4UIsession* session = 0; 96 #ifdef G4UI_USE_TCSH 97 session = new G4UIterminal(new G4UItcsh); 98 #else 99 session = new G4UIterminal(); 100 #endif 101 session->SessionStart(); 102 delete session; 103 } 104 else // Batch mode 105 { 87 G4UImanager* UI = G4UImanager::GetUIpointer(); 88 89 if (argc!=1) // batch mode 90 { 106 91 G4String command = "/control/execute "; 107 92 G4String fileName = argv[1]; 108 93 UI->ApplyCommand(command+fileName); 109 94 } 95 96 else //define visualization and UI terminal for interactive mode 97 { 98 #ifdef G4VIS_USE 99 G4VisManager* visManager = new G4VisExecutive; 100 visManager->Initialize(); 101 #endif 102 103 #ifdef G4UI_USE 104 G4UIExecutive * ui = new G4UIExecutive(argc,argv); 105 ui->SessionStart(); 106 delete ui; 107 #endif 108 109 #ifdef G4VIS_USE 110 delete visManager; 111 #endif 112 } 110 113 111 114 //job termination 112 #ifdef G4VIS_USE 113 delete visManager; 114 #endif 115 115 // 116 116 delete histo; 117 117 delete runManager; -
trunk/examples/extended/electromagnetic/TestEm11/TestEm11.in
r807 r1337 1 # $Id: TestEm11.in,v 1. 6 2006/10/26 10:18:59maire Exp $1 # $Id: TestEm11.in,v 1.8 2009/09/22 14:20:31 maire Exp $ 2 2 # 3 3 # limit the step size from histo 1 … … 10 10 /testem/det/setSizeYZ 1 mm 11 11 # 12 /testem/phys/addPhysics standard # em physics 13 ###/testem/phys/addPhysics penelope # em physics 12 /testem/phys/addPhysics local # em physics 14 13 # 15 14 /run/initialize … … 18 17 /gun/energy 500 keV 19 18 # 20 /testem/histo/setFileType hbook19 /testem/histo/setFileType root 21 20 /testem/histo/setFileName testem11 22 21 /testem/histo/setHisto 1 100 0. 1. mm #edep profile -
trunk/examples/extended/electromagnetic/TestEm11/TestEm11.out
r807 r1337 1 2 ############################################ 3 !!! WARNING - FPE detection is activated !!! 4 ############################################ 1 5 2 6 ************************************************************* 3 Geant4 version Name: g lobal-V09-00-03 (9-May-2008)7 Geant4 version Name: geant4-09-03-ref-06 (25-June-2010) 4 8 Copyright : Geant4 Collaboration 5 9 Reference : NIM A 506 (2003), 250-303 … … 8 12 9 13 10 ***** Table : Nb of materials = 1 3*****11 12 Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Imean: 166.000 eV14 ***** Table : Nb of materials = 14 ***** 15 16 Material: G4_Al density: 2.699 g/cm3 RadL: 8.896 cm Nucl.Int.Length: 38.894 cm Imean: 166.000 eV 13 17 ---> Element: Al (Al) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 14 18 15 Material: G4_Si density: 2.330 g/cm3 RadL: 9.366 cm Imean: 173.000 eV19 Material: G4_Si density: 2.330 g/cm3 RadL: 9.366 cm Nucl.Int.Length: 45.660 cm Imean: 173.000 eV 16 20 ---> Element: Si (Si) Z = 14.0 N = 28.1 A = 28.09 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 17 21 18 Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Imean: 286.000 eV22 Material: G4_Fe density: 7.874 g/cm3 RadL: 1.757 cm Nucl.Int.Length: 16.990 cm Imean: 286.000 eV 19 23 ---> Element: Fe (Fe) Z = 26.0 N = 55.8 A = 55.85 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 20 24 21 Material: G4_Ge density: 5.323 g/cm3 RadL: 2.301 cm Imean: 350.000 eV 25 Material: G4_Cu density: 8.960 g/cm3 RadL: 1.436 cm Nucl.Int.Length: 15.588 cm Imean: 322.000 eV 26 ---> Element: Cu (Cu) Z = 29.0 N = 63.5 A = 63.55 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 27 28 Material: G4_Ge density: 5.323 g/cm3 RadL: 2.301 cm Nucl.Int.Length: 27.431 cm Imean: 350.000 eV 22 29 ---> Element: Ge (Ge) Z = 32.0 N = 72.6 A = 72.61 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 23 30 24 Material: G4_Mo density: 10.220 g/cm3 RadL: 9.591 mm Imean: 424.000 eV31 Material: G4_Mo density: 10.220 g/cm3 RadL: 9.591 mm Nucl.Int.Length: 15.677 cm Imean: 424.000 eV 25 32 ---> Element: Mo (Mo) Z = 42.0 N = 95.9 A = 95.93 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 26 33 27 Material: G4_Ta density: 16.654 g/cm3 RadL: 4.094 mm Imean: 718.000 eV34 Material: G4_Ta density: 16.654 g/cm3 RadL: 4.094 mm Nucl.Int.Length: 11.887 cm Imean: 718.000 eV 28 35 ---> Element: Ta (Ta) Z = 73.0 N = 180.9 A = 180.95 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 29 36 30 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Imean: 727.000 eV37 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.312 cm Imean: 727.000 eV 31 38 ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 32 39 33 Material: G4_Au density: 19.320 g/cm3 RadL: 3.344 mm Imean: 790.000 eV40 Material: G4_Au density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.540 cm Imean: 790.000 eV 34 41 ---> Element: Au (Au) Z = 79.0 N = 197.0 A = 196.97 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 35 42 36 Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Imean: 823.000 eV43 Material: G4_Pb density: 11.350 g/cm3 RadL: 5.613 mm Nucl.Int.Length: 18.248 cm Imean: 823.000 eV 37 44 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 38 45 39 Material: G4_PbWO4 density: 8.280 g/cm3 RadL: 8.925 mm Imean: 541.543eV46 Material: G4_PbWO4 density: 8.280 g/cm3 RadL: 8.925 mm Nucl.Int.Length: 20.740 cm Imean: 542.741 eV 40 47 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 14.06 % ElmAbundance 66.67 % 41 48 ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ElmMassFraction: 45.54 % ElmAbundance 16.67 % 42 49 ---> Element: W (W) Z = 74.0 N = 183.8 A = 183.84 g/mole ElmMassFraction: 40.40 % ElmAbundance 16.67 % 43 50 44 Material: G4_SODIUM_IODIDE density: 3.667 g/cm3 RadL: 2.588 cm Imean: 452.000 eV51 Material: G4_SODIUM_IODIDE density: 3.667 g/cm3 RadL: 2.588 cm Nucl.Int.Length: 42.914 cm Imean: 452.000 eV 45 52 ---> Element: Na (Na) Z = 11.0 N = 23.0 A = 22.99 g/mole ElmMassFraction: 15.34 % ElmAbundance 50.00 % 46 53 ---> Element: I (I) Z = 53.0 N = 126.9 A = 126.90 g/mole ElmMassFraction: 84.66 % ElmAbundance 50.00 % 47 54 48 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm55 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.095 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm 49 56 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 0.01 % ElmAbundance 0.02 % 50 57 ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole ElmMassFraction: 75.53 % ElmAbundance 78.44 % … … 52 59 ---> Element: Ar (Ar) Z = 18.0 N = 39.9 A = 39.95 g/mole ElmMassFraction: 1.28 % ElmAbundance 0.47 % 53 60 54 Material: Water density: 1.000 g/cm3 RadL: 36.083 cm Imean: 75.000 eV61 Material: Water density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.423 cm Imean: 75.000 eV 55 62 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.19 % ElmAbundance 66.67 % 56 63 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.81 % ElmAbundance 33.33 % 57 64 58 59 Visualization Manager instantiating...60 Visualization Manager initialising...61 Registering graphics systems...62 63 You have successfully registered the following graphics systems.64 Current available graphics systems are:65 ASCIITree (ATree)66 DAWNFILE (DAWNFILE)67 G4HepRep (HepRepXML)68 G4HepRepFile (HepRepFile)69 RayTracer (RayTracer)70 VRML1FILE (VRML1FILE)71 VRML2FILE (VRML2FILE)72 OpenGLImmediateX (OGLIX)73 OpenGLStoredX (OGLSX)74 75 Registering model factories...76 77 You have successfully registered the following model factories.78 Registered model factories:79 generic80 drawByCharge81 drawByParticleID82 drawByOriginVolume83 drawByAttribute84 85 Registered filter factories:86 chargeFilter87 particleFilter88 originVolumeFilter89 attributeFilter90 65 91 66 /run/verbose 2 … … 95 70 /testem/det/setSizeYZ 1 mm 96 71 # 97 /testem/phys/addPhysics standard 98 PhysicsList::AddPhysicsList: <standard> 99 ###/testem/phys/addPhysics penelope # em physics 72 /testem/phys/addPhysics local 73 PhysicsList::AddPhysicsList: <local> 100 74 # 101 75 /run/initialize … … 109 83 physicsList->Construct() start. 110 84 physicsList->Construct() start. 85 physicsList->CheckParticleList() start. 111 86 physicsList->setCut() start. 112 87 PhysicsList::SetCuts:CutLength : 1 mm … … 115 90 /gun/energy 500 keV 116 91 # 117 /testem/histo/setFileType hbook92 /testem/histo/setFileType root 118 93 /testem/histo/setFileName testem11 119 /testem/histo/setHisto 1 100 0. 1. mm 94 /testem/histo/setHisto 1 100 0. 1. mm 120 95 ----> SetHisto 1: Edep (MeV/mm) along absorber (mm); 100 bins from 0 mm to 1 mm 121 /testem/histo/setHisto 2 102 0. 510. keV 96 /testem/histo/setHisto 2 102 0. 510. keV 122 97 ----> SetHisto 2: total Energy deposited in absorber (keV); 102 bins from 0 keV to 510 keV 123 /testem/histo/setHisto 3 100 0. 2. mm 98 /testem/histo/setHisto 3 100 0. 2. mm 124 99 ----> SetHisto 3: true track length of the primary particle (mm); 100 bins from 0 mm to 2 mm 125 /testem/histo/setHisto 4 100 0. 200. um 100 /testem/histo/setHisto 4 100 0. 200. um 126 101 ----> SetHisto 4: true step size of the primary particle (um); 100 bins from 0 um to 200 um 127 /testem/histo/setHisto 5 100 0. 1. mm 102 /testem/histo/setHisto 5 100 0. 1. mm 128 103 ----> SetHisto 5: projected range of the primary particle (mm); 100 bins from 0 mm to 1 mm 129 104 # … … 132 107 /run/beamOn 20000 133 108 134 phot: Total cross sections from Sandia parametrisation. 135 Sampling according PhotoElectric model 136 137 compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV 138 Sampling according Klein-Nishina model 139 tables are built for gamma 140 Lambda tables from 100 eV to 100 GeV in 90 bins. 141 142 conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; 143 sampling secondary e+e- according Bethe-Heitler model 144 tables are built for gamma 145 Lambda tables from 1.022 MeV to 100 GeV in 100 bins. 146 147 msc: Model variant of multiple scattering for e- 148 Lambda tables from 100 eV to 100 TeV in 120 bins. 149 LateralDisplacementFlag= 1 Skin= 0 150 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 151 152 eIoni: tables are built for e- 153 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 154 Lambda tables from threshold to 100 TeV in 120 bins. 155 Delta cross sections and sampling from MollerBhabha model 156 Good description from 1 KeV to 100 GeV. 157 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 158 CSDA range table up to 1 GeV in 70 bins. 159 160 eBrem: tables are built for e- 161 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 162 Lambda tables from threshold to 100 TeV in 120 bins. 163 Total cross sections and sampling from StandBrem model (based on the EEDL data library) 164 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 165 166 eIoni: tables are built for e+ 167 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 168 Lambda tables from threshold to 100 TeV in 120 bins. 169 Delta cross sections and sampling from MollerBhabha model 170 Good description from 1 KeV to 100 GeV. 171 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 172 CSDA range table up to 1 GeV in 70 bins. 173 174 eBrem: tables are built for e+ 175 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 176 Lambda tables from threshold to 100 TeV in 120 bins. 177 Total cross sections and sampling from StandBrem model (based on the EEDL data library) 178 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. LPM flag 1 179 180 annihil: Sampling according eplus2gg model 181 tables are built for e+ 182 Lambda tables from 100 eV to 100 TeV in 120 bins. 183 184 msc: Model variant of multiple scattering for proton 185 Lambda tables from 100 eV to 100 TeV in 120 bins. 186 LateralDisplacementFlag= 1 Skin= 0 187 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 188 189 hIoni: tables are built for proton 190 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 191 Lambda tables from threshold to 100 TeV in 120 bins. 192 Scaling relation is used from proton dE/dx and range. 193 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV 194 Parametrisation from Bragg for protons below. NuclearStopping= 1 195 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 196 CSDA range table up to 1 GeV in 70 bins. 197 198 msc: Model variant of multiple scattering for GenericIon 199 LateralDisplacementFlag= 0 Skin= 0 200 Boundary/stepping algorithm is active with RangeFactor= 0.2 Step limit type 1 201 202 ionIoni: tables are built for GenericIon 203 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 204 Lambda tables from threshold to 100 TeV in 120 bins. 205 Scaling relation is used from proton dE/dx and range. 206 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV 207 Parametrisation from BraggIon for protons below. NuclearStopping= 1 208 209 Stopping Power data for 8 ion/material pairs are used. 210 Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1, fluct: 1 211 CSDA range table up to 1 GeV in 70 bins. 212 213 hIoni: tables are built for anti_proton 214 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 215 Lambda tables from threshold to 100 TeV in 120 bins. 216 Scaling relation is used from proton dE/dx and range. 217 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV 218 Parametrisation from Bragg for protons below. NuclearStopping= 1 219 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 220 CSDA range table up to 1 GeV in 70 bins. 221 222 msc: Model variant of multiple scattering for mu+ 223 Lambda tables from 100 eV to 100 TeV in 120 bins. 224 LateralDisplacementFlag= 1 Skin= 0 225 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 226 227 muIoni: tables are built for mu+ 228 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 229 Lambda tables from threshold to 100 TeV in 120 bins. 230 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, 231 radiative corrections for E > 1 GeV 232 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 233 CSDA range table up to 1 GeV in 70 bins. 234 235 muBrems: tables are built for mu+ 236 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 237 Lambda tables from threshold to 100 TeV in 120 bins. 238 Parametrised model 239 240 muPairProd: tables are built for mu+ 241 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 242 Lambda tables from threshold to 100 TeV in 120 bins. 243 Parametrised model 244 245 muIoni: tables are built for mu- 246 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 247 Lambda tables from threshold to 100 TeV in 120 bins. 248 Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, 249 radiative corrections for E > 1 GeV 250 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 251 CSDA range table up to 1 GeV in 70 bins. 252 253 muBrems: tables are built for mu- 254 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 255 Lambda tables from threshold to 100 TeV in 120 bins. 256 Parametrised model 257 258 muPairProd: tables are built for mu- 259 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 260 Lambda tables from threshold to 100 TeV in 120 bins. 261 Parametrised model 262 263 hIoni: tables are built for pi+ 264 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 265 Lambda tables from threshold to 100 TeV in 120 bins. 266 Scaling relation is used from proton dE/dx and range. 267 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV 268 Parametrisation from Bragg for protons below. NuclearStopping= 1 269 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 270 CSDA range table up to 1 GeV in 70 bins. 271 272 msc: Model variant of multiple scattering for pi- 273 Lambda tables from 100 eV to 100 TeV in 120 bins. 274 LateralDisplacementFlag= 1 Skin= 0 275 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 1 276 277 hIoni: tables are built for pi- 278 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 279 Lambda tables from threshold to 100 TeV in 120 bins. 280 Scaling relation is used from proton dE/dx and range. 281 Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV 282 Parametrisation from Bragg for protons below. NuclearStopping= 1 283 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1 284 CSDA range table up to 1 GeV in 70 bins. 109 phot: for gamma SubType= 12 110 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 111 PhotoElectric : Emin= 0 eV Emax= 100 TeV 112 113 compt: for gamma SubType= 13 114 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 115 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 116 Klein-Nishina : Emin= 0 eV Emax= 100 TeV 117 118 conv: for gamma SubType= 14 119 Lambda tables from 1.022 MeV to 100 TeV in 240 bins, spline: 1 120 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 121 Bethe-Heitler : Emin= 0 eV Emax= 100 TeV 122 123 msc: for e- SubType= 10 124 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 125 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 126 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 127 UrbanMsc93 : Emin= 0 eV Emax= 100 TeV 128 129 eIoni: for e- SubType= 2 130 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 131 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 132 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 133 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 134 MollerBhabha : Emin= 0 eV Emax= 100 TeV 135 CSDA range table up to 1 GeV in 160 bins 136 137 eBrem: for e- SubType= 3 138 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 139 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 140 LPM flag: 1 for E > 1 GeV 141 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 142 eBrem : Emin= 0 eV Emax= 1 GeV 143 eBremRel : Emin= 1 GeV Emax= 100 TeV 144 145 eIoni: for e+ SubType= 2 146 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 147 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 148 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 149 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 150 MollerBhabha : Emin= 0 eV Emax= 100 TeV 151 CSDA range table up to 1 GeV in 160 bins 152 153 eBrem: for e+ SubType= 3 154 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 155 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 156 LPM flag: 1 for E > 1 GeV 157 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 158 eBrem : Emin= 0 eV Emax= 1 GeV 159 eBremRel : Emin= 1 GeV Emax= 100 TeV 160 161 annihil: for e+ SubType= 5 162 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 163 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 164 eplus2gg : Emin= 0 eV Emax= 100 TeV 165 166 msc: for proton SubType= 10 167 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 168 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 169 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 170 UrbanMsc90 : Emin= 0 eV Emax= 100 TeV 171 172 hIoni: for proton SubType= 2 173 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 174 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 175 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 176 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 177 Bragg : Emin= 0 eV Emax= 2 MeV 178 BetheBloch : Emin= 2 MeV Emax= 100 TeV 179 CSDA range table up to 1 GeV in 160 bins 180 181 hBrems: for proton SubType= 3 182 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 183 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 184 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 185 hBrem : Emin= 0 eV Emax= 100 TeV 186 187 hPairProd: for proton SubType= 4 188 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 189 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 190 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 191 hPairProd : Emin= 0 eV Emax= 100 TeV 192 193 msc: for GenericIon SubType= 10 194 RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0 195 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 196 UrbanMsc90 : Emin= 0 eV Emax= 100 TeV 197 198 ionIoni: for GenericIon SubType= 2 199 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 200 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 201 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 202 Stopping Power data for 17 ion/material pairs, nuclearStopping: 1 203 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 204 ParamICRU73 : Emin= 0 eV Emax= 100 TeV 205 CSDA range table up to 1 GeV in 160 bins 206 207 nuclearStopping: for GenericIon SubType= 1 208 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 209 ICRU49NucStopping : Emin= 0 eV Emax= 100 TeV 210 211 nuclearStopping: for alpha SubType= 1 212 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 213 ICRU49NucStopping : Emin= 0 eV Emax= 100 TeV 214 215 hIoni: for anti_proton SubType= 2 216 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 217 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 218 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 219 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 220 ICRU73QO : Emin= 0 eV Emax= 2 MeV 221 BetheBloch : Emin= 2 MeV Emax= 100 TeV 222 CSDA range table up to 1 GeV in 160 bins 223 224 hIoni: for kaon+ SubType= 2 225 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 226 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 227 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 228 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 229 Bragg : Emin= 0 eV Emax= 1.05231 MeV 230 BetheBloch : Emin= 1.05231 MeV Emax= 100 TeV 231 CSDA range table up to 1 GeV in 160 bins 232 233 hIoni: for kaon- SubType= 2 234 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 235 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 236 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 237 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 238 ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV 239 BetheBloch : Emin= 1.05231 MeV Emax= 100 TeV 240 CSDA range table up to 1 GeV in 160 bins 241 242 muMsc: for mu+ SubType= 10 243 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 244 RangeFactor= 0.04, step limit type: 2, lateralDisplacement: 1, polarAngleLimit(deg)= 0 245 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 246 UrbanMsc90 : Emin= 0 eV Emax= 100 TeV 247 248 muIoni: for mu+ SubType= 2 249 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 250 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 251 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 252 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 253 Bragg : Emin= 0 eV Emax= 200 keV 254 BetheBloch : Emin= 200 keV Emax= 1 GeV 255 MuBetheBloch : Emin= 1 GeV Emax= 100 TeV 256 CSDA range table up to 1 GeV in 160 bins 257 258 muBrems: for mu+ SubType= 3 259 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 260 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 261 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 262 MuBrem : Emin= 0 eV Emax= 100 TeV 263 264 muPairProd: for mu+ SubType= 4 265 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 266 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 267 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 268 muPairProd : Emin= 0 eV Emax= 100 TeV 269 270 muIoni: for mu- SubType= 2 271 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 272 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 273 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 274 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 275 ICRU73QO : Emin= 0 eV Emax= 200 keV 276 BetheBloch : Emin= 200 keV Emax= 1 GeV 277 MuBetheBloch : Emin= 1 GeV Emax= 100 TeV 278 CSDA range table up to 1 GeV in 160 bins 279 280 muBrems: for mu- SubType= 3 281 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 282 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 283 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 284 MuBrem : Emin= 0 eV Emax= 100 TeV 285 286 muPairProd: for mu- SubType= 4 287 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 288 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 289 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 290 muPairProd : Emin= 0 eV Emax= 100 TeV 291 292 hIoni: for pi+ SubType= 2 293 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 294 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 295 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 296 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 297 Bragg : Emin= 0 eV Emax= 297.504 keV 298 BetheBloch : Emin= 297.504 keV Emax= 100 TeV 299 CSDA range table up to 1 GeV in 160 bins 300 301 hBrems: for pi+ SubType= 3 302 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 303 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 304 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 305 hBrem : Emin= 0 eV Emax= 100 TeV 306 307 hPairProd: for pi+ SubType= 4 308 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 309 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 310 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 311 hPairProd : Emin= 0 eV Emax= 100 TeV 312 313 msc: for pi- SubType= 10 314 Lambda tables from 100 eV to 100 TeV in 240 bins, spline: 1 315 RangeFactor= 0.04, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 316 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 317 UrbanMsc90 : Emin= 0 eV Emax= 100 TeV 318 319 hIoni: for pi- SubType= 2 320 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 321 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 322 finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 323 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 324 ICRU73QO : Emin= 0 eV Emax= 297.504 keV 325 BetheBloch : Emin= 297.504 keV Emax= 100 TeV 326 CSDA range table up to 1 GeV in 160 bins 327 328 hBrems: for pi- SubType= 3 329 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 330 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 331 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 332 hBrem : Emin= 0 eV Emax= 100 TeV 333 334 hPairProd: for pi- SubType= 4 335 dE/dx and range tables from 100 eV to 100 TeV in 240 bins 336 Lambda tables from threshold to 100 TeV in 240 bins, spline: 1 337 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 338 hPairProd : Emin= 0 eV Emax= 100 TeV 285 339 286 340 Region <DefaultRegionForTheWorld> -- appears in <Absorber> world volume … … 288 342 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] 289 343 Materials : G4_Si 290 Production cuts : gamma 1 mm e- 1 mm e+ 1 mm344 Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 291 345 292 346 ========= Table of registered couples ============================== … … 294 348 Index : 0 used in the geometry : Yes recalculation needed : No 295 349 Material : G4_Si 296 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm297 Energy thresholds : gamma 6.88731 keV e- 540.718 keV e+ 521.113 keV350 Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm 351 Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 0 eV 298 352 Region(s) which use this couple : 299 353 DefaultRegionForTheWorld … … 322 376 Run Summary 323 377 Number of events processed : 20000 324 User= 22.26s Real=71.63s Sys=5.33s378 User=12.36s Real=13.63s Sys=0.52s 325 379 326 380 ======================== run summary ====================== … … 330 384 ============================================================ 331 385 332 Total Energy deposited = 44 9.403 keV +- 103.444keV333 334 Track length of primary track = 87 5.191 um +- 189.287um335 Range from EmCalculator = 9 37.365um (from full dE/dx)336 337 Projected range = 3 36.299 um +- 207.386um338 339 Nb of steps of primary track = 8 8.04 +- 18.91 Step size= 9.933 um +- 86.256nm340 341 absorbed = 7 8.50 % transmit = 5.42 % reflected = 16.08%386 Total Energy deposited = 443.666 keV +- 110.109 keV 387 388 Track length of primary track = 870.965 um +- 199.370 um 389 Range from EmCalculator = 943.428 um (from full dE/dx) 390 391 Projected range = 321.497 um +- 208.092 um 392 393 Nb of steps of primary track = 87.62 +- 19.93 Step size= 9.929 um +- 123.724 nm 394 395 absorbed = 76.69 % transmit = 5.10 % reflected = 18.20 % 342 396 343 397 --------- Ranecu engine status --------- 344 398 Initial seed (index) = 0 345 Current couple of seeds = 1 622151268, 626459451399 Current couple of seeds = 1325532186, 417920133 346 400 ---------------------------------------- 347 Graphics systems deleted.348 Visualization Manager deleting...349 401 UserDetectorConstruction deleted. 350 402 UserPhysicsList deleted. … … 353 405 G4 kernel has come to Quit state. 354 406 EventManager deleted. 355 Default detector region deleted.356 407 UImanager deleted. 357 408 Units table cleared. … … 359 410 RunManagerKernel is deleted. 360 411 RunManager is deleting. 412 -
trunk/examples/extended/electromagnetic/TestEm11/include/DetectorConstruction.hh
r1230 r1337 25 25 // 26 26 // $Id: DetectorConstruction.hh,v 1.2 2006/06/29 16:39:20 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/DetectorMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: DetectorMessenger.hh,v 1.2 2006/06/29 16:39:22 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/EventAction.hh
r1230 r1337 25 25 // 26 26 // $Id: EventAction.hh,v 1.2 2006/06/29 16:39:25 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/EventActionMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.hh,v 1.2 2006/06/29 16:39:27 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/HistoManager.hh
r1230 r1337 25 25 // 26 26 // $Id: HistoManager.hh,v 1.5 2007/11/07 17:22:15 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/HistoMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: HistoMessenger.hh,v 1.5 2008/01/11 13:01:57 jjacquem Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/PhysListEmStandard.hh
r1230 r1337 26 26 // 27 27 // $Id: PhysListEmStandard.hh,v 1.2 2006/06/29 16:39:41 gunter Exp $ 28 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 29 // 30 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/PhysicsList.hh
r1230 r1337 26 26 // 27 27 // $Id: PhysicsList.hh,v 1.3 2006/06/29 16:39:43 gunter Exp $ 28 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 29 // 30 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/PhysicsListMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.hh,v 1.2 2006/06/29 16:39:46 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/PrimaryGeneratorAction.hh
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorAction.hh,v 1.2 2006/06/29 16:39:49 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/PrimaryGeneratorMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorMessenger.hh,v 1.2 2006/06/29 16:39:52 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/RunAction.hh
r1230 r1337 25 25 // 26 26 // $Id: RunAction.hh,v 1.3 2007/08/19 20:52:53 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/StepMax.hh
r1230 r1337 25 25 // 26 26 // $Id: StepMax.hh,v 1.3 2006/06/29 16:39:57 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/StepMaxMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: StepMaxMessenger.hh,v 1.3 2006/06/29 16:40:00 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/SteppingAction.hh
r1230 r1337 25 25 // 26 26 // $Id: SteppingAction.hh,v 1.2 2006/06/29 16:40:03 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/SteppingVerbose.hh
r1230 r1337 25 25 // 26 26 // $Id: SteppingVerbose.hh,v 1.2 2006/06/29 16:40:06 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/include/TrackingAction.hh
r1230 r1337 25 25 // 26 26 // $Id: TrackingAction.hh,v 1.2 2006/06/29 16:40:08 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/sandia.mac
r1230 r1337 1 # $Id: sandia.mac,v 1. 9 2009/09/22 14:20:31 maire Exp $1 # $Id: sandia.mac,v 1.10 2010/05/20 13:02:41 maire Exp $ 2 2 # 3 3 # limit the step size from histo 8 … … 16 16 /testem/det/setSizeYZ 1 cm 17 17 # 18 /testem/phys/addPhysics local # em physics 18 /testem/phys/addPhysics emstandard_opt3 # em physics 19 ###/testem/phys/addPhysics standardSS # em physics 19 20 # 20 21 /run/initialize 21 22 # 23 /run/setCut 1 mm # will set cut for proton 24 # 22 25 /gun/particle e- 23 /gun/energy 314keV26 /gun/energy 521 keV 24 27 # 25 28 ###/testem/histo/setFileType hbook 26 29 /testem/histo/setFileType root 27 /testem/histo/setFileName sandia/ Al.314keV30 /testem/histo/setFileName sandia/ref4++/Al.521keV.opt3 28 31 /testem/histo/setHisto 8 100 0. 1.0 none #normalized Edep profile 29 32 ###/testem/histo/setHisto 8 100 0. 0.5 none #for Ta -
trunk/examples/extended/electromagnetic/TestEm11/src/DetectorConstruction.cc
r1230 r1337 25 25 // 26 26 // $Id: DetectorConstruction.cc,v 1.9 2009/11/16 15:09:10 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 29 29 // -
trunk/examples/extended/electromagnetic/TestEm11/src/DetectorMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: DetectorMessenger.cc,v 1.2 2006/06/29 16:40:12 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/EventAction.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: EventAction.cc,v 1. 2 2006/06/29 16:40:15 gunterExp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: EventAction.cc,v 1.4 2010/06/07 05:40:45 perl Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 37 37 38 38 #include "G4Event.hh" 39 #include "G4TrajectoryContainer.hh"40 #include "G4Trajectory.hh"41 #include "G4VVisManager.hh"42 39 43 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 72 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 73 70 74 void EventAction::EndOfEventAction(const G4Event* evt)71 void EventAction::EndOfEventAction(const G4Event*) 75 72 { 76 73 //plot energy deposited per event … … 80 77 histoManager->FillHisto(2,totalEdep); 81 78 } 82 83 //trajectories84 //85 if (G4VVisManager::GetConcreteInstance())86 {87 G4TrajectoryContainer* trajectoryContainer = evt->GetTrajectoryContainer();88 G4int n_trajectories = 0;89 if (trajectoryContainer) n_trajectories = trajectoryContainer->entries();90 for(G4int i=0; i<n_trajectories; i++)91 { G4Trajectory* trj = (G4Trajectory*)92 ((*(evt->GetTrajectoryContainer()))[i]);93 if (drawFlag == "all") trj->DrawTrajectory(1000);94 else if ((drawFlag == "charged")&&(trj->GetCharge() != 0.))95 trj->DrawTrajectory(1000);96 }97 }98 79 } 99 80 -
trunk/examples/extended/electromagnetic/TestEm11/src/EventActionMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.cc,v 1.2 2006/06/29 16:40:18 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/HistoManager.cc
r1230 r1337 25 25 // 26 26 // $Id: HistoManager.cc,v 1.10 2008/09/22 15:11:34 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/HistoMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: HistoMessenger.cc,v 1.6 2008/01/11 13:01:57 jjacquem Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/PhysListEmStandard.cc
r1230 r1337 25 25 // 26 26 // $Id: PhysListEmStandard.cc,v 1.11 2009/11/16 15:09:10 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/PhysicsList.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsList.cc,v 1. 7 2009/09/22 14:20:31 maire Exp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: PhysicsList.cc,v 1.8 2010/05/20 13:02:41 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 34 34 35 35 #include "PhysListEmStandard.hh" 36 #include "PhysListEmStandardSS.hh" 37 #include "PhysListEmStandardGS.hh" 38 #include "PhysListEmStandardWVI.hh" 36 39 37 40 #include "G4EmStandardPhysics.hh" … … 241 244 delete emPhysicsList; 242 245 emPhysicsList = new G4EmStandardPhysics_option3(); 243 246 247 } else if (name == "standardSS") { 248 249 emName = name; 250 delete emPhysicsList; 251 emPhysicsList = new PhysListEmStandardSS(name); 252 253 } else if (name == "standardGS") { 254 255 emName = name; 256 delete emPhysicsList; 257 emPhysicsList = new PhysListEmStandardGS(name); 258 259 } else if (name == "standardWVI") { 260 261 emName = name; 262 delete emPhysicsList; 263 emPhysicsList = new PhysListEmStandardWVI(name); 264 244 265 } else if (name == "emlivermore") { 245 266 emName = name; -
trunk/examples/extended/electromagnetic/TestEm11/src/PhysicsListMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.cc,v 1.2 2006/06/29 16:40:38 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/PrimaryGeneratorAction.cc
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorAction.cc,v 1.2 2006/06/29 16:40:40 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/PrimaryGeneratorMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorMessenger.cc,v 1.2 2006/06/29 16:40:43 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/RunAction.cc
r1230 r1337 25 25 // 26 26 // $Id: RunAction.cc,v 1.7 2007/08/19 20:52:53 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/StepMax.cc
r1230 r1337 25 25 // 26 26 // $Id: StepMax.cc,v 1.3 2006/06/29 16:40:47 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/StepMaxMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: StepMaxMessenger.cc,v 1.3 2006/06/29 16:40:53 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/SteppingAction.cc
r1230 r1337 25 25 // 26 26 // $Id: SteppingAction.cc,v 1.5 2007/08/19 20:52:53 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/SteppingVerbose.cc
r1230 r1337 25 25 // 26 26 // $Id: SteppingVerbose.cc,v 1.2 2006/06/29 16:41:54 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/src/TrackingAction.cc
r1230 r1337 25 25 // 26 26 // $Id: TrackingAction.cc,v 1.2 2006/06/29 16:41:56 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/electromagnetic/TestEm11/vis.mac
r807 r1337 10 10 /run/initialize 11 11 # 12 # Create empty scene ("world" is default)13 /vis/ scene/create12 # Use this open statement to create an OpenGL view: 13 /vis/open OGL 600x600-0+0 14 14 # 15 # Create a scene handler for a specific graphics system 16 # Edit the next line(s) to choose another graphic system 17 # 18 /vis/open OGLIX 19 # 15 # Use this open statement to create a .prim file suitable for 16 # viewing in DAWN: 20 17 #/vis/open DAWNFILE 21 18 # 22 # Draw scene 19 # Use this open statement to create a .heprep file suitable for 20 # viewing in HepRApp: 21 #/vis/open HepRepFile 22 # 23 # Use this open statement to create a .wrl file suitable for 24 # viewing in a VRML viewer: 25 #/vis/open VRML2FILE 26 # 27 # Disable auto refresh and quieten vis messages whilst scene and 28 # trajectories are established: 29 /vis/viewer/set/autoRefresh false 30 /vis/verbose errors 31 # 32 # Draw geometry: 33 /vis/drawVolume 34 # 35 # Specify view angle: 36 #/vis/viewer/set/viewpointThetaPhi 90. 0. 37 # 38 # Specify zoom value: 23 39 /vis/viewer/zoom 1.4 24 /vis/viewer/flush25 40 # 26 # for drawing the tracks 27 # if too many tracks cause core dump => storeTrajectory 0 28 /tracking/storeTrajectory 1 41 # Specify style (surface or wireframe): 42 #/vis/viewer/set/style wireframe 43 # 44 # Draw coordinate axes: 45 #/vis/scene/add/axes 0 0 0 1 m 46 # 47 # Draw smooth trajectories at end of event, showing trajectory points 48 # as markers 2 pixels wide: 49 /vis/scene/add/trajectories smooth 50 /vis/modeling/trajectories/create/drawByCharge 51 /vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true 52 /vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2 53 # (if too many tracks cause core dump => /tracking/storeTrajectory 0) 54 # 55 # Draw hits at end of event: 56 #/vis/scene/add/hits 57 # 58 # To draw only gammas: 59 #/vis/filtering/trajectories/create/particleFilter 60 #/vis/filtering/trajectories/particleFilter-0/add gamma 61 # 62 # To invert the above, drawing all particles except gammas, 63 # keep the above two lines but also add: 64 #/vis/filtering/trajectories/particleFilter-0/invert true 65 # 66 # Many other options are available with /vis/modeling and /vis/filtering. 67 # For example, to select colour by particle ID: 68 #/vis/modeling/trajectories/create/drawByParticleID 69 #/vis/modeling/trajectories/drawByParticleID-0/set e- blue 70 # 71 # To superimpose all of the events from a given run: 29 72 /vis/scene/endOfEventAction accumulate 30 73 # 31 /testem/event/drawTracks all 32 74 # Re-establish auto refreshing and verbosity: 75 /vis/viewer/set/autoRefresh true 76 /vis/verbose warnings 77 # 78 # For file-based drivers, use this to create an empty detector view: 79 #/vis/viewer/flush
Note: See TracChangeset
for help on using the changeset viewer.