File indexing completed on 2025-06-08 07:53:36
0001 #include "CalibSummary.h"
0002 #include "TFitResult.h"
0003 #include "TFitResultPtr.h"
0004
0005 ClassImp(CalibSummary);
0006
0007 bool CalibSummary::Fill(const TileCalib& tc){
0008 hLGped .Fill(tc.PedestalMeanL);
0009 hLGpedwidth .Fill(tc.PedestalSigL);
0010 hHGped .Fill(tc.PedestalMeanH);
0011 hHGpedwidth .Fill(tc.PedestalSigH);
0012 hLGscale .Fill(tc.ScaleL);
0013 if (tc.LGHGCorr != 0. ) hLGscaleCalc .Fill(tc.ScaleH/tc.LGHGCorr);
0014 else hLGscaleCalc .Fill(-1000);
0015 hLGscalewidth .Fill(tc.ScaleWidthL);
0016 hHGscale .Fill(tc.ScaleH);
0017 hHGscalewidth .Fill(tc.ScaleWidthH);
0018 hHGLGcorr .Fill(tc.HGLGCorr);
0019 hHGLGOffcorr .Fill(tc.HGLGCorrOff);
0020 hLGHGcorr .Fill(tc.LGHGCorr);
0021 hLGHGOffcorr .Fill(tc.LGHGCorrOff);
0022 return true;
0023 }
0024
0025 bool CalibSummary::Write(TFile* f){
0026 f->cd();
0027 hLGped .Write();
0028 hLGpedwidth .Write();
0029 hHGped .Write();
0030 hHGpedwidth .Write();
0031 hLGscale .Write();
0032 hLGscaleCalc .Write();
0033 hLGscalewidth .Write();
0034 hHGscale .Write();
0035 hHGscalewidth .Write();
0036 hHGLGcorr .Write();
0037 hHGLGOffcorr .Write();
0038 hLGHGcorr .Write();
0039 hLGHGOffcorr .Write();
0040 return true;
0041 }
0042
0043 bool CalibSummary::Analyse(){
0044 std::cout << "***********************************************************************************************************************" << std::endl;
0045 std::cout << "Run Nr.: "<< RunNr << "\t total entries: "<< hLGped.GetEntries() << std::endl;
0046 std::cout << Form("\t --> LG pedestal:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0047 hLGped.GetMean(), hLGped.GetRMS(),
0048 hLGped.GetBinContent(0)+hLGped.GetBinContent(hLGped.GetNbinsX()+1) )
0049 << std::endl;
0050 std::cout << Form("\t --> LG pedestal sigma:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0051 hLGpedwidth.GetMean(), hLGpedwidth.GetRMS(),
0052 hLGpedwidth.GetBinContent(0)+hLGpedwidth.GetBinContent(hLGpedwidth.GetNbinsX()+1) )
0053 << std::endl;
0054 std::cout << Form("\t --> HG pedestal:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0055 hHGped.GetMean(), hHGped.GetRMS(),
0056 hHGped.GetBinContent(0)+hHGped.GetBinContent(hHGped.GetNbinsX()+1) )
0057 << std::endl;
0058 std::cout << Form("\t --> HG pedestal sigma:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0059 hHGpedwidth.GetMean(), hHGpedwidth.GetRMS(),
0060 hHGpedwidth.GetBinContent(0)+hHGpedwidth.GetBinContent(hHGpedwidth.GetNbinsX()+1) )
0061 << std::endl;
0062 std::cout << Form("\t --> LG scale:\t\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0063 hLGscale.GetMean(), hLGscale.GetRMS(),
0064 hLGscale.GetBinContent(0)+hLGscale.GetBinContent(hLGscale.GetNbinsX()+1) )
0065 << std::endl;
0066 std::cout << Form("\t --> LG scale calc:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0067 hLGscaleCalc.GetMean(), hLGscaleCalc.GetRMS(),
0068 hLGscaleCalc.GetBinContent(0)+hLGscaleCalc.GetBinContent(hLGscaleCalc.GetNbinsX()+1) )
0069 << std::endl;
0070 std::cout << Form("\t --> LG scale width:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0071 hLGscalewidth.GetMean(), hLGscalewidth.GetRMS(),
0072 hLGscalewidth.GetBinContent(0)+hLGscalewidth.GetBinContent(hLGscalewidth.GetNbinsX()+1) )
0073 << std::endl;
0074 std::cout << Form("\t --> HG scale:\t\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0075 hHGscale.GetMean(), hHGscale.GetRMS(),
0076 hHGscale.GetBinContent(0)+hHGscale.GetBinContent(hHGscale.GetNbinsX()+1) )
0077 << std::endl;
0078 std::cout << Form("\t --> HG scale width:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0079 hHGscalewidth.GetMean(), hHGscalewidth.GetRMS(),
0080 hHGscalewidth.GetBinContent(0)+hHGscalewidth.GetBinContent(hHGscalewidth.GetNbinsX()+1) )
0081 << std::endl;
0082 std::cout << Form("\t --> LG-HG corr:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0083 hLGHGcorr.GetMean(), hLGHGcorr.GetRMS(),
0084 hLGHGcorr.GetBinContent(0)+hLGHGcorr.GetBinContent(hLGHGcorr.GetNbinsX()+1) )
0085 << std::endl;
0086 std::cout << Form("\t --> LG-HG corr offset:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0087 hLGHGOffcorr.GetMean(), hLGHGcorr.GetRMS(),
0088 hLGHGOffcorr.GetBinContent(0)+hLGHGOffcorr.GetBinContent(hLGHGOffcorr.GetNbinsX()+1) )
0089 << std::endl;
0090 std::cout << Form("\t --> HG-LG corr:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0091 hHGLGcorr.GetMean(), hHGLGcorr.GetRMS(),
0092 hHGLGcorr.GetBinContent(0)+hHGLGcorr.GetBinContent(hHGLGcorr.GetNbinsX()+1) )
0093 << std::endl;
0094 std::cout << Form("\t --> HG-LG corr offset:\t mean:\t%3.3f\t\tRMS:\t%3.3f\t\t%.0f\t out of bounds",
0095 hHGLGOffcorr.GetMean(), hHGLGOffcorr.GetRMS(),
0096 hHGLGOffcorr.GetBinContent(0)+hHGLGOffcorr.GetBinContent(hHGLGOffcorr.GetNbinsX()+1) )
0097 << std::endl;
0098 std::cout << "***********************************************************************************************************************" << std::endl;
0099 return true;
0100 }
