-//-----------------------------------------------------------------------------
-void AliD0toKpi::DrawPIDinTOF(TString pidScheme) const {
- // Draw parameterized PID probabilities in TOF
-
- const char *tofparampbpb = strstr(pidScheme.Data(),"TOFparam_PbPb");
- const char *tofparampp = strstr(pidScheme.Data(),"TOFparam_pp");
-
- TH2F* framePi = new TH2F("framePi","Tag probabilities for PIONS",2,0,2.5,2,0,1);
- framePi->SetXTitle("p [GeV/c]");
- framePi->SetStats(0);
- TH2F* frameK = new TH2F("frameK","Tag probabilities for KAONS",2,0,2.5,2,0,1);
- frameK->SetXTitle("p [GeV/c]");
- frameK->SetStats(0);
- TH2F* frameP = new TH2F("frameP","Tag probabilities for PROTONS",2,0,4.5,2,0,1);
- frameP->SetXTitle("p [GeV/c]");
- frameP->SetStats(0);
-
- TH1F* hPiPi = new TH1F("hPiPi","Tag probabilities for PIONS",kPiBins_PbPb,0,2.5);
- TH1F* hPiNid = new TH1F("hPiNid","Tag probabilities for PIONS",kPiBins_PbPb,0,2.5);
-
- TH1F* hKK = new TH1F("hKK","Tag probabilities for KAONS",kKBins_PbPb,0,2.5);
- TH1F* hKNid = new TH1F("hKNid","Tag probabilities for KAONS",kKBins_PbPb,0,2.5);
- TH1F* hKPi = new TH1F("hKPi","Tag probabilities for KAONS",kKBins_PbPb,0,2.5);
-
- TH1F* hPP = new TH1F("hPP","Tag probabilities for PROTONS",kPBins_PbPb,0,4.5);
- TH1F* hPNid = new TH1F("hPNid","Tag probabilities for PROTONS",kPBins_PbPb,0,4.5);
- TH1F* hPPi = new TH1F("hPPi","Tag probabilities for PROTONS",kPBins_PbPb,0,4.5);
-
-
- if(tofparampbpb) {
-
- for(Int_t i=1; i<=kPiBins_PbPb; i++) {
- hPiPi->SetBinContent(i,kPiTagPi_PbPb[i-1]);
- hPiNid->SetBinContent(i,kPiTagPi_PbPb[i-1]+kPiTagNid_PbPb[i-1]);
-
- hKK->SetBinContent(i,kKTagK_PbPb[i-1]);
- hKPi->SetBinContent(i,kKTagK_PbPb[i-1]+kKTagPi_PbPb[i-1]);
- hKNid->SetBinContent(i,kKTagK_PbPb[i-1]+kKTagPi_PbPb[i-1]+kKTagNid_PbPb[i-1]);
- }
- for(Int_t i=1; i<=kPBins_PbPb; i++) {
- hPP->SetBinContent(i,kPTagP_PbPb[i-1]);
- hPPi->SetBinContent(i,kPTagP_PbPb[i-1]+kPTagPi_PbPb[i-1]);
- hPNid->SetBinContent(i,kPTagP_PbPb[i-1]+kPTagPi_PbPb[i-1]+kPTagNid_PbPb[i-1]);
- }
-
- } else if(tofparampp) {
-
- for(Int_t i=1; i<=kPiBins_pp; i++) {
- hPiPi->SetBinContent(i,kPiTagPi_pp[i-1]);
- hPiNid->SetBinContent(i,kPiTagPi_pp[i-1]+kPiTagNid_pp[i-1]);
-
- hKK->SetBinContent(i,kKTagK_pp[i-1]);
- hKPi->SetBinContent(i,kKTagK_pp[i-1]+kKTagPi_pp[i-1]);
- hKNid->SetBinContent(i,kKTagK_pp[i-1]+kKTagPi_pp[i-1]+kKTagNid_pp[i-1]);
- }
- for(Int_t i=1; i<=kPBins_pp; i++) {
- hPP->SetBinContent(i,kPTagP_pp[i-1]);
- hPPi->SetBinContent(i,kPTagP_pp[i-1]+kPTagPi_pp[i-1]);
- hPNid->SetBinContent(i,kPTagP_pp[i-1]+kPTagPi_pp[i-1]+kPTagNid_pp[i-1]);
- }
-
- }
-
-
- TCanvas* c = new TCanvas("c","Parameterized PID in TOF",0,0,1000,400);
- c->Divide(3,1);
- c->cd(1);
- framePi->Draw();
- hPiNid->SetFillColor(18); hPiNid->Draw("same");
- hPiPi->SetFillColor(4); hPiPi->Draw("same");
- TPaveLabel* pav1 = new TPaveLabel(1,.2,1.4,.3,"#pi");
- pav1->SetBorderSize(0);
- pav1->Draw("same");
- TPaveLabel* pav2 = new TPaveLabel(1,.8,1.8,.9,"non-id");
- pav2->SetBorderSize(0);
- pav2->Draw("same");
-
- c->cd(2);
- frameK->Draw();
- hKNid->SetFillColor(18); hKNid->Draw("same");
- hKPi->SetFillColor(4); hKPi->Draw("same");
- hKK->SetFillColor(7); hKK->Draw("same");
- TPaveLabel* pav3 = new TPaveLabel(1,.2,1.5,.3,"K");
- pav3->SetBorderSize(0);
- pav3->Draw("same");
- TPaveLabel* pav4 = new TPaveLabel(1,.8,1.8,.9,"non-id");
- pav4->SetBorderSize(0);
- pav4->Draw("same");
- TPaveLabel* pav5 = new TPaveLabel(.4,.5,.8,.6,"#pi");
- pav5->SetBorderSize(0);
- pav5->Draw("same");
-
- c->cd(3);
- frameP->Draw();
- hPNid->SetFillColor(18); hPNid->Draw("same");
- hPPi->SetFillColor(4); hPPi->Draw("same");
- hPP->SetFillColor(3); hPP->Draw("same");
- TPaveLabel* pav6 = new TPaveLabel(1,.2,1.5,.3,"p");
- pav6->SetBorderSize(0);
- pav6->Draw("same");
- TPaveLabel* pav7 = new TPaveLabel(1,.8,2.6,.9,"non-id");
- pav7->SetBorderSize(0);
- pav7->Draw("same");
- TPaveLabel* pav8 = new TPaveLabel(.2,.5,1,.6,"#pi");
- pav8->SetBorderSize(0);
- pav8->Draw("same");
-
-
- return;
-}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff