// to analyse the backplane pulsing data per module
// changes by VH 210906: to eliminate hardcoded absolute paths, uses location of the analysis root files as additional input argument (accordingly changed in filename_creator.sh)

#include "TApplication.h"
#include "sls_detector_calibration/jungfrauCommonHeader.h"
#include "sls_detector_calibration/jungfrauCommonFunctions.h"

#include "sls_detector_calibration/jungfrauFile.C"
#include "sls_detector_calibration/jungfrauPedestal.C"

#include "TGraph.h"
#include "TGraphErrors.h"
#include "TF1.h"
#include "TFile.h"
#include "TPaveStats.h"
#include "TLegend.h"
#include "TPaveText.h"
#include "TCanvas.h"
#include "TFitResultPtr.h"
#include "TFitResult.h"

#include <algorithm>
#include <sys/stat.h>
#include <sstream>


TGraphErrors *grap_g0;
TGraphErrors *grap_g1;
TGraphErrors *grap_g2;

TF1 *fit_g0=0;
TF1 *fit_g1=0;
TF1 *fit_g2=0;

TGraphErrors *norm_g0;
TGraphErrors *norm_g1;
TGraphErrors *norm_g2;

TF1 *flat_g0;
TF1 *flat_g2;

TF1 *lin_g1_p02pc;
TF1 *lin_g1_p05pc;
TF1 *lin_g1_m02pc;
TF1 *lin_g1_m05pc;

TF1 *lin_g2_p02pc;
TF1 *lin_g2_p05pc;
TF1 *lin_g2_m02pc;
TF1 *lin_g2_m05pc;

Double_t par1;
Double_t par2;

double checkRangeMinForRCEffect(double range_min) {

  // check that the range minimum is at least 2 us, to avoid non-linearity of RC effect
     if (range_min < 2.0) {
         return 2.0;
     } else {
               return range_min;
     }
  
}
                 
 

