[154] | 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 | } |
---|