Skip to content
Snippets Groups Projects
Commit 31049362 authored by Maria Zurek's avatar Maria Zurek
Browse files

Draw clone before fitting

parent 34171c81
No related branches found
No related tags found
1 merge request!153Draft: Resolve "Add energy scan for Barrel Ecal"
...@@ -223,11 +223,12 @@ std::tuple <double, double, double, double, double, double, double, double> extr ...@@ -223,11 +223,12 @@ std::tuple <double, double, double, double, double, double, double, double> extr
hfsamImg->GetYaxis()->SetTitleOffset(1.4); hfsamImg->GetYaxis()->SetTitleOffset(1.4);
hfsamImg->SetLineWidth(2); hfsamImg->SetLineWidth(2);
hfsamImg->SetLineColor(kBlue); hfsamImg->SetLineColor(kBlue);
set_histo_range(hfsamImg);
hfsamImg->DrawClone();
double up_fit = hfsamImg->GetMean() + 5*hfsamImg->GetStdDev(); double up_fit = hfsamImg->GetMean() + 5*hfsamImg->GetStdDev();
double down_fit = hfsamImg->GetMean() - 5*hfsamImg->GetStdDev(); double down_fit = hfsamImg->GetMean() - 5*hfsamImg->GetStdDev();
if(down_fit <=0 ) down_fit = hfsamImg->GetXaxis()->GetBinUpEdge(1); if(down_fit <=0 ) down_fit = hfsamImg->GetXaxis()->GetBinUpEdge(1);
hfsamImg->Fit("gaus", "", "", down_fit, up_fit); hfsamImg->Fit("gaus", "", "", down_fit, up_fit);
set_histo_range(hfsamImg);
TF1 *gausImg = hfsamImg->GetFunction("gaus"); TF1 *gausImg = hfsamImg->GetFunction("gaus");
gausImg->SetLineWidth(2); gausImg->SetLineWidth(2);
gausImg->SetLineColor(kRed); gausImg->SetLineColor(kRed);
...@@ -235,7 +236,6 @@ std::tuple <double, double, double, double, double, double, double, double> extr ...@@ -235,7 +236,6 @@ std::tuple <double, double, double, double, double, double, double, double> extr
double sigmaImg = gausImg->GetParameter(2); double sigmaImg = gausImg->GetParameter(2);
double meanImg_err = gausImg->GetParError(1); double meanImg_err = gausImg->GetParError(1);
double sigmaImg_err = gausImg->GetParError(2); double sigmaImg_err = gausImg->GetParError(2);
hfsamImg->DrawClone();
save_canvas(c8, "fsamClusterImg", E_label, particle_label); save_canvas(c8, "fsamClusterImg", E_label, particle_label);
} }
{ {
...@@ -243,12 +243,13 @@ std::tuple <double, double, double, double, double, double, double, double> extr ...@@ -243,12 +243,13 @@ std::tuple <double, double, double, double, double, double, double, double> extr
hfsamScFi->GetYaxis()->SetTitleOffset(1.4); hfsamScFi->GetYaxis()->SetTitleOffset(1.4);
hfsamScFi->SetLineWidth(2); hfsamScFi->SetLineWidth(2);
hfsamScFi->SetLineColor(kBlue); hfsamScFi->SetLineColor(kBlue);
set_histo_range(hfsamScFi);
hfsamScFi->DrawClone();
hfsamScFi->Fit("gaus", "", "", 0.01, 0.1); hfsamScFi->Fit("gaus", "", "", 0.01, 0.1);
up_fit = hfsamScFi->GetMean() + 5*hfsamScFi->GetStdDev(); up_fit = hfsamScFi->GetMean() + 5*hfsamScFi->GetStdDev();
down_fit = hfsamScFi->GetMean() - 5*hfsamScFi->GetStdDev(); down_fit = hfsamScFi->GetMean() - 5*hfsamScFi->GetStdDev();
if(down_fit <=0 ) down_fit = hfsamScFi->GetXaxis()->GetBinUpEdge(1); if(down_fit <=0 ) down_fit = hfsamScFi->GetXaxis()->GetBinUpEdge(1);
hfsamScFi->Fit("gaus", "", "", down_fit, up_fit); hfsamScFi->Fit("gaus", "", "", down_fit, up_fit);
set_histo_range(hfsamScFi);
TF1 *gaus = hfsamScFi->GetFunction("gaus"); TF1 *gaus = hfsamScFi->GetFunction("gaus");
gaus->SetLineWidth(2); gaus->SetLineWidth(2);
gaus->SetLineColor(kRed); gaus->SetLineColor(kRed);
...@@ -256,7 +257,7 @@ std::tuple <double, double, double, double, double, double, double, double> extr ...@@ -256,7 +257,7 @@ std::tuple <double, double, double, double, double, double, double, double> extr
double sigmaScFi = gaus->GetParameter(2); double sigmaScFi = gaus->GetParameter(2);
double meanScFi_err = gaus->GetParError(1); double meanScFi_err = gaus->GetParError(1);
double sigmaScFi_err = gaus->GetParError(2); double sigmaScFi_err = gaus->GetParError(2);
hfsamScFi->DrawClone();
save_canvas(c7, "fsamClusterScFi", E_label, particle_label); save_canvas(c7, "fsamClusterScFi", E_label, particle_label);
return std::make_tuple(meanScFi, sigmaScFi, meanScFi_err, sigmaScFi_err, meanImg, sigmaImg, meanImg_err, sigmaImg_err); return std::make_tuple(meanScFi, sigmaScFi, meanScFi_err, sigmaScFi_err, meanImg, sigmaImg, meanImg_err, sigmaImg_err);
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment