2009

From New IAC Wiki
Revision as of 17:39, 23 October 2010 by Oborn (talk | contribs) (Created page with '=Documentation= Paste_painting_procedure TGEM_WasteHandling =1/23/09= == image== Image:Generic_SideView_GEMDetector.jpg 1.) Search the web for patent which coat...')
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Documentation

Paste_painting_procedure

TGEM_WasteHandling

1/23/09

image

Generic SideView GEMDetector.jpg

1.) Search the web for patent which coats GEM detector with neutron sensitive materials. I think it is for Thermal neutrons.


Materials of high neutron capture cross section are studied widely, an example is the following patent [[1]]


Most of the high neutron capture cross section materials were measured for different neutron energies. A comprehensive work is published in 2000. The project was supported by Korea Atomic Energy Research Institute and Brookhaven National Laboratory [[2]]

U-235 is the one of the best choices since it has a high neutron fission cross section and a long half life compared to other the other isotopes that may come under choice.

A second choice is Th-238 and U-238 which have fission cross section less than that of U-235 but still good for our experiment.[[3]]



A boron coated GEM foil is being made by the company below http://n-cdt.com/

Another method uses BF3.

We are interested in a fissionable material coated onto the copper foils that is thin enough to allow the fission fragments to escape the foil and ionize the gas in the detector.


2.) Search for companies which use either sputtering or coating technology to apply the above material to caopper PCboards with hole int them such that the material does not fill up the hole. Hole diameter = ?

The material sputter onto the copper would have thickness on the order of Angstroms.

The TGEM PCboard would have a surface area of 10 cm x 10 cm.

3.) Current neutron efficiency plots for several detector

NeutronDetectionEfficiency-vs-Energy He3 Tube PND INL.jpg

Media:NeutronDetectionEfficiency-vs-Energy_Ne-213_BaF.pdf

NeutronDetectionEfficiency-vs-Energy BC-408 Scintillator.jpg

NeutronDetectionEfficiency-vs-Energy GEM BoronCoatedFoil Neutron Efficiency.jpg

1/30/09

1.) Investigate if Thorium Oxide will be a good candidate for the fission chamber. You would use electrolysis to coat a TGEM board.

Ways to make thorium fission chamber

2.) Find reference for THGEM9, this was used to determine optimal THGEM design 2 years ago

GEM-copper plate has dimensions of 3X3cm or 10X10cm with thinkness 5 micrometer copper layers.It has holes with diameter 60-80 micrometer.(From Operation of a triple GEM detector with CsI photocathode in pure BF4)

2/6/09

Thin deposition ThO2 molecular plating


http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6TJM-4S03R9J-C&_user=489297&_rdoc=1&_fmt=&_orig=search&_sort=d&view=c&_acct=C000023739&_version=1&_urlVersion=0&_userid=489297&md5=cd4314aa4ce529c19369900da39fbbe1

2/13/09

Check if the people below can deliver Th coated Al or Cu sheets , coating thickness smaller than 5 micron, the thinner the better. http://www-wnt.gsi.de/tasca/

2/20/09

2.5 x 2.5 cm thorium coated TGEM cards coated by ISU chemistry.

chemistry Department cant not do the sputtering for safety purposes, they like to avoid radiation contamination.

Need a mask for the predrilled TGEM cards to prevent Thorium from entering holes


Radioactive waste procedure if we are allowed to sputter in chemistry.


Write report describing the process we want to do.

ThoriumPlatingTGEMproposal

Radioactiving of Thorium may be a stumbling blog because of waste generated.

If we go for non-radioactive materials look up fission X-sect for Bismuth and Dysprosium

Fission cross section for Bi starts to be effective when the neutron energy is more than the range of interest.

Neutron induced fission cross section for Bi.jpg


Dy-isotopes are very good for absorbing neutrons in the range between 0.01- 10 MeV.

Dy neutron cross sections.jpg

Media:Zaidi_RadChem_vol93_2005.pdf

2/27/09

1.) Dysprosium (Dy) makes a lot of gammas and maybe 100 less alphas


2.) Don't give up trying to make thorium coated materials

try to send e-mail to one of these authors

N. Takahashi, Zeitschrift für Physik A Hadrons and Nuclei Volume 353, Number 1 / March, 1995

3.) Fission cross-section n,f for Tb

3/13/09

1.) look for a company that does resistive evaporative coatings

This company sells the machine http://www.lesker.com/newweb/Deposition_Sources/ThermalEvaporationSources_Resistive.cfm

2.)(n,f) X-sect for Dysprosium (Dy) ?


3.) Thorium and thorium oxide thin films (19 to 61 nm thick) were RF-sputtered onto mirrors. RF sputtering onto copper plates?


The neutron fission cross sections of 92235U and 92238U between 0.3 and 12.5 MeV W W Osterhage et al 1978 J. Phys. G: Nucl. Phys. 4 587-595

http://www.inf.uu.se/Reports/publications.html

Neutron-induced fission cross sections of natPb and 197Au in the 45-180 MeV region, V.P. Eismont, A.V. Prokofiev, A.N. Smirnov, S.M. Soloviev, H. Condé, K. Elmgren, N. Olsson and P.-U. Renberg Conference Proc. ADTTA99, 1999, (in press).

Up-to-date status and problems of the experimental nucleon-induced fission cross section data base at intermediate energies, V.P. Eismont, A.V. Prokofiev, A.N. Smirnov, I.V. Ryzhov, G.A. Tutin, H. Cond, K. Elmgren, N. Olsson and P.-U. Renberg, Proc. ADTTA99, 1999, (in press).

Neutron-induced fission cross section ratios of 209Bi and 238U at 75 and 96 MeV, V.P. Eismont, A.V. Kireev, I.V. Ryzhov, S.M. Soloviev, G.A. Tutin, H. Condé, K. Elmgren, N. Olsson and P.-U. Renberg, Proc. ADTTA99, 1999, (in press).

Neutron-induced fission fragment angular distribution of 238U at 96 MeV, V.P. Eismont, A.V. Kireev, I.V. Ryzhov, S.M. Soloviev, G.A. Tutin, H. Condé, K. Elmgren, N. Olsson and P.-U. Renberg Proc. ADTTA99, 1999, (in press).

Measurements of neutron-induced fission cross sections of heavy nuclei in the intermediate energy region, V.P. Eismont, A.V. Prokofyev, A.N. Smirnov, K. Elmgren, J. Blomgren, H. Condé, J. Nilsson, N. Olsson and E. Ramström, Accelerator-Driven Transmutation Technologies and Applications, Kalmar, ed. H. Condé (Uppsala: Uppsala University, 1997) p 606-612.

4/2/09

Summury of events through the last two weeks:

1- Looking for a neutron fission cross section for Dyspromium.

2- Looking for other possile elements that can be sued beside our choices for Thorium and Dysprosium.
    these elemets should   have a high neutron cross for fission (n,f),(n,gamma),(n,p) or (n,alpha). 
3- Using Neutrons cross sections (by D. Hughes and R. Schwartz, 2nd edition, 1957) as comphensive reference for our choice.
  

Comparison between Thorium and Dysprosium

.

Thorium: chacterized by relatively high neutron fission cross for both of its isotopes (Th-232 nad Th-230) compared to the stable elements but it is one of low fission cross sections compared to radioctive actinides. A lot of efforts are spent even in coating or finding the appropriate place or group to do that, since most of people are totally disencouraged to coat radioctive elements.

Table of different neutron fission cross 11MeV - 13 MeV.jpg Th-232 fission cross section.jpg


Dysprosium: radioactively stable,it has 5 isotopes Dy-160, 161, 162, 163 and 164. It is one of the best elements for detecting the neutrons because of high neutron capture cross section in the energy range of interest.

Total and capture neutron cross sections for dy- 160 until 20 mev.jpg Total and capture neutron cross sections for dy- 161 until 20 mev.jpgTotal and capture neutron cross sections for dy- 162 until 20 mev.jpg Total and capture neutron cross sections for dy- 163 until 20 mev.jpgTotal and capture neutron cross sections for dy- 164 until 20 mev.jpg

Also dysprosium pulls out alpha paticles (σ(n,α)= 3.6 mb at 14.2 MeV).

Other Possible Elements

. Generally,the highest values for neutron cross sections (σ(n,f), σ(n,p), σ(n,γ) and σ(n,α)) are for radioactive isotopes, which are not desirable for coating. Fortunately, there are non- radioactive elements have a relatively high neutron cross sections, for example:

