Class Admin

Forest_SimPart_Syllabus

The ability to simulate nature using a computer has become a useful skill for physicists who work in disciplines ranging from basic research to video games. This class will teach you these fundamentals using a practical approach, based on a specific UNIX based programming environment, to simulate fundamental physics processes ranging from ionization and the photoelectric effect to producing anti-matter.

Homework

Homework is due at the beginning of class on the assigned day. If you have a documented excuse for your absence, then you will have 24 hours to hand in the homework after released by your doctor.

Class Policies

http://wiki.iac.isu.edu/index.php/Forest_Class_Policies

Instructional Objectives

Course Catalog Description: Simulations of Particle Interactions with Matter 3 credits. Lecture course with monte-carlo computation requirements. Topics include: Stopping power, interactions of electrons and photons with matter, hadronic interactions, and radiation detection devices.

Prequisites:Math 3360. Phys 3301 or 5561.

Course Description: A practical course applying theoretical descriptions of fundamental particle interactions such as the photoelectric effect, compton scattering and pair production to describe multiple interactions of particles with matter using the montecarlo method. A software package known as GEANT from CERN will be used that is freely available under the UNIX environment. The course assumes that the student has very limited experience with the UNIX environment and no experience with GEANT. Homework problems involve modifying and compiling example programs written in C++. A final project is required in which the student chooses a process to compare the predictions of GEANT with experimental data. A report is written in a format that would be publishable in a scientific journal. Publishing the report is not required but left as an option for the student.

Objectives and Outcomes

SPIM_ObjectiveNoutcomes

Homework Problems

HomeWork_Simulations_of_Particle_Interactions_with_Matter

Introduction

TF_SPIM_Intro

Energy Loss

TF_SPIM_StoppingPower

Ann. Phys. vol. 5, 325, (1930)

Interactions of Electrons and Photons with Matter

TF_SPIM_e-gamma

Hadronic Interactions

Proton Bremsstrahlung

= cross section for dipole radiation emitted at 90 degrees with respect to incident beam of particles scattered in a Coulomb field.

File:ProtonBrem Drell Huang PhysRev v99 n3 1955 pg686.pdf

Pluto event generator

A ROOT based Hadronic Simulation package based on Pluto

I installed Pluto V 5.14.1 on inca

I needed to set the environmental variables under tcsh

setenv ROOTSYS ~/src/ROOT/root
setenv PATH ${PATH}:${ROOTSYS}/bin
setenv LD_LIBRARY_PATH $ROOTSYS/lib

There is a subdirectory called "macros"

cd macros

Go to that subdirectory and type root, this will run the contents of the file "rootlogin.C"

cd macros

inca:~/src/Pluto/pluto_v5.14.1/macros> root
 *******************************************
 *                                         *
 *        W E L C O M E  to  R O O T       *
 *                                         *
 *   Version   5.17/03    30 August 2007   *
 *                                         *
 *  You are welcome to visit our Web site  *
 *          http://root.cern.ch            *
 *                                         *
 *******************************************
Compiled on 5 September 2007 for linux with thread support.
CINT/ROOT C/C++ Interpreter version 5.16.24, July 26, 2007
Type ? for help. Commands must be C++ statements.
Enclose multiple statements between { }.
 *********************************************************
 * The Pluto event generator                              
 * (C) HADES collaboration and all contributing AUTHORS   
 * www-hades.gsi.de/computing/pluto/html/PlutoIndex.html  
 * Version: 5.14.1
 * Compiled on 10 December 2008
 *********************************************************
Shared library Pluto.so loaded

to run a pp elastic model type

root [0] .x pp_elastic.C

a root ntuple is generate called "pp_elastic.root"

you can then analyze the data in the root file with

data->MakeClass();

the above command within root generates an analysis skeleton program.

using t.Show you can see the structure of the events within the ntuple. A few functions are also stored in the root tree which you can use. You can use the root file event to create an input file which GEANT4 can then use as its event generator. GEANT4 then reads the events in and propagates them through your geometry.