C program for neutron efficiency - Revision history

Oborn at 16:58, 13 May 2008

2008-05-13T16:58:36Z

Oborn at 19:12, 5 May 2008

2008-05-05T19:12:48Z

Oborn at 19:09, 5 May 2008

2008-05-05T19:09:35Z

Oborn at 19:05, 5 May 2008

2008-05-05T19:05:45Z

Oborn at 19:05, 5 May 2008

2008-05-05T19:05:24Z

Oborn at 19:04, 5 May 2008

2008-05-05T19:04:37Z

Oborn at 16:52, 5 May 2008

2008-05-05T16:52:03Z

Oborn at 16:23, 5 May 2008

2008-05-05T16:23:14Z

Oborn at 16:19, 5 May 2008

2008-05-05T16:19:16Z

New page

/**************************************************************************
* Author: Rebecca Bodily
* Date: 05/02/2008
* Organization: Idaho State University
*************************************************************************/

#include <iostream>
#include <cmath>
using namespace std;

double crossSectionHydrogenFunction();
double numberCarbonFunction(double thickness, const double
CUBICCENTIMETERCARBON);
double numberHydrogenFunction(double thickness, const double
CUBICCENTIMETERHYDROGEN);
double crossSectionCarbonFunction();
void efficiencyFunction(double numberHydrogen, double
numberCarbon, double
crossSectionHydrogen, double
crossSectionCarbon, double thickness);

/*************************************************************************
* This program is to find the
* efficiency of a neutron detector.
************************************************************************/
double main()
{
double neutronEnergy;
double numberHydrogen;
double numberCarbon;
double crossSectionHydrogen;
double crossSectionCarbon;
double thickness;
const double CUBICCENTIMETERCARBON = 5.23e22;
const double CUBICCENTIMETERHYDROGEN = 4.74e22;

cout << endl;
cout << "Enter thickness of detector (in centimeters): ";
cin >> thickness;

crossSectionHydrogen = crossSectionHydrogenFunction();
numberCarbon = numberCarbonFunction(thickness,
CUBICCENTIMETERCARBON);
numberHydrogen = numberHydrogenFunction(thickness,
CUBICCENTIMETERHYDROGEN);
crossSectionCarbon = crossSectionCarbonFunction();

efficiencyFunction(numberHydrogen, numberCarbon,
crossSectionHydrogen, crossSectionCarbon,
thickness);
}

/***********************************************************************
* Equation for cross section of hydrogen according to
* equation 15-9 from Radiation Detection and Measurement
* by Glenn Knoll
*********************************************************************/
double crossSectionHydrogenFunction()
{
double neutronEnergy;
double crossSectionHydrogen;

cout << endl;
cout << "Enter neutron energy (in MeV): ";
cin >> neutronEnergy;

crossSectionHydrogen = ((4.83 / sqrt(neutronEnergy))
- .578e-24);

return crossSectionHydrogen;
}

/**********************************************************************
* Equation for calculating the number of Carbons in
* the problem per cubic centimeter, dependent on what
* the user put in for the thickness of the detector.
*********************************************************************/
double numberCarbonFunction(double thickness, const
double CUBICCENTIMETERCARBON)
{
double numberCarbon;

numberCarbon = thickness * CUBICCENTIMETERCARBON;
return numberCarbon;
}

/*********************************************************************
* Equation for calculating the number of Hydrogens in
* the problem per cubic centimeter, dependent on what
* the user put in for the thickness of the detector.
********************************************************************/
double numberHydrogenFunction(double thickness, const
double CUBICCENTIMETERHYDROGEN)
{
double numberHydrogen;

numberHydrogen = thickness * CUBICCENTIMETERHYDROGEN;
return numberHydrogen;
}

/*********************************************************************
* Requesting the cross section from the Knoll book, and
* then calculating it from barns to centimeters.
********************************************************************/
double crossSectionCarbonFunction()
{
double crossSectionCarbon;

cout << endl;
cout << "Enter the cross section of Carbon by "
<< "looking in Radiation" << endl
<< "Detection and Measurement by Glenn Knoll, on"
<< " page 535, " << endl
<< "Figure 15-15b (in square centimeters): ";
cin >> crossSectionCarbon;

return crossSectionCarbon;
}

/********************************************************************
* Equation for efficiency according to equation 15-8b
* from Radiation Detection and Measurement
* by Glenn Knoll
******************************************************************/
void efficiencyFunction(double numberHydrogen, double
numberCarbon, double crossSectionHydrogen,
double crossSectionCarbon, double
thickness)
{
double efficiency;
double nHcH = numberHydrogen * crossSectionHydrogen;
double nCcC = numberCarbon * crossSectionCarbon;

efficiency = (((nHcH) / (nHcH + nCcC)) * (1-exp(-(nHcH
+ nCcC))
* thickness) * 100);
cout.precision(4);
cout << endl;
cout << "Therefore, the efficiency of our "
<< "system is: " << efficiency << "%" << endl;
cout << endl;
return;
}

@@ Line 131: / Line 131: @@
      cout << endl;
      return;
-   }--[[User:Daq|Daq]] 13:05, 5 May 2008 (MDT)

@@ Line 58: / Line 58: @@
      cin >> neutronEnergy;
-     crossSectionHydrogen = ((4.83 / sqrt(neutronEnergy))
+     crossSectionHydrogen = ((4.83e-24 / sqrt(neutronEnergy))
                                - .578e-24);
@@ Line 74: / Line 74: @@
      double numberCarbon;
-     numberCarbon = thickness * CUBICCENTIMETERCARBON;
+     numberCarbon = CUBICCENTIMETERCARBON;
      return numberCarbon;
+   }
@@ Line 88: / Line 88: @@
      double numberHydrogen;