σ(n,γ) at 0.025 MeV (in mb): Gd-158(710(70)), Sm-152 (670(100)), Br-81 (550(55)), Sm-154 (530(70)), Ce-142 (425(45)),

Hf-180(440(70)), Ru-96(320(60)), Ru-102 (390(40)), W-186 (300(40)), Zr-96 (240(40)), Hg-202(57(13))


σ(n,γ) at 1 MeV (in mb): Re-185(180), (390(40)), Ru-102 (30), Ru-104 (31), Hg-204(100), Br-81 (17)


σ(n,γ) at 14.5 MeV (in mb): Gd-160 (19(5)), Zr-96 (<4),

Element X-sect (mb)
Re-185 180
Ru-102 390

4/3/09

1.) Derive equations for voltages at each GEM stage and the net current for our voltage network and check them with measured values for the 4kV version of the voltage network.

2.) Reference [J.C. hadler, Radiation Measurements Vol 43 (2008) pf S334-S336] says

  1. R = 14 \mu m = mean range of fission fragment6s in U_3 O_8
  2. R= 12 um for UO_3
  3. 0.1 um thickness U( Th) will reduce mean length of fission tracks by 2%

3.) Plot N -vs- d using Equation 2 of Reference [J.C. hadler, Radiation Measurements Vol 43 (2008) pf S334 The number of particles per unit area is given by the equation:

 [math]N = 2/d *(1 - d/12.07) + 1/12.07 [/math]

d : the thickness in (um)

The range of the fragments in the emulsion is 12.07 um Fragment vs thickness (Hudler).jpg


4/17/09

A call conference with Dr.Wolfe, As result, the following is considered as next step for the project:

1- Check the melting temperature that a PC boeard can hold ( the size of the chamber can hold the 10 X 10 PC board ).
   
   The experiment is done, the temperature in average was aroubd 305oC, the PC board can hold that temperature easily, 

the PC board was just burned since the medium around is O2, but there is one thing I want ot be sure of, if the PC board has a very light layer of certain material covers its surface and can't hold up that temeperature?

  An experiment is done today (4/30/09) to check the PC board behavior at 350oC.
  A peace of copper sticked by normal used glue.

PC-board 1.jpg

   PC boeard can't hold up the temperature 350oC.

PC board 3.jpg



2- Checking from a vender for chunk bulk ThF4 with size min. 30 cc.
3- Checking the thermal properties of ThF4 specially the melting point and if Thorium is adhere to copper in that temperature.
4- A mask should be prepared from stainless steal (Al is not preferable),inaddtion to thinking of drilling process and who     
    is going to do it.


Laser machining can cut 1" stainless steel sheets

5- Dr. Forest is going to send the email after collecting all the information needed to Dr. Wolfe.

Coaing process is going to be by electron beam, thorium should be heated to 1750 C. (the melting point for pure Th).


Until this point the uniformity of the surface is not an important factor, but I think it would be important whenever we start thinking of the track and the direction for Th fission fragments.

4/23/09

1- Check the melting temperature that a PC board can hold ( the size of the chamber can hold the 10 X 10 PC board ).

Basically FR-4, FR-1, CEM-1 or CEM-3 PCB are made of polytetrafluoroethylene which has (327 °C (620.6 °F))as a melting point.

The experiment is done, the temperature in average was around 305 °C, the PC board can hold that temperature easily,the PC board was just burned since the medium contains O2, but there is one thing I want to be sure of, if the PC board has a very light layer of certain material covers its surface?


2- Checking from a vender for chunk bulk ThF4 with size min. 30 cc.

http://www.element-collection.com/RGB_Elements_OCT04_rev05.pdf

 Element-collection sells thorium with 190$/gram !

3- Will Thorium adhere to copper

Thorium flouride is used for optical purposes, according the technical applications they have a number of recommendations related to sputtering by E-beam, and other materials that helps in thorium flouride adhesion.[4]

Looks there is a little change, Dr. Forest is going to check a vendor for U-238, which has a better fission cross section (3 times higher compared to Th-232),in addition to , the high price for Th-232 and risk of having flakes after sputtering, Throium by itself does not adhere with the surfaces unless other materials are used through the sputtering process (as methioned in the case of coating glass with thorium fluoride).

5/1/09

1.) Results from PCboard heating test: Board melts at 350 but not 310 degrees Celcius (Documents suggest 327 as the melting point.)

Before heating to 350[math]^o C[/math] After Heating to 350[math]^o C[/math]
PC-board 2.jpg PC board 6.jpg
Before heating to 298[math]^o C[/math] After Heating to 298 [math]^o C[/math]
PC board before 298.jpg PC-board at 298.jpg

The experiment is done on the PC-board, it was heated for an hour under 298°C. The result is shown by the second photo, looks that PC-board should be replaced by another material that can hold this temeperature for longer time.

2.) Oak Ridge will give us 1 kg of metal Th-232. We pay for shipping and we need to do rad licenses. Dr. Forest e-mail Dr. Brey and the paper work is beginning.

Contact info

Lloyd J. Jollay Manager Nuclear Technology and Nonproliferation Y-12 National Security Complex P.O. Box 2009 Oak Ridge, TN 37831-8112 Office: 865-241-1872 Fax: 865-574-5169 Pager: 865-873-9146 Mobile: 865-206-9663

3.) Laser cutting can do up to 1" thick stainless steel sheets, we can make a mask! Need to learn CAD to create a drawing of the mask which will be uploaded to the laser cutting machine. make tolerances around 500 micron (1/2 mm).

Insert picture of our GEM PC board with a table of distances.

4.) Meting Point for U-238 =?

Melting Point: 1408 K (1135°C or 2075°F). Boiling Point: 4404 K (4131°C or 7468°F)

5/8/09

1.) It appears the FR-7 melts at 220 celsius so we may not be able to coat copper clad FR4 with U-28 or Th-232. This means we will most likely coat a copper sheet which is attached to a frame for tension and then laminate FR4 after we have coated the copper sheet with U-28 or The-232.

Determine melting point of our current PC boards.

Melting point for the PC board is 260 degrees Celsius,within this degree a change took place in the shape of the board which will make the possibility to use it for sputtering process is very little, the color of the PC-board starts to change at 200 degrees Celsius.

200 [math]^o C[/math] 210 [math]^o C[/math] 260 [math]^o C[/math]
PC 200.jpg PC board at 210 degrees.jpg PC 260.jpg, PC 260 top.jpg

2.) The distance between holes in the TGEM appear to be 1mm. A laser can cut through 1" of stainless steel. What accuracy can the laser have when cutting through the thinnest available SS we want to use for a mask.

 Find laser resolution (200 microns?)

The following company [5] can reach to a tolerance of 10 micrometer if the thickness of the plate is bigger than 300 micrometer, in our case the thickness of copper foil is 45 micrometer, and the hole is diameter is 1 mm, a tolerance of 200 micrometer is going to be good for the drilling the holes and doing the mask. an email from the company will confirm the expectations.

 Determine thinnest SS mask

still waiting for the email.


3.) HV distribution chain calculation. Low voltage version?

Order resistors based on power calculation, work with Tumuna on the order.

4.) Install GEANT4 in your Inca account, Start GEANT4 Fission model (CHIPS or GEISHA).

5/15/09

1.) Insert picture of PCboard at 220 and 260 and denote time spent at that temperature.


2.) The distance between holes in the TGEM appear to be 1mm and are staggered between rows. A laser can cut through 1" of stainless steel. We would like a laser with position accuracy to at least 100 microns (4 mils :1000 mils = 1 in) which can drill a hole whose diameter is accurate to 4 mils. What accuracy can the laser have when cutting through the thinnest available SS we want to use for a mask.

 Find laser resolution (10 microns?)
 Determine thinnest SS mask


3.) HV distribution chain calculation. Low voltage version?

Look at the following diagram :Media:GEM_HV_circuit.ps

[math] I = I_{11} + I_1^{\prime}[/math]


[math]I_{11} \times R_8 - I_1^{\prime} R_2 = 0 [/math]


[math] I_{11} \times R_3 + I_2^{\prime} R_4 = 0[/math]


[math] I_{11}= I_{22} + I_2^{\prime}[/math]


