00001 {
00002
00003 TCanvas *c1 = new TCanvas("Probability","",1000,1000);
00004
00005
00006
00007
00008
00009
00010
00011 TF1 *f1 = new TF1("f1","(3/4)*(1-x*x)",-1,1);
00012 TF1 *f2 = new TF1("f2","(3/8)*(1+x*x)",-1,1);
00013 TF1 *f3 = new TF1("f3",".5*(3/4)*(1-x*x)+.5*(3/8)*(1+x*x)",-1,1);
00014
00015
00016 c1->Divide(3,1);
00017 c1->cd(1);
00018 f1->GetXaxis()->SetTitle("cos(#theta)");
00019 f1->GetXaxis()->CenterTitle();
00020
00021 f1->GetHistogram()->SetTitle("Pure longitudinal angular distribution");
00022 f1->Draw("");
00023
00024 c1->cd(2);
00025
00026 f2->GetXaxis()->SetTitle("cos(#theta)");
00027 f2->GetXaxis()->CenterTitle();
00028
00029 f2->GetHistogram()->SetTitle("Pure transverse angular distribution");
00030 f2->Draw("");
00031
00032 c1->cd(3);
00033
00034 f3->GetXaxis()->SetTitle("cos(#theta)");
00035 f3->GetXaxis()->CenterTitle();
00036
00037 f3->GetHistogram()->SetTitle("#alpha = 0.5, #beta = 0.5");
00038 f3->Draw("");
00039
00040
00041 }