diff --git a/BP_analysis.cpp b/BP_analysis.cpp
index a1961b7..1d17ff3 100644
--- a/BP_analysis.cpp
+++ b/BP_analysis.cpp
@@ -659,29 +659,31 @@ int main(int argc, char* argv[]) {
       // get ratio measurements
       if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) {
 
-      	g0overg1_hist->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
-	g0overg1er_hist->Fill(sqrt(pow(fit_g0->GetParError(1)/fit_g0->GetParameter(1),2) + pow(fit_g1->GetParError(1)/fit_g1->GetParameter(1),2)));
+	double this_g0overg1 = fit_g0->GetParameter(1)/fit_g1->GetParameter(1);
+	double this_g0overg1er = sqrt(pow(fit_g0->GetParError(1)/fit_g0->GetParameter(1),2) + pow(fit_g1->GetParError(1)/fit_g1->GetParameter(1),2));
 
-      	g0overg1_map->Fill(i%NC,i/NC,fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
-	g0overg1er_map->Fill(i%NC,i/NC,sqrt(pow(fit_g0->GetParError(1)/fit_g0->GetParameter(1),2) + pow(fit_g1->GetParError(1)/fit_g1->GetParameter(1),2)));
+      	g0overg1_map->Fill(i%NC,i/NC,this_g0overg1);
+      	g0overg1_hist->Fill(this_g0overg1);
+	g0overg1er_map->Fill(i%NC,i/NC,this_g0overg1er);
+	g0overg1er_hist->Fill(this_g0overg1er);
 
       	if (isEdge(i)) {
-      	  g0overg1_hist_isEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_hist_isEdge->Fill(this_g0overg1);
       	}
       	if (isInnerEdge(i)) {
-      	  g0overg1_hist_isInnerEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_hist_isInnerEdge->Fill(this_g0overg1);
       	}
       	if (isDouble(i)) {
-      	  g0overg1_hist_isDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_hist_isDouble->Fill(this_g0overg1);
       	}
       	if (isNextToDouble(i)) {
-      	  g0overg1_hist_isNextToDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_hist_isNextToDouble->Fill(this_g0overg1);
       	}
       	if (isQuad(i)) {
-      	  g0overg1_hist_isQuad->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_hist_isQuad->Fill(this_g0overg1);
       	}
       	if (isBulk(i)) {
-      	  g0overg1_hist_isBulk->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_hist_isBulk->Fill(this_g0overg1);
       	}
       } else {
 	no_g0g1_calib_map->Fill(i%NC,i/NC,1);
diff --git a/CS_analysis.cpp b/CS_analysis.cpp
index e001575..4cbffef 100644
--- a/CS_analysis.cpp
+++ b/CS_analysis.cpp
@@ -1394,28 +1394,31 @@ int main(int argc, char* argv[]) {
       // get ratio measurements
       if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) {
 
-      	g0overg1hist->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
-      	g0overg1erhist->Fill(abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2)));
-      	g0overg1map->Fill(i%NC,i/NC,fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
-      	g0overg1ermap->Fill(i%NC,i/NC,abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2)));
+	double this_g0overg1 = fit_g0->GetParameter(1)/fit_g1->GetParameter(1);
+	double this_g0overg1er = abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2));
+
+      	g0overg1map->Fill(i%NC,i/NC,this_g0overg1);
+      	g0overg1hist->Fill(this_g0overg1);
+      	g0overg1ermap->Fill(i%NC,i/NC,this_g0overg1er);
+      	g0overg1erhist->Fill(this_g0overg1er);
 
       	if (isEdge(i)) {
-      	  g0overg1_isEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_isEdge->Fill(this_g0overg1);
       	}
       	if (isInnerEdge(i)) {
-      	  g0overg1_isInnerEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_isInnerEdge->Fill(this_g0overg1);
       	}
       	if (isDouble(i)) {
-      	  g0overg1_isDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_isDouble->Fill(this_g0overg1);
       	}
       	if (isNextToDouble(i)) {
-      	  g0overg1_isNextToDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_isNextToDouble->Fill(this_g0overg1);
       	}
       	if (isQuad(i)) {
-      	  g0overg1_isQuad->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_isQuad->Fill(this_g0overg1);
       	}
       	if (isBulk(i)) {
-      	  g0overg1_isBulk->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
+      	  g0overg1_isBulk->Fill(this_g0overg1);
       	}
       } else {
 	no_g0g1_calib_map->Fill(i%NC,i/NC,1);
@@ -1423,28 +1426,31 @@ int main(int argc, char* argv[]) {
 
       if (rangemax_g1 > rangemin_g1 && rangemax_g2 > rangemin_g2) {
 
-      	g1overg2hist->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
-      	g1overg2erhist->Fill(abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2)));
-      	g1overg2map->Fill(i%NC,i/NC,fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
-      	g1overg2ermap->Fill(i%NC,i/NC,abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2)));
+	double this_g1overg2 = fit_g1->GetParameter(1)/fit_g2->GetParameter(1);
+	double this_g1overg2er = abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2));
+
+      	g1overg2map->Fill(i%NC,i/NC,this_g1overg2);
+      	g1overg2hist->Fill(this_g1overg2);
+      	g1overg2ermap->Fill(i%NC,i/NC,this_g1overg2er);
+      	g1overg2erhist->Fill(this_g1overg2er);
 
       	if (isEdge(i)) {
-      	  g1overg2_isEdge->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
+      	  g1overg2_isEdge->Fill(this_g1overg2);
       	}
       	if (isInnerEdge(i)) {
-      	  g1overg2_isInnerEdge->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
+      	  g1overg2_isInnerEdge->Fill(this_g1overg2);
       	}
       	if (isDouble(i)) {
-      	  g1overg2_isDouble->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
+      	  g1overg2_isDouble->Fill(this_g1overg2);
       	}
       	if (isNextToDouble(i)) {
-      	  g1overg2_isNextToDouble->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
+      	  g1overg2_isNextToDouble->Fill(this_g1overg2);
       	}
       	if (isQuad(i)) {
-      	  g1overg2_isQuad->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
+      	  g1overg2_isQuad->Fill(this_g1overg2);
       	}
       	if (isBulk(i)) {
-      	  g1overg2_isBulk->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
+      	  g1overg2_isBulk->Fill(this_g1overg2);
       	}
       } else {
 	no_g1g2_calib_map->Fill(i%NC,i/NC,1);
diff --git a/DB_analysis.cpp b/DB_analysis.cpp
index 936b41a..fdbdf7c 100644
--- a/DB_analysis.cpp
+++ b/DB_analysis.cpp
@@ -735,28 +735,31 @@ int main(int argc, char* argv[]) {
       // get ratio measurements
       if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) {
 
-	g0overg1map->Fill(i%NC,i/NC,fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
-	g0overg1hist->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
-	g0overg1ermap->Fill(i%NC,i/NC,abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2)));
-	g0overg1erhist->Fill(abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2)));
+	double this_g0overg1 = fit_g0->GetParameter(1)/fit_g1->GetParameter(1);
+	double this_g0overg1er = abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2));
+
+	g0overg1map->Fill(i%NC,i/NC,this_g0overg1);
+	g0overg1hist->Fill(this_g0overg1);
+	g0overg1ermap->Fill(i%NC,i/NC,this_g0overg1er);
+	g0overg1erhist->Fill(this_g0overg1er);
 
 	if (isEdge(i)) {
-	  g0overg1hist_isEdge->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
+	  g0overg1hist_isEdge->Fill(this_g0overg1);
 	}
 	if (isInnerEdge(i)) {
-	  g0overg1hist_isInnerEdge->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
+	  g0overg1hist_isInnerEdge->Fill(this_g0overg1);
 	}
 	if (isDouble(i)) {
-	  g0overg1hist_isDouble->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
+	  g0overg1hist_isDouble->Fill(this_g0overg1);
 	}
 	if (isNextToDouble(i)) {
-	  g0overg1hist_isNextToDouble->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
+	  g0overg1hist_isNextToDouble->Fill(this_g0overg1);
 	}
 	if (isQuad(i)) {
-	  g0overg1hist_isQuad->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
+	  g0overg1hist_isQuad->Fill(this_g0overg1);
 	}
 	if (isBulk(i)) {
-	  g0overg1hist_isBulk->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
+	  g0overg1hist_isBulk->Fill(this_g0overg1);
 	}
       } else {
 	no_g0g1_calib_map->Fill(i%NC,i/NC,1);