[math] -I_2^{\prime} \times R_4 + I_{22} \times \frac{R_9 \times R_{10}}{ R_9 + R_{10}} = 0 [/math]


[math]I_{22} \times R_{11} - I_3^{\prime} R_{12} = 0 [/math]


[math] I_{22}= I_{33} + I_3{\prime}[/math]


[math]- I_3^{\prime} R_{12} + I_{33} \times \frac{R_{13} \times R_5}{ R_{13} + R_5} = 0 [/math]


[math] V_{in} - I_{11} \times (R_8 +R_3) - I_{22}\times\left( \frac{R_{10} \times R_9}{ R_{10} + R_9} +R_{11} \right) - I_{33} \times\left ( \frac{R_{13} \times R_5}{ R_{13} + R_5} +R_7 \right )= 0 [/math]


GEM HV Dist Net.jpg
Order resistors based on power calculation, work with Tumuna on the order

4.)GEANT4 is installed in Inca account, ExampleN02 was compiled and runs.


Add  GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program.  Goal is to find a model which reproduces the X-section data for Th-232 above.

5/22/09

1.) heat PC board for 3 hours at 210C


2.) Determine laser resolution for a 45 micron thick copper sheet


3.) HV distribution chain calculation. Do a sample calculation for [math]HV_{GEM} = 3200 V[/math] [math] \Delta V_{GEM1} \equiv V_{G1B}-V_{G1T}[/math] =?

4.) Going to order 25 resisters of each Ohm setting. Prepare a HV board for stuffing and assemble after resistors arrive.

5.)Add GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program. Goal is to find a model which reproduces the X-section data for Th-232 above.

6/9/09

1.) Write a brief description ( 2 paragraphs) of the experimental objectives and methods with justification for the specific radionuclides and quantities.

We propose constructing a fission chamber which contains copper PCboards coated with U-238 and Th-232. Our goal is to construct a fast neutron detector. A safe in the LDS will serve as a repository for the U-238 and Th-232 bulk material. A small fraction of the material (100 g) will be sent to a collaborator at another University who will use Electron Beam-Physical Vapor Deposition to coat a 10 cm x 10 cm PCboard with U-28 or Th-232. The radioactive material will be a 5 micron or less coating attached to the PCboard. The PCboard will reside inside a gaseous detector. The quantities of each bulk material are given in the table below. Oak Ridge national Lab will supply 1kg of Th-232 at no cost. Unfortunately, 1kg is the smallest size available.

U-238 Th-232
200 g 1000 g

U-238 has a atomic number 92, and mass number of 238.050785, it is mainly an alpha emitter. Th-232 has a atomic number 90, and mass number of 232.038051 , it is also an alpha emitter.


2) Describe the facility used to store the materials, diagram of room layout. location of safe.

A floorplan of the Laboratory for Detector Science is shown below.

Media:LDS dimension Model(1).pdf

The LDS has 2 safes available for storing radioactive materials. A solid steel safe manufactured by Bunker Hill with the outer dimensions of (25 cm)H X (35 cm) W X (25 cm) D and inner dimensions of 9-1/2” H X 13-1/2” W X 9-5/8” is controlled by a digital lock. The safe itself weighs 15 kg and is mounted to a wall in the LDS. A second iron safe, manufacturer unkown, uses a combination lock and has the outer dimensions of 62 H X57 W X 45 (cm)D.



3.) Describe radiation survey instruments available for monitoring.

We ask that the TSO provide us with all necessary monitoring devices.

4.) Description and estimate of radioactive waste being generated.

We do not expect to generate radioactive was.

6/12/09

1.) heat PC board for 3 hours at 210C


210 [math]^o C[/math]
PC board at 210 degrees.jpg


2.) Determine laser resolution for a 45 micron thick copper sheet

50 micron


3.) HV distribution chain calculation. Do a sample calculation for [math]HV_{GEM} = 3200 V[/math] [math] \Delta V_{GEM1} \equiv V_{G1B}-V_{G1T}[/math] =?

resisters have arrived currently soldering ready to measure HVs (remember, Digital voltmeter has 1 kV max)


4.)Add GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program. Goal is to find a model which reproduces the X-section data for Th-232 above.

6/19/09

1.) Th-232 mask due July 1,2009


2.) HV distribution chain calculation. Do a sample calculation for [math]HV_{GEM} = 3200 V[/math] [math] \Delta V_{GEM1} \equiv V_{G1B}-V_{G1T}[/math] =?

Abdel GEM HV distrib 6-19-09.png



[math] V_{in} - I \times (R_8 +R_3) - I \times \left( \frac{R_{10} \times R_9}{ R_{10} + R_9} +R_{11} \right) - I \times \left ( \frac{R_{13} \times R_5}{ R_{13} + R_5} +R_7 \right )= 0 [/math]



3.)Add GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program. Goal is to find a model which reproduces the X-section data for Th-232 above.

6/30/09

1.) Th-232 mask due July 1,2009

change design to be a series of lines covering the holes. Calculate the % area of the foil to be covered with Th-232.

2.) Insert HV distribution chain measurements in table form and compare to calculation.

Input Voltage = 1000 V

Current measured calculated
(uA) 262.3 220.7



From the figure above:

[math] R_{T11} = R_8 [/math]


[math] R_{T22} = \left( \frac{R_{10} \times R_9}{ R_{10} + R_9} \right )[/math]


[math] R_{T33} = \left( \frac{R_{13} \times R_5}{ R_{13} + R_5} \right )[/math]


3.)Add GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program. Goal is to find a model which reproduces the X-section data for Th-232 above.

Contact F.W. Jones, TRIUMF, 03-DEC-96 and ask for suggestion on simulating Th-232 (n,f) in GEANT4.  Which model to get X-sections right?

F.W. Jones sent by July, 7:

G4LFission is based on the GHEISHA model of fission, by H. Fesefeldt, from GEANT 3.21.

It is a parameterized model using empirical formulas, e.g. the one for average number of neutrons. The photofission part of the model is not implemented.

My role in this was mainly as a code translator, so unfortunately I can't give you an expert opinion on it. I would suggest that for your study, use one of the prepared physics lists which will include the optimal treatments of fission for various applications.

See e.g. http://geant4.cern.ch/support/proc_mod_catalog/physics_lists/referencePL.shtml

If you are in doubt about choosing a physics list, I would strongly recommend that you post a query on the Geant4 forums, describing your application and what you wish to study (indicate if you wish to study nuclear fragments or other yields).

The forum for hadronic physics is here:

http://hypernews.slac.stanford.edu/HyperNews/geant4/get/hadronprocess.html

This forum is monitored by experts who will be able to provide a more definitive answer.

Best regards, Fred Jones

7/10/09

1.) Th-232 mask due July 1,2009

Making the lines as wide as the holes consumes 50% of the surface area.

A 0.45 mm diameter circle covering the holes and connected by 0.3 mm thick lines will increase the area available for fissionable material to 64%.

Decision is to create a series of lines 0.4 mm thick to cover the holes and connect them with support material.

We need a design of lines 0.4 mm thick covering the holes and connecting them. The distance between line should be 0.3 mm and have a tolerance of 0.1 mm. This will be an effective area of 3/7 = 42% for fission material.


2.) Insert HV distribution chain measurements in table form and compare to calculation.

Redo with actual measured resistance values

3.)Add GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program. Goal is to find a model which reproduces the X-section data for Th-232 above.


Install model suggested by Forum below for U-232:


http://hypernews.slac.stanford.edu/HyperNews/geant4/get/hadronprocess.html

add the functions below to the physicslist. these are the fission models

G4HadronFissionProcess* Fission = new G4HadronFissionProcess();

G4ParaFissionModel* FissionModel = new G4ParaFissionModel();
thenFission->RegisterMe(FissionModel);
pManager->AddDiscreteProcess(Fission);

install this into detector construction


// U-235
  a = 235.01*g/mole;
  density = 19.050*g/cm3;
  U235 = new G4Material(name="U235", z= 92., a, density);


4.) Resistive Coating links

http://www.ellsworth.com/conformal.html?tab=Products

7/15/09

Got LLNL fission model with reference below

http://nuclear.llnl.gov/simulation/fission_usermanual.pdf


Created FissLib model in ExN02 using Uranium target.


The command below will turn on tracking printout

/tracking/verbose 1

The command below will send another neutron projectil in /run/beamOn 1

7/24/09

