Since g0 points are not fit, doesn't matter whether they exist or not.

Module number in all plots. Better stats size.
This commit is contained in:
redford_s
2017-12-01 11:35:37 +01:00
parent 4382f98227
commit cfc5c04fb0

View File

@@ -100,7 +100,7 @@ int main(int argc, char* argv[]) {
thisfile->rewind();
sprintf(savename,"plots/M%s/CurrentSource/pixelmask.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/pixelmask_M%s.png", module_str.c_str(), module_str.c_str());
pixelMaskObject->plotPixelMask(pixel_mask,savename);
cout << "after chip mask, n masked pixels is " << pixelMaskObject->getNMasked(pixel_mask) << endl;
@@ -153,21 +153,21 @@ int main(int argc, char* argv[]) {
pedestalsG0->GetYaxis()->SetTitle("Row");
pedestalsG0->GetYaxis()->SetTitleOffset(0.7);
pedestalsG0->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/pedeG0.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/pedeG0_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
pedestalsG1->GetXaxis()->SetTitle("Column");
pedestalsG1->GetYaxis()->SetTitle("Row");
pedestalsG1->GetYaxis()->SetTitleOffset(0.7);
pedestalsG1->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/pedeG1.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/pedeG1_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
pedestalsG2->GetXaxis()->SetTitle("Column");
pedestalsG2->GetYaxis()->SetTitle("Row");
pedestalsG2->GetYaxis()->SetTitleOffset(0.7);
pedestalsG2->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/pedeG2.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/pedeG2_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
pedeRMSG0->GetXaxis()->SetTitle("Column");
@@ -179,7 +179,7 @@ int main(int argc, char* argv[]) {
pedeRMSG0->GetZaxis()->SetRangeUser(0,25);
}
pedeRMSG0->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/pedeRMSG0.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/pedeRMSG0_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
pedeRMSG1->GetXaxis()->SetTitle("Column");
@@ -191,7 +191,7 @@ int main(int argc, char* argv[]) {
pedeRMSG1->GetZaxis()->SetRangeUser(0,15);
}
pedeRMSG1->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/pedeRMSG1.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/pedeRMSG1_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
pedeRMSG2->GetXaxis()->SetTitle("Column");
@@ -203,7 +203,7 @@ int main(int argc, char* argv[]) {
pedeRMSG2->GetZaxis()->SetRangeUser(0,15);
}
pedeRMSG2->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/pedeRMSG2.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/pedeRMSG2_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
TH2F *adcmap_avg_g0 [97];
@@ -343,26 +343,26 @@ int main(int argc, char* argv[]) {
adcmap_sel->GetYaxis()->SetTitle("Row");
adcmap_sel->GetYaxis()->SetTitleOffset(0.7);
adcmap_sel->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/adcmap_%i_sel.png", module_str.c_str(), j-1);
sprintf(savename,"plots/M%s/CurrentSource/adcmap_%i_sel_M%s.png", module_str.c_str(), j-1, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
adcmap_all->GetXaxis()->SetTitle("Column");
adcmap_all->GetYaxis()->SetTitle("Row");
adcmap_all->GetYaxis()->SetTitleOffset(0.7);
adcmap_all->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/adcmap_%i_all.png", module_str.c_str(), j-1);
sprintf(savename,"plots/M%s/CurrentSource/adcmap_%i_all_M%s.png", module_str.c_str(), j-1, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
gainmap_sel->GetXaxis()->SetTitle("Column");
gainmap_sel->GetYaxis()->SetTitle("Row");
gainmap_sel->GetYaxis()->SetTitleOffset(0.7);
gainmap_sel->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/gainmap_%i_sel.png", module_str.c_str(), j-1);
sprintf(savename,"plots/M%s/CurrentSource/gainmap_%i_sel_M%s.png", module_str.c_str(), j-1, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
gainmap_all->GetXaxis()->SetTitle("Column");
gainmap_all->GetYaxis()->SetTitle("Row");
gainmap_all->GetYaxis()->SetTitleOffset(0.7);
gainmap_all->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/gainmap_%i_all.png", module_str.c_str(), j-1);
sprintf(savename,"plots/M%s/CurrentSource/gainmap_%i_all_M%s.png", module_str.c_str(), j-1, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
adcmap_sel->Reset();
@@ -422,14 +422,14 @@ int main(int argc, char* argv[]) {
gain_avg->GetXaxis()->SetTitleOffset(1.1);
gain_avg->GetYaxis()->SetTitleOffset(1.5);
gain_avg->Draw();
sprintf(savename,"plots/M%s/CurrentSource/perStep/gain_avg_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/gain_avg_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
c1->SaveAs((const char *)(savename));
mapcanvas->cd();
gainmap_avg->GetXaxis()->SetTitle("Column");
gainmap_avg->GetYaxis()->SetTitle("Row");
gainmap_avg->GetYaxis()->SetTitleOffset(0.7);
gainmap_avg->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/perStep/gainmap_avg_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/gainmap_avg_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
c1->cd();
@@ -437,7 +437,7 @@ int main(int argc, char* argv[]) {
adc_avg_g0->GetXaxis()->SetTitleOffset(1.1);
adc_avg_g0->GetYaxis()->SetTitleOffset(1.5);
adc_avg_g0->Draw();
sprintf(savename,"plots/M%s/CurrentSource/perStep/adc_avg_g0_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adc_avg_g0_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
c1->SaveAs((const char *)(savename));
mapcanvas->cd();
adcmap_avg_g0[j/640-1]->GetXaxis()->SetTitle("Column");
@@ -445,14 +445,14 @@ int main(int argc, char* argv[]) {
adcmap_avg_g0[j/640-1]->GetYaxis()->SetTitleOffset(0.7);
adcmap_avg_g0[j/640-1]->Draw("colz");
mapcanvas->Update();
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcmap_avg_g0_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcmap_avg_g0_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
adcermap_avg_g0[j/640-1]->GetXaxis()->SetTitle("Column");
adcermap_avg_g0[j/640-1]->GetYaxis()->SetTitle("Row");
adcermap_avg_g0[j/640-1]->GetYaxis()->SetTitleOffset(0.7);
adcermap_avg_g0[j/640-1]->Draw("colz");
mapcanvas->Update();
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcermap_avg_g0_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcermap_avg_g0_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
c1->cd();
@@ -460,7 +460,7 @@ int main(int argc, char* argv[]) {
adc_avg_g1->GetXaxis()->SetTitleOffset(1.1);
adc_avg_g1->GetYaxis()->SetTitleOffset(1.5);
adc_avg_g1->Draw();
sprintf(savename,"plots/M%s/CurrentSource/perStep/adc_avg_g1_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adc_avg_g1_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
c1->SaveAs((const char *)(savename));
mapcanvas->cd();
adcmap_avg_g1[j/640-1]->GetXaxis()->SetTitle("Column");
@@ -468,14 +468,14 @@ int main(int argc, char* argv[]) {
adcmap_avg_g1[j/640-1]->GetYaxis()->SetTitleOffset(0.7);
adcmap_avg_g1[j/640-1]->Draw("colz");
mapcanvas->Update();
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcmap_avg_g1_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcmap_avg_g1_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
adcermap_avg_g1[j/640-1]->GetXaxis()->SetTitle("Column");
adcermap_avg_g1[j/640-1]->GetYaxis()->SetTitle("Row");
adcermap_avg_g1[j/640-1]->GetYaxis()->SetTitleOffset(0.7);
adcermap_avg_g1[j/640-1]->Draw("colz");
mapcanvas->Update();
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcermap_avg_g1_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcermap_avg_g1_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
c1->cd();
@@ -483,20 +483,20 @@ int main(int argc, char* argv[]) {
adc_avg_g2->GetXaxis()->SetTitleOffset(1.1);
adc_avg_g2->GetYaxis()->SetTitleOffset(1.5);
adc_avg_g2->Draw();
sprintf(savename,"plots/M%s/CurrentSource/perStep/adc_avg_g2_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adc_avg_g2_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
c1->SaveAs((const char *)(savename));
mapcanvas->cd();
adcmap_avg_g2[j/640-1]->GetXaxis()->SetTitle("Column");
adcmap_avg_g2[j/640-1]->GetYaxis()->SetTitle("Row");
adcmap_avg_g2[j/640-1]->GetYaxis()->SetTitleOffset(0.7);
adcmap_avg_g2[j/640-1]->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcmap_avg_g2_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcmap_avg_g2_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
adcermap_avg_g2[j/640-1]->GetXaxis()->SetTitle("Column");
adcermap_avg_g2[j/640-1]->GetYaxis()->SetTitle("Row");
adcermap_avg_g2[j/640-1]->GetYaxis()->SetTitleOffset(0.7);
adcermap_avg_g2[j/640-1]->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcermap_avg_g2_%i.png", module_str.c_str(), j/640);
sprintf(savename,"plots/M%s/CurrentSource/perStep/adcermap_avg_g2_%i_M%s.png", module_str.c_str(), j/640, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
gain_avg->Reset();
@@ -722,7 +722,7 @@ int main(int argc, char* argv[]) {
}
}
if (r0_adc.size() > 1 && r1_adc.size() > 1 && r3_adc.size() > 1) {
if (r1_adc.size() > 1 && r3_adc.size() > 1) {
TGraphErrors *grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0]));
TGraphErrors *grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0]));
@@ -896,7 +896,7 @@ int main(int argc, char* argv[]) {
fit2_e->Draw("same");
fit3_e->Draw("same");
mapcanvas->SetLogx();
sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d.png", module_str.c_str(), pixel_type.c_str(), i);
sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
mapcanvas->SetLogx(0);
@@ -925,7 +925,7 @@ int main(int argc, char* argv[]) {
grap1_norm->SetMaximum(100);
grap1_norm->Draw("AP");
flat_g1->Draw("same");
sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g1norm.png", module_str.c_str(), pixel_type.c_str(), i);
sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g1norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
grap2_norm->SetMarkerColor(kRed);
@@ -938,7 +938,7 @@ int main(int argc, char* argv[]) {
grap2_norm->SetMaximum(100);
grap2_norm->Draw("AP");
flat_g2->Draw("same");
sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g2norm.png", module_str.c_str(), pixel_type.c_str(), i);
sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g2norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
mapcanvas->SetLeftMargin(0.1);
@@ -970,46 +970,60 @@ int main(int argc, char* argv[]) {
g1hist->GetXaxis()->SetTitleOffset(1.1);
g1hist->GetYaxis()->SetTitleOffset(1.5);
g1hist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g1hist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1hist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
g2hist->GetXaxis()->SetTitle("G2 gradient");
g2hist->GetXaxis()->SetTitleOffset(1.1);
g2hist->GetYaxis()->SetTitleOffset(1.5);
g2hist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g2hist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g2hist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
g1overg2hist->GetXaxis()->SetTitle("G1 / G2");
g1overg2hist->GetXaxis()->SetTitleOffset(1.1);
g1overg2hist->GetYaxis()->SetTitleOffset(1.5);
g1overg2hist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g1overg2hist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1overg2hist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
g1overg2hist->Fit("gaus");
sprintf(savename,"plots/M%s/CurrentSource/g1overg2hist_fit.png", module_str.c_str());
g1overg2hist->Draw();
c1->Update();
TPaveText *pave2 = new TPaveText(0.6,0.8,0.94,0.94,"blNDC");
pave2->SetBorderSize(0);
pave2->SetFillStyle(0);
pave2->SetTextSize(0.04);
pave2->SetTextAlign(32);
TF1* g1overg2hist_gaus = g1overg2hist->GetFunction("gaus");
sprintf(savename,"Mean %0.4f #pm %0.4f", g1overg2hist_gaus->GetParameter(1), g1overg2hist_gaus->GetParError(1));
pave2->AddText((const char *)(savename));
sprintf(savename,"Sigma %0.4f #pm %0.4f", g1overg2hist_gaus->GetParameter(2), g1overg2hist_gaus->GetParError(2));
pave2->AddText((const char *)(savename));
pave2->Draw();
g1overg2hist->SetStats(kFALSE);
sprintf(savename,"plots/M%s/CurrentSource/g1overg2hist_fit_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
g1overg2erhist->GetXaxis()->SetTitle("G1 / G2 uncert");
g1overg2erhist->GetXaxis()->SetTitleOffset(1.1);
g1overg2erhist->GetYaxis()->SetTitleOffset(1.5);
g1overg2erhist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g1overg2erhist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1overg2erhist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
g1erhist->GetXaxis()->SetTitle("Fitted G1 uncert");
g1erhist->GetXaxis()->SetTitleOffset(1.1);
g1erhist->GetYaxis()->SetTitleOffset(1.5);
g1erhist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g1erhist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1erhist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
g2erhist->GetXaxis()->SetTitle("Fitted G2 uncert");
g2erhist->GetXaxis()->SetTitleOffset(1.1);
g2erhist->GetYaxis()->SetTitleOffset(1.5);
g2erhist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g2erhist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g2erhist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
diffFromLinG1hist->GetXaxis()->SetTitle("Normalised residual G1");
@@ -1017,7 +1031,7 @@ int main(int argc, char* argv[]) {
diffFromLinG1hist->GetXaxis()->SetTitleOffset(1.1);
diffFromLinG1hist->GetYaxis()->SetTitleOffset(1.5);
diffFromLinG1hist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG1hist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG1hist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
diffFromLinG2hist->GetXaxis()->SetTitle("Normalised residual G2");
@@ -1025,7 +1039,7 @@ int main(int argc, char* argv[]) {
diffFromLinG2hist->GetXaxis()->SetTitleOffset(1.1);
diffFromLinG2hist->GetYaxis()->SetTitleOffset(1.5);
diffFromLinG2hist->Draw();
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG2hist.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG2hist_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
cout << "G1 underflow " << diffFromLinG1hist->GetBinContent(0) << endl;
@@ -1064,7 +1078,7 @@ int main(int argc, char* argv[]) {
g1overg2_isNextToDouble->Draw("same");
g1overg2_isBulk->Draw("same");
leg->Draw("same");
sprintf(savename,"plots/M%s/CurrentSource/g1overg2_perType.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1overg2_perType_M%s.png", module_str.c_str(), module_str.c_str());
c1->SaveAs((const char *)(savename));
mapcanvas->cd();
@@ -1080,7 +1094,7 @@ int main(int argc, char* argv[]) {
g1map->Draw("colz");
pave->AddText("G1 [ADU/#mus]");
pave->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g1map.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1map_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
g2map->GetXaxis()->SetTitle("Column");
@@ -1095,7 +1109,7 @@ int main(int argc, char* argv[]) {
pave->Clear();
pave->AddText("G2 [ADU/#mus]");
pave->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g2map.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g2map_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
g1overg2map->GetXaxis()->SetTitle("Column");
@@ -1106,7 +1120,7 @@ int main(int argc, char* argv[]) {
pave->Clear();
pave->AddText("G1/G2");
pave->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g1overg2map.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1overg2map_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
g1overg2ermap->GetXaxis()->SetTitle("Column");
@@ -1117,7 +1131,7 @@ int main(int argc, char* argv[]) {
pave->Clear();
pave->AddText("#sigma_{G1/G2}");
pave->Draw();
sprintf(savename,"plots/M%s/CurrentSource/g1overg2ermap.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1overg2ermap_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
g1ermap->GetXaxis()->SetTitle("Column");
@@ -1128,7 +1142,7 @@ int main(int argc, char* argv[]) {
pave->AddText("#sigma_{G1}");
pave->Draw();
g1ermap->GetZaxis()->SetRangeUser(0,50);
sprintf(savename,"plots/M%s/CurrentSource/g1ermap.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g1ermap_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
g2ermap->GetXaxis()->SetTitle("Column");
@@ -1139,7 +1153,7 @@ int main(int argc, char* argv[]) {
pave->AddText("#sigma_{G2}");
pave->Draw();
g2ermap->GetZaxis()->SetRangeUser(0,5);
sprintf(savename,"plots/M%s/CurrentSource/g2ermap.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/g2ermap_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
diffFromLinG1map->GetXaxis()->SetTitle("Column");
@@ -1147,7 +1161,7 @@ int main(int argc, char* argv[]) {
diffFromLinG1map->GetYaxis()->SetTitleOffset(0.7);
diffFromLinG1map->Draw("colz");
diffFromLinG1map->GetZaxis()->SetRangeUser(-0.04,0.04);
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG1map.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG1map_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
diffFromLinG2map->GetXaxis()->SetTitle("Column");
@@ -1155,7 +1169,7 @@ int main(int argc, char* argv[]) {
diffFromLinG2map->GetYaxis()->SetTitleOffset(0.7);
diffFromLinG2map->Draw("colz");
diffFromLinG2map->GetZaxis()->SetRangeUser(-0.1,0.1);
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG2map.png", module_str.c_str());
sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG2map_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
sprintf(savename,"data/M%s/CS_ratio_M%s.root", module_str.c_str(), module_str.c_str());