This
website provides C++ source code and necessary input data files to simulate the
gamma ray energy spectrum of the radiative thermal neutron capture on
gadolinium-157 (157Gd) according to the "ANNRI-Gd model". The model was developed on the basis of
data that were taken with an enriched 157Gd sample placed inside a
neutron beam at the Material and Life Science Experimental Facility (MLF) of
J-PARC. The resulting gamma rays from thermal neutron captures were detected
with the germanium crystal arrays of the ANNRI detector.
Using
the provided data files, the source code generates the energies and directions
of gamma rays from the deexcitation of 158Gd after the thermal
neutron capture on 157Gd. All generated directions are totally
random, i.e., without angular correlations. The ANNRI-Gd
model bases on separate descriptions for the continuum component and the
discrete peaks in the gamma ray energy spectrum; they are finally combined to
form the total spectrum.
The
strength of the most discrete peaks in the model is tuned to match the measured
data. For the continuum component, the model employs a statistical description
based on a nuclear level density from the Hartree-Fock-Boguliubov
(HFB) method and a set of electric dipole (E1) photon strength functions (PSF),
which follows the Standard Lorentzian (SLO) model. Details on the ANNRI-Gd model and its performance compared to our own measured
data can be found in our publications.
We also include the model for 155Gd as well, which will
be published soon.
Both the source code and
the input data files are open for general use in any research, provided that
our papers are referred to in your publications.
ANNRI-Gd model is described the gamma-ray spectrum of the thermal
neutron capture on Gd. Note that reactions for neutrons of out
of the energy range are not considered and distinguished. For now, the
model only supports Geant4.9.6 version. The procedure for introducing the model
is shown below. Paths are provided assuming to be in the ANNRI-Gd installation
folder.
0. To define five global variables that the mode has.
The
variables are defined in ./exe_ANNRI_BGO.cc
Variable Name |
Note |
MODEL |
The Model Type |
Gd_CAPTURE |
The Type of Gd |
Gd_CASCADE |
The Type of Cascade |
Gd157_ROOTFile |
The Continuum Parameter for 157Gd |
Gd155_ROOTFile |
The Continuum Parameter for 155Gd |
G4int MODEL = 3:Default Model, 4:ANNRI-Gd Model
G4int Gd_CAPTURE = 1:Natural Gd, 2:Enriched
157Gd, 3:Enriched 155Gd
G4int Gd_CASCADE = 1:Discrete+Coninuum, 2:Only Discrete, 3:Only
Continuum
G4String
Gd157_ROOTFile = "Path of ROOT file of continuum parameter for 157Gd(n,
gamma)�h
G4String
Gd155_ROOTFile = "Path of ROOT file of continuum parameter for 155Gd(n,
gamma)�h
Download full package (61.2 MB)
Make file(./)
GNUmakefile
Main source file(./)
exe_ANNRI_BGO.cc
Source files (./src/)
ANNRIGd_156GdContinuumModelV2.cc
ANNRIGd_156GdDiscreteModel.cc
ANNRIGd_158GdContinuumModelV2.cc
ANNRIGd_158GdDiscreteModel.cc
ANNRIGd_Auxiliary.cc
ANNRIGd_DummyModel.cc
ANNRIGd_GdNCaptureGammaGenerator.cc
ANNRIGd_GeneratorConfigurator.cc
ANNRIGd_Model.cc
ANNRIGd_ModelType.cc
ANNRIGd_OutputConverter.cc
ANNRIGd_ReactionProduct.cc
BGOHit.cc
BGOSD.cc
CalHit.cc
CalorimeterSD.cc
DetectorConstruction.cc
DrawMessage.cc
EventAction.cc
GdCaptureGammas_ggarnet.cc
GdCaptureGammas_glg4sim.cc
GdNeutronHPCapture.cc
GdNeutronHPCaptureFS.cc
Physics_EMPhysics.cc
Physics_GeneralPhysics.cc
Physics_HadronPhysics.cc
Physics_MuonPhysics.cc
Physics_PhysicsList.cc
Physics_RadioactiveDecayPhysics.cc
PrimaryGeneratorAction.cc
RunAction.cc
TrackingAction.cc
Include files (./include/)
ANNRIGd_156GdContinuumModelV2.hh
ANNRIGd_156GdDiscreteModel.hh
ANNRIGd_158GdContinuumModelV2.hh
ANNRIGd_158GdDiscreteModel.hh
ANNRIGd_Auxiliary.hh
ANNRIGd_DummyModel.hh
ANNRIGd_GdNCaptureGammaGenerator.hh
ANNRIGd_GeneratorConfigurator.hh
ANNRIGd_Model.hh
ANNRIGd_ModelType.hh
ANNRIGd_OutputConverter.hh
ANNRIGd_Random.hh
ANNRIGd_ReactionProduct.hh
BGOHit.hh
BGOSD.hh
CalHit.hh
CalorimeterSD.hh
DetectorConstruction.hh
DrawMessage.hh
EventAction.hh
G4MultipleScattering.hh
GdCaptureGammas_ggarnet.hh
GdCaptureGammas_glg4sim.hh
GdNeutronHPCapture.hh
GdNeutronHPCaptureFS.hh
Physics_EMPhysics.hh
Physics_GeneralPhysics.hh
Physics_HadronPhysics.hh
Physics_MuonPhysics.hh
Physics_PhysicsList.hh
Physics_RadioactiveDecayPhysics.hh
PrimaryGeneratorAction.hh
RunAction.hh
TrackingAction.hh
Input data files (./cont_dat/)
156GdContTbl_E1SLO4_HFB.root
158GdContTbl_E1SLO4_HFB.root
Gd155.dat
Gd157_org.dat
Input data files (./macro/)
beam.mac
Before
compiling:
Compile:
After
compiling:
Default
output files:
·
Mandeep
Singh Reen (Okayama),
Email: Mandeepsingh@okayama-u.ac.jp
·
Takashi
Sudo (RCNP, Osaka University),
Email: stakasi@rcnp.osaka-u.ac.jp
·
Ajmi Ali (Padova)
Email: ajmi.ali@pd.infn.it
·
Enrico
Gazzola (Padova)
Email: gazzola.lann@gmail.com
·
Makoto
Sakuda (Okayama)
Email: sakuda-m@okayama-u.ac.jp