Incomplete codes for all 256 pixels

gain_measurement.cpp

@@ -45,7 +45,7 @@ TH1D *fit_1d(TH1D* &hist_255,int colors,int markers,int &amplitude,double fitmin
 
 
 
-void gain_6()
+std::vector<TH1D*> gain_measurement(int i, int j)
 {
   int x,y, nbins = 450;
   double xmin=850.0,xmax = 1750.0;
@@ -58,10 +58,14 @@ void gain_6()
 
   ifstream amp03,amp02,amp01;
 
-
-  amp03.open("amp03_255_255.dat");
-  amp02.open("amp02_255_255.dat");
-  amp01.open("amp01_255_255.dat");
+  stringstream sfilename;
+  sfilename << std::setw(3) << std::setfill('0') << i << "_"
+	    << std::setw(3) << std::setfill('0') << j <<".dat";
+  std::cout << sfilename.str() << std::endl;
+  amp03.open(std::string("amp03/" + sfilename.str()).c_str());
+  amp02.open(sfilename.str().c_str());
+  amp01.open(sfilename.str().c_str());
+    
 
   if(amp03.is_open())
     {
@@ -102,7 +106,7 @@ void gain_6()
     }
   if(amp01.is_open())
     {
-      cout<<"The data file 255_255 is opened..."<<std::endl;
+      cout<<"The data file " << i << "_" << j << " is opened..."<<std::endl;
       TH1D *h1 = new TH1D("h1","Threshold gain measurement",nbins,xmin,xmax);
       hist_255.push_back(h1);
       while(!amp01.eof())
@@ -134,5 +138,17 @@ void gain_6()
 
 
   cout<<"Done..."<<std::endl;
+  return hist_255;
+}
 
+void complete_plot()
+{
+  for(int row=0;row<256;row++)
+    {
+      for(int col=0;col<256;col++)
+	if(row==col)
+	  {
+	    auto histograms = gain_measurement(row,col);
+	  }
+    }
 }