This is a great help, thank you so much! All is running perfectly.
Forgive my ignorance, but I was wondering how you have made your gamma spectra above look like a spectrum. My histogram for kinetic energy of gamma rays just outputs individual points (below) - is there a way to draw this so it looks like a spectrum like yours do? This is probably a very basic question so sorry for that.
I uncomment the two lines in the .mac file (not needed if you are using the unmodified version of Hadr03.cc from the distribution):
/process/had/particle_hp/use_photo_evaporation true
/process/had/particle_hp/do_not_adjust_final_state true
I output a .ascii file with same name as the .root file (thatβs what the /analysis/h1/setAscii 2 line in the macro does). I then use a separate Octave program (it should be able run under Matlab too) to plot the spectrum. Iβve attached that program as Read_Hadr03Ascii_Spectum_Plot.m.txt. To use it, remove the .txt extension, move it to the same directory as the .ascii file, gzip the .ascii file, change the path to the ascii input file in the .m program to match where it is on your computer. Then run the program under Octave.
Hi Luis,
I have a LibreOffice spreadsheet for 35Cl(n,gamma), 36Ar(n,gamma), 40Ar(n,gamma), 58Ni(n,gamma) comparing Geant4 simulation energy and intensity results to the NNDC data base. However, this forum does not support .ods files for uploading. Please let me know your email address and I will send it to you.
Hi @John_McFee John, you can use manzanilla+lapp_in2p3_fr, with β+β->β@β and β_β β β.β
Could you please include the exact settings for your simulation?
Thanks in advance!