1.) Th-232 mask due July 1,2009

File:Th-232 Mask 7-24-09.pdf

Add alignment holes.

File:Copper foil 04 straight lines with adjustment.pdf


2.) Insert HV distribution chain measurements in table form and compare to calculation.

Redo with actual measured resistance values

Voltage (V) measured calculated
[math] R_{T11} [/math] 165 126.6
[math] R_{T22} [/math] 98.6 103.4
[math] R_{T33} [/math] 81.5 95.9


3.)Add GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program. Goal is to find a model which reproduces the X-section data for Th-232 above.


trying to simulate the experiment by geant4, once when the source is thorium and another when the source is uranium by hitting the target every time with 10^7 neutrons of energy of 7.5 MeV.

The result :

Until now I did not see any effect of fission processes, the thickness of the target is 0.01 cm. then I should make the target thinner.

I read the paper uploaded to the wiki, the editor explained that G4NDL3.10 has limaitions, under geant4 it just has data for U 7 isotopes.

Mostly all what I got in my simulation is gamma, pair production and Compton scattering.

I did search to see who is interested in tracking fission fragment: I found the following contact with "Michael Heffner" in Lawrence Livermore National Lab. who was engaged in tracking hte fission fragment for Cf-239 but unfortunately I did not get a reply for my email:

I am working on tracking the fission fragment of thorium-232 and uranium-238, part of my work is to use Geant4 for experiment simulations, I read about your work in tracking the fission fragment of Pu-239 (n,f). which model did you use for fission fragment tracking?

the link for their work is :

http://74.125.155.132/searchq=cache:SrXOI60kSs0J:www.nscl.msu.edu/~bickley/TPCWorkshop/Presentations/03_Heffner.pdf+tracking+fission+fragment+"Michael+Heffner"&cd=1&hl=en&ct=clnk&gl=us

Use Uranium target
 use 
http://wiki.iac.isu.edu/index.php/Simulations_of_Particle_Interactions_with_Matter#Changing_the_Random_number_seed_in_GEANT4
to save random number generator state for a fission event.
decrease target thickness until you see fission fragments

ExampleN02Detector construction.cc is edited, experiment simulated with 6.5 MeV neutron beam but different target thicknesses, the result is summarized as the following :

7/27/09

Zoom into interaction point to see all the particles (gamma and neutron).  
You will need to do this to see fission fragments anyway.
Beam Energy (MeV) Target thickness (um) the first fission event random number image
6.5 100 7458 N detector 6.5MeV 1mm.jpeg
6.5 10 17475
6.5 2 285078 100 px
6.5 1 2,691,504 N detector 6.5MeV 1um.jpeg

We can track the ions coming out of the fission process by defining the adding G4DiffuseElastic.cc to our physics list.

http://docs.google.com/gview?a=v&q=cache:hsRwLvrsdZMJ:www.oro.doe.gov/riaseb/wrkshop2003/papers/p-2-0-6.pdf+geant4+and+"+tracking+ions"+code&hl=en&gl=us&pli=1

8/3/09

1.) Th-232 mask due July 1,2009

File:Th-232 Mask 7-24-09.pdf

Add alignment holes.

File:Copper foil 04 straight lines with adjustment.pdf


Need to calculate gravitation deflection of stainless steel strip that is 0.4 mm and 11 cm long wide as a function of its thickness.

Plot mask thickness -vs- deflection due to gravity.


2.) Insert HV distribution chain measurements in table form and compare to calculation.

What is uncertainty in the current used to calculate the expected voltage drop.

Voltage (V) measured calculated
[math] R_{T11} [/math] 165 126.6
[math] R_{T22} [/math] 98.6 103.4
[math] R_{T33} [/math] 81.5 95.9


3.)Add GEANT4 Fission model (CHIPS or GEISHA) to ExampleN02 program. Goal is to find a model which reproduces the X-section data for Th-232 above.


Goal
Simulate the fission process using an incident neutron of energy 6.5 MeV and a layer of U-238 using Geant4.


Status
We are now seeing fission events for 1 micron thick sheets of U-238.
We do not see fission fragments yet


Suggestions:

1.) Move target to center of world and zoom in

2.) Add particle definition of ions

3.) add G4DiffuseElastic.cc process for charged pariticles > 0

8/10/09

1.)Need to calculate gravitation deflection of stainless steel strip that is 0.4 mm and 14 cm long wide as a function of its thickness.

Plot mask thickness -vs- deflection due to gravity.

Maximum deflection of plate of a uniform load, is given by the following formula:[6]

[math] \Delta = \frac{ 5 w L^4}{384 E I} [/math] ; [math] I = \frac{b h^3}{12} [/math]

where:

[math]\Delta[/math] : Deflection (m)

w : The weight of the plate per unit length (N/m)

b : The shorter side of the plate (m) (b = 0.4 mm)

E : Young Modulus (N/m^2) (for stainless steel = 200 GN/m^2)

h : The thickness of the plate (m)

L : The length of the plate (m) (d = 0.14 m)

I : The moment of inertia (m^4)

t (m) y_{max} (m)
1.00 10^-7 0.00003150
2.00 *10^-7 0.00000394
3.00 *10^-7 0.00000117
4.00 *10^-7 0.00000049
5.00 *10^-7 0.00000025
1.00 *10^-6 0.000000032
2.00 *10^-6 0.000000004



2.) What is uncertainty in the current used to calculate the expected voltage drop.

measured (Mohm) theoretical value (Mohm) Voltage (V)measured Voltage (V) calculated
[math] R_{T11} [/math] 0.56 [math]\pm[/math] 0.01 0.55 165 [math]\pm[/math] 1 126.6
[math] R_{T22} [/math] 0.47 0.46 98.6 103.4
[math] R_{T33} [/math] 0.42 0.41 81.5 95.9
[math] R_{T33} [/math] 0.56 0.55 244 220
[math] R_{T33} [/math] 0.47 0.46 200 220
[math] R_{T33} [/math] 0.42 0.41 212 220

the uncertainty in measuring the voltage is +_ 10, and for measuring the resistance is +_ 0.05 So the uncertainty in current theoritically is

3.)

Goal
Simulate the fission process using an incident neutron of energy 6.5 MeV and a layer of U-238 using Geant4.


Status
We are now seeing fission events for 1 micron thick sheets of U-2const
We do not see fission fragments yet


Suggestions

1.) Move target to center of world and zoom in

2.) Add particle definition of ions

3.) add G4DiffuseElastic.cc process for charged pariticles > 0

Uranium238 = new G4Isotope(name="U238", iz=92, n=238, a=238.03*g/mole);


G4IonTable *theIonTable =
G4ParticleTable::GetParticleTable()->GetIonTable();
G4RadioactiveDecay* theRadioactiveDecay =
new G4RadioactiveDecay();
for (G4int i=0; i<theIonTable->Entries(); i++)
{
G4String particleName =
theIonTable->GetParticle(i)->GetParticleName();
if (particleName == "GenericIon") {
G4ProcessManager* pmanager =
theIonTable->GetParticle(i)->GetProcessManager();
pmanager->AddProcess(theRadioactiveDecay);
}


// ions
  G4Deuteron::Deuteron();
  G4Triton::Triton();
  G4He3::He3();
  G4Alpha::Alpha();
  G4GenericIon::GenericIonDefinition();

} else if (particleName == "alpha" ||
               particleName == "He3" ||
               particleName == "GenericIon") {

      pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
      G4ionIonisation* ionIoni = new G4ionIonisation();
      ionIoni->SetStepFunction(0.1, 20*um);
      pmanager->AddProcess(ionIoni,                   -1, 2, 2);

    

8/17/09

1.)Convert table to a graph of "Thickness (mils)" on x-axis and "Deflection ([math]\mu[/math] m)" on y-axis. With horizontal line at 20 microns. Any thickness with deflection less than 20 microns may be used.

change units to micron and make y-axis logarithmic

Deflection thickness ss.jpeg


Thermal Expansion Coefficient

In addition to the probability of the mask deflection by its weight, thermal expansion is another effect that may decrease the efficiency of the mask in covering the holes. The thermal expansion coefficient for stainless steel is studied and equals to (um/[math]^o C[/math]): (according to [7])

Stainless steel type um/[math]^o C[/math] (10-6 in/inF)
Steel Stainless Austenitic (304) 17.3 9.6
Steel Stainless Austenitic (310) 14.4 8.0
Steel Stainless Austenitic (316) 16.0 8.9
Steel Stainless Ferritic (410) 9.9 5.5
Copper 16.6 9.3

