{
   gStyle->SetOptFit(111111);
   gStyle->SetOptStat(0);

   TCanvas *c0=new TCanvas("c0","plots Vana - VANA ",900,600);

   ifstream infile_chip1("IVsf.txt",ios::in);         

   // numb of events
   Int_t n=10;

   Int_t col1[n];

   //Variables
   Float_t col2[n];
   Float_t col3[n];
   Float_t col4[n];
   Float_t col5[n];
   Float_t col6[n];
   Float_t col7[n];
   Float_t col8[n];

   for (int i=0; i<n; i++){
      col1[i]=0;
      col2[i]=0.;
      col3[i]=0.;
      col4[i]=0.;
      col5[i]=0.;
      col6[i]=0.;
      col7[i]=0.;
      col8[i]=0.;
   }

   if (!infile_chip1){
      cout << "The data file is missing!!!" << endl;
      exit(1);
   }

   // Get the data
   while (!infile_chip1.eof()){
      for (int j=0; j<n ;j++){
         infile_chip1 >> col1[j] >> col2[j] >> col3[j]>>col4[j]>>col5[j]>>col6[j]>>col7[j]>>col8[j];   
      }
   }
  
   Double_t VSF[n];
   Double_t Vsf[n];
   Double_t Isf_0[n]; 
   Double_t Isf_50[n];
   Double_t Isf_100[n];
   Double_t Isf_150[n];
   Double_t Isf_200[n];
   Double_t Isf_250[n];

   for (int i=0; i<n; i++){    
      VSF[i]=0;
      Vsf[i]=0.; 
      Isf_0[i]=0.;
      Isf_50[i]=0.;
      Isf_100[i]=0.;
      Isf_150[i]=0.;
      Isf_200[i]=0.;
      Isf_250[i]=0.;
   } 

   for (int i=0; i<n; i++){
      VSF[i]=col1[i];
      Vsf[i]=col2[i];
      Isf_0[i]=col3[i];
      Isf_50[i]=col4[i];
      Isf_100[i]=col5[i];
      Isf_150[i]=col6[i];
      Isf_200[i]=col7[i];
      Isf_250[i]=col8[i];
   } 

   infile_chip1.close();            
 
   Double_t VSF_dac[n];
   Double_t vsf[n];
   Double_t isf_0[n];
   Double_t isf_50[n];
   Double_t isf_100[n];
   Double_t isf_150[n];
   Double_t isf_200[n];
   Double_t isf_250[n];

   TH2F *VSF_vsf= new TH2F("","Vsf - VSF DAC;VSF[DAC];Vsf [V]",256,0,256,200,1,2.5);
   TH2F *ivsf_vana_0= new TH2F("","IV sf curve: VANA;VSF[DAC];Isf[mA]",256,0,260,200,-1,90);
   TH2F *ivsf_vana_50= new TH2F("","IV sf curve: VANA;VSF[DAC];Iana[mA]",256,0,260,200,-1,90);
   TH2F *ivsf_vana_100= new TH2F("","IV sf curve: VANA;VSF[DAC];Iana[mA]",256,0,260,200,-1,90);
   TH2F *ivsf_vana_150= new TH2F("","IV sf curve: VANA;VSF[DAC];Iana[mA]",256,0,260,200,-1,90);
   TH2F *ivsf_vana_200= new TH2F("","IV sf curve: VANA;VSF[DAC];Iana[mA]",256,0,260,200,-1,90);
   TH2F *ivsf_vana_250= new TH2F("","IV sf curve: VANA;VSF[DAC];Iana[mA]",256,0,260,200,-1,90);

   Int_t n1=10;

   for(int i=0;i<n1;i++){
      VSF_dac[i]=VSF[i];
      vsf[i]=Vsf[i];
      cout<<VSF_dac[i]<<"\t"<<vsf[i]<<endl;
      VSF_vsf->Fill(VSF_dac[i],vsf[i]);
   }
 
   for(int i=0;i<n1;i++){
      VSF_dac[i]=VSF[i];
      isf_0[i]=Isf_0[i];
      cout<<VSF_dac[i]<<"\t"<<isf_0[i]<<endl;
      ivsf_vana_0->Fill(VSF_dac[i],isf_0[i]);
   }

   for(int i=0;i<n1;i++){
      VSF_dac[i]=VSF[i];
      isf_50[i]=Isf_50[i];
      cout<<VSF_dac[i]<<"\t"<<isf_50[i]<<endl;
      ivsf_vana_50->Fill(VSF_dac[i],isf_50[i]);
   }

 
   for(int i=0;i<n1;i++){
      VSF_dac[i]=VSF[i];
      isf_100[i]=Isf_100[i];
      cout<<VSF_dac[i]<<"\t"<<isf_100[i]<<endl;
      ivsf_vana_100->Fill(VSF_dac[i],isf_100[i]);
   }

 
   for(int i=0;i<n1;i++){
      VSF_dac[i]=VSF[i];
      isf_150[i]=Isf_150[i];
      cout<<VSF_dac[i]<<"\t"<<isf_150[i]<<endl;
      ivsf_vana_150->Fill(VSF_dac[i],isf_150[i]);
   }

   for(int i=0;i<n1;i++){
      VSF_dac[i]=VSF[i];
      isf_200[i]=Isf_200[i];
      cout<<VSF_dac[i]<<"\t"<<isf_200[i]<<endl;
      ivsf_vana_200->Fill(VSF_dac[i],isf_200[i]);
   }

 
   for(int i=0;i<n1;i++){
     VSF_dac[i]=VSF[i];
     isf_250[i]=Isf_250[i];
     cout<<VSF_dac[i]<<"\t"<<isf_250[i]<<endl;
     ivsf_vana_250->Fill(VSF_dac[i],isf_250[i]);
   }

   
   VSF_vsf->SetMarkerStyle(21);
   VSF_vsf->SetMarkerSize(0.5);
   VSF_vsf->SetMarkerColor(2);
   VSF_vsf->Draw();

   TF1 *fit1= new TF1("fit1","[0]+ x*[1]", 0,210);
   VSF_vsf->Fit("fit1","R");
}