Commit 3d7978c6 authored by Chad Lantz's avatar Chad Lantz
Browse files

Added total nPhotons histogram

parent bda8d19f
...@@ -65,6 +65,8 @@ class RPDAnalysis2021 : public Analysis{ ...@@ -65,6 +65,8 @@ class RPDAnalysis2021 : public Analysis{
TH2D *hPkVsDiffPk; TH2D *hPkVsDiffPk;
/** Charge sum histogram */ /** Charge sum histogram */
TH1D *hChargeSum; TH1D *hChargeSum;
/** Photon count histogram */
TH1D *hPhotons;
/** Peak height sum histogram */ /** Peak height sum histogram */
TH1D *hPeakSum; TH1D *hPeakSum;
/** Differential peak height sum histogram */ /** Differential peak height sum histogram */
......
...@@ -92,6 +92,7 @@ void RPDAnalysis2021::SetupHistograms( ){ ...@@ -92,6 +92,7 @@ void RPDAnalysis2021::SetupHistograms( ){
hChgVsPk = new TH2D("Charge vs Peak height","Q vs Peak; Q; Peak", 300, 0, 300000, 300, 0, 1100); hChgVsPk = new TH2D("Charge vs Peak height","Q vs Peak; Q; Peak", 300, 0, 300000, 300, 0, 1100);
hPkVsDiffPk = new TH2D("Peak height vs Differential peak height","Peak vs Diff Peak", 200, 0, 5000, 200, 0, 300); hPkVsDiffPk = new TH2D("Peak height vs Differential peak height","Peak vs Diff Peak", 200, 0, 5000, 200, 0, 300);
hChargeSum = new TH1D("RPD Integrated Signal","RPD Integrated Signal", 200, 0, 160000); hChargeSum = new TH1D("RPD Integrated Signal","RPD Integrated Signal", 200, 0, 160000);
hPhotons = new TH1D("RPD nPhotons", "Total number of photons", 100, 0, 50);
hPeakSum = new TH1D("RPD Peak Height Sum","RPD Peak Height Sum", 200, 0, 3000); hPeakSum = new TH1D("RPD Peak Height Sum","RPD Peak Height Sum", 200, 0, 3000);
hDiffPeakSum= new TH1D("RPD Differential Peak Height Sum","RPD Diff Peak Sum", 200, 0, 55000); hDiffPeakSum= new TH1D("RPD Differential Peak Height Sum","RPD Diff Peak Sum", 200, 0, 55000);
...@@ -163,6 +164,7 @@ void RPDAnalysis2021::AnalyzeEvent( ){ ...@@ -163,6 +164,7 @@ void RPDAnalysis2021::AnalyzeEvent( ){
PeakSum += rpd[row][col]->Peak_max[hit]; PeakSum += rpd[row][col]->Peak_max[hit];
DiffPeakSum += rpd[row][col]->Diff_max[hit]; DiffPeakSum += rpd[row][col]->Diff_max[hit];
if(rpd[row][col]->gain > 1) hPhotons->Fill(rpd[row][col]->nPhotons[hit]);
hChargeArr[row][col]->Fill(rpd[row][col]->Charge[hit]); hChargeArr[row][col]->Fill(rpd[row][col]->Charge[hit]);
hPeakArr[row][col]->Fill(rpd[row][col]->Peak_max[hit]); hPeakArr[row][col]->Fill(rpd[row][col]->Peak_max[hit]);
hDPeakArr[row][col]->Fill(rpd[row][col]->Diff_max[hit]); hDPeakArr[row][col]->Fill(rpd[row][col]->Diff_max[hit]);
...@@ -437,6 +439,7 @@ void RPDAnalysis2021::Finalize( ){ ...@@ -437,6 +439,7 @@ void RPDAnalysis2021::Finalize( ){
m_viz->DrawPlot(hChgVsPk,"RPD Q_{total}","RPD Peak_{sum}","RPD_TotalCharge.png",""); m_viz->DrawPlot(hChgVsPk,"RPD Q_{total}","RPD Peak_{sum}","RPD_TotalCharge.png","");
m_viz->DrawPlot(hPkVsDiffPk,"RPD Peak_{sum}","#frac{#partial V}{#partial t}_{max}","RPD_TotalCharge.png",""); m_viz->DrawPlot(hPkVsDiffPk,"RPD Peak_{sum}","#frac{#partial V}{#partial t}_{max}","RPD_TotalCharge.png","");
m_viz->DrawPlot(hChargeSum,"RPD Q_{total}","Counts","RPD_TotalCharge.png",""); m_viz->DrawPlot(hChargeSum,"RPD Q_{total}","Counts","RPD_TotalCharge.png","");
m_viz->DrawPlot(hPhotons,"RPD Photon Count","Counts","RPD_nPhotons.png","");
m_viz->DrawPlot(hPeakSum,"RPD Peak_{sum}","Counts","RPD_PeakSum.png",""); m_viz->DrawPlot(hPeakSum,"RPD Peak_{sum}","Counts","RPD_PeakSum.png","");
m_viz->DrawPlot(hDiffPeakSum,"#frac{#partial V}{#partial t}_{max}","Counts","RPD_DiffSum.png",""); m_viz->DrawPlot(hDiffPeakSum,"#frac{#partial V}{#partial t}_{max}","Counts","RPD_DiffSum.png","");
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment