]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGPP/macros/plotSparseL.C
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Split: removed dirs now in AliPhysics
[u/mrichter/AliRoot.git] / PWGPP / macros / plotSparseL.C
diff --git a/PWGPP/macros/plotSparseL.C b/PWGPP/macros/plotSparseL.C
deleted file mode 100644 (file)
index f3a0cc1..0000000
--- a/PWGPP/macros/plotSparseL.C
+++ /dev/null
@@ -1,444 +0,0 @@
-{
-
-//TFile f("output/pp-trunk010709-pdc1-00XXX.root")
-//TFile f("results_V0QA_v416rev08_090709.root");
-TFile f("output/pp-v417-release.root");
-
-THnSparseF* sparseL;
-
-sparseL=(THnSparseF*)(outputTRD->FindObject("sparseL"));
-
-sparseL->GetAxis(5)->SetRangeUser(0.,3.9);
-sparseL->GetAxis(9)->SetRangeUser(0.,3.9);
-
-
-TH1* ptgeantPPhys=sparseL->Projection(5);
-TH1* etageantPPhys=sparseL->Projection(6);
-TH1* phigeantPPhys=sparseL->Projection(7);
-
-TH1* ptResSinglePPhys=sparseL->Projection(13);
-TH1* etasinglePPhys=sparseL->Projection(14);
-TH1* phisinglePPhys=sparseL->Projection(15);
-
-
-TH1* ptgeantNPhys=sparseL->Projection(9);
-TH1* etageantNPhys=sparseL->Projection(10);
-TH1* phigeantNPhys=sparseL->Projection(11);
-
-TH1* ptResSingleNPhys=sparseL->Projection(21);
-TH1* etasingleNPhys=sparseL->Projection(22);
-TH1* phisingleNPhys=sparseL->Projection(23);
-
-
-TH1* efiEtaSingPPhys=(TH1*) etageantPPhys->Clone();
-efiEtaSingPPhys->Sumw2();
-efiEtaSingPPhys->Divide(etasinglePPhys,etageantPPhys);
-efiEtaSingPPhys->SetMinimum(0.);
-efiEtaSingPPhys->SetMaximum(1.2);
-efiEtaSingPPhys->SetAxisRange(-1.5,1.5);
-efiEtaSingPPhys->SetMarkerStyle(20);
-efiEtaSingPPhys->SetMarkerColor(2);
-
-
-TH1* efiEtaSingNPhys=(TH1*) etageantNPhys->Clone();
-efiEtaSingNPhys->Sumw2();
-efiEtaSingNPhys->Divide(etasingleNPhys,etageantNPhys);
-efiEtaSingNPhys->SetMinimum(0.);
-efiEtaSingNPhys->SetMaximum(1.2);
-efiEtaSingNPhys->SetAxisRange(-1.5,1.5);
-efiEtaSingNPhys->SetMarkerStyle(20);
-efiEtaSingNPhys->SetMarkerColor(2);
-
-
-TH1* efiPhiSingPPhys=(TH1*) phigeantPPhys->Clone();
-efiPhiSingPPhys->Sumw2();
-efiPhiSingPPhys->Divide(phisinglePPhys,phigeantPPhys);
-efiPhiSingPPhys->SetMinimum(0.);
-efiPhiSingPPhys->SetMaximum(1.2);
-efiPhiSingPPhys->SetAxisRange(-TMath::Pi(),TMath::Pi());
-efiPhiSingPPhys->SetMarkerStyle(20);
-efiPhiSingPPhys->SetMarkerColor(2);
-
-
-TH1* efiPhiSingNPhys=(TH1*) phigeantNPhys->Clone();
-efiPhiSingNPhys->Sumw2();
-efiPhiSingNPhys->Divide(phisingleNPhys,phigeantNPhys);
-efiPhiSingNPhys->SetMinimum(0.);
-efiPhiSingNPhys->SetMaximum(1.2);
-efiPhiSingNPhys->SetAxisRange(-TMath::Pi(),TMath::Pi());
-efiPhiSingNPhys->SetMarkerStyle(20);
-efiPhiSingNPhys->SetMarkerColor(2);
-
-// reconstructable tracks
-
-// track Lenght
-sparseL->GetAxis(8)->SetRangeUser(50,200);
-sparseL->GetAxis(12)->SetRangeUser(50,200);
-// momentum
-sparseL->GetAxis(5)->SetRangeUser(0.,3.9);
-sparseL->GetAxis(9)->SetRangeUser(0.,3.9);
-
-
-TH1* rgeantTrk=sparseL->Projection(3);
-TH1* ptgeantPTrk=sparseL->Projection(5);
-TH1* etageantPTrk=sparseL->Projection(6);
-TH1* phigeantPTrk=sparseL->Projection(7);
-
-
-TH1* ptgeantNTrk=sparseL->Projection(9);
-TH1* etageantNTrk=sparseL->Projection(10);
-TH1* phigeantNTrk=sparseL->Projection(11);
-
-sparseL->GetAxis(20)->SetRangeUser(0,2.);
-sparseL->GetAxis(28)->SetRangeUser(0,2.);
-
-TH1* rsingleTrk=sparseL->Projection(3);
-TH1* ptsinglePTrk=sparseL->Projection(5);
-TH1* etasinglePTrk=sparseL->Projection(6);
-TH1* phisinglePTrk=sparseL->Projection(7);
-
-
-
-
-TH1* ptsingleNTrk=sparseL->Projection(9);
-TH1* etasingleNTrk=sparseL->Projection(10);
-TH1* phisingleNTrk=sparseL->Projection(11);
-
-
-sparseL->GetAxis(32)->SetRangeUser(0,2.);
-sparseL->GetAxis(36)->SetRangeUser(0,2.);
-
-TH1* rV0Trk=sparseL->Projection(3);
-TH1* ptV0PTrk=sparseL->Projection(5);
-TH1* etaV0PTrk=sparseL->Projection(6);
-TH1* phiV0PTrk=sparseL->Projection(7);
-
-
-
-TH1* ptV0NTrk=sparseL->Projection(9);
-TH1* etaV0NTrk=sparseL->Projection(10);
-TH1* phiV0NTrk=sparseL->Projection(11);
-
-
-TH1* efiSingPPhys=(TH1*) ptgeantPPhys->Clone();
-efiSingPPhys->Sumw2();
-efiSingPPhys->Divide(ptgeantPTrk,ptgeantPPhys);
-efiSingPPhys->SetAxisRange(0.,5);
-efiSingPPhys->SetMinimum(0.);
-efiSingPPhys->SetMaximum(1.2);
-efiSingPPhys->SetMarkerStyle(20);
-efiSingPPhys->SetMarkerColor(2);
-
-
-TH1* efiSingNPhys=(TH1*) ptgeantNPhys->Clone();
-efiSingNPhys->Sumw2();
-efiSingNPhys->Divide(ptgeantNTrk,ptgeantNPhys);
-efiSingNPhys->SetAxisRange(0.,5);
-efiSingNPhys->SetMinimum(0.);
-efiSingNPhys->SetMaximum(1.2);
-efiSingNPhys->SetMarkerStyle(20);
-efiSingNPhys->SetMarkerColor(2);
-
-
-TH1* efiSingPTrk=(TH1*) ptgeantPTrk->Clone();
-efiSingPTrk->Sumw2();
-efiSingPTrk->Divide(ptsinglePTrk,ptgeantPTrk);
-efiSingPTrk->SetAxisRange(0.,5);
-efiSingPTrk->SetMinimum(0.);
-efiSingPTrk->SetMaximum(1.2);
-efiSingPTrk->SetMarkerStyle(21);
-efiSingPTrk->SetMarkerColor(4);
-
-TH1* efiSingNTrk=(TH1*) ptgeantNTrk->Clone();
-efiSingNTrk->Sumw2();
-efiSingNTrk->Divide(ptsingleNTrk,ptgeantNTrk);
-efiSingNTrk->SetAxisRange(0.,5);
-efiSingNTrk->SetMinimum(0.);
-efiSingNTrk->SetMaximum(1.2);
-efiSingNTrk->SetMarkerStyle(21);
-efiSingNTrk->SetMarkerColor(4);
-
-
-TH1* efiV0PTrk=(TH1*) ptsinglePTrk->Clone();
-efiV0PTrk->Sumw2();
-efiV0PTrk->Divide(ptV0PTrk,ptsinglePTrk);
-efiV0PTrk->SetAxisRange(0.,5);
-efiV0PTrk->SetMinimum(0.);
-efiV0PTrk->SetMaximum(1.2);
-efiV0PTrk->SetMarkerStyle(21);
-efiV0PTrk->SetMarkerColor(3);
-
-TH1* efiV0NTrk=(TH1*) ptsingleNTrk->Clone();
-efiV0NTrk->Sumw2();
-efiV0NTrk->Divide(ptV0NTrk,ptsingleNTrk);
-efiV0NTrk->SetAxisRange(0.,5);
-efiV0NTrk->SetMinimum(0.);
-efiV0NTrk->SetMaximum(1.2);
-efiV0NTrk->SetMarkerStyle(21);
-efiV0NTrk->SetMarkerColor(3);
-
-
-TH1* efiEtaSingPTrk=(TH1*) etageantPTrk->Clone();
-efiEtaSingPTrk->Sumw2();
-efiEtaSingPTrk->Divide(etasinglePTrk,etageantPTrk);
-efiEtaSingPTrk->SetMinimum(0.);
-efiEtaSingPTrk->SetMaximum(1.2);
-efiEtaSingPTrk->SetAxisRange(-1.5,1.5);
-efiEtaSingPTrk->SetMarkerStyle(21);
-efiEtaSingPTrk->SetMarkerColor(4);
-
-
-TH1* efiEtaSingNTrk=(TH1*) etageantNTrk->Clone();
-efiEtaSingNTrk->Sumw2();
-efiEtaSingNTrk->Divide(etasingleNTrk,etageantNTrk);
-efiEtaSingNTrk->SetMinimum(0.);
-efiEtaSingNTrk->SetMaximum(1.2);
-efiEtaSingNTrk->SetAxisRange(-1.5,1.5);
-efiEtaSingNTrk->SetMarkerStyle(21);
-efiEtaSingNTrk->SetMarkerColor(4);
-
-
-TH1* efiEtaV0PTrk=(TH1*) etasinglePTrk->Clone();
-efiEtaV0PTrk->Sumw2();
-efiEtaV0PTrk->Divide(etaV0PTrk,etasinglePTrk);
-efiEtaV0PTrk->SetMinimum(0.);
-efiEtaV0PTrk->SetMaximum(1.2);
-efiEtaV0PTrk->SetMarkerStyle(21);
-efiEtaV0PTrk->SetMarkerColor(3);
-
-
-TH1* efiEtaV0NTrk=(TH1*) etasingleNTrk->Clone();
-efiEtaV0NTrk->Sumw2();
-efiEtaV0NTrk->Divide(etaV0NTrk,etasingleNTrk);
-efiEtaV0NTrk->SetMinimum(0.);
-efiEtaV0NTrk->SetMaximum(1.2);
-efiEtaV0NTrk->SetMarkerStyle(21);
-efiEtaV0NTrk->SetMarkerColor(3);
-
-
-
-TH1* efiPhiSingPTrk=(TH1*) phigeantPTrk->Clone();
-efiPhiSingPTrk->Sumw2();
-efiPhiSingPTrk->Divide(phisinglePTrk,phigeantPTrk);
-efiPhiSingPTrk->SetMinimum(0.);
-efiPhiSingPTrk->SetMaximum(1.2);
-efiPhiSingPTrk->SetAxisRange(-TMath::Pi(),TMath::Pi());
-efiPhiSingPTrk->SetMarkerStyle(20);
-efiPhiSingPTrk->SetMarkerColor(4);
-
-
-TH1* efiPhiSingNTrk=(TH1*) phigeantNTrk->Clone();
-efiPhiSingNTrk->Sumw2();
-efiPhiSingNTrk->Divide(phisingleNTrk,phigeantNTrk);
-efiPhiSingNTrk->SetMinimum(0.);
-efiPhiSingNTrk->SetMaximum(1.2);
-efiPhiSingNTrk->SetAxisRange(-TMath::Pi(),TMath::Pi());
-efiPhiSingNTrk->SetMarkerStyle(20);
-efiPhiSingNTrk->SetMarkerColor(4);
-
-
-TH1* efiPhiV0PTrk=(TH1*) phisinglePTrk->Clone();
-efiPhiV0PTrk->Sumw2();
-efiPhiV0PTrk->Divide(phiV0PTrk,phisinglePTrk);
-efiPhiV0PTrk->SetMinimum(0.);
-efiPhiV0PTrk->SetMaximum(1.2);
-efiPhiV0PTrk->SetMarkerStyle(21);
-efiPhiV0PTrk->SetMarkerColor(3);
-
-
-TH1* efiPhiV0NTrk=(TH1*) phisingleNTrk->Clone();
-efiPhiV0NTrk->Sumw2();
-efiPhiV0NTrk->Divide(phiV0NTrk,phisingleNTrk);
-efiPhiV0NTrk->SetMinimum(0.);
-efiPhiV0NTrk->SetMaximum(1.2);
-efiPhiV0NTrk->SetMarkerStyle(21);
-efiPhiV0NTrk->SetMarkerColor(3);
-
-
-
-sparseL->GetAxis(32)->SetRangeUser(-2,0.);
-sparseL->GetAxis(36)->SetRangeUser(-2,0.);
-
-
-
-TH1* sigbXYZEPlusRecSingle=sparseL->Projection(17);
-TH1* sigbXYZEMinusRecSingle=sparseL->Projection(25);
-
-TH1* ptsinglePTrkNo=sparseL->Projection(5);
-TH1* etasinglePTrkNo=sparseL->Projection(6);
-TH1* phisinglePTrkNo=sparseL->Projection(7);
-
-TH1* ptsingleNTrkNo=sparseL->Projection(9);
-TH1* etasingleNTrkNo=sparseL->Projection(10);
-TH1* phisingleNTrkNo=sparseL->Projection(11);
-
-
-
-
-
-TH1* efiV0PTrkNo=(TH1*) ptsinglePTrkNo->Clone();
-efiV0PTrkNo->Sumw2();
-efiV0PTrkNo->Divide(ptsinglePTrkNo,ptsinglePTrk);
-efiV0PTrkNo->SetAxisRange(0.,5);
-efiV0PTrkNo->SetMinimum(0.);
-efiV0PTrkNo->SetMaximum(1.2);
-efiV0PTrkNo->SetMarkerStyle(21);
-efiV0PTrkNo->SetMarkerColor(6);
-
-
-
-TH1* efiV0NTrkNo=(TH1*) ptsingleNTrkNo->Clone();
-efiV0NTrkNo->Sumw2();
-efiV0NTrkNo->Divide(ptsingleNTrkNo,ptsingleNTrk);
-efiV0NTrkNo->SetAxisRange(0.,5);
-efiV0NTrkNo->SetMinimum(0.);
-efiV0NTrkNo->SetMaximum(1.2);
-efiV0NTrkNo->SetMarkerStyle(21);
-efiV0NTrkNo->SetMarkerColor(6);
-
-
-TH1* efiEtaV0PTrkNo=(TH1*) etasinglePTrkNo->Clone();
-efiEtaV0PTrkNo->Sumw2();
-efiEtaV0PTrkNo->Divide(etasinglePTrkNo,etasinglePTrk);
-efiEtaV0PTrkNo->SetAxisRange(-1.5,1.5);
-efiEtaV0PTrkNo->SetMinimum(0.);
-efiEtaV0PTrkNo->SetMaximum(1.2);
-efiEtaV0PTrkNo->SetMarkerStyle(21);
-efiEtaV0PTrkNo->SetMarkerColor(6);
-
-
-
-TH1* efiEtaV0NTrkNo=(TH1*) etasingleNTrkNo->Clone();
-efiEtaV0NTrkNo->Sumw2();
-efiEtaV0NTrkNo->Divide(etasingleNTrkNo,etasingleNTrk);
-efiEtaV0NTrkNo->SetAxisRange(-1.5,1.5);
-efiEtaV0NTrkNo->SetMinimum(0.);
-efiEtaV0NTrkNo->SetMaximum(1.2);
-efiEtaV0NTrkNo->SetMarkerStyle(21);
-efiEtaV0NTrkNo->SetMarkerColor(6);
-
-
-TH1* efiPhiV0PTrkNo=(TH1*) phisinglePTrkNo->Clone();
-efiPhiV0PTrkNo->Sumw2();
-efiPhiV0PTrkNo->Divide(phisinglePTrkNo,phisinglePTrk);
-efiPhiV0PTrkNo->SetAxisRange(-3.5,3.5);
-efiPhiV0PTrkNo->SetMinimum(0.);
-efiPhiV0PTrkNo->SetMaximum(1.2);
-efiPhiV0PTrkNo->SetMarkerStyle(21);
-efiPhiV0PTrkNo->SetMarkerColor(6);
-
-
-
-TH1* efiPhiV0NTrkNo=(TH1*) phisingleNTrkNo->Clone();
-efiPhiV0NTrkNo->Sumw2();
-efiPhiV0NTrkNo->Divide(phisingleNTrkNo,phisingleNTrk);
-efiPhiV0NTrkNo->SetAxisRange(-3.5,3.5);
-efiPhiV0NTrkNo->SetMinimum(0.);
-efiPhiV0NTrkNo->SetMaximum(1.2);
-efiPhiV0NTrkNo->SetMarkerStyle(21);
-efiPhiV0NTrkNo->SetMarkerColor(6);
-
-
-
-
-
- TCanvas* c1 = new TCanvas("c1","MC Info",200,10,700,900);
-  c1->Divide(2,4);
-  c1->SetFillColor(0);
-  c1->GetFrame()->SetFillColor(0);
-  c1->SetBorderMode(0);
- TPad *c1_1 = (TPad*)(c1->GetPrimitive("c1_1"));
- TPad *c1_2 = (TPad*)(c1->GetPrimitive("c1_2"));
- TPad *c1_3 = (TPad*)(c1->GetPrimitive("c1_3"));
- TPad *c1_4 = (TPad*)(c1->GetPrimitive("c1_4"));
- TPad *c1_5 = (TPad*)(c1->GetPrimitive("c1_5"));
- TPad *c1_6 = (TPad*)(c1->GetPrimitive("c1_6"));
-
-
- c1_1->SetGridx(1);
- c1_3->SetGridx(1);
- c1_2->SetGridx(1);
- c1_4->SetGridx(1);
- c1_5->SetGridx(1);
- c1_6->SetGridx(1);
-
- c1_1->SetGridy(1);
- c1_3->SetGridy(1);
- c1_2->SetGridy(1);
- c1_4->SetGridy(1);
- c1_5->SetGridy(1);
- c1_6->SetGridy(1);
-
-
-c1->cd(1);
-efiSingPTrk->Draw();
-efiSingPTrk->SetXTitle("1/sqrt(pt) pos (sqrt(GeV))");
-efiSingPPhys->Draw("same");
-efiV0PTrk->Draw("same");
-efiV0PTrkNo->Draw("same");
-
-TLegend *leg1 = new TLegend(0.6,0.75,0.9,0.9);
-leg1->AddEntry(efiSingPTrk,"Single tracking eff");
-leg1->AddEntry(efiV0PTrk,"V0 tracking eff");
-leg1->AddEntry(efiV0PTrkNo,"V0 tracking ineff");
-leg1->Draw();
-c1->cd(2);
-efiSingNTrk->Draw();
-efiSingNTrk->SetXTitle("1/sqrt(pt) neg (sqrt(GeV))");
-efiSingNPhys->Draw("same");
-efiV0NTrk->Draw("same");
-efiV0NTrkNo->Draw("same");
-
-c1->cd(3);
-efiEtaSingPTrk->Draw();
-//efiEtaSingPPhys->Draw("same");
-efiEtaV0PTrk->Draw("same");
-efiEtaV0PTrkNo->Draw("same");
-
-c1->cd(4);
-efiEtaSingNTrk->Draw();
-//efiEtaSingNPhys->Draw("same");
-efiEtaV0NTrk->Draw("same");
-efiEtaV0NTrkNo->Draw("same");
-
-c1->cd(5);
-efiPhiSingPTrk->Draw();
-//efiPhiSingPPhys->Draw("same");
-efiPhiV0PTrk->Draw("same");
-efiPhiV0PTrkNo->Draw("same");
-
-c1->cd(6);
-efiPhiSingNTrk->Draw();
-//efiPhiSingNPhys->Draw("same");
-efiPhiV0NTrk->Draw("same");
-efiPhiV0NTrkNo->Draw("same");
-
-c1->cd(7);
-TH1* efiRSingTrk=(TH1*) rgeantTrk->Clone();
-efiRSingTrk->Sumw2();
-efiRSingTrk->Divide(rsingleTrk,rgeantTrk);
-
-
-efiRSingTrk->SetMinimum(0.);
-efiRSingTrk->SetMaximum(1.2);
-efiRSingTrk->SetMarkerStyle(21);
-efiRSingTrk->SetMarkerColor(4);
-
-efiRSingTrk->Draw();
-
-TH1* efiRV0Trk=(TH1*) rgeantTrk->Clone();
-efiRV0Trk->Sumw2();
-efiRV0Trk->Divide(rV0Trk,rsingleTrk);
-
-
-efiRV0Trk->SetMinimum(0.);
-efiRV0Trk->SetMaximum(1.2);
-efiRV0Trk->SetMarkerStyle(21);
-efiRV0Trk->SetMarkerColor(3);
-efiRV0Trk->Draw("same");
-
-
-//c1->Print("Effi-L-norefit-v416rev08-090709.ps"); 
-//c1->Print("Effi-L-norefit-trunk010709.ps"); 
-c1->Print("Effi-L-norefit-v417release-110709.ps"); 
-}