Difference between revisions of "Simulations and Backgrounds"
Jump to navigation
Jump to search
Line 9: | Line 9: | ||
= '''GEANT4 Simulation''' = | = '''GEANT4 Simulation''' = | ||
− | * GEANT4 Photon-Inelastic Model: CHIPS Gamma Nuclear has this cross section for our reaction: | + | * GEANT4 Photon-Inelastic Model: CHIPS Gamma Nuclear has this (wrong) cross section for our reaction: |
[[file:G4_CrossSection.gif||500px|]] | [[file:G4_CrossSection.gif||500px|]] |
Revision as of 11:37, 4 November 2015
GEANT4 vs. Experimental Data
- Benchmarking of monte carlo simulation of bremsstrahlung from thick targets at radiotherapy energies
B. A. Faddegon, M. Asai, J. Perl, C. Ross, J. Sempau, J. Tinslay, and F. Salvat. Med. Phys. 35 (2008) 4308 [1] media:MedPhys.35.4308.pdf
GEANT4 Simulation
- GEANT4 Photon-Inelastic Model: CHIPS Gamma Nuclear has this (wrong) cross section for our reaction:
First Attempt
- Quick start for the simulation: media:G4_BubbleEx_README.txt. There is also a README file that comes with the example.
- Bubble Chamber GEANT4 Example (change .txt to .tar): media:BubbleEx96.txt
Second Attempt
- Using github.com to manage the simulation software
- The current geometry for the Bubble Chamber is in gdml and is shown below:
Using the ifarm or jlabl4 nodes
- Login remotely to the JLab gateway: $ ssh login.jlab.org
- Connect to either ifarm or Jlabl4:
- $ ssh ifarm
- $ ssh jlabl4
- Load Geant4 and ROOT Environments by running the following commands: (needed every login) (see [2] for more details)
- $ source /site/12gev_phys/production.csh 1.1
- Should see a short list of programs such as ROOT and Geant4 listed:
> Common Environment Version: <1.1> (Mon, 20 Oct 2014) > Running as suleiman on ifarm1401 > OS Release: Linux_CentOS6.5-x86_64-gcc4.4.7 > JLAB_ROOT set to: /site/12gev_phys > JLAB_SOFTWARE set to: /site/12gev_phys/1.1/Linux_CentOS6.5-x86_64-gcc4.4.7 > CLHEP version: 2.1.3.1 > XERCES version: 3.1.1 > QTDIR version: system > Geant4 version: 4.9.6.p03 > ROOT version: 5.34.13 > GEMC version: 2.1 > JANA version: 0.7.2 > Build version: 1.1 > EVIO version: 4.0 > Banks version: 1.0
Compiling the simulation
- Simulation is available at https://github.com/JeffersonLab/BubbleEx
- Create directory for simulation (on JLab nodes can use /group/bubble/ directory):
- $ mkdir /group/bubble/<user>
- $ cd /group/bubble/<user>
- $ git clone https://github.com/JeffersonLab/BubbleEx.git
- $ cd BubbleEx
- $ git pull (to update files from github)
- $ mkdir build
- $ cd build
- $ use cmake (only if using JLab nodes)
- $ cmake ..
- $ make
- $ cd ..
Running the simulation
- In the install directory
- $ build/BubbleEx (remember to run the source command before running)
- Terminal will list what has been loaded and then display output command box
- In session run:
- $ /vis/open OGL
- $ /run/beamOn 100
- Initial run of the simulation is shown:
To-do List
- Geometry is in gdml format, need to get correct Radiator, Collimator, Beamline, PhotonDump geometry.
- Add sweep magnet and electrondump geometry in gdml format.
- Complete: Include the Bubble Chamber geometry using the Bubble Chamber STEP file (change .txt to .stp): media:BubbleChamber.txt.
- Geometry has been created and is on github, BubbleMother.gdml and all daughter files
- The beamline STEP file will be ready very soon.
- GEANT4 uses models that calculate wrong photo–nuclear cross sections (see above). GEANT4 does not allow for user’s cross sections. What to do?
- Use GEANT4 to produce the photon spectra impinging on the superheated liquid.
- Fold the above photon spectra with our cross sections in stand-alone codes.
- Could we not create a correct physicslist.cc in geant4 instead?
Relevant Cross Sections
- Photo-Nuclear Cross Sections:
- 16O(γ,α)12C media:cs16Oga.txt (Eg (MeV), Cross Section (b))
- 17O(γ,α)13C media:cs17Oga.txt
- 18O(γ,α)14C media:cs18Oga.txt
- 14N(γ,p)13C media:cs14Ngp.txt
- Neutron Elastic Cross Sections:
- 16O(n,n) media:cs16Onn.txt (En (MeV), Cross Section (b))
- 14N(n,n) media:cs14Nnn.txt
- Fluorine Cross Section:
- 19F(γ,α)15N media:cs19Fga.txt (Eg (MeV), Cross Section (b))