HomeWork Simulations of Particle Interactions with Matter - Revision history

Foretony: /* Homework 6 */

2025-03-19T16:47:56Z

Homework 6

Foretony: /* Homework 4 */

2025-02-24T18:43:03Z

Homework 4

Foretony: /* Homework 5 */

2025-02-24T18:41:33Z

Homework 5

Foretony: /* Homework 4 */

2025-02-24T17:47:56Z

Homework 4

Foretony: /* 3.) Histograms using ROOT */

2025-01-22T19:39:15Z

3.) Histograms using ROOT

Foretony: /* 3.) Histograms using ROOT */

2025-01-18T19:06:13Z

3.) Histograms using ROOT

Foretony: /* 3.) Histograms using ROOT */

2025-01-18T19:02:28Z

3.) Histograms using ROOT

Foretony: /* 3.) Histograms using ROOT */

2025-01-18T19:01:46Z

3.) Histograms using ROOT

Foretony: /* 2.) MC calculation of Pi */

2025-01-18T18:58:24Z

2.) MC calculation of Pi

Foretony: /* Homework 7 */

2025-01-15T19:15:11Z

Homework 7

← Older revision		Revision as of 16:47, 19 March 2025
Line 166:		Line 166:

	[[Media:SPIM_BremE-Spectrum-Tantalum.pdf]]		[[Media:SPIM_BremE-Spectrum-Tantalum.pdf]]
		+
		+	sudo apt-get install texlive-latex-base

	[[File:SPIM_LaTex_TemplateFile.txt]]		[[File:SPIM_LaTex_TemplateFile.txt]]

@@ Line 138: / Line 138: @@
 :<math>m_e</math> = mass of target electron initially at rest
-.) Use GEANT4 to determine the Range of the particle chosen in Homework 3 through liquid hydrogen as a function of at least three of the Energies used in Homework 3.  Use the equation from question 1 above to calculate the max energy loss for this relativistic particle and compare it to the values in the PDF file for range below.
+Use the above equation to calculate the max energy loss (transfer) for a 10 GeV proton.
 [[Image:RangeInLH2.pdf]]

← Older revision		Revision as of 18:41, 24 February 2025
Line 155:		Line 155:


−	2.) Alter GEANT4 example ~~N02~~ to check your answer for problem 1 above. I expect you to hand in a screen shot showing GEANT4 tracking the proton from 600 MeV to 400 MeV.	+	2.) Alter GEANT4 example B2a to check your answer for problem 1 above. I expect you to hand in a screen shot showing GEANT4 tracking the proton from 600 MeV to 400 MeV.

	3.) Find <math>\frac{\sigma_R}{R}</math> using GEANT4 for a 600 MeV proton traveling through a slab of copper. You will need to make the copper thick enough to stop the proton. Then output the stopping distance to a file which you can read into ROOT using some of the software we used for Homework 1's RNG problem.		3.) Find <math>\frac{\sigma_R}{R}</math> using GEANT4 for a 600 MeV proton traveling through a slab of copper. You will need to make the copper thick enough to stop the proton. Then output the stopping distance to a file which you can read into ROOT using some of the software we used for Homework 1's RNG problem.

@@ Line 138: / Line 138: @@
 :<math>m_e</math> = mass of target electron initially at rest
-.) Use GEANT4 to determine the Range of the particle chosen in Homework 3 through liquid hydrogen as a function of at least three of the Energies used in Homework 3.
+.) Use GEANT4 to determine the Range of the particle chosen in Homework 3 through liquid hydrogen as a function of at least three of the Energies used in Homework 3.  Use the equation from question 1 above to calculate the max energy loss for this relativistic particle and compare it to the values in the PDF file for range below.
 [[Image:RangeInLH2.pdf]]

@@ Line 88: / Line 88: @@
 <pre>
-rns->Draw("x*4");
+Sim->Draw("evt.x*4");
-rns->Draw("x:y","x>0.5");
+Sim->Draw("evt.x:evt.y","evt.x>0.5");
-rns->Draw("asin(x)");
+Sim->Draw("asin(evt.x)");
-rns->Draw("asin(x-y)");   Does this look like a Normal/Gaussian Distributi
+Sim->Draw("asin(evt.x-evt.y)");   Does this look like a Normal/Gaussian Distributi
 on?
 </pre>

