gwm17/Mask
1
0
Fork 0
mirror of https://github.com/gwm17/Mask.git synced 2024-11-22 10:18:50 -05:00
Code Issues Projects Releases Wiki Activity

Compare commits

base: gwm17:29a136a14b1ab80fdad95901e67d4ee8870c6077
Branches Tags
gwm17:main
gwm17:devel
...
compare: gwm17:cf22fb6384de0c5bb1e00bc90b8c39c061e7b853
Branches Tags
gwm17:main
gwm17:devel

2 Commits
29a136a14b ... cf22fb6384

Author SHA1 Message Date
Gordon McCann cf22fb6384 Update gitignore 2023-05-01 09:27:13 -04:00
Gordon McCann 6815d3d2e7 Bugfixes for ThreeStepSystem, add some utility plots to RootPlot 2023-05-01 09:26:46 -04:00
4 changed files with 7 additions and 3 deletions
Show Stats Expand all files Collapse all files

1
.gitignore vendored
View File

@ -18,6 +18,7 @@ Makefile
build/
lib/
bin/
scripts/
###Keep this file###

2
src/Detectors/DetectorApp.cpp
View File

@ -141,5 +141,7 @@ void DetectorApp::Run()
else if(m_fileWriter.Write())
++count;
}
std::cout << std::endl;
}

4
src/Mask/ThreeStepSystem.cpp
View File

@ -56,8 +56,8 @@ namespace Mask {
m_nuclei[3] = CreateNucleus(zr, ar); //residual
m_nuclei[4] = CreateNucleus(step2Params.Z[1], step2Params.A[1]); //breakup1
m_nuclei[5] = CreateNucleus(zb2, ab2); //breakup2
m_nuclei[5] = CreateNucleus(step3Params.Z[1], step3Params.A[1]); //breakup3
m_nuclei[5] = CreateNucleus(zb4, ab4); //breakup4
m_nuclei[6] = CreateNucleus(step3Params.Z[1], step3Params.A[1]); //breakup3
m_nuclei[7] = CreateNucleus(zb4, ab4); //breakup4
m_step1.BindNuclei(&(m_nuclei[0]), &(m_nuclei[1]), &(m_nuclei[2]), &(m_nuclei[3]));
m_step2.BindNuclei(&(m_nuclei[3]), nullptr, &(m_nuclei[4]), &(m_nuclei[5]));

3
src/Plotters/RootPlotter.cpp
View File

@ -92,7 +92,8 @@ void RootPlotter::FillData(const Mask::Nucleus& nuc)
MyFill(th_vs_ph_name, th_vs_ph_title, nuc.vec4.Theta()*s_rad2deg, FullPhi(nuc.vec4.Phi())*s_rad2deg, 2);
MyFill(ex_name, ex_title, 260, -1.0, 25, nuc.GetExcitationEnergy());
MyFill(angdist_name, angdist_title, 20, -1.0, 1.0, std::cos(nuc.thetaCM));
if(nuc.isDetected)
MyFill(hist_ke_th_name, hist_ke_th_title, 180, 0.0, 180.0, 400, 0.0, 20.0, nuc.vec4.Theta()*s_rad2deg, nuc.GetKE());
if(nuc.isDetected && nuc.detectedKE > 0.25)
{
MyFill(ke_vs_th_name_det, ke_vs_th_title_det, nuc.vec4.Theta()*s_rad2deg, nuc.detectedKE, 2);
MyFill(ke_vs_ph_name_det, ke_vs_ph_title_det, FullPhi(nuc.vec4.Phi())*s_rad2deg, nuc.detectedKE, 4);