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

Last change on this file since 154 was 154, checked in by campagne, 18 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.