1 / 15

Getting started with musrSim

Getting started with musrSim. November 3, 2009 Tsutomu Mibe (KEK). Preparation. unzip musrSim.zip cd musrSim make Add low energy EM data files Goto http:// geant4.web.cern.ch/geant4/support/download.shtml Click “download” su cd /opt/ heplib /Geant4/ gtar zxf G4EMLow6.2.tar.gz

kieran-travis
Download Presentation

Getting started with musrSim

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. Getting started with musrSim November 3, 2009 Tsutomu Mibe (KEK)

  2. Preparation • unzip musrSim.zip • cdmusrSim • make • Add low energy EM data files • Gotohttp://geant4.web.cern.ch/geant4/support/download.shtml • Click “download” • su • cd /opt/heplib/Geant4/ • gtarzxf G4EMLow6.2.tar.gz • exit • export G4LEDATA=/opt/heplib/Geant4/data/G4EMLOW6.2 • cd run • copy example01.mac to this directory.

  3. Example 01 • example01.mac • Target: • vacuum • Detector: • No detector • Beam: • Stopped mu+ at (z,y,z)=(0,0,0)

  4. Running example01.mac (1)

  5. Running example01.mac (2) • Close the VRML View Pro window • A mu+ generates and decays. • A new window with decay particles pops-up. Geant session ends after this window is closed.

  6. Turning off visualization • The visualization of the geometry and particles can be turned off by modifying example01.mac as shown on the right.

  7. Increasing number of events • Number of events can be changed by changing the number after /run/beamOn.

  8. Let’s run again • No window appears at this time. • A new data file is generated under ./data/ • cd data • root musr_0.root • “t1” is a name of the output root tree. • A command t1->Print() will provide list of contents

  9. positron momentum • x-component of the positoron momentum vector • t1->Draw(“posIniMomX”)

  10. positron momentum • t1->Draw(“sqrt(posIniMomX^2+posIniMomY^2+posIniMomZ^2)”)

  11. positron angle • Polar angle (angle between z-axis and positron) • t1->Draw(“posIniMomZ/sqrt(posIniMomX^2+posIniMomY^2+posIniMomZ^2)”) cosq

  12. End root session

  13. Exercises • Generate 10000 events. • Draw following distributions • positron angle vs momentum • muon decay positions (muDecayPos*) (*=X,Y,Z) • muon polarization at decay (muDecayPol*) • muon decay time (muDecayTime) • Fit the exponential function to the muon decay time distribution and obtain lifetime • Check positron angler distribution and momentum distribution is consistent with the expectation from the weak interaction theory.

  14. appendix

  15. My G4 environmental variables

More Related