10 #include "TGraphErrors.h" 14 #include "rootlogon.C" 20 TFile*
file =
new TFile(
"/nova/ana/users/kalra/tmva/SystOptimisation/Unfolding/RootFilesFor_UnfoldingMetrics_UseInMakePlotsMacro/ForPi0KE/UnfoldingMetricsEnergy.root");
29 float fCoverage[iIter];
31 float fChiSqErr[iIter];
33 float fBiasErr[iIter];
34 float fCoverageErr[iIter];
37 TGraph* hChiSqCov[iSyst];
38 TGraph* hAverageCor[iSyst];
40 TGraph* hCoverage[iSyst];
43 for(
int i = 0;
i < iSyst;
i++){
45 char coveragename[50];
47 char correlationname[50];
48 sprintf(biasname,
"%s_%i_bias",
"unf_RecopiKE",
i);
49 sprintf(coveragename,
"%s_%i_coverage",
"unf_RecopiKE",
i);
50 sprintf(chi2name,
"%s_%i_chi",
"unf_RecopiKE",
i);
51 sprintf(correlationname,
"%s_%i_correlation",
"unf_RecopiKE",
i);
58 hChiSqCov[
i] = (TGraph*)file->Get(chi2name);
59 hAverageCor[
i] = (TGraph*)file->Get(correlationname);
60 hBias[
i] = (TGraph*)file->Get(biasname);
61 hCoverage[
i] = (TGraph*)file->Get(coveragename);
64 for(
int j = 0;
j < iIter;
j++){
65 TH1F* hChi2Holder =
new TH1F(
"hChi2Holder",
";#chi^{2};Count",100,0,200);
66 TH1F* hCorrel =
new TH1F(
"hCorrel",
";Correlation;Count",100,0,1);
67 TH1F* hBiasHolder =
new TH1F(
"hBiasHolder",
";Bias;Count",100,-10,10);
68 TH1F* hCoverHolder =
new TH1F(
"hCoverHolder",
";Coverage;Count",100,0,1);
70 for(
int i = 0;
i < iSyst;
i++){
74 hChiSqCov[
i]->GetPoint(
j,x,y);
77 hAverageCor[
i]->GetPoint(
j,x,y);
81 hBias[
i]->GetPoint(
j,x,y);
84 hCoverage[
i]->GetPoint(
j,x,y);
85 hCoverHolder->Fill(y);
88 fChiSq[
j] = (
float)hChi2Holder->GetMean();
89 fCor[
j] = (
float)hCorrel->GetMean();
90 fBias[
j] = (
float)hBiasHolder->GetMean();
91 fCoverage[
j] = (
float)hCoverHolder->GetMean();
93 fChiSqErr[
j] = (
float)hChi2Holder->GetRMS();
94 fCorErr[
j] = (
float)hCorrel->GetRMS();
95 fBiasErr[
j] = (
float)hBiasHolder->GetRMS();
96 fCoverageErr[
j] = (
float)hCoverHolder->GetRMS();
100 float xArray[100] = {0};
101 for(
int i = 0;
i < iIter;
i++){
104 float xArrayErr[100] = {0};
106 for(
int i = 0;
i < iIter;
i++){
108 << fCor[
i] <<
"\t" << fBias[
i] <<
"\t" 111 TGraphErrors* hCoverageFinal =
112 new TGraphErrors(100, xArray, fCoverage, xArrayErr, fCoverageErr);
114 TGraphErrors* hChiSqFinal =
115 new TGraphErrors(100, xArray, fChiSq, xArrayErr, fChiSqErr);
117 TGraphErrors* hCorrelationFinal =
118 new TGraphErrors(100, xArray, fCor, xArrayErr, fCorErr);
120 TGraphErrors* hBiasFinal =
121 new TGraphErrors(100, xArray, fBias, xArrayErr, fBiasErr);
125 TH1F* hFrame1 =
new TH1F(
"hFrame1",
";Iterations; #chi^{2}",
128 hFrame1->SetBinContent(
ibin,-100);
131 TH1F* hFrame2 =
new TH1F(
"hFrame2",
";Iterations; Average Correlation",
134 hFrame2->SetBinContent(
ibin,-100);
136 TH1F* hFrame3 =
new TH1F(
"hFrame3",
";Iterations; Average Bias",
139 hFrame3->SetBinContent(
ibin,-100);
141 TH1F* hFrame4 =
new TH1F(
"hFrame4",
";Iterations; Average Coverage",
144 hFrame4->SetBinContent(
ibin,-100);
147 hFrame1->GetXaxis()->CenterTitle();
148 hFrame1->GetYaxis()->CenterTitle();
150 hFrame2->GetXaxis()->CenterTitle();
151 hFrame2->GetYaxis()->CenterTitle();
152 hFrame3->GetXaxis()->CenterTitle();
153 hFrame3->GetYaxis()->CenterTitle();
154 hFrame4->GetXaxis()->CenterTitle();
155 hFrame4->GetYaxis()->CenterTitle();
158 TCanvas*
c1 =
new TCanvas(
"c1",
"c1");
159 hChiSqFinal->SetMarkerStyle(7);
160 hFrame1->GetYaxis()->SetRangeUser(0,110);
162 hChiSqFinal->Draw(
"PL same");
163 c1->SaveAs(
"unf_RecopiKETuned_chisq_averaged.pdf");
165 TCanvas*
c2 =
new TCanvas(
"c2",
"c2");
166 hBiasFinal->SetMarkerStyle(7);
167 hFrame3->GetYaxis()->SetRangeUser(0,0.03);
169 hBiasFinal->Draw(
"PL same");
170 c2->SaveAs(
"unf_RecopiKETuned_bias_averaged.pdf");
172 TCanvas*
c3 =
new TCanvas(
"c3",
"c3");
173 hCorrelationFinal->SetMarkerStyle(7);
174 hFrame2->GetYaxis()->SetRangeUser(0,1.0);
176 hCorrelationFinal->Draw(
"PL same");
177 c3->SaveAs(
"unf_RecopiKETuned_correlation_averaged.pdf");
179 TCanvas*
c4 =
new TCanvas(
"c4",
"c4");
180 hCoverageFinal->SetMarkerStyle(7);
181 hFrame4->GetYaxis()->SetRangeUser(0,1.0);
183 hCoverageFinal->Draw(
"PL same");
184 c4->SaveAs(
"unf_RecopiKETuned_coverage_averaged.pdf");