bool isHGX=false;
int main(int argc, char* argv[]) {

 //nonblock(NB_ENABLE);
 cout <<"opening the rootapp" <<endl;
 TApplication rootapp("example",&argc, argv);
  
  jungfrauStyle();
  //gROOT->SetBatch(1);
  gStyle->SetOptFit(11);

  /*
  if (argc != 3) {
    cout << "Correct usage:" << endl;
    cout << "arg 1: specify module number" << endl;
    cout << "arg 2: specify data location" << endl;
    cout << "arg 3: specify column (x)" << endl;
    cout << "arg 4: specify row (y)" << endl;
    exit(1);
  }
  */ //uncomment for SR

  if (argc != 4) {  
    cout << "Correct usage:" << endl;
    cout << "arg 1: specify module number" << endl;
    cout << "arg 2: specify pixel x position from 0 to 1023" << endl;
    cout << "arg 3: specify pixel y position from 0 to 511" << endl;
    exit(1);
  } //uncomment for VH 210906

  string module_str = argv[1];

  string str2 =("HG0G1G2");
  string str3 =("HGOG1G2"); //filename creator had this bug 
   
  int column = atoi(argv[2]); 
  int row = atoi(argv[3]);
  
  int pixel = column+row*1024;

  char data_loc[256];
  sprintf(data_loc,"/mnt/sls_det_storage/jungfrau_calib/jungfrau_ana_sophie/M%s_CalibAna/", module_str.c_str()); 
  cout << data_loc << endl;
  
  std::string folder_path(data_loc);
  if (folder_path.find(str2) != string::npos) isHGX=true;  
  if (folder_path.find(str3) != string::npos) isHGX=true; 


  // cout << data_loc.find(str2)<<" "  << string::npos << " " << str2 << " " << data_loc <<endl;


  if (isHGX) {
    cout << " HG0->HG1->HG2 sequence - dynamicHG0" <<endl;
   // plotfolder_str="BackplanePulsing_HG0G1G2"; 
  }
  else {
    cout << " G0->G1->G2 sequence - dynamicG0" <<endl;
   // plotfolder_str="BackplanePulsing"; 
  }

 
 // Define vector with module name for the different exposure time scan

  vector<string> scan_v1;
  scan_v1.push_back("002");
  scan_v1.push_back("003");

  vector<string> scan_v2;
  scan_v2.push_back("032");

  vector<string> scan_v3;

  vector<string> scan_v4;
  scan_v4.push_back("040");
  scan_v4.push_back("044");

  vector<string> scan_v5;
  scan_v5.push_back("006");
  scan_v5.push_back("022");
  scan_v5.push_back("055");
  scan_v5.push_back("065");
  scan_v5.push_back("066");
  scan_v5.push_back("068");
  scan_v5.push_back("069");
  scan_v5.push_back("088");
  scan_v5.push_back("094");

  vector<string> scan_v6;
  scan_v6.push_back("008");
  scan_v6.push_back("011");
  scan_v6.push_back("021");
  scan_v6.push_back("023");
  scan_v6.push_back("025");
  scan_v6.push_back("029");
  scan_v6.push_back("035");
  scan_v6.push_back("039");
  scan_v6.push_back("043");
  scan_v6.push_back("045");
  scan_v6.push_back("046");
  scan_v6.push_back("048");
  scan_v6.push_back("049");
  scan_v6.push_back("050");
  scan_v6.push_back("051");
  scan_v6.push_back("053");
  scan_v6.push_back("054");
  scan_v6.push_back("057");
  scan_v6.push_back("058");
  scan_v6.push_back("059");
  scan_v6.push_back("060");
  scan_v6.push_back("061");
  scan_v6.push_back("062");
  scan_v6.push_back("063");
  scan_v6.push_back("064");
  scan_v6.push_back("067");
  scan_v6.push_back("070");
  scan_v6.push_back("071");
  scan_v6.push_back("072");
  scan_v6.push_back("073");
  scan_v6.push_back("074");
  scan_v6.push_back("075");
  scan_v6.push_back("076");
  scan_v6.push_back("079");
  scan_v6.push_back("080");
  scan_v6.push_back("086");
  scan_v6.push_back("089");
  scan_v6.push_back("090");
  scan_v6.push_back("091");
  scan_v6.push_back("092");
  scan_v6.push_back("093");
  scan_v6.push_back("095");
  scan_v6.push_back("097");
  scan_v6.push_back("098");
  scan_v6.push_back("099");
  scan_v6.push_back("100");
  scan_v6.push_back("102");
  scan_v6.push_back("103");
  scan_v6.push_back("104");
  scan_v6.push_back("106");
  scan_v6.push_back("107");
  scan_v6.push_back("108");
  scan_v6.push_back("109");
  scan_v6.push_back("110");
  scan_v6.push_back("111");
  scan_v6.push_back("112");
  scan_v6.push_back("113");
  scan_v6.push_back("114");
  scan_v6.push_back("115");
  scan_v6.push_back("116");
  scan_v6.push_back("117");
  scan_v6.push_back("118");
  scan_v6.push_back("119");
  scan_v6.push_back("120");
  scan_v6.push_back("121");
  scan_v6.push_back("122");
  scan_v6.push_back("123");
  scan_v6.push_back("124");
  scan_v6.push_back("125");
  scan_v6.push_back("126");
  scan_v6.push_back("129");
  scan_v6.push_back("130");
  scan_v6.push_back("131");
  scan_v6.push_back("132");
  scan_v6.push_back("133");
  scan_v6.push_back("134");
  scan_v6.push_back("135");
  scan_v6.push_back("200");
  scan_v6.push_back("201");
  scan_v6.push_back("202");
  scan_v6.push_back("203");
  scan_v6.push_back("204");
  scan_v6.push_back("205");
  scan_v6.push_back("206");
  scan_v6.push_back("208");
  scan_v6.push_back("209");
  scan_v6.push_back("210");
  scan_v6.push_back("212");
  scan_v6.push_back("214");
  scan_v6.push_back("215");
  scan_v6.push_back("216");
  scan_v6.push_back("217");
  scan_v6.push_back("218");
  scan_v6.push_back("219");
  scan_v6.push_back("221");
  scan_v6.push_back("222");
  scan_v6.push_back("223");
  scan_v6.push_back("224");
  scan_v6.push_back("225");
  scan_v6.push_back("226");
  scan_v6.push_back("227");
  scan_v6.push_back("228");
  scan_v6.push_back("229");
  scan_v6.push_back("232");
  scan_v6.push_back("233");
  scan_v6.push_back("234");
  scan_v6.push_back("235");
  scan_v6.push_back("236");
  scan_v6.push_back("242");
  scan_v6.push_back("244");
  scan_v6.push_back("245");
  scan_v6.push_back("246");
  scan_v6.push_back("248");
  scan_v6.push_back("249");
  scan_v6.push_back("250");
  scan_v6.push_back("251");
  scan_v6.push_back("252");
  scan_v6.push_back("253");
  scan_v6.push_back("254");
  scan_v6.push_back("255");
  scan_v6.push_back("260");
  scan_v6.push_back("261");
  scan_v6.push_back("263");
  scan_v6.push_back("264");
  scan_v6.push_back("265");
  scan_v6.push_back("266");
  scan_v6.push_back("267");
  scan_v6.push_back("269");
  scan_v6.push_back("270");
  scan_v6.push_back("272");
  scan_v6.push_back("273");
  scan_v6.push_back("274");
  scan_v6.push_back("275");
  scan_v6.push_back("276");
  scan_v6.push_back("277");
  scan_v6.push_back("278");
  scan_v6.push_back("279");
  scan_v6.push_back("281");
  scan_v6.push_back("282");
  scan_v6.push_back("283");
  scan_v6.push_back("285");
  scan_v6.push_back("286");
  scan_v6.push_back("287");
  scan_v6.push_back("288");
  scan_v6.push_back("289");
  scan_v6.push_back("290");
  scan_v6.push_back("291");
  scan_v6.push_back("292");
  scan_v6.push_back("294");
  scan_v6.push_back("296");
  scan_v6.push_back("298");
  scan_v6.push_back("299");
  scan_v6.push_back("301");
  scan_v6.push_back("307");
  scan_v6.push_back("312");
  scan_v6.push_back("351");
  scan_v6.push_back("352");
  scan_v6.push_back("354");
  scan_v6.push_back("360");
  scan_v6.push_back("373");
  scan_v6.push_back("380");
  scan_v6.push_back("304");
  scan_v6.push_back("321");
  scan_v6.push_back("346");
  scan_v6.push_back("375");
  scan_v6.push_back("381");
  scan_v6.push_back("383");
  scan_v6.push_back("310");
  scan_v6.push_back("313");
  scan_v6.push_back("314");
  scan_v6.push_back("334");
  scan_v6.push_back("337");
  scan_v6.push_back("347");
  scan_v6.push_back("361");
  scan_v6.push_back("362");
  scan_v6.push_back("355");
  scan_v6.push_back("364");
  scan_v6.push_back("366");
  scan_v6.push_back("243");
  scan_v6.push_back("339");
  scan_v6.push_back("341");
  scan_v6.push_back("368");
  scan_v6.push_back("376");
  scan_v6.push_back("379");

  scan_v6.push_back("359");
  scan_v6.push_back("363");
  scan_v6.push_back("365");
  scan_v6.push_back("347");
  scan_v6.push_back("366");

  scan_v6.push_back("268");
  scan_v6.push_back("300");
  scan_v6.push_back("302");
  scan_v6.push_back("303");
  scan_v6.push_back("309");
  scan_v6.push_back("326");

  scan_v6.push_back("308");
  scan_v6.push_back("322");
  scan_v6.push_back("328");
  scan_v6.push_back("332");
  scan_v6.push_back("335");
  scan_v6.push_back("353");
  scan_v6.push_back("335_FW2_2");
  scan_v6.push_back("353_FW2_2");

  scan_v6.push_back("262_FWcheck_1_1_1");
  scan_v6.push_back("262_FWcheck_1_1_1_test");
  scan_v6.push_back("262_FWcheck_1_2_rc210831");
  scan_v6.push_back("262_FWcheck_1_2_defdacs");

  scan_v6.push_back("311");
  scan_v6.push_back("316");
  scan_v6.push_back("318");
  scan_v6.push_back("319");
  scan_v6.push_back("320");
  scan_v6.push_back("324");

  scan_v6.push_back("272_fullspeedtest");
  scan_v6.push_back("272_halfspeedtest");

  scan_v6.push_back("352_fullspeed"); //jf4mb2-00 (SLS)
  scan_v6.push_back("307_fullspeed"); //jf4mb2-01 (SLS)
  scan_v6.push_back("264_fullspeed"); //jf4mb2-02 (SLS)
  scan_v6.push_back("253_fullspeed"); //jf4mb2-03 (SLS)
  scan_v6.push_back("351_fullspeed"); //jf4mb2-04 (SLS)
  scan_v6.push_back("312_fullspeed"); //jf4mb2-05 (SLS)
  scan_v6.push_back("261_fullspeed"); //jf4mb2-06 (SLS)
  scan_v6.push_back("310_fullspeed"); //jf4mb2-07 (SLS)

  scan_v6.push_back("384"); //jfcalib-5
  scan_v6.push_back("510"); //jfcalib-0
  scan_v6.push_back("511"); //jfcalib-1
  scan_v6.push_back("512"); //jfcalib-2
  scan_v6.push_back("556"); //jfcalib-3
  scan_v6.push_back("572"); //jfcalib-4

  scan_v6.push_back("329"); //jfcalib-0
  scan_v6.push_back("331"); //jfcalib-1
  scan_v6.push_back("333"); //jfcalib-2
  scan_v6.push_back("336"); //jfcalib-3
  scan_v6.push_back("356"); //jfcalib-4
  scan_v6.push_back("377"); //jfcalib-5

  scan_v6.push_back("387"); //jfcalib-0
  scan_v6.push_back("388"); //jfcalib-1
  scan_v6.push_back("389"); //jfcalib-2
  scan_v6.push_back("396"); //jfcalib-3
  scan_v6.push_back("397"); //jfcalib-4
  scan_v6.push_back("398"); //jfcalib-5

  scan_v6.push_back("390"); //jfcalib-0
  scan_v6.push_back("391"); //jfcalib-1
  scan_v6.push_back("393"); //jfcalib-2
  scan_v6.push_back("399"); //jfcalib-3
  scan_v6.push_back("516"); //jfcalib-4
  scan_v6.push_back("517"); //jfcalib-5

  scan_v6.push_back("519"); //jfcalib-0
  scan_v6.push_back("520"); //jfcalib-1
  scan_v6.push_back("532"); //jfcalib-2
  scan_v6.push_back("539"); //jfcalib-3
  scan_v6.push_back("555"); //jfcalib-4

  //Calib for VELA beamtime in March 2022
  scan_v6.push_back("262"); //jfcalib-0
  scan_v6.push_back("231"); //jfcalib-1
  scan_v6.push_back("305"); //jfcalib-2
  
  scan_v6.push_back("392"); //jfcalib-1
  scan_v6.push_back("394"); //jfcalib-2
  scan_v6.push_back("513"); //jfcalib-3
  scan_v6.push_back("514"); //jfcalib-4
  scan_v6.push_back("515"); //jfcalib-5
  
  scan_v6.push_back("239"); //jfcalib-0  //bad
  scan_v6.push_back("064"); //jfcalib-1  //bad //already added above
  scan_v6.push_back("023"); //jfcalib-2  //bad //already added above
  scan_v6.push_back("349"); //jfcalib-3  //bad
  scan_v6.push_back("518"); //jfcalib-4
  scan_v6.push_back("521"); //jfcalib-5
  
  scan_v6.push_back("523"); //jfcalib-0
  scan_v6.push_back("526"); //jfcalib-1
  scan_v6.push_back("529"); //jfcalib-2
  scan_v6.push_back("530"); //jfcalib-3
  scan_v6.push_back("531"); //jfcalib-4
  scan_v6.push_back("533"); //jfcalib-5
  
  scan_v6.push_back("524"); //1M for ESRF
  scan_v6.push_back("525"); //1M for ESRF
  
  scan_v6.push_back("306"); //jfcalib-0 (strixel)
  scan_v6.push_back("350"); //jfcalib-1 (strixel)
  
  scan_v6.push_back("534"); //jfcalib-0
  scan_v6.push_back("537"); //jfcalib-1
  scan_v6.push_back("538"); //jfcalib-2
  scan_v6.push_back("542"); //jfcalib-3
  scan_v6.push_back("543"); //jfcalib-4
  scan_v6.push_back("544"); //jfcalib-5

  scan_v6.push_back("402"); //jfcalib-0
  scan_v6.push_back("404"); //jfcalib-1
  scan_v6.push_back("407"); //jfcalib-2
  scan_v6.push_back("409"); //jfcalib-3
  scan_v6.push_back("410"); //jfcalib-4
  scan_v6.push_back("411"); //jfcalib-5

  //scan_v6.push_back("008"); //jfcalib-0 //already in the list from previous calibration
  scan_v6.push_back("433"); //jfcalib-1
  scan_v6.push_back("550"); //jfcalib-2
  scan_v6.push_back("545"); //jfcalib-3
  scan_v6.push_back("546"); //jfcalib-4
  scan_v6.push_back("547"); //jfcalib-5

  //scan_v6.push_back("416"); //jfcalib-0 //destroyed, machanical sample
  scan_v6.push_back("417"); //jfcalib-1
  scan_v6.push_back("418"); //jfcalib-2
  scan_v6.push_back("420"); //jfcalib-3
  scan_v6.push_back("421"); //jfcalib-4
  scan_v6.push_back("430"); //jfcalib-5

  scan_v6.push_back("417_fullspeed"); //jfcalib-1
  scan_v6.push_back("418_fullspeed"); //jfcalib-2
  scan_v6.push_back("420_fullspeed"); //jfcalib-3
  scan_v6.push_back("421_fullspeed"); //jfcalib-4
  scan_v6.push_back("430_fullspeed"); //jfcalib-5

  scan_v6.push_back("434"); //jfcalib-0
  scan_v6.push_back("437"); //jfcalib-1
  scan_v6.push_back("480"); //jfcalib-2
  scan_v6.push_back("464"); //jfcalib-3
  scan_v6.push_back("473"); //jfcalib-4
  scan_v6.push_back("474"); //jfcalib-5

  scan_v6.push_back("076"); //jfcalib-0
  scan_v6.push_back("089"); //jfcalib-3
  scan_v6.push_back("207"); //jfcalib-5

  scan_v6.push_back("244"); //jfcalib-3

  scan_v6.push_back("474_fullspeed"); //jfcalib-0
  scan_v6.push_back("437_fullspeed"); //jfcalib-1
  scan_v6.push_back("480_fullspeed"); //jfcalib-2
  scan_v6.push_back("482_fullspeed"); //jfcalib-3
  scan_v6.push_back("487_fullspeed"); //jfcalib-5

  scan_v6.push_back("482"); //jfcalib-3
  scan_v6.push_back("487"); //jfcalib-5

  scan_v6.push_back("499_fullspeed"); //jfcalib-0
  scan_v6.push_back("548_fullspeed"); //jfcalib-1
  scan_v6.push_back("549_fullspeed"); //jfcalib-2
  scan_v6.push_back("427_fullspeed"); //jfcalib-3

  scan_v6.push_back("499"); //jfcalib-0
  scan_v6.push_back("548"); //jfcalib-1
  scan_v6.push_back("549"); //jfcalib-2
  scan_v6.push_back("427"); //jfcalib-3
  scan_v6.push_back("435"); //jfcalib-5

  scan_v6.push_back("403"); //jfcalib-0
  scan_v6.push_back("406"); //jfcalib-1
  scan_v6.push_back("423"); //jfcalib-2
  scan_v6.push_back("436"); //jfcalib-3
  scan_v6.push_back("479"); //jfcalib-5

  // Calibration 450 um sensor for ESRF
 
  scan_v6.push_back("583"); //jfcalib-0
  scan_v6.push_back("593"); //jfcalib-1
  scan_v6.push_back("601"); //jfcalib-2
  scan_v6.push_back("604"); //jfcalib-3
  scan_v6.push_back("605"); //jfcalib-5

  // Calibration 320 um sensors

  scan_v6.push_back("441"); //jfcalib-0
  scan_v6.push_back("460"); //jfcalib-1
  scan_v6.push_back("463"); //jfcalib-2
  scan_v6.push_back("476"); //jfcalib-3
  scan_v6.push_back("371"); //jfcalib-4

  scan_v6.push_back("429"); //jfcalib-0
  scan_v6.push_back("442"); //jfcalib-1
  scan_v6.push_back("444"); //jfcalib-2
  scan_v6.push_back("446"); //jfcalib-3
  scan_v6.push_back("447"); //jfcalib-4

  scan_v6.push_back("428"); //jfcalib-0
  scan_v6.push_back("439"); //jfcalib-1
  scan_v6.push_back("443"); //jfcalib-2
  scan_v6.push_back("445"); //jfcalib-3

  scan_v6.push_back("448"); //jfcalib-0
  scan_v6.push_back("454"); //jfcalib-1
  scan_v6.push_back("455"); //jfcalib-2
  scan_v6.push_back("456"); //jfcalib-3
  scan_v6.push_back("466"); //jfcalib-4
  scan_v6.push_back("471"); //jfcalib-5

  scan_v6.push_back("481"); //jfcalib-0
  scan_v6.push_back("488"); //jfcalib-1
  scan_v6.push_back("490"); //jfcalib-2
  scan_v6.push_back("491"); //jfcalib-3
  scan_v6.push_back("492"); //jfcalib-4
  scan_v6.push_back("493"); //jfcalib-5

  scan_v6.push_back("494"); //jfcalib-0
  scan_v6.push_back("495"); //jfcalib-1
  scan_v6.push_back("497"); //jfcalib-2
  scan_v6.push_back("498"); //jfcalib-3
  scan_v6.push_back("501"); //jfcalib-4
  scan_v6.push_back("502"); //jfcalib-5

  scan_v6.push_back("036"); //jfcalib-0
  scan_v6.push_back("036_fullspeed"); //jfcalib-0
  scan_v6.push_back("560"); //jfcalib-1
  scan_v6.push_back("536"); //jfcalib-2
  scan_v6.push_back("535"); //jfcalib-3
  scan_v6.push_back("527"); //jfcalib-4
  scan_v6.push_back("522"); //jfcalib-5

  scan_v6.push_back("571"); //jfcalib-0
  scan_v6.push_back("571_fullspeed"); //jfcalib-0
  scan_v6.push_back("573"); //jfcalib-1
  scan_v6.push_back("573_fullspeed"); //jfcalib-1
  scan_v6.push_back("574"); //jfcalib-2
  scan_v6.push_back("574_fullspeed"); //jfcalib-2
  scan_v6.push_back("575"); //jfcalib-3
  scan_v6.push_back("575_fullspeed"); //jfcalib-3
  scan_v6.push_back("576"); //jfcalib-4
  scan_v6.push_back("576_fullspeed"); //jfcalib-4
  scan_v6.push_back("577"); //jfcalib-5
  scan_v6.push_back("577_fullspeed"); //jfcalib-5

  scan_v6.push_back("563"); //jfcalib-0
  scan_v6.push_back("563_fullspeed"); //jfcalib-0
  scan_v6.push_back("564"); //jfcalib-1
  scan_v6.push_back("564_fullspeed"); //jfcalib-1
  scan_v6.push_back("566"); //jfcalib-2
  scan_v6.push_back("566_fullspeed"); //jfcalib-2
  scan_v6.push_back("567"); //jfcalib-3
  scan_v6.push_back("567_fullspeed"); //jfcalib-3
  scan_v6.push_back("568"); //jfcalib-4
  scan_v6.push_back("568_fullspeed"); //jfcalib-4
  scan_v6.push_back("569"); //jfcalib-5
  scan_v6.push_back("569_fullspeed"); //jfcalib-5
  
  scan_v6.push_back("579"); //QUAD thin
  scan_v6.push_back("579_fullspeed"); //QUAD thin

  scan_v6.push_back("454_fullspeed"); //jfcalib-0
  scan_v6.push_back("466_fullspeed"); //jfcalib-1
  scan_v6.push_back("536_fullspeed"); //jfcalib-2
  scan_v6.push_back("561"); //jfcalib-4
  scan_v6.push_back("561_fullspeed"); //jfcalib-3
  scan_v6.push_back("562"); //jfcalib-4
  scan_v6.push_back("562_fullspeed"); //jfcalib-4
  scan_v6.push_back("038"); //jfcalib-5
  scan_v6.push_back("038_fullspeed"); //jfcalib-5

  scan_v6.push_back("431_JF12");


  scan_v6.push_back("428_fullspeed"); //jfcalib-0
  scan_v6.push_back("443_fullspeed"); //jfcalib-1
  scan_v6.push_back("444_fullspeed"); //jfcalib-2
  scan_v6.push_back("455_fullspeed"); //jfcalib-3
  scan_v6.push_back("471_fullspeed"); //jfcalib-4
  scan_v6.push_back("493_fullspeed"); //jfcalib-5

  scan_v6.push_back("448_fullspeed"); //jfcalib-0
  scan_v6.push_back("497_fullspeed"); //jfcalib-1
  scan_v6.push_back("498_fullspeed"); //jfcalib-2
  scan_v6.push_back("501_fullspeed"); //jfcalib-3
  scan_v6.push_back("502_fullspeed"); //jfcalib-4
  scan_v6.push_back("527_fullspeed"); //jfcalib-5

  scan_v6.push_back("580_fullspeed"); //jfcalib-5
  scan_v6.push_back("580"); //jfcalib-5

  scan_v6.push_back("596"); //jfcalib-0
  scan_v6.push_back("596_fullspeed"); //jfcalib-0
  scan_v6.push_back("597"); //jfcalib-1
  scan_v6.push_back("597_fullspeed"); //jfcalib-1
  scan_v6.push_back("598"); //jfcalib-2
  scan_v6.push_back("598_fullspeed"); //jfcalib-2
  scan_v6.push_back("599"); //jfcalib-3
  scan_v6.push_back("599_fullspeed"); //jfcalib-3
  scan_v6.push_back("602"); //jfcalib-4
  scan_v6.push_back("602_fullspeed"); //jfcalib-4
  scan_v6.push_back("603"); //jfcalib-5
  scan_v6.push_back("603_fullspeed"); //jfcalib-5

<<<<<<< HEAD
=======

>>>>>>> refs/remotes/origin/master
  scan_v6.push_back("578"); //jfcalib-0
  scan_v6.push_back("578_fullspeed"); //jfcalib-0
  scan_v6.push_back("588"); //jfcalib-2
  scan_v6.push_back("588_fullspeed"); //jfcalib-2
  scan_v6.push_back("592"); //jfcalib-3
  scan_v6.push_back("592_fullspeed"); //jfcalib-3
  scan_v6.push_back("590"); //jfcalib-4
  scan_v6.push_back("590_fullspeed"); //jfcalib-4
  scan_v6.push_back("591"); //jfcalib-5
  scan_v6.push_back("591_fullspeed"); //jfcalib-5

  scan_v6.push_back("594"); //jfcalib-0
  scan_v6.push_back("594_fullspeed"); //jfcalib-0
  scan_v6.push_back("595"); //jfcalib-1
  scan_v6.push_back("595_fullspeed"); //jfcalib-1
  scan_v6.push_back("606"); //jfcalib-3
  scan_v6.push_back("606_fullspeed"); //jfcalib-3
  scan_v6.push_back("608"); //jfcalib-4
  scan_v6.push_back("608_fullspeed"); //jfcalib-4
  scan_v6.push_back("589"); //jfcalib-5
  scan_v6.push_back("589_fullspeed"); //jfcalib-5

<<<<<<< HEAD
=======

>>>>>>> refs/remotes/origin/master
  scan_v6.push_back("586"); //jfcalib-0
  scan_v6.push_back("586_fullspeed"); //jfcalib-0
  scan_v6.push_back("607"); //jfcalib-1
  scan_v6.push_back("607_fullspeed"); //jfcalib-1
  scan_v6.push_back("609"); //jfcalib-2
  scan_v6.push_back("609_fullspeed"); //jfcalib-2
  scan_v6.push_back("625"); //jfcalib-3
  scan_v6.push_back("625_fullspeed"); //jfcalib-3
<<<<<<< HEAD
  scan_v6.push_back("626"); //jfcalib-4
  scan_v6.push_back("626_fullspeed"); //jfcalib-4

  scan_v6.push_back("607"); //jfcalib-1
  scan_v6.push_back("607_fullspeed"); //jfcalib-1
  scan_v6.push_back("609"); //jfcalib-2
  scan_v6.push_back("609_fullspeed"); //jfcalib-2
  scan_v6.push_back("367"); //jfcalib-4
  scan_v6.push_back("367_fullspeed"); //jfcalib-4
  scan_v6.push_back("604"); //jfcalib-5
  scan_v6.push_back("604_fullspeed"); //jfcalib-5


  scan_v6.push_back("367_Kirsty"); //jfcalib-4
  scan_v6.push_back("367_Kirsty_fullspeed"); //jfcalib-4
  scan_v6.push_back("604_Kirsty"); //jfcalib-5
  scan_v6.push_back("604_Kirsty_fullspeed"); //jfcalib-5


  scan_v6.push_back("628"); //jfcalib-0
  scan_v6.push_back("628_fullspeed"); //jfcalib-0
  scan_v6.push_back("627"); //jfcalib-1
  scan_v6.push_back("627_fullspeed"); //jfcalib-1
  scan_v6.push_back("629"); //jfcalib-5
  scan_v6.push_back("629_fullspeed"); //jfcalib-5
  scan_v6.push_back("642"); //jfcalib-3
  scan_v6.push_back("642_fullspeed"); //jfcalib-3
  scan_v6.push_back("604"); //jfcalib-4
  scan_v6.push_back("604_fullspeed"); //jfcalib-4
  scan_v6.push_back("701_Vadym");
=======
  scan_v6.push_back("526"); //jfcalib-4
  scan_v6.push_back("526_fullspeed"); //jfcalib-4

>>>>>>> refs/remotes/origin/master



   // filter out points at zero and points with 0 error (mean taken from one entry in hist)
    int m_max = 0;
    if (find(scan_v2.begin(), scan_v2.end(), module_str) != scan_v2.end()) {
      m_max = 27;
    } else if (find(scan_v4.begin(), scan_v4.end(), module_str) != scan_v4.end()) {
      m_max = 31;
    } else if (find(scan_v3.begin(), scan_v3.end(), module_str) != scan_v3.end()) {
      m_max = 35;
    } else if (find(scan_v5.begin(), scan_v5.end(), module_str) != scan_v5.end()) {
      m_max = 51;
    } else if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
      m_max = 97;
    } else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
      m_max = 90;
    }
 
  std::vector<double> G0_pixel(97, 0);
  std::vector<double> G0_err_pixel(97, 0);
  std::vector<double> G1_pixel(97, 0);
  std::vector<double> G1_err_pixel(97, 0);
  std::vector<double> G2_pixel(97, 0);  
  std::vector<double> G2_err_pixel(97, 0);
  
  
  char filename[256]; //uncomment for VH 210902


