Difference between revisions of "Simulations and Backgrounds"

From Ciswikidb
Jump to navigation Jump to search
Line 57: Line 57:
 
## $ '''git config --global user.name "Your Name"'''
 
## $ '''git config --global user.name "Your Name"'''
 
## $ '''git config --global user.email your_git_email@example.com'''
 
## $ '''git config --global user.email your_git_email@example.com'''
−
# Done once at the beginning to add the master to your disk
+
# Done once at the beginning to add the master to your disk:
−
## $ '''cd detectors/injector/''
+
## $ '''cd detectors/injector/'''
 
## $ '''git remote add upstream https://github.com/gemc/source.git'''
 
## $ '''git remote add upstream https://github.com/gemc/source.git'''
 
## $ '''git fetch upstream'''
 
## $ '''git fetch upstream'''
Line 65: Line 65:
  
 
* How to work with '''git''':
 
* How to work with '''git''':
−
# to update your fork at github.com:
+
# To update your fork at github.com:
 
## $ '''git add ''newfile'' ''' (if any)
 
## $ '''git add ''newfile'' ''' (if any)
 
## $ '''git commit -a -m"comment"'''
 
## $ '''git commit -a -m"comment"'''

Revision as of 11:45, 8 December 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

G4 CrossSection.gif



gemc - Now (Fall 2015)

  • We will use GEant4 Monte Carlo GEMC: [2]


  • To try:
  1. $ ssh ifarm
  2. $ mkdir /group/bubble/<user>
  3. $ cd /group/bubble/<user>
  4. $ source /site/12gev_phys/production.csh 1.3 (for more details, see [3])
  5. $ gemc


  • Simple Geometry and Material Tutorial: [4]


  • Drawings used to generate gemc Bubble Chamber Model:
  1. Beamline Drawings:
    1. Copper Radiator: media:JL0015733_CU_RADIATOR_WINDOW_ASSY.pdf
    2. Copper Collimator: media:JL0015777_CU_COLLIMATOR.pdf
    3. Photon Dump: media:JL0009386_ALUMINUM_PHOTON_BEAM_DUMP.pdf
  2. Bubble Chamber Drawings:
    1. High Pressure Vessel: media:NPD-20-001-R1.pdf
    2. Pressure Vessel Blank Flange: media:NPD-20-002.pdf
    3. Gamma Beam Window: media:NPD-20-011.pdf
    4. Blank Flange MOD - Gamma Port: media:NPD-20-014.pdf
    5. Glass Cell - Metal Adapter: media:NPD-21-106.pdf
    6. Bubble Chamber Dimensions: media:Bubble_Chamber_Reference_Dimensions.pdf


  • How to start with git:
  1. Done once at the beginning on web:
    1. Create an account (user-name and password) at https://github.com/
    2. Fork gemc/detectors
  2. Done once at the beginning in your work directory:
    1. $ ssh ifarm
    2. $ git clone https://github.com/your_git_user_name/detectors.git
    3. $ git config --global user.name "Your Name"
    4. $ git config --global user.email your_git_email@example.com
  3. Done once at the beginning to add the master to your disk:
    1. $ cd detectors/injector/
    2. $ git remote add upstream https://github.com/gemc/source.git
    3. $ git fetch upstream
    4. $ git checkout master


  • How to work with git:
  1. To update your fork at github.com:
    1. $ git add newfile (if any)
    2. $ git commit -a -m"comment"
    3. $ git push
  2. To sync your files at disk with the master:
    1. $ git merge upstream/master (update disk copy from master: https://github.com/gemc/detectors.git)


  • To ask the project manager to incorporate your changes into the “master":
  1. Create pull request: done with website, navigate to your fork - create pull request by pushing green button


  • How to run gemc:
  1. $ cd detectors/injector
  2. $ source /site/12gev_phys/production.csh 1.3
  3. $ ./geometry.pl config.dat
  4. $ gemc bubble.gcard


  • To check overlaps:
  1. $ gemc bubble.gcard -USE_GUI=0 -CHECK_OVERLAPS=10


  • To check overlaps:
  1. $ gemc bubble.gcard -USE_GUI=0 -CHECK_OVERLAPS=10


  • To draw axes: /vis/scene/add/axes 0 0 0 0.5 m (Draws axes at (x0, y0, z0) of given length and colour – x, y and z will be red, green and blue)


  • To print physics summary: gemc bubble.gcard -PHYS_VERBOSITY=10


  • To run gemc in batch mode:
  1. $ gemc bubble.gcard -USE_GUI=0 -N=1000000 -PRINT_EVENT=100000


  • To convert to root:
  1. $ $BANKS/bin/evio2root -INPUTF=out.ev


  • The current geometry for the Bubble Chamber is in gemc and is shown below:


Bubble gemc.png





First Attempt - Spring 2014



Second Attempt - Summer 2014

  • Use github.com to manage the simulation software
  • The current geometry for the Bubble Chamber is in gdml and is shown below:

BCjul7.jpeg


BC1jul7.jpeg


Using ifarm

  • Login remotely to the JLab gateway: $ ssh login.jlab.org
  • Connect to either ifarm or Jlabl4:
  1. $ ssh ifarm
  2. $ ssh jlabl4
  • Load Geant4 and ROOT Environments by running the following commands: (needed every login) (see [5] for more details)
  1. $ 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 Simulation

  1. $ mkdir /group/bubble/<user>
  2. $ cd /group/bubble/<user>
  3. $ git clone https://github.com/JeffersonLab/BubbleEx.git
  4. $ cd BubbleEx
  5. $ git pull (to update files from github)
  6. $ mkdir build
  7. $ cd build
  8. $ use cmake (only if using JLab nodes)
  9. $ cmake ..
  10. $ make
  11. $ cd ..


Running Simulation

  • In the install directory
  1. $ 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:
  1. $ /vis/open OGL
  2. $ /run/beamOn 100
  • Initial run of the simulation is shown:

BubbleEx Output.jpg



To-do List

  1. Geometry is in gdml format, need to get correct Radiator, Collimator, Beamline, PhotonDump geometry.
  2. Add sweep magnet and electrondump geometry in gdml format.
  3. Complete: Include the Bubble Chamber geometry using the Bubble Chamber STEP file (change .txt to .stp): media:BubbleChamber.txt.
    1. Geometry has been created and is on github, BubbleMother.gdml and all daughter files
  4. The beamline STEP file will be ready very soon.
  5. GEANT4 uses models that calculate wrong photo–nuclear cross sections (see above). GEANT4 does not allow for user’s cross sections. What to do?
    1. Use GEANT4 to produce the photon spectra impinging on the superheated liquid.
    2. Fold the above photon spectra with our cross sections in stand-alone codes.
    3. Could we not create a correct physicslist.cc in geant4 instead?



Relevant Theoretical Cross Sections

  • Photo-Nuclear Cross Sections:
CrossSections.gif
  1. 16O(γ,α)12C media:cs16Oga.txt (Eg (MeV), Cross Section (b))
  2. 17O(γ,α)13C media:cs17Oga.txt
  3. 18O(γ,α)14C media:cs18Oga.txt
  4. 14N(γ,p)13C media:cs14Ngp.txt


  • Neutron Elastic Cross Sections:
  1. 16O(n,n) media:cs16Onn.txt (En (MeV), Cross Section (b))
  2. 14N(n,n) media:cs14Nnn.txt
Cs16Onn.gif
Cs14Nnn.gif


  • Fluorine Cross Section:
  1. 19F(γ,α)15N media:cs19Fga.txt (Eg (MeV), Cross Section (b))
Cs19Fga.gif
Retrieved from "https://wiki.jlab.org/ciswiki/index.php?title=Simulations_and_Backgrounds&oldid=7632"