1 | { |
---|
2 | TFile* fn = new TFile("Proba.root"); |
---|
3 | TTree* t1 = (TTree*)fn->Get("ProbaOscill"); |
---|
4 | |
---|
5 | |
---|
6 | TCanvas* c1 = new TCanvas("c1","Current",100,100,600,600); |
---|
7 | c1->SetGridx(); |
---|
8 | c1->SetGridy(); |
---|
9 | t1->Draw("Pmumu:LovE"); |
---|
10 | TGraph* gmumu = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gmumu")); |
---|
11 | t1->Draw("Pmue:LovE"); |
---|
12 | TGraph* gmue = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gmue")); |
---|
13 | |
---|
14 | t1->Draw("PMmumu:LovE"); |
---|
15 | TGraph* gMmumu = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gMmumu")); |
---|
16 | t1->Draw("PMmue:LovE"); |
---|
17 | TGraph* gMmue = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gMmue")); |
---|
18 | |
---|
19 | t1->Draw("Pbmumu:LovE"); |
---|
20 | TGraph* gbmumu = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gbmumu")); |
---|
21 | t1->Draw("Pbmue:LovE"); |
---|
22 | TGraph* gbmue = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gbmue")); |
---|
23 | |
---|
24 | t1->Draw("PMbmumu:LovE"); |
---|
25 | TGraph* gMbmumu = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gMbmumu")); |
---|
26 | t1->Draw("PMbmue:LovE"); |
---|
27 | TGraph* gMbmue = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gMbmue")); |
---|
28 | |
---|
29 | t1->Draw("(Pbmue-Pmue)/(Pbmue+Pmue):LovE"); |
---|
30 | TGraph* gAsymMuE = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gAsymMuE")); |
---|
31 | t1->Draw("(PMbmue-PMmue)/(PMbmue+PMmue):LovE"); |
---|
32 | TGraph* gMAsymMuE = (TGraph*)(gPad->GetPrimitive("Graph")->Clone("gMAsymMuE")); |
---|
33 | |
---|
34 | |
---|
35 | |
---|
36 | //color: 1:black, 2: red, 3:green, 4:blue, 5: yellow, 6:cyan, 7:light blue |
---|
37 | //Line Style 1=solid, 2=dash, 3=dot-dot, 4=dash-dot |
---|
38 | //Width in pixel |
---|
39 | |
---|
40 | //P mu->mu alone |
---|
41 | TCanvas* cmumu = new TCanvas("cmumu","Pmumu",100,100,600,600); |
---|
42 | cmumu->SetGridx(); |
---|
43 | cmumu->SetGridy(); |
---|
44 | |
---|
45 | gmumu->SetLineColor(4); //blue |
---|
46 | gmumu->SetTitle("P(#nu_{#mu}#rightarrow #nu_{#mu})"); |
---|
47 | gmumu->GetXaxis()->SetTitle("L/E [km/GeV]"); |
---|
48 | gmumu->Draw("AL"); |
---|
49 | gMmumu->SetLineColor(2); //red |
---|
50 | gMmumu->SetLineStyle(2); |
---|
51 | gMmumu->Draw("Csame"); |
---|
52 | cmumu->Update(); |
---|
53 | |
---|
54 | //P mub->mub alone |
---|
55 | TCanvas* cbmumu = new TCanvas("cbmumu","Pbmumu",100,100,600,600); |
---|
56 | cbmumu->SetGridx(); |
---|
57 | cbmumu->SetGridy(); |
---|
58 | |
---|
59 | gbmumu->SetLineColor(4); //blue |
---|
60 | gbmumu->SetTitle("P(#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{#mu})"); |
---|
61 | gbmumu->GetXaxis()->SetTitle("L/E [km/GeV]"); |
---|
62 | gbmumu->Draw("AL"); |
---|
63 | gMbmumu->SetLineColor(2); //red |
---|
64 | gMbmumu->SetLineStyle(2); |
---|
65 | gMbmumu->Draw("Csame"); |
---|
66 | cbmumu->Update(); |
---|
67 | |
---|
68 | //P mu->e alone |
---|
69 | TCanvas* cmue = new TCanvas("cmue","Pmue",100,100,600,600); |
---|
70 | cmue->SetGridx(); |
---|
71 | cmue->SetGridy(); |
---|
72 | |
---|
73 | gmue->SetLineColor(4); //blue |
---|
74 | gmue->SetTitle("P(#nu_{#mu}#rightarrow #nu_{e})"); |
---|
75 | gmue->GetXaxis()->SetTitle("L/E [km/GeV]"); |
---|
76 | gmue->Draw("AL"); |
---|
77 | gMmue->SetLineColor(2); //red |
---|
78 | gMmue->SetLineStyle(2); |
---|
79 | gMmue->Draw("Csame"); |
---|
80 | cmue->Update(); |
---|
81 | |
---|
82 | //P mub->eb alone |
---|
83 | TCanvas* cbmue = new TCanvas("cbmue","Pbmue",100,100,600,600); |
---|
84 | cbmue->SetGridx(); |
---|
85 | cbmue->SetGridy(); |
---|
86 | |
---|
87 | gbmue->SetLineColor(4); //blue |
---|
88 | gbmue->SetTitle("P(#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e})"); |
---|
89 | gbmue->GetXaxis()->SetTitle("L/E [km/GeV]"); |
---|
90 | gbmue->Draw("AL"); |
---|
91 | gMbmue->SetLineColor(2); //red |
---|
92 | gMbmue->SetLineStyle(2); |
---|
93 | gMbmue->Draw("Csame"); |
---|
94 | cbmue->Update(); |
---|
95 | |
---|
96 | //Assymetry |
---|
97 | //P mub->eb alone |
---|
98 | TCanvas* cAsymMuE = new TCanvas("cAsymMuE","PAsymMuE",100,100,600,600); |
---|
99 | cAsymMuE->SetGridx(); |
---|
100 | cAsymMuE->SetGridy(); |
---|
101 | |
---|
102 | gAsymMuE->SetLineColor(4); //blue |
---|
103 | gAsymMuE->SetTitle("P(#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e})-P(#nu_{#mu} #rightarrow #nu_{e})"); |
---|
104 | gAsymMuE->GetXaxis()->SetTitle("L/E [km/GeV]"); |
---|
105 | gAsymMuE->Draw("AL"); |
---|
106 | gMAsymMuE->SetLineColor(2); //red |
---|
107 | gMAsymMuE->SetLineStyle(2); |
---|
108 | gMAsymMuE->Draw("Csame"); |
---|
109 | cAsymMuE->Update(); |
---|
110 | |
---|
111 | |
---|
112 | } |
---|