double filter[97];
  if (find(scan_v4.begin(), scan_v4.end(), module_str) != scan_v4.end()) {
    for (int i = 0; i < 5; i++) {
      filter[i] = 0.05+(i*0.05);
    }
    for (int i = 0; i < 9; i++) {
      filter[i+5] = 0.5+(i*0.25);
    }
    for (int i = 0; i < 17; i++) {
      filter[i+5+9] = 5.+(i*2.5);
    }
  } else if (find(scan_v5.begin(), scan_v5.end(), module_str) != scan_v5.end()) {
    for (int i = 0; i < 9; i++) {
      filter[i] = 0.05+(i*0.05);
    }
    for (int i = 0; i < 17; i++) {
      filter[i+9] = 0.5+(i*0.25);
    }
    for (int i = 0; i < 17; i++) {
      filter[i+9+17] = 5.+(i*2.5);
    }
    for (int i = 0; i < 8; i++) {
      filter[i+9+17+17] = 50.+(i*10.);
    }
  } else if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
    // loop 1
    for (int i = 0; i < 9; i++) {
      filter[i] = 0.05+(i*0.05);
    }
    // loop 2a
    for (int i = 0; i < 11; i++) {
      filter[i+9] = 0.5+(i*0.05);
    }
    // loop 2b
    for (int i = 0; i < 10; i++) {
      filter[i+9+11] = 1.10+(i*0.10);
    }
    // loop 2c
    for (int i = 0; i < 10; i++) {
      filter[i+9+11+10] = 2.25+(i*0.25);
    }
    // loop 3a
    for (int i = 0; i < 11; i++) {
      filter[i+9+11+10+10] = 5.+(i*0.5);
    }
    // loop 3b
    for (int i = 0; i < 10; i++) {
      filter[i+9+11+10+10+11] = 11.+(i*1.0);
    }
    // loop 3c
    for (int i = 0; i < 10; i++) {
      filter[i+9+11+10+10+11+10] = 22.5+(i*2.5);
    }
    // loop 4a
    for (int i = 0; i < 11; i++) {
      filter[i+9+11+10+10+11+10+10] = 50.+(i*5.);
    }
    // loop 4b
    for (int i = 0; i < 15; i++) {
      filter[i+9+11+10+10+11+10+10+11] = 110.+(i*10.);
    }
  } else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
    for (int i = 0; i < 10; i++) {
      filter[i] = 0.2+(i*0.2);
    }
    for (int i = 0; i < 8; i++) {
      filter[i+10] = 3+(i*1);
    }
    for (int i = 0; i < 25; i++) {
      filter[i+10+8] = 12+(i*2);
    }
    for (int i = 0; i < 47; i++) {
      filter[i+10+8+25] = 80+(i*20);
    }
  } else if (find(scan_v3.begin(), scan_v3.end(), module_str) != scan_v3.end()) {
    for (int i = 0; i < 9; i++) {
      filter[i] = 0.05+(i*0.025);
    }
    for (int i = 0; i < 9; i++) {
      filter[i+9] = 0.5+(i*0.25);
    }
    for (int i = 0; i < 17; i++) {
      filter[i+9+9] = 5.+(i*2.5);
    }
  } else if (find(scan_v2.begin(), scan_v2.end(), module_str) != scan_v2.end()) {
    for (int i = 0; i < 9; i++) {
      filter[i] = 0.05+(i*0.025);
    }
    for (int i = 0; i < 9; i++) {
      filter[i+9] = 0.5+(i*0.25);
    }
    for (int i = 0; i < 9; i++) {
      filter[i+9+9] = 5.+(i*2.5);
    }
  }