In our case, the plate expansion toward the width is negligible at 260 [math]^o C[/math] (assuming the room temperature 20 [math]^o C[/math]), but heating stainless steel (304) plate to the same temperature will cause (0.58mm) change in length! we should be careful as we are welding the ends, or it is better if there is a way to clip the edges to avoid the deflection by length expansion.

Areal thermal expansion

At room temperature, a 0.4 mm length square of stainless steel is going to cover a hole in a copper plate of diameter 0.4 mm. But, as soon as both plates expose to (300 [math]^o C[/math]) furnace, the area will increase by the following formula:

[math] \frac{ \delta A}{A} = 2 \alpha (\delta T) [/math]

So the change in area for SS-plate is 7.6 *10^-4 mm^2, but for the copper is 2.3 * 10^-3 mm^2, a wider mask line should be used.i.e 0.4mm SS-plate width will not cover completely the hole at 300 [math]^o C[/math].


2.) What is uncertainty in the current used to calculate the expected voltage drop.

the uncertainty in voltage can be given by the following formula:


Insert expression for [math]I_{11}[/math] and calculate it using the measured resistances below also propagate error of measured resistances

measured (Mohm) Voltage (V)measured I measured (mA)
[math] R_{T11} [/math] 0.56 [math]\pm[/math] 0.01 165.0 [math]\pm[/math] 0.1 375 [math]\pm[/math] 6.7
[math] R_{T22} [/math] 0.47 [math]\pm[/math] 0.01 98.6 [math]\pm[/math] 0.1 253 [math]\pm[/math] 6.7
[math] R_{T33} [/math] 0.42 [math]\pm[/math] 0.01 81.5 [math]\pm[/math] 0.1 259 [math]\pm[/math] 6.7
Why aren't the currents the same?
Measure the voltage and resistance at each junction , compute current and compare

8/19/2009

measured (Mohm) ( [math]\pm[/math] 0.01) Voltage (V)measured ( [math]\pm[/math] 0.1) I (uA)
[math] R_{8} [/math] 0.56 165.0 294 [math]\pm[/math] 5
[math] R_{3} [/math] 1.01 244.0 241 [math]\pm[/math] 2
[math] R_{9} [/math] 0.51 99.0 194 [math]\pm[/math] 4
[math] R_{10} [/math] 5.51 99.0 16 [math]\pm[/math] 0.03
[math] R_{11} [/math] 1.02 198.6 195 [math]\pm[/math] 2
[math] R_{13} [/math] 0.51 81.4 159 [math]\pm[/math] 0.2
[math] R_{5} [/math] 2.32 81.4 33 [math]\pm[/math] 0.15
[math] R_{7} [/math] 1.04 211.9 203 [math]\pm[/math] 2


8/28/2009

1.)Final version of deflection plot

change units to micron and make y-axis logarithmic

2.) Construct table of measured [math]\Delta[/math] V for each GEM foil preamp and the predicted [math]\Delta[/math] V when [math]V_{GEM}[/math] is 1000, 500, and 200 volts (also try 2 and 4 kV). Use the current measure by CAEN power supply as input to calculation. Measure the voltage difference using a voltmeter.

Perhaps our previous inconsistencies were due to bad current measurements?

3.) Dr. Forest installed GEANT 4.9.2 on Inca. But still no ion tracking.

9/8/2009

1.)Final version of deflection plot

Insert caption and increase font size of labels.  Try bold lines.  Make 20 micron horizontal line dashed.  Make theory curve line black and thicker.

Deflection vs thickness 20umline.jpeg


Insert paragraph describing the results in the curve


2.) Construct table of measured [math]\Delta[/math] V for each GEM foil preamp and the predicted [math]\Delta[/math] V when [math]V_{GEM}[/math] is 1000, 500, and 200 volts (also try 2 and 4 kV). Use the current measure by CAEN power supply as input to calculation. Measure the voltage difference using a voltmeter.

Perhaps our previous inconsistencies were due to bad current measurements?

Measured 200V [math] V_{expected} [/math] 500V [math] V_{expected} [/math] 1000V [math] V_{expected} [/math] 2000V [math] V_{expected} [/math] 4000V [math] V_{expected} [/math]
[math] V_{T11} [/math] (V) [math]\pm[/math] 0.1 24.3 24 [math]\pm[/math]3 59.8 61 [math]\pm[/math]6 118.6 123 [math]\pm[/math]12 238.0 251.4 [math]\pm[/math]23 481.0 530.9 [math]\pm[/math]48
[math] V_{T22} [/math] (V) [math]\pm[/math] 0.1 20.3 20 [math]\pm[/math]3 50.2 51 [math]\pm[/math]6 99.9 103 [math]\pm[/math]11 201.9 211 [math]\pm[/math]23 421.0 445 [math]\pm[/math]48
[math] V_{T33} [/math] (V) [math]\pm[/math] 0.1 18.4 18 [math]\pm[/math]3 45.3 46 [math]\pm[/math]6 90.0 92.8 182.0 188 [math]\pm[/math]22 381.0 398.2 [math]\pm[/math]48
[math] I_{tot} [/math] (uA) [math]\pm[/math] 1 43 44 110 111 221 211 449 443 948 885 [math]\pm[/math]60*
  • considering resistance error accumulates for 6 resistances connected in series (each one has 5%).

The measurements show that we have better accuracy in measuring the current and the voltage within the expected range of error.

Now do a table of Vin , and [math]\Delta[/math] V Gem 1 ,2,3

3.) Dr. Forest installed GEANT 4.9.2 on Inca. But still no ion tracking.

Koi, Tatsumi Phone: (650) 926-4816 E-Mail: tkoi@SLAC.Stanford.EDU

suggests adding the physics list below


source/physics_lists/lists/src/HadronPhysicsQGSP_BIC.cc

9/10/09 Ion Tracking

You can track CF-232 if you add the line

G4GenericIon::GenericIonDefinition();

to the PhysicsList baryon contructor function as shown below

void ExN02PhysicsList::ConstructBaryons()
{
  //  barions                                                                     
  G4Proton::ProtonDefinition();
  G4AntiProton::AntiProtonDefinition();
  G4GenericIon::GenericIonDefinition();
  G4Neutron::NeutronDefinition();
  G4AntiNeutron::AntiNeutronDefinition();
}

9/15/09

1.)Final version of deflection plot

Insert caption and increase font size of labels.  Try bold lines.  Make 20 micron horizontal line dashed.  Make theory curve line black and thicker.

Deflection vs thickness 20umline.jpeg

Insert paragraph describing the results in the curve.

The curve represents the relationship between the thickness and the deflection for a plate of stainless steel of width 0.4 mm and length 140mm when the weight is the only force that acts on it. Please look at [8] for more details about the factors that affects the value of the deflection. The dashed line in the graph intersects with the with deflection line to determine the needed thickness for the plate to have a deflection of 20 um.

Now do a table of Vin , an, and [math]\Delta[/math] V Gem 1 ,2,3

Measured 200V 500V 1000V 2000V 4000V ?
[math] \Delta V_{T11} [/math] (V) [math]\pm[/math] 0.1 11.9 29.8 59.1 119.2 241.7 3000
[math] \Delta V_{T22} [/math] (V) [math]\pm[/math] 0.1 10.2 25.2 50.8 103.2 215.5 2800
[math] \Delta V_{T33} [/math] (V) [math]\pm[/math] 0.1 9.21 23.1 46.1 93.2 195.5 2500

[math] \Delta V_{T33} [/math] is the potential difference between the protection resistance (10 Mohm) and the resistance which is mainly responsible for providing the voltage for the THGEM foil. Please look at the figure [9]. For instance, [math] \Delta V_{T33} [/math] is the voltage between LP3 and LP4.

9/22/09

1.) deflections

According to the figure I contacted a prof. in civilian engineering, He recommended another formula looks simpler and appropriate for a beam fixed from both sides. I am working on the figure and the information that we may extract from formula.

Enter formula with reference

I also talked to him about anything we should consider affects the mask job when it is heated to 300 degrees Celsius, He suggested to buy a Stainless steel wire (beam) that has the same dimension to the one that will cover a row of holes, then heat it and check what is going to happen to the beam. In other words, build a one beam mask and study what will happen to it (we are considering only weight deflection!)


