From e5bc631cdba6a1dfe1648a62a87b413bfcf8f6f8 Mon Sep 17 00:00:00 2001
From: Andreas Redelbach <redelbach@cbmfles01.gsi.de>
Date: Fri, 25 Jun 2021 16:16:42 +0200
Subject: [PATCH] Formatting of CbmMcbmCheckTimingAlgo.cxx
---
.../mcbm2018/tasks/CbmMcbmCheckTimingAlgo.cxx | 148 ++++++++----------
1 file changed, 68 insertions(+), 80 deletions(-)
diff --git a/fles/mcbm2018/tasks/CbmMcbmCheckTimingAlgo.cxx b/fles/mcbm2018/tasks/CbmMcbmCheckTimingAlgo.cxx
index c2104af446..d0ad5c5321 100644
--- a/fles/mcbm2018/tasks/CbmMcbmCheckTimingAlgo.cxx
+++ b/fles/mcbm2018/tasks/CbmMcbmCheckTimingAlgo.cxx
@@ -437,91 +437,79 @@ void CbmMcbmCheckTimingAlgo::WriteHistos()
}
break;
}
-
- case ECbmModuleId::kMuch: {
- if (DetAverageSingle > 0) {
- TF1* gs_much =
- new TF1("gs_much", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fMuchPeakWidthNs, DetPeakPosSingle + 2 * fMuchPeakWidthNs);
- gs_much->SetParameters(DetAverageSingle, DetPeakPosSingle, fMuchPeakWidthNs, DetAverageSingle);
- fvhDetToRefDiff[uDetIdx]->Fit("gs_much", "R");
- TF1* fitresult_much = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_much");
- LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_much->GetParameter(0) << ", "
- << fitresult_much->GetParameter(1) << ", " << fitresult_much->GetParameter(2);
- }
- break;
- }
-
- case ECbmModuleId::kTrd: {
- if (DetAverageSingle > 0) {
- TF1* gs_trd =
- new TF1("gs_trd", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fTrdPeakWidthNs, DetPeakPosSingle + 2 * fTrdPeakWidthNs);
- gs_trd->SetParameters(0.7 * DetAverageSingle, DetPeakPosSingle, fTrdPeakWidthNs, DetAverageSingle);
- fvhDetToRefDiff[uDetIdx]->Fit("gs_trd", "R");
- TF1* fitresult_trd = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_trd");
- LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_trd->GetParameter(0) << ", "
+ case ECbmModuleId::kMuch: {
+ if (DetAverageSingle > 0) {
+ TF1* gs_much = new TF1("gs_much", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fMuchPeakWidthNs, DetPeakPosSingle + 2 * fMuchPeakWidthNs);
+ gs_much->SetParameters(DetAverageSingle, DetPeakPosSingle, fMuchPeakWidthNs, DetAverageSingle);
+ fvhDetToRefDiff[uDetIdx]->Fit("gs_much", "R");
+ TF1* fitresult_much = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_much");
+ LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_much->GetParameter(0) << ", "
+ << fitresult_much->GetParameter(1) << ", " << fitresult_much->GetParameter(2);
+ }
+ break;
+ }
+ case ECbmModuleId::kTrd: {
+ if (DetAverageSingle > 0) {
+ TF1* gs_trd = new TF1("gs_trd", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fTrdPeakWidthNs, DetPeakPosSingle + 2 * fTrdPeakWidthNs);
+ gs_trd->SetParameters(0.7 * DetAverageSingle, DetPeakPosSingle, fTrdPeakWidthNs, DetAverageSingle);
+ fvhDetToRefDiff[uDetIdx]->Fit("gs_trd", "R");
+ TF1* fitresult_trd = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_trd");
+ LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_trd->GetParameter(0) << ", "
<< fitresult_trd->GetParameter(1) << ", " << fitresult_trd->GetParameter(2);
- }
- break;
- }
-
- case ECbmModuleId::kT0: {
- if (DetAverageSingle > 0) {
- TF1* gs_tof =
- new TF1("gs_tof", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fTofPeakWidthNs, DetPeakPosSingle + 2 * fTofPeakWidthNs);
- gs_tof->SetParameters(DetAverageSingle, DetPeakPosSingle, fTofPeakWidthNs, DetAverageSingle);
- fvhDetToRefDiff[uDetIdx]->Fit("gs_tof", "R");
- TF1* fitresult_tof = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_tof");
- LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_tof->GetParameter(0) << ", "
+ }
+ break;
+ }
+ case ECbmModuleId::kT0: {
+ if (DetAverageSingle > 0) {
+ TF1* gs_tof =
+ new TF1("gs_tof", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fTofPeakWidthNs, DetPeakPosSingle + 2 * fTofPeakWidthNs);
+ gs_tof->SetParameters(DetAverageSingle, DetPeakPosSingle, fTofPeakWidthNs, DetAverageSingle);
+ fvhDetToRefDiff[uDetIdx]->Fit("gs_tof", "R");
+ TF1* fitresult_tof = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_tof");
+ LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_tof->GetParameter(0) << ", "
<< fitresult_tof->GetParameter(1) << ", " << fitresult_tof->GetParameter(2);
- }
- break;
- }
-
- case ECbmModuleId::kTof: {
- if (DetAverageSingle > 0) {
- TF1* gs_tof =
- new TF1("gs_tof", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fTofPeakWidthNs, DetPeakPosSingle + 2 * fTofPeakWidthNs);
- gs_tof->SetParameters(DetAverageSingle, DetPeakPosSingle, fTofPeakWidthNs, DetAverageSingle);
- fvhDetToRefDiff[uDetIdx]->Fit("gs_tof", "R");
- TF1* fitresult_tof = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_tof");
- LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_tof->GetParameter(0) << ", "
- << fitresult_tof->GetParameter(1) << ", " << fitresult_tof->GetParameter(2);
- }
- break;
- }
-
- case ECbmModuleId::kRich: {
- if (DetAverageSingle > 0) {
- TF1* gs_rich =
- new TF1("gs_rich", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fRichPeakWidthNs, DetPeakPosSingle + 2 * fRichPeakWidthNs);
- gs_rich->SetParameters(0.5 * DetAverageSingle, DetPeakPosSingle, fRichPeakWidthNs, DetAverageSingle);
- fvhDetToRefDiff[uDetIdx]->Fit("gs_rich", "R");
- TF1* fitresult_rich = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_rich");
- LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_rich->GetParameter(0) << ", "
- << fitresult_rich->GetParameter(1) << ", " << fitresult_rich->GetParameter(2);
- }
- break;
- }
-
- case ECbmModuleId::kPsd: {
- if (DetAverageSingle > 0) {
- TF1* gs_psd =
- new TF1("gs_psd", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fPsdPeakWidthNs, DetPeakPosSingle + 2 * fPsdPeakWidthNs);
- gs_psd->SetParameters(DetAverageSingle, DetPeakPosSingle, fPsdPeakWidthNs, DetAverageSingle);
- fvhDetToRefDiff[uDetIdx]->Fit("gs_psd", "R");
- TF1* fitresult_psd = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_psd");
- LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_psd->GetParameter(0) << ", "
- << fitresult_psd->GetParameter(1) << ", " << fitresult_psd->GetParameter(2);
- }
+ }
break;
+ }
+ case ECbmModuleId::kTof: {
+ if (DetAverageSingle > 0) {
+ TF1* gs_tof =
+ new TF1("gs_tof", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fTofPeakWidthNs, DetPeakPosSingle + 2 * fTofPeakWidthNs);
+ gs_tof->SetParameters(DetAverageSingle, DetPeakPosSingle, fTofPeakWidthNs, DetAverageSingle);
+ fvhDetToRefDiff[uDetIdx]->Fit("gs_tof", "R");
+ TF1* fitresult_tof = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_tof");
+ LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_tof->GetParameter(0) << ", "
+ << fitresult_tof->GetParameter(1) << ", " << fitresult_tof->GetParameter(2);
+ }
+ break;
+ }
+ case ECbmModuleId::kRich: {
+ if (DetAverageSingle > 0) {
+ TF1* gs_rich = new TF1("gs_rich", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fRichPeakWidthNs, DetPeakPosSingle + 2 * fRichPeakWidthNs);
+ gs_rich->SetParameters(0.5 * DetAverageSingle, DetPeakPosSingle, fRichPeakWidthNs, DetAverageSingle);
+ fvhDetToRefDiff[uDetIdx]->Fit("gs_rich", "R");
+ TF1* fitresult_rich = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_rich");
+ LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_rich->GetParameter(0) << ", "
+ << fitresult_rich->GetParameter(1) << ", " << fitresult_rich->GetParameter(2);
}
-
- default: {
- LOG(info)<< "Detector ID for fitting is not valid.";
- break;
+ break;
+ }
+ case ECbmModuleId::kPsd: {
+ if (DetAverageSingle > 0) {
+ TF1* gs_psd =
+ new TF1("gs_psd", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fPsdPeakWidthNs, DetPeakPosSingle + 2 * fPsdPeakWidthNs);
+ gs_psd->SetParameters(DetAverageSingle, DetPeakPosSingle, fPsdPeakWidthNs, DetAverageSingle);
+ fvhDetToRefDiff[uDetIdx]->Fit("gs_psd", "R");
+ TF1* fitresult_psd = fvhDetToRefDiff[uDetIdx]->GetFunction("gs_psd");
+ LOG(debug)<< fvDets[uDetIdx].sName << " parameters from Gauss fit = " << fitresult_psd->GetParameter(0) << ", "
+ << fitresult_psd->GetParameter(1) << ", " << fitresult_psd->GetParameter(2);
}
-
-
+ break;
+ }
+ default: {
+ LOG(info)<< "Detector ID for fitting is not valid.";
+ break;
+ }
}
fvhDetToRefDiff[uDetIdx]->Write(); //At the end in order to include fitting results in histos
--
GitLab