Skip to content
Snippets Groups Projects
Commit e5bc631c authored by Andreas Redelbach's avatar Andreas Redelbach
Browse files

Formatting of CbmMcbmCheckTimingAlgo.cxx

parent 703107b2
No related branches found
No related tags found
1 merge request!394Add peak fitting to CbmMcbmCheckTimingAlgo
......@@ -437,11 +437,9 @@ 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);
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");
......@@ -450,11 +448,9 @@ void CbmMcbmCheckTimingAlgo::WriteHistos()
}
break;
}
case ECbmModuleId::kTrd: {
if (DetAverageSingle > 0) {
TF1* gs_trd =
new TF1("gs_trd", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fTrdPeakWidthNs, DetPeakPosSingle + 2 * fTrdPeakWidthNs);
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");
......@@ -463,7 +459,6 @@ void CbmMcbmCheckTimingAlgo::WriteHistos()
}
break;
}
case ECbmModuleId::kT0: {
if (DetAverageSingle > 0) {
TF1* gs_tof =
......@@ -476,7 +471,6 @@ void CbmMcbmCheckTimingAlgo::WriteHistos()
}
break;
}
case ECbmModuleId::kTof: {
if (DetAverageSingle > 0) {
TF1* gs_tof =
......@@ -489,11 +483,9 @@ void CbmMcbmCheckTimingAlgo::WriteHistos()
}
break;
}
case ECbmModuleId::kRich: {
if (DetAverageSingle > 0) {
TF1* gs_rich =
new TF1("gs_rich", "gaus(0)+pol0(3)", DetPeakPosSingle - 2 * fRichPeakWidthNs, DetPeakPosSingle + 2 * fRichPeakWidthNs);
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");
......@@ -502,7 +494,6 @@ void CbmMcbmCheckTimingAlgo::WriteHistos()
}
break;
}
case ECbmModuleId::kPsd: {
if (DetAverageSingle > 0) {
TF1* gs_psd =
......@@ -515,13 +506,10 @@ void CbmMcbmCheckTimingAlgo::WriteHistos()
}
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
......
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment