diff --git a/.gitignore b/.gitignore
index c5ad14d..1246fa2 100644
--- a/.gitignore
+++ b/.gitignore
@@ -4,8 +4,8 @@ JFMC_CuFluoPeak
 JFMC_CuFluoPeakFit
 JFMC_CurrentSourceScan
 JFMC_DirectBeamScan
-dummyRootFileMaker
 BP_DB_Comp
+JFMC_ModuleComp
 
 data/
 plots/
diff --git a/JFMC_ModuleComp.cpp b/JFMC_ModuleComp.cpp
new file mode 100644
index 0000000..4ef69cd
--- /dev/null
+++ b/JFMC_ModuleComp.cpp
@@ -0,0 +1,159 @@
+// to extract the calibration constants of all modules from the gain maps and compare
+
+#include "/afs/psi.ch/project/mythen/sophie/sls_detector_calibration/jungfrauCommonHeader.h"
+#include "/afs/psi.ch/project/mythen/sophie/sls_detector_calibration/jungfrauCommonFunctions.h"
+
+#include <fstream> 
+
+#include "TH1F.h"
+#include "TCanvas.h"
+#include "TF1.h"
+#include "TGraphErrors.h"
+
+struct GainMaps
+{
+  double g0VALs[NCH];
+  double g1VALs[NCH];
+  double g2VALs[NCH];
+  double hg0VALs[NCH];
+};
+
+
+int main(int argc, char* argv[]) {
+
+  jungfrauStyle();
+
+  GainMaps *mapsObject = new GainMaps;
+
+  vector<int> modules;
+  modules.push_back(6);
+  modules.push_back(8);
+  modules.push_back(21);
+  modules.push_back(22);
+  modules.push_back(27);
+  modules.push_back(28);
+  modules.push_back(32);
+  modules.push_back(40);
+  modules.push_back(44);
+  modules.push_back(45);
+  modules.push_back(48);
+  modules.push_back(50);
+  modules.push_back(55);
+  modules.push_back(57);
+  modules.push_back(59);
+  modules.push_back(65);
+  modules.push_back(66);
+  modules.push_back(68);
+  modules.push_back(69);
+  modules.push_back(79);
+  modules.push_back(80);
+  modules.push_back(88);
+  modules.push_back(94);
+  modules.push_back(102);
+  modules.push_back(119);
+  modules.push_back(126);
+  
+  vector<double> modules_d;
+  vector<double> module_errs;
+  vector<double> g0_means;
+  vector<double> g0_rmss;
+  vector<double> g1_means;
+  vector<double> g1_rmss;
+  vector<double> g2_means;
+  vector<double> g2_rmss;
+  vector<double> hg0_means;
+  vector<double> hg0_rmss;
+
+  char savename[128];
+
+  TH1F* g0hist = new TH1F("g0hist","",100,30,50);
+  TH1F* g1hist = new TH1F("g1hist","",100,-3,0);
+  TH1F* g2hist = new TH1F("g2hist","",100,-0.2,0);
+  TH1F* hg0hist = new TH1F("hg0hist","",100,80,120);
+
+  for (size_t j = 0; j < modules.size(); j++) {
+
+    // clear gain maps
+    for (int i = 0; i < NCH; i++) {
+      mapsObject->g0VALs[i] = 0;
+      mapsObject->g1VALs[i] = 0;
+      mapsObject->g2VALs[i] = 0;
+      mapsObject->hg0VALs[i] = 0;
+    }
+
+    // get gain maps
+    fstream infile;
+    sprintf(savename,"../JFModuleCalib_CuFBPCS/data/M%3.3d/gainMaps_M%3.3d.bin", modules[j], modules[j]);
+    infile.open((char*)savename, ios::in | ios::binary);
+    
+    infile.read((char*)mapsObject->g0VALs, sizeof(mapsObject->g0VALs));
+    infile.read((char*)mapsObject->g1VALs, sizeof(mapsObject->g1VALs));
+    infile.read((char*)mapsObject->g2VALs, sizeof(mapsObject->g2VALs));
+    infile.read((char*)mapsObject->hg0VALs, sizeof(mapsObject->hg0VALs));
+
+    infile.close();
+
+    for (int i = 0; i < NCH; i++) {
+      g0hist->Fill(mapsObject->g0VALs[i]);
+      g1hist->Fill(mapsObject->g1VALs[i]);
+      g2hist->Fill(mapsObject->g2VALs[i]);
+      hg0hist->Fill(mapsObject->hg0VALs[i]);
+    }
+
+    modules_d.push_back(modules[j]);
+    module_errs.push_back(0.);
+    g0_means.push_back(g0hist->GetMean());
+    g0_rmss.push_back(g0hist->GetRMS());
+    g1_means.push_back(g1hist->GetMean());
+    g1_rmss.push_back(g1hist->GetRMS());
+    g2_means.push_back(g2hist->GetMean());
+    g2_rmss.push_back(g2hist->GetRMS());
+    hg0_means.push_back(hg0hist->GetMean());
+    hg0_rmss.push_back(hg0hist->GetRMS());
+
+    g0hist->Reset();
+    g1hist->Reset();
+    g2hist->Reset();
+    hg0hist->Reset();
+    
+  }
+
+  TCanvas *c1 = new TCanvas("c1","",150,10,800,400);
+  c1->SetLeftMargin(0.12);
+  c1->SetRightMargin(0.03);
+  c1->SetTopMargin(0.08);
+  c1->SetBottomMargin(0.15);
+
+  TGraphErrors* grap_g0 = new TGraphErrors(modules_d.size(), &(modules_d[0]), &(g0_means[0]), &(module_errs[0]), &(g0_rmss[0]));
+  grap_g0->SetMarkerStyle(20);
+  grap_g0->GetXaxis()->SetTitle("Module");
+  grap_g0->GetYaxis()->SetTitle("Mean G0 [ADU/keV]");
+  grap_g0->GetYaxis()->SetTitleOffset(0.9);
+  grap_g0->Draw("AP");
+  c1->SaveAs("plots/ModuleSummary_g0.png");
+
+  TGraphErrors* grap_g1 = new TGraphErrors(modules_d.size(), &(modules_d[0]), &(g1_means[0]), &(module_errs[0]), &(g1_rmss[0]));
+  grap_g1->SetMarkerStyle(20);
+  grap_g1->GetXaxis()->SetTitle("Module");
+  grap_g1->GetYaxis()->SetTitle("Mean G1 [ADU/keV]");
+  grap_g1->GetYaxis()->SetTitleOffset(0.9);
+  grap_g1->Draw("AP");
+  c1->SaveAs("plots/ModuleSummary_g1.png");
+
+  TGraphErrors* grap_g2 = new TGraphErrors(modules_d.size(), &(modules_d[0]), &(g2_means[0]), &(module_errs[0]), &(g2_rmss[0]));
+  grap_g2->SetMarkerStyle(20);
+  grap_g2->GetXaxis()->SetTitle("Module");
+  grap_g2->GetYaxis()->SetTitle("Mean G2 [ADU/keV]");
+  grap_g2->GetYaxis()->SetTitleOffset(0.9);
+  grap_g2->Draw("AP");
+  c1->SaveAs("plots/ModuleSummary_g2.png");
+
+  TGraphErrors* grap_hg0 = new TGraphErrors(modules_d.size(), &(modules_d[0]), &(hg0_means[0]), &(module_errs[0]), &(hg0_rmss[0]));
+  grap_hg0->SetMarkerStyle(20);
+  grap_hg0->GetXaxis()->SetTitle("Module");
+  grap_hg0->GetYaxis()->SetTitle("Mean HG0 [ADU/keV]");
+  grap_hg0->GetYaxis()->SetTitleOffset(0.9);
+  grap_hg0->Draw("AP");
+  c1->SaveAs("plots/ModuleSummary_hg0.png");
+  
+}
diff --git a/makefile b/makefile
index d4f671c..afb570e 100644
--- a/makefile
+++ b/makefile
@@ -14,5 +14,11 @@ JFMC_CurrentSourceScan: JFMC_CurrentSourceScan.cpp
 JFMC_CalibWriter: JFMC_CalibWriter.cpp
 	g++ -Wall -O3 -m64 -I$(ROOTSYS)/include -L$(ROOTSYS)/lib -lGui -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d -lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic JFMC_CalibWriter.cpp -o JFMC_CalibWriter
 
-dummyRootFileMaker: dummyRootFileMaker.cpp
-	g++ -Wall -O3 -m64 -I$(ROOTSYS)/include -L$(ROOTSYS)/lib -lGui -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d -lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic dummyRootFileMaker.cpp -o dummyRootFileMaker
\ No newline at end of file
+CS_BaselineTest: CS_BaselineTest.cpp
+	g++ -Wall -O3 -m64 -I$(ROOTSYS)/include -L$(ROOTSYS)/lib -lGui -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d -lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic CS_BaselineTest.cpp -o CS_BaselineTest
+
+BP_DB_Comp: BP_DB_Comp.cpp
+	g++ -Wall -O3 -m64 -I$(ROOTSYS)/include -L$(ROOTSYS)/lib -lGui -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d -lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic BP_DB_Comp.cpp -o BP_DB_Comp
+
+JFMC_ModuleComp: JFMC_ModuleComp.cpp
+	g++ -Wall -O3 -m64 -I$(ROOTSYS)/include -L$(ROOTSYS)/lib -lGui -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d -lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic JFMC_ModuleComp.cpp -o JFMC_ModuleComp
\ No newline at end of file