@@ Line 39: / Line 39: @@
 Create histograms of the random numbers stored a the ROOT ntuple .
-First step is to create a file with the 2 columns of random numbers in it.
+First step is to create a file with the 2 columns of random numbers in it using the PI simulation program
   ./PI

@@ Line 43: / Line 43: @@
   ./PI
-now load the program "ascii2root.C" into ROOT.
+now load the program "asci2root.C" into ROOT.
 If root is not available you can try to add the program ROOT to your UNIX path with the command below
@@ Line 56: / Line 56: @@
 you can run the asci2root program in ROOT with the command
-  root [0] .x ascii2root.C
+  root [0] .x asci2root.C
 You now have a ROOT file called "sim.root"

@@ Line 41: / Line 41: @@
 First step is to create a file with the 2 columns of random numbers in it.
-  ./PI_MC 1000 > temp.dat
+  ./PI
 now load the program "ascii2root.C" into ROOT.
-Add the program ROOT to your UNIX path with the command below
+If root is not available you can try to add the program ROOT to your UNIX path with the command below
 source ~foretony/src/ROOT/root-6.06.02/bin/thisroot.sh
@@ Line 54: / Line 54: @@
   root -l
-now load the asci2root program into ROOT
+you can run the asci2root program in ROOT with the command
-  root [0] .L ascii2root.C
+  root [0] .x ascii2root.C
-now run the program so it reads the data in from the file "temp.dat" created above
+You now have a ROOT file called "sim.root"
 for simplicity lets quite root and restart it
@@ Line 71: / Line 67: @@
 now restart it and tell it to load the root file "rns.root" .
-  root -l rns.root
+  root -l sim.root
 You don't have to do it this way. You could use the browser to load it "new TBrowser()".
@@ Line 79: / Line 75: @@
 Look under the "ROOT Files" menu and you will see and icon labeled "rns.root".
-Double click on the icon labeled "rns.root"
+Double click on the icon labeled "sim.root"
-You now see a folder icon named "rns"
+You now see a folder icon named "Sim"
 Double click again.
-You now see a leaf icon with the names "rnd1" and "rnd2".  These should contain the numbers from the file temp.dat.
+You now see a leaf icon with the names "x" and "y".  These should contain the numbers from the file sim.dat.
 double click on one of the leaf icons and a histogram will appear.
@@ Line 92: / Line 88: @@
 <pre>
-rns->Draw("rnd1*4");
+rns->Draw("x*4");
-rns->Draw("rnd1:rnd1","rnd1>0.5");
+rns->Draw("x:y","x>0.5");
-rns->Draw("asin(rnd1)");
+rns->Draw("asin(x)");
-rns->Draw("asin(rnd1-rnd2)");   Does this look like a Normal/Gaussian Distributi
+rns->Draw("asin(x-y)");   Does this look like a Normal/Gaussian Distributi
 on?
 </pre>
 =Homework 2=

← Older revision		Revision as of 18:58, 18 January 2025
Line 35:		Line 35:

	You may use the program I gave you in the lab for this but you must change the program to indicate that you were able to understand its structure and you were able to recompile it.		You may use the program I gave you in the lab for this but you must change the program to indicate that you were able to understand its structure and you were able to recompile it.
−
−
−	~~Xtra help: http://physics.isu.edu/~tforest/Classes/NucSim/Day1/PI_MC/~~

	==3.) Histograms using ROOT==		==3.) Histograms using ROOT==

← Older revision		Revision as of 19:15, 15 January 2025
Line 187:		Line 187:

	https://indico.cern.ch/event/629841/contributions/2712690/attachments/1518832/2371867/PE_models_G4GV.pdf		https://indico.cern.ch/event/629841/contributions/2712690/attachments/1518832/2371867/PE_models_G4GV.pdf
		+
		+	https://geant4-userdoc.web.cern.ch/UsersGuides/PhysicsReferenceManual/html/index.html

	= Homework 8 =		= Homework 8 =