Changeset 154 for GLBFrejus


Ignore:
Timestamp:
May 30, 2006, 11:32:37 AM (18 years ago)
Author:
campagne
Message:

new version

Location:
GLBFrejus/HEAD/run
Files:
12 added
7 edited

Legend:

Unmodified
Added
Removed

  • GLBFrejus/HEAD/run/discoCP.C

    r2 r154  
    99  Double_t logSin2EdgeMin = logSin2Min - logSin2BinWidth/2.;
    1010
    11   Int_t nDeltaCP  = 21;
     11  Int_t nDeltaCP  = 41;
    1212  Double_t deltaMin = 0.; //rad/pi
    13   Double_t deltaMax =  1.;
     13  Double_t deltaMax =  2.;
    1414  Double_t deltaBinWidth = (deltaMax - deltaMin)/double(nDeltaCP-1.);
    1515  Double_t deltaEdgeMax = deltaMax + deltaBinWidth/2.;
    … …  
    4646
    4747
    48   TFile f1("./SPLCPDiscovery2p100-h0-o0.root");
     48  TFile f1("./SPLCPDisco-h0-o0.root");
    4949  TTree* tree1 = (TTree*)f1.Get("SplGlb");
    5050  TH2D* histo1 = new TH2D("histo1"," Theta vs Delta",
    … …  
    6060  overlay->cd();
    6161  histo1->Draw("CONT3 SAME");
    62   leg->AddEntry(histo1,"syst. 2\%","l");
     62  leg->AddEntry(histo1,"H^{t} Us","l");
    6363
    6464
    65   TPaveLabel* text = new TPaveLabel(-3.6,1.,-2.5,1.1,"true sin^{2}(#theta_{23})=0.5");
     65  TPaveLabel* text = new TPaveLabel(-3.0,1.,-2.5,1.1,"true sin^{2}#theta_{23}=0.5");
    6666  text->SetBorderSize(0);
    6767  text->SetFillColor(10);
    … …  
    7171
    7272
    73   TFile f2("./SPLCPDiscovery3p100-h0-o0.root");
     73  TFile f2("./SPLCPDiscoRefereeA-h0-o0.root");
    7474  TTree* tree2 = (TTree*)f2.Get("SplGlb");
    7575  TH2D* histo2 = new TH2D("histo2"," Theta vs Delta",
    … …  
    8181  histo2->SetContour(1);
    8282  histo2->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
    83   histo2->SetLineColor(5); //yellow
     83  histo2->SetLineColor(2); //red
    8484  overlay->cd();
    8585  histo2->Draw("CONT3 SAME");
    86   leg->AddEntry(histo2,"syst. 3\%","l");
     86  leg->AddEntry(histo2,"H^{t} Ref.A","l");
    8787
    8888
    89   TFile f3("./SPLCPDiscovery5p100-h0-o0.root");
     89  TFile f3("./SPLCPDisco-h1-o0.root");
    9090  TTree* tree3 = (TTree*)f3.Get("SplGlb");
    9191  TH2D* histo3 = new TH2D("histo3"," Theta vs Delta",
    … …  
    9797  histo3->SetContour(1);
    9898  histo3->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
    99   //  histo3->SetLineStyle(2); //dashed
    100   histo3->SetLineColor(3); //green
     99  histo3->SetLineStyle(2); //dashed
     100  histo3->SetLineColor(4); //blue
    101101  overlay->cd();
    102102  histo3->Draw("CONT3 SAME");
    103   leg->AddEntry(histo3,"syst. 5%","l");
     103  leg->AddEntry(histo3,"H^{w} Us","l");
    104104
    105105
    106   TFile f4("./SPLCPDiscovery10p100-h0-o0.root");
     106  TFile f4("./SPLCPDiscoRefereeA-h1-o0.root");
    107107  TTree* tree4 = (TTree*)f4.Get("SplGlb");
    108108  TH2D* histo4 = new TH2D("histo4"," Theta vs Delta",
    … …  
    114114  histo4->SetContour(1);
    115115  histo4->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
    116   //  histo4->SetLineStyle(2); //dashed
     116  histo4->SetLineStyle(2); //dashed
    117117  histo4->SetLineColor(2); //red
    118118  overlay->cd();
    119119  histo4->Draw("CONT3 SAME");
    120   leg->AddEntry(histo4,"syst. 10%","l");
     120  leg->AddEntry(histo4,"H^{w} Ref. A","l");
    121121
    122   TFile f5("./SPLCPDiscoveryFake-h1-o0.root");
    123   TTree* tree5 = (TTree*)f5.Get("SplGlb");
    124   TH2D* histo5 = new TH2D("histo5"," Theta vs Delta",
    125                           nLogSin22q13,logSin2EdgeMin,logSin2EdgeMax,
    126                           nDeltaCP,deltaEdgeMin,deltaEdgeMax
    127                           );
    128   tree5->Project(histo5->GetName(),"deltaCP/pi:TMath::Log10(pow(sin(2*theta13),2))","chi2");
     122//   TFile f5("./SPLCPDiscoveryFake-h1-o0.root");
     123//   TTree* tree5 = (TTree*)f5.Get("SplGlb");
     124//   TH2D* histo5 = new TH2D("histo5"," Theta vs Delta",
     125//                        nLogSin22q13,logSin2EdgeMin,logSin2EdgeMax,
     126//                        nDeltaCP,deltaEdgeMin,deltaEdgeMax
     127//                        );
     128//   tree5->Project(histo5->GetName(),"deltaCP/pi:TMath::Log10(pow(sin(2*theta13),2))","chi2");
    129129
    130   histo5->SetContour(1);
    131   histo5->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
    132   histo5->SetLineStyle(2); //dashed
    133   histo5->SetLineColor(4); //red
    134   overlay->cd();
    135   histo5->Draw("CONT3 SAME");
    136   leg->AddEntry(histo5,"syst 2%, Fake Hier","l");
     130//   histo5->SetContour(1);
     131//   histo5->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
     132//   histo5->SetLineStyle(2); //dashed
     133//   histo5->SetLineColor(4); //red
     134//   overlay->cd();
     135//   histo5->Draw("CONT3 SAME");
     136//   leg->AddEntry(histo5,"syst 2%, Fake Hier","l");
    137137
    138138

  • GLBFrejus/HEAD/run/discoCP.data

    r2 r154  
    11WRONG-OCTANT-OPTION     0
    22WRONG-HIERAR-OPTION     1
    3 THETA12-RADIAN          0.594816801
    4 THETA23-RADIAN          0.684719203
     3THETA12-RADIAN          0.579639740
     4THETA23-RADIAN          0.785398163
    55DMQ21-EV2               7.9e-5
    66DMQ31-EV2               2.4e-3
    7 ROOT-OUTPUT-FILE        SPLCPDiscoveryFake
    8 NBIN-DELTA              21
     7ROOT-OUTPUT-FILE        SPLCPDisco
     8NBIN-DELTA              41
    99NBIN-THETA13            21
    1010MIN-LOG10(SIN2(2T13))   -3.3

  • GLBFrejus/HEAD/run/discoTheta13.C

    r2 r154  
    4747
    4848
    49   TFile f1("./SPLThetaDiscovery-h0-o0.root");
     49  TFile f1("./SPLThetaDiscoTest-h0-o0.root");
    5050  TTree* tree1 = (TTree*)f1.Get("SplGlb");
    5151  TH2D* histo1 = new TH2D("histo1"," Theta vs Delta",
    … …  
    6262  histo1->Draw("CONT3 SAME");
    6363  //  leg->AddEntry(histo1,"H^{tr} O^{tr}","l");
    64   leg->AddEntry(histo1,"True Octant","l");
     64  leg->AddEntry(histo1,"440kT","l");
    6565
    6666
    … …  
    9494
    9595
    96   TFile f3("./SPLThetaDiscovery-h0-o1.root");
    97   TTree* tree3 = (TTree*)f3.Get("SplGlb");
    98   TH2D* histo3 = new TH2D("histo3"," Theta vs Delta",
    99                           nLogSin22q13,logSin2EdgeMin,logSin2EdgeMax,
    100                           nDeltaCP,deltaEdgeMin,deltaEdgeMax
    101                           );
    102   tree3->Project(histo3->GetName(),"deltaCP/pi:TMath::Log10(pow(sin(2*theta13),2))","chi2");
     96//   TFile f3("./SPLThetaDiscovery220kT-h0-o0.root");
     97//   TTree* tree3 = (TTree*)f3.Get("SplGlb");
     98//   TH2D* histo3 = new TH2D("histo3"," Theta vs Delta",
     99//                        nLogSin22q13,logSin2EdgeMin,logSin2EdgeMax,
     100//                        nDeltaCP,deltaEdgeMin,deltaEdgeMax
     101//                        );
     102//   tree3->Project(histo3->GetName(),"deltaCP/pi:TMath::Log10(pow(sin(2*theta13),2))","chi2");
    103103
    104   histo3->SetContour(1);
    105   histo3->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
    106   histo3->SetLineStyle(2); //dashed 
    107   histo3->SetLineColor(4); //blue
    108   overlay->cd();
    109   histo3->Draw("CONT3 SAME");
    110   leg->AddEntry(histo3,"Fake Octant","l");
    111   //  leg->AddEntry(histo3,"H^{tr} O^{wr}","l");
     104//   histo3->SetContour(1);
     105//   histo3->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
     106//   histo3->SetLineStyle(2); //dashed 
     107//   histo3->SetLineColor(4); //blue
     108//   overlay->cd();
     109//   histo3->Draw("CONT3 SAME");
     110//   leg->AddEntry(histo3,"220kT","l");
     111//   //  leg->AddEntry(histo3,"H^{tr} O^{wr}","l");
    112112
    113113

  • GLBFrejus/HEAD/run/discoTheta13.data

    r2 r154  
    1 WRONG-OCTANT-OPTION     1
     1WRONG-OCTANT-OPTION     0
    22WRONG-HIERAR-OPTION     0
    3 THETA12-RADIAN          0.594816801
     3THETA12-RADIAN          0.579639740
    44THETA23-RADIAN          0.684719203
    55DMQ21-EV2               7.9e-5
    66DMQ31-EV2               2.4e-3
    7 ROOT-OUTPUT-FILE        SPLThetaDiscovery
     7ROOT-OUTPUT-FILE        SPLThetaDiscoTest
    88NBIN-DELTA              31
    99NBIN-THETA13            21

  • GLBFrejus/HEAD/run/dm31th23.C

    r2 r154  
    9090
    9191  histo1->SetContour(1);
    92   histo1->SetContourLevel(0,1.);  //Chi2 for 1 dof and 1sigma
    93   //histo1->SetContourLevel(0,4); //Chi2 for 1 dof and 2sigma
    94   //histo1->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
     92  histo1->SetContourLevel(0,9.2);  //Chi2 for 2 dof and 99%
    9593  histo1->SetLineColor(2); //red
    9694  overlay->cd();
    … …  
    9997
    10098
    101 //   TH2D* h1bis = histo1->Clone("h1bis");
    102 //   h1bis->SetContour(1);
    103 //   h1bis->SetContourLevel(0,9); //Chi2 for 1 dof and 3sigma
    104 //   overlay->cd();
    105 //   h1bis->SetLineColor(4); //blue
    106 //   h1bis->Draw("CONT3 SAME");
    107 //   leg->AddEntry(h1bis,"1#sigma 5yrs (+)","l");
    10899
    109100

  • GLBFrejus/HEAD/run/dm31th23.data

    r2 r154  
    1 THETA12-RADIAN          0.594816801
    2 THETA23-RADIAN          0.725253222
     1THETA12-RADIAN          0.579639740
     2THETA23-RADIAN          0.785398163
    33THETA13-RADIAN          0.
    44DELTA-RADIAN            0.
    55DMQ21-EV2               7.9e-5
    6 DMQ31-EV2               2.4e-3
    7 ROOT-OUTPUT-FILE        SPLDmTheta23MeasZoom2
    8 NBIN-DM31               21
    9 NBIN-THETA23            31
    10 MIN-DM31                2.3e-3
    11 MAX-DM31                2.5e-3
    12 MIN-SIN2TH23            0.41
    13 MAX-SIN2TH23            0.48
     6DMQ31-EV2               2.2e-3
     7ROOT-OUTPUT-FILE        SPLDmTheta23Meas35GeV
     8NBIN-DM31               41
     9NBIN-THETA23            51
     10MIN-DM31                2.0e-3
     11MAX-DM31                3.0e-3
     12MIN-SIN2TH23            0.3
     13MAX-SIN2TH23            0.7

  • GLBFrejus/HEAD/run/sensi.C

    r2 r154  
    11{
    22
    3   gROOT->LoadMacro("TDSyst.C++");
     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
    424  TCanvas* c1 = new TCanvas("c1","Contour",100,100,600,600);
    5   c1->SetGridx();
    6   c1->SetGridy();
     25  //  c1->SetGridx();
     26  // c1->SetGridy();
    727  c1->SetLogy();
    828
    9   TH2D* hdum = new TH2D("hdum","Theta vs Delta",41,-180,180,31,1e-4,1e-1);
    10   hdum->SetXTitle("#delta_{CP} (deg)");
     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.)");
    1135  hdum->SetYTitle("sin^{2}(2#theta_{13})");
    1236  hdum->Draw();
    … …  
    1943  overlay->SetFillStyle(4000);
    2044  overlay->Draw();
    21   overlay->Range(-180.,-4,180.,-1);
     45  overlay->Range(deltaMin,logSin2Min,deltaMax,logSin2Max);
    2246
    23   //Original
    24   TFile* file1 = new TFile("./SPLDeltaTheta2+8OldFluxCorr.root");
    25   TTree* tree1 = (TTree*)file1->Get("SplGlb");
    26   TDSyst* comp1 = new TDSyst(tree1,"h1");
    27   comp1->Loop();
     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                          );
    2853
    29   TH2D* h1 = comp1->GetHisto();
    30   h1->SetContour(1);
    31   h1->SetContourLevel(0,4.605); //Chi2 for 2 dof and 90% CL
    32   //  h1->SetLineStyle(2);
    33   h1->SetLineColor(4); //blue
     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
    3458  overlay->cd();
    35   h1->Draw("CONT3 SAME");
    36   overlay->Update();
     59  histo1->Draw("CONT3 SAME");
     60  leg->AddEntry(histo1,"440kT","l");
    3761
    38   //Comp
    39   TFile* file2 = new TFile("./OptiFlux/SPLDeltaThetaSens.root");
    40   TTree* tree2 = (TTree*)file2->Get("SplGlb");
    41   TDSyst* comp2 = new TDSyst(tree2,"h2");
    42   comp2->Loop();
    4362
    44   TH2D* h2 = comp2->GetHisto();
    45   h2->SetContour(1);
    46   h2->SetContourLevel(0,4.605); //Chi2 for 2 dof and 90% CL
    47   h2->SetLineColor(2); //red
     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
    4875  overlay->cd();
    49   h2->Draw("CONT3 SAME");
    50   overlay->Update();
     76  histo2->Draw("CONT3 SAME");
     77  leg->AddEntry(histo2,"330kT","l");
    5178
    52   //Comp
    53   TFile* file3 = new TFile("./OptiFlux/SPLDeltaThetaRateCorr.root");
    54   TTree* tree3 = (TTree*)file3->Get("SplGlb");
    55   TDSyst* comp3 = new TDSyst(tree3,"h3");
    56   comp3->Loop();
    5779
    58   TH2D* h3 = comp3->GetHisto();
    59   h3->SetContour(1);
    60   h3->SetContourLevel(0,4.605); //Chi2 for 2 dof and 90% CL
    61   h3->SetLineColor(6); //blue
     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
    6292  overlay->cd();
    63   h3->Draw("CONT3 SAME");
     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");
    64118  overlay->Update();
    65119
Note: See TracChangeset for help on using the changeset viewer.