-
Notifications
You must be signed in to change notification settings - Fork 0
frankliuao/plot_gbeam_matlab
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Folders and files
| Name | Name | Last commit message | Last commit date | |
|---|---|---|---|---|
Repository files navigation
<head> <title> How to use plot_gbeam_matlab </title> </head> <body style="margin-left: 80px; margin-right: 80px"> <h2 align="center" style="color:green"> Author: Ao Liu (a.k.a. frankliuao) </h2> <h1 align="center"> How to use the script </h1> <p align="center"> First of all and obviously, you need to have a valid MATLAB on your computer. Then, download the .m file and move it to a directory you like. I recommend that you do "Set Path" in MATLAB and add the path to this directory to the MATLAB search path. Then, in your MATLAB Command Window, type in </p> <p align="center" style="font-weight:bold"> plot_gbeam_matlab </p> <p align="center"> Then, a window will pop out asking you for the beam distribution in the G4Beamline ASCII format. Because as a personal preference, I'd like to give beam files with the default G4BL ASCII format an extension of <span style="font-weight:bold">.gbeam</span>, plot_gbeam_matlab expects a "<span style="font-weight:bold">.gbeam</span>" file, but you may choose "All files" to read in any extension. </p> <p align="center"> Then you need to choose what to plot. If you want to create a 1-D histogram, choose which property you'd like to plot; if you want to create a 2-D histogram (scattered plot), choose "Number of particles". <p align="center"> <img src="./pics/pic0.png" width="300pt"> </p> </p> <p align="center"> Next is to choose the base, or the coordinate system to plot the property on. For example, if you chose "Number of particles" at the previous step, opt to choose from the following bases. <p align="center"> <img src="./pics/pic1.png" width="300pt"> </p> </p> <p align="center"> The next step is to choose the number of bins in your histogram. Notice that, the # of bins corresponds to each dimension, i.e. if you are plotting a 2-D histogram, the number of bins (grids) will be n-by-n. <p align="center"> <img src="./pics/pic2.png" width="300pt"> </p> </p> <p align="center"> The next window is to ask you for a momentum cut. Eventually I want to turn this into an arbitrary cut. But it is what it is now. <p align="center"> <img src="./pics/pic3.png" width="300pt"> </p> </p> <p align="center"> Then, choose which kind(s) of particles you need to plot, by their PID. Notice, that the particles won't be mixed together - even though you choose multiple particle IDs, they will be plotted separately. <p align="center"> <img src="./pics/pic4.png" width="300pt"> </p> </p> <p align="center"> Finally, you got what you want, e.g. <p align="center"> <img src="./pics/pic5.png" width="500pt"> </p> </p> <h1 align="center"> Stay tuned </h1> <p align="center"> https://github.com/frankliuao/plot_gbeam_matlab </p> <h1 align="center"> Known issues </h1> <p align="left"> 1. The windows can not be closed when you reach the 2nd selection window. The best way is to kill the run, or to just randomly select what you want to plot and proceed to the plotting step. </body>
About
No description, website, or topics provided.
Resources
Stars
Watchers
Forks
Releases
No releases published
Packages 0
No packages published