00001 {
00002
00003 TCanvas *c1 = new TCanvas("Polarization","",800,800);
00004
00005
00006 TF1 *f1 = new TF1("f1","(3/4)*(1-x*x)",-1,1);
00007 TF1 *f2 = new TF1("f2","(3/8)*(1+x*x)",-1,1);
00008
00009 c1->Divide(2,1);
00010 c1->cd(1);
00011 f1->GetXaxis()->SetTitle("cos(#theta)");
00012 f1->GetXaxis()->CenterTitle();
00013
00014 f1->GetHistogram()->SetTitle("Pure longitudinal angular distribution");
00015 f1->Draw("C");
00016
00017 c1->cd(2);
00018 f2->GetXaxis()->SetTitle("cos(#theta)");
00019 f2->GetXaxis()->CenterTitle();
00020 f2->GetHistogram()->SetTitle("Pure transverse angular distribution");
00021 f2->Draw("C");
00022
00023 }