1 / 17

GAMOS tutorial Plug-in’s Exercises

GAMOS tutorial Plug-in’s Exercises. Pedro Arce Dubois CIEMAT http://fismed.ciemat.es/GAMOS. This set of exercises are meant for a primer user to practice with each of the most common components used in the simulation of detectors Geometry Visualisation Primary generator

cheche
Download Presentation

GAMOS tutorial Plug-in’s Exercises

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. GAMOS tutorialPlug-in’sExercises Pedro Arce Dubois CIEMAT http://fismed.ciemat.es/GAMOS Introducción a GEANT4

  2. This set of exercises are meant for a primer user to practice with each of the most common components used in the simulation of detectors Geometry Visualisation Primary generator Electromagnetic/hadronic physics Production cuts and user limits Scoring Sensitive detectors They represent examples on How to do simulation using GAMOS commands How to extend GAMOS adding new C++ code They are sequential Use the commands of the previous exercise (only change what indicated) Some examples use code that is compiled in the previous examples Introduction Introducción a GEANT4

  3. Ex. 1: Basic example Ex. 2: More difficult geometry Generator Ex. 3: Physics/cuts Ex. 4: User actions Ex. 5: Scoring Ex. 6: Sensitive detectors RT simulation Exercises Introducción a GEANT4

  4. Geometry: Tube radius 100cm length 200cm of Water inside Box 30cm X 30cm X 40cm of Aluminun Standard electromagnetic physics Primary particle: gamma 10 MeV at (0,0,0) in random direction Exercise 1: simple geometry Introducción a GEANT4

  5. Visualize geometry and tracks OGLIX VRML2FILE DAWNFILE Exercise 1b: visualisation Introducción a GEANT4

  6. Use low energy electromagnetic Use Penelope (Make statistics on processes/particles) Exercise 1c: use different e.m. physics Introducción a GEANT4

  7. “Calorimeter” 1cm X 1cm X 5cm LSO crystals, placed along X axis LSO: Lu 76.4016%, Si 6.1328%, O 17.4656% Block of 4X4 crystals Place three crystal blocks at (5,-5,0), (5,0,0), (5,5,0) cm “Tracker” Box 40cm X 20cm X 20cm of CO2 (50%) Al (50%) Place at (50,0,0) cm Exercise 2: more difficult geometry Introducción a GEANT4

  8. Position of primary gammas at (-200,0,0) mm Direction along X axis Exercise 2b: change generator position and direction Introducción a GEANT4

  9. Create position distribution: Along a semicircumpherence of radius 1 cm in YZ plane, centered in (0,0,0) Change radius at the input script Make it a plug-in and select it at the input script Exercise 2c: create C++ generator position distribution Introducción a GEANT4

  10. Production cut 0.1 mm in crystals, 0.01 mm in tracker Limit step to 1 mm in crystals Exercise 3: Cuts by region Introducción a GEANT4

  11. Primary particle: neutron 10 MeV Use QGSP_BERT_HP physics list Make it a plug-in Exercise 3b: Hadronic physics Introducción a GEANT4

  12. Kill event if primary particle does not leave all its energy in crystal If it exits the crystals (it might exit and enter back and leave some energy in the meantime…) Exercise 4a: User actions Introducción a GEANT4

  13. For each event count the energy deposited in the crystal when the process that defined the step is ionisation Make an histogram of it Exercise 4b: User actions Introducción a GEANT4

  14. Count number of secondaries with energy > 100 keV produced in the crystals per event Count energy deposited in the tracker per event Exercise 5: Scoring Introducción a GEANT4

  15. Make the crystals and the tracker sensitive detectors Plot hits variables with GmHitsHistosUA Save hits in a file Exercise 6: Sensitive detector / hits Introducción a GEANT4

  16. Use as primary particle F18 at position (50,0,0)mm with 1 miliCurie Make energy resolution of crystals 20% and tracker 2% Make measuring time 1 milisecond for both SD Make dead time 1 milisecond for both SD (look for number of hits and good hits in hits.root/hits.csv) Make dead time 10 nanoseconds for both SD Make energy resolution of crystals 0% and tracker 0% (look at hits output file) Exercise 6b: Study detector effects Introducción a GEANT4

  17. Assign Sensitive Detector to crystals Create hits with: Track ID Deposited energy Position Print number of hits in each event Print Track ID, deposited energy and position of each hit Exercise 6c: Create SD with C++ Introducción a GEANT4

More Related