void compare_color(){


  double x_ref[7]    =  {2, 4, 5, 9, 10, 11, 12};
  double amp_ref[7]  =  {2872, 2568.3, 2564.6, 988.4, 707.6, 861.3, 701};
  double res_ref[7]  =  {4.17, 4.26, 3.92, 4.50, 4.85, 4.54, 4.63};
  
  
  double amp_col[7]  =  {2848.5, 2234.0, 2047.5, 864.5, 695.7, 623.5, 534.3};
  //  double res_col[7]  =  {3.9, 4.57, 4.79, 4.82, 4.64, 5.63, 5.37};
  double res_col[7]  =  {4., 4.57, 4.5, 4.82, 4.64, 5.63, 5.37};


  double amp_test3[7]  =  {2808.67, 2127.94, 2245.1,  919.802, 723.213, 780.707, 673.918};
  double res_test3[7]  =  {3.83721, 4.39402, 4.26688, 4.9626, 4.55655,  4.88945, 4.65384};

  double amp_test4[7]  =  {2838.55, 2237.19, 2317.88,  849.364, 680.937, 783.534, 681.358};
  double res_test4[7]  =  {3.9, 4.4, 4.72,  4.5, 4.81211, 4.57362, 5.0};

  double amp_test5[7]  =  {2469.38, 2132.05, 2168.52, 790.596, 668.812, 733.728, 632.361};
  double res_test5[7]  =  {3.79504, 3.98778, 3.99093, 4.48286, 3.73557, 4.29 , 4.50209};

  double amp_test6[7]  =  {2689.1, 2273.23, 2519.57, 923.645, 680.283, 761.797, 701.635};
  double res_test6[7]  =  {4.1, 3.90489, 3.75641, 4.52417, 4.61937, 4.53736, 4.65003};

  double amp_test7[7]  =  {2722.16, 2367.72, 2242.15, 0, 685.778, 794.98, 666.679 };
  double res_test7[7]  =  {4.32867, 4.2914, 4.44, 0, 4.58676, 4.63805, 4.82898};

  double amp_test8[7]  =  {2841.99,  2334.67, 2401.59, 919.152, 720.675, 869.219, 692.408};
  double res_test8[7]  =  {3.94167, 4.53, 4.33016, 4.67, 4.54614, 4.45315, 4.74};


  double ex_ref[7];
  
  memset(ex_ref,0,sizeof(ex_ref));


  double rel_amp[7];

  double rel_amp_test3[7];
  double rel_amp_test4[7];
  double rel_amp_test5[7];
  double rel_amp_test6[7];
  double rel_amp_test7[7];
  double rel_amp_test8[7];

  for(int ii = 0; ii < 7; ii++){

    rel_amp[ii] = (amp_col[ii] - amp_ref[ii])/amp_ref[ii];
    rel_amp_test3[ii] = (amp_test3[ii] - amp_ref[ii])/amp_ref[ii];
    rel_amp_test4[ii] = (amp_test4[ii] - amp_ref[ii])/amp_ref[ii];

    rel_amp_test5[ii] = (amp_test5[ii] - amp_ref[ii])/amp_ref[ii];
    rel_amp_test6[ii] = (amp_test6[ii] - amp_ref[ii])/amp_ref[ii];
    rel_amp_test7[ii] = (amp_test7[ii] - amp_ref[ii])/amp_ref[ii];

    rel_amp_test8[ii] = (amp_test8[ii] - amp_ref[ii])/amp_ref[ii];

  }



  TGraphErrors *gr10 = new TGraphErrors(7, x_ref, rel_amp, ex_ref, ex_ref);

  TGraphErrors *gr11 = new TGraphErrors(7, x_ref, rel_amp_test3, ex_ref, ex_ref);

  TGraphErrors *gr12 = new TGraphErrors(7, x_ref, rel_amp_test4, ex_ref, ex_ref);

  TGraphErrors *gr15 = new TGraphErrors(7, x_ref, rel_amp_test5, ex_ref, ex_ref);

  TGraphErrors *gr16 = new TGraphErrors(7, x_ref, rel_amp_test6, ex_ref, ex_ref);

  TGraphErrors *gr17 = new TGraphErrors(7, x_ref, rel_amp_test7, ex_ref, ex_ref);

  TGraphErrors *gr18 = new TGraphErrors(7, x_ref, rel_amp_test8, ex_ref, ex_ref);



  TCanvas *c1 = new TCanvas("c1","c1", 400, 10, 600, 400);

  gPad->SetGrid(1);

  gr10->SetMarkerStyle(20);
  gr10->SetMarkerColor(4);

  gr10->SetMaximum(0.1);
  gr10->SetMinimum(-0.3);
  
  gr10->GetXaxis()->SetTitle("Module ID");
  gr10->GetYaxis()->SetTitle("(A_{COL} - A_{REF})/A_{REF}");


  gr10->SetTitle("");

  gr10->Draw("AP");

  gr11->SetMarkerStyle(21);
  gr11->SetMarkerColor(2);
 
  gr11->Draw("Psame");

  gr12->SetMarkerStyle(22);
  gr12->SetMarkerColor(6);

  gr12->Draw("Psame");

  gr15->SetMarkerStyle(23);
  gr15->SetMarkerColor(1);

  gr15->Draw("Psame");

  gr16->SetMarkerStyle(24);
  gr16->SetMarkerColor(7);

  gr16->Draw("Psame");

  gr17->SetMarkerStyle(25);
  gr17->SetMarkerColor(8);

  gr17->Draw("Psame");

  gr18->SetMarkerStyle(26);
  gr18->SetMarkerColor(9);

  gr18->Draw("Psame");


  leg1 = new TLegend(0.2,0.65,0.5,0.85);
  leg1->SetFillColor(0);
  leg1->SetTextSize(0.035);
  
  leg1->AddEntry(gr10,"Colored crystals ","p");
  leg1->AddEntry(gr11,"Good crystals, Test 3","p");
  leg1->AddEntry(gr12,"Good crystals, Test 4","p");
  leg1->AddEntry(gr15,"CRYTUR, Test 5","p");

  leg1->AddEntry(gr16,"CRYTUR, Test 6","p");
  leg1->AddEntry(gr17,"CRYTUR, Test 7","p");

  leg1->AddEntry(gr18,"CRYTUR, Test 8","p");

  leg1->Draw();



  TGraphErrors *gr20 = new TGraphErrors(7, x_ref, res_ref, ex_ref, ex_ref);
  TGraphErrors *gr21 = new TGraphErrors(7, x_ref, res_col, ex_ref, ex_ref);
  TGraphErrors *gr22 = new TGraphErrors(7, x_ref, res_test3, ex_ref, ex_ref);
  TGraphErrors *gr23 = new TGraphErrors(7, x_ref, res_test4, ex_ref, ex_ref);

  TGraphErrors *gr25 = new TGraphErrors(7, x_ref, res_test5, ex_ref, ex_ref);
  TGraphErrors *gr26 = new TGraphErrors(7, x_ref, res_test6, ex_ref, ex_ref);
  TGraphErrors *gr27 = new TGraphErrors(7, x_ref, res_test7, ex_ref, ex_ref);

  TGraphErrors *gr28 = new TGraphErrors(7, x_ref, res_test8, ex_ref, ex_ref);

  TCanvas *c2 = new TCanvas("c2","c2", 400, 10, 600, 400);

  gPad->SetGrid(1);

  gr20->SetMarkerStyle(22);
  gr20->SetMarkerColor(6);

  gr20->SetMaximum(6.);
  gr20->SetMinimum(3.);
  
  gr20->GetXaxis()->SetTitle("Module ID");
  gr20->GetYaxis()->SetTitle("#sigma_{Amp}/A   (%)");


  gr20->SetTitle("");

  gr20->Draw("AP");

  gr21->SetMarkerStyle(20);
  gr21->SetMarkerColor(4);
  
  gr21->Draw("Psame");

  gr22->SetMarkerStyle(21);
  gr22->SetMarkerColor(2);
  
  gr22->Draw("Psame");


  gr23->SetMarkerStyle(22);
  gr23->SetMarkerColor(6);
  
  gr23->Draw("Psame");

  gr25->SetMarkerStyle(23);
  gr25->SetMarkerColor(1);
  
  gr25->Draw("Psame");

  gr26->SetMarkerStyle(24);
  gr26->SetMarkerColor(7);
  
  gr26->Draw("Psame");

  gr27->SetMarkerStyle(25);
  gr27->SetMarkerColor(8);
 
  gr27->Draw("Psame");

  gr28->SetMarkerStyle(26);
  gr28->SetMarkerColor(9);
 
  gr28->Draw("Psame");

  leg = new TLegend(0.2,0.65,0.5,0.85);
  leg->SetFillColor(0);
  leg->SetTextSize(0.035);
  
  leg->AddEntry(gr20,"Good crystals, Test 1  ","p");
  leg->AddEntry(gr21,"Colored crystals ","p");
  leg->AddEntry(gr22,"Good crystals, Test 3 ","p");
  leg->AddEntry(gr23,"Good crystals, Test 4 ","p");
  leg->AddEntry(gr25,"CRYTUR, Test 5 ","p");
  leg->AddEntry(gr26,"CRYTUR, Test 6 ","p");
  leg->AddEntry(gr27,"CRYTUR, Test 7 ","p");
  leg->AddEntry(gr28,"CRYTUR, Test 8 ","p");

  leg->Draw();

}