Visulaization macro for ExB effect (Marian)
[u/mrichter/AliRoot.git] / TPC / macros / AliTPCExBdraw.C
1
2 void Draw(Double_t theta){
3   //
4   // Draw  "primary track distortion"
5   //
6   //Double_t theta=0.5;
7   TF1 *fdrfi_rpi0=new TF1("fdrfi_rpi0",Form("AliTPCExB::GetDrphi(x,0*pi/2,%f*x)",theta),20,250);
8   TF1 *fdrfi_rpi1=new TF1("fdrfi_rpi1",Form("AliTPCExB::GetDrphi(x,1*pi/2,%f*x)",theta),20,250);
9   TF1 *fdrfi_rpi2=new TF1("fdrfi_rpi2",Form("AliTPCExB::GetDrphi(x,2*pi/2,%f*x)",theta),20,250);
10   TF1 *fdrfi_rpi3=new TF1("fdrfi_rpi3",Form("AliTPCExB::GetDrphi(x,3*pi/2,%f*x)",theta),20,250);
11
12   fdrfi_rpi0->GetXaxis()->SetTitle("r (cm)");
13   fdrfi_rpi0->GetYaxis()->SetTitle("drd#phi (cm)");
14   
15   fdrfi_rpi0->SetLineColor(2);
16   fdrfi_rpi1->SetLineColor(3);
17   fdrfi_rpi2->SetLineColor(4);
18   fdrfi_rpi3->SetLineColor(5);
19   
20   fdrfi_rpi0->Draw();
21   fdrfi_rpi1->Draw("same");
22   fdrfi_rpi2->Draw("same");
23   fdrfi_rpi3->Draw("same");
24
25   TLegend *legend = new TLegend(0.15,0.70,0.6,0.85, Form("ExB distortion alog track (tan(#theta)=%f)",theta));
26   legend->AddEntry(fdrfi_rpi0, "#phi=0", "l");
27   legend->AddEntry(fdrfi_rpi1, "#phi=#pi/2", "l");
28   legend->AddEntry(fdrfi_rpi2, "#phi=pi", "l");
29   legend->AddEntry(fdrfi_rpi3, "#phi=3*#pi/2", "l");
30   legend->Draw();
31 }