2.) TGEM HV

Applied Potential (V) [math] V_{T11} [/math] (V) [math] \Delta V_{T22} [/math] (V) [math] V_{T22} [/math] (V) [math] \Delta V_{T22} [/math] (V) [math] V_{T33} [/math] (V) [math] \Delta V_{T33} [/math] (V)
200 24.3 11.9 20.3 10.2 18.4 9.2
500 59.8 29.8 50.2 25.2 45.3 23.1
1000 118.6 59.1 103.0 50.8 90.0 46.1
2000 238.0 119.2 211.0 103.2 182.0 93.2
4000 481.0 241.7 421.0 215.5 381.0 195.5

According to the measurements taken, we conclude :

    [math] \ V_{T} [/math] = (approximately) 2 [math] \Delta V_{T} [/math]

So if [math] \Delta V_{T} [/math] is needed to be as the the table below then the value for [math] \ V_{T} [/math] can be determined depending on the previous conclusion:


THGEM foil resistance (Mohm) [math] \Delta V_{T} [/math](V) [math] V_{T} [/math] (V) I (mA) P (W)
[math] R_{11} [/math] = 0.56 3000 6000 11 64
[math] R_{22} [/math] = 0.46 2800 5600 12 66
[math] R_{33} [/math] = 0.42 2500 5000 12 61

From the table, the power supply potential should be 16,600 V.


Study the 12 Volt power system in the links below and determine how to copy


3.) Resistive paste I called the paste company [10], there were in process to make the paste in their labs (because of that it took them that long), looks they succeeded, but until now there is not any information about the shipping which will probably be available on Friday. So I will call (bug) them again on that day.

The paste was shipped on Friday, we are expecting to have it by Wednesday Sept. 30, the company should call me again for tracking number! but I will call them tomorrow.


9/29/09

Paste

The paste was shipped on Friday, we are expecting to have it by Wednesday Sept. 30, the company should call me again for tracking number! but I looks I should call again!


Deflection vs Thickness

Depending on the equation [11] :

Deflection vs thickness 20umline.jpeg


HV

M. Cortesi used a CAEN N471A to individually supply each THGEM stage. The HV supply was set to trip at 50 nA and a 15-20 M Ohm resister was attached in series with the Power supply to limit discharge currents.

Circuit:

The paper gives an idea to how to use the active divider to avoid using more that high voltage power supply, the active divider circuit.

There is a little problem in using that circuit: the maximum voltage provided to the THGEM foil is maximum 1800 V, and that will create a current 22 mA.

Most of the experiments, they starts to run separate HV-power supply connected to a series of high resistances (15- 20 Mohm).

10/6/09

1.) Resistive Paste:

Describe past mixture procedure and curing.

technical data sheetFile:ED7100 Series paste.pdf

safety data sheet File:ED7100 OSHA MSDS.pdf

2.) Change units on deflection plot and insert new plot

3.) Insert link for resistor to use in HV network

references (HV circuit)




4.) Simulation

a.) Zoom in to see target and fission fragments

b.) implement Fission Model G4BertiniEvaporation in Physics list

10/7/09

(n,f) event

I added the following to the physics list

The header files are put into ExN02PhysicsList.cc

#include "G4HENeutronInelastic.hh"
#include "G4NeutronInelasticProcess.hh"
#include "G4CascadeInterface.hh"
#include "G4BertiniEvaporation.hh"

   else if (particleName == "neutron") {
 
      G4CascadeInterface* bertiniModel = new G4CascadeInterface();
      G4NeutronInelasticProcess* inelProcess = new G4NeutronInelasticProcess();
      inelProcess->RegisterMe(bertiniModel);
      pmanager->AddDiscreteProcess(inelProcess);
}


And I saw the event



*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    6.5 MeV     0 eV      0 fm      0 fm       Target    initStep
    1      0 fm      0 fm    -10 cm      0 eV      0 eV   14.3 nm   14.3 nm       Target  NeutronInelastic

*********************************************************************************************************
* G4Track Information:   Particle = U238[0.0],   Track ID = 8,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   25.7 keV     0 eV      0 fm      0 fm       Target    initStep
    1 -2.14e+03 fm -7.54e+03 fm    -10 cm      0 eV   25.7 keV   101 nm    101 nm       Target       hIoni

By oct. 8 the follwing event is tracked by the GEANT4 with U-238 target with 7 MeV neutron beam:

:*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm      7 MeV     0 eV      0 fm      0 fm       Target    initStep
    1      0 fm      0 fm    -10 cm      7 MeV     0 eV      5 um      5 um       Target  Transportation
    2      0 fm      0 fm  -9.85 cm      0 eV      0 eV   1.51 mm   1.51 mm      Tracker  NeutronInelastic

*********************************************************************************************************
* G4Track Information:   Particle = N14[0.0],   Track ID = 5,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -9.85 cm   1.45 MeV     0 eV      0 fm      0 fm      Tracker    initStep
    1  -63.2 um   33.9 um  -9.69 cm    177 keV  1.27 MeV  1.56 mm   1.56 mm      Tracker       hIoni
    2   -181 um     63 um  -9.59 cm      0 eV    177 keV  1.13 mm   2.69 mm      Tracker       hIoni

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 4,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -9.85 cm    382 keV     0 eV      0 fm      0 fm      Tracker    initStep
    1  -2.37 mm  -1.45 mm    -10 cm    382 keV     0 eV   3.16 mm   3.16 mm      Tracker  Transportation
    2  -2.38 mm  -1.45 mm    -10 cm    382 keV     0 eV   10.5 um   3.17 mm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -9.85 cm   1.66 MeV     0 eV      0 fm      0 fm      Tracker    initStep
    1     10 cm   1.17 cm  -2.11 cm   1.66 MeV     0 eV   12.7 cm   12.7 cm      Tracker  Transportation
    2     10 cm   1.17 cm  -2.11 cm   1.66 MeV     0 eV   6.35 um   12.7 cm   OutOfWorld  Transportation

I thick we still need to improve the stepping verbose, since whenever the beam is run, there is a limited number of events are detected and usually are the neutrons that just pass, I run the beam only with 100 neutrons only to track the useful ones.

By oct. 12 , fission event is detected after I turned off the inelastic process code using U-238 as target , the result :


### Run 2 start.

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    6.5 MeV     0 eV      0 fm      0 fm       Target    initStep
    1      0 fm      0 fm    -10 cm      0 eV      0 eV    131 nm    131 nm       Target    nFission

*********************************************************************************************************
* G4Track Information:   Particle = U238[0.0],   Track ID = 8,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   25.7 keV     0 eV      0 fm      0 fm       Target    initStep
    1 -2.14e+03 fm -7.54e+03 fm    -10 cm      0 eV   25.7 keV   101 nm    101 nm       Target       hIoni

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 7,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   96.7 keV     0 eV      0 fm      0 fm       Target    initStep
    1   -759 nm   -940 nm    -10 cm   96.7 keV     0 eV   3.58 um   3.58 um       Target  Transportation
    2   -4.5 cm  -5.58 cm     10 cm   96.7 keV     0 eV   21.2 cm   21.2 cm      Tracker  Transportation
    3   -4.5 cm  -5.58 cm     10 cm   96.7 keV     0 eV   3.72 um   21.2 cm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 6,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    343 keV     0 eV      0 fm      0 fm       Target    initStep
    1   1.49 um    3.5 um    -10 cm    343 keV     0 eV   4.07 um   4.07 um       Target  Transportation
    2   3.46 um   8.14 um    -10 cm    343 keV     0 eV    5.4 um   9.47 um        World  Transportation
    3   4.25 cm     10 cm  -5.86 cm    343 keV     0 eV   11.6 cm   11.6 cm      Tracker  Transportation
    4   4.25 cm     10 cm  -5.86 cm    343 keV     0 eV   4.07 um   11.6 cm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 5,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    597 keV     0 eV      0 fm      0 fm       Target    initStep
    1  -2.93 um   -3.5 um    -10 cm    597 keV     0 eV   4.65 um   4.65 um       Target  Transportation
    2  -11.6 um  -13.8 um    -10 cm    597 keV     0 eV   13.7 um   18.4 um        World  Transportation
    3  -8.38 cm    -10 cm  -7.56 cm    597 keV     0 eV   13.3 cm   13.3 cm      Tracker  Transportation
    4  -8.38 cm    -10 cm  -7.56 cm    597 keV     0 eV   4.65 um   13.3 cm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 4,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   1.07 MeV     0 eV      0 fm      0 fm       Target    initStep
    1    210 nm  -1.63 um    -10 cm   1.07 MeV     0 eV   3.99 um   3.99 um   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   1.68 MeV     0 eV      0 fm      0 fm       Target    initStep
    1  -3.06 um    3.5 um    -10 cm   1.68 MeV     0 eV   5.24 um   5.24 um   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 2,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   2.67 MeV     0 eV      0 fm      0 fm       Target    initStep
    1  -1.44 um    3.5 um    -10 cm   2.67 MeV     0 eV   3.84 um   3.84 um   OutOfWorld  Transportation

