SHMS GEANT4 Simulations
A geant4 simulation toolkit for SHMS using GEANT 4.9.6 .--Buddhini 15:05, 17 December 2013 (UTC)
The source code can be found in hallcgit repository : https://hallcweb.jlab.org/git/?p=buddhini/shms_geant.git;a=summary Follow the instructions in the README file to compile the code.
- Physics List - Geant4 Physics List simulation engine: QGSP_BERT_HP 3.0
- The default cuts used in the physics list have a cut off value of 0.7*CLHEP::mm.
- Geometry - Geometry information is read in via GDML.
- The geometry uses a coordinate system centered at the bender parallel to the optics coordinate system.
- Currently we are only using this tool for radiation heat studies on the horizontal bender coils. Therefore the only geometry included in the source directory (shms_geom.gdml) is the horizontal bender magnet (coils, bore, cryostat and the york), the beamline, beam dump vacuum window, beam dump helium vessel and the dump face. All of these have rough geometric shapes based on the real dimensions of the components. But the program is setup such that it is flexible to be used with any geometric setup written in GDML.
- Primary Generator - The primary generator is defined in the SHMSPrimaryGeneratorAction class.
- The default kinematics is a 11 GeV electrons randomly distributed in a 2mm x2 mm square in the (X,Y) plane.
- Magnetic Field - Magnetic field values are read in from a table generated by TOSCA (Michael Moor's work). Currently we are only considering the field generated by the horizontal bender.
- The field can be turned On or OFF via a Boolean flag used in the SHMSDetectorConstruction class.
- The field tables are not included in the source directory due to their large size. Please contact me or Michael to get a copy of the tables.
- Outputs - ntuples and histograms of selected quantities for sensitive detectors are stored in rootfiles.
- ntuples (step by step information) - local and global position of the hits (cm), initial kinetic energy(MeV), final kinetic energy(MeV), energy deposition (MeV), particle ID and parent ID
- histograms (event by event information) energy deposited in the volume (MeV)
- text output with the total energy deposited on all the sensitive volumes during the run.