+A  Click here to enlarge/reduce to/back from full screen 
Paul Scherrer Institut PSI LMU : Laboratory for Muon Spin Spectroscopy

PSI Home NUM Home LMU Home Contact Feedback



Updated:
12.01.2010
E-Mail: thomas.prokscha@psi.ch


Printer Friendly Printout without Logo or Navigation Elements available here... just click and print

Simulations in Geant4 (musrSim)

GEANT4 as a simulation framework in µSR - musrSim


Numerical simulations in µSR

The increased use of numerical simulations in recent years has demonstrated the potential of Monte Carlo methods also for the µSR technique, where it has opened up new application possibilities by offering:
  • New insights on the basics of µSR.
  • A better understanding of the existing spectrometers, allowing a targeted performance improvement.
  • The design and optimization of new instruments.
sw_arch.png
Software architecture of the GEANT4-based musrSim and musrSimAna simulation tool.


GEANT4 as a simulation framework

GEANT4 is a Monte Carlo radiation transport toolkit which includes a complete range of functionalities required to build flexible simulation frameworks. The main features of the musrSim program:
  • Great user friendliness:
    • Only one program (executable) for different muSR instruments (different studies).
    • The geometry and other parameters defined in the macro file.
    • Ideally, no C++ knowledge required.
  • Any geometry and material:
    • Simple geometry defined in the macro text file.
    • Complex shapes (e.g. volumes with holes) can be predefined in the source, and finaly defined in the macro text file.
  • Any electromagnetic field:
    • From field map or constant field.
    • 2D, 3D and axi-symmetric.
    • Whatever field superposition.
  • Root format outputs

a. New insights on the basics of µSR

The versatility of the new tool has permitted the modelling of the incoming muon beam, the study of the outgoing positrons behaviour, the investigation of the geometrical effects, etc.
beam_spot_meas_sim (9 kB)
Measured (left) and simulated (right) beam spot of a low-energy muon beam at z = 0.

b. Simulation-based instrument design

The new, high magnetic field instruments being built at RAL and PSI, whose design is entirely based on realistic Monte Carlo simulations, represent another example of the usefulness of numerical methods for the µSR technique.
Hi_field_PSI (13 kB) LEM2_view (34 kB)
Left: High-field spectrometer for TF-µSR measurements being built at PSI. Right: Simulation of the low-energy µSR spectrometer.

Future uses of the developed platform could include also:
  • Feasibility analysis
  • Experiment planning and optimization
  • Use as an interactive didactic tool

For more information on G4 Monte Carlo simulations in µSR one can refer to the following poster and the related paper.

This project was supported by the EC under the 6th Framework Programme through the Key Action: Strengthening the European Research Area, Research Infrastructures. Contract no.: RII3-CT-2003-505925.

Source code of musrSim

The program musrSim is intended to be a general tool for the whole µSR community. One of the main improvements to the older versions of musrSim should be a general analysis program musrSimAna (a subdirectory in the musrSim directory), which could help the users to easily analyse the output of the simulation in the case of a muSR analysis (i.e. to create muSR spectra and other histograms). To compile and run musrSim, the user has to have Geant4 and Root installed. At the moment, musrSim was tested with Geant version 4.9.5 and Root version 5.32. Root version should not be critical, but Geant should be 4.9.5 or 4.9.4 (4.9.3 is probably fine too).


References

  1. GEANT4 home page
    S. Agostinelli, et al., (GEANT4 Collaboration), Nucl. Instr. and Meth. A 506 (2003) 250.
    J. Allison, et al., IEEE Trans. Nucl. Sci. 53 (2006) 270.
  2. T. Lancaster, et al., Nucl. Instr. and Meth. A 580 (2007) 1578.
  3. T. Shiroka, et al., Nucl. Instr. and Meth. A 591 (2008) 306.
  4. K. Sedlak, et al., Physica B 404 (2009) 970.
  5. K. Sedlak, et al., Physica B 404 (2009) 974.
  6. T. Shiroka, et al., Physica B 404 (2009) 966.
  7. Z. Salman, et al., Physica B 404 (2009) 978
  8. T. Paraïso, et al., Physica B 374 (2006) 498.