TCanvas  *A2 = new TCanvas("A2","Plot scan",150,10,800,400);
TCanvas  *A3 = new TCanvas("A3","Plot G1 fit",150,10,800,400);
TCanvas  *A4 = new TCanvas("A4","Plot G1 residuals",150,10,800,400);
TCanvas  *A5 = new TCanvas("A5","Plot G2 fit",150,10,800,400);
TCanvas  *A6 = new TCanvas("A6","Plot G2 residuals",150,10,800,400);
TCanvas  *A7 = new TCanvas("A7","Plot G0 fit",150,10,800,400);
TCanvas  *A8 = new TCanvas("A8","Plot G0 residuals",150,10,800,400);
  
  

  //file name
  sprintf(filename,"%sCS_histos_M%s.root",folder_path.c_str(), module_str.c_str());
  cout << "Loading file " << filename << endl;
  TFile* f = new TFile((const char *)(filename),"READ");
   

  for (int j = 0; j < m_max; j++) {
    TH2F* hist0=(TH2F*)f->Get(Form("adcmap_avg_g0_%d",j));
    G0_pixel[j]=hist0->GetBinContent((column+1),(row+1));
    TH2F* hist0er=(TH2F*)f->Get(Form("adcermap_avg_g0_%d",j));
    G0_err_pixel[j]=hist0er->GetBinContent((column+1),(row+1));
    TH2F* hist1=(TH2F*)f->Get(Form("adcmap_avg_g1_%d",j));
    G1_pixel[j]=hist1->GetBinContent((column+1),(row+1));
    TH2F* hist1er=(TH2F*)f->Get(Form("adcermap_avg_g1_%d",j));    
    G1_err_pixel[j]=hist1er->GetBinContent((column+1),(row+1));
    TH2F* hist2=(TH2F*)f->Get(Form("adcmap_avg_g2_%d",j));
    G2_pixel[j]=hist2->GetBinContent((column+1),(row+1));
    TH2F* hist2er=(TH2F*)f->Get(Form("adcermap_avg_g2_%d",j));    
    G2_err_pixel[j]=hist2er->GetBinContent((column+1),(row+1));
    //cout << "Data for pixel "<< pixel << "is loaded" << endl;
    }  

  
  A2->SetLeftMargin(0.13);
  A2->SetRightMargin(0.05);

    vector<double> r0_adc;
    vector<double> r0_filter;

    vector<double> r0_adcerr;
    vector<double> r0_ferr;

    vector<double> r1_adc;
    vector<double> r1_filter;
    vector<double> r1_adcerr;
    vector<double> r1_ferr;

    vector<double> r2_adc;
    vector<double> r2_filter;
    vector<double> r2_adcerr;
    vector<double> r2_ferr;

    for (int j = 0; j < m_max; j++) {

      if (G0_pixel[j] != 0) {
	    r0_filter.push_back(filter[j]);
	    r0_ferr.push_back(0.);
	    r0_adc.push_back(G0_pixel[j]);
	    r0_adcerr.push_back(G0_err_pixel[j]);
      }
      if (G1_pixel[j] != 0) {
	    r1_filter.push_back(filter[j]);
	    r1_ferr.push_back(0.);
	    r1_adc.push_back(G1_pixel[j]);
	    r1_adcerr.push_back(G1_err_pixel[j]);
      }
      if (G2_pixel[j] != 0) {
	    r2_filter.push_back(filter[j]);
	    r2_ferr.push_back(0.);
	    r2_adc.push_back(G2_pixel[j]);
	    r2_adcerr.push_back(G2_err_pixel[j]);
      }   
    }
    cout << "Filter array is created" << endl;

    double rangemin_g0 = 0;
    double rangemax_g0 = 0;
    double rangemin_g1 = 0;
    double rangemax_g1 = 0;
    double rangemin_g2 = 0;
    double rangemax_g2 = 0;
     
    A2->cd(); 

    // define graphs
    if (r0_adc.size() > 1) {
      grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0]));
      grap_g0->SetMarkerStyle(20);
      grap_g0->SetMarkerColor(kBlue);
      grap_g0->SetLineColor(kBlue);

    }

    if (r1_adc.size() > 1) {
      grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0]));
      grap_g1->SetMarkerStyle(20);
      grap_g1->SetMarkerColor(kGreen+2);
      grap_g1->SetLineColor(kGreen+2);
    }

    if (r2_adc.size() > 1) {
      grap_g2 = new TGraphErrors(r2_adc.size(),&(r2_filter[0]),&(r2_adc[0]),&(r2_ferr[0]),&(r2_adcerr[0]));
      grap_g2->SetMarkerStyle(20);
      grap_g2->SetMarkerColor(kRed);
      grap_g2->SetLineColor(kRed);
    }

    // plot the datapoints
    if (r1_adc.size() > 1) {

    if (find(scan_v6.begin(), scan_v6.end(), module_str) != scan_v6.end()) {
	  grap_g1->GetXaxis()->SetLimits(0.009,5E2);
	} else if (find(scan_v1.begin(), scan_v1.end(), module_str) != scan_v1.end()) {
	  grap_g1->GetXaxis()->SetLimits(0.09,1.1E3);
	} else {
	  grap_g1->GetXaxis()->SetLimits(0.009,1.1E2);
	}

    grap_g1->GetXaxis()->SetTitle("Integration time [#mus]");
    grap_g1->GetYaxis()->SetTitle("ADC [ADU]");
    grap_g1->GetYaxis()->SetTitleOffset(0.9);
	grap_g1->SetMinimum(0);
	grap_g1->SetMaximum(16000);    
    grap_g1->Draw("AP");
	if (r2_adc.size() > 1) {
	  grap_g2->Draw("P");
	}
	if (r0_adc.size() > 1) {
	  grap_g0->Draw("P");
	}
    A2->SetLogx();
    A2->Modified();
    A2->Update();
    
      }
    
    // define fit ranges and fit
    if (r0_adc.size() > 1) {
      A7->cd();
      
      rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end());
      rangemin_g0 = checkRangeMinForRCEffect(rangemin_g0);
      rangemax_g0 = highestPointBeforeSwitching(r0_filter, r1_filter);

     if (rangemax_g0 > rangemin_g0) {
      
	fit_g0 = new TF1("fit_g0","[0]+[1]*x",rangemin_g0, rangemax_g0);
	fit_g0->SetParameter(0, 3000.);
	fit_g0->SetParameter(1, 30000.);
	fit_g0->SetLineColor(kBlue);
	fit_g0->SetParName(0,"G0 const");
	fit_g0->SetParName(1,"G0 grad");
	grap_g0->Fit(fit_g0,"QR+","");
      
	grap_g0->SetMinimum(0);
	grap_g0->SetMaximum(18000);
	grap_g0->GetXaxis()->SetTitle("Integration time [#mus]");
	grap_g0->GetYaxis()->SetTitle("ADC [ADU]");
	grap_g0->GetYaxis()->SetTitleOffset(0.9);
	grap_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
	grap_g0->Draw("AP");
	fit_g0->Draw("same");
	A7->Update();
	if (TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats")) {
	    st0->SetX1NDC(0.2);
	    st0->SetX2NDC(0.54);
	    st0->SetY1NDC(0.71);
	    st0->SetY2NDC(0.90);
	    st0->SetBorderSize(0);
	    st0->SetTextColor(kBlue);
	    A7->SetLogx();
        A7->Modified();
        A7->Update();
      
    //A7->cd();
    }
	vector<double> r0_adc_norm;
	for (size_t j = 0; j < r0_adc.size(); j++) {
	    r0_adc_norm.push_back(r0_adc[j] - fit_g0->Eval(r0_filter[j]));
	}

    A8->cd();

	  TGraphErrors *norm_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_ferr[0]),&(r0_adcerr[0]));
	  norm_g0->SetMarkerColor(kBlue);
	  norm_g0->SetLineColor(kBlue);

	  TF1* flat_g0 = new TF1("flat_g0","0",rangemin_g0, rangemax_g0);
	  flat_g0->SetLineColor(kBlue);

	  norm_g0->GetXaxis()->SetTitle("Integration time [#mus]");
	  norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]");
	  norm_g0->GetYaxis()->SetTitleOffset(0.9);
	  norm_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
	  norm_g0->SetMinimum(-100);
	  norm_g0->SetMaximum(100);
	  norm_g0->Draw("AP");
	  flat_g0->Draw("same");

      A8->Update();

	    //sprintf(savename,"plots/M%s/%s/pixel_%s_%d_g0norm_M%s.png", module_str.c_str(), plotfolder_str.c_str(),pixel_type.c_str(), i,  module_str.c_str());
	    //mapcanvas->SaveAs((const char *)(savename));

	    //delete norm_g0;
	  	    
    }
   }
   A3->cd();
    if (r1_adc.size() > 1) {

      rangemin_g1 = lowestPointAfterSwitching(r1_filter, r0_filter);
      //cout << rangemin_g1 << endl;
      rangemin_g1 = checkRangeMinForRCEffect(rangemin_g1);
      rangemax_g1 = highestPointBeforeSwitching(r1_filter, r2_filter);

      //cout << rangemin_g1 << endl;
      //cout << rangemax_g1 << endl;

      if (rangemax_g1 > rangemin_g1) {

	fit_g1 = new TF1("fit_g1","[0]+[1]*x",rangemin_g1,rangemax_g1);
	fit_g1->SetParameter(0, 10000.);
	fit_g1->SetParameter(1, -0.1);
	fit_g1->SetLineColor(kGreen+2);
	fit_g1->SetParName(0,"G1 const");
	fit_g1->SetParName(1,"G1 grad");
	grap_g1->Fit(fit_g1,"QR","");
	TFitResultPtr r1 = grap_g1->Fit(fit_g1,"S");  
    par1 = r1->Parameter(1);      
	grap_g1->SetMinimum(0);
	grap_g1->SetMaximum(16000);
	grap_g1->GetXaxis()->SetTitle("Integration time [#mus]");
	grap_g1->GetYaxis()->SetTitle("ADC [ADU]");
	grap_g1->GetYaxis()->SetTitleOffset(0.9);
	grap_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
	grap_g1->Draw("AP");
	fit_g1->Draw("same");
	A3->Update();
	TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats");
	st0->SetX1NDC(0.2);
	st0->SetX2NDC(0.54);
	st0->SetY1NDC(0.17);
	st0->SetY2NDC(0.36);
	st0->SetBorderSize(0);
	st0->SetTextColor(kGreen+2);
	A3->SetLogx();
    A3->Modified();
    A3->Update();


	vector<double> r1_adc_norm;
	for (size_t j = 0; j < r1_adc.size(); j++) {
	r1_adc_norm.push_back(r1_adc[j] - fit_g1->Eval(r1_filter[j]));
	}

    A4->cd();
	TGraphErrors *norm_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcerr[0]));
	norm_g1->SetMarkerColor(kGreen+2);
	norm_g1->SetLineColor(kGreen+2);

	TF1* flat_g1 = new TF1("flat_g1","0",rangemin_g1,rangemax_g1);
	flat_g1->SetLineColor(kGreen+2);

	TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
	lin_g1_p02pc->SetParameter(0,fit_g1->GetParameter(0)/500.);
	lin_g1_p02pc->SetParameter(1,fit_g1->GetParameter(1)/500.);
	lin_g1_p02pc->SetLineColor(kOrange);
      
    TF1* lin_g1_p05pc = new TF1("lin_g1_p05pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
	lin_g1_p05pc->SetParameter(0,fit_g1->GetParameter(0)/200.);
    lin_g1_p05pc->SetParameter(1,fit_g1->GetParameter(1)/200.);
	lin_g1_p05pc->SetLineColor(kOrange+1);

	TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
	lin_g1_m02pc->SetParameter(0,fit_g1->GetParameter(0)/-500.);
	lin_g1_m02pc->SetParameter(1,fit_g1->GetParameter(1)/-500.);
	lin_g1_m02pc->SetLineColor(kOrange);

	TF1* lin_g1_m05pc = new TF1("lin_g1_m05pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
	lin_g1_m05pc->SetParameter(0,fit_g1->GetParameter(0)/-200.);
	lin_g1_m05pc->SetParameter(1,fit_g1->GetParameter(1)/-200.);
	lin_g1_m05pc->SetLineColor(kOrange+1);

	norm_g1->GetXaxis()->SetTitle("Integration time [#mus]");
	norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]");
	norm_g1->GetYaxis()->SetTitleOffset(0.9);
	norm_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
	norm_g1->SetMinimum(1.5*lin_g1_m05pc->Eval(rangemin_g1));
	norm_g1->SetMaximum(1.5*lin_g1_p05pc->Eval(rangemin_g1));
	norm_g1->Draw("AP");
	flat_g1->Draw("same");
	lin_g1_p02pc->Draw("same");
	lin_g1_p05pc->Draw("same");
	lin_g1_m02pc->Draw("same");
	lin_g1_m05pc->Draw("same");
	norm_g1->Draw("P");
    A4->Modified();
    A4->Update();

   }
}

