{
 #define DisplayTitle 0

 #define Constrain 1
 Int_t n=7;
 TCanvas *c=new TCanvas("c","",10,10,1000,600);
 TFile *f=new TFile("RootFiles/Response_ecal.root","RECREATE");   
 gStyle->SetOptFit();
 gStyle->SetOptStat(0);
 Double_t E[n]={5,10,15,20,50,75,100}, par[3];
 Int_t N[n]={77497,45361,28876,49716,13794,12370};

 Double_t E_Mean1[n]={125.68,214.17,356.12,496.4,1352.61,2067.54,3068.15};
 Double_t errE_Mean1[n]={1.95,1.7,2.53,2.93,12.08,52.08,52.08};
 Double_t ChiSq1[n]={2.691,0.7392,0.308,0.9878,0.7756,0.204,0.03};

 Double_t E_Mean2[n]={118.43,215.47,359.58,500.81,1349.5,2019.48,3052.16};
 Double_t errE_Mean2[n]={0.3,1.97,2.24,2.86,10.78,37.96,47.94};
 Double_t ChiSq2[n]={4.3226,0.759,0.42,0.6718,0.4543,0.4475,0.2775};

 Double_t Err1[n],Err2[n];
 
 Double_t errWeight=0.00;
 Double_t errE[n]={0,0,0,0,0,0,0};
 Double_t errFitted[n];

   gStyle->SetOptFit(111);
   c->SetFillColor(kWhite);
   c->SetGrid();
   c->GetFrame()->SetBorderSize(1);



 for(Int_t i=0;i<n;i++) {
    Err1[i] = errE_Mean1[i] * TMath::Sqrt(ChiSq1[i]);
    Err1[i] *= 0.001;
//    E_Mean1[i] *= 0.001 / E[i];
    E_Mean1[i] *= 0.001;
    Err2[i] = errE_Mean2[i] * TMath::Sqrt(ChiSq2[i]);
    Err2[i] *= 0.001;
//    E_Mean2[i] *= 0.001 / E[i];
    E_Mean2[i] *= 0.001;
    cout << E[i] << "\t" << E_Mean1[i] << "\t"<< Err1[i] << "\t" << E_Mean2[i] << "\t"<< Err2[i] << endl;
 }

/* for(Int_t i=0;i<n;i++) {
  E_Mean[i]=E_Mean2[i];
  errE_Mean1[i]=errE_Mean1[i]*sqrt(ChiSq1[i]);
  errE_Mean2[i]=errE_Mean2[i]*sqrt(ChiSq2[i]);
  Err[i]=(E_Mean1[i]-E_Mean2[i])*(E_Mean1[i]-E_Mean2[i])+errE_Mean1[i]*errE_Mean1[i]+errE_Mean2[i]*errE_Mean2[i]+errWeight*errWeight;  
  Err[i]=sqrt(Err[i]);
  printf("Mean %lf (MeV) Err %lf (MeV)\n" ,E_Mean1[i],Err[i]);
  Err[i]=0.001*Err[i];  
  E_Mean[i]=E_Mean[i]*0.001;   
 }*/

 TGraph *gr=new TGraphErrors(7,E,E_Mean1,errE,Err1);
 gr->Draw("P");
 gr->SetLineColor(2);
 gr->SetMarkerStyle(23);
 gr->SetMarkerColor(2);
 gr->SetMarkerSize(1.0);

 TGraph *gr1=new TGraphErrors(7,E,E_Mean2,errE,Err2);
 gr1->Draw("P");
 gr1->SetLineColor(4);
 gr1->SetMarkerStyle(22);
 gr1->SetMarkerColor(4);
 gr1->SetMarkerSize(1.0);
 
   TMultiGraph *mg = new TMultiGraph();
   mg->Add(gr);
   mg->Add(gr1);
   mg->Draw("AP");
   mg->SetTitle("");
   mg->GetXaxis()->SetTitle("E_{#pi^{-}} GeV");
   mg->GetXaxis()->CenterTitle(kFALSE);
   mg->GetXaxis()->SetLabelSize(0.02);
//   mg->GetYaxis()->SetTitle("#LTE_{live}#GT #/ E_{#pi^{-}}");
   mg->GetYaxis()->SetTitle("#LTE_{live}#GT GeV");
   mg->GetYaxis()->CenterTitle(kFALSE);
   mg->GetYaxis()->SetLabelSize(0.02);

   gr->Fit("pol1");
   TF1 *fline1=gr->GetFunction("pol1");
   if(Constrain)fline1->SetParameter(0,0);
   par[0] = fline1->GetParameter(0);
   par[1] = fline1->GetParameter(1);
   cout << par[0] << " " << par[1] << endl;
   fline1->SetLineColor(kRed);

   gr1->Fit("pol1");
   TF1 *fline2=gr1->GetFunction("pol1");
   if(Constrain)fline2->SetParameter(0,0);
   par[0] = fline2->GetParameter(0);
   par[1] = fline2->GetParameter(1);
   cout << par[0] << " " << par[1] << endl;
   fline2->SetLineColor(kBlue);

//  c->Print("./GIF/Response_ecal.gif");
//  c->Write();
}