-     numberHydrogen = thickness * CUBICCENTIMETERHYDROGEN;
+     numberHydrogen = CUBICCENTIMETERHYDROGEN;
      return numberHydrogen;
+   }
@@ Line 105: / Line 105: @@
           << "Detection and Measurement by Glenn Knoll, on"
           << " page 535, " << endl
-          << "Figure 15-15b (in square centimeters): ";
+          << "Figure 15-15b (in barns): ";
      cin >> crossSectionCarbon;
      return crossSectionCarbon;
@@ Line 126: / Line 127: @@
      efficiency = (((nHcH) / (nHcH + nCcC)) * (1-exp(-(nHcH
-                                                    + nCcC))
+                                                    + nCcC)
-                                              * thickness) * 100);
+                                              * thickness)));
      cout.precision(4);
      cout << endl;
      cout << "Therefore, the efficiency of our "
-          << "system is: " << efficiency << "%" << endl;
+          << "system is: " << efficiency *100 << "%" << endl;
      cout << endl;
      return;
+   }

@@ Line 4: / Line 4: @@
   double crossSectionHydrogenFunction();
-  double numberCarbonFunction(double thickness, const double CUBICCENTIMETERCARBON);
+  double numberCarbonFunction(double thickness,
-  double numberHydrogenFunction(double thickness, const double CUBICCENTIMETERHYDROGEN);
+                             const double CUBICCENTIMETERCARBON);
   double crossSectionCarbonFunction();
-  void efficiencyFunction(double numberHydrogen, double numberCarbon, double crossSectionHydrogen,
+  void efficiencyFunction(double numberHydrogen, double numberCarbon,
-                         double crossSectionCarbon, double thickness);
+                         double crossSectionHydrogen, double crossSectionCarbon,
   <code>/************************************************************************
@@ Line 41: / Line 44: @@
+  }
-  /*********************************************************************
+  /**********************************************************************
    * Equation for cross section of hydrogen according to
    * equation 15-9 from Radiation Detection and Measurement
@@ Line 61: / Line 64: @@
+   }
-  /*********************************************************************
+  /**********************************************************************
    * Equation for calculating the number of Carbons in
    * the problem per cubic centimeter, dependent on what
@@ Line 75: / Line 78: @@
+   }
-  /********************************************************************
+  /*********************************************************************
    * Equation for calculating the number of Hydrogens in
    * the problem per cubic centimeter, dependent on what
@@ Line 89: / Line 92: @@
+   }
-  /********************************************************************
+  /*********************************************************************
    * Requesting the cross section from the Knoll book, and
    * then calculating it from barns to centimeters.
@@ Line 108: / Line 111: @@
+   }
-  /******************************************************************
+  /*******************************************************************
    * Equation for efficiency according to equation 15-8b
    * from Radiation Detection and Measurement

@@ Line 1: / Line 1: @@
   #include <iostream>
   #include <cmath>
@@ Line 136: / Line 131: @@
      cout << endl;
      return;
+   }--[[User:Daq|Daq]] 13:05, 5 May 2008 (MDT)

@@ Line 1: / Line 1: @@
-**************************************************************************
+ /**************************************************************************
   * Author: Rebecca Bodily
   * Date: 05/02/2008
   * Organization: Idaho State University
   *************************************************************************/
-#include <iostream>
+ double crossSectionHydrogenFunction();
-#include <cmath>
+ double numberCarbonFunction(double thickness, const double CUBICCENTIMETERCARBON);
-using namespace std;
+ double numberHydrogenFunction(double thickness, const double CUBICCENTIMETERHYDROGEN);
-double crossSectionHydrogenFunction();
+ <code>/*************************************************************************
   * This program is to find the
   * efficiency of a neutron detector.
   ************************************************************************/
-double  main()
+ double  main()
+ {</code>
-  double neutronEnergy;
+   double neutronEnergy;
-  double numberHydrogen;
+   double numberHydrogen;
-  double numberCarbon;
+   double numberCarbon;
-  double crossSectionHydrogen;
+   double crossSectionHydrogen;
-  double crossSectionCarbon;
+   double crossSectionCarbon;
-  double thickness;
+   double thickness;
-  const double CUBICCENTIMETERCARBON = 5.23e22;
+   const double CUBICCENTIMETERCARBON = 5.23e22;
-  const double CUBICCENTIMETERHYDROGEN = 4.74e22;
+   const double CUBICCENTIMETERHYDROGEN = 4.74e22;
-  cout << endl;
+   cout << endl;
-  cout << "Enter thickness of detector (in centimeters): ";
+   cout << "Enter thickness of detector (in centimeters): ";
-  cin >> thickness;
+   cin >> thickness;
-  crossSectionHydrogen = crossSectionHydrogenFunction();
+   crossSectionHydrogen = crossSectionHydrogenFunction();
-  numberCarbon = numberCarbonFunction(thickness,
+   numberCarbon = numberCarbonFunction(thickness,
-                                      CUBICCENTIMETERCARBON);
+                                       CUBICCENTIMETERCARBON);
-  numberHydrogen = numberHydrogenFunction(thickness,
+   numberHydrogen = numberHydrogenFunction(thickness,
-                                          CUBICCENTIMETERHYDROGEN);
+                                           CUBICCENTIMETERHYDROGEN);
-  crossSectionCarbon = crossSectionCarbonFunction();
+   crossSectionCarbon = crossSectionCarbonFunction();
-  efficiencyFunction(numberHydrogen, numberCarbon,
+   efficiencyFunction(numberHydrogen, numberCarbon,
-                     crossSectionHydrogen, crossSectionCarbon,
+                      crossSectionHydrogen, crossSectionCarbon,
-                     thickness);
+                      thickness);
-−
+ }
   /**********************************************************************