Fission event is also detected after I turned off the inelastic process code using Th-232 as target , the result :

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    6.5 MeV     0 eV      0 fm      0 fm       Target    initStep
    1      0 fm      0 fm    -10 cm      0 eV      0 eV    400 nm    400 nm       Target    nFission

*********************************************************************************************************
* G4Track Information:   Particle = Th232[0.0],   Track ID = 8,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm     39 keV     0 eV      0 fm      0 fm       Target    initStep
    1 -2.21e+04 fm -2.88e+04 fm    -10 cm      0 eV     39 keV   177 nm    177 nm       Target       hIoni

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 7,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    143 keV     0 eV      0 fm      0 fm       Target    initStep
    1    3.5 um  -92.2 nm    -10 cm    143 keV     0 eV   4.75 um   4.75 um   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 6,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    843 keV     0 eV      0 fm      0 fm       Target    initStep
    1   -1.1 um    165 nm    -10 cm    843 keV     0 eV   3.29 um   3.29 um       Target  Transportation
    2  -7.09 cm   1.07 cm     10 cm    843 keV     0 eV   21.2 cm   21.2 cm      Tracker  Transportation
    3  -7.09 cm   1.07 cm     10 cm    843 keV     0 eV   3.72 um   21.2 cm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 5,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    956 keV     0 eV      0 fm      0 fm       Target    initStep
    1  -1.85 um  -2.71 um    -10 cm    956 keV     0 eV    5.1 um    5.1 um   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 4,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   1.42 MeV     0 eV      0 fm      0 fm       Target    initStep
    1   -3.5 um   -300 nm    -10 cm   1.42 MeV     0 eV   4.41 um   4.41 um   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   1.43 MeV     0 eV      0 fm      0 fm       Target    initStep
    1   2.67 um    3.5 um    -10 cm   1.43 MeV     0 eV   4.46 um   4.46 um   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = gamma,   Track ID = 2,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   1.66 MeV     0 eV      0 fm      0 fm       Target    initStep
    1    393 nm   -3.5 um    -10 cm   1.66 MeV     0 eV   3.54 um   3.54 um       Target  Transportation
    2    494 nm  -4.39 um    -10 cm   1.66 MeV     0 eV    904 nm   4.44 um   OutOfWorld  Transportation
>>> Event 0
    8 trajectories stored in this event.


10/20/09

1.) Redo resistor network below to reflect one in paper. Find resistance required to have 2000 Volt difference between two the top and bottom of a GEM card and have correct power rating.


Thick_GEM_versus_thin_GEM_in_two_phase_argon_avalanche_detectors (HV circuit)[12]

The new HV-circuit is represented by the figure:


Circuit 1
Circuit 2
Circuit 3


Resistor value 0.3[13] 0.62 [14] 1 [15] 1.25 [16] 1.5 [17] 2 [18] 3[19] 5 [20] 10[21] price
[math] R_{1} [/math] 8.6 1 1 1 1 10.13
[math] R_{2} [/math] 4.3 1 2 11.96
[math] R_{3} [/math] 10 1 5.59
[math] R_{4} [/math] 3.3 1 1 1 6.3
[math] R_{5} [/math] 4.8 1 3.95
total 2 1 2 4 2 1



according to the suggestion to add new line of resistances as represented in the last figure, this change the the old circuit properties. the following table summarizes these differences:

2_lines _THGEM_circuit 3l_line_THGEM_circuit
Total Resistance (Mohm) 48.45 41.7
Current (mA) 0.31 0.36
Max. Voltage Per THGEM Foil (V) (1.43Mohm)444.3 (1.23Mohm)442.3
Power (W) 0.14 0.16


2.) reproduce Figure 3 in paper below

http://www.helsinki.fi/~miheikki/system/refs/heikkinen/chep09geant4.pdf

10/27/09 (Paste processing procedure)

Paste

Paste_painting_procedure

HV Circuit

Update the figure for the high voltage circuit.[22]

Deflection

Figure is updated for the deflection figure.[23]

Simulation

Chips interface is used but unfortunately the result still the same, a neutron hits then two of them leave without any fragment tracking.



G4StringChipsInterface* theCascade = new G4StringChipsInterface();
G4HadronFissionProcess* fissionProcess = new G4HadronFissionProcess();
fissionProcess->RegisterMe(theCascade);
pmanager->AddDiscreteProcess(fissionProcess);


*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    6.5 MeV     0 eV      0 fm      0 fm       Target    initStep
    1      0 fm      0 fm    -10 cm      0 eV      0 eV   4.31 um   4.31 um       Target    nFission

*********************************************************************************************************
* G4Track Information:   Particle = U237[0.0],   Track ID = 4,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm   25.2 keV     0 eV      0 fm      0 fm       Target    initStep
    1 -2.35e+03 fm 4.81e+03 fm    -10 cm      0 eV   25.2 keV  99.7 nm   99.7 nm       Target       hIoni

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    165 keV     0 eV      0 fm      0 fm       Target    initStep
    1   1.44 um    -77 nm    -10 cm    165 keV     0 eV    1.6 um    1.6 um       Target  Transportation
    2     10 cm  -5.35 mm  -5.22 cm    165 keV     0 eV   11.1 cm   11.1 cm      Tracker  Transportation
    3     10 cm  -5.35 mm  -5.22 cm    165 keV     0 eV   5.55 um   11.1 cm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 2,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -10 cm    156 keV     0 eV      0 fm      0 fm       Target    initStep
    1   -2.7 um     -5 um    -10 cm    156 keV     0 eV   5.69 um   5.69 um       Target  Transportation
    2  -41.9 um  -77.6 um    -10 cm    156 keV     0 eV   82.6 um   88.3 um   OutOfWorld  Transportation

10/30/09

Resistors for THGEM HV-circuit
5 MOHM [24] and [25] 10 MOHM [26]
quantity 6 12

11/3/09

Paste

need to finish procedure in link below. Insert info from vendor describing reducing the paster for an ink jet printer.

Paste_painting_procedure

Give a copy of paste MSDS to Roy dunker and ask for permission to use their oven located in the LDS hood

HV Circuit

Update the figure for the high voltage circuit.[27]

Arrange with Tamuna some time to use the power supply to check power on the circuit design.

Deflection

Redo deflection figure. [28]

Simulation

reproduce Figure 3 in paper below

File:Abla incl.pdf

File:Abla incl models.pdf

11/11/09

Paste

please follow the link for the update: [29]


HV Circuit

The circuit were tested as was suggested the figure [30] , very low current was passing through the circuit which relatively decreased the dissipated power through the resistors.

Resistor arrangements and values Current(P.S.) (uA)[math]\pm[/math] 1 Current through each line (uA)[math]\pm[/math] 1 (2 lines only) [math] V_{THGEM1} [/math] (V) [math]\pm[/math] 1 [math] V_{THGEM2} [/math] (V)[math]\pm[/math] 1 dissipated power (W)[math]\pm[/math] 0.1%
Third circuit [31] 50 25 105 85 0.025
removing [math] R_{2}[/math] 50 25 136 136 0.006
[math] R_{1}[/math] = 20 Mohm 41 20.5 177 181 0.008
[math] R_{1}[/math] = 30 Mohm (with one protection resistor) 38 19 195 186 0.11
[math] R_{1}[/math] = 30 Mohm (with all figure resistors) 38 19 159 122 0.11

please note that all the previous measurements with [math] HV+[/math] = 1 kV


Result: highly recommended to increase the main resistor(s) which will provide the THGEM foil with needed potential to 30 Mohm which make the current less that 19 uA in each one, but at the same time this will provide each foil with around 200 V when the the power supply is providing 1 kV.So using 15 kV power supply will easily provide each foil with 3000 V (which is our goal).


