Changeset 1337 for trunk/examples/extended/optical/LXe
- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- Location:
- trunk/examples/extended/optical/LXe
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/optical/LXe/History
r807 r1337 1 $Id: History,v 1. 9 2007/09/30 22:33:42gum Exp $1 $Id: History,v 1.14 2010/06/03 21:23:13 gum Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 03 June 2010 Peter Gumplinger (LXe-V09-03-00) 19 - modify LXeTrajectory class to compile without warnings against 20 tracking-V09-03-03 by implementing a DrawTrajectory() method. 21 22 09 November 2009 Peter Gumplinger (LXe-V09-02-01) 23 - exhibit the usage of a complex index of refraction for the photo cathode 24 show usage of G4OpticalPhysics builder/messenger (instead of LXeOpticalPhysics) 25 26 30 October 2009 John Allison (LXe-V09-02-00) 27 - Introduced G4UIExecutive. 28 29 20 November 2008 Peter Gumplinger (LXe-V09-01-01) 30 - update code in LXeOpticalPhysics.cc 31 32 16 July 2008 Peter Gumplinger (LXe-V09-01-00) 33 - add Birks Constant for the LXe scintillator, Polystyrene scintillator 34 modified: LXeDetectorConstruction.cc, LXeOpticalPhysics.cc 17 35 18 36 30 Sept 2007 Peter Gumplinger (LXe-V09-00-00) -
trunk/examples/extended/optical/LXe/LXe.in
r807 r1337 8 8 /LXe/eventVerbose 1 9 9 #/tracking/verbose 1 10 /run/initialize 11 #/optics_engine/selectOpProcess OpBoundary 12 #/optics_engine/setOpProcessVerbose 1 10 13 /run/beamOn 3 -
trunk/examples/extended/optical/LXe/LXe.out
r807 r1337 1 ********************************************** 2 Geant4 version $Name: $ 3 (12-December-2003) 4 Copyright : Geant4 Collaboration 5 ********************************************** 6 Visualization Manager instantiating... 1 2 ############################################ 3 !!! WARNING - FPE detection is activated !!! 4 ############################################ 5 6 ************************************************************* 7 Geant4 version Name: geant4-09-03-ref-06 (25-June-2010) 8 Copyright : Geant4 Collaboration 9 Reference : NIM A 506 (2003), 250-303 10 WWW : http://cern.ch/geant4 11 ************************************************************* 12 13 Visualization Manager instantiating with verbosity "warnings (3)"... 7 14 Visualization Manager initialising... 8 15 Registering graphics systems... 9 16 10 You have successfully chosen to usethe following graphics systems.17 You have successfully registered the following graphics systems. 11 18 Current available graphics systems are: 12 19 ASCIITree (ATree) 13 20 DAWNFILE (DAWNFILE) 14 G AGTree (GAGTree)21 G4HepRep (HepRepXML) 15 22 G4HepRepFile (HepRepFile) 16 G4HepRep (HepRepXML)17 23 RayTracer (RayTracer) 18 24 VRML1FILE (VRML1FILE) 19 25 VRML2FILE (VRML2FILE) 26 gMocrenFile (gMocrenFile) 27 OpenGLStoredX (OGL) 28 OpenGLImmediateX (OGLI) 29 OpenGLStoredX (OGLS) 20 30 OpenGLImmediateX (OGLIX) 21 31 OpenGLStoredX (OGLSX) 22 OpenGLImmediateXm (OGLIXm) 23 OpenGLStoredXm (OGLSXm) 24 25 phot: Total cross sections from Sandia parametrisation. 26 MuonMinusCaptureAtRest is created 27 28 compt: Total cross sections from a parametrisation. Good description from 10 KeV to (100/Z) GeV. 29 Scattered gamma energy according Klein-Nishina. 30 PhysicsTables from 1 keV to 100 GeV in 80 bins. 31 32 conv: Total cross sections from a parametrisation. Good description from 1.5 MeV to 100 GeV for all Z. 33 e+e- energies according Bethe-Heitler 34 PhysicsTables from 1.022 MeV to 100 GeV in 100 bins. 35 36 msc: Model variant of multiple scattering for e- 37 Lambda tables from 100 eV to 100 TeV in 120 bins. 38 Boundary algorithm is active with facrange= 0.199 39 40 eIoni: tables are built for e- 41 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 42 Lambda tables from threshold to 100 TeV in 120 bins. 43 Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 1 44 Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. 45 46 eBrem: tables are built for e- 47 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 48 Lambda tables from threshold to 100 TeV in 120 bins. 49 Total cross sections from a parametrisation based on the EEDL data library. 50 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. 51 52 eIoni: tables are built for e+ 53 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 54 Lambda tables from threshold to 100 TeV in 120 bins. 55 Step function: finalRange(mm)= 1, dRoverRange= 1, integral: 1 56 Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV. 57 58 eBrem: tables are built for e+ 59 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 60 Lambda tables from threshold to 100 TeV in 120 bins. 61 Total cross sections from a parametrisation based on the EEDL data library. 62 Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. 63 64 annihil: Total cross section from Heilter formula(annihilation into 2 photons). 65 gamma energies sampled according Heitler 66 PhysicsTables from 10 keV to 10 TeV in 100 bins. 67 68 msc: Model variant of multiple scattering for mu+ 69 Lambda tables from 100 eV to 100 TeV in 120 bins. 70 Boundary algorithm is active with facrange= 0.199 71 72 MuIoni: tables are built for mu+ 73 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 74 Lambda tables from threshold to 100 TeV in 120 bins. 75 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 0 76 Bether-Bloch model for E > 0.2 MeV parametrisation of Bragg peak below. 77 78 MuBrems: tables are built for mu+ 79 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 80 Lambda tables from threshold to 100 TeV in 120 bins. 81 Parametrised model 82 83 MuPairProd: tables are built for mu+ 84 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 85 Lambda tables from threshold to 100 TeV in 120 bins. 86 Parametrised model 87 ### All dEdx and Range tables are built ##### 88 89 MuIoni: tables are built for mu- 90 dE/dx and range tables from 100 eV to 100 TeV in 120 bins. 91 Lambda tables from threshold to 100 TeV in 120 bins. 92 Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 0 93 Bether-Bloch model for E > 0.2 MeV parametrisation of Bragg peak below. 94 Energy weighted position of hits in LXe : (5.18102,30.2388,-98.8831) 32 33 Registering model factories... 34 35 You have successfully registered the following model factories. 36 Registered model factories: 37 generic 38 drawByCharge 39 drawByParticleID 40 drawByOriginVolume 41 drawByAttribute 42 43 Registered filter factories: 44 chargeFilter 45 particleFilter 46 originVolumeFilter 47 attributeFilter 48 49 G4Cerenkov::G4Cerenkov constructor 50 NOTE: this is now a G4VProcess! 51 Required change in UserPhysicsList: 52 change: pmanager->AddContinuousProcess(theCerenkovProcess); 53 to: pmanager->AddProcess(theCerenkovProcess); 54 pmanager->SetProcessOrdering(theCerenkovProcess,idxPostStep); 55 56 phot: for gamma SubType= 12 57 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 58 PhotoElectric : Emin= 0 eV Emax= 10 TeV 59 60 compt: for gamma SubType= 13 61 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 62 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 63 Klein-Nishina : Emin= 0 eV Emax= 10 TeV 64 65 conv: for gamma SubType= 14 66 Lambda tables from 1.022 MeV to 10 TeV in 77 bins, spline: 1 67 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 68 Bethe-Heitler : Emin= 0 eV Emax= 10 TeV 69 70 msc: for e- SubType= 10 71 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 72 RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1 73 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 74 UrbanMsc93 : Emin= 0 eV Emax= 10 TeV 75 76 eIoni: for e- SubType= 2 77 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 78 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 79 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 80 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 81 MollerBhabha : Emin= 0 eV Emax= 10 TeV 82 83 eBrem: for e- SubType= 3 84 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 85 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 86 LPM flag: 1 for E > 1 GeV 87 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 88 eBrem : Emin= 0 eV Emax= 1 GeV 89 eBremRel : Emin= 1 GeV Emax= 10 TeV 90 91 eIoni: for e+ SubType= 2 92 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 93 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 94 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 95 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 96 MollerBhabha : Emin= 0 eV Emax= 10 TeV 97 98 eBrem: for e+ SubType= 3 99 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 100 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 101 LPM flag: 1 for E > 1 GeV 102 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 103 eBrem : Emin= 0 eV Emax= 1 GeV 104 eBremRel : Emin= 1 GeV Emax= 10 TeV 105 106 annihil: for e+ SubType= 5 107 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 108 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 109 eplus2gg : Emin= 0 eV Emax= 10 TeV 110 111 muMsc: for mu+ SubType= 10 112 Lambda tables from 100 eV to 10 TeV in 77 bins, spline: 1 113 RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 0 114 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 115 UrbanMsc90 : Emin= 0 eV Emax= 10 TeV 116 117 muIoni: for mu+ SubType= 2 118 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 119 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 120 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 121 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 122 Bragg : Emin= 0 eV Emax= 200 keV 123 BetheBloch : Emin= 200 keV Emax= 1 GeV 124 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 125 126 muBrems: for mu+ SubType= 3 127 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 128 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 129 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 130 MuBrem : Emin= 0 eV Emax= 10 TeV 131 132 muPairProd: for mu+ SubType= 4 133 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 134 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 135 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 136 muPairProd : Emin= 0 eV Emax= 10 TeV 137 138 muIoni: for mu- SubType= 2 139 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 140 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 141 finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 142 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 143 ICRU73QO : Emin= 0 eV Emax= 200 keV 144 BetheBloch : Emin= 200 keV Emax= 1 GeV 145 MuBetheBloch : Emin= 1 GeV Emax= 10 TeV 146 147 muBrems: for mu- SubType= 3 148 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 149 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 150 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 151 MuBrem : Emin= 0 eV Emax= 10 TeV 152 153 muPairProd: for mu- SubType= 4 154 dE/dx and range tables from 100 eV to 10 TeV in 77 bins 155 Lambda tables from threshold to 10 TeV in 77 bins, spline: 1 156 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 157 muPairProd : Emin= 0 eV Emax= 10 TeV 158 ============================================================================================ 159 HADRONIC PROCESSES SUMMARY (verbose level 1) 160 161 Hadronic Processes for <mu-> 162 muMinusCaptureAtRest 163 ============================================================================================ 164 ### G4EmSaturation::FindBirksCoefficient Birks coefficient for LXe 0.126 mm/MeV 165 Energy weighted position of hits in LXe : (0.106357,-0.0981109,-103.428) 95 166 Total energy deposition in scintillator : 511 (keV) 96 Reconstructed position of hits in LXe : (3.87256,8.56373,-27.017) 97 Number of photons that hit PMTs in this event : 2160 167 Reconstructed position of hits in LXe : (4.05096,0.730982,-14.686) 168 WARNING: G4VisManager::IsValidView(): Attempt to draw when no graphics system 169 has been instantiated. Use "/vis/open" or "/vis/sceneHandler/create". 170 Alternatively, to avoid this message, suppress instantiation of vis 171 manager (G4VisExecutive), possibly by setting G4VIS_NONE, and ensure 172 drawing code is executed only if G4VVisManager::GetConcreteInstance() 173 is non-zero. 174 Number of photons that hit PMTs in this event : 1398 98 175 Number of PMTs above threshold(1) : 32 99 Number of photons produced by scintillation in this event : 6040100 Number of photons produced by cerenkov in this event : 1101 Number of photons absorbed (OpAbsorption) in this event : 3881176 Number of photons produced by scintillation in this event : 5625 177 Number of photons produced by cerenkov in this event : 2 178 Number of photons absorbed (OpAbsorption) in this event : 4229 102 179 Number of photons absorbed at boundaries (OpBoundary) in this event : 0 103 180 Unacounted for photons in this event : 0 104 Energy weighted position of hits in LXe : ( -28.7618,-1.33268,-9.54847)181 Energy weighted position of hits in LXe : (0.148663,0.0658336,-97.3519) 105 182 Total energy deposition in scintillator : 511 (keV) 106 Reconstructed position of hits in LXe : ( -5.96696,-2.2751,-2.77019)107 Number of photons that hit PMTs in this event : 2176183 Reconstructed position of hits in LXe : (1.15832,-1.1588,-14.6319) 184 Number of photons that hit PMTs in this event : 1349 108 185 Number of PMTs above threshold(1) : 32 109 Number of photons produced by scintillation in this event : 6167186 Number of photons produced by scintillation in this event : 5502 110 187 Number of photons produced by cerenkov in this event : 1 111 Number of photons absorbed (OpAbsorption) in this event : 3992188 Number of photons absorbed (OpAbsorption) in this event : 4154 112 189 Number of photons absorbed at boundaries (OpBoundary) in this event : 0 113 190 Unacounted for photons in this event : 0 114 Energy weighted position of hits in LXe : ( 1.6474,-12.2959,-71.5899)191 Energy weighted position of hits in LXe : (-0.148027,0.14522,-99.4255) 115 192 Total energy deposition in scintillator : 511 (keV) 116 Reconstructed position of hits in LXe : ( 1.21798,-3.59905,-26.1855)117 Number of photons that hit PMTs in this event : 2230193 Reconstructed position of hits in LXe : (-1.74264,-1.85574,-16.0139) 194 Number of photons that hit PMTs in this event : 1484 118 195 Number of PMTs above threshold(1) : 32 119 Number of photons produced by scintillation in this event : 6217196 Number of photons produced by scintillation in this event : 5629 120 197 Number of photons produced by cerenkov in this event : 1 121 Number of photons absorbed (OpAbsorption) in this event : 3988198 Number of photons absorbed (OpAbsorption) in this event : 4146 122 199 Number of photons absorbed at boundaries (OpBoundary) in this event : 0 123 200 Unacounted for photons in this event : 0 124 201 Graphics systems deleted. 125 202 Visualization Manager deleting... 203 -
trunk/examples/extended/optical/LXe/WALKTHROUGH
r807 r1337 16 16 #3) Run a basic event 17 17 18 #3a)Turn on verbose output at end of event 18 #3a) Run initialize - do this only once 19 /run/initialize 20 21 #3b)Turn on verbose output at end of event 19 22 # We'll leave it on for the rest of the events in the walkthrough too 20 23 /LXe/eventVerbose 1 21 24 22 #3 b)Run25 #3c)Run 23 26 /run/beamOn 24 27 -
trunk/examples/extended/optical/LXe/include/LXeTrajectory.hh
r807 r1337 45 45 LXeTrajectory(LXeTrajectory &); 46 46 virtual ~LXeTrajectory(); 47 47 48 virtual void DrawTrajectory() const; 48 49 virtual void DrawTrajectory(G4int i_mode=0) const; 49 50 -
trunk/examples/extended/optical/LXe/src/LXeTrajectory.cc
r807 r1337 64 64 { 65 65 66 } 67 68 //_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_ 69 void LXeTrajectory::DrawTrajectory() const 70 { 71 // Invoke the default implementation in G4VTrajectory... 72 G4VTrajectory::DrawTrajectory(); 73 // ... or override with your own code here. 66 74 } 67 75
Note: See TracChangeset
for help on using the changeset viewer.