Skip to content

Code to calculate the geometric efficiency of circular and gaussian sources on a circular detector surface

Notifications You must be signed in to change notification settings

MichaelHeines/Geometric_efficiency

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

34 Commits
 
 
 
 
 
 
 
 

Repository files navigation

This github directory contains code to determine the geometric efficiency of a circular/annular detector with circular uniform/gaussian sources, considering isotropic emission.

The program uses monte carlo methods to pick a starting location and initial direction, which are than used the extrapolate the trajectories. For each detector-trajectory pair, it is determined whether or not the emission ens up in the detector. The errors can be estimated straightforwardly with Poisson statistics.

After cloning the git, execute the command "chmod +x ./build.sh" once to set up permission to use the shell script build.sh in the emission folder.

To run, from main path: "./build.isotropic.exe 'source' 'detector'" Where 'source' can be 'uniform' or 'gaussian'; 'detector can be 'circular' or 'annular' The program will ask for some parameters: zmin/rd: the minimal distance for which the geometric efficiency will be calculated (in detector radius units; for annular = outer radius); zmax/rd: the maximal distance for which the geometric efficiency will be calculated (in detector radius units; for annular = outer radius); number of points: number of linspace points in which the geometric efficiency is calculated; source/rd: source spread (in detector radius units; for annular = outer radius); for circular source = source radius, for gaussian source = sigma; Power: 10^x monte carlo points used per distance; Detector outer/inner: ratio of outer radius to inner radius (only for annular detector); Filename: name of output Filename;

The ouput file is of csv type with columns showing: distance_from_source(detector radius units) point_source_approximation model_value relative_uncertainty(%)

For more information about the code: contact '[email protected]'

About

Code to calculate the geometric efficiency of circular and gaussian sources on a circular detector surface

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published