Simulation
geant4.9.2.p01/examples/extended/radioactivedecay/exrdm/src

Abla model has a little problem to be defined, I tried to use the following statement in ExN02PhysicsList but unfortunately is not working.

    G4InclCascadeInterface* bertiniModel = new G4InclCascadeInterface();
    G4AblaFission* fissionProcess = new G4AblaFission(); 
    fissionProcess->RegisterMe(bertiniModel);
    pmanager->AddDiscreteProcess(fissionProcess);


Depending on a reference used in the same paper File:Abla incl models.pdf, Binary and Bertini models are used for neutron or proton interactions


G4BinaryCascade* bcModel = new G4BinaryCascade();
G4HadronFissionProcess* fissionProcess = new G4HadronFissionProcess();
fissionProcess->RegisterMe(bcModel);
pmanager->AddDiscreteProcess(fissionProcess);


(same as ABLA and the results can be compared). So i added Binary Cascade model to ExN02PhysicsList But unfortunately the result is weird!!! the following error message is not convenient !!

Computed tolerance = 2.0001e-09 mm
Target is 0.001 cm of Carbon
There are 1 chambers in the tracker region. The chambers are 50 mm of ArgonGas
 The distance between chamber is 10 cm
Segmentation fault

12/8/09

      else if (particleName == "proton") { 
      //protons
        G4InclAblaCascadeInterface *theModel  = new G4InclAblaCascadeInterface();
        theModel->SetMinEnergy(0.0 * GeV);
        theModel->SetMaxEnergy(3.0 * GeV);
        G4ProtonInelasticProcess *protonInelasticProcess = new G4ProtonInelasticProcess();
        protonInelasticProcess->AddDataSet(new G4ProtonInelasticCrossSection());
        protonInelasticProcess->RegisterMe(theModel);
        pmanager->AddDiscreteProcess(protonInelasticProcess);
      }
      else if (particleName == "neutron") { 
      //neutrons
        G4InclAblaCascadeInterface *theModel  = new G4InclAblaCascadeInterface();
        theModel->SetMinEnergy(0.0 * GeV);
        theModel->SetMaxEnergy(3.0 * GeV);
        G4NeutronInelasticProcess *neutronInelasticProcess = new G4NeutronInelasticProcess();
        neutronInelasticProcess->AddDataSet(new G4NeutronHPFissionData());
        neutronInelasticProcess->RegisterMe(theModel);
        pmanager->AddDiscreteProcess(neutronInelasticProcess);

Below is a sample event which appears to be the reaction

[math]n + {330 \atop 92 }U_{238} \rightarrow {173 \atop 49 }In_{124} + {152 \atop 43 }Tc_{109} + 6n[/math]


 *********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm    -80 cm     10 MeV     0 eV      0 fm      0 fm        World    initStep
    1      0 fm      0 fm     -5 mm     10 MeV     0 eV   79.5 cm   79.5 cm        World  Transportation
    2      0 fm      0 fm  -1.94 mm      0 eV      0 eV   3.06 mm   79.8 cm       Target  NeutronInelastic

*********************************************************************************************************
* G4Track Information:   Particle = In124[16830.9],   Track ID = 9,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm   84.1 MeV     0 eV      0 fm      0 fm       Target    initStep
    1   -686 um   1.79 mm     -5 mm   84.1 MeV     0 eV   3.61 mm   3.61 mm       Target  Transportation
    2  -17.9 cm   46.7 cm    -80 cm   84.1 MeV     0 eV   93.8 cm   94.2 cm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 8,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm    215 keV     0 eV      0 fm      0 fm       Target    initStep
    1  -3.93 mm     -5 mm  -3.86 mm    215 keV     0 eV   6.64 mm   6.64 mm       Target  Transportation
    2  -62.9 cm    -80 cm  -30.8 cm    215 keV     0 eV   1.06 m    1.06 m    OutOfWorld  Transportation

*********************************************************************************************************
*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 7,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm   2.42 MeV     0 eV      0 fm      0 fm       Target    initStep
    1  -2.86 mm    271 um     -5 mm   2.42 MeV     0 eV    4.2 mm    4.2 mm       Target  Transportation
    2  -74.8 cm   7.07 cm    -80 cm   2.42 MeV     0 eV   1.09 m     1.1 m    OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = Tc109[16830.9],   Track ID = 6,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm   90.1 MeV     0 eV      0 fm      0 fm       Target    initStep
    1  -2.79 mm   1.98 mm     -5 mm   90.1 MeV     0 eV   4.59 mm   4.59 mm       Target  Transportation
    2    -73 cm   51.8 cm    -80 cm   90.1 MeV     0 eV   1.19 m     1.2 m    OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 5,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm   2.85 MeV     0 eV      0 fm      0 fm       Target    initStep
    1  -2.26 mm      5 mm  -3.15 mm   2.85 MeV     0 eV   5.62 mm   5.62 mm       Target  Transportation
    2  -36.2 cm     80 cm  -19.5 cm   2.85 MeV     0 eV   89.4 cm   89.9 cm   OutOfWorld  Transportation
*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 4,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm   1.36 MeV     0 eV      0 fm      0 fm       Target    initStep
    1  -4.91 mm      5 mm   -4.1 mm   1.36 MeV     0 eV   7.33 mm   7.33 mm       Target  Transportation
    2  -78.6 cm     80 cm  -34.6 cm   1.36 MeV     0 eV   1.17 m    1.17 m    OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm   7.03 MeV     0 eV      0 fm      0 fm       Target    initStep
    1   -2.6 mm      5 mm  -2.38 mm   7.03 MeV     0 eV   5.65 mm   5.65 mm       Target  Transportation
    2  -41.6 cm     80 cm  -7.21 cm   7.03 MeV     0 eV   89.9 cm   90.4 cm   OutOfWorld  Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 2,   Parent ID = 1
*********************************************************************************************************

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm      0 fm  -1.94 mm    287 keV     0 eV      0 fm      0 fm       Target    initStep
    1   1.54 mm  -3.74 mm     -5 mm    287 keV     0 eV   5.07 mm   5.07 mm       Target  Transportation
    2     33 cm    -80 cm  -65.6 cm    287 keV     0 eV   1.08 m    1.08 m    OutOfWorld  Transportation
>>> Event 4
    9 trajectories stored in this event.
5 events have been kept for refreshing and/or reviewing.


12/14/09

1.) For class project overlap GEANT4 fission atomic mass distribution with data locate in NIM article: n U-235 measurements : Diiorio NIM vol B147 (1977)487

2.) Create Xmgrace file to reproduce graph below. Upload original measurements and make a table in the wiki.

File:Tingjin IntConfNucDataforScienceTech 2007.pdf

3.) HV measurements


Pb Proton-fission

The mass number is tracked for the fragment instead of he name by

1-Edit ExN02SteppingVerbose.cc in for-loop, replace ParticleName by AtomicMass.

2- Add the the following to for-loopif:


 (fTrack->GetDefinition()->GetPDGEncoding()==2212 &&fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName()=="NeutronInelastic" && fTrack->GetVolume()->GetName() =="Target" && (*fSeco$)  [lp1]->GetDefinition()->GetAtomicMass()>1)
                   outfile
                   << (*fSecondary)[lp1]->GetDefinition()->GetAtomicMass()


n U-235 measurements : Diiorio NIM vol B147 (1977)487

File:Tingjin IntConfNucDataforScienceTech 2007.pdf


Th-232 Fission CrossSection

Th-232 fission cross section.jpg

12/15/2009

Thermal Neutron and U_235 fission
File:Diiorio NIM vol B147 (1977)487.pdf


1) For thermal neutron (0.025 MeV) we can not use the same class mentioned in [32] th Simulation were running for a long time but not event.

2) Different Neutron data sets were used also:

  neutronInelasticProcess->AddDataSet(new G4HadronFissionDataSet()); (does not give any event)
  neutronInelasticProcess->AddDataSet(new G4NeutronHPorLFissionData());(compiles perfectly but GEANT4 crashes after 1st neutron fired)

3) Recently a new presentation done by describes the the thermal neutron interactions, he is replacing the Inelastic process to Elastic process even registered by other classes or used with neutron data sets. File:THermal neutrons interactions.pdf