source: GLBFrejus/HEAD/run/newComp.C @ 154

{
  TCanvas* c1 = new TCanvas("c1","Contour",100,100,600,600);
  c1->SetGridx();
  c1->SetGridy();
  
  //Reference plot
  TFile* file1 = new TFile("./Mauro/Original/SPLHierTOctTv2.root");
  TTree* tree1 = (TTree*)file1->Get("SplGlb");
  tree1->Draw("deltaCP:sin(2*theta13)^2>>h1(41,-0.001,0.081,31,-1.025,0.525)","chi2","SURF");
  TH2F *h1 = (TH2F*)gDirectory->Get("h1");
  h1->SetXTitle("sin^{2}(2#theta_{13})");
  h1->SetYTitle("#delta_{CP}/#pi");
  h1->SetContour(1);
  h1->SetContourLevel(0,4.605);
  h1->SetLineColor(2); //red
  h1->Draw("CONT3");

//   TH2D h1bis; 
//   h1->Copy(h1bis);
//   h1bis.SetContour(1);
//   h1bis.SetContourLevel(0,6.0); //Chi2 for 2 dof and 95% CL
//   h1bis.SetLineColor(2); //red
//   h1bis.SetLineStyle(2); //dashed
//   h1bis.Draw("CONT3 SAME");


  //Second to be compared
  TFile* file2 = new TFile("./SPLHierTOctFv2.root");
  TTree* tree2 = (TTree*)file2->Get("SplGlb");
  tree2->Draw("deltaCP:sin(2*theta13)^2>>h2(41,-0.001,0.081,31,-1.025,0.525)","chi2","SURF");
  TH2F *h2 = (TH2F*)gDirectory->Get("h2");
  h2->SetContour(1);
  h2->SetContourLevel(0,4.605);
  h2->SetLineStyle(1); //full
  h2->SetLineColor(4); //blue
  h2->Draw("CONT3 SAME");

//   TH2D h2bis; 
//   h2->Copy(h2bis);
//   h2bis.SetContour(1);
//   h2bis.SetContourLevel(0,6.0); //Chi2 for 2 dof and 95% CL
//   h2bis.SetLineColor(4); //blue
//   h2bis.SetLineStyle(2); //dashed
//   h2bis.Draw("CONT3 SAME");


}