| 1 | {
|
|---|
| 2 |
|
|---|
| 3 | Int_t nLogSin22q13 = 21;
|
|---|
| 4 | Double_t logSin2Min = -4;
|
|---|
| 5 | Double_t logSin2Max = -2;
|
|---|
| 6 | Double_t logSin2BinWidth =
|
|---|
| 7 | (logSin2Max - logSin2Min)/double(nLogSin22q13-1.);
|
|---|
| 8 | Double_t logSin2EdgeMax = logSin2Max + logSin2BinWidth/2.;
|
|---|
| 9 | Double_t logSin2EdgeMin = logSin2Min - logSin2BinWidth/2.;
|
|---|
| 10 |
|
|---|
| 11 | Int_t nDeltaCP = 31;
|
|---|
| 12 | Double_t deltaMin = -180.; //degree
|
|---|
| 13 | Double_t deltaMax = 180.;
|
|---|
| 14 | Double_t deltaBinWidth = (deltaMax - deltaMin)/double(nDeltaCP-1.);
|
|---|
| 15 | Double_t deltaEdgeMax = deltaMax + deltaBinWidth/2.;
|
|---|
| 16 | Double_t deltaEdgeMin = deltaMin - deltaBinWidth/2.;
|
|---|
| 17 |
|
|---|
| 18 | TLegend* leg = new TLegend(0.6,0.7,0.9,0.9);
|
|---|
| 19 | leg->SetTextFont(gStyle->GetTextFont());
|
|---|
| 20 | leg->SetTextSize(0.07);
|
|---|
| 21 | leg->SetFillColor(10);
|
|---|
| 22 |
|
|---|
| 23 |
|
|---|
| 24 | TCanvas* c1 = new TCanvas("c1","Contour",100,100,600,600);
|
|---|
| 25 | // c1->SetGridx();
|
|---|
| 26 | // c1->SetGridy();
|
|---|
| 27 | c1->SetLogy();
|
|---|
| 28 |
|
|---|
| 29 | TH2D* hdum = new TH2D("hdum","",
|
|---|
| 30 | nDeltaCP,deltaMin,deltaMax,
|
|---|
| 31 | nLogSin22q13,TMath::Power(10.,logSin2Min),
|
|---|
| 32 | TMath::Power(10.,logSin2Max)
|
|---|
| 33 | );
|
|---|
| 34 | hdum->SetXTitle("#delta_{CP} (deg.)");
|
|---|
| 35 | hdum->SetYTitle("sin^{2}(2#theta_{13})");
|
|---|
| 36 | hdum->Draw();
|
|---|
| 37 |
|
|---|
| 38 | TPad *overlay = new TPad("overlay","",
|
|---|
| 39 | gStyle->GetPadLeftMargin(),
|
|---|
| 40 | gStyle->GetPadBottomMargin(),
|
|---|
| 41 | 1-gStyle->GetPadRightMargin(),
|
|---|
| 42 | 1-gStyle->GetPadTopMargin());
|
|---|
| 43 | overlay->SetFillStyle(4000);
|
|---|
| 44 | overlay->Draw();
|
|---|
| 45 | overlay->Range(deltaMin,logSin2Min,deltaMax,logSin2Max);
|
|---|
| 46 |
|
|---|
| 47 | TFile f1("./SPLSensi440kT.root");
|
|---|
| 48 | TTree* tree1 = (TTree*)f1.Get("SplGlb");
|
|---|
| 49 | TH2D* histo1 = new TH2D("histo1"," Theta vs Delta",
|
|---|
| 50 | nDeltaCP,deltaEdgeMin,deltaEdgeMax,
|
|---|
| 51 | nLogSin22q13,logSin2EdgeMin,logSin2EdgeMax
|
|---|
| 52 | );
|
|---|
| 53 |
|
|---|
| 54 | tree1->Project(histo1->GetName(),"TMath::Log10(pow(sin(2*theta13),2)):deltaCP/pi*180.","chi2");
|
|---|
| 55 | histo1->SetContour(1);
|
|---|
| 56 | histo1->SetContourLevel(0,4.605); //Chi2 for 2 dof and 90% CL
|
|---|
| 57 | histo1->SetLineColor(2); //red
|
|---|
| 58 | overlay->cd();
|
|---|
| 59 | histo1->Draw("CONT3 SAME");
|
|---|
| 60 | leg->AddEntry(histo1,"440kT","l");
|
|---|
| 61 |
|
|---|
| 62 |
|
|---|
| 63 | TFile f2("./SPLSensi330kT.root");
|
|---|
| 64 | TTree* tree2 = (TTree*)f2.Get("SplGlb");
|
|---|
| 65 | TH2D* histo2 = new TH2D("histo2"," Theta vs Delta",
|
|---|
| 66 | nDeltaCP,deltaEdgeMin,deltaEdgeMax,
|
|---|
| 67 | nLogSin22q13,logSin2EdgeMin,logSin2EdgeMax
|
|---|
| 68 | );
|
|---|
| 69 |
|
|---|
| 70 | tree2->Project(histo2->GetName(),"TMath::Log10(pow(sin(2*theta13),2)):deltaCP/pi*180.","chi2");
|
|---|
| 71 | histo2->SetContour(1);
|
|---|
| 72 | histo2->SetContourLevel(0,4.605); //Chi2 for 2 dof and 90% CL
|
|---|
| 73 | histo2->SetLineStyle(2); //dashed
|
|---|
| 74 | histo2->SetLineColor(2); //red
|
|---|
| 75 | overlay->cd();
|
|---|
| 76 | histo2->Draw("CONT3 SAME");
|
|---|
| 77 | leg->AddEntry(histo2,"330kT","l");
|
|---|
| 78 |
|
|---|
| 79 |
|
|---|
| 80 | TFile f3("./SPLSensi220kT.root");
|
|---|
| 81 | TTree* tree3 = (TTree*)f3.Get("SplGlb");
|
|---|
| 82 | TH2D* histo3 = new TH2D("histo3"," Theta vs Delta",
|
|---|
| 83 | nDeltaCP,deltaEdgeMin,deltaEdgeMax,
|
|---|
| 84 | nLogSin22q13,logSin2EdgeMin,logSin2EdgeMax
|
|---|
| 85 | );
|
|---|
| 86 |
|
|---|
| 87 | tree3->Project(histo3->GetName(),"TMath::Log10(pow(sin(2*theta13),2)):deltaCP/pi*180.","chi2");
|
|---|
| 88 | histo3->SetContour(1);
|
|---|
| 89 | histo3->SetContourLevel(0,4.605); //Chi2 for 2 dof and 90% CL
|
|---|
| 90 | histo3->SetLineStyle(3); //dotted
|
|---|
| 91 | histo3->SetLineColor(2); //red
|
|---|
| 92 | overlay->cd();
|
|---|
| 93 | histo3->Draw("CONT3 SAME");
|
|---|
| 94 | leg->AddEntry(histo3,"220kT","l");
|
|---|
| 95 |
|
|---|
| 96 |
|
|---|
| 97 |
|
|---|
| 98 |
|
|---|
| 99 |
|
|---|
| 100 | //Add text
|
|---|
| 101 |
|
|---|
| 102 | TLatex* txt1 = new TLatex(-150,-3.7,"90% CL");
|
|---|
| 103 | txt1->SetTextSize(0.07);
|
|---|
| 104 | txt1->Draw();
|
|---|
| 105 |
|
|---|
| 106 | TLatex* txt2 = new TLatex(-150,-2.3,"2y(+) & 8y(-)");
|
|---|
| 107 | txt2->SetTextSize(0.07);
|
|---|
| 108 | txt2->Draw();
|
|---|
| 109 |
|
|---|
| 110 |
|
|---|
| 111 | TLatex* txt3 = new TLatex(60,-3.7,"#delta^{tr}_{CP}=#theta^{tr}_{13}=0");
|
|---|
| 112 | txt3->SetTextSize(0.07);
|
|---|
| 113 | txt3->Draw();
|
|---|
| 114 |
|
|---|
| 115 |
|
|---|
| 116 |
|
|---|
| 117 | leg->Draw("SAME");
|
|---|
| 118 | overlay->Update();
|
|---|
| 119 |
|
|---|
| 120 | }
|
|---|