A5->cd();

    if (r2_adc.size() > 1) {

      rangemin_g2 = lowestPointAfterSwitching(r2_filter, r1_filter);
      rangemin_g2 = checkRangeMinForRCEffect(rangemin_g2);

      rangemax_g2 = *max_element(r2_filter.begin(),r2_filter.end());

      if (rangemax_g2 > rangemin_g2) {

	    fit_g2 = new TF1("fit_g2","[0]+[1]*x",rangemin_g2, rangemax_g2);
	    fit_g2->SetParameter(0, 10000.);
	    fit_g2->SetParameter(1, -0.01);
	    fit_g2->SetLineColor(kRed);
	    fit_g2->SetParName(0,"G2 const");
	    fit_g2->SetParName(1,"G2 grad");
	    grap_g2->Fit(fit_g2,"QR+","");
      TFitResultPtr r = grap_g2->Fit(fit_g2,"S");  
      par2 = r->Parameter(1);     
	  grap_g2->SetMinimum(0);
	  grap_g2->SetMaximum(16000);
	  grap_g2->GetXaxis()->SetTitle("Integration time [#mus]");
	  grap_g2->GetYaxis()->SetTitle("ADC [ADU]");
	  grap_g2->GetYaxis()->SetTitleOffset(0.9);
	  grap_g2->GetXaxis()->SetLimits(*min_element(r2_filter.begin(),r2_filter.end()),*max_element(r2_filter.begin(),r2_filter.end()));
	  grap_g2->Draw("AP");
	  fit_g2->Draw("same");
	  A5->Update();
	  TPaveStats *st0 = (TPaveStats*)grap_g2->FindObject("stats");
	  st0->SetX1NDC(0.2);
	  st0->SetX2NDC(0.54);
	  st0->SetY1NDC(0.17);
	  st0->SetY2NDC(0.36);
	  st0->SetBorderSize(0);
	  st0->SetTextColor(kRed);
	  A5->SetLogx();
      A5->Modified();
      A5->Update();

	  vector<double> r2_adc_norm;
	  for (size_t j = 0; j < r2_adc.size(); j++) {
	    r2_adc_norm.push_back(r2_adc[j] - fit_g2->Eval(r2_filter[j]));
	  }

      A6->cd();

	  TGraphErrors *norm_g2 = new TGraphErrors(r2_adc.size(),&(r2_filter[0]),&(r2_adc_norm[0]),&(r2_ferr[0]),&(r2_adcerr[0]));
	  norm_g2->SetMarkerColor(kRed);
	  norm_g2->SetLineColor(kRed);
	  
	  TF1* flat_g2 = new TF1("flat_g2","0",rangemin_g2, rangemax_g2);	  
	  flat_g2->SetLineColor(kRed);

	  TF1* lin_g2_p02pc = new TF1("lin_g2_p02pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
	  lin_g2_p02pc->SetParameter(0,fit_g2->GetParameter(0)/500.);
	  lin_g2_p02pc->SetParameter(1,fit_g2->GetParameter(1)/500.);
	  lin_g2_p02pc->SetLineColor(kOrange);

	  TF1* lin_g2_p05pc = new TF1("lin_g2_p05pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
	  lin_g2_p05pc->SetParameter(0,fit_g2->GetParameter(0)/200.);
	  lin_g2_p05pc->SetParameter(1,fit_g2->GetParameter(1)/200.);
	  lin_g2_p05pc->SetLineColor(kOrange+1);

	  TF1* lin_g2_m02pc = new TF1("lin_g2_m02pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
	  lin_g2_m02pc->SetParameter(0,fit_g2->GetParameter(0)/-500.);
	  lin_g2_m02pc->SetParameter(1,fit_g2->GetParameter(1)/-500.);
	  lin_g2_m02pc->SetLineColor(kOrange);

	  TF1* lin_g2_m05pc = new TF1("lin_g2_m05pc","[0]+[1]*x",rangemin_g2,rangemax_g2);
	  lin_g2_m05pc->SetParameter(0,fit_g2->GetParameter(0)/-200.);
	  lin_g2_m05pc->SetParameter(1,fit_g2->GetParameter(1)/-200.);
	  lin_g2_m05pc->SetLineColor(kOrange+1);

	  norm_g2->GetXaxis()->SetTitle("Integration time [#mus]");
	  norm_g2->GetYaxis()->SetTitle("Normalised ADC [ADU]");
	  norm_g2->GetYaxis()->SetTitleOffset(0.9);
	  norm_g2->GetXaxis()->SetLimits(*min_element(r2_filter.begin(),r2_filter.end()),*max_element(r2_filter.begin(),r2_filter.end()));
	  norm_g2->SetMinimum(1.5*lin_g2_m05pc->Eval(rangemin_g2));
	  norm_g2->SetMaximum(1.5*lin_g2_p05pc->Eval(rangemin_g2));
	  norm_g2->Draw("AP");
	  flat_g2->Draw("same");
	  lin_g2_p02pc->Draw("same");
	  lin_g2_p05pc->Draw("same");
	  lin_g2_m02pc->Draw("same");
	  lin_g2_m05pc->Draw("same");
	  norm_g2->Draw("P");
      A6->Modified();
      A6->Update();
     }
    }
    // cout << par1/par2 << endl; 
    rootapp.Run();
    return 0;

}