source: GLBFrejus/HEAD/run/Proba.C@ 700

Last change on this file since 700 was 154, checked in by campagne, 19 years ago

new version
File size: 3.4 KB
Line 
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}
Note: See TracBrowser for help on using the repository browser.